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

Exploding h/rho in rhoSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 17, 2010, 11:41
Default Exploding h/rho in rhoSimpleFoam
  #1
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
Dear Solver-Experts,

I am running a case in rhoSimpleFoam and I get exploding values for temperature and especially density after few iterations, but no divergence. This even happens for turbulence turned off. I have to blame a few cells in my prism layer (low-Re, layer next to wall) for this.
I am trying to bound rho and h in some way, but without success so far. I tried:

Code:
interpolationSchemes
{
      interpolate(rho) limitedLimitedLinear phi 1.0 1.0 1.25;

...
}
to bound rho between 1.0 and 1.25 but this does not work at all. I get higer densities after 4 iterations... How can these values be bounded? Thanks.

Regards

Last edited by bastil; February 17, 2010 at 11:59.
bastil is offline   Reply With Quote

Old   February 18, 2010, 06:05
Default
  #2
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: http://olesenm.github.io/
Posts: 777
Rep Power: 18
olesen will become famous soon enough
Quote:
Originally Posted by bastil View Post
Dear Solver-Experts,

I am running a case in rhoSimpleFoam and I get exploding values for temperature and especially density after few iterations
Have you tried severely underrelaxing rho? For example, a value of 0.05 might get you over the rough spots.
olesen is offline   Reply With Quote

Old   February 18, 2010, 06:53
Default
  #3
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 471
Rep Power: 11
bastil is on a distinguished road
Quote:
Originally Posted by olesen View Post
Have you tried severely underrelaxing rho? For example, a value of 0.05 might get you over the rough spots.
Thanks Mark. Yes I do this. I still have problems in some critical areas. I am also wondering about BCs at walls. We use for walls:
fixedValue, 1e-20 for k and epsilon (for LaunderSharma-k-epsilon)
zeroGradient for temperature
I am not sure for temperature... The only similar tutorial case is NacaAirfoil from turbFoam. This uses zeroGradient for temperature but also for k and epsilon which is obviously wrong?
Any other ideas? Thanks.

Regards BastiL
bastil is offline   Reply With Quote

Reply

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
TRANSONIC FLOW in RHOSIMPLEFOAM dinonettis OpenFOAM Running, Solving & CFD 9 October 21, 2009 12:02
rhoSimpleFoam OF1.6 sheintz OpenFOAM 4 September 3, 2009 07:41
Monitor viscosity - RhoSimpleFoam ABouvier OpenFOAM Post-Processing 5 April 24, 2009 18:37
RhoSimpleFoam FoamX spv24 OpenFOAM Running, Solving & CFD 1 July 21, 2008 10:29
Stability startup problems with rhoSimpleFoam olesen OpenFOAM Running, Solving & CFD 1 July 18, 2006 08:09


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