[Edoardo1993]

Hello,
I am trying to simulate a flammable gas flow in atmosphere with reactingFoam. You can find the whole case directory attached at the end.

If I do not ignite the gas (/constant/combustionProperties --> active off and /constant/chemistryProperties --> chemistry off) the gas leak is totally diffusive and results are good.
However, if I try to ignite the flammable gas (/constant/combustionProperties --> active on and /constant/chemistryProperties --> chemistry on) the reaction does not occur. Therefore the results are exactly the same as withouth ignition.

The reaction is a single-step reaction, implemented in /constant/reactionGRI.

reactionGRI

elements
(
)
;

species
4
(
H2
O2
H2O
N2
)
;


reactions
{
    un-named-reaction-2
    {
        type            irreversibleinfiniteReaction; 
        reaction        "2H2 + O2 + 3.79N2 = 2H2O + 3.79N2";
    }
}

chemistryProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      chemistryProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

chemistryType
{
    chemistrySolver   ode;
    chemistryThermo   psi;
    TDAC              off;
}

initialChemicalTimeStep 1e-07;

chemistry       on;

odeCoeffs
{
    solver          seulex;  // Rosenbrock34, seulex or rodas23
    absTol          1e-12;
    relTol           1e-7;
}

// ************************************************************************* //

combustionProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      combustionProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

combustionModel  EDC;

active  on;

EDCCoeffs
{
    version v2005;
}

// ************************************************************************* //

thermophysicalProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{

    type            hePsiThermo;
    mixture        reactingMixture; //singleStepReactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

inertSpecie N2;

chemistryReader foamChemistryReader;
foamChemistryFile "<constant>/reactionsGRI";
foamChemistryThermoFile "<constant>/thermo.compressibleGasGRI";

dpdt false;

 // ************************************************************************* //

I have the following question:


1) Why does the reaction not occur?

2) Is irreversibleinfiniteReaction in /constant/reactionGRI correct? Maybe is it not fast enough to show something in a short time?

3) In /constant/thermophysicalProperties should be selected mixture reactingMixture or singleStepReactingMixture?

4) Could be the problem the lack of ignition?




Really thank for your help.