tutorials/incompressible/simpleFoam/pitzDaily: Setup to run with range of turbulence models

This commit is contained in:
Henry 2015-01-28 18:36:03 +00:00
parent 06e4f4db6e
commit acbd73cb5b
7 changed files with 166 additions and 16 deletions

View File

@ -0,0 +1,49 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object v2;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
upperWall
{
type fWallFunction;
value uniform 0;
}
lowerWall
{
type fWallFunction;
value uniform 0;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

View File

@ -0,0 +1,50 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 39.6;
boundaryField
{
inlet
{
type fixedValue;
value uniform 39.6;
}
outlet
{
type zeroGradient;
}
upperWall
{
type omegaWallFunction;
value uniform 39.6;
}
lowerWall
{
type omegaWallFunction;
value uniform 39.6;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

View File

@ -0,0 +1,50 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object v2;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.25;
boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type zeroGradient;
}
upperWall
{
type v2WallFunction;
value $internalField;
}
lowerWall
{
type v2WallFunction;
value $internalField;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

View File

@ -19,7 +19,7 @@ simulationType RAS;
RAS
{
RASModel kEpsilon;
RASModel v2f; //kEpsilon; //realizableKE; //kOmega; //kOmegaSST;
turbulence on;

View File

@ -17,19 +17,19 @@ FoamFile
application simpleFoam;
startFrom latestTime;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 1000;
endTime 2000;
deltaT 1;
writeControl timeStep;
writeInterval 50;
writeInterval 100;
purgeWrite 0;

View File

@ -28,12 +28,11 @@ gradSchemes
divSchemes
{
default none;
div(phi,U) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,R) bounded Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) bounded Gauss upwind;
div(phi,U) bounded Gauss linearUpwind grad(U);
div(phi,k) bounded Gauss limitedLinear 1;
div(phi,epsilon) bounded Gauss limitedLinear 1;
div(phi,omega) bounded Gauss limitedLinear 1;
div(phi,v2) bounded Gauss limitedLinear 1;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}
@ -58,5 +57,10 @@ fluxRequired
p ;
}
wallDist
{
method meshWave;
}
// ************************************************************************* //

View File

@ -31,7 +31,7 @@ solvers
mergeLevels 1;
}
"(U|k|epsilon|R|nuTilda)"
"(U|k|epsilon|omega|f|v2)"
{
solver smoothSolver;
smoother symGaussSeidel;
@ -48,7 +48,7 @@ SIMPLE
{
p 1e-2;
U 1e-3;
"(k|epsilon|omega)" 1e-3;
"(k|epsilon|omega|f|v2)" 1e-3;
}
}
@ -61,10 +61,7 @@ relaxationFactors
equations
{
U 0.7;
k 0.7;
epsilon 0.7;
R 0.7;
nuTilda 0.7;
".*" 0.5;
}
}