[harshad88]

Hello all,

I am still a newbie with OpenFoam. I want to know what does the boolean statement useReactionRate TRUE achieve when using the reactingFoam solver ??

Any help would be greatly appreciated. Thanks!

[harshad88]

Anyone who can help?

[alexeym]

Hi,

Assuming that you mean useReactionRate in PaSRCoeffs dictionary in constant/combustionProperties file (though it was non-trivial to figure this out), it has no direct connection with reactingFoam but it's a property of combustoin model.

If you go to src/combustionModels/PaSR and look as PaSR.C, you'll find there:

template<class Type>
void Foam::combustionModels::PaSR<Type>::correct()
{
    if (this->active())
    {
        ...
        if (!useReactionRate_)
        {
            this->chemistryPtr_->solve(t - dt, dt);
        }
        else
        {
            this->chemistryPtr_->calculate();
        }
...
}

there type is defined by the line "combustionModel PaSR<psiChemistryCombustion>" in combustionProperties. Properties of psiChemistryCombustion are defined in chemistryProperties dictionary. If you then go to src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel and open basicChemistryModel.H file, you'll find the following code and comments:

            // Chemistry solution
                
                //- Calculates the reaction rates
                virtual void calculate() = 0;

                //- Solve the reaction system for the given start time and
                //  timestep and return the characteristic time
                virtual scalar solve(const scalar t0, const scalar deltaT) = 0;

Now you can use Doxygen generated documentation and sources to find concrete implementations of these methods.

[Howard]

Originally Posted by alexeym

Hi,

Assuming that you mean useReactionRate in PaSRCoeffs dictionary in constant/combustionProperties file (though it was non-trivial to figure this out), it has no direct connection with reactingFoam but it's a property of combustoin model.

If you go to src/combustionModels/PaSR and look as PaSR.C, you'll find there:

template<class Type>
void Foam::combustionModels::PaSR<Type>::correct()
{
    if (this->active())
    {
        ...
        if (!useReactionRate_)
        {
            this->chemistryPtr_->solve(t - dt, dt);
        }
        else
        {
            this->chemistryPtr_->calculate();
        }
...
}

there type is defined by the line "combustionModel PaSR<psiChemistryCombustion>" in combustionProperties. Properties of psiChemistryCombustion are defined in chemistryProperties dictionary. If you then go to src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel and open basicChemistryModel.H file, you'll find the following code and comments:

            // Chemistry solution
 
                //- Calculates the reaction rates
                virtual void calculate() = 0;
 
                //- Solve the reaction system for the given start time and
                //  timestep and return the characteristic time
                virtual scalar solve(const scalar t0, const scalar deltaT) = 0;

Now you can use Doxygen generated documentation and sources to find concrete implementations of these methods.

Thank you very much for your answer. I am new user of OpenFoam and I also try to understand the parametersn and the model uses in the reactiongFoam. Could you explain how to see the detail thing by 'Doxygen generated documentation and sources', I am not quite clear.

[amin_jalalian]

Hi
the "useReactionRate" has been removed in OF301. It is not in src's files or other folder's, but it is in some tutorial files which is not true and their programmers forgot to remove that phrase from tutorials.
regards.

[amin_jalalian]

Originally Posted by amin_jalalian

Hi
the "useReactionRate" has been removed in OF301. It is not in src's files or other folder's, but it is in some tutorial files which is not true and their programmers forgot to remove that phrase from tutorials.
regards.

better:

Hi
the "useReactionRate" in has been removed in OF301. It is not in src's files or other folder's, but it is still in 3 tutorial files (file: combustionProperties) which is not right and it seems its programmers forgot to remove "useReactionRate" from tutorials.
regards.

[lx882211]

Originally Posted by alexeym

Hi,

Assuming that you mean useReactionRate in PaSRCoeffs dictionary in constant/combustionProperties file (though it was non-trivial to figure this out), it has no direct connection with reactingFoam but it's a property of combustoin model.

If you go to src/combustionModels/PaSR and look as PaSR.C, you'll find there:

template<class Type>
void Foam::combustionModels::PaSR<Type>::correct()
{
    if (this->active())
    {
        ...
        if (!useReactionRate_)
        {
            this->chemistryPtr_->solve(t - dt, dt);
        }
        else
        {
            this->chemistryPtr_->calculate();
        }
...
}

there type is defined by the line "combustionModel PaSR<psiChemistryCombustion>" in combustionProperties. Properties of psiChemistryCombustion are defined in chemistryProperties dictionary. If you then go to src/thermophysicalModels/chemistryModel/chemistryModel/basicChemistryModel and open basicChemistryModel.H file, you'll find the following code and comments:

            // Chemistry solution
                
                //- Calculates the reaction rates
                virtual void calculate() = 0;

                //- Solve the reaction system for the given start time and
                //  timestep and return the characteristic time
                virtual scalar solve(const scalar t0, const scalar deltaT) = 0;

Now you can use Doxygen generated documentation and sources to find concrete implementations of these methods.

I meet a problem in Openfoam 2.3.0, reactingFoam, i need the data of CH4 reaction rate, but i failed to autowrite it like U or Yi.

First i added this code in the CreatField.H,

volScalarField Rrate
(
IOobject
(
"Rrate",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("Rrate", dimMass/dimVolume/dimTime, 0.0)
);

Then, i added this code in YEqn,

forAll(Y, i)
{
if (Y[i].name() != "CH4")
RR = reaction->R(Yi);
}

After wmake, it shows

YEqn.H:26:14: error: no match for ‘operator=’ (operand types are ‘Foam::volScalarField {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>}’ and ‘Foam::tmp<Foam::fvMatrix<double> >’)
Rrate = reaction->R(Yi);

Can you tell me how to do that? Thank you very much.