From e4f8804a372638ab81fdc7ca975398326574d816 Mon Sep 17 00:00:00 2001 From: Henry Date: Thu, 22 Aug 2013 16:53:02 +0100 Subject: [PATCH] counterFlowFlame2D: make tutorials consistent --- .../reactingFoam/ras/counterFlowFlame2D/0/CO2 | 48 ++++++ .../reactingFoam/ras/counterFlowFlame2D/0/H2O | 48 ++++++ .../reactingFoam/ras/counterFlowFlame2D/0/N2 | 13 +- .../reactingFoam/ras/counterFlowFlame2D/0/O2 | 6 +- .../ras/counterFlowFlame2D/0/alphat | 6 +- .../ras/counterFlowFlame2D/0/epsilon | 4 +- .../reactingFoam/ras/counterFlowFlame2D/0/k | 6 +- .../reactingFoam/ras/counterFlowFlame2D/0/mut | 6 +- .../reactingFoam/ras/counterFlowFlame2D/0/p | 1 - .../reactingFoam/ras/counterFlowFlame2D/0/rho | 46 ++++++ .../constant/combustionProperties | 2 +- .../ras/counterFlowFlame2D/constant/reactions | 10 +- .../ras/counterFlowFlame2D/system/fvSchemes | 11 +- .../counterFlowFlame2D/0/N2 | 1 - .../counterFlowFlame2D/0/O2 | 1 - .../counterFlowFlame2D/Allrun | 11 -- .../constant/combustionProperties | 2 +- .../constant/thermo.compressibleGas | 129 ++++++++++++++++ .../constant/thermo.icoPoly8 | 139 ------------------ .../constant/thermophysicalProperties | 8 +- .../counterFlowFlame2D/system/fvSchemes | 21 +-- .../counterFlowFlame2D/system/fvSolution | 4 +- 22 files changed, 317 insertions(+), 206 deletions(-) create mode 100644 tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/CO2 create mode 100644 tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/H2O create mode 100644 tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/rho delete mode 100755 tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun create mode 100644 tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.compressibleGas delete mode 100644 tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/CO2 b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/CO2 new file mode 100644 index 0000000000..2fef379ec7 --- /dev/null +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/CO2 @@ -0,0 +1,48 @@ +/*--------------------------------*- 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/H2O b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/H2O new file mode 100644 index 0000000000..ca22735630 --- /dev/null +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/H2O @@ -0,0 +1,48 @@ +/*--------------------------------*- 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/N2 b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/N2 index ac188ad505..8b2b2e8b5d 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/N2 +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/N2 @@ -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 { diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/O2 b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/O2 index d7ab36547c..95649fdac4 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/O2 +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/O2 @@ -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 { diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/alphat b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/alphat index f2b4a35b53..c74018bbdb 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/alphat +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/alphat @@ -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 { diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/epsilon b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/epsilon index e536911162..ad65f7ea0f 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/epsilon +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/epsilon @@ -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 { diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/k b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/k index ea9dd15921..3dcedaa2da 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/k +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/k @@ -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 { diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/mut b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/mut index 1cf2a1efca..47d65d567d 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/mut +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/mut @@ -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 { diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/p b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/p index 52393be19b..79f071ce31 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/p +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/p @@ -15,7 +15,6 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - dimensions [1 -1 -2 0 0 0 0]; internalField uniform 1e5; diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/rho b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/rho new file mode 100644 index 0000000000..d1620aaeb0 --- /dev/null +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/0/rho @@ -0,0 +1,46 @@ +/*--------------------------------*- 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/combustionProperties b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/combustionProperties index 6c1f0e9cf3..9c283d5933 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/combustionProperties +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/combustionProperties @@ -22,7 +22,7 @@ active true; PaSRCoeffs { Cmix 0.1; - turbulentReaction on; + turbulentReaction off; } diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/reactions b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/reactions index 839b5e4e45..e548bae90f 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/reactions +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/constant/reactions @@ -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; } } diff --git a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/system/fvSchemes b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/system/fvSchemes index 6d4a329ff2..3731c816a3 100644 --- a/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/system/fvSchemes +++ b/tutorials/combustion/reactingFoam/ras/counterFlowFlame2D/system/fvSchemes @@ -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 diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/N2 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/N2 index 3c48882346..8b2b2e8b5d 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/N2 +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/N2 @@ -34,7 +34,6 @@ boundaryField outlet { type inletOutlet; -// inletValue uniform 0.77; inletValue uniform 1; value uniform 1; diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/O2 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/O2 index 29b38fbb21..95649fdac4 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/O2 +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/O2 @@ -34,7 +34,6 @@ boundaryField outlet { type inletOutlet; -// inletValue uniform 0.23; inletValue uniform 0; value uniform 0; } diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun deleted file mode 100755 index b729092026..0000000000 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun +++ /dev/null @@ -1,11 +0,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` diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties index 754f22dad2..06715e90f1 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/combustionProperties @@ -22,7 +22,7 @@ active true; PaSRCoeffs { Cmix 0.1; - turbulentReaction on; + turbulentReaction off; useReactionRate true; } diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.compressibleGas b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.compressibleGas new file mode 100644 index 0000000000..ff00e560b4 --- /dev/null +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.compressibleGas @@ -0,0 +1,129 @@ +/*--------------------------------*- 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 dictionary; + location "constant"; + object thermo.compressibleGas; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +O2 +{ + specie + { + nMoles 1; + molWeight 31.9988; + } + thermodynamics + { + 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 + { + As 1.67212e-06; + Ts 170.672; + } +} + +H2O +{ + specie + { + nMoles 1; + molWeight 18.0153; + } + thermodynamics + { + 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 + { + As 1.67212e-06; + Ts 170.672; + } +} + +CH4 +{ + specie + { + nMoles 1; + molWeight 16.0428; + } + thermodynamics + { + 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 + { + As 1.67212e-06; + Ts 170.672; + } +} + +CO2 +{ + specie + { + nMoles 1; + molWeight 44.01; + } + thermodynamics + { + 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 + { + As 1.67212e-06; + Ts 170.672; + } +} + +N2 +{ + specie + { + nMoles 1; + molWeight 28.0134; + } + thermodynamics + { + 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 + { + As 1.67212e-06; + Ts 170.672; + } +} + + +// ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 deleted file mode 100644 index b25328716a..0000000000 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 +++ /dev/null @@ -1,139 +0,0 @@ -/*--------------------------------*- 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 dictionary; - location "constant"; - object thermo.icoPoly8; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -O2 -{ - specie - { - 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 ); - } - 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 ); - } -} - -H2O -{ - specie - { - 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 ); - } - 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 ); - } -} - -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 ); - } - 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 ); - } - 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 ); - } -} - -CO2 -{ - specie - { - 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 ); - } - 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 ); - } -} - -N2 -{ - specie - { - 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 ); - } - 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 ); - } -} - - -// ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties index 8b36a40f5d..4c54af61d5 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties @@ -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"; // ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes index 58be93e54d..dd00d8dd67 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes @@ -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 diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution index f1e4568f0e..552f948a84 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution @@ -86,9 +86,9 @@ PIMPLE rhoMin 0.1; rhoMax 1.5; - maxCo 0.5; + maxCo 0.1; rDeltaTSmoothingCoeff 0.1; - maxDeltaT 1; + maxDeltaT 2e-4; alphaTemp 0.005; }