CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   reactingFoam vs rhoReactingFoam (https://www.cfd-online.com/Forums/openfoam-solving/86464-reactingfoam-vs-rhoreactingfoam.html)

Scot March 23, 2011 10:31

reactingFoam vs rhoReactingFoam
 
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 April 5, 2011 18:08

I'm also interested.....

linch May 20, 2011 08:18

same for me

Edison_Ge May 22, 2011 08:29

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

ni-openfoam-user July 9, 2014 09:21

Incompressible ReactingFoam
 
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 April 5, 2016 09:24

anyone able to answer in 2016 ? :confused:

wyldckat April 16, 2016 14:25

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

      Solver for combustion with chemical reactions.
    2. Description in the source code for rhoReactingFoam:
      Code:

      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:
    Code:

    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 April 18, 2016 04:03

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

lx882211 June 2, 2016 12:28

Quote:

Originally Posted by wyldckat (Post 595352)
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:
      Code:

      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:
    Code:

    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.


All times are GMT -4. The time now is 19:07.