Physics:First law of thermodynamics (fluid mechanics)

In physics, the first law of thermodynamics is an expression of the conservation of total energy of a system. The increase of the energy of a system is equal to the sum of work done on the system and the heat added to that system:

[math]\displaystyle{ dE_t=dQ + dW }[/math]

where

• [math]\displaystyle{ E_t }[/math] is the total energy of a system.
• [math]\displaystyle{ W }[/math] is the work done on it.
• [math]\displaystyle{ Q }[/math] is the heat added to that system.

In fluid mechanics, the first law of thermodynamics takes the following form:

[math]\displaystyle{ \frac{D E_t}{D t}=\frac{D W}{D t} + \frac{D Q}{D t} \to \frac{D E_t}{D t} = \nabla\cdot({\mathbf \sigma\cdot v}) - \nabla\cdot{\mathbf q} }[/math]

where

• [math]\displaystyle{ \mathbf\sigma }[/math] is the Cauchy stress tensor.
• [math]\displaystyle{ \mathbf v }[/math] is the flow velocity.
• and [math]\displaystyle{ \mathbf q }[/math] is the heat flux vector.

Because it expresses conservation of total energy, this is sometimes referred to as the energy balance equation of continuous media. The first law is used to derive the non-conservation form of the Navier–Stokes equations.

Note

[math]\displaystyle{ {\mathbf\sigma} = -p{\mathbf I} + {\mathbf T} }[/math]

Where

• [math]\displaystyle{ p }[/math] is the pressure
• [math]\displaystyle{ \mathbf I }[/math] is the identity matrix
• [math]\displaystyle{ \mathbf T }[/math] is the deviatoric stress tensor

That is, pulling is positive stress and pushing is negative stress.

Compressible fluid

For a compressible fluid the left hand side of equation becomes:

[math]\displaystyle{ \frac{D E_t}{D t} = \frac{\partial E}{\partial t}+\nabla\cdot(E\mathbf v) }[/math]

because in general

[math]\displaystyle{ \nabla\cdot\mathbf v\neq 0. }[/math]

Integral form

[math]\displaystyle{ \int_V \frac{\partial E}{\partial t} \, dV =-\oint_{\partial V}E{\mathbf v}\cdot d{\mathbf A}+\oint_{\partial V}({\mathbf \sigma\cdot v})\cdot d{\mathbf A}-\oint_{\partial V}{\mathbf q}\cdot d{\mathbf A} }[/math]

That is, the change in the internal energy of the substance within a volume is the negative of the amount carried out of the volume by the flow of material across the boundary plus the work done compressing the material on the boundary minus the flow of heat out through the boundary. More generally, it is possible to incorporate source terms.

Alternative representation

[math]\displaystyle{ \rho\frac{D h}{D t}=\frac{D p}{D t}+ \nabla\cdot(k \, \nabla T)+\Phi }[/math]

where [math]\displaystyle{ h }[/math] is specific enthalpy, [math]\displaystyle{ \Phi={\mathbf \tau}:\nabla{\mathbf v} }[/math] is dissipation function and [math]\displaystyle{ T }[/math] is temperature. And where

• [math]\displaystyle{ E_t=\rho(e+\frac{1}{2}v^2-\mathbf{g\cdot r}) }[/math]

i.e. internal energy per unit volume equals mass density times the sum of: proper energy per unit mass, kinetic energy per unit mass, and gravitational potential energy per unit mass.

• [math]\displaystyle{ -\nabla\cdot{\mathbf q}=+\nabla\cdot(k\,\nabla T) }[/math]

i.e. change in heat per unit volume (negative divergence of heat flow) equals the divergence of heat conductivity times the gradient of the temperature.

• [math]\displaystyle{ \nabla\cdot({\mathbf \sigma\cdot v})={\mathbf v\cdot\nabla\cdot\sigma}+\sigma :\nabla \mathbf v }[/math]

i.e. divergence of work done against stress equals flow of material times divergence of stress plus stress times divergence of material flow.

• [math]\displaystyle{ \sigma :\nabla {\mathbf v}=\Phi - p\,\nabla\cdot{\mathbf v} }[/math]

i.e. stress times divergence of material flow equals deviatoric stress tensor times divergence of material flow minus pressure times material flow.

• [math]\displaystyle{ h=e+\frac{p}{\rho} }[/math]

i.e. enthalpy per unit mass equals proper energy per unit mass plus pressure times volume per unit mass (reciprocal of mass density).

• [math]\displaystyle{ \nabla\cdot\sigma=\frac{DJ}{Dt}-\mathbf f }[/math]
• [math]\displaystyle{ -p\,\nabla\cdot{\mathbf v}=\frac{Dp}{Dt}-\rho\frac{D}{Dt} \left( \frac{p}{\rho}\right) }[/math]

Alternative form data

• [math]\displaystyle{ \nabla\cdot{\mathbf\sigma} = \frac{D{\mathbf J}}{D t} - \mathbf f }[/math] left hand side of Navier–Stokes equations minus body force (per unit volume) acting on fluid.
• [math]\displaystyle{ -p\nabla\cdot{\mathbf v}=\frac{D p}{D t}-\rho\frac{D}{D t} \left( \frac{p}{\rho} \right) }[/math] this relation is derived using this relationship [math]\displaystyle{ \rho\nabla\cdot{\mathbf v}=-\frac{D \rho}{D t} }[/math] which is alternative form of continuity equation [math]\displaystyle{ \frac{\partial\rho}{\partial t}+\nabla\cdot{\mathbf J} = 0 }[/math]

See also

• Clausius–Duhem inequality
• Continuum mechanics
• First law of thermodynamics
• Material derivative
• Incompressible flow