Understanding energy initialization

chriss85 · August 7, 2014, 08:47

Hello,

I'm currently trying to understand how the energy fields are initialized from the initial temperature field.

In heThermo.C, init() the fields are calculated. I don't understand why the cells are initialized with HE() ( -> J/kg) while the boundary is initialized with he() ( -> J/kmol).

Can anyone clarify this?

Background: I'm using tabulated material data for enthalpy and I'm trying to get rid of the root finding method in the temperature calculation, given a function h(T, p), by directly using a T(h, p) function. I hope to get increased solver stability and possibly speed.

chriss85 · August 7, 2014, 08:53

I just noticed that the called he() function is actually not the one from thermo.H but rather from heThermo.H, and this one uses J/kg. So it's actually just an unlucky naming convention

.

August 7, 2014, 08:47	Understanding energy initialization	#1
chriss85 Senior Member Join Date: Oct 2013 Posts: 397 Rep Power: 18	Hello, I'm currently trying to understand how the energy fields are initialized from the initial temperature field. In heThermo.C, init() the fields are calculated. I don't understand why the cells are initialized with HE() ( -> J/kg) while the boundary is initialized with he() ( -> J/kmol). Can anyone clarify this? Background: I'm using tabulated material data for enthalpy and I'm trying to get rid of the root finding method in the temperature calculation, given a function h(T, p), by directly using a T(h, p) function. I hope to get increased solver stability and possibly speed.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
energy in sonicFoam	joern	OpenFOAM Running, Solving & CFD	1	September 24, 2019 00:15
energy spectra from LES simulations	kumar	OpenFOAM Post-Processing	2	February 23, 2016 12:25
Natural convection, problems with total energy balance	MdoNascimento	STAR-CCM+	0	March 13, 2014 12:15
ATTENTION! Reliability problems in CFX 5.7	Joseph	CFX	14	April 20, 2010 15:45
SIMPLE and energy equation convergence	Fabio	Main CFD Forum	0	June 1, 2007 06:06

August 7, 2014, 08:53		#2
chriss85 Senior Member Join Date: Oct 2013 Posts: 397 Rep Power: 18	I just noticed that the called he() function is actually not the one from thermo.H but rather from heThermo.H, and this one uses J/kg. So it's actually just an unlucky naming convention .