Commit e4f8804a authored by Henry's avatar Henry
Browse files

counterFlowFlame2D: make tutorials consistent

parent c3a425ae
/*--------------------------------*- 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 CO2;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
fuel
{
type fixedValue;
value uniform 0;
}
air
{
type fixedValue;
value uniform 0;
}
outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //
/*--------------------------------*- 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 H2O;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
fuel
{
type fixedValue;
value uniform 0;
}
air
{
type fixedValue;
value uniform 0;
}
outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //
......@@ -17,21 +17,26 @@ FoamFile
dimensions [0 0 0 0 0 0 0];
internalField uniform 0.77;
internalField uniform 1;
boundaryField
{
fuel
{
type calculated;
type fixedValue;
value uniform 0.0;
}
air
{
type calculated;
type fixedValue;
value uniform 0.77;
}
outlet
{
type calculated;
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
frontAndBack
{
......
......@@ -17,7 +17,7 @@ FoamFile
dimensions [0 0 0 0 0 0 0];
internalField uniform 0.23;
internalField uniform 0;
boundaryField
{
......@@ -34,8 +34,8 @@ boundaryField
outlet
{
type inletOutlet;
inletValue uniform 0.23;
value uniform 0.23;
inletValue uniform 0;
value uniform 0;
}
frontAndBack
{
......
......@@ -17,19 +17,19 @@ FoamFile
dimensions [1 -1 -1 0 0 0 0];
internalField uniform 0.0;
internalField uniform 0;
boundaryField
{
fuel
{
type fixedValue;
value uniform 0.0;
value uniform 0;
}
air
{
type fixedValue;
value uniform 0.0;
value uniform 0;
}
outlet
{
......
......@@ -24,12 +24,12 @@ boundaryField
fuel
{
type fixedValue;
value uniform 1.e-5;
value uniform 1e-05;
}
air
{
type fixedValue;
value uniform 1.e-5;
value uniform 1e-05;
}
outlet
{
......
......@@ -17,19 +17,19 @@ FoamFile
dimensions [0 2 -2 0 0 0 0];
internalField uniform 3.75e-5;
internalField uniform 3.75e-05;
boundaryField
{
fuel
{
type fixedValue;
value uniform 3.75e-5;
value uniform 3.75e-05;
}
air
{
type fixedValue;
value uniform 3.75e-5;
value uniform 3.75e-05;
}
outlet
{
......
......@@ -17,19 +17,19 @@ FoamFile
dimensions [1 -1 -1 0 0 0 0];
internalField uniform 0.0;
internalField uniform 0;
boundaryField
{
fuel
{
type fixedValue;
value uniform 0.0;
value uniform 0;
}
air
{
type fixedValue;
value uniform 0.0;
value uniform 0;
}
outlet
{
......
......@@ -15,7 +15,6 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 1e5;
......
/*--------------------------------*- 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 rho;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -3 0 0 0 0 0];
internalField uniform 1.18001;
boundaryField
{
fuel
{
type calculated;
value uniform 0.243121;
}
air
{
type calculated;
value uniform 0.438564;
}
outlet
{
type calculated;
value uniform 1.18001;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //
......@@ -22,7 +22,7 @@ active true;
PaSRCoeffs
{
Cmix 0.1;
turbulentReaction on;
turbulentReaction off;
}
......
......@@ -11,10 +11,10 @@ reactions
{
methaneReaction
{
type irreversibleArrheniusReaction;
reaction "CH4 + 2O2 = CO2 + 2H2O";
A 5.2e16;
beta 0;
Ta 14906;
type irreversibleArrheniusReaction;
reaction "CH4 + 2O2 = CO2 + 2H2O";
A 5.2e16;
beta 0;
Ta 14906;
}
}
......@@ -23,7 +23,6 @@ ddtSchemes
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}
divSchemes
......@@ -42,13 +41,7 @@ divSchemes
laplacianSchemes
{
default Gauss linear uncorrected;
laplacian(muEff,U) Gauss linear uncorrected;
laplacian(mut,U) Gauss linear uncorrected;
laplacian(DkEff,k) Gauss linear uncorrected;
laplacian(DepsilonEff,epsilon) Gauss linear uncorrected;
laplacian((rho*(1|A(U))),p) Gauss linear uncorrected;
laplacian(alphaEff,h) Gauss linear uncorrected;
default Gauss linear orthogonal;
}
interpolationSchemes
......@@ -58,7 +51,7 @@ interpolationSchemes
snGradSchemes
{
default uncorrected;
default orthogonal;
}
fluxRequired
......
......@@ -34,7 +34,6 @@ boundaryField
outlet
{
type inletOutlet;
// inletValue uniform 0.77;
inletValue uniform 1;
value uniform 1;
......
......@@ -34,7 +34,6 @@ boundaryField
outlet
{
type inletOutlet;
// inletValue uniform 0.23;
inletValue uniform 0;
value uniform 0;
}
......
#!/bin/sh
cd ${0%/*} || exit 1 # run from this directory
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
# create mesh
runApplication blockMesh
# run the solver
runApplication `getApplication`
......@@ -22,7 +22,7 @@ active true;
PaSRCoeffs
{
Cmix 0.1;
turbulentReaction on;
turbulentReaction off;
useReactionRate true;
}
......
......@@ -11,7 +11,7 @@ FoamFile
format ascii;
class dictionary;
location "constant";
object thermo.icoPoly8;
object thermo.compressibleGas;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
......@@ -22,20 +22,18 @@ O2
nMoles 1;
molWeight 31.9988;
}
equationOfState
{
rhoCoeffs<8> ( 4.1388 -0.016603 3.3039e-05 -3.6455e-08 2.3302e-11 -8.5664e-15 1.6783e-18 -1.356e-22 );
}
thermodynamics
{
Hf 0;
Sf 6407.93;
CpCoeffs<8> ( 904.2 -0.35136 0.0019339 -2.6639e-06 1.8222e-09 -6.75e-13 1.2962e-16 -1.0131e-20 );
Tlow 200;
Thigh 5000;
Tcommon 1000;
highCpCoeffs ( 3.69758 0.00061352 -1.25884e-07 1.77528e-11 -1.13644e-15 -1233.93 3.18917 );
lowCpCoeffs ( 3.21294 0.00112749 -5.75615e-07 1.31388e-09 -8.76855e-13 -1005.25 6.03474 );
}
transport
{
muCoeffs<8> ( 1.2346e-06 6.5097e-08 -2.8362e-11 8.301e-15 -9.7467e-19 0 0 0 );
kappaCoeffs<8> ( -0.0017889 9.5996e-05 -3.129e-08 7.9407e-12 -8.1996e-16 0 0 0 );
As 1.67212e-06;
Ts 170.672;
}
}
......@@ -46,20 +44,18 @@ H2O
nMoles 1;
molWeight 18.0153;
}
equationOfState
{
rhoCoeffs<8> ( 2.3302 -0.0093477 1.8601e-05 -2.0524e-08 1.3119e-11 -4.8229e-15 9.4489e-19 -7.6342e-23 );
}
thermodynamics
{
Hf -13424600;
Sf 10475.4;
CpCoeffs<8> ( 1570.5 1.6886 -0.0039382 6.2482e-06 -4.9342e-09 2.0513e-12 -4.3361e-16 3.6812e-20 );
Tlow 200;
Thigh 5000;
Tcommon 1000;
highCpCoeffs ( 2.67215 0.00305629 -8.73026e-07 1.201e-10 -6.39162e-15 -29899.2 6.86282 );
lowCpCoeffs ( 3.38684 0.00347498 -6.3547e-06 6.96858e-09 -2.50659e-12 -30208.1 2.59023 );
}
transport
{
muCoeffs<8> ( 1.2346e-06 6.5097e-08 -2.8362e-11 8.301e-15 -9.7467e-19 0 0 0 );
kappaCoeffs<8> ( 0.0069975 0.00013726 8.6579e-09 -7.1514e-12 8.441e-16 0 0 0 );
As 1.67212e-06;
Ts 170.672;
}
}
......@@ -68,22 +64,20 @@ CH4
specie
{
nMoles 1;
molWeight 16.043;
}
equationOfState
{
rhoCoeffs<8> ( 2.0751 -0.0083243 1.6564e-05 -1.8277e-08 1.1683e-11 -4.2949e-15 8.4144e-19 -6.7984e-23 );
molWeight 16.0428;
}
thermodynamics
{
Hf -4.66836e+06;
Sf 11597.7;
CpCoeffs<8> ( 996.92 4.0854 -0.00092999 1.957e-06 -2.9777e-09 1.7459e-12 -4.5594e-16 4.4774e-20 );
Tlow 200;
Thigh 6000;
Tcommon 1000;
highCpCoeffs ( 1.63543 0.0100844 -3.36924e-06 5.34973e-10 -3.15528e-14 -10005.6 9.9937 );
lowCpCoeffs ( 5.14988 -0.013671 4.91801e-05 -4.84744e-08 1.66694e-11 -10246.6 -4.64132 );
}
transport
{
muCoeffs<8> ( 1.2346e-06 6.5097e-08 -2.8362e-11 8.301e-15 -9.7467e-19 0 0 0 );
kappaCoeffs<8> ( -0.018501 0.00021332 1.4678e-07 -7.7129e-11 1.0936e-14 0 0 0 );
As 1.67212e-06;
Ts 170.672;
}
}
......@@ -94,20 +88,18 @@ CO2
nMoles 1;
molWeight 44.01;
}
equationOfState
{
rhoCoeffs<8> ( 5.6924 -0.022836 4.544e-05 -5.014e-08 3.2048e-11 -1.1782e-14 2.3083e-18 -1.865e-22 );
}
thermodynamics
{
Hf -8.9423e+06;
Sf 4856.6;
CpCoeffs<8> ( 462.23 1.5803 -0.0010515 7.4618e-08 3.216e-10 -1.9646e-13 4.7679e-17 -4.3027e-21 );
Tlow 200;
Thigh 5000;
Tcommon 1000;
highCpCoeffs ( 4.45362 0.00314017 -1.27841e-06 2.394e-10 -1.66903e-14 -48967 -0.955396 );
lowCpCoeffs ( 2.27572 0.00992207 -1.04091e-05 6.86669e-09 -2.11728e-12 -48373.1 10.1885 );
}
transport
{
muCoeffs<8> ( 1.2346e-06 6.5097e-08 -2.8362e-11 8.301e-15 -9.7467e-19 0 0 0 );
kappaCoeffs<8> ( -0.0076857 0.00010773 -2.7043e-08 4.4967e-12 -3.0218e-16 0 0 0 );
As 1.67212e-06;
Ts 170.672;
}
}
......@@ -118,20 +110,18 @@ N2
nMoles 1;
molWeight 28.0134;
}
equationOfState
{
rhoCoeffs<8> ( 3.6234 -0.014535 2.8924e-05 -3.1915e-08 2.0399e-11 -7.4994e-15 1.4693e-18 -1.1871e-22 );
}
thermodynamics
{
Hf 51.2816;
Sf 6836.49;
CpCoeffs<8> ( 1049.8 -0.16537 0.00036105 2.0824e-07 -5.3568e-10 3.2274e-13 -8.3335e-17 8.054e-21 );
Tlow 200;
Thigh 5000;
Tcommon 1000;
highCpCoeffs ( 2.92664 0.00148798 -5.68476e-07 1.0097e-10 -6.75335e-15 -922.798 5.98053 );
lowCpCoeffs ( 3.29868 0.00140824 -3.96322e-06 5.64152e-09 -2.44486e-12 -1020.9 3.95037 );
}
transport
{
muCoeffs<8> ( 1.2346e-06 6.5097e-08 -2.8362e-11 8.301e-15 -9.7467e-19 0 0 0 );
kappaCoeffs<8> ( 0.0028213 8.6668e-05 -1.5148e-08 1.2499e-12 3.1675e-17 0 0 0 );
As 1.67212e-06;
Ts 170.672;
}
}
......
......@@ -19,10 +19,10 @@ thermoType
{
type heRhoThermo;
mixture reactingMixture;
transport polynomial;
thermo hPolynomial;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState icoPolynomial;
equationOfState perfectGas;
specie specie;
}
......@@ -32,7 +32,7 @@ chemistryReader foamChemistryReader;
foamChemistryFile "$FOAM_CASE/constant/reactions";
foamChemistryThermoFile "$FOAM_CASE/constant/thermo.icoPoly8";
foamChemistryThermoFile "$FOAM_CASE/constant/thermo.compressibleGas";
// ************************************************************************* //
......@@ -23,31 +23,26 @@ ddtSchemes
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,Yi_h) Gauss upwind;
div(phi,h) Gauss upwind;
div(phi,K) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,U) Gauss limitedLinearV 1;
div(phi,Yi_h) Gauss limitedLinear01 1;
div(phi,h) Gauss limitedLinear 1;
div(phi,K) Gauss limitedLinear 1;
div(phid,p) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,k) Gauss limitedLinear 1;
div((muEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear orthogonal;
laplacian(muEff,U) Gauss linear orthogonal;
laplacian(mut,U) Gauss linear orthogonal;
laplacian(DkEff,k) Gauss linear orthogonal;
laplacian(DepsilonEff,epsilon) Gauss linear orthogonal;
laplacian((rho*(1|A(U))),p) Gauss linear orthogonal;
laplacian(alphaEff,h) Gauss linear orthogonal;
}
interpolationSchemes
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment