zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

TUT: add kEpsilonPhitF RAS model example to bump2D

Browse Source 
- STYLE: remove redundant whitespaces
This commit is contained in:
Kutalmis Bercin
2019-12-17 11:44:55 +00:00
parent e38e5d6717
commit e6238a5c62
9 changed files with 195 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

60
tutorials/incompressible/simpleFoam/bump2D/0/epsilon Normal file
View File

@ -0,0 +1,60 @@
/*--------------------------------*- 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 volScalarField;
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -3 0 0 0 0 ];
internalField uniform 469.872; // computed from omega and k
boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}
top
{
type symmetryPlane;
}
"(symUp|symDown)"
{
type symmetryPlane;
}
bump
{
type epsilonWallFunction;
value $internalField;
lowReCorrection on;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

56
tutorials/incompressible/simpleFoam/bump2D/0/f Normal file
View File

@ -0,0 +1,56 @@
/*--------------------------------*- 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 volScalarField;
object f;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 0 -1 0 0 0 0 ];
internalField uniform 1;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
top
{
type symmetryPlane;
}
"(symUp|symDown)"
{
type symmetryPlane;
}
bump
{
type fixedValue;
value uniform 0;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

2
tutorials/incompressible/simpleFoam/bump2D/0/nuTilda
View File

@ -45,7 +45,7 @@ boundaryField
bump
{
type fixedValue;
type fixedValue;
value uniform 0;
}

2
tutorials/incompressible/simpleFoam/bump2D/0/p
View File

@ -43,7 +43,7 @@ boundaryField
bump
{
type zeroGradient;
type zeroGradient;
}
frontAndBack

56
tutorials/incompressible/simpleFoam/bump2D/0/phit Normal file
View File

@ -0,0 +1,56 @@
/*--------------------------------*- 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 volScalarField;
object phit;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 0 0 0 0 0 0 ];
internalField uniform 0.66;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
top
{
type symmetryPlane;
}
"(symUp|symDown)"
{
type symmetryPlane;
}
bump
{
type fixedValue;
value uniform 1e-10;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

1
tutorials/incompressible/simpleFoam/bump2D/Allrun
View File

@ -10,6 +10,7 @@ runApplication blockMesh
models="
kOmegaSST
SpalartAllmaras
kEpsilonPhitF
"
# Extract a value (Eg, from boundaryField/bump/value)

2
tutorials/incompressible/simpleFoam/bump2D/plot
View File

@ -4,7 +4,7 @@ cd "${0%/*}" || exit # Run from this directory
#------------------------------------------------------------------------------
# Turbulence closure models
models=('kOmegaSST' 'SpalartAllmaras')
models=('kOmegaSST' 'SpalartAllmaras' 'kEpsilonPhitF')
# Note: CFL3D data available from:
# https://turbmodels.larc.nasa.gov/bump_sa.html

10
tutorials/incompressible/simpleFoam/bump2D/system/fvSchemes
View File

@ -27,12 +27,14 @@ gradSchemes
divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwind grad(U);
default none;
div(phi,U) bounded Gauss linearUpwind grad(U);
div(phi,epsilon) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,omega) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,nuTilda) bounded Gauss upwind;
div(phi,phit) bounded Gauss upwind;
div(phi,f) bounded Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

16
tutorials/incompressible/simpleFoam/bump2D/system/fvSolution
View File

@ -25,13 +25,21 @@ solvers
relTol 0.1;
}
"(U|k|epsilon|omega|nuTilda)"
"(U|k|epsilon|omega|nuTilda|phit)"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-8;
relTol 0;
}
f
{
solver PBiCGStab;
preconditioner DIC;
tolerance 1e-8;
relTol 0;
}
}
SIMPLE
@ -42,7 +50,7 @@ SIMPLE
{
p 1e-5;
U 1e-5;
"(k|epsilon|omega|nuTilda)" 1e-5;
"(k|epsilon|omega|nuTilda|phit|f)" 1e-5;
}
}
@ -51,7 +59,7 @@ relaxationFactors
equations
{
U 0.9;
"(k|epsilon|omega|nuTilda)" 0.7;
"(k|epsilon|omega|nuTilda|phit|f)" 0.7;
}
}
@ -59,6 +67,8 @@ relaxationFactors
cache
{
grad(U);
grad(k);
grad(phit);
}
// ************************************************************************* //