Pressure-Poisson formulation and projection methods
by Yi Zhang
Take divergence of the momentum equation, one will finally have
This is the pressure-Poisson formulation, taking divergence-free condition of velocity into consideration. This equation presents the explicit form of pressure, instead of implicit in “standard” momentum formulations. Pressure-Poisson equation is widely used in projection methods, in which the velocity is decomposed into solenoidal part and irroational part by Hodge or Helmholtz decomposition. Chorin & Marsden’s this book contains extensive introduction for the latter, together with many other insightful comments.
Basically, in projection the velocity if splitted as
Or one can say that the velocity space is the direct product of two subspaces. The solenoidal nature of a flow field is projected to one of the subspaces: . This projection, denoted by operator , function as
The stability of projection operator follows the fact that its norm is less or equal than one, which can be intuitively seen from the action of “projection”.
Generally, in a typical predictor-corrector method using pressure to correct the velocity pseudo-step, the velocity in predictor sub-step is obtained using information from time step , i.e. , and (explicit methods) or together with information at time step (implicit methods). Since does not satisfy divergence-free condition, divergence operator is applied to the equation connecting and , in order to obtain the pressue required to satisfy solenoidal and , and this is the pressue-corrector sub-step. Usually, the second sub-step would give a Poisson equation of pressue, as the example by Idelsohn et al I cited in this post.