A split-step finite-element method for incompressible Navier-Stokes equations with high-order accuracy up-to the boundary

Abstract

An efficient and accurate finite-element algorithm is described for the numerical solution of the incompressible Navier-Stokes (INS) equations. The new algorithm that solves the INS equations in a velocity-pressure reformulation is based on a split-step scheme in conjunction with the standard finite-element method. The split-step scheme employed for the temporal discretization of our algorithm completely separates the pressure updates from the solution of velocity variables. When the pressure equation is formed explicitly, the algorithm avoids solving a saddle-point problem; therefore, our algorithm has more flexibility in choosing finite-element spaces. For efficiency and robustness, Lagrange finite elements of equal order for both velocity and pressure are used. Motivated by a post-processing technique that calculates derivatives of a finite element solution with super-convergent error estimates, an alternative numerical boundary condition is proposed for the pressure equation at the discrete level. The new numerical pressure boundary condition that can be regarded as a better implementation of the compatibility boundary condition improves the boundary-layer errors of the pressure solution. Normal-mode analysis is performed using a simplified model problem on a uniform mesh to demonstrate the numerical properties of our methods. Convergence study using ${\mathbb{P}}_1$ elements confirms the analytical results and demonstrates that our algorithm with the new numerical boundary condition achieves the optimal second-order accuracy for both velocity and pressure up-to the boundary. Benchmark problems are also computed and carefully compared with existing studies. Finally, as an example to illustrate that our approach can be easily adapted for higher-order finite elements, we solve the classical flow-past-a-cylinder problem using ${\mathbb{P}}_n$ finite elements with $n\ge 1$.

References in corpus (6)

• A partitioned scheme for fluid-composite structure interaction problems
• A stable partitioned FSI algorithm for incompressible flow and deforming beams
• A stable partitioned FSI algorithm for rigid bodies and incompressible flow in three dimensions
• A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis
• A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part II: General formulation
• A stable added-mass partitioned (AMP) algorithm for elastic solids and incompressible flow

Cited by in corpus (5)

• A POD-Galerkin reduced order model for a LES filtering approach
• Efficient split-step schemes for fluid-structure interaction involving incompressible generalised Newtonian flows
• High-order Methods for a Pressure Poisson Equation Reformulation of the Navier-Stokes Equations with Electric Boundary Conditions
• Improving stability of moving particle semi-implicit method by source terms based on time-scale correction of particle-level impulses
• Border mapping multi-resolution (BMMR) technique for incompressible projection-based particle methods