Reduced Order Model (ROM)

Category: Analysis | Integrated 2026-04-06
reduced-order-model-theory
Theory & Physics

Reduced Order Model (ROM): Theoretical Foundations

Overview

๐Ÿง‘โ€๐ŸŽ“

Teacher! Today's topic is about Reduced Order Models (ROM), right? What are they?


๐ŸŽ“

It's a method that drastically reduces the degrees of freedom of high-dimensional CAE models using Proper Orthogonal Decomposition (POD) or Dynamic Mode Decomposition (DMD). It enables real-time simulation and parametric analysis.



Governing Equations


๐ŸŽ“

Expressing this mathematically, it looks like this.


$$\mathbf{u}(\mathbf{x},t) \approx \bar{\mathbf{u}}(\mathbf{x}) + \sum_{i=1}^{r} a_i(t) \boldsymbol{\phi}_i(\mathbf{x})$$

๐Ÿง‘โ€๐ŸŽ“

Hmm, just the equation alone doesn't really click... What does it represent?


๐ŸŽ“

Construction of POD basis:



$$\mathbf{C}\boldsymbol{\phi}_i = \lambda_i \boldsymbol{\phi}_i, \quad \mathbf{C} = \frac{1}{N}\mathbf{X}\mathbf{X}^T$$

Theoretical Foundation

๐Ÿง‘โ€๐ŸŽ“

I've heard of "theoretical foundation," but I might not fully understand it...


๐ŸŽ“

Reduced Order Models (ROM) are an important technique aiming for the fusion of data-driven approaches and physics-based modeling. While computational cost is a major bottleneck in traditional CAE analysis, introducing Reduced Order Models (ROM) can significantly improve the trade-off between computational efficiency and prediction accuracy. The mathematical foundation of this method is based on function approximation theory and statistical learning theory, with theoretical research topics including guarantees of generalization performance and rigorous analysis of convergence. Particularly, dealing with the "curse of dimensionality" when the input dimension is high is a key practical challenge, and approaches like dimensionality reduction and leveraging sparsity are important.


๐Ÿง‘โ€๐ŸŽ“

I see... ROMs seem simple at first glance, but they're actually very profound.


Details of Mathematical Formulation

๐Ÿง‘โ€๐ŸŽ“

Next is "Details of Mathematical Formulation"! What is this about?


๐ŸŽ“

It shows the basic mathematical framework for applying machine learning models to CAE.



Loss Function Composition

๐Ÿง‘โ€๐ŸŽ“

What does "loss function composition" mean specifically?


๐ŸŽ“

The loss function in AIร—CAE is composed as a weighted sum of a data-driven term and a physics constraint term:



$$ \mathcal{L} = \lambda_d \mathcal{L}_{\text{data}} + \lambda_p \mathcal{L}_{\text{physics}} + \lambda_r \mathcal{L}_{\text{reg}} $$


๐ŸŽ“

Here, $\mathcal{L}_{\text{data}}$ is the squared error with observed data, $\mathcal{L}_{\text{physics}}$ is the residual of the governing equation, and $\mathcal{L}_{\text{reg}}$ is the regularization term. Adjusting the weight parameters $\lambda$ greatly affects learning stability and accuracy.




Generalization Performance and Extrapolation Problem

๐Ÿง‘โ€๐ŸŽ“

Please tell me about "Generalization Performance and the Extrapolation Problem"!


๐ŸŽ“

The biggest challenge for surrogate models is prediction accuracy outside the range of training data (extrapolation region). Incorporating physical laws can improve extrapolation performance, but complete guarantees are difficult.




Curse of Dimensionality

๐Ÿง‘โ€๐ŸŽ“

Please tell me about the "Curse of Dimensionality"!


๐ŸŽ“

When the dimension of the input parameter space is high, the required number of samples increases exponentially. Efficient sample placement through Active Learning or Latin Hypercube Sampling (LHS) is super important.



$$ N_{\text{samples}} \propto d^{\alpha}, \quad \alpha \geq 1 $$

Assumptions and Applicability Limits

๐Ÿง‘โ€๐ŸŽ“

Is this formula not universal? When can't it be used?


๐ŸŽ“
  • The training data must sufficiently represent the physics of the analysis target.
  • The relationship between input parameters and output must be smooth (if there are discontinuities, domain partitioning is necessary).
  • Reducing computational cost is the main purpose; conventional solvers should be used in conjunction for final verification requiring high accuracy.
  • If the quality of training data (mesh-converged, V&V completed) is insufficient, model reliability decreases.

๐Ÿง‘โ€๐ŸŽ“

Ah, I see! So that's how the training data being the analysis target works.


Dimensionless Parameters and Dominant Scales

๐Ÿง‘โ€๐ŸŽ“

Teacher, please tell me about "Dimensionless Parameters and Dominant Scales"!


๐ŸŽ“

Understanding the dimensionless parameters governing the physical phenomenon being analyzed forms the basis for appropriate model selection and parameter setting.


๐ŸŽ“
  • Peclet Number Pe: Relative importance of convection and diffusion. Pe >> 1 indicates convection dominance (stabilization techniques required).
  • Reynolds Number Re: Ratio of inertial forces to viscous forces. Fundamental parameter for fluid problems.
  • Biot Number Bi: Ratio of internal conduction to surface convection. For Bi < 0.1, the lumped capacitance method is applicable.
  • Courant Number CFL: Indicator of numerical stability. For explicit methods, CFL โ‰ค 1 is required.

๐Ÿง‘โ€๐ŸŽ“

Ah, I see! So that's how the physics of the analysis target works.



Verification by Dimensional Analysis

๐Ÿง‘โ€๐ŸŽ“

Please tell me about "Verification by Dimensional Analysis"!


๐ŸŽ“

Dimensional analysis based on Buckingham's ฮ  theorem is effective for order-of-magnitude estimation of analysis results. Using characteristic length $L$, characteristic velocity $U$, and characteristic time $T = L/U$, the order of each physical quantity is estimated beforehand to confirm the validity of the analysis results.


๐Ÿง‘โ€๐ŸŽ“

I see. So if the physics of the analysis target is understood, it's basically okay to start?


Classification of Boundary Conditions and Mathematical Characteristics

๐Ÿง‘โ€๐ŸŽ“

I've heard that if you get the boundary conditions wrong, everything fails...


TypeMathematical ExpressionPhysical MeaningExample
Dirichlet Condition$u = u_0$ on $\Gamma_D$Specification of variable valueFixed wall, specified temperature
Neumann Condition$\partial u/\partial n = g$ on $\Gamma_N$Specification of gradient (flux)Heat flux, force
Robin Condition$\alpha u + \beta \partial u/\partial n = h$Linear combination of variable and gradientConvective heat transfer
Periodic Boundary Condition$u(x) = u(x+L)$Spatial periodicityUnit cell analysis
๐ŸŽ“

Choosing appropriate boundary conditions directly affects solution uniqueness and physical validity. Insufficient boundary conditions lead to an ill-posed problem, while excessive ones cause contradictions.



๐Ÿง‘โ€๐ŸŽ“

I've grasped the overall picture of Reduced Order Models (ROM)! I'll try to be mindful of it in my practical work from tomorrow.


๐ŸŽ“

Yeah, you're doing great! Actually trying things out is the best way to learn. If you have any questions, feel free to ask anytime.


Coffee Break Casual Talk

Mathematical Foundation of ROMโ€”POD (Fluid Version of Principal Component Analysis) and Galerkin Projection

The most classical method for Reduced Order Models (ROM) is the POD-Galerkin method. It performs Singular Value Decomposition (SVD) on numerous FEM snapshot solutions to extract "directions with large solution variation" (POD modes), and uses Galerkin projection to reduce the original equation system to a low-dimensional system of ordinary differential equations. A fluid problem with 1 million degrees of freedom can be transformed into an ODE with 10-50 modes, reducing computational cost by over 1000 times. However, Galerkin projection often has stability issues, and for nonlinear problems, the "closure problem" (the effect of omitted modes) causes accuracy degradation.

Computational Methods for Reduced Order Model (ROM)

๐ŸŽ“

Explains numerical methods and algorithms for implementing Reduced Order Models (ROM).



Discretization and Calculation Procedure

๐Ÿง‘โ€๐ŸŽ“

How do you actually solve this equation on a computer?


๐ŸŽ“

As data preprocessing, normalization/standardization of input features is important. Since CAE data have vastly different scales for each physical quantity, appropriate selection of Min-Max normalization or Z-score normalization is necessary. For learning algorithm selection, appropriate methods should be chosen according to data volume, dimensionality, and degree of nonlinearity.



Implementation Considerations

๐Ÿง‘โ€๐ŸŽ“

What is the most important thing to be careful about when using Reduced Order Models (ROM) in practical work?


๐ŸŽ“

Implementation using the Python ecosystem (scikit-learn, PyTorch, TensorFlow) is common. Key implementation aspects include learning acceleration via GPU parallelization, automatic hyperparameter tuning, and preventing overfitting through cross-validation. Using the HDF5 format is recommended for efficient I/O processing of large-scale CAE data.



Verification Methods

๐Ÿง‘โ€๐ŸŽ“

Teacher, please tell me about "Verification Methods"!


๐ŸŽ“

It's important to use k-fold cross-validation, Leave-One-Out method, and holdout method appropriately according to the purpose, and to evaluate prediction performance comprehensively using coefficient of determination Rยฒ, RMSE, MAE, and maximum error.


๐Ÿง‘โ€๐ŸŽ“

Now I understand what my senior meant when they said, "At least do cross-validation properly."


Code Quality and Reproducibility

๐Ÿง‘โ€๐ŸŽ“

What is the most important thing to be careful about when using Reduced Order Models (ROM) in practical work?


๐ŸŽ“

Implementation using the Python ecosystem (scikit-learn, PyTorch, TensorFlow) is common. Key implementation aspects include learning acceleration via GPU parallelization, automatic hyperparameter tuning, and preventing overfitting through cross-validation. Using the HDF5 format is recommended for efficient I/O processing of large-scale CAE data.


Related fields

Structural AnalysisFluid AnalysisV&V ยท Quality Assurance
Rate this article
Thank you for your feedback!
Helpful
More details
Report error
Helpful
0
More details
0
Report error
0
Written by NovaSolver Contributors
Anonymous Engineers & AI โ€” Sitemap
About the Authors