Full-order Model Research Articles

Enhanced Drag Force Estimation in Automotive Design: A Surrogate Model Leveraging Limited Full-Order Model Drag Data and Comprehensive Physical Field Integration

In this paper, a novel surrogate model for shape-parametrized vehicle drag force prediction is proposed. It is assumed that only a limited dataset of high-fidelity CFD results is available, typically less than ten high-fidelity CFD solutions for different shape samples. The idea is to take advantage not only of the drag coefficients but also physical fields such as velocity, pressure, and kinetic energy evaluated on a cutting plane in the wake of the vehicle and perpendicular to the road. This additional “augmented” information provides a more accurate and robust prediction of the drag force compared to a standard surface response methodology. As a first step, an original reparametrization of the shape based on combination coefficients of shape principal components is proposed, leading to a low-dimensional representation of the shape space. The second step consists in determining principal components of the x-direction momentum flux through a cutting plane behind the car. The final step is to find the mapping between the reduced shape description and the momentum flux formula to achieve an accurate drag estimation. The resulting surrogate model is a space-parameter separated representation with shape principal component coefficients and spatial modes dedicated to drag-force evaluation. The algorithm can deal with shapes of variable mesh by using an optimal transport procedure that interpolates the fields on a shared reference mesh. The Machine Learning algorithm is challenged on a car concept with a three-dimensional shape design space. With only two well-chosen samples, the numerical algorithm is able to return a drag surrogate model with reasonable uniform error over the validation dataset. An incremental learning approach involving additional high-fidelity computations is also proposed. The leading algorithm is shown to improve the model accuracy. The study also shows the sensitivity of the results with respect to the initial experimental design. As feedback, we discuss and suggest what appear to be the correct choices of experimental designs for the best results.

A projection-based time-segmented reduced order model for fluid-structure interactions

In this paper, a type of novel projection-based, time-segmented reduced order model (ROM) is proposed for dynamic fluid-structure interaction (FSI) problems based upon the arbitrary Lagrangian–Eulerian (ALE)-finite element method (FEM) in a monolithic frame, where spatially, each variable is separated from others in terms of their attribution (fluid/structure), category (velocity/pressure) and component (two/three dimension) while temporally, the proper orthogonal decomposition (POD) bases are constructed in some deliberately partitioned time segments tailored through extensive numerical trials. By the combination of spatial and temporal decompositions, the developed ROM approach enables prolonged simulations under prescribed accuracy thresholds. Numerical experiments are carried out to compare numerical performances of the proposed ROM with corresponding full-order model (FOM) by solving a two-dimensional FSI benchmark problem that involves a vibrating elastic beam in the fluid, where the performance of offline ROM on perturbed physical parameters in the online phase is investigated as well. Extensive numerical results demonstrate that the proposed ROM has a comparable accuracy to while much higher efficiency than the FOM. The developed ROM approach is dimension-independent and can be seamlessly extended to solve high dimensional FSI problems.

Evaluating Reduced-Order Urban Wind Models for Simulating Flight Dynamics of Advanced Aerial Mobility Aircraft

Advanced Aerial Mobility (AAM) platforms are poised to begin high-density operations in urban areas nationwide. This new category of aviation platforms spans a broad range of sizes, from small package delivery drones to passenger-carrying vehicles. Unlike traditional aircraft, AAM vehicles operate within the urban boundary layer, where large structures, such as buildings, interrupt the flow. This study examines the response of a package delivery drone, a general aviation aircraft, and a passenger-carrying urban air mobility aircraft through an urban wind field generated using Large Eddy Simulations (LES). Since it is burdensome to simulate flight dynamics in real-time using the full-order solution, reduced-order wind models are created. Comparing trajectories for each aircraft platform using full-order or reduced-order solutions reveals little difference; reduced-order wind representations appear sufficient to replicate trajectories as long as the spatiotemporal wind field is represented. However, examining control usage statistics and time histories creates a stark difference between the wind fields, especially for the lower wing-loading package delivery drone where control saturation was encountered. The control saturation occurrences were inconsistent across the full-order and reduced-order winds, advising caution when using reduced-order models for lightly wing-loaded aircraft. The results presented demonstrate the effectiveness of using a simulation environment to evaluate reduced-order models by directly comparing their trajectories and control activity metrics with the full-order model. This evaluation provides designers valuable insights for making informed decisions for disturbance rejection systems. Additionally, the results indicate that using Reynolds-averaged Navier–Stokes (RANS) solutions to represent urban wind fields is inappropriate. It was observed that the mean wind field trajectories fall outside the 95% confidence intervals, a finding consistent with the authors’ previous research.

Nonlinear unsteady aerodynamic forces prediction and aeroelastic analysis of wind-induced bridge response at multiple wind speeds: A deep learning-based reduced-order model

Machine learning-based aerodynamic reduced-order models (ROMs) combine high accuracy with extremely low computational costs, making them highly effective in predicting nonlinear and unsteady bridge aerodynamic forces. Although several machine learning-based nonlinear aerodynamic models have been developed, the majority are built on a single wind speed parameter. However, in nonlinear aerodynamic prediction and aeroelastic analysis of bridges, the variability in incoming wind speed significantly influences the computed results. A ROM relying solely on a single wind speed lacks the ability to accurately forecast the intricate dynamic behaviors arising from changes in wind speed. When the incoming wind speed changes, the model's prediction accuracy significantly decreases. Usually, it is necessary to establish a new database and train a new model, which not only increases time and cost but also greatly reduces the convenience of the ROM. Addressing this challenge, this study proposes a multiple-wind-speed (MWS) nonlinear unsteady bridge aerodynamic model based on the LSTM deep neural network. Taking the Taohuayu Yellow River Bridge in the Henan Province of China as an example, the modeling process of the proposed MWS-ROM is demonstrated, along with non-linear aerodynamic predictions of the deck under various conditions and aerodynamic-elastic analysis of the deck under different wind speeds. The research results show that the trained LSTM network can accurately predict the nonlinear aerodynamic forces of bridges under single and double degrees of freedom vibration conditions. The MWS-ROM performed well in predicting convergent vibrations at low wind speeds and limits cycle oscillations at high wind speeds, aligning closely with results from the CFD full-order model. Compared to CFD, the aerodynamic ROM based on the LSTM network significantly enhances computational efficiency, consequently boosting the convenience and efficiency of bridge flutter analysis. Additionally, the methodology proposed herein can be extended for wind-induced vibration control and response prediction in other types of deck sections.

Impact of Phase Angle Jump on a Doubly Fed Induction Generator under Low-Voltage Ride-Through Based on Transfer Function Decomposition

During the fault period, a phase angle jump may occur at the stator or the point of common coupling, which will deteriorate the low-voltage ride-through (LVRT) characteristics of a doubly fed induction generator (DFIG). The existing LVRT studies focus on the impact of a voltage drop on DFIGs but often ignore that of a phase angle jump. The time-domain simulation is accurate in describing the response of a DFIG during the LVRT process, but it is time-consuming for a DFIG with the full-order model. In this paper, by using the voltage magnitude and phase angle of the stator or the point of common coupling as the inputs, and the state variables as the outputs, the transfer function of a DFIG is derived to analyze its response and find the LVRT measures against the voltage drop and, especially, the phase angle jump. Firstly, the differential-algebraic equations of the DFIG are linearized to propose their transfer function model. Secondly, considering its high-order characteristic, a model reduction method for the transfer function of the DFIG using the Schur decomposition is proposed, and the analytical expression of the output variables of the DFIG with the phase angle jump is derived by the inverse Laplace transformation to judge the necessity of the LVRT measures. Finally, the simulation results of the DFIG are provided to verify the accuracy of the transfer function model and its reduced-order form and validate the feasibility of the LVRT against the phase angle jump with the proposed models.

Research on Model Reduction of AUV Underwater Support Platform Based on Digital Twin

Digital twin technology, as a data-driven and model-driven innovation means, plays a crucial role in the process of digital transformation and intelligent upgrading of the marine industry, helping the industry to move towards a new stage of more intelligent and efficient development. In order to solve the defects of the Autonomous Underwater Vehicle (AUV) underwater support platform structure deformation field, digital twin technology and model reduction technology are applied to an AUV underwater support platform, and a five-dimensional digital twin model of the AUV underwater support platform is studied, including five dimensions: physical world, digital world, twin data center, service application, and data connection. The digital twin of the subsea support platform is established by using the digital twin modeling technology. The POD method is used to calculate the deformation field matrix of the support structure of the subsea support platform under the 0–5 sea state, and the corresponding eigenvalues and eigenvectors are obtained. By intercepting the eigenvectors corresponding to the eigenvalues of the high energy proportion, the low-order equation is constructed, and the reduced-order model under each sea state can be quickly solved. The experimental results show that the model reduction technology can greatly shorten the model solving time, and the calculated results are highly consistent with the simulation results of the finite element full-order model, which can realize the rapid analysis of the deformation response of the subsea support platform structure, and provide a theoretical basis and technical support for the subsequent simulation, state evaluation, visual monitoring, and predictive maintenance.

Dynamic Effects Analysis in Fractional Memristor-Based Rulkov Neuron Model.

Mathematical models such as Fitzhugh-Nagoma and Hodgkin-Huxley models have been used to understand complex nervous systems. Still, due to their complexity, these models have made it challenging to analyze neural function. The discrete Rulkov model allows the analysis of neural function to facilitate the investigation of neuronal dynamics or others. This paper introduces a fractional memristor Rulkov neuron model and analyzes its dynamic effects, investigating how to improve neuron models by combining discrete memristors and fractional derivatives. These improvements include the more accurate generation of heritable properties compared to full-order models, the treatment of dynamic firing activity at multiple time scales for a single neuron, and the better performance of firing frequency responses in fractional designs compared to integer models. Initially, we combined a Rulkov neuron model with a memristor and evaluated all system parameters using bifurcation diagrams and the 0-1 chaos test. Subsequently, we applied a discrete fractional-order approach to the Rulkov memristor map. We investigated the impact of all parameters and the fractional order on the model and observed that the system exhibited various behaviors, including tonic firing, periodic firing, and chaotic firing. We also found that the more I tend towards the correct order, the more chaotic modes in the range of parameters. Following this, we coupled the proposed model with a similar one and assessed how the fractional order influences synchronization. Our results demonstrated that the fractional order significantly improves synchronization. The results of this research emphasize that the combination of memristor and discrete neurons provides an effective tool for modeling and estimating biophysical effects in neurons and artificial neural networks.

An efficient reduced order model for nonlinear transient porous media flow with time-varying injection rates

An intrusive Reduced Order Model (ROM) is developed for nonlinear porous media flow problems with transient and time-discontinuous fluid injection rates. The proposed ROM is significantly more computationally efficient than the Full Order Model (FOM). The training regime is generated using the FOM with constant injection rates during the offline stage. The trained ROM exhibits high accuracy for complex pumping schedules (rate vs time) simulated online. The proposed ROM uses the combination of Proper Orthogonal Decomposition and Discrete Empirical Interpolation Method (POD-DEIM), which is compared with the classical POD-Galerkin. The use of an approximated column-reduced Jacobian is shown to be vital to achieving a substantial speedup of ROM vs FOM run-times. An analysis of the trade-off between accuracy and run-time is conducted for ROMs of different sizes and hyper-parameters. The impact of the training regime on the performance of the presented ROM is assessed. The performance of the ROM is studied in the context of a two-dimensional analysis of time-varying injection into a two-well system in a layered porous media reservoir. The accuracy and efficiency of POD-DEIM motivate its potential use as a surrogate model in the real-time control and monitoring of fluid injection processes.

A Simplified 4-DOF Dynamic Model of a Series-Parallel Hybrid Electric Vehicle

To research the dynamic response of a new type of dedicated transmission for a hybrid electric vehicle, a detailed dynamics model should be built. However, a model with too many degrees of freedom has a negative effect on controller design, which means the detailed model should be simplified. In this paper, two dynamic models are established. One is an original and detailed powertrain dynamics model (ODPDM), which can capture the transient response, and it is validated that the ODPDM can be used to accurately describe the real vehicle in some specific operating conditions. The other is a simplified torsional vibration dynamics model to study the torsional vibration characteristics of the hybrid electric vehicle. Compared with the full-order model, which is based on the ODPDM, the simplified model has a very similar vibration in low frequency. This study provides a basis for further vibration control of the hybrid powertrain during the process of a driving-mode switch.

Fast model calibration for predicting the response of breast cancer to chemotherapy using proper orthogonal decomposition

Constructing digital twins for predictive tumor treatment response models can have a high computational demand that presents a practical barrier for their clinical adoption. In this work, we demonstrate that proper orthogonal decomposition, by which a low-dimensional representation of the full model is constructed, can be used to dramatically reduce the computational time required to calibrate a partial differential equation model to magnetic resonance imaging (MRI) data for rapid predictions of tumor growth and response to chemotherapy. In the proposed formulation, the reduction basis is based on each patient’s own MRI data and controls the overall size of the “reduced order model”. Using the full model as the reference, we validate that the reduced order mathematical model can accurately predict response in 50 triple negative breast cancer patients receiving standard of care neoadjuvant chemotherapy. The concordance correlation coefficient between the full and reduced order models was 0.986 ± 0.012 (mean ± standard deviation) for predicting changes in both tumor volume and cellularity across the entire model family, with a corresponding median local error (inter-quartile range) of 4.36 % (1.22 %, 15.04 %). The total time to estimate parameters and to predict response dramatically improves with the reduced framework. Specifically, the reduced order model accelerates our calibration by a factor of (mean ± standard deviation) 378.4 ± 279.8 when compared to the full order model for a non-mechanically coupled model. This enormous reduction in computational time can directly help realize the practical construction of digital twins when the access to computational resources is limited.

Second-order Arnoldi accelerated boundary element method for two-dimensional broadband acoustic shape sensitivity analysis

This paper proposes a novel approach for broadband acoustic shape sensitivity analysis based on the direct differentiation approach. Since the system matrices of the boundary element method (BEM) for the analysis of acoustic state and acoustic sensitivity have frequency dependence, repeated calculations are needed at different frequencies. This is very time-consuming, especially for sensitivity calculations used in shape optimization design. The Taylor series expansion of the Hankel function is carried out to separate the frequency-dependent and frequency-independent terms in the acoustic shape sensitivity boundary integral equation to construct a frequency-independent system matrix. In addition, due to the formation of asymmetric full-coefficient matrices in acoustic shape sensitivity equations based on the BEM, repeatedly solving system equations is also extremely time-consuming at broadband frequencies for large scale issues. The second-order Arnoldi approach was employed to create a reduced-order model that maintains the key features of the initial full-order model. The strong singular and supersingular integrals within the sensitivity equations can be calculated directly utilizing the singularity elimination technique. Finally, several numerical examples confirm the accuracy and efficiency of the proposed algorithm.

Efficient aerodynamic design using BézierGAN and model order reduction: A computational study

Bézier Generative Adversarial Networks (BézierGANs) is one of the most widely used Convolutional Neural Network (CNN) technique-based algorithms for generating smooth curves to analyze and optimize complex shapes and flow patterns, particularly aerodynamics. Airfoils, the basic cross-sectional shape to design air wings, are generated here by selecting the vital parameters wisely after alternating the typical BézierGANs algorithm instead of relying on pre-existing models. When these predicted airfoils are used to develop the physical structure of the adjacent air layer to observe the phenomena, they become the subject of computationally expensive simulations of turbulent airflow. Simulations of turbulent airflow that are computationally expensive are required when these anticipated airfoils are utilized to create the physical structure of the adjacent air layer to monitor the physical events. To prevail over this computational burden, the Model Order Reduction (MOR) technique is applied to develop the best-fitted lower-order model of the original turbulent air layer model around the airfoils by filtering the most effective system's singular values while preserving all the essential characteristics of the original large-scale model. As evidence of the effectiveness of the Reduced Order Model (ROM) of the turbulent air layer regarding the truncation of the design time and computational costs, the comparative analysis of the computational time is occupied by the original and prescribed system to ensure about 99.7660% of less time consumption by the ROM. The regeneration of the reduced-order physical model is also covered in this paper, which analyzes the physical behaviors regarding the air pressure and the velocity patterns of the airflow around the airfoils and compares the similarities to the physical characteristics of the Full Order Model (FOM). It was found that the ROM can be one of the best alternatives to FOM in real-life implementation.

Parameter identification and uncertainty propagation of hydrogel coupled diffusion-deformation using POD-based reduced-order modeling

AbstractThis study explores reduced-order modeling for analyzing time-dependent diffusion-deformation of hydrogels. The full-order model describing hydrogel transient behavior consists of a coupled system of partial differential equations in which the chemical potential and displacements are coupled. This system is formulated in a monolithic fashion and solved using the finite element method. We employ proper orthogonal decomposition as a model order reduction approach. The reduced-order model performance is tested through a benchmark problem on hydrogel swelling and a case study simulating co-axial printing. Then, we embed the reduced-order model into an optimization loop to efficiently identify the coupled problem’s material parameters using full-field data. Finally, a study is conducted on the uncertainty propagation of the material parameter.

Reduced-order modeling of unsteady fluid flow using neural network ensembles

The use of deep learning has become increasingly popular in reduced-order models (ROMs) to obtain low-dimensional representations of full-order models. Convolutional autoencoders (CAEs) are often used to this end as they are adept at handling data that are spatially distributed, including solutions to partial differential equations. When applied to unsteady physics problems, ROMs also require a model for time-series prediction of the low-dimensional latent variables. Long short-term memory (LSTM) networks, a type of recurrent neural network useful for modeling sequential data, are frequently employed in data-driven ROMs for autoregressive time-series prediction. When making predictions at unseen design points over long time horizons, error propagation is a frequently encountered issue, where errors made early on can compound over time and lead to large inaccuracies. In this work, we propose using bagging, a commonly used ensemble learning technique, to develop a fully data-driven ROM framework referred to as the CAE-eLSTM ROM that uses CAEs for spatial reconstruction of the full-order model and LSTM ensembles for time-series prediction. When applied to two unsteady fluid dynamics problems, our results show that the presented framework effectively reduces error propagation and leads to more accurate time-series prediction of latent variables at unseen points.

Projection-based reduced-order modelling of time-periodic problems, with application to flow past flapping hydrofoils

Simulating forced time-periodic flows in industrial applications presents significant computational challenges, partly due to the need to overcome costly transients before achieving time-periodicity. Reduced-order modelling emerges as a promising method to speed-up computations. We extend upon the work of Lotz et al. (2024) where a time-periodic space–time model is introduced. We present a time-periodic reduced-order model that directly finds the time-periodic solution without requiring extensive time integration. The reduced-order model gives a reduction in variables in both space and time. Our approach involves a POD-Galerkin reduced-order model based on a time-periodic full-order model that employs isogeometric analysis, residual-based variational multiscale turbulence modelling and weak boundary conditions. The projection-based reduced-order model inherits these features. We evaluate the reduced-order model with numerical experiments on moving hydrofoils. The motion is known a priori and we restrict ourselves to two spatial dimensions. In these experiments we vary the Strouhal and Reynolds numbers, and the motion profile respectively. Reduced-order model solutions agree well with those of the full-order model. The errors over the entire time period of thrust and lift forces are less than 0.2%. This includes complex scenarios such as the transition from drag to thrust production with increasing Strouhal number. Our time-periodic reduced-order model offers speed-ups ranging from O(102) to O(103) compared to the full-order model, depending upon the basis size. This makes it an appealing solution for prescribed time-periodic problems, with potential for additional speedup through nonlinear reduction techniques such as hyper-reduction.

Non-intrusive reduced order models for partitioned fluid–structure interactions

The main goal of this work is to develop a data-driven Reduced Order Model (ROM) strategy from high-fidelity simulation result data of a Full Order Model (FOM). The goal is to predict at lower computational cost the time evolution of solutions of Fluid–Structure Interaction (FSI) problems. For some FSI applications, the elastic solid FOM (often chosen as quasi-static) can take far more computational time than the fluid one. In this context, for the sake of performance one could only derive a ROM for the structure and try to achieve a partitioned FOM fluid solver coupled with a ROM solid one. In this paper, we present a data-driven partitioned ROM on two study cases: (i) a simplified 1D-1D FSI problem representing an axisymmetric elastic model of an arterial vessel, coupled with an incompressible fluid flow; (ii) an incompressible 2D wake flow over a cylinder facing an elastic solid with two flaps. We evaluate the accuracy and performance of the proposed ROM-FOM strategy on these cases while investigating the effects of the model’s hyperparameters. We demonstrate a high prediction accuracy and significant speedup achievements using this strategy.

A reduced order model formulation for left atrium flow: an atrial fibrillation case

A data-driven reduced order model (ROM) based on a proper orthogonal decomposition-radial basis function (POD-RBF) approach is adopted in this paper for the analysis of blood flow dynamics in a patient-specific case of atrial fibrillation (AF). The full order model (FOM) is represented by incompressible Navier–Stokes equations, discretized with a finite volume (FV) approach. Both the Newtonian and the Casson’s constitutive laws are employed. The aim is to build a computational tool able to efficiently and accurately reconstruct the patterns of relevant hemodynamics indices related to the stasis of the blood in a physical parametrization framework including the cardiac output in the Newtonian case and also the plasma viscosity and the hematocrit in the non-Newtonian one. Many FOM-ROM comparisons are shown to analyze the performance of our approach as regards errors and computational speed-up.

Data-Driven Nonintrusive Model-Order Reduction for Aerodynamic Design Optimization

Fast and accurate evaluation of aerodynamic characteristics is essential for aerodynamic design optimization because aircraft programs require many years of design and optimization. Therefore, it is imperative to develop sufficiently fast, robust, and accurate computational tools for industry routine analysis. This paper presents a nonintrusive machine-learning method for building reduced-order models (ROMs) using an autoencoder neural network architecture. An optimization framework was developed to identify the optimal solution by exploring the low-dimensional subspace generated by the trained autoencoder. To demonstrate the convergence, stability, and reliability of the ROM, a subsonic inverse design problem and a transonic drag minimization problem of the airfoil were studied and validated using two different parameterization strategies. The robustness and accuracy demonstrated by the method suggest that it is valuable in parametric studies, such as aerodynamic design and optimization, and requires only a small fraction of the cost of full-order modeling.

Uncertainty qualification in seismic analysis of concrete dams based on model order reduction accelerated stochastic SBFEM

A pioneering framework for uncertainty quantification in seismic analysis of concrete dams rooted in the stochastic scaled boundary finite element method (SSBFEM) is proposed. For seismic analysis of dams, octree decomposition is used for mesh generation. SSBFEM is developed to predict the structural responses with randomly distributed material properties. The script in this paper can integrate octree mesh and SSBFEM into ABAQUS user element (UEL). Implementing octree SSBFEM on commercial software provides accurate and convenient full order models for uncertainty quantification. Monte Carlo simulation (MCs) is deployed to calculate the statistical characteristics of structural responses under random variables. Additionally, singular value decomposition (SVD) and radial basis function (RBF) are leveraged to refine traditional MCs. The solution space of MCs is decomposed into lower-order subspaces. Any structural responses can be rapidly evaluated from linear combinations of subspaces. Particularly, the method has been successfully applied to seismic analysis of concrete dams with different random material properties. Finally, several illustrative examples are provided to validate the accuracy and efficacy of the algorithm.

Geometrically parametrised reduced order models for studying the hysteresis of the Coanda effect in finite element-based incompressible fluid dynamics

This article presents a general reduced order model (ROM) framework for addressing fluid dynamics problems involving time-dependent geometric parametrisations. The framework integrates Proper Orthogonal Decomposition (POD) and Empirical Cubature Method (ECM) hyper-reduction techniques to effectively approximate incompressible computational fluid dynamics simulations. To demonstrate the applicability of this framework, we investigate the behaviour of a planar contraction-expansion channel geometry exhibiting bifurcating solutions known as the Coanda effect. By introducing time-dependent deformations to the channel geometry, we observe hysteresis phenomena in the solution.The paper provides a detailed formulation of the framework, including the stabilised finite elements full order model (FOM) and ROM, with a particular focus on the considerations related to geometric parametrisation. Subsequently, we present the results obtained from the simulations, analysing the solution behaviour in a phase space for the fluid velocity at a probe point, considered as the Quantity of Interest (QoI). Through qualitative and quantitative evaluations of the ROMs and hyper-reduced order models (HROMs), we demonstrate their ability to accurately reproduce the complete solution field and the QoI.While HROMs offer significant computational speedup, enabling efficient simulations, they do exhibit some errors, particularly for testing trajectories. However, their value lies in applications where the detection of the Coanda effect holds paramount importance, even if the selected bifurcation branch is incorrect. Alternatively, for more precise results, HROMs with lower speedups can be employed.