[bany]

When I study the standardChemistryModel.C, the function which calculates dc/dt called omega is different from that in openfoam v7. There is an additional label variable "li" in v9 compared to that in v7.

in openfoam v9:

        virtual void omega
        (
            const scalar p,
            const scalar T,
            const scalarField& c,
            const label li,
            scalarField& dcdt
        ) const;

in openfoam v7:

        virtual void omega
        (
            const scalarField& c,
            const scalar T,
            const scalar p,
            scalarField& dcdt
        ) const;

However, when I went deep into the function omega, the label variable of the function omega in Reaction.C is used in code :

    omegaf = this->kf(p, clippedT, c, li);
    omegar = this->kr(omegaf, p, clippedT, c, li);

For reversible reactions, kf if defined in ReversibleReaction.C:

template<class ReactionThermo, class ReactionRate>
Foam::scalar Foam::ReversibleReaction<ReactionThermo, ReactionRate>::kf
(
    const scalar p,
    const scalar T,
    const scalarField& c,
    const label li
) const
{
    return k_(p, T, c, li);
}

For Arrhenius reaction rate, above k_ can invoke function operator() in ArrheniusReactionRateI.H:

inline Foam::scalar Foam::ArrheniusReactionRate::operator()
(
    const scalar p,
    const scalar T,
    const scalarField&,
    const label
) const
{
    scalar ak = A_;

    if (mag(beta_) > vSmall)
    {
        ak *= pow(T, beta_);
    }

    if (mag(Ta_) > vSmall)
    {
        ak *= exp(-Ta_/T);
    }

    return ak;
}

At last, i could not find where is the label variable li used.

What dose the meaning of the existence of the useless label variable?

Any help are appreciated!