# Energy equation

### From Thermal-FluidsPedia

The integral energy equation is

The surface integrals on the right-hand side of the above equation can be rewritten as volume integrals using

from transformation formula, i.e.,

Substituting eqs. (1) – (3) into the integral energy equation and considering yields

For eq. (4) to be true for any arbitrary control volume, the integrand must be zero. The result is the general differential form of the energy equation:

The left-hand side of eq. (5) can be rewritten as

According to the continuity eq. from continuity equation, the first bracketed term on the left-hand side of eq. (6) is zero. The second term on the left-hand side may be written more simply in substantial derivative form, i.e.,

which is the total energy (including thermal and mechanical energies) balance equation.
It would be more convenient to remove the mechanical energy terms from eq. (7). The mechanical energy balance equation can be obtained by forming a dot (scalar) product of the momentum equation from momentum equation with the velocity vector *V*_{rel}, i.e.,

which can be rearranged to obtain

Subtracting eq. (9) from eq. (7) yields the following thermal energy equation:

where the left-hand side represents the rate of gain of energy per unit volume. The terms in the right-hand side are the rate of energy input by heat transfer per unit volume, the internal heat generation per unit volume, the reversible rate of energy increase per unit volume by compression, and the irreversible rate of energy increase per unit volume by viscous dissipation, respectively; the viscous dissipation is the contraction [see eq. (G.35)] of two tensors of rank two: and , i.e.,

Equation (10) is the energy equation expressed in terms of internal energy. In order to obtain an equation that contains enthalpy, the definition of enthalpy is employed:

Substituting eq. (12) into eq. (10), and considering continuity equation from continuity equation, the energy equation in terms of enthalpy and temperature is obtained:

For a multicomponent system, the enthalpy can be expressed as

where ω_{i} and *h*_{i} are the mass fraction and specific enthalpy of the *i*^{th} component. Substituting eq. (14) into eq. (13), the energy equation for a multicomponent system becomes

The heat flux vector can be obtained from Fourier's law, :

Substituting eq. (16) into eq. (15), the energy equation becomes

where the second and third terms on the right hand side, which represent the contributions of interdiffusional convection and Dufour effect, are both negligible for most applications.

For a pure substance, eq. (17) is reduced to:

The enthalpy of a pure substance (single component) can be expressed as a function of temperature and pressure, *h* = *f*(*T*,*p*), i.e.,

Thermodynamic relations give us

where

is the coefficient of thermal expansion. Therefore, eq. (19) becomes

Substituting eq. (22) into eq. (18), the energy equation becomes

which can be simplified for different substances as demonstrated below.

For ideal gases, substituting the equation of state for ideal gas [ρ = *p* / (*R*_{g}*T*)] into eq. (21) yields β = 1 / *T* and eq. (23) becomes

For incompressible fluids, the constant density yields β = 0 and eq. (23) becomes

For solids, the density is constant and the velocity is zero, and the energy equation becomes

which is a heat conduction equation.

## References

Faghri, A., and Zhang, Y., 2006, *Transport Phenomena in Multiphase Systems,* Elsevier, Burlington, MA

Faghri, A., Zhang, Y., and Howell, J. R., 2010, *Advanced Heat and Mass Transfer,* Global Digital Press, Columbia, MO.