CFD Online Discussion Forums - Turbulence modelling in chrmultiregionsimplefoam

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Turbulence modelling in chrmultiregionsimplefoam (https://www.cfd-online.com/Forums/openfoam-solving/184983-turbulence-modelling-chrmultiregionsimplefoam.html)

Turbulence modelling in chrmultiregionsimplefoam

HI All,

I have been trying to solve a CHT problem, with air, water and solid. It is a test case before solving the main problem. It is a simple box with two tubes. Box is of steel and two tubes contain water and gas respectively.

I have solved the case with Laminar case and results are reasonable. But as soon as I twitched on turbulence model with "turbulence true;" in the turbulence properties dict, the run crashes.

I have tried numerous combinations. The run goes well when "turbulence false;"

Please help !

Turbulence= true?

Upload your turbulence files and the log of the error you're getting.

What turbulence model are you using?

Code:

FoamFile

{

    version    2.0;

    class    dictionary;

    format    ascii;

    location    "constant";

    object    turbulenceProperties;

}

simulationType    RAS;

RAS

{

    printCoeffs    on;

    RASModel    kEpsilon;

    turbulence    true;

    kEpsilonCoeffs

    {

        C3    -0.33;

        Cmu    0.09;

        sigmak    1.0;

        C1    1.44;

        sigmaEps    1.3;

        C2    1.92;

    }

}

Code:

Time = 7





Solving for fluid region gas

DILUPBiCG:  Solving for Ux, Initial residual = 0.9998339, Final residual = 6.225594e-005, No Iterations 1

DILUPBiCG:  Solving for Uy, Initial residual = 0.9997901, Final residual = 9.356765e-005, No Iterations 1

DILUPBiCG:  Solving for Uz, Initial residual = 0.9999247, Final residual = 6.460889e-005, No Iterations 1

DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 4.503303e-005, No Iterations 1

Min/max T:300 650

DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 9.244747e-005, No Iterations 259

DICPCG:  Solving for p_rgh, Initial residual = 2.396295e-005, Final residual = 8.989188e-007, No Iterations 20

DICPCG:  Solving for p_rgh, Initial residual = 8.997803e-007, Final residual = 8.997803e-007, No Iterations 0

DICPCG:  Solving for p_rgh, Initial residual = 8.997803e-007, Final residual = 8.997803e-007, No Iterations 0

time step continuity errors : sum local = 1.405708e+047, global = 6.975396e+042, cumulative = 6.975396e+042

Min/max rho:0.5365993 1.16163

DILUPBiCGStab:  Solving for epsilon, Initial residual = 1, Final residual = 4.706257e-008, No Iterations 6

bounding epsilon, min: -3.629814e+084 max: 2.95134e+086 average: 6.291088e+081

DILUPBiCGStab:  Solving for k, Initial residual = 1, Final residual = 5.365471e-008, No Iterations 5

bounding k, min: -4.599895e+068 max: 3.286911e+090 average: 2.401271e+066



Solving for fluid region water

DILUPBiCG:  Solving for Ux, Initial residual = 0.9992455, Final residual = 8.898344e-005, No Iterations 5

DILUPBiCG:  Solving for Uy, Initial residual = 0.9928763, Final residual = 8.119978e-005, No Iterations 5

DILUPBiCG:  Solving for Uz, Initial residual = 0.9997006, Final residual = 7.40432e-005, No Iterations 5

DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = nan, No Iterations 1001

Min/max T:300 650

DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = nan, No Iterations 501

DICPCG:  Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 501





--> FOAM FATAL IO ERROR: 

wrong token type - expected Scalar, found on line 0 the word 'nan'



file: E:/water/data.solverPerformance.p_rgh at line 0.



--> FOAM Warning : 

    From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)

    in file lnInclude/Scalar.C at line 107.



FOAM exiting

I have attached the turbulence model details and error messages... It crashes after few iterations..

Are these warnings anything to do with this problem?
I am using windows version.

Code:

--> FOAM Warning :

    From function dlOpen(const fileName&, const bool)

    in file MSwindows.C at line 1313

\   dlopen error : "Error 126: The specified module could not be found.

"

--> FOAM Warning :

    From function dlLibraryTable::open(const fileName&, const bool)

    in file db/dynamicLibrary/dlLibraryTable/dlLibraryTable.C at line 99

    could not load "libturbulenceModelSchemes.dll"

From your log, it seems like there is something wrong with this file:

"E:/Speed/08_Work/SimTest/Test/Test/system/water/data.solverPerformance.p_rgh"

I'm not sure if this is a file from the set up (it seems like there is some divergence in the simulation) but if it is, please upload it so I can take a look. Or better yet, upload your case

Thanks for the reply. I do not see any such file. How can upload the case? Give me some idea.. It is 11 MB