turbulentDFSEMInlet for LES imulation for naca0012 aerofoil
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 :