|
[Sponsors] |
Question regarding Energy equation in reactingFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 3, 2014, 18:36 |
Question regarding Energy equation in reactingFoam
|
#1 |
Member
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12 |
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! |
|
July 4, 2014, 17:52 |
|
#2 |
Member
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12 |
Foamers,
Anyone who can help? would really appreciate it |
|
August 1, 2014, 11:27 |
|
#3 | |
Senior Member
Armin
Join Date: Feb 2011
Location: Helsinki, Finland
Posts: 156
Rep Power: 19 |
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:
Code:
applications/solvers/combustion/reactingFoam/ 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 |
||
August 2, 2014, 16:07 |
|
#4 |
Member
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12 |
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! |
|
August 3, 2014, 01:46 |
|
#5 |
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18 |
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 |
|
August 5, 2014, 00:08 |
|
#6 |
Member
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 12 |
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? |
|
August 5, 2014, 01:22 |
|
#7 |
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |