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

Question regarding Energy equation in reactingFoam

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

Like Tree1Likes
  • 1 Post By Cyp

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 3, 2014, 18:36
Default Question regarding Energy equation in reactingFoam
  #1
Member
 
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12
Dan1788 is on a distinguished road
Dear all,

I am using reactingFoam for simulating a 2d axisymmetric jet flame. The following in the default energy equation in reactingFoam:

{
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;
}

What does the term "fvc::ddt(rho, K) + fvc::div(phi, K)" signify in the above, specifically what is the variable K in this equation?

Also what does the term "fvOptions(rho, he)" mean?

I would be really glad if anyone could help me out. Thanks!
Dan1788 is offline   Reply With Quote

Old   July 4, 2014, 17:52
Default
  #2
Member
 
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12
Dan1788 is on a distinguished road
Foamers,

Anyone who can help? would really appreciate it
Dan1788 is offline   Reply With Quote

Old   August 1, 2014, 11:27
Default
  #3
Senior Member
 
dkxls's Avatar
 
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 156
Rep Power: 19
dkxls will become famous soon enough
Hi Daniel,

maybe at first a general advise: Please use the CODE tags, makes it much easier to read your post!
See also here: http://www.cfd-online.com/Forums/ope...tml#post424508

Quote:
Originally Posted by Dan1788 View Post
What does the term "fvc::ddt(rho, K) + fvc::div(phi, K)" signify in the above, specifically what is the variable K in this equation?
K is the kinetic energy; see the files createFields.H and pEqn.H in
Code:
applications/solvers/combustion/reactingFoam/
This "new" formulation of the energy equation got introduced with version 2.1, see OpenFOAM 2.0 for an alternative formulation of the energy equation (or a good text book will also do ).

Quote:
Originally Posted by Dan1788 View Post
Also what does the term "fvOptions(rho, he)" mean?
If you don't know what it is you can forget about it at this point. See here for a more detailed explanation:
http://www.openfoam.org/version2.2.0/fvOptions.php
dkxls is offline   Reply With Quote

Old   August 2, 2014, 16:07
Default
  #4
Member
 
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12
Dan1788 is on a distinguished road
Thanks Armin, this info certainly helped

I had another question regarding the momentum predictor statement in OpenFoam.

In the default version of reactingFoam (the solver I am using), the momentum predictor is turned off. Does this mean that the momentum equation is not solved ?? How is the velocity then computed ??

Please let me know if you have some knowledge about this. Thanks again!
Dan1788 is offline   Reply With Quote

Old   August 3, 2014, 01:46
Default
  #5
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
Hi Daniel,

When your momentum predictor is set to off, indeed the matrix UEqn is not solved. However, you will notice in the source code that the matrix in created and updated with the flux of the previous time step. The variable H and rAU depend on this updated matrix and they are used to solve p.

U is then evaluated during the corrector step, at the end of pEqn file.

Best,
Cyprien
Cyp is offline   Reply With Quote

Old   August 5, 2014, 00:08
Default
  #6
Member
 
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12
Dan1788 is on a distinguished road
Thanks Cyprien,

Can you suggest me some references or books where I can read about the pressure velocity coupling and the nomenclature that OpenFoam uses?
Dan1788 is offline   Reply With Quote

Old   August 5, 2014, 01:22
Default
  #7
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
Sure,

The reference for me remains the book of Patankar "Numerical Heat Transfer And Fluid Flow"

Another good reference, Ferziger and Peric "Computational Method for Fluid Dynamics"

The original paper of PISO is also instersting: Issa "Solution of the implicitly discretised fluid flow equations by operator-splitting"

More specific to OpenFOAM, there is the PhD dissertation of Hrvoje Jasak one of the main contributor to OpenFOAM "Error Analysis and Estimation for the Finite Volume Method with Application to Fluid Flows".

And a personal contribution: http://www.scribd.com/doc/181588911/...ion-of-icoFoam.

Cheers,
Cyprien
Dan1788 likes this.
Cyp 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
Source Term due to evaporation in energy transport equation styleworker OpenFOAM Programming & Development 3 September 7, 2022 04:09
Energy equation for compressible flow majkl OpenFOAM Running, Solving & CFD 5 February 10, 2021 16:17
Energy equation in bouyant*Foam solvers makaveli_lcf OpenFOAM Running, Solving & CFD 0 December 21, 2009 05:09
Viscosity and the Energy Equation Rich Main CFD Forum 0 December 16, 2009 15:01
Energy equation Suresh Balasubramanian FLUENT 1 April 6, 2003 03:54


All times are GMT -4. The time now is 09:47.