[zeinelserfy]

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)

inlet condition

    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.17;
        value           uniform (0 0 0);
    }

the turbulenceProperties file

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\	  /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:	www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
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;
}
}
// ************************************************************************* //

After running the simulation for the first time step, the only 5 eddy is generated at the inlet

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

Does anyone have used this boundary condition for external flow case ?

[Haoran]

Hello,

Maybe it is too late. I met the same problem. And I found that even though the OpenFOAM website says we can set constant values for U, L and R, but it does not work. For my case, U and L can be set as constant, but for R, I have to use the mapping method, otherwise it will stop at the first iter. And the deltaT is extremly small.