|
[Sponsors] |
simpleFoam: Reynolds Stress, turbulenceProperties:R, no output |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 26, 2019, 08:57 |
simpleFoam: Reynolds Stress, turbulenceProperties:R, no output
|
#1 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
Hello,
I am having troubles with extracting the Reynolds stress from a simulation. Wether I use Code:
"simpleFoam -postProcess -func R" Code:
functions { turbulenceFields1 { type turbulenceFields; libs ("libfieldFunctionObjects.so"); field R; writeControl onEnd; } } Code:
No MRF models present No finite volume options present functionObjects::turbulenceFields R writing field: turbulenceProperties:R 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) 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? |
|
May 9, 2019, 03:46 |
|
#2 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
Still no solutions in hand, but more information gathered:
When running either Code:
postProcess -func R Code:
postProcess -funcs '(turbulenceFields)' 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 |
|
May 15, 2019, 13:55 |
|
#3 |
New Member
Join Date: Sep 2018
Posts: 27
Rep Power: 8 |
Try this maybe :
Code:
YourSolverHere -postProcess -func R |
|
May 16, 2019, 02:41 |
|
#4 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
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 |
|
May 16, 2019, 04:23 |
|
#5 | |
New Member
Join Date: Sep 2018
Posts: 27
Rep Power: 8 |
Quote:
Check your time folders for a file named "turbulenceProperties:R", in folder "450" for example. |
||
May 16, 2019, 05:24 |
|
#6 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
No such file nor anything similar.
If I postprocess for yplus for example, I do get a "yPlus"-file with the results. |
|
May 16, 2019, 05:56 |
|
#7 |
New Member
Join Date: Sep 2018
Posts: 27
Rep Power: 8 |
That is strange, what case are you running?
|
|
May 16, 2019, 07:03 |
|
#8 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
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. |
|
May 16, 2019, 09:40 |
|
#9 |
New Member
Join Date: Sep 2018
Posts: 27
Rep Power: 8 |
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; } |
|
May 16, 2019, 10:25 |
|
#10 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
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? |
|
May 16, 2019, 10:30 |
|
#11 | |
New Member
Join Date: Sep 2018
Posts: 27
Rep Power: 8 |
Quote:
|
||
May 16, 2019, 11:40 |
|
#12 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
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. |
|
May 16, 2019, 15:25 |
|
#13 | |
New Member
Join Date: Sep 2018
Posts: 27
Rep Power: 8 |
Quote:
|
||
May 17, 2019, 02:20 |
|
#14 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
||
May 17, 2019, 02:45 |
|
#15 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
New error type..
Using: Code:
postProcess -func 'components(turbulenceProperties:R)' Code:
Time = 0 Reading fields: Executing functionObjects --> FOAM Warning : functionObjects::components components(turbulenceProperties:R) failed to execute. |
|
May 19, 2019, 14:37 |
|
#16 |
New Member
Join Date: Sep 2018
Posts: 27
Rep Power: 8 |
I use OpenFOAM 5, and yes works everytime.
|
|
October 10, 2019, 07:24 |
|
#17 |
New Member
Join Date: Mar 2019
Posts: 12
Rep Power: 7 |
Switched OF version from the windows binary docker version to ubuntu via windows, and everything works like a charm now.
|
|
Tags |
reynolds stress, v1812 |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |