Dirichlet and Neumann boundary conditions for the pressure poisson equation of incompressible flow

Abstract

SUMMARYIn a recent paper Gresho and Sani' showed that Dirichlet and Neumann boundary conditions for the pressure Poisson equation give the same solution. The purpose of this paper is to confirm this (for one case at least) by numerically solving the pressure equation with Dirichlet and Neumann boundary conditions for the inviscid stagnation point flow problem. The Dirichlet boundary condition is obtained by integrating the tangential component of the momentum equation along the boundary. The Neumann boundary condition is obtained by applying the normal component of the momentum equation at the boundary. In this work solutions for the Neumann problem exist only if a compatibility condition is satisfied. A consistent finite difference procedure which satisfies this condition on non-staggered grids is used for the solution of the pressure equation with Neumann conditions. Two test cases are computed. In the first case the velocity field is given from the analytical solution and the pressure is recovered from the solution of the associated Poisson equation. The computed results are identical for both Dirichlet and Neumann boundary conditions. However, the Dirichlet problem converges faster than the Neumann case. In the second test case the velocity field is computed from the momentum equations, which are solved iteratively with the pressure Poisson equation. In this case the Neumann problem converges faster than the Dirichlet problem.

Leray weak solutions of the incompressible Navier Stokes system on exterior domains via the artificial compressibility method

“…The issue of the correct boundary condition to assign to the pressure was also studied by many other authors. Among them Gresho and Sani [22], (see also [39], [1]) showed that Dirichlet and Neumann boundary conditions for the pressure give the same solution. In fact they recover appropriate Dirichlet boundary conditions for the pressure by taking into account the normal derivative of the pressure associated with the Neumann boundary condition and by using Green's functions.…”

Computational relaxation method for modeling electrostatic systems with non-trivial geometries

“…We demonstrated with examples the core principles of this method, and how it can be applied to specific problems in physics. Furthermore, this model can be optimized, and adapted for a wide range of systems (see [1] and [2]). This, coupled with the easy access to the needed resources, makes the relaxation method and its applications an excellent topic to introduce numerical methods for solving physical problems.…”

“…which means that the average value of the solution over a sphere is equal to the function evaluated at its center 1 .…”

“…The relaxation method is used in many sub-disciplines of physics, from fluid mechanics [1], to astrophysics [2]. Its versatility and relative simplicity make it a perfect introduction to computational methods, and its application to a familiar problem in electrostatics is sure to be illuminating for an eager learner.…”

Direct primitive variable incompressible solutions for two and threedimensional viscous flows

“…For method (i), all the boundary conditions are satisfied through the groupings described in section (3); no additional numerical boundary conditions are required. In method (ii), the boundary conditions are satisfied through similar groupings, as portrayed in Fig.3.…”