A novel convergence enhancement method based on Online Dimension Reduction Optimization

11/29/2022
by   Wenbo Cao, et al.
0

Iterative steady-state solvers are widely used in computational fluid dynamics. Unfortunately, it is difficult to obtain steady-state solution for unstable problem caused by physical instability and numerical instability. Optimization is a better choice for solving unstable problem because steady-state solution is always the extreme point of optimization regardless of whether the problem is unstable or ill-conditioned, but it is difficult to solve partial differential equations (PDEs) due to too many optimization variables. In this study, we propose an Online Dimension Reduction Optimization (ODRO) method to enhance the convergence of the traditional iterative method to obtain the steady-state solution of unstable problem. This method performs proper orthogonal decomposition (POD) on the snapshots collected from a few iteration steps, optimizes PDE residual in the POD subspace to get a solution with lower residual, and then continues to iterate with the optimized solution as the initial value, repeating the above three steps until the residual converges. Several typical cases show that the proposed method can efficiently calculate the steady-state solution of unstable problem with both the high efficiency and robustness of the iterative method and the good convergence of the optimization method. In addition, this method is easy to implement in almost any iterative solver with minimal code modification.

READ FULL TEXT

page 6

page 7

page 10

research
08/02/2022

Neural Basis Functions for Accelerating Solutions to High Mach Euler Equations

We propose an approach to solving partial differential equations (PDEs) ...
research
02/11/2022

Formal verification of iterative convergence of numerical algorithms

Physical systems are usually modeled by differential equations, but solv...
research
09/08/2020

L1-based reduced over collocation and hyper reduction for steady state and time-dependent nonlinear equations

The task of repeatedly solving parametrized partial differential equatio...
research
11/11/2020

Continuation of Probability Density Functions using a Generalized Lyapunov Approach

Techniques from numerical bifurcation theory are very useful to study tr...
research
08/04/2021

Partial Reuse AMG Setup Cost Amortization Strategy for the Solution of Non-Steady State Problems

The partial reuse algebraic multigrid (AMG) setup cost amortization stra...
research
09/21/2017

Convergence characteristics of the generalized residual cutting method

The residual cutting (RC) method has been proposed for efficiently solvi...
research
06/27/2012

An Iterative Locally Linear Embedding Algorithm

Local Linear embedding (LLE) is a popular dimension reduction method. In...

Please sign up or login with your details

Forgot password? Click here to reset