Update tutorials

This commit is contained in:
Henry
2013-01-25 12:43:07 +00:00
parent 2eb05ebb87
commit 8bcb8132bf
18 changed files with 151 additions and 2379 deletions

View File

@ -107,26 +107,37 @@ edges
( (
); );
patches boundary
( (
// is there no way of defining all my 'defaultFaces' to be 'wall'? // is there no way of defining all my 'defaultFaces' to be 'wall'?
wall front front
{
type wall;
faces
( (
// inlet block // inlet block
frontQuad(in1, join1, join2, in2) frontQuad(in1, join1, join2, in2)
// outlet block // outlet block
frontQuad(poro1, out1, out2, poro2) frontQuad(poro1, out1, out2, poro2)
) );
}
wall back back
{
type wall;
faces
( (
// inlet block // inlet block
backQuad(in1, join1, join2, in2) backQuad(in1, join1, join2, in2)
// outlet block // outlet block
backQuad(poro1, out1, out2, poro2) backQuad(poro1, out1, out2, poro2)
) );
}
wall wall wall
{
type wall;
faces
( (
// inlet block // inlet block
quad2D(in1, join1) quad2D(in1, join1)
@ -134,9 +145,13 @@ patches
// outlet block // outlet block
quad2D(poro1, out1) quad2D(poro1, out1)
quad2D(out2, poro2) quad2D(out2, poro2)
) );
}
wall porosityWall porosityWall
{
type wall;
faces
( (
// porosity block // porosity block
frontQuad(join1, poro1, poro2, join2) frontQuad(join1, poro1, poro2, join2)
@ -145,17 +160,26 @@ patches
// porosity block // porosity block
quad2D(join1, poro1) quad2D(join1, poro1)
quad2D(poro2, join2) quad2D(poro2, join2)
) );
}
patch inlet inlet
{
type patch;
faces
( (
quad2D(in2, in1) quad2D(in2, in1)
) );
}
patch outlet outlet
{
type patch;
faces
( (
quad2D(out2, out1) quad2D(out2, out1)
) );
}
); );
mergePatchPairs mergePatchPairs

View File

@ -28,23 +28,17 @@ porosity1
DarcyForchheimerCoeffs DarcyForchheimerCoeffs
{ {
d d [0 -2 0 0 0 0 0] (5e7 -1000 -1000); d d [0 -2 0 0 0 0 0] (7e5 -1000 -1000);
f f [0 -1 0 0 0 0 0] (0 0 0); f f [0 -1 0 0 0 0 0] (0 0 0);
coordinateSystem coordinateSystem
{ {
type cartesian;
origin (0 0 0);
coordinateRotation
{
type axesRotation;
e1 (0.70710678 0.70710678 0); e1 (0.70710678 0.70710678 0);
e2 (0 0 1); e3 (0 0 1);
}
} }
} }
} }
} }
// ************************************************************************* // //************************************************************************ //

View File

@ -21,14 +21,14 @@ solvers
{ {
solver PCG; solver PCG;
preconditioner DIC; preconditioner DIC;
tolerance 1e-06; tolerance 1e-07;
relTol 0.01; relTol 0.01;
} }
pFinal pFinal
{ {
$p; $p;
tolerance 1e-06; tolerance 1e-07;
relTol 0; relTol 0;
} }
@ -36,14 +36,14 @@ solvers
{ {
solver PBiCG; solver PBiCG;
preconditioner DILU; preconditioner DILU;
tolerance 1e-05; tolerance 1e-06;
relTol 0.1; relTol 0.1;
} }
"(rho|U|h|k|epsilon|omega)Final" "(rho|U|h|k|epsilon|omega)Final"
{ {
$U; $U;
tolerance 1e-05; tolerance 1e-06;
relTol 0; relTol 0;
} }
@ -56,8 +56,8 @@ PIMPLE
nOuterCorrectors 50; nOuterCorrectors 50;
nCorrectors 1; nCorrectors 1;
nNonOrthogonalCorrectors 0; nNonOrthogonalCorrectors 0;
rhoMin rhoMin [ 1 -3 0 0 0 ] 0.1; rhoMin rhoMin [ 1 -3 0 0 0 ] 0.5;
rhoMax rhoMax [ 1 -3 0 0 0 ] 3.0; rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0;
residualControl residualControl
{ {
@ -75,13 +75,12 @@ relaxationFactors
{ {
fields fields
{ {
"p.*" 1; "p.*" 0.3;
"rho.*" 1; "rho.*" 1;
} }
equations equations
{ {
"(U|h|k|epsilon|omega).*" 0.85; "(U|h|k|epsilon|omega).*" 0.7;
"p.*" 1;
} }
} }

View File

@ -64,7 +64,7 @@ SIMPLE
{ {
p 1e-2; p 1e-2;
U 1e-4; U 1e-4;
T 1e-3; e 1e-3;
// possibly check turbulence fields // possibly check turbulence fields
"(k|epsilon|omega)" 1e-3; "(k|epsilon|omega)" 1e-3;
@ -76,7 +76,7 @@ relaxationFactors
fields fields
{ {
p 0.3; p 0.3;
rho 0.05; rho 0.03;
} }
equations equations
{ {

View File

@ -58,7 +58,7 @@ SIMPLE
{ {
p 1e-3; p 1e-3;
U 1e-4; U 1e-4;
T 1e-3; e 1e-3;
// possibly check turbulence fields // possibly check turbulence fields
"(k|epsilon|omega)" 1e-3; "(k|epsilon|omega)" 1e-3;

View File

@ -64,7 +64,7 @@ SIMPLE
{ {
p 1e-2; p 1e-2;
U 1e-4; U 1e-4;
T 1e-3; e 1e-3;
// possibly check turbulence fields // possibly check turbulence fields
"(k|epsilon|omega)" 1e-3; "(k|epsilon|omega)" 1e-3;
@ -76,7 +76,7 @@ relaxationFactors
fields fields
{ {
p 0.3; p 0.3;
rho 0.02; rho 0.05;
} }
equations equations
{ {

View File

@ -54,6 +54,16 @@ SIMPLE
rhoMin rhoMin [1 -3 0 0 0] 0.1; rhoMin rhoMin [1 -3 0 0 0] 0.1;
rhoMax rhoMax [1 -3 0 0 0] 1.0; rhoMax rhoMax [1 -3 0 0 0] 1.0;
transonic yes; transonic yes;
residualControl
{
p 1e-3;
U 1e-4;
e 1e-3;
// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
} }
relaxationFactors relaxationFactors

View File

@ -28,7 +28,7 @@ gradSchemes
divSchemes divSchemes
{ {
default none; default none;
div(rhoFlux,rho) Gauss limitedLinear 1; div(rhoFlux,rho) Gauss upwind;
} }
laplacianSchemes laplacianSchemes

View File

@ -10,36 +10,37 @@ FoamFile
version 2.0; version 2.0;
format ascii; format ascii;
class volScalarField; class volScalarField;
location "0";
object k; object k;
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0]; dimensions [0 2 -2 0 0 0 0];
internalField uniform 10.0; internalField uniform 10;
boundaryField boundaryField
{ {
inlet inlet
{ {
type fixedValue; type turbulentIntensityKineticEnergyInlet;
value uniform 0.5; intensity 0.05;
value uniform 0.05;
} }
outlet outlet
{ {
type zeroGradient; type zeroGradient;
} }
walls walls
{ {
type zeroGradient; type kqRWallFunction;
value uniform 10;
} }
frontBack frontBack
{ {
type empty; type empty;
} }
} }
// ************************************************************************* // // ************************************************************************* //

View File

@ -10,27 +10,30 @@ FoamFile
version 2.0; version 2.0;
format ascii; format ascii;
class volScalarField; class volScalarField;
location "0";
object nuSgs; object nuSgs;
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0]; dimensions [0 2 -1 0 0 0 0];
internalField uniform 1e-11; internalField uniform 0;
boundaryField boundaryField
{ {
inlet inlet
{ {
type zeroGradient; type calculated;
value uniform 0;
} }
outlet outlet
{ {
type zeroGradient; type calculated;
value uniform 0;
} }
walls walls
{ {
type fixedValue; type nutUSpaldingWallFunction;
value uniform 0; value uniform 0;
} }
frontBack frontBack
@ -39,4 +42,5 @@ boundaryField
} }
} }
// ************************************************************************* // // ************************************************************************* //

View File

@ -28,9 +28,9 @@ interpolationSchemes
divSchemes divSchemes
{ {
default none; default none;
div(phiv,rho) Gauss limitedLinear 0.2; div(phiv,rho) Gauss vanLeer;
div(phi,U) Gauss filteredLinear2V 0.2 0; div(phi,U) Gauss LUST grad(U);
div(phiv,k) Gauss filteredLinear2 0.2 0; div(phiv,k) Gauss LUST grad(k);
div((muEff*dev(T(grad(U))))) Gauss linear; div((muEff*dev(T(grad(U))))) Gauss linear;
} }

View File

@ -16,14 +16,15 @@ FoamFile
dimensions [0 2 -2 0 0]; dimensions [0 2 -2 0 0];
internalField uniform 10.0; internalField uniform 10;
boundaryField boundaryField
{ {
inlet inlet
{ {
type fixedValue; type turbulentIntensityKineticEnergyInlet;
value uniform 0.5; intensity 0.05;
value uniform 0.05;
} }
outlet outlet
@ -33,7 +34,8 @@ boundaryField
walls walls
{ {
type zeroGradient; type kqRWallFunction;
value uniform 10;
} }
frontBack frontBack

View File

@ -10,34 +10,38 @@ FoamFile
version 2.0; version 2.0;
format ascii; format ascii;
class volScalarField; class volScalarField;
location "0";
object nuSgs; object nuSgs;
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0]; dimensions [0 2 -1 0 0 0 0];
internalField uniform 1e-11; internalField uniform 0;
boundaryField boundaryField
{ {
inlet inlet
{ {
type zeroGradient; type calculated;
value uniform 0;
} }
outlet outlet
{ {
type zeroGradient; type calculated;
value uniform 0;
} }
walls walls
{ {
type fixedValue; type nutUSpaldingWallFunction;
value uniform 0; value uniform 0;
} }
frontBack frontBack
{ {
type fixedValue; type nutUSpaldingWallFunction;
value uniform 0; value uniform 0;
} }
} }
// ************************************************************************* // // ************************************************************************* //

View File

@ -16,14 +16,15 @@ FoamFile
dimensions [0 2 -2 0 0]; dimensions [0 2 -2 0 0];
internalField uniform 10.0; internalField uniform 10;
boundaryField boundaryField
{ {
inlet inlet
{ {
type fixedValue; type turbulentIntensityKineticEnergyInlet;
value uniform 0.5; intensity 0.05;
value uniform 0.05;
} }
outlet outlet
@ -33,7 +34,8 @@ boundaryField
walls walls
{ {
type zeroGradient; type kqRWallFunction;
value uniform 10;
} }
frontBack frontBack

View File

@ -10,34 +10,38 @@ FoamFile
version 2.0; version 2.0;
format ascii; format ascii;
class volScalarField; class volScalarField;
location "0";
object nuSgs; object nuSgs;
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0]; dimensions [0 2 -1 0 0 0 0];
internalField uniform 1e-11; internalField uniform 0;
boundaryField boundaryField
{ {
inlet inlet
{ {
type zeroGradient; type calculated;
value uniform 0;
} }
outlet outlet
{ {
type zeroGradient; type calculated;
value uniform 0;
} }
walls walls
{ {
type fixedValue; type nutUSpaldingWallFunction;
value uniform 0; value uniform 0;
} }
frontBack frontBack
{ {
type fixedValue; type nutUSpaldingWallFunction;
value uniform 0; value uniform 0;
} }
} }
// ************************************************************************* // // ************************************************************************* //

View File

@ -29,8 +29,8 @@ cp -r 0.org 0
runApplication blockMesh runApplication blockMesh
refineMeshByCellSet 1 2 3 refineMeshByCellSet 1 2 3
echo "mapping fields from 2D throttle case" #echo "mapping fields from 2D throttle case"
runApplication mapFields ../throttle -sourceTime latestTime #runApplication mapFields ../throttle -sourceTime latestTime
runApplication decomposePar runApplication decomposePar
runParallel $application 4 runParallel $application 4

View File

@ -28,9 +28,9 @@ interpolationSchemes
divSchemes divSchemes
{ {
default none; default none;
div(phiv,rho) Gauss limitedLinear 0.2; div(phiv,rho) Gauss vanLeer;
div(phi,U) Gauss filteredLinear2V 0.2 0; div(phi,U) Gauss LUST grad(U);
div(phiv,k) Gauss filteredLinear2 0.2 0; div(phiv,k) Gauss LUST grad(k);
div((muEff*dev(T(grad(U))))) Gauss linear; div((muEff*dev(T(grad(U))))) Gauss linear;
} }