Flame propagation (2D)

Update: January 29, 2018
OpenFOAM 4.x

Case directory

$FOAM_TUTORIALS/combustion/XiFoam/ras/moriyoshiHomogeneous

Summary

We calculate the flame when a premixed reactant is ignited in two-dimensional space. We use the b-Ξ two-equation model.

Model geometry Model geometry

Symmetry conditions have been set to the 4 boundaries, left, top, right, and bottom, of the computational domain. Also, in the "ignitionSites" of file "constant/combustionProperties, we set the ignition position to the lower left point of the region (0, 0, 0.0005) as shown below.

ignite          yes;

ignitionSites
(
    {
        location        (0 0 0.0005);
        diameter        0.003;
        start           0;
        duration        0.003;
        strength        2;
    }
);

ignitionSphereFraction 1;
ignitionThickness ignitionThickness [0 1 0 0 0 0 0] 0.001;
ignitionCircleFraction 0.5;
ignitionKernelArea ignitionKernelArea [0 2 0 0 0 0 0] 0.001;		

The reactants and products of combustion are defined in the "reactants" and "products" of the files "constant/thermophysicalProperties", respectively.

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

Mesh Meshes

The calculation result is as follows.

Temperature at 0.001 s (T) Temperature at 0.001 s (T)
Temperature at 0.002 s (T) Temperature at 0.002 s (T)
Temperature at 0.003 s (T) Temperature at 0.003 s (T)
Temperature at 0.004 s (T) Temperature at 0.004 s (T)
Temperature at 0.005 s (T) Temperature at 0.005 s (T)

Commands

cp -r $FOAM_TUTORIALS/combustion/XiFoam/ras/moriyoshiHomogeneous moriyoshiHomogeneous
cd moriyoshiHomogeneous

blockMesh
XiFoam

paraFoam

Calculation time

29.76 seconds *Single, Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz

Reference