CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

turbulentDFSEMInlet for LES imulation for naca0012 aerofoil

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 22, 2019, 23:03
Default turbulentDFSEMInlet for LES imulation for naca0012 aerofoil
  #1
New Member
 
zein elserfy
Join Date: May 2018
Posts: 25
Rep Power: 7
zeinelserfy is on a distinguished road
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 ?
Attached Images
File Type: jpg Capture.JPG (21.3 KB, 28 views)
zeinelserfy is offline   Reply With Quote

Old   December 4, 2020, 11:50
Default
  #2
New Member
 
Haoran Shi
Join Date: Jun 2020
Posts: 1
Rep Power: 0
Haoran is on a distinguished road
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.
Haoran is offline   Reply With Quote

Reply

Tags
les inlet conditions, openfoam 1906, turbulentdfseminlet

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
2D Aerofoil with rotating cylinders - Help & Advice! Bluejay OpenFOAM 1 December 31, 2012 05:40


All times are GMT -4. The time now is 10:59.