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

CHANNEL FLOW: a question and a request

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

LinkBack Thread Tools Display Modes
Old   August 19, 2002, 10:53
Default CHANNEL FLOW: a question and a request
Posts: n/a

I have one "stupid" question and one demand:

Question number one:

I have created a code for the steady Navier-Stokes equations, and it seems to me that it gives me good results for the velocity, but for the pressure, he gives me the value zero everywhere. I left in the exit "free outflow" Boundary condition (Pressure gradient +pression) *normal = 0., and for the inlet, the usual parabolic profile. I solve the discretized problem with an iterative method (of type Schur -Complement). The question is: for this type of Boundary conditions, which type of pressure one has to find out? What type of "predicted" pressure is well adapted as initial value in the iterative process?


Question number two: Has anybody a FORTRAN77 file with an implementation of the GMRES solver method for a System Ax=B? Because I have one but I have a suspect of their performances

....Okay, Thank You very Much

  Reply With Quote

Old   August 19, 2002, 23:20
Default Re: CHANNEL FLOW: a question and a request
Posts: n/a
Why not give prescribed pressure at outlet?
  Reply With Quote

Old   August 19, 2002, 23:32
Default Re: CHANNEL FLOW: a question and a request
Posts: n/a
Hi, Carlos,

1) If the velocity is correct, check your pressure to

make sure if it is zero or not.

2) A simple boundary condition (gradient of pressure = 0 at the outlet, prescribed velocity profile at inlet) should give you good result.

3) Use a few simple linear systems to test you GMRES module. You can also get the subroutine from Netlib or elsewhere. If it works for such system, it should also works for your NSe.


  Reply With Quote

Old   August 20, 2002, 15:22
Default Re: CHANNEL FLOW: a question and a request
Posts: n/a
Thank you by your answers. Indeed, for the canal plan, I impose parabolic profile for the entry, zero for walls, but NOTHING for the pressure, and as result, I obtain parabolic profile everywhere, but for the pressure, a " zero " everywhere. One said to me that the fact of putting as condition ( Deformation+pression)* normal=0 gives exactly the idea of an outside pressure (and eventually zero), but what me surprise it is the fact of obtaining pressure zero everywhere, being given that in my calculations for the speed I use the speed as "slave" later to resolve equations for the pressure " Bt*inverse (S) B*p=F...


  Reply With Quote

Old   August 23, 2002, 05:55
Default Re: CHANNEL FLOW: a question and a request
Andrea Panizza
Posts: n/a
Hello carlos, fortran77 code GMRES can be found in the libraries dftemplates.tgz (double precision) and sftemplates.tgz (single) at The files are gmres.f and gmresrevcom.f. For what it concerns your code, you are solving incompressible steady NVs in a plane channel with constant area duct, aren't you? Then pressure has a constant negative pressure gradient, and cannot be uniformly zero, but decreases linearly from inflow to outflow. Actually, this is strictly true for infinitely long ducts or equivalently ducts in which periodic b.c. are imposed on velocity only <=> Poiseuille flow, but it should hold in your case, provided that H/L <<1 where L is the plane channel length and H its width. So I am surprised than you can get the right velocity, even though the pressure is wrong: you should have a velocity field whose divergence is non zero. Therefore, I would suggest you to check that the assigned boundary velocities are correct, and also that, _at_convergence_ , mass is conserved both globally and locally down to your prescribed tolerance (e.g., check differences between mass flow at inlet and mass flow at outlet, and checkl divergence in any cell). I am also very interested in your choice of a Schur Complement method for the treatment of pressure: could you please provide further details and references on your implementation? I guess that you start with a trial pressure field and 1) Solve momentum equations for velocity, considering p as a right hand side term, by GMRES and not CG, since, due to convection terms, the coefficient matrix is not symmetrical. 2) With the new velocity, you calculate the divergence of the velocity field and add it to the old pressure, multiplied by some factor beta. 3) With the new pressure, you go back to point 1)and iterate to convergence. Am I correct?

Goodbye, Andrea
  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
Open Channel Flow ElanMorin FLUENT 4 February 25, 2015 17:26
channel flow question keith Main CFD Forum 5 September 14, 2012 05:14
Flow Split Setup Question ed FLUENT 1 February 16, 2007 17:55
A question about the Modelling of a Channel selen CD-adapco 1 August 24, 2003 12:06
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 21:31

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