Updated transonic tutorial

This commit is contained in:
Hrvoje Jasak 2014-09-02 12:39:34 +01:00 committed by Dominik Christ
parent 0c1ad0a78e
commit fd2c3bc914
24 changed files with 98 additions and 836 deletions

View file

@ -1,79 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 300;
boundaryField
{
INLE1
{
type fixedValue;
value $internalField;
// type waveTransmissiveInlet;
// phi phi;
// rho rho;
// psi psi;
// U U;
// gamma 1.4;
// inletOutlet false;
// correctSupercritical false;
// lInf 0.01;//HJ
// fieldInf 300;
// value $internalField;
}
PRES2
{
type zeroGradient;
// type waveTransmissive;
// phi phi;
// rho rho;
// psi psi;
// U U;
// gamma 1.4;
// inletOutlet true;
// correctSupercritical false;
// lInf 0.0;
// fieldInf 300;
// value $internalField;
}
WALL3
{
type zeroGradient;
}
SYMP5
{
type empty;
}
SYMP6
{
type empty;
}
WALL4
{
type zeroGradient;
}
}
// ************************************************************************* //

View file

@ -1,86 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
(-1.5 0 -0.1)
(-0.5 0 -0.1)
( 0.5 0 -0.1)
( 1.5 0 -0.1)
(-1.5 1 -0.1)
(-0.5 1 -0.1)
( 0.5 1 -0.1)
( 1.5 1 -0.1)
(-1.5 0 0.1)
(-0.5 0 0.1)
( 0.5 0 0.1)
( 1.5 0 0.1)
(-1.5 1 0.1)
(-0.5 1 0.1)
( 0.5 1 0.1)
( 1.5 1 0.1)
);
blocks
(
hex (0 1 5 4 8 9 13 12) (160 160 1) simpleGrading (1 1 1)
hex (1 2 6 5 9 10 14 13) (160 160 1) simpleGrading (1 1 1)
hex (2 3 7 6 10 11 15 14) (160 160 1) simpleGrading (1 1 1)
);
edges
(
arc 1 2 (0 0.1 -0.1)
arc 9 10 (0 0.1 0.1)
);
patches
(
patch INLE1
(
(0 8 12 4)
)
patch PRES2
(
(3 7 15 11)
)
wall WALL3
(
(4 12 13 5)
(5 13 14 6)
(6 14 15 7)
)
wall WALL4
(
(0 1 9 8)
(1 2 10 9)
(2 3 11 10)
)
);
mergePatchPairs
(
);
// ************************************************************************* //

View file

@ -1,76 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver BiCGStab;
preconditioner DILU;
minIter 0;
maxIter 1000;
tolerance 1e-8;
relTol 0.0;
};
U
{
solver BiCGStab;
preconditioner DILU;
minIter 0;
maxIter 1000;
tolerance 1e-8;
relTol 0.0;
};
h
{
solver BiCGStab;
preconditioner DILU;
minIter 0;
maxIter 1000;
tolerance 1e-8;
relTol 0.0;
};
}
PIMPLE
{
nOuterCorrectors 1;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
}
relaxationFactors
{
// Note: under-relaxation factors used in wave-transmissive schemes
U 0.4;
p 0.2;
h 0.5;
rho 0.5;
T 1;
}
fieldBounds
{
// With bounding
p 50 1e6;
T 20 3000;
U 1000;
}
// ************************************************************************* //

View file

@ -1,74 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 300;
boundaryField
{
INLE1
{
type fixedValue;
value $internalField;
// type waveTransmissiveInlet;
// phi phi;
// rho rho;
// psi psi;
// U U;
// gamma 1.4;
// inletOutlet false;
// correctSupercritical false;
// lInf 0.01;//HJ
// fieldInf 300;
// value $internalField;
}
PRES2
{
// type zeroGradient;
type waveTransmissive;
phi phi;
rho rho;
psi psi;
U U;
gamma 1.4;
inletOutlet true;
correctSupercritical false;
lInf 0.5;
fieldInf 300;
value $internalField;
}
WALL3
{
type zeroGradient;
}
WALL4
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //

View file

@ -1,65 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (173 0 0);
boundaryField
{
INLE1
{
type fixedValue;
value uniform (173 0 0);
}
PRES2
{
// type inletOutlet;
// inletValue uniform (0 0 0);
type waveTransmissive;
phi phi;
rho rho;
psi psi;
U U;
gamma 1.4;
inletOutlet true;
correctSupercritical false;
lInf 0.5;
fieldInf (0 0 0);
value $internalField;
}
WALL3
{
type fixedValue;
value uniform (0 0 0);
}
WALL4
{
type fixedValue;
value uniform (0 0 0);
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //

View file

@ -1,74 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 101325;
boundaryField
{
INLE1
{
type zeroGradient;
// type waveTransmissiveInlet;
// phi phi;
// rho rho;
// psi psi;
// U U;
// gamma 1.4;
// inletOutlet false;
// correctSupercritical false;
// lInf 0.0;
// fieldInf 101325;
// value uniform 101325;
}
PRES2
{
// type fixedValue;
// value uniform 101325;
type waveTransmissive;
phi phi;
rho rho;
psi psi;
U U;
gamma 1.4;
inletOutlet false;
correctSupercritical true;
lInf 0.5;
fieldInf 101325;
value uniform 101325;
}
WALL3
{
type zeroGradient;
}
WALL4
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //

View file

@ -1,7 +0,0 @@
#!/bin/sh
# Source tutorial clean functions
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
cleanCase
\rm -rf VTK

View file

@ -1,8 +0,0 @@
#!/bin/sh
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
application="steadyCompressibleFoam"
runApplication blockMesh
runApplication $application

View file

@ -1,32 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
RASModel laminar;
turbulence on;
printCoeffs on;
laminarCoeffs
{}
wallFunctionCoeffs
{
kappa 0.4187;
E 9;
}
// ************************************************************************* //

View file

@ -1,22 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType hPsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>;
// name, nMoles, mol weight, CP, Hf, mu, Pr;
mixture air 1 28.9 1007 0 0 0.7;
// ************************************************************************* //

View file

@ -1,106 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 10000;
deltaT 1;
writeControl timeStep;
writeInterval 50;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
graphFormat raw;
runTimeModifiable yes;
functions
(
flux
{
type divFlux;
phiName phi;
functionObjectLibs ("libutilityFunctionObjects.so");
}
Mach
{
type MachNumber;
UName U;
functionObjectLibs ("libutilityFunctionObjects.so");
}
minMaxU
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name U;
}
minMaxP
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name p;
}
minMaxRho
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name rho;
}
minMaxT
{
type minMaxField;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldFunctionObjects.so");
name T;
}
);
// ************************************************************************* //

View file

@ -1,83 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default none;
// ddt(rho,U) steadyInertial phi rho 1;
ddt(rho,U) steadyState;
// ddt(rho,h) steadyState;
ddt(rho,h) steadyInertial phi rho 1;
ddt(psi,p) steadyInertial phi rho 1;
U steadyState;
T steadyState;
p steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
// div(phi,U) Gauss linearUpwind faceLimited Gauss linear 1.0;
// div(phi,U) Gauss blended 0.9;
// div(phi,U) Gauss vanLeerDC;
// div(phi,U) Gauss SuperBeeDC;
// div(phi,U) Gauss GammaVDC 0.5;
// div(phi,h) Gauss upwind;
// div(phi,h) Gauss GammaDC 0.5;
div(phi,h) Gauss vanLeerDC;
// div(phid,p) Gauss upwind;
div(phid,p) Gauss vanLeerDC;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div((muEff*dev2(grad(U).T()))) Gauss linear;
div(U,p) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
}
// ************************************************************************* //

View file

@ -9,44 +9,48 @@ FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
5
(
dimensions [0 0 0 1 0 0 0];
internalField uniform 300;
boundaryField
{
INLE1
{
type patch;
nFaces 160;
startFace 152960;
type isentropicTotalTemperature;
phi phi;
rho none;
psi psi;
p p;
T0 uniform 327.335857;
value uniform 300;
}
PRES2
{
type patch;
nFaces 160;
startFace 153120;
type zeroGradient;
}
WALL3
{
type wall;
nFaces 480;
startFace 153280;
type zeroGradient;
}
WALL4
{
type wall;
nFaces 480;
startFace 153760;
type zeroGradient;
}
defaultFaces
{
type empty;
nFaces 153600;
startFace 154240;
}
)
}
// ************************************************************************* //

View file

@ -16,32 +16,20 @@ FoamFile
dimensions [0 1 -1 0 0 0 0];
internalField uniform (234 0 0);
internalField uniform (234.636777 0 0);
boundaryField
{
INLE1
{
type fixedValue;
value uniform (234 0 0);
value uniform (234.636777 0 0);
}
PRES2
{
type inletOutlet;
inletValue uniform (0 0 0);
// type waveTransmissive;
// phi phi;
// rho rho;
// psi psi;
// U U;
// gamma 1.4;
// inletOutlet true;
// correctSupercritical false;
// lInf 0;
// fieldInf (0 0 0);
// value $internalField;
}
WALL3
@ -50,21 +38,17 @@ boundaryField
value uniform (0 0 0);
}
SYMP5
{
type empty;
}
SYMP6
{
type empty;
}
WALL4
{
type fixedValue;
value uniform (0 0 0);
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //

View file

@ -22,37 +22,16 @@ boundaryField
{
INLE1
{
type zeroGradient;
// type waveTransmissiveInlet;
// phi phi;
// rho rho;
// psi psi;
// U U;
// gamma 1.4;
// inletOutlet false;
// correctSupercritical false;
// lInf 0.0;
// fieldInf 101325;
// value uniform 101325;
type isentropicTotalPressure;
U U;
p0 uniform 137491.986;
value uniform 101325;
}
PRES2
{
type fixedValue;
value uniform 101325;
// type waveTransmissive;
// phi phi;
// rho rho;
// psi psi;
// U U;
// gamma 1.4;
// inletOutlet false;
// correctSupercritical true;
// lInf 0.1;
// fieldInf 101325;
// value uniform 101325;
}
WALL3
@ -60,20 +39,16 @@ boundaryField
type zeroGradient;
}
SYMP5
{
type empty;
}
SYMP6
{
type empty;
}
WALL4
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //

View file

@ -23,27 +23,27 @@ vertices
( 0.5 0 -0.1)
( 1.5 0 -0.1)
(-1.5 1 -0.1)
(-0.5 1 -0.1)
( 0.5 1 -0.1)
( 1.5 1 -0.1)
(-1.5 1.5 -0.1)
(-0.5 1.5 -0.1)
( 0.5 1.5 -0.1)
( 1.5 1.5 -0.1)
(-1.5 0 0.1)
(-0.5 0 0.1)
( 0.5 0 0.1)
( 1.5 0 0.1)
(-1.5 1 0.1)
(-0.5 1 0.1)
( 0.5 1 0.1)
( 1.5 1 0.1)
(-1.5 1.5 0.1)
(-0.5 1.5 0.1)
( 0.5 1.5 0.1)
( 1.5 1.5 0.1)
);
blocks
(
hex (0 1 5 4 8 9 13 12) (110 160 1) simpleGrading (1 1 1)
hex (1 2 6 5 9 10 14 13) (110 160 1) simpleGrading (1 1 1)
hex (2 3 7 6 10 11 15 14) (110 160 1) simpleGrading (1 1 1)
hex (0 1 5 4 8 9 13 12) (50 75 1) simpleGrading (1 1 1)
hex (1 2 6 5 9 10 14 13) (50 75 1) simpleGrading (1 1 1)
hex (2 3 7 6 10 11 15 14) (50 75 1) simpleGrading (1 1 1)
);
edges

View file

@ -20,32 +20,32 @@ FoamFile
INLE1
{
type patch;
nFaces 160;
startFace 105110;
nFaces 75;
startFace 22275;
}
PRES2
{
type patch;
nFaces 160;
startFace 105270;
nFaces 75;
startFace 22350;
}
WALL3
{
type wall;
nFaces 330;
startFace 105430;
nFaces 150;
startFace 22425;
}
WALL4
{
type wall;
nFaces 330;
startFace 105760;
nFaces 150;
startFace 22575;
}
defaultFaces
{
type empty;
nFaces 105600;
startFace 106090;
nFaces 22500;
startFace 22725;
}
)

View file

@ -20,7 +20,7 @@ startTime 0;
stopAt endTime;
endTime 2000;
endTime 3000;
deltaT 1;
@ -34,7 +34,7 @@ writeFormat ascii;
writePrecision 6;
writeCompression uncompressed;
writeCompression compressed;
timeFormat general;

View file

@ -18,11 +18,9 @@ ddtSchemes
{
default none;
// ddt(rho,U) steadyInertial phi rho 1;
ddt(rho,U) steadyState;
// ddt(rho,h) steadyState;
ddt(rho,h) steadyInertial phi rho 1;
ddt(psi,p) steadyInertial phi rho 1;
ddt(rho,h) steadyState;
ddt(psi,p) steadyInertial phi rho 0.5;
U steadyState;
T steadyState;
@ -40,23 +38,21 @@ divSchemes
div(phi,U) Gauss upwind;
// div(phi,U) Gauss linearUpwind faceLimited Gauss linear 1.0;
// div(phi,U) Gauss blended 0.9;
// div(phi,U) Gauss vanLeerDC;
// div(phi,U) Gauss SuperBeeDC;
// div(phi,U) Gauss GammaVDC 0.5;
// div(phi,h) Gauss upwind;
// div(phi,h) Gauss GammaDC 0.5;
div(phi,h) Gauss vanLeerDC;
div(phi,h) Gauss upwind;
// div(phi,h) Gauss vanLeerDC;
// div(phid,p) Gauss upwind;
div(phid,p) Gauss vanLeerDC;
div(phid,p) Gauss upwind;
// div(phid,p) Gauss vanLeerDC;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div((muEff*dev2(grad(U).T()))) Gauss linear;
div(U,p) Gauss linear;
div(U) Gauss linear;
}
laplacianSchemes

View file

@ -16,6 +16,20 @@ FoamFile
solvers
{
// Coupled
Up
{
solver BiCGStab;
preconditioner Cholesky;
tolerance 1e-09;
relTol 0.0;
minIter 1;
maxIter 500;
}
// Segregated
p
{
solver BiCGStab;
@ -25,7 +39,7 @@ solvers
maxIter 1000;
tolerance 1e-8;
relTol 0.0;
};
}
U
{
solver BiCGStab;
@ -35,7 +49,8 @@ solvers
maxIter 1000;
tolerance 1e-8;
relTol 0.0;
};
}
h
{
solver BiCGStab;
@ -45,7 +60,7 @@ solvers
maxIter 1000;
tolerance 1e-8;
relTol 0.0;
};
}
}
PIMPLE
@ -58,19 +73,19 @@ PIMPLE
relaxationFactors
{
// Note: under-relaxation factors used in wave-transmissive schemes
U 0.4;
p 0.2;
h 0.5;
U 0.1;
p 0.25;
h 0.1;
rho 0.5;
T 1;
T 0.5;
}
fieldBounds
{
// With bounding
p 50 1e6;
T 20 3000;
U 1000;
p 5e4 3e5;
T 230 500;
U 500;
}
// ************************************************************************* //