CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Low-Re turbulence models (

grtabor September 1, 2009 11:34

Low-Re turbulence models
Hi. I'm trying to run simpleFoam using the low-Re turbulence models. I've run a preliminary calculation using the standard k-e 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.45357e-05, Final residual = 9.97272e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.0012608, Final residual = 2.48792e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.000676913, Final residual = 7.49901e-07, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.0957274, Final residual = 9.47486e-07, No Iterations 318
DICPCG: Solving for p, Initial residual = 0.00598807, Final residual = 9.84517e-07, No Iterations 207
time step continuity errors : sum local = 3.38296e-08, global = -5.26895e-10, cumulative = -5.26895e-10

request for volScalarField RASModel::G from objectRegistry region0 failed
available objects of type volScalarField are


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.


gocarts September 1, 2009 19:20

Replace *WallFunction in k and epsilon fields
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 high-Re to low-Re. Something to check.

ancsa February 15, 2010 10:52

Hi Richard,
by replacing you mean deleting it and simply adding a small value like in epsilon and k file e. g. :
value uniform 1e-10;

gocarts February 15, 2010 11:10

You have to provide meaningful epsilon and k values, for more see:
but otherwise yes.

If you are not sure about the difference between high-Re and low-Re turbulence models then I'd suggest sticking with high-Re, i.e., turbulence models that use wall functions.

Hope this helps.

ancsa February 15, 2010 11:25

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.

All times are GMT -4. The time now is 16:48.