Artificial pressure for pressure-linked equation

Abstract The discretized pressure-linked equation can be proven to possess a singular coefficient matrix. This implies that the pressure solution corresponding to a given velocity would not exist unless this velocity precisely satisfies the continuity equation. In the present investigation, an artificial source term is added to the pressure-linked equation. This additional term will generate an extra pressure called the ‘artificial pressure’ for each updated velocity to compensate for their nonzero dilation. The use of the artificial pressure is equivalent to creating an extra mass for the need of the updated velocity in satisfying the continuity equation. This treatment guarantees the existence and uniqueness of the pressure solution such that the pressure can be directly solved without recourse to the conventional pressure correction equation. Based on the concept of artificial pressure, the APPLE (Artificial Pressure for Pressure-Linked Equation) and the NAPPLE (Nonstaggered APPLE) algorithms are developed for incompressible flows. Through two well-known examples, the APPLE algorithm is found to produce essentially the same results with only about 40% CPU time as compared with the SIMPLER algorithm. Due to its simplicity, the NAPPLE algorithm has a potential use for problems with arbitrarily shaped domain as long as the grid size is not large.

[1]  S. Biringen,et al.  On pressure boundary conditions for the incompressible Navier-Stokes equations using nonstaggered grids , 1988 .

[2]  F. W. Schmidt,et al.  USE OF A PRESSURE-WEIGHTED INTERPOLATION METHOD FOR THE SOLUTION OF THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON A NONSTAGGERED GRID SYSTEM , 1988 .

[3]  S. Abdallah,et al.  Numerical solutions for the incompressible Navier-Stokes equations in primitive variables using a non-staggered grid, 11 , 1987 .

[4]  W. L. Hankey,et al.  Use of Primitive Variables in the Solution of Incompressible Navier-Stokes Equations , 1979 .

[5]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[6]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[7]  Suhas V. Patankar,et al.  A Calculation Procedure for Two-Dimensional Elliptic Situations , 1981 .

[8]  Suhas V. Patankar,et al.  Recent Developments in Computational Heat Transfer , 1988 .

[9]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[10]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[11]  L. Shong-Leih.  Weighting function scheme and its application on multidimensional conservation equations , 1989 .

[12]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

[13]  Andrew Pollard,et al.  COMPARISON OF PRESSURE-VELOCITY COUPLING SOLUTION ALGORITHMS , 1985 .

[14]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[15]  S. Majumdar.  Role of underrelaxation in momentum interpolation for calculation of flow with nonstaggered grids , 1988 .

[16]  M. Peric.  A finite volume method for the prediction of three-dimensional fluid flow in complex ducts , 1985 .

[17]  Shong-Leih Lee.  A STRONGLY IMPLICIT SOLVER FOR TWO-DIMENSIONAL ELLIPTIC DIFFERENTIAL EQUATIONS , 1990 .

[18]  S. Abdallah.  Numerical solutions for the pressure poisson equation with Neumann boundary conditions using a non-staggered grid, 1 , 1987 .

[19]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.