CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

How to add Arrhenius Reaction Rate to energy equation?

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By Zack

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 3, 2017, 14:10
Default How to add Arrhenius Reaction Rate to energy equation?
  #1
New Member
 
Join Date: Sep 2012
Posts: 23
Rep Power: 10
Zack is on a distinguished road
Hi All

I want to modify the below EEqn.H in order to add Arrhenius Reaction Rate to it. I am wondering whether I need to call it or add its equation separately to the file?

Current EEqn:
Code:
{
    volScalarField& he = thermo.he();

    fvScalarMatrix EEqn
    (
        fvm::ddt(rho, he) + mvConvection->fvmDiv(phi, he)
      + fvc::ddt(rho, K) + fvc::div(phi, K)
      + (
            he.name() == "e"
          ? fvc::div
            (
                fvc::absolute(phi/fvc::interpolate(rho), U),
                p,
                "div(phiv,p)"
            )
          : -dpdt
        )
      - fvm::laplacian(turbulence->alphaEff(), he)
//    - fvm::laplacian(turbulence->muEff(), he)  // unit lewis no.
     ==
        reaction->Sh()
      + fvOptions(rho, he)
    );

    EEqn.relax();

    fvOptions.constrain(EEqn);

    EEqn.solve();

    fvOptions.correct(he);

    thermo.correct();

    Info<< "min/max(T) = "
        << min(T).value() << ", " << max(T).value() << endl;
}
My energy equation should be like this:

Code:
 ∂(ρE)/∂t + ∂[u(ρE + p)]/∂x + ∂[v(ρE + p)]/∂y = 

∂(uσxx + vσxy − qx)/∂x + ∂(vσyy + uσyx − qy)/∂y + qω
where:

ω == Arrhenius Reaction Rate
Zack is offline   Reply With Quote

Old   January 3, 2017, 14:14
Default
  #2
New Member
 
Join Date: Sep 2012
Posts: 23
Rep Power: 10
Zack is on a distinguished road
Arrhenius Reaction Rate is available at :

src/thermophysicalModels/specie/reaction/reactionRate/ArrheniusReactionRate/
khaledhmz likes this.
Zack is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Porous Modeling of Energy equation in OpenFOAM mohammad_kordo OpenFOAM Running, Solving & CFD 9 November 22, 2020 07:18
Setting the height of the stream in the free channel kevinmccartin CFX 10 July 9, 2015 21:36
error message cuteapathy CFX 14 March 20, 2012 06:45
Segmentation fault in running alternateSteadyReactingFoam,why? NewKid OpenFOAM 18 January 20, 2011 16:55
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 15:45


All times are GMT -4. The time now is 23:56.