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

penalty-projection method for incompressible flows

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

LinkBack Thread Tools Display Modes
Old   May 8, 2006, 03:49
Default penalty-projection method for incompressible flows
Posts: n/a
Dear CFD Forum, I have some problems when I try to implement a scheme to solve 3D Navier Stokes equations in a domain which has a boundary as an inflow and another one as an open boundary. Indeed, it is a rectangular tube whose dimensions are (0.002x0.04x0.002)=(x dimension ,y dimension,z dimension). As an inflow (in y=0.04) I put a constant inflow through (u,v,w)=(0,5.3,0) because the flow goes down the y axis. On the walls I put a non-slip boundary conditions through (u,v,w)=(0,0,0) and, finally, on the output I put Neumann boundary condition that says (viscosity* grad v- grad p)=0.

The scheme I am using is the one proposed by M. Jobelin, their interesting article called "A finite element penalty-projection method for incompressible flows". I am using FEM, first order in time and Q2/Q1 as elements for velocity and pressure respectively.

I solve the two systems consecutively: first velocity, second pressure. Later, I correct velocity components and pressure. I must say I control the outflow in each iteration so it is the same as the inflow, but when convergence is reached, the difference between inflow and outflow is of 5.3%, which is a lot since I don't stop until a residual of 1e-12 is got in both, velocity and pressure. Besides, on the output boundary I obtain a nearly Poiseuille flow, but not a perfect one because the velocity component has a small slope towards "y axis" and "x axis".

I don't know where the problem is, because I think I am following all the steps pointed out in the article; it is a penalty projection method.

I would be really grateful if anyone could give some advice.

Thanks for your time and best regards!

  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
SIMPLE pressure correction in compressible flows Maarten de Jong Main CFD Forum 7 April 30, 2012 02:23
About flowfield-dependent variation(FDV) method? Jinwon Main CFD Forum 1 December 4, 2007 22:13
How is penalty method sound for unsteady flows? TT Main CFD Forum 1 August 17, 2006 10:42
PPE in Projection Method Woo Tae Kim Main CFD Forum 2 June 29, 2003 01:46
PDF Method for turbulent flows Bipin Lokhande Main CFD Forum 1 November 8, 1998 14:07

All times are GMT -4. The time now is 03:05.