Model reduction of fractional impedance spectra for time–frequency analysis of batteries, fuel cells, and supercapacitors

Joint time–frequency analysis is an emerging method for interpreting the underlying physics in fuel cells,batteries,and supercapacitors.To increase the reliability of time–frequency analysis,a theoretical correlation between frequency-domain stationary analysis and time-domain transient analysis is urgently required.The present work formularizes a thorough model reduction of fractional impedance spectra for electrochemical energy devices involving not only the model reduction from fractional-order models to integer-order models and from high-to low-order RC circuits but also insight into the evolution of the characteristic time constants during the whole reduction process.The following work has been carried out:(i)the model-reduction theory is addressed for typical Warburg elements and RC circuits based on the continued fraction expansion theory and the response error minimization technique,respectively;(ii)the order effect on the model reduction of typical Warburg elements is quantitatively evaluated by time–frequency analysis;(iii)the results of time–frequency analysis are confirmed to be useful to determine the reduction order in terms of the kinetic information needed to be captured;and(iv)the results of time–frequency analysis are validated for the model reduction of fractional impedance spectra for lithium-ion batteries,supercapacitors,and solid oxide fuel cells.In turn,the numerical validation has demonstrated the powerful function of the joint time–frequency analysis.The thorough model reduction of fractional impedance spectra addressed in the present work not only clarifies the relationship between time-domain transient analysis and frequency-domain stationary analysis but also enhances the reliability of the joint time–frequency analysis for electrochemical energy devices.

Statistical detection of structural damage based on model reduction

This paper proposes a statistical method for damage detection based on the finite element （FE） model reduction technique that utilizes measured modal data with a limited number of sensors. A deterministic damage detection process is formulated based on the model reduction technique. The probabilistie process is integrated into the deterministic damage detection process using a perturbation technique, resulting in a statistical structural damage detection method. This is achieved by deriving the first- and second-order partial derivatives of uncertain parameters, such as elasticity of the damaged member, with respect to the measurement noise, which allows expectation and covariance matrix of the uncertain parameters to be calculated. Besides the theoretical development, this paper reports numerical verification of the proposed method using a portal frame example and Monte Carlo simulation.

Model reduction and active control of flexible beam using internal balance technique

The internal balance technique is effective for the model reduction in flexible structures, especially the ones with dense frequencies. However, due to the difficulty in extracting the internal balance modal coordinates from the physical sensor readings, research on this topic has been mostly theoretical so far, and little has been done in experiments or engineering applications. This paper studies the internal balance method theoretically as well as experimentally and designs an active controller based on the reduction model. The research works on a digital signal processor （DSP） TMS320F2812- based experiment system with a flexible beam and proposes an approximate approach to access the internal balance modal coordinates. The simulation and test results have shown that the proposed approach is feasible and effective, and the designed controller is successful in restraining the beam vibration.

Stable Routh-Padé-type approximation in model reduction of interval systems

An interval Pade-type approximation is introduced and then Routh-Pade-type method （IRPTM） is presented to model reduction in interval systems. The denominator in reduced model is obtained from the stable Routh table, and its numerator is constructed by the interval Pade-type definition. Compared to the existing Routh-Pade method, IRPTM does not need to solve linear interval equations theoretical analysis shows that IRPTM has example is given to illustrate our method. Hence, we do not have to compute smaller computational cost than that interval division in the process. Moreover, of Routh-Pade method. A typical numerical

Static output feedback stabilization for second-order singular systems using model reduction methods

In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.

Model reduction using the genetic algorithmand routh approxi mations

A new method of model reduction combining the genetic algorithm（GA） with the Routh approximation method is presented. It is suggested that a high-order system can be approximated by a low-order model with a time delay. The denominator parameters of the reduced-order model are determined by the Routh approximation method, then the numerator parameters and time delay are identified by the GAL. The reduced-order models obtained by the proposed method will always be stable if the original system is stable and produce a good approximation to the original system in both the frequency domain and time domain. Two numerical examples show that the method is cornputationally simple and efficient.

A matrix Padé-type-Routh model reduction method for multivariable linear systems

In the present paper, the formulae for matrix Padé-type approximation were improved. The mixed model reduction method of matrix Padé-type-Routh for the multivariable linear systems was presented. A well-known example was given to illustrate that the mixed method is efficient.

Model Reduction Methods for High-Rise Buildings with Active Mass Damper Control Systems

To improve the performance of an active mass damper control system,the controller should be designed based on a reduced-order model. An improved method based on balanced truncation method was proposed to reduce the dimension of high-rise buildings,and was compared with other widely used reduction methods by using a framework with ten floors. This optimized method has improvement of reduction process and choice of the order. Based on the reduced-order model obtained by the improved method and pole-assignment algorithm,a controller was designed. Finally,a comparative analysis of structural responses,transfer functions,and poles was conducted on an actual high-rise building. The results show the effectiveness of the improved method.

EMTP-type Program Realization of Krylov Subspace Based Model Reduction Methods for Large-Scale Active Distribution Network

Effective model reduction methods are required to deal with new challenges in active distribution network simulations that are on a large scale and have complicated structures.In the development of advanced electromagnetic transient simulation programs,automated model reduction plays an important role.This paper proposes an automated realization algorithm for the Krylov subspace based model reduction methods of an active distribution network with which the reduced model can be automatically established according to a given threshold of reduction error.The combined state-space nodal analysis framework is employed to apply the automated model reduction algorithm in popular EMTP-type simulation programs.Simulations are performed using PSCAD and a self-developed program to show the feasibility and validity of the proposed methods.

MODEL REDUCTION AND ACTIVE CONTROL FOR A FLEXIBLE PLATE

The internal balance technique is effective for model reduction in flexible structures, especially those with dense frequencies. However, due to the difficulty in extracting the internal balance modal coordinates from the physical sensor readings, research so far on this topic has been mostly theoretic and little on experiment or engineering applications. This paper, by working on a DSP TMS320F2812-based experiment system with a flexible plate and bringing forward an approximating approach to accessing the internal balance modal coordinates, studies the internal balance method theoretically as well as experimentally, and further designs an active controller based on the reduced model. Simulation and test results have proven the proposed approximating approach feasible and effective, and the designed controller successful in restraining the plate vibration.

Hilbert-Schmidt-Hankel norm model reduction for matrix second-order linear systems

This paper considers the optimal model reduction problem of matrix second-order linear systems in the sense of Hilbert-Schmidt-Hankel norm, with the reduced order systems preserving the structure of the original systems. The expressions of the error function and its gradient are derived. Two numerical examples are given to illustrate the presented model reduction technique.

Model Reduction and Controller Design for a Nonlinear Heat Conduction Problem Using Finite Element Method

The mathematical models for dynamic distributed parameter systems are given by systems of partial differential equations. With nonlinear material properties, the corresponding finite element (FE) models are large systems of nonlinear ordinary differential equations. However, in most cases, the actual dynamics of interest involve only a few of the variables, for which model reduction strategies based on system theoretical concepts can be immensely useful. This paper considers the problem of controlling a three dimensional profile on nontrivial geometries. Dynamic model is obtained by discretizing the domain using FE method. A nonlinear control law is proposed which transfers any arbitrary initial temperature profile to another arbitrary desired one. The large dynamic model is reduced using proper orthogonal decomposition (POD). Finally, the stability of the control law is proved through Lyapunov analysis. Results of numerical implementation are presented and possible further extensions are identified.

Correlation analysis of PCB and comparison of test-analysis model reduction methods

The validity of correlation analysis between finite element model（FEM） and modal test data is strongly affected by three factors, i.e., quality of excitation and measurement points in modal test,FEM reduction methods, and correlation check techniques. A new criterion based on modified mode participation（MMP） for choosing the best excitation point is presented. Comparison between this new criterion and mode participation（MP） criterion is made by using Case 1 with a simple printed circuit board（PCB）. The result indicates that this new criterion produces better results. In Case 2, 35 measurement points are selected to perform modal test and correlation analysis while 9 selected in Case 3.System equivalent reduction expansion process（SEREP）, modal assurance criteria（MAC）, coordinate modal assurance criteria（CoMAC）, pseudo orthogonality check（POC） and coordinate orthogonality check（CORTHOG） are used to show the error introduced by modal test in Cases 2 and 3. Case 2 shows that additional errors which cannot be identified by using CoMAC can be found by using CORTHOG.In both Cases 2 and 3, Guyan reduction, improved reduced system（IRS） method, SEREP and Hybrid reduction are compared for accuracy and robustness. The results suggest that the quality of the reduction process is problem dependent. However, the IRS method is an improvement over the Guyan reduction, and the Hybrid reduction is an improvement over the SEREP reduction.

Development of stability-preserving time-limited model reduction framework for 2-D and 1-D models with error bound

Due to their complex structure,2-D models are challenging to work with;additionally,simulation,analysis,design,and control get increasingly difficult as the order of the model grows.Moreover,in particular time intervals,Gawronski and Juang’s time-limited model reduction schemes produce an unstable reduced-order model for the 2-D and 1-D models.Researchers revealed some stability preservation solutions to address this key flaw which ensure the stability of 1-D reduced-order systems;nevertheless,these strategies result in large approximation errors.However,to the best of the authors’knowledge,there is no literature available for the stability preserving time-limited-interval Gramian-based model reduction framework for the 2-D discrete-time systems.In this article,2-D models are decomposed into two separate sub-models(i.e.,two cascaded 1-D models)using the condition of minimal rank-decomposition.Model reduction procedures are conducted on these obtained two 1-D sub-models using limited-time Gramian.The suggested methodology works for both 2-D and 1-D models.Moreover,the suggested methodology gives the stability of the reduced model as well as a priori error-bound expressions for the 2-D and 1-D models.Numerical results and comparisons between existing and suggested methodologies are provided to demonstrate the effectiveness of the suggested methodology.

An Improved Model Reduction Method on AIMs for N-S Equations Using Multilevel Finite Element Method and Hierarchical Basis

A numerical method is proposed to approach the Approximate Inertial Man-ifolds(AIMs)in unsteady incompressible Navier-Stokes equations,using multilevel fi-nite element method with hierarchical basis functions.Following AIMS,the unknown variables,velocity and pressure in the governing equations,are divided into two com-ponents,namely low modes and high modes.Then,the couplings between low modes and high modes,which are not accounted by standard Galerkin method,are consid-ered by AIMs,to improve the accuracy of the numerical results.Further,the multilevel finite element method with hierarchical basis functions is introduced to approach low modes and high modes in an efficient way.As an example,the flow around airfoil NACA0012 at different angles of attack has been simulated by the method presented,and the comparisons show that there is a good agreement between the present method and experimental results.In particular,the proposed method takes less computing time than the traditional method.As a conclusion,the present method is efficient in numer-ical analysis of fluid dynamics,especially in computing time.

Frequency-domain generelaized singular peruturbation method for relative error model order reduction

A new mixed method for relative error model order reduction is proposed. In the proposed method the frequency domain balanced stochastic truncation method is improved by applying the generalized singular perturbation method to the frequency domain balanced system in the reduction procedure. The frequency domain balanced stochastic truncation method, which was proposed in [15] and [17] by the author, is based on two recently developed methods, namely frequency domain balanced truncation within a desired frequency bound and inner-outer factorization techniques. The proposed method in ttiis paper is a carry over of the frequency-domain balanced stochastic truncation and is of interest for practical model order reduction because in this context it shows to keep the accuracy of the approximation as high as possible without sacrificing the computational efficiency and important system properties. It is shown that some important properties of the frequency domain stochastic balanced reduction technique are extended to the proposed reduction method by using the concept and properties of the reciprocal systems. Numerical results show the accuracy, simplicity and flexibility enhancement of the method.

A new dynamic bacterial foraging optimization and its application on model reduction

Inspired by the foraging behavior of E.coli bacteria,bacterial foraging optimization(BFO)has emerged as a powerful technique for solving optimization problems.However,BFO shows poor performance on complex and high-dimensional optimization problems.In order to improve the performance of BFO,a new dynamic bacterial foraging optimization based on clonal selection(DBFO-CS)is proposed.Instead of fixed step size in the chemotaxis operator,a new piecewise strategy adjusts the step size dynamically by regulatory factor in order to balance between exploration and exploitation during optimization process,which can improve convergence speed.Furthermore,reproduction operator based on clonal selection can add excellent genes to bacterial populations in order to improve bacterial natural selection and help good individuals to be protected,which can enhance convergence precision.Then,a set of benchmark functions have been used to test the proposed algorithm.The results show that DBFO-CS offers significant improvements than BFO on convergence,accuracy and robustness.A complex optimization problem of model reduction on stable and unstable linear systems based on DBFO-CS is presented.Results show that the proposed algorithm can efficiently approximate the systems.

Comparison of Invariant Manifolds for Model Reduction in Chemical Kinetics

A modern approach to model reduction in chemical kinetics is often based on the notion of slow invariant manifold.The goal of this paper is to give a comparison of various methods of construction of slow invariant manifolds using a simple Michaelis-Menten catalytic reaction.We explore a recently introduced Method of Invariant Grids(MIG)for iteratively solving the invariance equation.Various initial approximations for the grid are considered such as Quasi Equilibrium Manifold,Spectral Quasi Equilibrium Manifold,Intrinsic Low Dimensional Manifold and Symmetric Entropic Intrinsic Low Dimensional Manifold.Slow invariant manifold was also computed using the Computational Singular Perturbation(CSP)method.A comparison between MIG and CSP is also reported.

The Role of Thermodynamics in Model Reduction when Using Invariant Grids

In the present work,we develop in detail the process leading to reduction of models in chemical kinetics when using the Method of Invariant Grids(MIG).To this end,reduced models(invariant grids)are obtained by refining initial approximations of slow invariant manifolds,and used for integrating smaller and less stiff systems of equations capable to recover the detailed description with high accuracy.Moreover,we clarify the role played by thermodynamics in model reduction,and carry out a comparison between detailed and reduced solutions for a model hydrogen oxidation reaction.

Multi-Objective Adaptive Optimization Model Predictive Control:Decreasing Carbon Emissions from a Zinc Oxide Rotary Kiln

The zinc oxide rotary kiln,as an essential piece of equipment in the zinc smelting industrial process,is presenting new challenges in process control.China’s strategy of achieving a carbon peak and carbon neutrality is putting new demands on the industry,including green production and the use of fewer resources;thus,traditional stability control is no longer suitable for multi-objective control tasks.Although researchers have revealed the principle of the rotary kiln and set up computational fluid dynamics(CFD)simulation models to study its dynamics,these models cannot be directly applied to process control due to their high computational complexity.To address these issues,this paper proposes a multi-objective adaptive optimization model predictive control(MAO-MPC)method based on sparse identification.More specifically,with a large amount of data collected from a CFD model,a sparse regression problem is first formulated and solved to obtain a reduction model.Then,a two-layered control framework including real-time optimization(RTO)and model predictive control(MPC)is designed.In the RTO layer,an optimization problem with the goal of achieving optimal operation performance and the lowest possible resource consumption is set up.By solving the optimization problem in real time,a suitable setting value is sent to the MPC layer to ensure that the zinc oxide rotary kiln always functions in an optimal state.Our experiments show the strength and reliability of the proposed method,which reduces the usage of coal while maintaining high profits.