CFD Online Discussion Forums - How to use DESModelRegions function object

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

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

- - How to use DESModelRegions function object (https://www.cfd-online.com/Forums/openfoam-solving/117083-how-use-desmodelregions-function-object.html)

hakonbar

May 1, 2013 18:37

How to use DESModelRegions function object

Hi! I'm doing some DES simulations, and I would like to see which parts of my flow field are in RANS mode and which are in LES mode. I see the new OF release (2.2.0) introduced a new function object called DESModelRegions, which outputs this data. The thing is, I don't know which lines to put in the controlDict in order to activate this object, and which library to load.

Have any of you used this function object before, and do you know which libs to load etc?

best regards,
Håkon Line

wyldckat

May 2, 2013 08:48

Greetings Håkon,

Without an example case, I can't test any further than the following instructions:

Edit the file "system/controlDict".
Append/add the following block:

Code:

functions { desField { type DESModelRegions; functionObjectLibs ("libutilityFunctionObjects.so"); log true; //region "region0"; enabled on; storeFilter on; timeStart 0.0; timeEnd 10.0; outputControl timeStep; outputInterval 5000; } }

Note: it writes and outputs information for each 5000 steps in this example! And if you do not define the "timeStart" and "timeEnd", it will run for the whole simulation.
It created the file "postProcessing/desField/0/DESModelRegions.dat", but I didn't use DES, so I only got this header:

Code:

# DES model region coverage (% volume) # time LES RAS
And since I used "log true;", I only got this message:

Code:

DESModelRegions output: No DES turbulence model found in database
And it creates the field "DES::LESRegion".

Like I wrote at the beginning, without an example, I can't test this properly.

edit: I do have an example case now and is used here: http://www.cfd-online.com/Forums/ope...tml#post424797 - post #10

Best regards,
Bruno

hakonbar

May 3, 2013 08:15

It works perfectly, thanks! Now I've learned the general syntax for function objects as well. =)

potentialFoam

July 5, 2016 07:26

Dear Foamers,

can you help me to use 'DESModelRegions' with OF30X?

I fear, it is not implemented in OF30X. But how could you determine the RANS content of a DES?

Using an older version (like OF231) for postprocessing does not work, 'cos the constant/turbulenceProperties-file changed since this version.

Even this thread

HTML Code:

http://www.cfd-online.com/Forums/openfoam-solving/104123-les-content-detached-eddy-simulation-openfoam.html

does not provide an answer for OF30X.

I would be happy to hear your ideas!
Regards,
Peter

EDIT:
According to a hint I installed OpenFOAM-v1606+ from openfoam.com
and this perfectly works for postprocessing cases simulated with OF301 :-)

Bazinga

December 2, 2016 03:38

Quote:

Originally Posted by potentialFoam (Post 608048)

HTML Code:

http://www.cfd-online.com/Forums/openfoam-solving/104123-les-content-detached-eddy-simulation-openfoam.html

Dear Peter,

I am also trying to find my LES and RANS regions for my DES simulations. Can you guide me how you used 1606+ to find those regions. Thanks

potentialFoam

December 2, 2016 04:01

Sure, inside
system/controlDict.functions
you need to add

Code:

functions

{

 //   Determine LES-content

    desField

    {

      type            DESModelRegions;

      functionObjectLibs ("libutilityFunctionObjects.so");

      log             true;

      enabled on;

      storeFilter on;

      outputControl   outputTime;

    }

}

And afterwards, you can just execute

Code:

execFlowFunctionObjects -time 0.4 | tee log.desField

for the time you want.

Good luck!

Bazinga

December 2, 2016 04:09

Quote:

Originally Posted by potentialFoam (Post 627911)

Sure, inside
system/controlDict.functions
you need to add

Code:

functions

{

 //   Determine LES-content

    desField

    {

      type            DESModelRegions;

      functionObjectLibs ("libutilityFunctionObjects.so");

      log             true;

      enabled on;

      storeFilter on;

      outputControl   outputTime;

    }

}

And afterwards, you can just execute

Code:

execFlowFunctionObjects -time 0.4 | tee log.desField

for the time you want.

Good luck!

Thank you so much for the fast reply. Will try it now. Best regards :)

Bazinga

December 2, 2016 09:39

Hmmm, I installed v1606+ but it says "execFlowFunctionObjects: command not found".

My OpenFOAM 2.4.0 installation can be used for "execFlowFunctionObjects" but for this the folder and file structure has changed too much since, so it doesn't work.

Anything I am missing here?

wyldckat

December 3, 2016 15:18

Quote:

Originally Posted by Bazinga (Post 627970)

Hmmm, I installed v1606+ but it says "execFlowFunctionObjects: command not found".

Quick answer: Something might have went wrong with your installation of v1606+, because execFlowFunctionObjects is built by default in that version as well.

Or maybe you had OpenFOAM 4.0 or 4.1 active in your shell? Because it was replaced by in 4.0 by postProcess, along with (almost) all solvers having their one "-postProcess" argument as well.

Bazinga

December 4, 2016 01:55

Thanks for the quick reply. I think something went wrong with the installation. Will try again :)

edit: So, you were right. Not everything was installed during the 1606+ installation. I was able to run wmake for the execFlowFunctionObjects myself and now it is working :)

hanis_rda

August 14, 2018 06:18

Dear Foamer,

I have done IDDES simulation in OpenFoam 5x.
But when I try to use DESModelRegions function for post-processing in OpenFoam 5x. It is not recognised.
Anyone can help me to find the region that used RANS and LES using OpenFOAM 5x?

Thank you.

wyldckat

August 17, 2018 16:11

Quick answer: https://bugs.openfoam.org/view.php?id=2763

hanis_rda

August 17, 2018 16:25

Quote:

Originally Posted by wyldckat (Post 703033)

Quick answer: https://bugs.openfoam.org/view.php?id=2763

Thank you for the information.
Is there any other function that can differentiate between RANS and LES region?

wyldckat

August 20, 2018 20:45

Greetings hanis_rda,

I was curious about this and was inspired to work on this. I've managed to create the function object for OpenFOAM 5 and 6, but I have not fully tested it.

The repository is provided here: https://github.com/wyldckat/DESModelRegions
And you can find download and build instructions on that page as well.

Please report back here if this works or not and with which OpenFOAM version it worked and with which DES-type turbulence model!

Best regards,
Bruno

hanis_rda

August 21, 2018 09:49

Quote:

Originally Posted by wyldckat (Post 703323)

Thank you very much for your help.
I run for IDDES turbulence model in OpenFOAM5x.
But I got this error message:

Code:



DESModelRegions desField write:

    writing field DES::LESRegion

--> FOAM Warning : 

    From function Foam::OFstream& Foam::functionObjects::logFiles::file()

    in file db/functionObjects/logFiles/logFiles.C at line 131

    Requested single file, but multiple files are present

--------------------------------------------------------------------------

A process has executed an operation involving a call to the

"fork()" system call to create a child process.  Open MPI is currently

operating in a condition that could result in memory corruption or

other system errors; your job may hang, crash, or produce silent

data corruption.  The use of fork() (or system() or other calls that

create child processes) is strongly discouraged.



The process that invoked fork was:



  Local host:          [[16608,0],0] (PID 1396)



If you are *absolutely sure* that your application will successfully

and correctly survive a call to fork(), you may disable this warning

by setting the mpi_warn_on_fork MCA parameter to 0.

--------------------------------------------------------------------------

[0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?

[0] #1  Foam::sigSegv::sigHandler(int) at ??:?

[0] #2  ? in "/usr/lib64/libc.so.6"

[0] #3  Foam::functionObjects::logFiles::file() at ??:?

[0] #4  Foam::functionObjects::DESModelRegions::write() at ??:?

[0] #5  Foam::functionObjects::timeControl::write() at ??:?

[0] #6  Foam::functionObjectList::execute() at ??:?

[0] #7  Foam::Time::loop() at ??:?

[0] #8  ? at ??:?

[0] #9  __libc_start_main in "/usr/lib64/libc.so.6"

[0] #10  ? at ??:?

wyldckat

August 21, 2018 10:18

Quick answer: Sorry about that, I forgot to properly re-test after adapting the missing features that were in the original function object.

I've now fixed it and is available at the repository. To update your build:

If you downloaded by using git, then go into the folder "DESModelRegions" and run:

Code:

git pull wmake wclean
If you downloaded from ZIP, then you must delete the old ZIP file and the folder "DESModelRegions-master". Then do once again the steps for downloading the ZIP file and so on, as instructed on the repository.

hanis_rda

August 21, 2018 10:37

Now it is running for IDDES simulation in OpenFOAM5x. Thank you very much.

# DES model region coverage (% volume)
# Time LES RAS
1.4001 5.134440e+01 4.865560e+01
1.4002 5.134438e+01 4.865562e+01
1.4003 5.134434e+01 4.865566e+01

S123	September 18, 2018 17:27

Receiving error: Attempt to cast type cubeRootVol to type IDDESDelta

Hello, I am trying to run a SpalartAllmarasIDDES simulation in OpenFoam v5.0. However, I receive the error "Attempt to cast type cubeRootVol to type IDDESDelta". Since you were able to successfully run IDDES in v5, please could you provide some assistance. cubeRootVol is the delta notation in my LESProperties file.

This is how my LESProperties file looks:

Code:

simulationType  LES;



LES

{

        LESModel            SpalartAllmarasIDDES;



        delta               cubeRootVol;



        printCoeffs         on;



    turbulence          on;        



        cubeRootVolCoeffs

{

    deltaCoeff      1;

}





laminarCoeffs

{

}



oneEqEddyCoeffs

{

    ck               0.07;

    ce               1.05;

}



dynOneEqEddyCoeffs

{

    ce               1.05;

    filter          simple;

}



locDynOneEqEddyCoeffs

{

    ce               1.05;

    filter          simple;

}



SmagorinskyCoeffs

{

    ce               1.05;

    ck               0.07;

}



Smagorinsky2Coeffs

{

    ce               1.05;

    ck               0.07;

    cD2              0.02;

}



spectEddyViscCoeffs

{

    ce               1.05;

    cB               8.22;

    cK1              0.83;

    cK2              1.03;

    cK3              4.75;

    cK4              2.55;

}



dynSmagorinskyCoeffs

{

    ce               1.05;

    filter          simple;

}



mixedSmagorinskyCoeffs

{

    ce               1.05;

    ck               0.07;

    filter          simple;

}



dynMixedSmagorinskyCoeffs

{

    ce               1.05;

    filter          simple;

}



LRRDiffStressCoeffs

{

    ce               1.05;

    ck               0.09;

    c1               1.8;

    c2               0.6;

}



DeardorffDiffStressCoeffs

{

    ce               1.05;

    ck               0.09;

    cm               4.13;

}



SpalartAllmarasCoeffs

{

    alphaNut         1.5;

    Cb1              0.1355;

    Cb2              0.622;

    Cw2              0.3;

    Cw3              2;

    Cv1              7.1;

    Cv2              5.0;

    CDES             0.65;

    ck               0.07;

}



SpalartAllmarasDDESCoeffs

{

    alphaNut         1.5;

    Cb1              0.1355;

    Cb2              0.622;

    Cw2              0.3;

    Cw3              2.0;

    Cv1              7.1;

    Cv2              5.0;

    CDES             0.65;

    ck               0.07;

}



SpalartAllmarasIDDESCoeffs

{

     delta  IDDESDelta;



     IDDESDeltaCoeffs

     {

         deltaCoeff       1;



         //according to src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C

         cw               0.15;

     }



     alphaNut         1.5;

     kappa            0.4187;

     Cb1              0.1355;

     Cb2              0.622;

     Cw2              0.3;

     Cw3              2.0;

     Cv1              7.1;

     Cv2              5.0;

     CDES             0.65;

     ck               0.07;



//These are the default values

     sigmaNut        0.66666;

     kappa           0.41;

     Cb1             0.1355;

     Cb2             0.622;

     Cv1             7.1;

     Cv2             5;

     CDES            0.65;

     ck              0.07;

     Cw2             0.3;

     Cw3             2;

     fwStar          0.424;

     cl              3.55;

     ct              1.63;



}



cubeRootVolCoeffs

{

    deltaCoeff      1;

}



PrandtlCoeffs

{

    delta           cubeRootVol;

    cubeRootVolCoeffs

    {

        deltaCoeff      1;

    }

    smoothCoeffs

    {

        delta           cubeRootVol;

        cubeRootVolCoeffs

        {

            deltaCoeff      1;

        }

        maxDeltaRatio   1.1;

    }

    Cdelta           0.158;

}



vanDriestCoeffs

{

    delta           cubeRootVol;

    cubeRootVolCoeffs

    {

        deltaCoeff      1;

    }

    smoothCoeffs

    {

        delta           cubeRootVol;

        cubeRootVolCoeffs

        {

            deltaCoeff      1;

        }

        maxDeltaRatio   1.1;

    }

    Aplus            26;

    Cdelta           0.158;

}



smoothCoeffs

{

    delta           cubeRootVol;

    cubeRootVolCoeffs

    {

        deltaCoeff      1;

    }

    maxDeltaRatio   1.1;

}



kappa            0.4187;



wallFunctionCoeffs

{

    E                9;

}

Many thanks.

Quote:

Originally Posted by hanis_rda (Post 703398)

hanis_rda

September 20, 2018 04:53

Quote:

Originally Posted by S123 (Post 706750)

Sorry for late reply.

I have changed all 'cubeRootVol' to 'IDDESDelta' in my LESProperties. Hope it will help you.

S123	September 21, 2018 05:52

Thank you very much, Hanis. This worked!! Many thanks & best wishes.

Quote:

Originally Posted by hanis_rda (Post 706893)

Sorry for late reply.

I have changed all 'cubeRootVol' to 'IDDESDelta' in my LESProperties. Hope it will help you.

ch_jaehnel

September 26, 2018 03:24

Quote:

Originally Posted by potentialFoam (Post 627911)

Sure, inside
system/controlDict.functions
you need to add

Code:

functions

{

 //   Determine LES-content

    desField

    {

      type            DESModelRegions;

      functionObjectLibs ("libutilityFunctionObjects.so");

      log             true;

      enabled on;

      storeFilter on;

      outputControl   outputTime;

    }

}

And afterwards, you can just execute

Code:

execFlowFunctionObjects -time 0.4 | tee log.desField

for the time you want.

Good luck!

Thank you for that hint. I did a simulation with interfoam and now get the following error after executing "execFlowFunctionObjects -time 1.5 | tee log.desField":

Code:

Create time



Create mesh for time = 1.5



Time = 1.5

    Reading phi

    Reading U

    Reading p

No finite volume options present



--> FOAM Warning :

--> FOAM FATAL IO ERROR:

keyword transportModel is undefined in dictionary "/lustre/scratch2/s2665038/FFW2.5_n24_350l/constant/transportProperties"



file: /lustre/scratch2/s2665038/FFW2.5_n24_350l/constant/transportProperties from line 18 to line 66.



    From function dictionary::lookupEntry(const word&, bool, bool) const

    in file db/dictionary/dictionary.C at line 437.



End

Here are my transportProperties:

Code:

phases (water air);



water

{

    transportModel  Newtonian;

    nu              nu [ 0 2 -1 0 0 0 0 ]  1e-06;

    rho             rho [ 1 -3 0 0 0 0 0 ] 1000;

    

    CrossPowerLawCoeffs    

    {

        nu0            nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;

        nuInf        nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;

        m            m [ 0 0 1 0 0 0 0 ] 1;

        n            n [ 0 0 0 0 0 0 0 ] 0;

    }

    

    BirdCarreauCoeffs   

    {

        nu0            nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;

        nuInf        nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;

        k            k [ 0 0 1 0 0 0 0 ] 99.6;

        n            n [ 0 0 0 0 0 0 0 ] 0.1003;

    }

}



air

{

    transportModel  Newtonian;

    nu              nu [ 0 2 -1 0 0 0 0 ]  1.48e-05;

    rho             rho [ 1 -3 0 0 0 0 0 ] 1;

    

    CrossPowerLawCoeffs    

    {

        nu0            nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;

        nuInf        nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;

        m            m [ 0 0 1 0 0 0 0 ] 1;

        n            n [ 0 0 0 0 0 0 0 ] 0;

    }

    

    BirdCarreauCoeffs    

    {

        nu0            nu0 [ 0 2 -1 0 0 0 0 ] 0.0142515;

        nuInf        nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;

        k            k [ 0 0 1 0 0 0 0 ] 99.6;

        n            n [ 0 0 0 0 0 0 0 ] 0.1003;

    }

}



sigma           sigma [ 1 0 -2 0 0 0 0 ] 0.07;

so the transportModel is obviously defined, but for every phase. How can I get around this error?
Thanks a lot!

wyldckat

September 26, 2018 16:51

Quick question @ch_jaehnel: Which OpenFOAM version are you using?

ch_jaehnel

September 28, 2018 06:35

Quote:

Originally Posted by wyldckat (Post 707729)

Quick question @ch_jaehnel: Which OpenFOAM version are you using?

Sorry forgot to mention: I use OF 2.3.1

I have also acces to OF5.0, the result is the same.

saiguruprasad

December 4, 2018 02:58

Quote:

Originally Posted by wyldckat (Post 703323)

I am trying to use the DESModelRegions for the k-Omega SST DES model, but I don't think it is supported with this code.

Is there any way to make this code work for the k-Omega SST DES model?

I got this error when I used it for k-Omega SST DES : "No DES turbulence model found in database"

wyldckat

December 4, 2018 16:41

Quick answers:

@ch_jaehnel: Sorry for the very late reply:

Quote:

Originally Posted by ch_jaehnel (Post 708099)

Sorry forgot to mention: I use OF 2.3.1

I have also acces to OF5.0, the result is the same.

The problem is that execFlowFunctionObjects does not load the necessary dependencies for transport properties, turbulence and so on.
With OpenFOAM 5, use the "-postProcess" option with the solver, which will load in the necessary properties and fields. See the User Guide for more details.

----------

@saiguruprasad:

Quote:

Originally Posted by saiguruprasad (Post 717845)

I... I don't even know what to tell you exactly... please read this thread: https://www.cfd-online.com/Forums/op...-get-help.html - hopefully you will understand why I don't even know what answer to give you :(

saiguruprasad

December 4, 2018 22:19

I am sorry for not being verbose about the problem before.

I am trying to simulate Couette flow at Re_w=8600 using k- $\omega$ SST DES model in OF-5.x. I used perturbUChannel to set up the initial conditions for the flow to sustain turbulence. Running a LES simulation, the instantaneous and mean velocity should be different. But I was not able to achieve this using k- $\omega$ SST DES.

I wanted to use DESModelRegions to check how much of the area was being solved using LES and RANS. I downloaded the code and compiled it without errors. I added the folowing line to my controlDict:

Code:

functions

{

desField

    {

        type            DESModelRegions;

        functionObjectLibs ("libDESModelRegions.so");

        writeControl writeTime;

    }

}

In the log file, I got the following output:

Quote:

DESModelRegions desField write:
No DES turbulence model found in database

I would like to know if I can use DESModelRegions with k- $\omega$ SST DES model in OF-5.x.

wyldckat

December 5, 2018 05:46

Quick question: What is the exact line you used to define the turbulence model in the file "turbulenceProperties"?

saiguruprasad

December 5, 2018 22:39

Code:



simulationType  LES;



LES

{

    LESModel        kOmegaSSTDES;



    turbulence      on;



    printCoeffs     on;



    delta               cubeRootVol;



    printCoeffs         on;



    turbulence          on;



    cubeRootVolCoeffs

    {

        deltaCoeff      1;

    }



    PrandtlCoeffs

    {

        delta           cubeRootVol;

        cubeRootVolCoeffs

        {

            deltaCoeff      1;

        }

        smoothCoeffs

        {

            delta           cubeRootVol;

            cubeRootVolCoeffs

            {

                deltaCoeff      1;

            }

            maxDeltaRatio   1.1;

        }

        Cdelta           0.158;

    }



    vanDriestCoeffs

    {

        delta           cubeRootVol;

        cubeRootVolCoeffs

        {

            deltaCoeff      1;

        }

        smoothCoeffs

        {

            delta           cubeRootVol;

            cubeRootVolCoeffs

            {

                deltaCoeff      1;

            }

            maxDeltaRatio   1.1;

        }

        Aplus            26;

        Cdelta           0.158;

    }



    smoothCoeffs

    {

        delta           cubeRootVol;

        cubeRootVolCoeffs

        {

            deltaCoeff      1;

        }

        maxDeltaRatio   1.1;

    }

}

This is my turbulenceProperties file

wyldckat

December 6, 2018 12:38

Quick answer @saiguruprasad: Many thanks for all of the details, because I'm not very familiar with DES and how it's implemented.

After taking a better look, I roughly remembered things and then noticed what I wrote back when I created the repository and I quote: https://github.com/wyldckat/DESModel...#how-to-use-it

Quote:

You must use the turbulence model SpalartAllmarasDES or one of its derivatives, as shown in the diagram here: https://cpp.openfoam.org/v5/classFoa...lmarasDES.html

In other words, "kOmegaSSTDES" is not supported because it does not provide the necessary calculations for differentiate between RAS and LES :(

In more detail, with the SpalartAllmarasDES models, the calculation of the model regions field is calculated as indicated here: https://cpp.openfoam.org/v5/SpalartA...ce.html#l00379
But the problem is that it's not clear how it should be calculated for "kOmegaSSTDES", not even back in OpenFOAM 2.3.

If you do some research on the topic and figure out how it should be calculated, we can then work on implementing that calculation.

sibo	February 21, 2019 10:44

DESModelRegions Usage

1 Attachment(s)

Quote:

Originally Posted by wyldckat (Post 703323)

Hi Bruno,

Firstly thanks a lot for providing this utility for OF 5.
I want to test this function in a tutorial case pitzDaily using SpalartAllmarasDES model. But the results show that

Code:

DESModelRegions desField write:

    No DES turbulence model found in database

I was following the instruction https://github.com/wyldckat/DESModelRegions to build it.

I attached my case here. Could you help me to test it?
Thanks!

wyldckat

February 24, 2019 18:10

Quick answer:

Quote:

Originally Posted by sibo (Post 725554)

I attached my case here. Could you help me to test it?

Many thanks for the feedback and test case! I had made a mistake in how the model is constructed for compressible flow.

If you downloaded the ZIP file, then please do another download and build again, because I have committed the bug fix just now.

If you downloaded using Git, then go into the folder where you've placed the "DESModelRegions" source code and run:

Code:

git pull

wmake

Then try running the case once again! It should now work as intended.

sibo	February 24, 2019 18:31

Thanks a lot!

It works fine now.

CarlosGRR

April 19, 2019 19:12

Implementation for kwSST-DES

Quote:

Originally Posted by wyldckat (Post 703323)

Thank you for this function for the latest OpenFOAM versions. As you mentioned in your GitHub documentation, this function only supports Spalart-Allmaras-DES models. It would be great that this were uptated to work with kwSST-DES hybrid models, available in OpenFOAM too.

Thanks again for this contribution,

Carlos

wyldckat

April 21, 2019 06:17

Quick answer:

Even though I have been informed and given permission to bring in adapted code to the repo I'm maintaining, I haven't had time to work on it: https://github.com/wyldckat/DESModelRegions/issues/1
On that report I was informed that the versions at OpenFOAM.com do have that capability, namely they do have the DESModelRegions for the k-omega variant with DES.

All times are GMT -4. The time now is 16:43.