[Scot]

I just downloaded OpenFOAM (1.7.1) (SUSE linux) and I got reactingFoam but did not get rhoReactingFoam. Would someone be so kind as to compare and contrast the two? Why use one versus the other? I see reactingFoam uses hCombustionThermo.H and psiChemistryModel.H while rhoReactingFoam uses hReactingThermo.H and rhoChemistryModel.H. What am I missing by not having rhoReactingFoam? Thanks!

[Hìr0]

I'm also interested.....

[linch]

same for me

[Edison_Ge]

same here
And I'm interested in a incompressible combustion solver, any suggestion?

[ni-openfoam-user]

Dear colleagues,

Did you ever get an answer to your question or resolve this issue?

I am attempting to simulate a high pressure hydrogen release from a pipe using ReactingFoam (release only = no chemistry and no reaction) and wish to run it as an incompressible flow (simplification)

Is it possible to modify the ReactingFoam solver from compressible to incompressible?

Many thanks,

James

[adrieno]

anyone able to answer in 2016 ?

[wyldckat]

Quick answers:

1. Reason why the two solvers exist:
   1. Description in the source code for reactingFoam:
      

      Solver for combustion with chemical reactions.

   2. Description in the source code for rhoReactingFoam:
      

      Solver for combustion with chemical reactions using density based thermodynamics package.

2. You can do incompressible flow by configuring the file "constant/thermophysicalProperties" as an incompressible fluid. Run the following command to look for examples:
   

   find $FOAM_TUTORIALS -name thermophysicalProperties

3. As to why an explicit implementation of incompressible combustion solver doesn't exist, it is most likely because the developers of OpenFOAM aim for generic implementations, to minimize maintenance. For example, see this bug report: http://www.openfoam.org/mantisbt/view.php?id=1856

[adrieno]

thank you for your reply wyldckat,
Sorry for my ignorance but how did you see that ? I mean in the source code ?

[lx882211]

Originally Posted by wyldckat

Quick answers:

1. Reason why the two solvers exist:
   1. Description in the source code for reactingFoam:
   2. Description in the source code for rhoReactingFoam:
      

      Solver for combustion with chemical reactions using density based thermodynamics package.

2. You can do incompressible flow by configuring the file "constant/thermophysicalProperties" as an incompressible fluid. Run the following command to look for examples:
   

   find $FOAM_TUTORIALS -name thermophysicalProperties

3. As to why an explicit implementation of incompressible combustion solver doesn't exist, it is most likely because the developers of OpenFOAM aim for generic implementations, to minimize maintenance. For example, see this bug report: http://www.openfoam.org/mantisbt/view.php?id=1856

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.