Skip to content
Permalink
Browse files

chtMultiRegionFoam: Improvements to reverseBurner tutorial

The tutorial has been converted from two-dimensions to a wedge and the
flow has been swirl stabilised. The turbulence parameters have been made
physical. The transport schemes have been increased to second order. The
reaction mechanism has been changed to one from a publically accessible
reference. The gas thermodynamics have been made incompressible, and the
pressure offset around zero, which improves the behaviour of the
pressure solution.
  • Loading branch information...
Will Bainbridge
Will Bainbridge committed Jan 24, 2019
1 parent faf6e24 commit a0850290b315ea09e2193e97dab71e7cf0172c34
Showing with 179 additions and 94 deletions.
  1. +0 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/N2
  2. +0 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/O2
  3. +3 −3 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/U
  4. +3 −2 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/epsilon
  5. +3 −2 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/k
  6. +1 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p
  7. +1 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/gas/p_rgh
  8. +2 −2 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/0/solid/p
  9. +2 −0 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/Allrun
  10. +1 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/chemistryProperties
  11. +1 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/combustionProperties
  12. +4 −4 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/fvOptions
  13. +1 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/g
  14. +1 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/radiationProperties
  15. +20 −3 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/reactions
  16. +21 −1 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermo.compressibleGas
  17. +2 −2 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/gas/thermophysicalProperties
  18. +7 −10 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/constant/solid/thermophysicalProperties
  19. +49 −29 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/blockMeshDict
  20. +5 −4 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/controlDict
  21. +7 −7 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSchemes
  22. +6 −17 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/fvSolution
  23. +39 −0 tutorials/heatTransfer/chtMultiRegionFoam/reverseBurner/system/gas/setFieldsDict
@@ -38,7 +38,6 @@ boundaryField
type inletOutlet;
inletValue $internalField;
value $internalField;

}
gas_to_solid
{
@@ -38,7 +38,6 @@ boundaryField
type inletOutlet;
inletValue $internalField;
value $internalField;

}
gas_to_solid
{
@@ -26,12 +26,12 @@ boundaryField
inletFuel
{
type fixedValue;
value uniform (0 0.5 0);
value uniform (0 0.15 0.15);
}
inletAir
{
type fixedValue;
value uniform (0 1.5 0);
value uniform (0 0.15 -0.15);
}
outlet
{
@@ -40,7 +40,7 @@ boundaryField
}
gas_to_solid
{
type noSlip;
type slip;
}
}

@@ -17,15 +17,16 @@ FoamFile

dimensions [ 0 2 -3 0 0 0 0 ];

internalField uniform 0.01;
internalField uniform 1.27351e-06;

boundaryField
{
#includeEtc "caseDicts/setConstraintTypes"

"inlet.*"
{
type fixedValue;
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.1;
value $internalField;
}
outlet
@@ -17,15 +17,16 @@ FoamFile

dimensions [ 0 2 -2 0 0 0 0 ];

internalField uniform 0.1;
internalField uniform 8.4375e-05;

boundaryField
{
#includeEtc "caseDicts/setConstraintTypes"

"inlet.*"
{
type fixedValue;
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value $internalField;
}
outlet
@@ -17,7 +17,7 @@ FoamFile

dimensions [ 1 -1 -2 0 0 0 0 ];

internalField uniform 1e5;
internalField uniform 0;

boundaryField
{
@@ -17,7 +17,7 @@ FoamFile

dimensions [ 1 -1 -2 0 0 0 0 ];

internalField uniform 1e5;
internalField uniform 0;

boundaryField
{
@@ -15,9 +15,9 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];
dimensions [ 1 -1 -2 0 0 0 0 ];

internalField uniform 1e5;
internalField uniform 0;

boundaryField
{
@@ -15,6 +15,8 @@ paraFoam -region solid -touch

runApplication decomposePar -allRegions

runParallel setFields -region gas

runParallel $(getApplication)

runApplication reconstructPar -allRegions -newTimes
@@ -10,7 +10,7 @@ FoamFile
version 2.0;
format ascii;
class dictionary;
location "constant/gas";
location "constant";
object chemistryProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -10,7 +10,7 @@ FoamFile
version 2.0;
format ascii;
class dictionary;
location "constant/gas";
location "constant";
object combustionProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -10,7 +10,7 @@ FoamFile
version 2.0;
format ascii;
class dictionary;
location "constant/gas";
location "constant";
object fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -21,15 +21,15 @@ options
{
type scalarSemiImplicitSource;

timeStart 0.2;
duration 0.3;
timeStart 0;
duration 1;
selectionMode cellSet;
cellSet ignition;
volumeMode specific;

injectionRateSuSp
{
h (2e7 0); // kg/m/s^3
e (5e7 0); // kg/m/s^3
}
}
}
@@ -16,6 +16,6 @@ FoamFile

dimensions [0 1 -2 0 0 0 0];

value (0 -9.81 0);
value (0 0 0);

// ************************************************************************* //
@@ -10,7 +10,7 @@ FoamFile
version 2.0;
format ascii;
class dictionary;
location "constant/gas";
location "constant";
object radiationProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -4,6 +4,23 @@
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
-------------------------------------------------------------------------------
Reference:
Westbrook, C. K., & Dryer, F. L. (1981).
Simplified reaction mechanisms for the oxidation of hydrocarbon fuels in
flames.
Combustion science and technology, 27(1-2), 31-43.
Notes:
The reference gives parameters in a unit system of cm-sec-mole-kcal-Kelvin,
whilst OpenFOAM uses m-sec-kmol-J-Kelvin, both the pre-exponential factors
and the activation temperature/energy require unit conversion.
The mechanism used is set 3 from table 2. This should be the most accurate
single-step methane mechanism without any negative exponents. Set 2 is not
usable because the negative exponent on the methane concentration causes
divergence in regions without any fuel.
\*---------------------------------------------------------------------------*/
FoamFile
{
@@ -29,10 +46,10 @@ reactions
methaneReaction
{
type irreversibleArrheniusReaction;
reaction "CH4 + 2O2 = CO2 + 2H2O";
A 5.2e16;
reaction "CH4^0.2 + 2O2^1.3 = CO2 + 2H2O";
A 2.11873e+11;
beta 0;
Ta 14906;
Ta 24357;
}
}

@@ -10,7 +10,7 @@ FoamFile
version 2.0;
format ascii;
class dictionary;
location "constant/gas";
location "constant";
object thermo.compressibleGas;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -21,6 +21,10 @@ O2
{
molWeight 31.9988;
}
equationOfState
{
pRef 1e6;
}
thermodynamics
{
Tlow 200;
@@ -42,6 +46,10 @@ H2O
{
molWeight 18.0153;
}
equationOfState
{
pRef 1e6;
}
thermodynamics
{
Tlow 200;
@@ -63,6 +71,10 @@ CH4
{
molWeight 16.0428;
}
equationOfState
{
pRef 1e6;
}
thermodynamics
{
Tlow 200;
@@ -84,6 +96,10 @@ CO2
{
molWeight 44.01;
}
equationOfState
{
pRef 1e6;
}
thermodynamics
{
Tlow 200;
@@ -105,6 +121,10 @@ N2
{
molWeight 28.0134;
}
equationOfState
{
pRef 1e6;
}
thermodynamics
{
Tlow 200;
@@ -21,9 +21,9 @@ thermoType
mixture reactingMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
equationOfState incompressiblePerfectGas;
specie specie;
energy sensibleEnthalpy;
energy sensibleInternalEnergy;
}

inertSpecie N2;
@@ -29,23 +29,20 @@ mixture
{
specie
{
molWeight 50;
molWeight 50;
}

transport
equationOfState
{
kappa 80;
rho 8000;
}

thermodynamics
{
Hf 0;
Cp 450;
Hf 0;
Cp 450;
}

equationOfState
transport
{
rho 8000;
kappa 80;
}
}

0 comments on commit a085029

Please sign in to comment.
You can’t perform that action at this time.