[Anne Lincke]

I have the same in my computation, meaning to use cubeRootVol for delta and IDDESDelta inside SpalartAllmaras model.
This works in any case.
Would be interesting to know if it is possible to set the general delta to "IDDESDelta".
So far, this does not work with OF-2.1.1.

[MaryBau]

It runs perfect! Thank you ingcorra.

But I might have spoken to fast. Anne has a point.

So to check if IDDESDelta is being used and not cubeRoot I simply added a line in the IDDESDelta.C file:

Info<<delta_.internalField()<<endl;

after the calculation of delta.

It prints the value of delta before the first iteration, but not in each iteration.

Does delta is something that is calculated once? (which will mean the IDDESDelta works fine) or it is calculated every iteration? (and there is something wrong in the code)

[Joanna Huang]

Originally Posted by Anne Lincke

Hey Foamers,

I started with OpenFOAM®+DES, SpalartAllmarasIDDES.
I would like to determine the LES-Content of my computation.
Does someone know an indicator or a way to determine the content?

It would be desirable to see the averaged LES content...

Thanks
Anne

Hello,
Anne, I notice that you are familar with IDDES, I am currently doing some simulation relating to this, but everytime it has errors like this:


request for RASModel RASProperties from objectRegistry region0 failed
available objects of type RASModel are
0
(
)

I attached the case file below(0,constant,system), could you please kindly help me with this, I will really appreciate that since I am suffering this now.

[hakonbar]

Hi everyone! I was just wondering, what's the consensus on the delta for IDDES? The way I understand the issue is that in old releases of OpenFOAM (pre 2.1) it was possible to define the delta as "IDDESDelta" in the root of the "LESProperties" file (as per post #11), but that in the newer versions, the delta has to be defined as something else first (like in post #20).

Is this a bug in the new releases, or has it been done on purpose? If it's on purpose, is the first definition of delta just a "dummy definition" which isn't really used in the simulation?

I get the following output when I start an IDDES simulation, suggesting that the solver switches back and forth between different definitions of delta.

Selecting turbulence model type LESModel
Selecting LES turbulence model SpalartAllmarasIDDES
Selecting LES delta type maxDeltaxyz
Selecting LES delta type maxDeltaxyz
Selecting LES delta type IDDESDelta
Selecting LES delta type maxDeltaxyz

When I select a regular LES-model, there's no such switching back and forth:

Selecting turbulence model type LESModel
Selecting LES turbulence model oneEqEddy
Selecting LES delta type maxDeltaxyz

So is this a bug or a feature, and how do we find out?

[Anne Lincke]

The computation of IDDES delta is based on a general delta. You can find in IDDESDelta.C the formula of IDDESDelta

delta_.internalField() =
            deltaCoeff_
           *min
            (
                max
                (
                    max(cw_*wallDist(mesh()).y(), cw_*hmax),
                    faceToFacenMax()
                ),
                hmax
            );

where "hmax" needs a general delta as you can see in SpalartAllmarasIDDES.C

hmax_
    (
        LESdelta::New
        (
            "hmax",
            mesh_,
            *this
        )
    ),

The hmax equals to the maximum local grid spacing. When using structured meshes the defintion of maxDeltaxyz is fine, when using unstructured meshes the cubeRootVolDelta is more appropriate.
In the paper of Shur, Spalart, Strelets and Travin, International Jorunal of Heat and Fluid Flow 29 (2008) 1638-1649 you can read more.

Hope this helped.

Kind Regards,
Anne

[hakonbar]

Hi Anne, and thanks for the thorough explanation! So in a structured grid, like one would make in blockMesh, all the cells are aligned with the principal axes of the domain, and the "maxDeltaxyz" definition is fine, but since this might not be the case for an unstructured grid, it would be safer to use the "cubeRootVol" definition..is this correct?

Shur et al. do talk a lot about how the "max"-definition is preferable to the cube root, but the choice of delta definition is maybe not as significant for DDES and IDDES, which have the shielding term?

cheers,
Håkon

[vonboett]

Hi Håkon,

The LeMoS extensions for OpenFOAM provide an URANS-LES Detatched Eddy approach where the regions of URANS and LES are decomposed dynamically in every timestep dependent on flow properties and local mesh resolution. They have very good validation results using ship hydrodynamic test cases. So if you want to use LES with an unstructured grid, this might be helpful.

[Anne Lincke]

Hey Ha°kon,

So in a structured grid, like one would make in blockMesh, all the cells are aligned with the principal axes of the domain, and the "maxDeltaxyz" definition is fine, but since this might not be the case for an unstructured grid, it would be safer to use the "cubeRootVol" definition..is this correct?

Yes, this is correct.

Kind Regards
Anne

[baklanton]

Hi Anne,

Could you please post or attach the post-processing tool that you created, the one that writes \tilde{fd}*(1+fe) and 1-\tilde{fd} .

Many Thanks

Anton

[Anne Lincke]

Dear Anton,

sorry for this late response. I found the tool which I implemented a while ago. If you are still interested, I could send it to you via e-mail.

Kind Regards,
Anne

[Chanterz]

(Sorry, ignore me)

[Chanterz]

For anyone still following this thread, the correction Corrado posted in response to Anne's LESProperties post worked for me,

Thanks

[Alhasan]

Hi All,

My two cents on the intial question,

This should be added to the functions in the controlDict, to get the LES and DES regions during the simulations so it can be visualised after the simulation.

desField
{
   type DESModelRegions;
   functionObjectLibs ("libutilityFunctionObjects.so");
   log true;
   enabled on;
   storeFilter on;
   timeStart 0.15;
   outputControl outputTime;
}

Regards,
Hasan K.J