Analysis of velocity-components decoupled projection method for the incompressible Navier–Stokes equations. (April 2016)

• Record Type: Journal Article
• Authors: Pan, Xiaomin; Lee, Changhoon; Kim, Kyoungyoun; Choi, Jung-Il
• Is Part Of: Computers & mathematics with applications. Volume 71:issue 8(2016)
• Page Start: 1722
• Page End: 1743
• Publication Date: 2016-04
• Subjects: Projection method -- Velocity-components decoupling -- Second-order temporal accuracy -- Energy estimation -- von Neumann analysis -- Numerical stability; Electronic data processing -- Periodicals; Mathematics -- Data processing -- Periodicals; 510.28541
• DOI: 10.1016/j.camwa.2016.03.002
• Languages: English
• ISSNs: 0898-1221

Abstract

Abstract: We study the temporal accuracy and stability of the velocity-components decoupled projection method (VDPM) for fully discrete incompressible Navier–Stokes equations. In particular, we investigate the effect of three formulations of the nonlinear convection term, which include the advective, skew-symmetric, and divergence forms, on the temporal accuracy and stability. Second-order temporal accuracy of the VDPM for both velocity and pressure is verified by establishing global error estimates in terms of a discrete l 2 -norm. Considering the energy evolution, we demonstrate that the VDPM is stable when the time step is less than or equal to a constant. Stability diagrams, which display the distributions of the maximum magnitude of the eigenvalues of the corresponding amplification matrices, are obtained using von Neumann analysis. These diagrams indicate that the advective form is more stable than the other formulations of the nonlinear convection term. Numerical tests are performed in order to support the mathematical findings involving temporal accuracy and stability, and the effects of the formulations of the nonlinear convection term are analyzed. Overall, our results indicate that the VDPM along with an advective discrete convection operator is almost unconditionally stable, second-order accurate in time, and computationally efficient because of the non-iterative solution procedure in solving the decoupled momentum equations.