penalty-projection method for incompressible flows

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

 May 8, 2006, 03:49 penalty-projection method for incompressible flows #1 Isabel Guest   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, C.Lapuerta...in 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! Isa

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Maarten de Jong Main CFD Forum 7 April 30, 2012 02:23 Jinwon Main CFD Forum 1 December 4, 2007 22:13 TT Main CFD Forum 1 August 17, 2006 10:42 Woo Tae Kim Main CFD Forum 2 June 29, 2003 01:46 Bipin Lokhande Main CFD Forum 1 November 8, 1998 14:07

All times are GMT -4. The time now is 08:39.