[ebtedaei]

Dear All, I want to use Reynolds Stress Model (for example LRR) at OpenFOAM but I do not know how to use and define at the software?

Please help me.

[ebtedaei]

Where is someone to help me???
Thanks in advance

[Ardali]

Please give more details,
LRR needs epsilon and R file to start,
Ardalan

[ebtedaei]

Hi, Thanks for your reply,
My case is a CYCLONE. I want to use Reynolds Stress Model (LRR) at OpenFOAM but I do not know how to use and define at the software.
To begin, I found a file R and changed it. My boundary conditions are inlet, outlet and walls.

File U is as fallow:
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type noSlip;
}
inlet
{
type pressureInletOutletVelocity;
value uniform (-10 0 0);
}
outlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
}

File R is as fallows:
dimensions [0 2 -2 0 0 0 0];
internalField uniform (0 0 0 0 0 0);
boundaryField
{
wallS
{
type kqRWallFunction;
value uniform (0 0 0 0 0 0);
}

inlet
{
type pressureInletOutletVelocity;
value uniform (0.0015325 0 0 0.0015325 0 0.00306504);
}

outlet
{
type zeroGradient;
}
}

According to the file U, What changes do I need on R?
Have you a file R to send for me?
I need the fvSolution and fvSchemes files for model LRR.

Best,
Ali

[ebtedaei]

--> FOAM FATAL IO ERROR:
Cannot find patchField entry for walls

file: /home/dadehnegar/OpenFOAM/dadehnegar-4.0/run/project/MYcyclone/air+water/rsmLRR/0/R.boundaryField from line 25 to line 37.

According to the files U and R in previous comment, help me to solve this error please.
Thanks
Ali

[Ardali]

Please check the spelling for walls in R file.
and let me know?

[ebtedaei]

Hi, the error was solved. Except the file R, What other changes are needed for model LRR?
Now I think that the fvSolution file and the fvSchemes must be changed based on file R. Do you agree? Please send me a example of these files for model LRR.





Thanks,

Ali

[Ardali]

Report the error, then I will help you.
It seems you have big y+ close to the wall.
I recommend to use y+~1 and start with LaunderSharmaKe.
fvSchmes and fvSolution for LRR is not something special. Add an alforithem for R in fvSolution and appropriate shcmes in fvScheme.
Ardalan

[ebtedaei]

Originally Posted by Ardali

Report the error, then I will help you.
It seems you have big y+ close to the wall.
I recommend to use y+~1 and start with LaunderSharmaKe.
fvSchmes and fvSolution for LRR is not something special. Add an alforithem for R in fvSolution and appropriate shcmes in fvScheme.
Ardalan

Thanks for your reply.

--> FOAM FATAL IO ERROR:
[4] attempt to read beyond EOF
[4]
[4] file: IOstream.divSchemes.div(phi,R) at line 0.
[4]
[4] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[4] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[4]
FOAM parallel run exiting
[4]
[3] file: IOstream.divSchemes.div(phi,R)[7]
at line 0.
[3]
[3] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[3] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[3]
FOAM parallel run exiting

-----------------------------------------------------------------------------------------------------------------
fvSchemes:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(phi,R) bounded Gauss limitedLinear;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}


// ************************************************** *********************** //

[ebtedaei]

Originally Posted by ebtedaei

Thanks for your reply.

--> FOAM FATAL IO ERROR:
[4] attempt to read beyond EOF
[4]
[4] file: IOstream.divSchemes.div(phi,R) at line 0.
[4]
[4] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[4] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[4]
FOAM parallel run exiting
[4]
[3] file: IOstream.divSchemes.div(phi,R)[7]
at line 0.
[3]
[3] From function virtual Foam::Istream& Foam::ITstream::read(Foam::token&)
[3] in file db/IOstreams/Tstreams/ITstream.C at line 82.
[3]
FOAM parallel run exiting

-----------------------------------------------------------------------------------------------------------------
fvSchemes:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(phi,R) bounded Gauss limitedLinear;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}


// ************************************************** *********************** //

What changes I have to do in fvSchemes and fvSulotion files?

[Ardali]

Change all schemes to upwind.
Put all your boundary conditions to something same. You had noSlip. Put fixedValue.
Try it!
Ardalan

[ebtedaei]

Hi, in the first stage, fvSchemes file is as fallow:

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

the error file is attached for above file.

[Ardali]

you have not defined div(phi,R).
Simply define:
div(phi,R) Gauss upwind;

Most likely you will receive similar errors for other schemes.
Add them and put upwind.
Ardalan

[ebtedaei]

Thanks. I solved the errors at the fvSchemes file. Now I got a another error at the fvSolution file is as fallow:

[0] --> FOAM FATAL IO ERROR:
[0] keyword RFinal is undefined in dictionary "/home/dadehnegar/OpenFOAM/dadehnegar-4.0/run/project/MYcyclone/air+water/rsmLRR/system/fvSolution.solvers"
[0]
[0] file: /home/dadehnegar/OpenFOAM/dadehnegar-4.0/run/project/MYcyclone/air+water/rsmLRR/system/fvSolution.solvers from line 22 to line 63.


fvSolution:
solvers
{
"alpha.water.*"
{
nAlphaCorr 2;
nAlphaSubCycles 1;
cAlpha 1;

MULESCorr yes;
nLimiterIter 3;

solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-8;
relTol 0;
}

pcorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-5;
relTol 0;
}

p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.05;
}

p_rghFinal
{
$p_rgh;
relTol 0;
}

"(U|k|epsilon).*"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-06;
relTol 0;
minIter 1;
}
}

PIMPLE
{
momentumPredictor no;
nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
equations
{
".*" 1;
}
}


Please help me.
Best,
Ali

[Ardali]

Please read the tutorials, these are basics in OpenFoam.
You use pimple algorithem and it need RFinal in fvSolution.
Introduce it.
Have fun!
Ardalan

[ebtedaei]

Hi, How do I add this option (RFinal)? Where? Can you introduce a fvSolution file that it has a Simple algorithem?

Thanks in advance,
Ali

[Ardali]

put this in you fvSolution

"(R|RFinal)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-9;
relTol 0.01;
nSweeps 1;
};

see if it works!
Ardalan

[ebtedaei]

Originally Posted by Ardali

put this in you fvSolution

"(R|RFinal)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-9;
relTol 0.01;
nSweeps 1;
};

see if it works!
Ardalan

I used it but solution was divergent!

[ebtedaei]

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.000120212, Final residual = 2.15417e-10, No Iterations 2
Phase-1 volume fraction = 0.961255 Min(alpha.water) = -4.52004e-13 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.961255 Min(alpha.water) = -3.40531e-07 Max(alpha.water) = 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0476108, No Iterations 632
time step continuity errors : sum local = 508582, global = -1554.68, cumulative = -1554.68
DICPCG: Solving for p_rgh, Initial residual = 1.52076e-11, Final residual = 1.52076e-11, No Iterations 0
time step continuity errors : sum local = 9.09126e+13, global = 402788, cumulative = 401233
DICPCG: Solving for p_rgh, Initial residual = 1.57267e-11, Final residual = 1.57267e-11, No Iterations 0
time step continuity errors : sum local = 9.40157e+13, global = 1.78437e+12, cumulative = 1.78437e+12
[1] [5] #0 Foam::error:rintStack(Foam::Ostream&)#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[5] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #2 ? at ??:?
[5] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
at ??:?
[1] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const[5] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[1] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[5] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[1] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[5] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[5] #7 at ??:?
[1] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::diction ary const&) at ??:?
[5] #8 at ??:?
[1] #8 Foam::fvMatrix<double>::solve() at ??:?
[5] #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&)Foam::fvMatrix<double>::solve() at ??:?
[5] #10 Foam::RASModels::LRR<Foam::IncompressibleTurbulenc eModel<Foam::transportModel> >::correct() at ??:?
[1] #9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
[5] #11 at ??:?
[1] #10 Foam::RASModels::LRR<Foam::IncompressibleTurbulenc eModel<Foam::transportModel> >::correct()? at ??:?
[1] #11 at ??:?
[5] #12 __libc_start_main? in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #13 at ??:?
[1] #12 __libc_start_main? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #13 at ??:?
[DNCO:05536] *** Process received signal ***
[DNCO:05536] Signal: Floating point exception (8)
[DNCO:05536] Signal code: (-6)
[DNCO:05536] Failing at address: 0x3e8000015a0
?[DNCO:05536] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) [0x7f2942ef8cb0]
[DNCO:05536] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f2942ef8c37]
[DNCO:05536] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) [0x7f2942ef8cb0]
[DNCO:05536] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6sm oothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5 _RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduI nterfaceFieldEEEhi+0x2f5) [0x7f29440b4d35]
[DNCO:05536] [ 4] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6s moothERNS_5FieldIdEERKS2_hi+0x2d) [0x7f29440b510d]
[DNCO:05536] [ 5] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5 FieldIdEERKS2_h+0x3a7) [0x7f29440afc97]
[DNCO:05536] [ 6] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegr egatedERKNS_10dictionaryE+0x11e) [0x7f29463c7afe]
[DNCO:05536] [ 7] interFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictio naryE+0x108) [0x49a528]
[DNCO:05536] [ 8] interFoam(_ZN4Foam8fvMatrixIdE5solveEv+0xcb) [0x49a7cb]
[DNCO:05536] [ 9] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam5solv eIdEENS_17SolverPerformanceIT_EERKNS_3tmpINS_8fvMa trixIS2_EEEE+0x27) [0x7f29479d31b7]
[DNCO:05536] [10] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASM odels3LRRINS_29IncompressibleTurbulenceModelINS_14 transportModelEEEE7correctEv+0x533) [0x7f2947a2a343]
[DNCO:05536] [11] interFoam() [0x439aa9]
[DNCO:05536] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f2942ee3f45]
[DNCO:05536] [13] interFoam() [0x43e479]
[DNCO:05536] *** End of error message ***
at ??:?
[DNCO:05532] *** Process received signal ***
[DNCO:05532] Signal: Floating point exception (8)
[DNCO:05532] Signal code: (-6)
[DNCO:05532] Failing at address: 0x3e80000159c
[DNCO:05532] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) [0x7fe1059ebcb0]
[DNCO:05532] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7fe1059ebc37]
[DNCO:05532] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36cb0) [0x7fe1059ebcb0]
[DNCO:05532] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6sm oothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5 _RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduI nterfaceFieldEEEhi+0x2f5) [0x7fe106ba7d35]
[DNCO:05532] [ 4] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6s moothERNS_5FieldIdEERKS2_hi+0x2d) [0x7fe106ba810d]
[DNCO:05532] [ 5] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5 FieldIdEERKS2_h+0x3a7) [0x7fe106ba2c97]
[DNCO:05532] [ 6] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegr egatedERKNS_10dictionaryE+0x11e) [0x7fe108ebaafe]
[DNCO:05532] [ 7] interFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictio naryE+0x108) [0x49a528]
[DNCO:05532] [ 8] interFoam(_ZN4Foam8fvMatrixIdE5solveEv+0xcb) [0x49a7cb]
[DNCO:05532] [ 9] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam5solv eIdEENS_17SolverPerformanceIT_EERKNS_3tmpINS_8fvMa trixIS2_EEEE+0x27) [0x7fe10a4c61b7]
[DNCO:05532] [10] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASM odels3LRRINS_29IncompressibleTurbulenceModelINS_14 transportModelEEEE7correctEv+0x533) [0x7fe10a51d343]
[DNCO:05532] [11] interFoam() [0x439aa9]
[DNCO:05532] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7fe1059d6f45]
[DNCO:05532] [13] interFoam() [0x43e479]
[DNCO:05532] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 5532 on node DNCO exited on signal 8 (Floating point exception).

[ebtedaei]

Please send me your email address.

Thanks
Ali