A compact finite differences exact projection method for the Navier–Stokes equations on a staggered grid with fourth-order spatial precision. (2nd September 2015)

• Record Type: Journal Article
• Authors: Reis, G.A.; Tasso, I.V.M.; Souza, L.F.; Cuminato, J.A.
• Is Part Of: Computers & fluids. Volume 118(2015)
• Page Start: 19
• Page End: 31
• Publication Date: 2015-09-02
• Subjects: Navier–Stokes equations -- Compact finite differences -- Exact projection -- High-order methods; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2015.06.015
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: The method uses primitive variables (velocity and pressure) on a staggered grid. Temporal discretization is based on a second-order Adams–Bashforth scheme. Spatial discretizations are based solely on fourth-order compact finite differences. Discrete incompressibility constraint is imposed exactly, up to machine precision. Formal accuracy is verified for both velocity and pressure. Abstract: An exact projection method for the numerical solution of the incompressible Navier–Stokes equations is devised. In all spatial discretizations, fourth-order compact finite differences are used, including domain boundaries and the Poisson equation that arises from the projection method. The integration in time is carried out by a second-order Adams–Bashforth scheme. The discrete incompressibility constraint is imposed exactly (up to machine precision) by a simple and efficient discretization of the Poisson equation. Spatial and temporal accuracies, for both velocity and pressure, are verified through the use of analytical and manufactured solutions. The results show that the method converges with fourth-order accuracy in space and second-order accuracy in time, for both velocity and pressure. Additionally, two popular benchmark problems, the flow over a backward facing step and the lid-driven cavity flow, are used to demonstrate the robustness and correctness of the code.