Combustion of methane

Update: January 24, 2018
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/combustion/fireFoam/les/smallPoolFire3D

Summary

Methane flows in and burns from the region "inlet", which is a part of the floor of the air-filled space. The regions other than the floor are assumed to be open.

Model geometry Model geometry

The chemical species and combustion reactions with methane and air are defined in the file "reactions" in the directory "constant" as follows. The properties of each chemical species are defined in the file "thermo.compressibleGas" in the same directory.

species
(
    O2
    H2O
    CH4
    CO2
    N2
);
reactions
{
    propaneReaction
    {
        type         irreversibleinfiniteReaction;
        reaction      "CH4 + 2O2 + 7.5N2 = CO2 + 2H2O + 7.5N2";
    }
}

These two file "reactions", "thermo.compressibleGas", are imported in the file "thermophysicalProperties" in the directory "constant" to read the settings.

The meshes are as follows, and the number of mesh is 216000.

Mesh Meshes

The calculation result is as follows. The figures are rendered with volume rendering.

CO2 at 0.1 sec (CO2) CO2 at 0.1 sec (CO2)
CO2 at 0.5 sec (CO2) CO2 at 0.5 sec (CO2)
CO2 at 0.8 sec (CO2) CO2 at 0.8 sec (CO2)
CO2 at 1 sec (CO2) CO2 at 1 sec (CO2)
CO2 at 2 sec (CO2) CO2 at 2 sec (CO2)
CO2 at 4 sec (CO2) CO2 at 4 sec (CO2)

Commands

cp -r $FOAM_TUTORIALS/combustion/fireFoam/les/oppositeBurningPanels oppositeBurningPanels
cd oppositeBurningPanels

blockMesh
topoSet
createPatch -overwrite

cp 0/ph_rgh.orig 0/ph_rgh

decomposePar -force
mpirun -np 4 fireFoam -parallel
reconstructPar

paraFoam

Calculation time

1 hour 27 minutes 49.01 seconds *4 parallel, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz