GitLab

Summary

I am trying to use turbulentDFSEMInlet at the inlet for simulation flow over aerofoil (The domain is attached below domain size ( 5c upstream, 10 c downstream and spanwise 0.1 c ) c is the chord length (c=0.3048 m) and in the log file at the first time step only 5 eddies are generated at the inlet as shown below

`starting time loop

Courant Number mean: 0.000449944 max: 0.0468918 deltaT = 1.2e-07 Time = 1.2e-07

PIMPLE: iteration 1

Turbulent DFSEM patch: inlet seeded 5 eddies with total volume 0.0519412

smoothSolver: Solving for Ux, Initial residual = 3.70186e-06, Final residual = 5.58285e-09, No Iterations 1

smoothSolver: Solving for Uy, Initial residual = 2.02752e-06, Final residual = 7.27509e-08, No Iterations 1

smoothSolver: Solving for Uz, Initial residual = 0.000726498, Final residual = 8.2538e-08, No Iterations 1

GAMG: Solving for p, Initial residual = 0.999923, Final residual = 1.86291, No Iterations 1000

GAMG: Solving for p, Initial residual = 0.0407024, Final residual = 0.0142608, No Iterations 1000

time step continuity errors : sum local = 9.09561e-07, global = -1.32729e-08, cumulative = -1.32729e-08

GAMG: Solving for p, Initial residual = 0.0149751, Final residual = 0.00549751, No Iterations 1000

GAMG: Solving for p, Initial residual = 0.0055067, Final residual = 0.00168609, No Iterations 1000

time step continuity errors : sum local = 1.49201e-07, global = -4.97451e-09, cumulative = -1.82474e-08`

Steps to reproduce

boundary conditions

`/--------------------------------- C++ -----------------------------------\

| ========= | |

| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |

| \ / O peration | Version: v1906 |

| \ / A nd | Web: www.OpenFOAM.com |

| \/ M anipulation | |

*---------------------------------------------------------------------------*/

FoamFile

{

version     2.0;

format      ascii;

class       volVectorField;

object      U;

} // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [ 0 1 -1 0 0 0 0 ];

internalField uniform ( 71.3 0 0 );

boundaryField

{

aerofoil


{

    type            fixedValue;

    value           uniform ( 0 0 0 );

}

top

{

    type            symmetryPlane;

}

bottom

{
    type            symmetryPlane;

}

front

{
    type            cyclic;

}

back

{

    type            cyclic;

}


inlet
{
    type            turbulentDFSEMInlet;
    delta           2;
    nCellPerEddy    3;
    mapMethod       nearestCell;
    R               uniform (0.5 0 0 0.5 0 0.5);
    U               uniform (71.3 0 0);
    L               uniform 0.055;
    value           $internalField;
}
outlet    
{
    type            freestream;
    freestreamValue uniform (71.3 0 0);
    value           uniform (71.3 0 0);
}

}

// ************************************************************************* // FoamFile

{

version     2.0;
format	ascii;
class	dictionary;
location    "constant";
object	turbulenceProperties;

} // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType LES;

LES { LESModel kEqn; turbulence on; printCoeffs on;

delta vanDriest;

vanDriestCoeffs { delta cubeRootVol; cubeRootVolCoeffs { deltaCoeff 1; }

Aplus           26;
Cdelta          0.158;

} } // ************************************************************************* //`

Example case

What is the current bug behaviour?

only 5 eddies are generated at the inlet and inlet velocity higher than average specified velocity as shown in the attached picture (velocity inlet=113 although the Umean=71.3 m/s)

What is the expected correct behavior?

<expected higher numbers of eddies are generated at the inlet with umean=71.3 m/s

Relevant logs and/or images

Paste any relevant logs - please use code blocks (```) to format console output, logs, and code as it's very hard to read otherwise. -->

Environment information

• OpenFOAM version :1906
• Operating system :ubuntu 15.04
• Hardware info :
• Compiler :

Possible fixes

## Reattaching the author to the issue ticket: @zein.elserfy ##