CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Pressure Poisson Equation

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   February 23, 2003, 18:00
Default Pressure Poisson Equation
Posts: n/a
I am trying to implement a finite element for the incompressible navier-stokes and i have trouble with the Pressure Poisson equation. the discrete laplace operator is as follows CT*ML-1*C , where C: gradient matrix and ML is the lumped mass, CT:divergence matrix the problem is, the boundary conditions coming from the velocity. I can form this matrix in two forms. first one: reducing the rows and columns of the known velocity boundary condition from the C and CT matrices and second without reducing these entries. In two cases the dimension of my laplace operator will be the same but they won't be the same. Which one should I implement?
  Reply With Quote

Old   March 4, 2003, 04:14
Default Re: Pressure Poisson Equation
Posts: n/a
I am not familar with projection FEM.

but according to my recent re-look at projection FDM, in the case of with pressure poisson equation, the projection procedure could easily cause the FD to lose one of its good properties: strong form.

however, this could be not a problem at all for FEM, since classical FEM is already in weak form and has nothing to lose.

Acusim, seeming the first FEM software using Galerkin/Least Square in a projection procedure, may have something for you. please check

and there are a big bulk of papers on projection FEM.

and if you still have the smae question in summer, you can email me at, because i'll implement the projection version of my current mixed-formulation sometime in May.
  Reply With Quote

Old   March 4, 2003, 14:41
Default Re: Pressure Poisson Equation
Posts: n/a
The projection / fractional step method can be interpretated with the generalized block LU decomposition (see JB Perot, JCPH (1993), 108:51-58).

In fact, the pressure poisson equation is obtained using the LU decomposition. In this sense, it should contain the velocity BC from the momentum eqs. You may choose not to as in your 2nd case. The difference will show in a numerical boundary layer - reduced accuracy for pressure in time. You may find tons of paper out there for this issue.

Read the mentioned reference first.

Hope this helps.

  Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Implicit elliptic pressure field equation. TheBoyce Main CFD Forum 2 June 3, 2011 13:08
Strcture of pressure correction equation in PISO stable over time ? sebi Main CFD Forum 1 October 26, 2010 22:19
Two-Phase Buoyant Flow Issue Miguel Baritto CFX 4 August 31, 2006 12:02
newbie: simple questions about pressure poisson eq cfd-newbie Main CFD Forum 19 February 4, 2006 07:48
FEM pressure poisson equation (Implicit) cfd-beginner Main CFD Forum 0 August 9, 2005 13:32

All times are GMT -4. The time now is 07:41.