CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

buoyantSimpleFoam: pressure instabilities

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 15, 2011, 06:43
Default buoyantSimpleFoam: pressure instabilities
  #1
New Member
 
Massimo Marro
Join Date: Mar 2011
Location: Lyon
Posts: 2
Rep Power: 0
massimo.marro is on a distinguished road
Hello,

I'm trying to simulate the velocity field for a heated cavity with buoyantSimpleFoam and a k-epsilon turbulence model. I'm using OpenFOAM-2.0.1.
I copied the characteristics of the solver from the directory "tutorial/heatTransfer/buoyantSimpleFoma/hotroom/system".

At the inlet a turbulent boundary layer is imposed, with fixed profiles of velocity, epsilon and k.
At the outlet zeroGradient for velocity, k, epsilon. On the walls, the boundary conditions are no slip condition for the velocity and compressible::wallFunction for k and epsilon. At the top fixed values of velocity, k and epsilon are imposed.
The boundary condition for the pressure is buoyantPressure everywhere.

If the inlet values of k and epsilon are low, pressure instabilities rise at the inlet and the simulation crashes.

Thanks for your help,
Massimo
massimo.marro is offline   Reply With Quote

Old   November 16, 2011, 08:53
Default
  #2
New Member
 
Alex
Join Date: Apr 2011
Location: München
Posts: 13
Rep Power: 6
al_pr is on a distinguished road
Hello Massimo,

can you post your boundary condition files for k, epsilon and pressure? I think this would help us to find an answer. Could you also post your fvSchemes and fvSolution file?

I found some threads of this topic in this forum. Many try to increase simply epsilon to make the simulation more stable.


Regards,
Alex
al_pr is offline   Reply With Quote

Old   November 16, 2011, 09:47
Default buoyantSimpleFoam
  #3
New Member
 
Massimo Marro
Join Date: Mar 2011
Location: Lyon
Posts: 2
Rep Power: 0
massimo.marro is on a distinguished road
Hello Alex,

thank you for your answer.
I have posted the boundary condition files for k, epsilon and p_rgh. The inlet condition for k is an experimental profile of a boundary layer.
The boundary condition for p are $internalField.

I have also posted the fvSchemes and fvSolution files.

If the inlet and top values of epsilon are "too" low, large pressure instabilities rise at the inlet.

Regards,
Massimo

epsilon.gz

k.gz

p_rgh.gz

fvSchemes.gz

fvSolution.gz
massimo.marro is offline   Reply With Quote

Reply

Tags
buoyantpressure, buoyantsimplefoam, heated cavity, pressure instabilities

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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Pressure instabilities with interDyMFoam for the floatingObject case nbadano OpenFOAM Running, Solving & CFD 14 March 8, 2011 09:19
Using pressure outlet BC Shreyansh FLUENT 0 November 5, 2010 15:57
Help ASAP! pressure inlet & outlet engahmed FLUENT 0 June 13, 2010 15:33
Setting pressure and velocity in inlet Asghari FLUENT 5 September 22, 2006 13:23


All times are GMT -4. The time now is 04:38.