CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   turbulentDFSEMInlet for LES imulation for naca0012 aerofoil (https://www.cfd-online.com/Forums/openfoam-solving/221606-turbulentdfseminlet-les-imulation-naca0012-aerofoil.html)

zeinelserfy October 22, 2019 22:03

turbulentDFSEMInlet for LES imulation for naca0012 aerofoil
 
1 Attachment(s)
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
Code:

    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
Code:

/*--------------------------------*- 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
Code:

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 December 4, 2020 10:50

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.


All times are GMT -4. The time now is 17:47.