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.
翻译:在计算流体动态中广泛使用稳定状态稳定解答器。 不幸的是, 很难找到稳定状态的方法来解决由物理不稳定和数字不稳定造成的不稳定问题。 优化是解决不稳定问题的更好选择, 因为稳定状态解决方案总是最优化的极点, 不论问题不稳定还是条件不当, 但是由于优化变量太多, 很难解决部分差异方程式( PDEs ) 。 在此研究中, 我们提议了一种在线降低维度优化化( ODRO) 方法, 以加强传统迭代方法的趋同, 以获得不稳定问题的稳定状态解决方案。 这种方法在从几个迭代步骤中采集的近似镜片上进行适当的 orthogoal 解析( POD), 优化 PDE 在 POD 子空间的剩余部分以较低残余值获得解决方案, 然后继续使用优化的解决方案作为初始值, 重复以上三个步骤, 直至剩余部分趋同。 几个典型的案例表明, 拟议的方法可以有效地计算出不稳定问题的稳定状态解决方案, 以高效和稳健的稳健性方法 。 在迭代法中, 最易采用的方法是最佳的方法。