An improved numerical scheme for the SIMPLER method on nonorthogonal curvilinear coordinates: SIMPLERM

In this article, an improved numerical algorithm named SIMPLERM is proposed for incompressible fluid flow computations on the nonstaggered and nonorthogonal curvilinear grid system. In the proposed algorithm, the contravariant velocities are chosen as the cell face velocities and the Cartesian components as the primary variables. The velocity under-relaxation factor is incorporated into the momentum interpolation, and special treatment is adopted to avoid the underrelaxation factor dependence of the velocity solution. In addition, a 1-delta pressure difference is introduced into the interfacial contravariant velocity determination. Compared with the existing implementation methods of the SIMPLE family on non-staggered and nonorthogonal grids, the SIMPLERM algorithm can guarantee the coupling between velocity and pressure, underrelaxation independence of the solution, and satisfaction of the conservation law, while still possessing sufficient robustness.