CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Pre-Processing

Numerical problem in buoyantBoussinesqPimpleFoam

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

LinkBack Thread Tools Display Modes
Old   July 3, 2012, 08:43
Default Numerical problem in buoyantBoussinesqPimpleFoam
New Member
Vitor Campos
Join Date: Jun 2012
Posts: 1
Rep Power: 0
vtcampos is on a distinguished road

I am trying to simulate heat transfer in a small vertical cilindrical pipe. The pipe is established with its axis in the y direction, which is the direction of the acceleration g.
Air comes from the inlet which is the lower part of the pipe, flows through it and exits at the upper part. It looks similar to this:

| |
| | walls
| |
| |

I am using the buoyantBoussinesqPimpleFoam with kOmegaSST turbulence model.

The boundary conditions for time 0 are as following:
U: inlet fixedValue (0 8 0); outlet zeroGradient; wall fixedValue (0 0 0) (no slip condition)
T: inled fixedValue, 300; outlet zeroGradient; wall of the pipe fixedValue 450;
p: calculated uniform 0 in all the boundaries;
p_rgh: buoyantPressure uniform 0 in all the boundaries;
alphat: inlet fixedValue 0; outlet zeroGradient; wall alphatWallFunction uniform 0
kappat: inlet calculated uniform 0; outlet zeroGradient; wall kappatJayatillekeWallFunction value=uniform 0 and Prt 0 0.85;
nut: inlet calculated uniform 0; outlet zeroGradient; wall nutkWallFunction uniform 0;
k: inlet fixedValue uniform 0.3267; outlet zeroGradient: wall kqRWallFunction uniform 0.3267
omega: inlet fixedValue uniform 3760; outlet zeroGradient: wall omegaWallFunction uniform 3760

The problem appears when I try to run the solver and it asks for a reference point and value for the field p_rgh. I have set the reference value as 0, since my conditions are null pressure at outlet. Therefore I decided to set my reference point at some cell on the outlet boundary. However, no matter which cell or reference point I choose, in the outlet or not, it comes to a numerical error with unphysical velocity and pressure values (for example U = 500 m/s) at that exact point.
Can any of you help me finding my error on the simulation?
vtcampos is offline   Reply With Quote

Old   July 4, 2012, 05:12
Senior Member
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 254
Rep Power: 11
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Dear Vitor,

Try to put a

fixedValue uniform 0 for p_rgh on the outlet instead of buoyantpressure there, I think the solver then also does not need the reference cell or reference value, which it should only need in a closed volume.

tomf is offline   Reply With Quote

Old   July 7, 2013, 07:55
New Member
Korichi Abdelkader
Join Date: Jan 2013
Posts: 11
Rep Power: 4
AbdelkaderDZ is on a distinguished road
Dear Tom,
I have a problem to simulate the natural convection in channel, I can't found the appropriate boundary conditions for p and p_rgh at the inlet and outlat channel.
AbdelkaderDZ is offline   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
conduction problem venkataramana OpenFOAM 3 December 1, 2013 08:30
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43
Problem Importing Geometry ProE to CFX fatb0y CFX 3 January 14, 2012 20:42
Can I solve this problem by Fluent? Kai_kc FLUENT 1 October 27, 2010 05:29
Numerical solution to the rotating disk problem? johny Main CFD Forum 7 September 5, 2005 05:53

All times are GMT -4. The time now is 21:35.