[Inquisitor101]

I was trying some simple tutorial about a backward-facing step in steady state via different turbulence models in RAS using the simpleFoam solver. One of which was the LRR (Reynolds stress model). The tutorials themselves were done so in another version of OpenFOAM and the one I am using was the latest version: 4.1.

All else worked splendidly, however in the case of the LRR, I keep getting an error returned from the fvSchemes file. Basically, the steps in this method that I am following are: after copying the kEpsilon file (by file I mean the whole kEpsilon simulation folder -- i.e. the directory containing: 0/ constant/ system/ ), I execute the command:

simpleFoam -postProcess -func R

I get, instead of the expected R file, a file named turbulenceProperties:R in the 0/ directory. Moving on, I rename it to R.

Next, I change the turbulence model from kEpsilon to LRR as intended. Afterwards, I edit the fvSchemes and include the following for the divergence terms:

divSchemes
{
    default             none;
    div(phi,U)         bounded Gauss linearUpwind grad(U);
    div(phi,k)          bounded Gauss linearUpwind grad(k);
    div(phi,epsilon) bounded Gauss linearUpwind grad(epsilon);
    div(phi,omega)  bounded Gauss limitedLinear 1;
    div(phi,v2)        bounded Gauss limitedLinear 1;
    div(phi,R)         bounded Gauss upwind;
    div(R)              Gauss linear;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
    div(nonlinearStress)             Gauss linear;
}

Then, after running the command simpleFoam, I get the following error:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 0.01
    field U     tolerance 0.001
    field "(k|epsilon|omega|f|v2|R)"     tolerance 0.001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model LRR


--> FOAM FATAL IO ERROR: 
cannot find file

file: /home/username/OpenFOAM/username-4.1/run/tutorials/incompressible/simpleFoam/bfs/BFS_baseCase/0/R at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 72.

FOAM exiting

Then, after trying different stuff, I came up with the following and it worked. I simply change the default condition (in the divSchemes) from a none into some other scheme (e.g. Gauss linear) and all worked splendid.


My question hence is why did this happen? is there anything wrong with my implementation or is it a bug ?

Thanks !

[cryabroad]

Take a look at this thread

fvSchemes

According to Dr. Jasak, "It does not recognise it because this is a straightforward div(thing) term and not convection div(phi, thing). Therefore you cannot do upwinding, since there's no flux to decide the direction from."

[Giacomodb]

Thank you very much!! You fixed my problem.


I saw browsing from the internet but I didn't realized it was different from the one:

div((nuEff*dev2(T(grad(U))))) Gauss linear;

Later I got another error:

keyword wallDist is undefined in dictionary

I googled it and i added at the end of fvSolutions the following definition:

wallDist
{
    method meshWave;
}

Then it worked fine!


Thank you very much again!!


Giacomo

[lpz456]

Originally Posted by Inquisitor101

I was trying some simple tutorial about a backward-facing step in steady state via different turbulence models in RAS using the simpleFoam solver. One of which was the LRR (Reynolds stress model). The tutorials themselves were done so in another version of OpenFOAM and the one I am using was the latest version: 4.1.

All else worked splendidly, however in the case of the LRR, I keep getting an error returned from the fvSchemes file. Basically, the steps in this method that I am following are: after copying the kEpsilon file (by file I mean the whole kEpsilon simulation folder -- i.e. the directory containing: 0/ constant/ system/ ), I execute the command:

simpleFoam -postProcess -func R

I get, instead of the expected R file, a file named turbulenceProperties:R in the 0/ directory. Moving on, I rename it to R.

Next, I change the turbulence model from kEpsilon to LRR as intended. Afterwards, I edit the fvSchemes and include the following for the divergence terms:

divSchemes
{
    default             none;
    div(phi,U)         bounded Gauss linearUpwind grad(U);
    div(phi,k)          bounded Gauss linearUpwind grad(k);
    div(phi,epsilon) bounded Gauss linearUpwind grad(epsilon);
    div(phi,omega)  bounded Gauss limitedLinear 1;
    div(phi,v2)        bounded Gauss limitedLinear 1;
    div(phi,R)         bounded Gauss upwind;
    div(R)              Gauss linear;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
    div(nonlinearStress)             Gauss linear;
}

Then, after running the command simpleFoam, I get the following error:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 0.01
    field U     tolerance 0.001
    field "(k|epsilon|omega|f|v2|R)"     tolerance 0.001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model LRR


--> FOAM FATAL IO ERROR: 
cannot find file

file: /home/username/OpenFOAM/username-4.1/run/tutorials/incompressible/simpleFoam/bfs/BFS_baseCase/0/R at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 72.

FOAM exiting

Then, after trying different stuff, I came up with the following and it worked. I simply change the default condition (in the divSchemes) from a none into some other scheme (e.g. Gauss linear) and all worked splendid.


My question hence is why did this happen? is there anything wrong with my implementation or is it a bug ?

Thanks !

Hi，I was using LRR model to simulate pitzDaily by simpleFoam solver, but the results are unsteay, have you meet the same question?

[Ellen Caroline]

Originally Posted by Inquisitor101

I was trying some simple tutorial about a backward-facing step in steady state via different turbulence models in RAS using the simpleFoam solver. One of which was the LRR (Reynolds stress model). The tutorials themselves were done so in another version of OpenFOAM and the one I am using was the latest version: 4.1.

All else worked splendidly, however in the case of the LRR, I keep getting an error returned from the fvSchemes file. Basically, the steps in this method that I am following are: after copying the kEpsilon file (by file I mean the whole kEpsilon simulation folder -- i.e. the directory containing: 0/ constant/ system/ ), I execute the command:

simpleFoam -postProcess -func R

I get, instead of the expected R file, a file named turbulenceProperties:R in the 0/ directory. Moving on, I rename it to R.

Next, I change the turbulence model from kEpsilon to LRR as intended. Afterwards, I edit the fvSchemes and include the following for the divergence terms:

divSchemes
{
    default             none;
    div(phi,U)         bounded Gauss linearUpwind grad(U);
    div(phi,k)          bounded Gauss linearUpwind grad(k);
    div(phi,epsilon) bounded Gauss linearUpwind grad(epsilon);
    div(phi,omega)  bounded Gauss limitedLinear 1;
    div(phi,v2)        bounded Gauss limitedLinear 1;
    div(phi,R)         bounded Gauss upwind;
    div(R)              Gauss linear;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
    div(nonlinearStress)             Gauss linear;
}

Then, after running the command simpleFoam, I get the following error:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 0.01
    field U     tolerance 0.001
    field "(k|epsilon|omega|f|v2|R)"     tolerance 0.001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model LRR


--> FOAM FATAL IO ERROR: 
cannot find file

file: /home/username/OpenFOAM/username-4.1/run/tutorials/incompressible/simpleFoam/bfs/BFS_baseCase/0/R at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 72.

FOAM exiting

Then, after trying different stuff, I came up with the following and it worked. I simply change the default condition (in the divSchemes) from a none into some other scheme (e.g. Gauss linear) and all worked splendid.


My question hence is why did this happen? is there anything wrong with my implementation or is it a bug ?

Thanks !

Is it possible to create a R file if I didn't use the command "- postProcess - func R"? My simulation is already done and a I just typed "twoPhaseEulerFoam" in the terminal to initialize it.