Bridging Proper Orthogonal Decomposition methods and augmented Newton-Krylov algorithms: an adaptive model order
P Kerfriden1, P Gosselet2, S Adhikari3
1Cardiff University, Queen's Buildings, The Parade, Cardiff CF24 3AA, Wales, UK.
Summary
This study introduces a novel algorithm that enhances reduced-order modeling for complex, nonlinear problems with topological changes. The corrected hyperreduction method improves model accuracy efficiently, even when dealing with damage initiation.
Area of Science:
- Computational mechanics
- Numerical analysis
- Model order reduction
Background:
- Model order reduction (MOR) techniques, such as Proper Orthogonal Decomposition (POD), are essential for efficiently simulating complex systems.
- Classical Newton/Krylov solvers are widely used for solving nonlinear problems but can be computationally expensive for large-scale simulations.
- Simulating highly nonlinear problems with strong topological changes, like damage initiation, poses significant challenges for traditional MOR and solvers.
Purpose of the Study:
- To develop an efficient algorithm that bridges POD-based model order reduction and Newton/Krylov solvers.
- To enable real-time correction of reduced-order models for highly nonlinear problems with significant topological changes.
- To improve the accuracy and relevance of reduced-order models in scenarios involving damage initiation.
Main Methods:
- A novel algorithm is derived by connecting POD-based model order reduction with Newton/Krylov solvers.
- The algorithm implements a corrected hyperreduction method to address problems with strong topological changes.
- The approach focuses on 'on-the-fly' correction of the reduced-order model during simulation.
Main Results:
- The proposed algorithm effectively corrects reduced-order models for highly nonlinear problems.
- Significant improvements in the relevancy of the reduced-order model are achieved.
- The method demonstrates efficacy even when dealing with strong topological changes, such as those in damage initiation problems.
- The accuracy enhancements come with reasonable additional computational costs.
Conclusions:
- The developed bridge between POD-based MOR and Newton/Krylov solvers offers an efficient solution for complex nonlinear simulations.
- The corrected hyperreduction method significantly enhances the performance of reduced-order models in the presence of topological changes.
- This approach provides a valuable tool for simulating phenomena like damage initiation with improved accuracy and efficiency.
Related Concept Videos
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
393
Mechanistic models play a crucial role in algorithms for numerical problem-solving, particularly in nonlinear mixed effects modeling (NMEM). These models aim to minimize specific objective functions by evaluating various parameter estimates, leading to the development of systematic algorithms. In some cases, linearization techniques approximate the model using linear equations.
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
393
Linear Approximation in Time Domain
394
Nonlinear systems often require sophisticated approaches for accurate modeling and analysis, with state-space representation being particularly effective. This method is especially useful for systems where variables and parameters vary with time or operating conditions, such as in a simple pendulum or a translational mechanical system with nonlinear springs.
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
394
Newton’s Method
121
Newton’s Method is a powerful iterative technique for approximating the roots of real-valued, differentiable functions, particularly when analytical solutions are impractical. This approach is widely used in scientific computing, engineering, and finance, where equations may be too complex for traditional algebraic methods to handle. The method relies on an iterative process that refines an initial estimate using the function’s derivative to approach the true solution progressively.
121
Mechanical Systems
830
Mechanical systems are analogous to to electrical networks where springs and masses play similar roles to inductors and capacitors, respectively. A viscous damper in mechanical systems functions similarly to a resistor in electrical networks, dissipating energy. The forces acting on a mass in such systems include an applied force in the direction of motion, counteracted by forces from the spring, a viscous damper, and the mass's acceleration. This interplay of forces is mathematically...
830
Methods of Medium Optimization
15
Optimizing growth media enhances microbial proliferation and maximizes product yield. Statistical experimental design methodologies provide structured and reproducible approaches, offering progressively higher levels of robustness and efficiency.The One-Factor-at-a-Time (OFAT) MethodThe One-Factor-at-a-Time (OFAT) method involves adjusting a single variable while keeping all others constant. However, it cannot detect interactions between variables, often leading to suboptimal outcomes when...
15
Unsymmetric Loading of Thin-Walled Members: Problem Solving
582
The shear center of a channel section with uniform thickness, height, and width, is determined by computing the shear force in the member and calculating the moments of inertia of the sections.
To compute the shear forces, find the shear flow at a specific distance from the endpoint using the vertical shear and the moment of inertia values. The total shear force on the flange is calculated by integrating the shear flow from one end of the flange to the other.
Next, calculate the moments of...
To compute the shear forces, find the shear flow at a specific distance from the endpoint using the vertical shear and the moment of inertia values. The total shear force on the flange is calculated by integrating the shear flow from one end of the flange to the other.
Next, calculate the moments of...
582


