Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational...Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.展开更多
This paper describes a method proposed for modeling large deflection of aircraft in nonlinear aeroelastic analysis by developing reduced order model（ROM）.The method is applied for solving the static aeroelastic and ...This paper describes a method proposed for modeling large deflection of aircraft in nonlinear aeroelastic analysis by developing reduced order model（ROM）.The method is applied for solving the static aeroelastic and static aeroelastic trim problems of flexible aircraft containing geometric nonlinearities;meanwhile,the non-planar effects of aerodynamics and follower force effect have been considered.ROMs are computational inexpensive mathematical representations compared to traditional nonlinear finite element method（FEM） especially in aeroelastic solutions.The approach for structure modeling presented here is on the basis of combined modal/finite element（MFE） method that characterizes the stiffness nonlinearities and we apply that structure modeling method as ROM to aeroelastic analysis.Moreover,the non-planar aerodynamic force is computed by the non-planar vortex lattice method（VLM）.Structure and aerodynamics can be coupled with the surface spline method.The results show that both of the static aeroelastic analysis and trim analysis of aircraft based on structure ROM can achieve a good agreement compared to analysis based on the FEM and experimental result.展开更多
Based on Recursive Radial Basis Function(RRBF)neural network,the Reduced Order Model(ROM)of compressor cascade was established to meet the urgent demand of highly efficient prediction of unsteady aerodynamics performa...Based on Recursive Radial Basis Function(RRBF)neural network,the Reduced Order Model(ROM)of compressor cascade was established to meet the urgent demand of highly efficient prediction of unsteady aerodynamics performance of turbomachinery.One novel ROM called ASA-RRBF model based on Adaptive Simulated Annealing(ASA)algorithm was developed to enhance the generalization ability of the unsteady ROM.The ROM was verified by predicting the unsteady aerodynamics performance of a highly-loaded compressor cascade.The results show that the RRBF model has higher accuracy in identification of the dimensionless total pressure and dimensionless static pressure of compressor cascade under nonlinear and unsteady conditions,and the model behaves higher stability and computational efficiency.However,for the strong nonlinear characteristics of aerodynamic parameters,the RRBF model presents lower accuracy.Additionally,the RRBF model predicts with a large error in the identification of aerodynamic parameters under linear and unsteady conditions.For ASA-RRBF,by introducing a small-amplitude and highfrequency sinusoidal signal as validation sample,the width of the basis function of the RRBF model is optimized to improve the generalization ability of the ROM under linear unsteady conditions.Besides,this model improves the predicting accuracy of dimensionless static pressure which has strong nonlinear characteristics.The ASA-RRBF model has higher prediction accuracy than RRBF model without significantly increasing the total time consumption.This novel model can predict the linear hysteresis of dimensionless static pressure happened in the harmonic condition,but it cannot accurately predict the beat frequency of dimensionless total pressure.展开更多
Active stability augmentation system is an attractive and promising technology to suppress flutter and limit cycle oscillation (LCO). In order to design a good active control law, the control plant model with low orde...Active stability augmentation system is an attractive and promising technology to suppress flutter and limit cycle oscillation (LCO). In order to design a good active control law, the control plant model with low order and high accuracy must be provided, which is one of the most important key points. The traditional model is based on low fidelity aerodynamics model such as panel method, which is unsuitable for transonic flight regime. The physics-based high fidelity tools, reduced order model (ROM) and CFD/CSD coupled aeroservoelastic solver are used to design the active control law. The Volterra/ROM is applied to constructing the low order state space model for the nonlinear unsteady aerodynamics and static output feedback method is used to active control law design. The detail of the new method is demonstrated by the Goland+ wing/store system. The simulation results show that the effectiveness of the designed active augmentation system, which can suppress the flutter and LCO successfully.展开更多
This paper presents a design of continuous-time sliding mode control for the higher order systems via reduced order model. It is shown that a continuous-time sliding mode control designed for the reduced order model g...This paper presents a design of continuous-time sliding mode control for the higher order systems via reduced order model. It is shown that a continuous-time sliding mode control designed for the reduced order model gives similar performance for thc higher order system. The method is illustrated by numerical examples. The paper also introduces a technique for design of a sliding surface such that the system satisfies a cost-optimality condition when on the sliding surface.展开更多
A Reduced Order Model(ROM)based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper.The unsteady aerodynamic model is established by a system identification technique combined wi...A Reduced Order Model(ROM)based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper.The unsteady aerodynamic model is established by a system identification technique combined with a set of Aerodynamic Influence Coefficients(AIC).Subsequently,the aerodynamic model is encoded into the state space and then coupled with the structural dynamic equations,resulting in a ROM of the cascade aeroelasticity.The cascade flutter can be determined by solving the eigenvalues of the ROM.Bending-torsional coupled mode flutter analysis for the Standard Configuration Eleven(SC11)cascade is used to validate the proposed method.展开更多
This paper considers the problem of simulating the humidity distributions of a grain storage system. The distributions are described by partial differential equations(PDE). It is quite difficult to obtain the humidity...This paper considers the problem of simulating the humidity distributions of a grain storage system. The distributions are described by partial differential equations(PDE). It is quite difficult to obtain the humidity profiles from the PDE model. Hence, a discretization method is applied to obtain an equivalent ordinary differential equation model. However, after applying the discretization technique, the cost of solving the system increases as the size increases to a few thousands. It may be noted that after discretization,the degree of freedom of the system remain the same while the order increases. The large dynamic model is reduced using a proper orthogonal decomposition based technique and an equivalent model but of much reduced size is obtained. A controller based on optimal control theory is designed to obtain an input such that the output humidity reaches a desired profile and also its stability is analyzed.Numerical results are presented to show the validity of the reduced model and possible further extensions are identified.展开更多
Aerothermoelasticity is one of the key technologies for hypersonic vehicles. Accurate and efficient computation of the aerothermodynamics is one of the primary challenges for hypersonic aerothermoelastic analysis. Aim...Aerothermoelasticity is one of the key technologies for hypersonic vehicles. Accurate and efficient computation of the aerothermodynamics is one of the primary challenges for hypersonic aerothermoelastic analysis. Aimed at solving the shortcomings of engineering calculation, compu- tation fluid dynamics （CFD） and experimental investigation, a reduced order modeling （ROM） framework for aerothermodynamics based on CFD predictions using an enhanced algorithm of fast maximin Latin hypercube design is developed. Both proper orthogonal decomposition （POD） and surrogate are considered and compared to construct ROMs. Two surrogate approaches named Kriging and optimized radial basis function （ORBF） are utilized to construct ROMs. Furthermore, an enhanced algorithm of fast maximin Latin hypercube design is proposed, which proves to be helpful to improve the precisions of ROMs. Test results for the three-dimensional aerothermody- namic over a hypersonic surface indicate that： the ROMs precision based on Kriging is better than that by ORBF, ROMs based on Kriging are marginally more accurate than ROMs based on POD- Kriging. In a word, the ROM framework for hypersonic aerothermodynamics has good precision and efficiency.展开更多
Reduced order models(ROMs) based on the snapshots on the CFD high-fidelity simulations have been paid great attention recently due to their capability of capturing the features of the complex geometries and flow confi...Reduced order models(ROMs) based on the snapshots on the CFD high-fidelity simulations have been paid great attention recently due to their capability of capturing the features of the complex geometries and flow configurations. To improve the efficiency and precision of the ROMs, it is indispensable to add extra sampling points to the initial snapshots, since the number of sampling points to achieve an adequately accurate ROM is generally unknown in prior, but a large number of initial sampling points reduces the parsimony of the ROMs. A fuzzy-clustering-based adding-point strategy is proposed and the fuzzy clustering acts an indicator of the region in which the precision of ROMs is relatively low. The proposed method is applied to construct the ROMs for the benchmark mathematical examples and a numerical example of hypersonic aerothermodynamics prediction for a typical control surface. The proposed method can achieve a 34.5% improvement on the efficiency than the estimated mean squared error prediction algorithm and shows same-level prediction accuracy.展开更多
Flow around a cavity is characterized by a self-sustained mechanism in which the shear layer impinges on the downstream edge of the cavity resulting in a feedback mechanism.Direct Numerical Simulations of the flow at ...Flow around a cavity is characterized by a self-sustained mechanism in which the shear layer impinges on the downstream edge of the cavity resulting in a feedback mechanism.Direct Numerical Simulations of the flow at low Reynolds number has been carried out to get pressure and velocity fluctuations,for the case of un-actuated and multi frequency actuation.A Reduced Order Model for the isentropic compressible equations based on the method of Proper Orthogonal Decomposition has been constructed.The model has been extended to include the effect of control.The Reduced Order dynamical system shows a divergence in time integration.A method of calibration based on the minimization of a linear functional of error,to the sensitivity of the modes,is proposed.The calibrated low order model is used to design a feedback control of cavity flows based on an observer design.For the experimental implementation of the controller,a state estimate based on the observed pressure measurements is obtained through a linear stochastic estimation.Finally the obtained control is introduced into the Direct Numerical Simulation to obtain a decrease in spectra of the cavity acoustic mode.展开更多
This paper discusses a physics-informed methodology aimed at reconstructing efficiently the fluid state of a system.Herein,the generation of an accurate reduced order model of twodimensional unsteady flows from data l...This paper discusses a physics-informed methodology aimed at reconstructing efficiently the fluid state of a system.Herein,the generation of an accurate reduced order model of twodimensional unsteady flows from data leverages on sparsity-promoting statistical learning techniques.The cornerstone of the approach is l_(1) regularised regression,resulting in sparselyconnected models where only the important quadratic interactions between modes are retained.The original dynamical behaviour is reproduced at low computational costs,as few quadratic interactions need to be evaluated.The approach has two key features.First,interactions are selected systematically as a solution of a convex optimisation problem and no a priori assumptions on the physics of the flow are required.Second,the presence of a regularisation term improves the predictive performance of the original model,generally affected by noise and poor data quality.Test cases are for two-dimensional lid-driven cavity flows,at three values of the Reynolds number for which the motion is chaotic and energy interactions are scattered across the spectrum.It is found that:(A)the sparsification generates models maintaining the original accuracy level but with a lower number of active coefficients;this becomes more pronounced for increasing Reynolds numbers suggesting that extension of these techniques to real-life flow configurations is possible;(B)sparse models maintain a good temporal stability for predictions.The methodology is ready for more complex applications without modifications of the underlying theory,and the integration into a cyberphysical model is feasible.展开更多
A computationally efficient two-surface plasticity model is assessed against crystal plasticity. Focus is laid on the mechanical behavior of magnesium alloys in the presence of ductility-limiting defects, such as void...A computationally efficient two-surface plasticity model is assessed against crystal plasticity. Focus is laid on the mechanical behavior of magnesium alloys in the presence of ductility-limiting defects, such as voids. The two surfaces separately account for slip and twinning such that the constitutive formulation captures the evolving plastic anisotropy and evolving tension-compression asymmetry. For model identification, a procedure is proposed whereby the initial guess is based on a combination of experimental data and computationally intensive polycrystal calculations from the literature. In drawing direct comparisons with crystal plasticity, of which the proposed model constitutes a heuristically derived reduced-order model, the available crystal plasticity simulations are grouped in two datasets. A calibration set contains minimal data for both pristine and porous material subjected to one loading path. Then the two-surface model is assessed against a broader set of crystal plasticity simulations for voided unit cells under various stress states and two loading orientations. The assessment also includes microstructure evolution(rate of growth of porosity and void distortion). The ability of the two-surface model to capture essential features of crystal plasticity is analyzed along with an evaluation of computational cost. The prospects of using the model in guiding the development of physically sound damage models in Mg alloys are put forth in the context of high-throughput simulations.展开更多
The fast and accurate reduced-order modeling of fluidized beds is a challenging task in the field of fluid dynamics,owing to their high dimensionality and nonlinear dynamic behavior.In this study,a nonintrusive reduce...The fast and accurate reduced-order modeling of fluidized beds is a challenging task in the field of fluid dynamics,owing to their high dimensionality and nonlinear dynamic behavior.In this study,a nonintrusive reduced order modeling method,the reduced order model based on principal component analysis and bidirectional long short-term memory networks(PBLSTM ROM),was developed to capture complex spatio-temporal dynamics of fluidized beds.By combining principal component analysis and Bidirectional long-short-term memory networks,the PBLSTM ROM effectively extracted dynamic evolution information without any prior knowledge of governing equations,enabling reduced-order modeling of unsteady flow fields.The PBLSTM ROM was validated using the solid volume fraction and gas velocity flow fields of a fluidized bed with immersed tubes,showing superior performance over both the PLSTM and PANN ROMs in accurately capturing temporal changes in the fluidization fields,especially in the region near immersed tubes where severe fluctuations appear.Moreover,the PBLSTM ROM improved the simulation speed by five orders of magnitude compared to traditional computational fluid dynamics simulations.These findings suggest that the PBLSTM ROM presents a promising approach for analyzing the complex fluid flows in engineering practice.展开更多
Accurate and efficient online parameter identification and state estimation are crucial for leveraging digital twin simulations to optimize the operation of near-carbon-free nuclear energy systems.In previous studies,...Accurate and efficient online parameter identification and state estimation are crucial for leveraging digital twin simulations to optimize the operation of near-carbon-free nuclear energy systems.In previous studies,we developed a reactor operation digital twin(RODT).However,non-differentiabilities and discontinuities arise when employing machine learning-based surrogate forward models,challenging traditional gradient-based inverse methods and their variants.This study investigated deterministic and metaheuristic algorithms and developed hybrid algorithms to address these issues.An efficient modular RODT software framework that incorporates these methods into its post-evaluation module is presented for comprehensive comparison.The methods were rigorously assessed based on convergence profiles,stability with respect to noise,and computational performance.The numerical results show that the hybrid KNNLHS algorithm excels in real-time online applications,balancing accuracy and efficiency with a prediction error rate of only 1%and processing times of less than 0.1 s.Contrastingly,algorithms such as FSA,DE,and ADE,although slightly slower(approximately 1 s),demonstrated higher accuracy with a 0.3%relative L_2 error,which advances RODT methodologies to harness machine learning and system modeling for improved reactor monitoring,systematic diagnosis of off-normal events,and lifetime management strategies.The developed modular software and novel optimization methods presented offer pathways to realize the full potential of RODT for transforming energy engineering practices.展开更多
Dynamic mode decomposition(DMD) aims at extracting intrinsic mechanisms in a time sequence via linear recurrence relation of its observables, thereby predicting later terms in the sequence. Stability is a major concer...Dynamic mode decomposition(DMD) aims at extracting intrinsic mechanisms in a time sequence via linear recurrence relation of its observables, thereby predicting later terms in the sequence. Stability is a major concern in DMD predictions. We adopt a regularized form and propose a Regularized DMD(Re DMD) algorithm to determine the regularization parameter. This leverages stability and accuracy. Numerical tests for Burgers' equation demonstrate that Re DMD effectively stabilizes the DMD prediction while maintaining accuracy. Comparisons are made with the truncated DMD algorithm.展开更多
Gust alleviation is very important to a large flexible aircraft.A nonlinear low-order aerodynamic state space model is required to model the nonlinear aeroelastic responses due to gust.Based on the proper orthogonal d...Gust alleviation is very important to a large flexible aircraft.A nonlinear low-order aerodynamic state space model is required to model the nonlinear aeroelastic responses due to gust.Based on the proper orthogonal decomposition method,a reduced order modeling of gust loads was proposed.And then the open-loop and closed-loop reduced order state space model for the transonic aeroelastic system was developed.The static output feed back control scheme was used to design a simple multiple-in multiple-out(MIMO)gust alleviation control law.The control law was demonstrated with the Goland+wing model with four control surfaces.The simulation results of different discrete gusts show the capability and good performance of the designed MIMO controller in transonic gust alleviation.展开更多
The improved line sampling （LS） technique, an effective numerical simulation method, is employed to analyze the probabilistic characteristics and reliability sensitivity of flutter with random structural parameter i...The improved line sampling （LS） technique, an effective numerical simulation method, is employed to analyze the probabilistic characteristics and reliability sensitivity of flutter with random structural parameter in transonic flow. The improved LS technique is a novel methodology for reliability and sensitivity analysis of high dimensionality and low probability problem with implicit limit state function, and it does not require any approximating surrogate of the implicit limit state equation. The improved LS is used to estimate the flutter reliability and the sensitivity of a two-dimensional wing, in which some structural properties, such as frequency, parameters of gravity center and mass ratio, are considered as random variables. Computational fluid dynamics （CFD） based unsteady aerodynamic reduced order model （ROM） method is used to construct the aerodynamic state equations. Coupling structural state equations with aerodynamic state equations, the safety margin of flutter is founded by using the critical velocity of flutter. The results show that the improved LS technique can effectively decrease the computational cost in the random uncertainty analysis of flutter. The reliability sensitivity, defined by the partial derivative of the failure probability with respect to the distribution parameter of random variable, can help to identify the important parameters and guide the structural optimization design.展开更多
The hydraulic fracturing is a nonlinear,fluid-solid coupling and transient problem,in most cases it is always time-consuming to simulate this process numerically.In recent years,although many numerical methods were pr...The hydraulic fracturing is a nonlinear,fluid-solid coupling and transient problem,in most cases it is always time-consuming to simulate this process numerically.In recent years,although many numerical methods were proposed to settle this problem,most of them still require a large amount of computer resources.Thus it is a high demand to develop more efficient numerical approaches to achieve the real-time monitoring of the fracture geometry during the hydraulic fracturing treatment.In this study,a reduced order modeling technique namely Proper Generalized Decomposition(PGD),is applied to accelerate the simulations of the transient,non-linear coupled system of hydraulic fracturing problem,to match this extremely tight response time constraint.The separability of the solution in space and time dimensions is studied for a simplified model problem.The solid and fluid equations are coupled explicitly by inverting the solid discrete problem,and a simple iterative procedure to handle the non-linear characteristic of the hydraulic fracturing problem is proposed in this work.Numeral validation illustrates that the results of PGD match well with these of standard finite element method in terms o f fracture opening and fluid pressure in the hydro-fracture.Moreover,after the off-line calculations,the numerical results can be obtained in real time.展开更多
The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model(ROM) framework f...The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model(ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit(GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation(LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.展开更多
A projection-based reduced order model(ROM)based on the Fourier collocation method is proposed for compressible flows.The incorporation of localized artificial viscosity model and filtering is pursued to enhance the r...A projection-based reduced order model(ROM)based on the Fourier collocation method is proposed for compressible flows.The incorporation of localized artificial viscosity model and filtering is pursued to enhance the robustness and accuracy of the ROM for shock-dominated flows.Furthermore,for Euler systems,ROMs built on the conservative and the skew-symmetric forms of the governing equation are compared.To ensure efficiency,the discrete empirical interpolation method(DEIM)is employed.An alternative reduction approach,exploring the sparsity of viscosity is also investigated for the viscous terms.A number of one-and two-dimensional benchmark cases are considered to test the performance of the proposed models.Results show that stable computations for shock-dominated cases can be achieved with ROMs built on both the conservative and the skew-symmetric forms without additional stabilization components other than the viscosity model and filtering.Under the same parameters,the skew-symmetric form shows better robustness and accuracy than its conservative counterpart,while the conservative form is superior in terms of efficiency.展开更多
基金supported by the National Key R&D Program of China under Grant No.2021ZD0110400.
文摘Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.
文摘This paper describes a method proposed for modeling large deflection of aircraft in nonlinear aeroelastic analysis by developing reduced order model（ROM）.The method is applied for solving the static aeroelastic and static aeroelastic trim problems of flexible aircraft containing geometric nonlinearities;meanwhile,the non-planar effects of aerodynamics and follower force effect have been considered.ROMs are computational inexpensive mathematical representations compared to traditional nonlinear finite element method（FEM） especially in aeroelastic solutions.The approach for structure modeling presented here is on the basis of combined modal/finite element（MFE） method that characterizes the stiffness nonlinearities and we apply that structure modeling method as ROM to aeroelastic analysis.Moreover,the non-planar aerodynamic force is computed by the non-planar vortex lattice method（VLM）.Structure and aerodynamics can be coupled with the surface spline method.The results show that both of the static aeroelastic analysis and trim analysis of aircraft based on structure ROM can achieve a good agreement compared to analysis based on the FEM and experimental result.
基金co-National Science and Technology Major Project(No.2017-II-0009-0023)Innovation Guidance Support Project for Taicang Top Research Institutes(No.TC2019DYDS09)。
文摘Based on Recursive Radial Basis Function(RRBF)neural network,the Reduced Order Model(ROM)of compressor cascade was established to meet the urgent demand of highly efficient prediction of unsteady aerodynamics performance of turbomachinery.One novel ROM called ASA-RRBF model based on Adaptive Simulated Annealing(ASA)algorithm was developed to enhance the generalization ability of the unsteady ROM.The ROM was verified by predicting the unsteady aerodynamics performance of a highly-loaded compressor cascade.The results show that the RRBF model has higher accuracy in identification of the dimensionless total pressure and dimensionless static pressure of compressor cascade under nonlinear and unsteady conditions,and the model behaves higher stability and computational efficiency.However,for the strong nonlinear characteristics of aerodynamic parameters,the RRBF model presents lower accuracy.Additionally,the RRBF model predicts with a large error in the identification of aerodynamic parameters under linear and unsteady conditions.For ASA-RRBF,by introducing a small-amplitude and highfrequency sinusoidal signal as validation sample,the width of the basis function of the RRBF model is optimized to improve the generalization ability of the ROM under linear unsteady conditions.Besides,this model improves the predicting accuracy of dimensionless static pressure which has strong nonlinear characteristics.The ASA-RRBF model has higher prediction accuracy than RRBF model without significantly increasing the total time consumption.This novel model can predict the linear hysteresis of dimensionless static pressure happened in the harmonic condition,but it cannot accurately predict the beat frequency of dimensionless total pressure.
基金National Natural Science Foundation of China (10902082)New Faculty Research Foundation of XJTUthe Fundamental Research Funds for the Central Universities (xjj20100126)
文摘Active stability augmentation system is an attractive and promising technology to suppress flutter and limit cycle oscillation (LCO). In order to design a good active control law, the control plant model with low order and high accuracy must be provided, which is one of the most important key points. The traditional model is based on low fidelity aerodynamics model such as panel method, which is unsuitable for transonic flight regime. The physics-based high fidelity tools, reduced order model (ROM) and CFD/CSD coupled aeroservoelastic solver are used to design the active control law. The Volterra/ROM is applied to constructing the low order state space model for the nonlinear unsteady aerodynamics and static output feedback method is used to active control law design. The detail of the new method is demonstrated by the Goland+ wing/store system. The simulation results show that the effectiveness of the designed active augmentation system, which can suppress the flutter and LCO successfully.
文摘This paper presents a design of continuous-time sliding mode control for the higher order systems via reduced order model. It is shown that a continuous-time sliding mode control designed for the reduced order model gives similar performance for thc higher order system. The method is illustrated by numerical examples. The paper also introduces a technique for design of a sliding surface such that the system satisfies a cost-optimality condition when on the sliding surface.
基金supported by the National Science and Technology Major Project, China (No. 2017-II-0009-0023)the Aeronautical Science Foundation of China(No. 2020Z039053004)the Fundamental Research Funds for the Central Universities, China (No. 3102019OQD701)
文摘A Reduced Order Model(ROM)based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper.The unsteady aerodynamic model is established by a system identification technique combined with a set of Aerodynamic Influence Coefficients(AIC).Subsequently,the aerodynamic model is encoded into the state space and then coupled with the structural dynamic equations,resulting in a ROM of the cascade aeroelasticity.The cascade flutter can be determined by solving the eigenvalues of the ROM.Bending-torsional coupled mode flutter analysis for the Standard Configuration Eleven(SC11)cascade is used to validate the proposed method.
文摘This paper considers the problem of simulating the humidity distributions of a grain storage system. The distributions are described by partial differential equations(PDE). It is quite difficult to obtain the humidity profiles from the PDE model. Hence, a discretization method is applied to obtain an equivalent ordinary differential equation model. However, after applying the discretization technique, the cost of solving the system increases as the size increases to a few thousands. It may be noted that after discretization,the degree of freedom of the system remain the same while the order increases. The large dynamic model is reduced using a proper orthogonal decomposition based technique and an equivalent model but of much reduced size is obtained. A controller based on optimal control theory is designed to obtain an input such that the output humidity reaches a desired profile and also its stability is analyzed.Numerical results are presented to show the validity of the reduced model and possible further extensions are identified.
基金supported by the National Natural Science Foundation of China (Nos. 11372036, 50875024)Excellent Young Scholars Research Fund of Beijing Institute of Technology of China (No. 2010Y0102)
文摘Aerothermoelasticity is one of the key technologies for hypersonic vehicles. Accurate and efficient computation of the aerothermodynamics is one of the primary challenges for hypersonic aerothermoelastic analysis. Aimed at solving the shortcomings of engineering calculation, compu- tation fluid dynamics （CFD） and experimental investigation, a reduced order modeling （ROM） framework for aerothermodynamics based on CFD predictions using an enhanced algorithm of fast maximin Latin hypercube design is developed. Both proper orthogonal decomposition （POD） and surrogate are considered and compared to construct ROMs. Two surrogate approaches named Kriging and optimized radial basis function （ORBF） are utilized to construct ROMs. Furthermore, an enhanced algorithm of fast maximin Latin hypercube design is proposed, which proves to be helpful to improve the precisions of ROMs. Test results for the three-dimensional aerothermody- namic over a hypersonic surface indicate that： the ROMs precision based on Kriging is better than that by ORBF, ROMs based on Kriging are marginally more accurate than ROMs based on POD- Kriging. In a word, the ROM framework for hypersonic aerothermodynamics has good precision and efficiency.
基金Supported by National Natural Science Foundation of China(Grant No.11372036)
文摘Reduced order models(ROMs) based on the snapshots on the CFD high-fidelity simulations have been paid great attention recently due to their capability of capturing the features of the complex geometries and flow configurations. To improve the efficiency and precision of the ROMs, it is indispensable to add extra sampling points to the initial snapshots, since the number of sampling points to achieve an adequately accurate ROM is generally unknown in prior, but a large number of initial sampling points reduces the parsimony of the ROMs. A fuzzy-clustering-based adding-point strategy is proposed and the fuzzy clustering acts an indicator of the region in which the precision of ROMs is relatively low. The proposed method is applied to construct the ROMs for the benchmark mathematical examples and a numerical example of hypersonic aerothermodynamics prediction for a typical control surface. The proposed method can achieve a 34.5% improvement on the efficiency than the estimated mean squared error prediction algorithm and shows same-level prediction accuracy.
基金supported by a Marie Curie Early Stage Research Training Fellowship of the European Community’s Sixth Framework Programme under contract number MEST CT 2005020301.
文摘Flow around a cavity is characterized by a self-sustained mechanism in which the shear layer impinges on the downstream edge of the cavity resulting in a feedback mechanism.Direct Numerical Simulations of the flow at low Reynolds number has been carried out to get pressure and velocity fluctuations,for the case of un-actuated and multi frequency actuation.A Reduced Order Model for the isentropic compressible equations based on the method of Proper Orthogonal Decomposition has been constructed.The model has been extended to include the effect of control.The Reduced Order dynamical system shows a divergence in time integration.A method of calibration based on the minimization of a linear functional of error,to the sensitivity of the modes,is proposed.The calibrated low order model is used to design a feedback control of cavity flows based on an observer design.For the experimental implementation of the controller,a state estimate based on the observed pressure measurements is obtained through a linear stochastic estimation.Finally the obtained control is introduced into the Direct Numerical Simulation to obtain a decrease in spectra of the cavity acoustic mode.
文摘This paper discusses a physics-informed methodology aimed at reconstructing efficiently the fluid state of a system.Herein,the generation of an accurate reduced order model of twodimensional unsteady flows from data leverages on sparsity-promoting statistical learning techniques.The cornerstone of the approach is l_(1) regularised regression,resulting in sparselyconnected models where only the important quadratic interactions between modes are retained.The original dynamical behaviour is reproduced at low computational costs,as few quadratic interactions need to be evaluated.The approach has two key features.First,interactions are selected systematically as a solution of a convex optimisation problem and no a priori assumptions on the physics of the flow are required.Second,the presence of a regularisation term improves the predictive performance of the original model,generally affected by noise and poor data quality.Test cases are for two-dimensional lid-driven cavity flows,at three values of the Reynolds number for which the motion is chaotic and energy interactions are scattered across the spectrum.It is found that:(A)the sparsification generates models maintaining the original accuracy level but with a lower number of active coefficients;this becomes more pronounced for increasing Reynolds numbers suggesting that extension of these techniques to real-life flow configurations is possible;(B)sparse models maintain a good temporal stability for predictions.The methodology is ready for more complex applications without modifications of the underlying theory,and the integration into a cyberphysical model is feasible.
基金support of this work by the National Science Foundation (CMMI Award no.1932975)。
文摘A computationally efficient two-surface plasticity model is assessed against crystal plasticity. Focus is laid on the mechanical behavior of magnesium alloys in the presence of ductility-limiting defects, such as voids. The two surfaces separately account for slip and twinning such that the constitutive formulation captures the evolving plastic anisotropy and evolving tension-compression asymmetry. For model identification, a procedure is proposed whereby the initial guess is based on a combination of experimental data and computationally intensive polycrystal calculations from the literature. In drawing direct comparisons with crystal plasticity, of which the proposed model constitutes a heuristically derived reduced-order model, the available crystal plasticity simulations are grouped in two datasets. A calibration set contains minimal data for both pristine and porous material subjected to one loading path. Then the two-surface model is assessed against a broader set of crystal plasticity simulations for voided unit cells under various stress states and two loading orientations. The assessment also includes microstructure evolution(rate of growth of porosity and void distortion). The ability of the two-surface model to capture essential features of crystal plasticity is analyzed along with an evaluation of computational cost. The prospects of using the model in guiding the development of physically sound damage models in Mg alloys are put forth in the context of high-throughput simulations.
基金supported by the National Key R&D Program of China(grant No.2021YFF0500400)Key Research Program of Shaanxi Province(grant No.2022GXLH-01-08)+2 种基金National Key R&D Program of China(grant No.2018YFB1501003)Shaanxi Province Qin Chuangyuan“Scientist+Engineer”Team(grant No.2022KXJ-179)Targeted Funding Program of Power Construction Corporation of China(grant No.DJ-PTZX-2021-03).
文摘The fast and accurate reduced-order modeling of fluidized beds is a challenging task in the field of fluid dynamics,owing to their high dimensionality and nonlinear dynamic behavior.In this study,a nonintrusive reduced order modeling method,the reduced order model based on principal component analysis and bidirectional long short-term memory networks(PBLSTM ROM),was developed to capture complex spatio-temporal dynamics of fluidized beds.By combining principal component analysis and Bidirectional long-short-term memory networks,the PBLSTM ROM effectively extracted dynamic evolution information without any prior knowledge of governing equations,enabling reduced-order modeling of unsteady flow fields.The PBLSTM ROM was validated using the solid volume fraction and gas velocity flow fields of a fluidized bed with immersed tubes,showing superior performance over both the PLSTM and PANN ROMs in accurately capturing temporal changes in the fluidization fields,especially in the region near immersed tubes where severe fluctuations appear.Moreover,the PBLSTM ROM improved the simulation speed by five orders of magnitude compared to traditional computational fluid dynamics simulations.These findings suggest that the PBLSTM ROM presents a promising approach for analyzing the complex fluid flows in engineering practice.
基金supported by the Natural Science Foundation of Shanghai(No.23ZR1429300)Innovation Funds of CNNC(Lingchuang Fund,Contract No.CNNC-LCKY-202234)the Project of the Nuclear Power Technology Innovation Center of Science Technology and Industry(No.HDLCXZX-2023-HD-039-02)。
文摘Accurate and efficient online parameter identification and state estimation are crucial for leveraging digital twin simulations to optimize the operation of near-carbon-free nuclear energy systems.In previous studies,we developed a reactor operation digital twin(RODT).However,non-differentiabilities and discontinuities arise when employing machine learning-based surrogate forward models,challenging traditional gradient-based inverse methods and their variants.This study investigated deterministic and metaheuristic algorithms and developed hybrid algorithms to address these issues.An efficient modular RODT software framework that incorporates these methods into its post-evaluation module is presented for comprehensive comparison.The methods were rigorously assessed based on convergence profiles,stability with respect to noise,and computational performance.The numerical results show that the hybrid KNNLHS algorithm excels in real-time online applications,balancing accuracy and efficiency with a prediction error rate of only 1%and processing times of less than 0.1 s.Contrastingly,algorithms such as FSA,DE,and ADE,although slightly slower(approximately 1 s),demonstrated higher accuracy with a 0.3%relative L_2 error,which advances RODT methodologies to harness machine learning and system modeling for improved reactor monitoring,systematic diagnosis of off-normal events,and lifetime management strategies.The developed modular software and novel optimization methods presented offer pathways to realize the full potential of RODT for transforming energy engineering practices.
基金supported by the National Nature Science Foundation of China (Grant No.11988102)National Undergraduate Training Program for Innovation and Entrepreneurship。
文摘Dynamic mode decomposition(DMD) aims at extracting intrinsic mechanisms in a time sequence via linear recurrence relation of its observables, thereby predicting later terms in the sequence. Stability is a major concern in DMD predictions. We adopt a regularized form and propose a Regularized DMD(Re DMD) algorithm to determine the regularization parameter. This leverages stability and accuracy. Numerical tests for Burgers' equation demonstrate that Re DMD effectively stabilizes the DMD prediction while maintaining accuracy. Comparisons are made with the truncated DMD algorithm.
基金supported by the National Natural Science Foundation of China(Grant Nos.11272005,10902082,91016008)
文摘Gust alleviation is very important to a large flexible aircraft.A nonlinear low-order aerodynamic state space model is required to model the nonlinear aeroelastic responses due to gust.Based on the proper orthogonal decomposition method,a reduced order modeling of gust loads was proposed.And then the open-loop and closed-loop reduced order state space model for the transonic aeroelastic system was developed.The static output feed back control scheme was used to design a simple multiple-in multiple-out(MIMO)gust alleviation control law.The control law was demonstrated with the Goland+wing model with four control surfaces.The simulation results of different discrete gusts show the capability and good performance of the designed MIMO controller in transonic gust alleviation.
基金Foundation items： National Natural Science Foundation of China （NSFC 10572117, 10802063, 50875213） National High-tech Research and Development Program （2007AA04Z401）+2 种基金 Aeronautical Science Foundation of China （2007ZA53012） New Century Program For Excellent Talents of Ministry of Education of China （NCET-05-0868） Ph.D. Program Foundation of Northwestern Polytechnical University （CX200801）.
文摘The improved line sampling （LS） technique, an effective numerical simulation method, is employed to analyze the probabilistic characteristics and reliability sensitivity of flutter with random structural parameter in transonic flow. The improved LS technique is a novel methodology for reliability and sensitivity analysis of high dimensionality and low probability problem with implicit limit state function, and it does not require any approximating surrogate of the implicit limit state equation. The improved LS is used to estimate the flutter reliability and the sensitivity of a two-dimensional wing, in which some structural properties, such as frequency, parameters of gravity center and mass ratio, are considered as random variables. Computational fluid dynamics （CFD） based unsteady aerodynamic reduced order model （ROM） method is used to construct the aerodynamic state equations. Coupling structural state equations with aerodynamic state equations, the safety margin of flutter is founded by using the critical velocity of flutter. The results show that the improved LS technique can effectively decrease the computational cost in the random uncertainty analysis of flutter. The reliability sensitivity, defined by the partial derivative of the failure probability with respect to the distribution parameter of random variable, can help to identify the important parameters and guide the structural optimization design.
基金the National Science Foundation of China(Grant Nos.51804033 and 51936001)China Postdoctoral Science and Foundation(Grant No.2018M641254)+3 种基金Beijing Postdoctoral Research Foundation(2018-ZZ-045)the Project of Construction of Innovative Teams and Teacher Career Development for Universities and Colleges Under Beijing Municipality(Grant No.IDHT20170507)Program of Great Wall Scholar(Grant No.CIT&TCD20180313)Jointly Projects of Beijing Natural Science Foundation and Beijing Municipal Education Commission(Grant No.KZ201810017023).
文摘The hydraulic fracturing is a nonlinear,fluid-solid coupling and transient problem,in most cases it is always time-consuming to simulate this process numerically.In recent years,although many numerical methods were proposed to settle this problem,most of them still require a large amount of computer resources.Thus it is a high demand to develop more efficient numerical approaches to achieve the real-time monitoring of the fracture geometry during the hydraulic fracturing treatment.In this study,a reduced order modeling technique namely Proper Generalized Decomposition(PGD),is applied to accelerate the simulations of the transient,non-linear coupled system of hydraulic fracturing problem,to match this extremely tight response time constraint.The separability of the solution in space and time dimensions is studied for a simplified model problem.The solid and fluid equations are coupled explicitly by inverting the solid discrete problem,and a simple iterative procedure to handle the non-linear characteristic of the hydraulic fracturing problem is proposed in this work.Numeral validation illustrates that the results of PGD match well with these of standard finite element method in terms o f fracture opening and fluid pressure in the hydro-fracture.Moreover,after the off-line calculations,the numerical results can be obtained in real time.
基金supported by the National Natural Science Foundation of China(No.11672018).
文摘The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model(ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit(GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation(LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.
基金supported by the National Natural Science Foundation of China(grant numbers 11972064 and 11721202).
文摘A projection-based reduced order model(ROM)based on the Fourier collocation method is proposed for compressible flows.The incorporation of localized artificial viscosity model and filtering is pursued to enhance the robustness and accuracy of the ROM for shock-dominated flows.Furthermore,for Euler systems,ROMs built on the conservative and the skew-symmetric forms of the governing equation are compared.To ensure efficiency,the discrete empirical interpolation method(DEIM)is employed.An alternative reduction approach,exploring the sparsity of viscosity is also investigated for the viscous terms.A number of one-and two-dimensional benchmark cases are considered to test the performance of the proposed models.Results show that stable computations for shock-dominated cases can be achieved with ROMs built on both the conservative and the skew-symmetric forms without additional stabilization components other than the viscosity model and filtering.Under the same parameters,the skew-symmetric form shows better robustness and accuracy than its conservative counterpart,while the conservative form is superior in terms of efficiency.
