
[Sponsors] 
September 1, 2009, 11:34 
LowRe turbulence models

#1 
Senior Member
Gavin Tabor
Join Date: Mar 2009
Posts: 181
Rep Power: 8 
Hi. I'm trying to run simpleFoam using the lowRe turbulence models. I've run a preliminary calculation using the standard ke model, then switched models. Each time on startup I get the following output;
Starting time loop Time = 1001 DILUPBiCG: Solving for Ux, Initial residual = 9.45357e05, Final residual = 9.97272e08, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.0012608, Final residual = 2.48792e06, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.000676913, Final residual = 7.49901e07, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.0957274, Final residual = 9.47486e07, No Iterations 318 DICPCG: Solving for p, Initial residual = 0.00598807, Final residual = 9.84517e07, No Iterations 207 time step continuity errors : sum local = 3.38296e08, global = 5.26895e10, cumulative = 5.26895e10 request for volScalarField RASModel::G from objectRegistry region0 failed available objects of type volScalarField are 14 ( pPrevIter (1(0.3*exp(sqr((sqr(k)(nu*epsilon)))))) k ((((Cmu*(1exp((Am*(((sqrt(k)*y)nu)+1e15)))))((1exp((Aepsilon*(((sqrt(k)*y)nu)+1e15))))+1e15))*sqr(k))(epsilon+epsilonSmall)) nu p (symm(grad(U))&&grad(U)) ((((Cmu*((1exp((Am*(((sqrt(k)*y)nu)+1e15))))((1exp((Aepsilon*(((sqrt(k)*y)nu)+1e15))))+1e15)))*sqr(k))epsilon)*(symm(grad(U))&&grad(U))) ((1exp((Am*(((sqrt(k)*y)nu)+1e15))))((1exp((Aepsilon*(((sqrt(k)*y)nu)+1e15))))+1e15)) (sqr(k)(nu*epsilon)) (((sqrt(k)*y)nu)+1e15) (((C2*(1(0.3*exp(sqr((sqr(k)(nu*epsilon)))))))*epsilon)k) epsilon y ) followed by a stack trace. What am I doing wrong? I've tried this with both the LaunderSharmaKE and the LienLeschzinerLowRe models with the same result. From an initial inspection of the code I can't see that G is being read in (in fact I think its calculated) and its not complaining about that. Any suggestions gratefully received. Gavin 

September 1, 2009, 19:20 
Replace *WallFunction in k and epsilon fields

#2 
Senior Member
Richard Smith
Join Date: Mar 2009
Location: Enfield, NH, USA
Posts: 138
Blog Entries: 4
Rep Power: 8 
I came across the same problem  but it was my fault.
I hadn't replaced the epsilonWallFunction and kqRWallFunction conditions in the epsilon and k fields (case time directory) when I switched turbulence models highRe to lowRe. Something to check.
__________________
Symscape, Computational Fluid Dynamics for all 

February 15, 2010, 10:52 

#3 
New Member
Aniko Rakai
Join Date: Oct 2009
Location: Budapest
Posts: 29
Rep Power: 7 
Hi Richard,
by replacing you mean deleting it and simply adding a small value like in epsilon and k file e. g. : SL6 { value uniform 1e10; } ? 

February 15, 2010, 11:10 

#4 
Senior Member
Richard Smith
Join Date: Mar 2009
Location: Enfield, NH, USA
Posts: 138
Blog Entries: 4
Rep Power: 8 
You have to provide meaningful epsilon and k values, for more see:
http://www.cfdonline.com/Wiki/Turbu...ary_conditions but otherwise yes. If you are not sure about the difference between highRe and lowRe turbulence models then I'd suggest sticking with highRe, i.e., turbulence models that use wall functions. Hope this helps.
__________________
Symscape, Computational Fluid Dynamics for all 

February 15, 2010, 11:25 

#5 
New Member
Aniko Rakai
Join Date: Oct 2009
Location: Budapest
Posts: 29
Rep Power: 7 
Thanks for fast reply. I will try to go on as you suggested, I know how to calculate the values.
I'm calculating my eps and k profile from measurement data, only I have no points near the wall, I was using wall functions hence, I just wanted to see the difference in the results if I change and I read in the forums to use small figures but not zero. I will read more about it. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Turbulence model for low Reynolds number flow?  Nokadu  Main CFD Forum  3  May 26, 2013 11:42 
Discussion: Reason of Turbulence!!  Wen Long  Main CFD Forum  3  May 15, 2009 09:52 
Turbulence Models for very low Reynolds Numbers  Sudhi  CFX  1  March 2, 2007 19:03 
Multicomponent fluid  Andrea  CFX  2  October 11, 2004 05:12 