CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

simpleFoam: Reynolds Stress, turbulenceProperties:R, no output

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 26, 2019, 08:57
Post simpleFoam: Reynolds Stress, turbulenceProperties:R, no output
  #1
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
Hello,

I am having troubles with extracting the Reynolds stress from a simulation. Wether I use

Code:
"simpleFoam -postProcess -func R"
or (in the controlDict)

Code:
 functions
{
    turbulenceFields1
    {
        type        turbulenceFields;
        libs        ("libfieldFunctionObjects.so");
        field		R;
        writeControl        onEnd;
    }
}
The result entire simulation executes well, but the extracted Reynolds stresses is the same:

Code:
No MRF models present

No finite volume options present
    functionObjects::turbulenceFields R writing field: turbulenceProperties:R
Which seems fine, except there is not output. Nothing to be found neither in the folder nor when viewing results via Paraview.

Running more postproccessing together with the Reynolds stress via the controlDict, for instance wallshear stress, works fine and produces:

Code:
    functionObjects::turbulenceFields turbulenceFields1 writing field: turbulenceProperties:R
wallShearStress wallShearStress1 write:
    writing field wallShearStress
    min/max(ground) = (-0.119511319929 -1.42561710357e-18 -2.02731918065e-05), (-0.115200085848 2.42927606607e-18 1.24279923199e-05)
But again, no Reynolds stresses. Which I need.

I am using the "Windows Binary v1812" version of OpenFOAM.

Does anyone have suggestions on what to do, to get the Reynolds stress?
Or a way of fixing my problem?
FoamingSimon is offline   Reply With Quote

Old   May 9, 2019, 03:46
Unhappy
  #2
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
Still no solutions in hand, but more information gathered:

When running either

Code:
postProcess -func R
or

Code:
postProcess -funcs '(turbulenceFields)'
The result is:

Code:
Executing functionObjects


--> FOAM FATAL ERROR:
Turbulence model not found in database, deactivating

    From function bool Foam::functionObjects::turbulenceFields::compressible()
    in file turbulenceFields/turbulenceFields.C at line 107.

FOAM exiting
FoamingSimon is offline   Reply With Quote

Old   May 15, 2019, 13:55
Default
  #3
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 3
blackpingu is on a distinguished road
Try this maybe :

Code:
 YourSolverHere -postProcess -func R
blackpingu is offline   Reply With Quote

Old   May 16, 2019, 02:41
Default
  #4
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
Hello Blackpingu

Unfortunately, it does not make any difference. An excerpt of the result from running above code is:

Code:
Time = 450
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 kEpsilon
kEpsilonCoeffs
{
    Cmu             0.03;
    C1              1.21;
    C2              1.92;
    sigmaEps        1.3;
    C3              0;
    sigmak          1;
}

No MRF models present

No finite volume options present
    functionObjects::turbulenceFields R writing field: turbulenceProperties:R
Also, still no "R" in the postprocess folder
FoamingSimon is offline   Reply With Quote

Old   May 16, 2019, 04:23
Default
  #5
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 3
blackpingu is on a distinguished road
Quote:
Originally Posted by FoamingSimon View Post
Hello Blackpingu

Unfortunately, it does not make any difference. An excerpt of the result from running above code is:

Code:
Time = 450
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 kEpsilon
kEpsilonCoeffs
{
    Cmu             0.03;
    C1              1.21;
    C2              1.92;
    sigmaEps        1.3;
    C3              0;
    sigmak          1;
}

No MRF models present

No finite volume options present
    functionObjects::turbulenceFields R writing field: turbulenceProperties:R
Also, still no "R" in the postprocess folder


Check your time folders for a file named "turbulenceProperties:R", in folder "450" for example.
blackpingu is offline   Reply With Quote

Old   May 16, 2019, 05:24
Default
  #6
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
No such file nor anything similar.
If I postprocess for yplus for example, I do get a "yPlus"-file with the results.
FoamingSimon is offline   Reply With Quote

Old   May 16, 2019, 05:56
Default
  #7
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 3
blackpingu is on a distinguished road
That is strange, what case are you running?
blackpingu is offline   Reply With Quote

Old   May 16, 2019, 07:03
Default
  #8
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
Indeed.

It does not matter what I run, it seems.
But currently, I run a modified "turbine siting tutorial". Basically a square domain, with ABL inflow.

I have attached a 2D-ABL flow study I have made, where I cannot extract the Reynold's stresses as well.
Attached Files
File Type: zip 2D_ABL_flow.zip (15.3 KB, 2 views)
FoamingSimon is offline   Reply With Quote

Old   May 16, 2019, 09:40
Default
  #9
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 3
blackpingu is on a distinguished road
Try to change to this in the controldict file instead of "onEnd"

Code:
turbulenceFields1
    {
        type        turbulenceFields;
        libs        ("libfieldFunctionObjects.so");
        field		R;
        writeControl        writeTime;
    }

or

Code:
turbulenceFields1
    {
        type        turbulenceFields;
        libs        ("libfieldFunctionObjects.so");
        field		R;
        writeControl        timeStep;
    }
blackpingu is offline   Reply With Quote

Old   May 16, 2019, 10:25
Default
  #10
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
Thank you for your perseverance, Blackpingu, but no - changing the writeControl to either writeTime or timeStep did not result in any Reynolds stress outputs. Neither in the postProcess folder, nor in any of the arbitrary timestep folders.

Although, the writeControl "timeStep" did result in a lot of empty folders - which makes sense, since an output is written for every time increment.

Do you have any other ideas?
FoamingSimon is offline   Reply With Quote

Old   May 16, 2019, 10:30
Default
  #11
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 3
blackpingu is on a distinguished road
Quote:
Originally Posted by FoamingSimon View Post
Thank you for your perseverance, Blackpingu, but no - changing the writeControl to either writeTime or timeStep did not result in any Reynolds stress outputs. Neither in the postProcess folder, nor in any of the arbitrary timestep folders.

Although, the writeControl "timeStep" did result in a lot of empty folders - which makes sense, since an output is written for every time increment.

Do you have any other ideas?
Not that i can think of, but I wasn't able to run your case properly, didn't get files for U or P or anything, do you get decent results for the velocity? Which OF version are you running?
blackpingu is offline   Reply With Quote

Old   May 16, 2019, 11:40
Default
  #12
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
That sounds weird, as I get reasonable results, as can be seen on the attached .png-file, showing the velocity profiles through the domain.

- Did you run it via the Allrun-file?
Otherwise, you would have to copy the 0.orig and rename it 0 before running it.

I am using the "Windows Binary v1812" version of OpenFOAM.
Attached Images
File Type: png 4_again_VelocityProfiles.png (38.7 KB, 8 views)
FoamingSimon is offline   Reply With Quote

Old   May 16, 2019, 15:25
Default
  #13
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 3
blackpingu is on a distinguished road
Quote:
Originally Posted by FoamingSimon View Post
That sounds weird, as I get reasonable results, as can be seen on the attached .png-file, showing the velocity profiles through the domain.

- Did you run it via the Allrun-file?
Otherwise, you would have to copy the 0.orig and rename it 0 before running it.

I am using the "Windows Binary v1812" version of OpenFOAM.
Yep I did but no luck
blackpingu is offline   Reply With Quote

Old   May 17, 2019, 02:20
Default
  #14
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
Quote:
Originally Posted by blackpingu View Post
Yep I did but no luck
Maybe this is due to some differences in versions of OF.
What version do you use?
And is it normally possible for you to extract the Reynolds stresses?
FoamingSimon is offline   Reply With Quote

Old   May 17, 2019, 02:45
Default
  #15
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
New error type..
Using:
Code:
postProcess -func 'components(turbulenceProperties:R)'
The result is:
Code:
Time = 0

Reading fields:

Executing functionObjects
--> FOAM Warning :     functionObjects::components components(turbulenceProperties:R) failed to execute.
For each time..
FoamingSimon is offline   Reply With Quote

Old   May 19, 2019, 14:37
Default
  #16
New Member
 
Join Date: Sep 2018
Posts: 27
Rep Power: 3
blackpingu is on a distinguished road
I use OpenFOAM 5, and yes works everytime.
blackpingu is offline   Reply With Quote

Old   October 10, 2019, 07:24
Default
  #17
New Member
 
Join Date: Mar 2019
Posts: 10
Rep Power: 2
FoamingSimon is on a distinguished road
Switched OF version from the windows binary docker version to ubuntu via windows, and everything works like a charm now.
FoamingSimon is offline   Reply With Quote

Reply

Tags
reynolds stress, v1812

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Turbulent DFSEM reynolds stress calculation ssa_cfd OpenFOAM Running, Solving & CFD 0 February 1, 2019 08:23
[ANSYS Meshing] Error in Workbench "Unable to attach to geometry file" Jeremie84 ANSYS Meshing & Geometry 70 October 25, 2017 16:03
New implemented algebraic Reynolds stress model volker OpenFOAM 13 October 21, 2013 00:01
Reynolds Stress Models Jade M Main CFD Forum 0 April 21, 2010 16:38
Reynolds Stress Models Roued Main CFD Forum 20 February 8, 2000 02:58


All times are GMT -4. The time now is 11:06.