CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

SimpleFoam boundary conditions

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 6, 2005, 07:01
Default I am trying to evaluate OpenFO
  #1
Senior Member
 
Håkan Nilsson
Join Date: Mar 2009
Location: Gothenburg, Sweden
Posts: 203
Rep Power: 18
hani is on a distinguished road
I am trying to evaluate OpenFOAM for a workshop in hydraulic water turbine drafttube computations (www.turbine-99.org). One of the cases in the workshop is to generate a steady solution of the flow in the drafttube, so I use simpleFoam. Niklas Nordin has helped me generate the appropriate inlet boundary conditions for velocity and turbulent quantities. The grid is highly non-orthogonal (specified at the workshop homepage), so I understand that I must have at least one non-orthogonal correctors in the SIMPLE algorithm. Is that correct?

Of course I have convergence problems, and I see at this discussion page that another user was recommended to use the inletOutlet boundary condition when there is a risk for recirculation at the outlet. In the drafttube case I am sure that I will get recirculation at the outlet. Right now I use zeroGradient for U, k and epsilon, and fixedValue, uniform 0 for pressure. Then I get the same kind of problem with divergence of the continuity as the other user had. What boundary condition is recommended to use, and how do I specify it for U, p, k, epsilon?

In our in-house code a Neumann condition on all variables works fine even if there is recirculation. Of course the inflow values are not correct, but it does not seem to influence the results much. We don't specify a constant pressure at the outlet, but we set a reference pressure anywhere in the domain.

Is there any other reason why the continuity diverges although the prime variables seem to converge?

Håkan.
hani is offline   Reply With Quote

Old   January 10, 2007, 03:08
Default Dear Håkan, I an new user.
  #2
New Member
 
Ram Poudel
Join Date: Mar 2009
Location: Flagstaff, arizona, USA
Posts: 9
Rep Power: 17
rcpoudel is on a distinguished road
Dear Håkan,

I an new user. Just wondering how did you people (including Niklas Nordin) generate the appropriate inlet boundary conditions for velocity and turbulent quantities.
Appriciate your detail answer/hints.

Thanks in advance,
Ram
rcpoudel is offline   Reply With Quote

Old   January 10, 2007, 03:44
Default Ram, I have previously written
  #3
Senior Member
 
Håkan Nilsson
Join Date: Mar 2009
Location: Gothenburg, Sweden
Posts: 203
Rep Power: 18
hani is on a distinguished road
Ram, I have previously written a pretty detailed instruction of how to generate inlet profiles in the forum. However, I now also post an example (warning, it is an ugly piece of code and should be improved):

setCombBC.tgz

This example is for a circular inlet facing the y-direction with axisymmetric flow in the negative y-direction. The inlet boundary condition is linearly interpolated from profiles that are determined from the boundaryProperties and
inletProperties.dat, which you should place in your constant directory. The example is a pre-processing code, i.e. it should be run as a separate application just to set the boundary condition.

Read through it and adapt it to your application, maybe you will learn something.

Håkan.
hani is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Own boundary condition modified simpleFoam erorr in parallel execution sponiar OpenFOAM Running, Solving & CFD 1 August 27, 2008 10:16
Define non uniform TKE boundary condition in simpleFOAM qtian OpenFOAM Running, Solving & CFD 9 July 30, 2007 15:54
SimpleFoam boundary conditions changed in OF 14 adorean OpenFOAM Running, Solving & CFD 5 June 22, 2007 08:50
Rotating boundary simplefoam gabriel OpenFOAM Running, Solving & CFD 4 June 4, 2007 04:41
Integral boundary conditions turbulent intensitylength boundary conditions olesen OpenFOAM Running, Solving & CFD 0 July 27, 2006 08:18


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