fvSchemes for DES

saeed jamshidi · November 7, 2023, 12:01

Hello Foamers,

I'm going to use DES for my simulation. I have read that in DES the scheme is specified using:

Code:

divSchemes
{
    default         none;
    div(phi,U)      Gauss DEShybrid
        linear                        // scheme 1
        linearUpwind grad(U)          // scheme 2
        0.65                          // DES coefficient, typically = 0.65
        30                            // Reference velocity scale
        2                             // Reference length scale
        0                             // Minimum sigma limit (0-1)
        1                             // Maximum sigma limit (0-1)
        1e-3;                         // Limiter of B function, typically 1e-03
}

Besides, in order to use the scheme, libs ("libturbulenceModelSchemes.so") must be added to controlDict.

Althouth I'v done this procedure, I,v come across fatal error:

Code:

file: stream.divSchemes.div(phi,U) at line 0.
[2]
[2]     From static Foam::tmp<Foam::surfaceInterpolationScheme<Type> > Foam::surfaceInterpolationScheme<Type>::New(const Foam::fvMesh&, const surfaceScalarField&, Foam::Istream&) [with Type = Foam::Vector<double>; Foam::surfaceScalarField = Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>]
[2]     in file lnInclude/surfaceInterpolationScheme.C at line 114.
[2]
FOAM parallel run exiting
[2]
[3]
[3]
[3] --> FOAM FATAL IO ERROR: (openfoam-2112 patch=220610)
[3] Unknown discretisation type DEShybrid

Valid discretisation types :

63
(
CoBlended
Gamma
GammaV
LUST
MUSCL
MUSCLV
Minmod
MinmodV
OSPRE
OSPREV
Phi
QUICK
QUICKV
SFCD
SFCDV
SuperBee
SuperBeeV
UMIST
UMISTV
biLinearFit
blended
cellCoBlended
clippedLinear
cubic
cubicUpwindFit
deferredCorrection
downwind
filteredLinear
filteredLinear2
filteredLinear2V
filteredLinear3
filteredLinear3V
fixedBlended
limitWith
limitedCubic
limitedCubicV
limitedLinear
limitedLinearV
limiterBlended
linear
linearFit
linearPureUpwindFit
linearUpwind
linearUpwindV
localBlended
localMax
localMin
midPoint
outletStabilised
pointLinear
quadraticFit
quadraticLinearFit
quadraticLinearUpwindFit
quadraticUpwindFit
reverseLinear
skewCorrected
upwind
vanAlbada
vanAlbadaV
vanLeer
vanLeerV
weighted
weightedFlux
)

How can I handle this problem?

thanks.

dlahaye · November 7, 2023, 12:37

What do tutorials do?

saeed jamshidi · November 7, 2023, 12:40

Dear dlahaye, thank you for the prompt reply.
Which tutorial, do you know which tutorial is the best fit?

dlahaye · November 7, 2023, 12:57

Any tutorial using DES on a sufficiently small mesh, I imagine.

Unix is your best friend.

Code:

tut 
cd incompressible 
for jfile in `find . -name turbulenceProperties` ; do echo $jfile; grep -i des $jfile; done

saeed jamshidi · November 7, 2023, 13:53

I cloudn't find the problem

Here are my fvschemes & fvsolution:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         backward;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss DEShybrid
        linear                        // scheme 1
        linearUpwind grad(U)          // scheme 2
        0.65                          // DES coefficient, typically = 0.65
        1                            // Reference velocity scale
        0.04                             // Reference length scale
        0                             // Minimum sigma limit (0-1)
        1                             // Maximum sigma limit (0-1)
        1e-3;                         // Limiter of B function, typically 1e-03
		
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,omega)  Gauss limitedLinear 1;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method          meshWave;
    nRequired       yes;
}


// ************************************************************************* //

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          GAMG;
        smoother        DICGaussSeidel;
        tolerance       1e-06;
        relTol          0.05;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    "(U|k|omega)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0.1;
    }

    "(U|k|omega)Final"
    {
        $U;
        relTol          0;
    }
	/*
    cellDisplacement
    {
        solver          GAMG;
        tolerance       1e-6;
        relTol          0;
        smoother        GaussSeidel;
    }
	*/
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       0;
}

PIMPLE
{
    nOuterCorrectors 1;
    nCorrectors      3;
    nNonOrthogonalCorrectors 1;
    pRefCell        0;
    pRefValue       0;
}

relaxationFactors
{
    fields
    {
        p               0.3;
    }
    equations
    {
        "(U|k|omega)"   0.7;
        "(U|k|omega)Final" 1.0;
    }
}

cache
{
    grad(U);
}

// ************************************************************************* //

I will appreciate any assistance.

dlahaye · November 8, 2023, 01:54

I do suggest running a tutorial that uses DES first.

Use Unix tools to find such a tutorial.

saeed jamshidi · November 8, 2023, 05:16

Dear dlahaye, thank you again for the reply.

I looked at all of my tutorials and there was not any DES turbulence model!!
Furtheremore, by reading below thread I'm going to accept that maybe my openfoam v2112 does not support DEShybrid !!

k-omega SST IDDES Setup

Yann · November 8, 2023, 05:41

Hello Saeed,

As Domenico suggested, you can use Unix tools to search for things in the tutorials.
For instance, this should lead you to tutorials using DEShydrid:

Code:

grep -r "DEShybrid" $FOAM_TUTORIALS

Yann

saeed jamshidi · November 8, 2023, 06:05

Thank all of you.

I saw them, and for the periodic hill case SpalartAllmarasIDDES turbulent model is used with the fvschem of:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         backward;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
}

divSchemes
{
    default         none;

    div(phi,U)      Gauss DEShybrid
        linear                    // scheme 1
        linearUpwind grad(U)      // scheme 2
        hmax
        0.65                      // DES coefficient, typically = 0.65
        1                         // Reference velocity scale
        0.028                     // Reference length scale
        0                         // Minimum sigma limit (0-1)
        1                         // Maximum sigma limit (0-1)
        1; // 1.0e-03;                  // Limiter of B function, typically 1e-03

    div(phi,k)      Gauss limitedLinear 1;
    div(phi,nuTilda) Gauss limitedLinear 1;

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

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method          meshWave;
    nRequired       yes;
}


// ************************************************************************* //

and I'm going to use kOmegaSSTIDDES. However, my fvscheme is:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         backward;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss DEShybrid
        linear                        // scheme 1
        linearUpwind grad(U)          // scheme 2
        0.65                          // DES coefficient, typically = 0.65
        1                            // Reference velocity scale
        0.04                             // Reference length scale
        0                             // Minimum sigma limit (0-1)
        1                             // Maximum sigma limit (0-1)
        1e-3;                         // Limiter of B function, typically 1e-03
		
    div(phi,k)      Gauss limitedLinear 1;
    div(phi,omega)  Gauss limitedLinear 1;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method          meshWave;
    nRequired       yes;
}


// ************************************************************************* //

major difference is in divSchemes part, where I should consider omega instead of nutilda (SpalartAllmaras model).

I am wondering about this really!!!!

dlahaye · November 8, 2023, 07:20

It should suffice to extend

Code:

div(phi,k)      Gauss limitedLinear 1;

to

Code:

div(phi,k)              Gauss limitedLinear 1;
div(phi,omega)      Gauss limitedLinear 1;

saeed jamshidi · November 9, 2023, 04:09

Finally, I found the problem, and I want to share it for readers of this thread.

Please, refer to here:

k-omega SST IDDES Setup

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Error in fvSchemes	TGS	OpenFOAM	0	May 10, 2023 11:20
fvschemes for laminar particles of fluidisedBed	gryphaea1635	OpenFOAM Pre-Processing	1	April 26, 2020 16:10
First order in fvSchemes does not seem to work	gerritgroot	OpenFOAM Running, Solving & CFD	0	September 30, 2015 20:06
Unrealistic values in cells at edge (high pressure & velocity) - fvSchemes?!	matzbanni	OpenFOAM Running, Solving & CFD	8	August 24, 2015 16:46
fvschemes and fvsolutions in MRFSimpleFoam	renyun0511	OpenFOAM Running, Solving & CFD	23	August 3, 2011 04:07

November 7, 2023, 12:37		#2
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 723 Blog Entries: 1 Rep Power: 17	What do tutorials do?

November 7, 2023, 12:40		#3
saeed jamshidi Senior Member Saeed Jamshidi Join Date: Aug 2019 Posts: 214 Rep Power: 7	Dear dlahaye, thank you for the prompt reply. Which tutorial, do you know which tutorial is the best fit?

November 7, 2023, 12:57		#4
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 723 Blog Entries: 1 Rep Power: 17	Any tutorial using DES on a sufficiently small mesh, I imagine. Unix is your best friend. Code: tut cd incompressible for jfile in `find . -name turbulenceProperties` ; do echo $jfile; grep -i des $jfile; done

November 8, 2023, 01:54		#6
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 723 Blog Entries: 1 Rep Power: 17	I do suggest running a tutorial that uses DES first. Use Unix tools to find such a tutorial.

November 8, 2023, 05:16		#7
saeed jamshidi Senior Member Saeed Jamshidi Join Date: Aug 2019 Posts: 214 Rep Power: 7	Dear dlahaye, thank you again for the reply. I looked at all of my tutorials and there was not any DES turbulence model!! Furtheremore, by reading below thread I'm going to accept that maybe my openfoam v2112 does not support DEShybrid !! k-omega SST IDDES Setup

November 8, 2023, 05:41		#8
Yann Senior Member Yann Join Date: Apr 2012 Location: France Posts: 1,067 Rep Power: 26	Hello Saeed, As Domenico suggested, you can use Unix tools to search for things in the tutorials. For instance, this should lead you to tutorials using DEShydrid: Code: grep -r "DEShybrid" $FOAM_TUTORIALS Yann

November 8, 2023, 07:20		#10
dlahaye Senior Member Domenico Lahaye Join Date: Dec 2013 Posts: 723 Blog Entries: 1 Rep Power: 17	It should suffice to extend Code: div(phi,k) Gauss limitedLinear 1; to Code: div(phi,k) Gauss limitedLinear 1; div(phi,omega) Gauss limitedLinear 1;

November 9, 2023, 04:09		#11
saeed jamshidi Senior Member Saeed Jamshidi Join Date: Aug 2019 Posts: 214 Rep Power: 7	Finally, I found the problem, and I want to share it for readers of this thread. Please, refer to here: k-omega SST IDDES Setup