[Mr.Hack]

Hi,

I am simulating a combustion problem with the standard EBU model. Everything works fine and I am almost done. Now I wanted to plot the reaction rate, but I think there is no field function existing (only the chemistry heat release rate, which is "the sum of reaction rate multiplied by formation enthalpy for each species"). I would like to plot the reaction rate for just one species.

Is there anybody who can help me with this issue?

Maybe it would already help me to see the definition of the "chemistry heat release function", but I do not know, if this is possible?

Thank you in advance for your help.

[Mr.Hack]

Ok, I hopefully solved it already by my own. I still need to check if the results from my field function are correct.

(By the way; maybe cd-adapco should revise the user guide...)

Best regards...

[TomLboro]

Hi,

There is a function inbuilt in Star that shows the reaction rate for the reactions that you have set. It is a primitive function but could be a starting point.

It is in the field functions node and would be have the same name as the reaction.

[lx882211]

Originally Posted by Mr.Hack

Ok, I hopefully solved it already by my own. I still need to check if the results from my field function are correct.

(By the way; maybe cd-adapco should revise the user guide...)

Best regards...

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")
Rrate = 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.