Conservation of Internal Energy

Conservation of internal energy is used to model temperature dependent material behavior. It also allows an energy balance evaluation.

However, internal energy is only calculated if it is turned on, to reduce computation time in problems not involving heat transfer.

The conservation of energy is given by:(1)

\frac{\partial ρ e}{\partial t} - ((w_{i} - v_{i}) \cdot \frac{\partial ρ e}{\partial x_{i}}) + (ρ e + P) \frac{\partial v_{K}}{\partial x_{K}} = 0

Where,

Applying a Galerkin variation form for the solution gives:(2)

\int_{}^{V} ψ (\frac{\partial ρ e}{\partial t} - ((w_{i} - v_{i}) \cdot \frac{\partial ρ e}{\partial x_{i}}) + (ρ e + P) \frac{\partial v_{K}}{\partial x_{K}}) = 0

Making the following assumptions:

$ψ$ =1

$ρ_{e}$ = constant over control volume $V$

Equation 2 reduces to:(3)

\int_{V}^{} \frac{\partial ρ e}{\partial t} d V + \int_{V}^{} (ρ e + P) \frac{\partial v_{K}}{\partial x_{K}} d V = 0

Applying the divergence theorem gives:(4)

\int_{}^{V} \frac{\partial ρ e}{\partial t} d V + \int_{}^{S} (ρ e (v_{j} \cdot n_{j}) d S) + \int_{}^{V} P \frac{\partial v_{K}}{\partial x_{K}} d V = 0

Hence:(5)

\frac{d}{d t} \int_{V}^{} ρ e d V = \int_{S}^{} ρ e (w_{j} - v_{j}) n_{j} d S - \int_{V}^{} P \frac{\partial v_{K}}{\partial x_{K}} d V

This formula is still valid if $e$ is not assumed uniform over volume $V$ .

New formulation based on total energy conservation are under investigation.