Optimization of the nonlinear convergence rate of the SIMPLER Picard algorithm for steady incompressible internal flows through enforcement of continuity projections
We discuss the performance optimization of the Semi-Implicit Method for Pressure-Linked Equations-Revised (SIMPLER) Picard algorithm for steady incompressible internal flows. We discuss the nonlinear convergence of the Picard iteration as a function of the pressure and scalar potential continuity projections stemming from the SIMPLER algorithm, for three example problems. In particular, we discuss the choice of under-relaxation method, and choice of under-relaxation factors; the choice of the projection algorithm; and the required tolerance for the linear solve of the generalized Poisson equations for the pressure and scalar potential equations that arise from the projection operations. We conclude that the convergence of the nonlinear Picard iteration can be effectively controlled by the optimal enforcement of the continuity projections. Copyright (c) 2005 John Wiley & Sons, Ltd.
convergence rate simpler picard iterations projection methods navier-stokes equations solvers