
[Sponsors] 
September 20, 2010, 10:50 
Viscous terms in rhoCentralFoam

#1  
Member
Florian Ettner
Join Date: Mar 2009
Location: Munich, Germany
Posts: 40
Rep Power: 8 
Dear Foamers,
I have a few questions on the viscous terms in rhoCentralFoam. In the source code, the energy equation looks like this: Quote:
2. I think I understand the heat conduction in line 200. What kind of conduction of internal energy is in lines 198199 ? 3. What is the use of the repeated use of fvm and fvc discretisation of the same expressions in lines 197199? Why is it not applied to the heat conduction in line 200? I'd be glad if anybody could help me understanding this issue. It would be quite important to me, as I currently work on extending the solver to reactive flows. Thank you in advance! 

March 14, 2012, 05:47 

#2 
Member
Alexander
Join Date: Mar 2009
Posts: 39
Rep Power: 8 
Good day to everyone!
I also use rhoCentralFoam and I am interesting in questions about viscous dissipation. And about your questions: 2. "I think I understand the heat conduction in line 200. What kind of conduction of internal energy is in lines 198199 ?" I want to know answer too =) Probably, it is turbulent heat transfer. In new version of rhoCentralFoam it is written as Code:
 fvm::laplacian(turbulence>alphaEff(), e) + fvc::laplacian(turbulence>alpha(), e) Answer is simple: there is solving of "corrector" step, because "fvm" means that for this term the matrix is prepared and the system wil be solved ("m" for "matrix"). Term with "fvc" means that term will be calculated in explicit way (the right part of matrix equation; "c" for "calculate"). So, lines 197200 is the corrector which include the influence of thermal conductivity part in energy equation. Why turbulence is there? I do not know. I hope, alpha is zero for laminar flow =) But the main question: is the dissipation correctly taken into account? Because in my laminar transsonic tests (with low Reynolds number about 100) heating due to dissipation is absent (although it should be meaningful) =( 

March 21, 2012, 04:05 

#3 
Member
Alexander
Join Date: Mar 2009
Posts: 39
Rep Power: 8 
I've just found that Prandtl number in rhoCentralFoam is always 1 (independent of what you have set in constant/thermophysicalProperties).
See bug http://www.openfoam.org/mantisbt/view.php?id=475 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Implicit treatment of advection terms and pressure correction  nikosb  Main CFD Forum  0  January 17, 2010 17:07 
Question in definition of terms in solve  titio  OpenFOAM Running, Solving & CFD  0  March 19, 2009 17:02 
Limited Gradients for Viscous fluxes and Sources  andy  Main CFD Forum  7  June 16, 2006 11:38 
Implicitisation of viscous terms  Michael HallsMoore  Main CFD Forum  0  April 3, 2006 07:52 
KEpsilon model?  Brindaban Ghosh  Main CFD Forum  2  June 24, 2000 04:22 