CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Error in heat conduction rate calculation in OpenFOAM while using rhoCentralFOAM

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

Like Tree1Likes
  • 1 Post By Daedalus

LinkBack Thread Tools Search this Thread Display Modes
Old   January 29, 2019, 06:53
Default Error in heat conduction rate calculation in OpenFOAM while using rhoCentralFOAM
New Member
Gaurav Kumar
Join Date: Jul 2015
Location: India
Posts: 1
Rep Power: 0
Daedalus is on a distinguished road

This chunk of code represents viscous correction for energy equation in rhoCentralFoam solver in OpenFOAM-dev.
if (!inviscid)

fvm::ddt(rho, e) - fvc::ddt(rho, e)
  - fvm::laplacian(turbulence->alphaEff(), e)  

rhoE = rho*(e + 0.5*magSqr(U));
Here, turbulence->alphaEff() is defined as kappaEff/Cp and using energy definition as sensibleInternalEnergy in constant/thermophysicalProperties file, e=Cv*T.
So, above equation effectively translates into

\frac{\partial}{\partial x_j}(\frac{\kappa_{Eff}}{\gamma}\frac{\partial T}{\partial x_j})

But, it should be only

\frac{\partial}{\partial x_j}(\kappa_{Eff}\frac{\partial T}{\partial x_j}).

Since, it would not make much difference if \gamma is constant. But when using Janaf tables at high temperature, these two gives different solutions because \gamma varies with temperature. Which one is correct? Can somebody help?

I have same question regarding calculation of wallHeatFlux. Should it be

calcWallHeatFlux(turbModel.alphaEff(), turbModel.transport().he())
calcWallHeatFlux(turbModel.kappaEff(), turbModel.transport().T())

in file at line 220

febriyan91 likes this.

Last edited by Daedalus; January 29, 2019 at 06:57. Reason: formating error
Daedalus is offline   Reply With Quote


janafthermo, rhocentralfoam, wallheattransfer

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
OpenFOAM Heat Insulation Boundary Condition for Heat Conduction in Solids abkahraman OpenFOAM Running, Solving & CFD 3 August 15, 2020 11:11
Heat Release vs. Apparent Heat Release // Comparison with GT Power MFGT CONVERGE 10 August 17, 2016 03:38
Problem with total heat transfer rate aswathy_raghu FLUENT 0 July 26, 2016 08:39
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 07:28
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 08:00

All times are GMT -4. The time now is 10:27.