Governing equations on differential form

The Differential Equations on Conservation Form

Conservation of Mass

Apply Gauss's divergence theorem on the surface integral in Eqn. \ref{eq:governing:cont:int} gives

$$\underset{\partial \Omega }{\iint }\rho 𝐯\cdot 𝐧dS=\underset{\Omega }{\iiint }\mathrm{\nabla }\cdot (\rho 𝐯)dV$$

Also, if $\Omega$ is a fixed control volume

$$\frac{d}{dt}\underset{\Omega }{\iiint }\rho dV=\underset{\Omega }{\iiint }\frac{\partial \rho }{\partial t}dV$$

The continuity equation can now be written as a single volume integral.

$$\underset{\Omega }{\iiint }[\frac{\partial \rho }{\partial t}+\mathrm{\nabla }\cdot (\rho 𝐯)]dV=0$$

$\Omega$ is an arbitrary control volume and thus

$$\frac{\partial \rho }{\partial t}+\mathrm{\nabla }\cdot (\rho 𝐯)=0$$

which is the continuity equation on partial differential form.

Conservation of Momentum

As for the continuity equation, the surface integral terms are rewritten as volume integrals using Gauss's divergence theorem.

$$\underset{\partial \Omega }{\iint }(\rho 𝐯\cdot 𝐧)𝐯dS=\underset{\Omega }{\iiint }\mathrm{\nabla }\cdot (\rho 𝐯𝐯)dV$$

$$\underset{\partial \Omega }{\iint }p𝐧dS=\underset{\Omega }{\iiint }\mathrm{\nabla }pdV$$

Also, if $\Omega$ is a fixed control volume

$$\frac{d}{dt}\underset{\Omega }{\iiint }\rho 𝐯dV=\underset{\Omega }{\iiint }\frac{\partial }{\partial t}(\rho 𝐯)dV$$

The momentum equation can now be written as one single volume integral

$$\underset{\Omega }{\iiint }[\frac{\partial }{\partial t}(\rho 𝐯)+\mathrm{\nabla }\cdot (\rho 𝐯𝐯)+\mathrm{\nabla }p-\rho 𝐟]dV=0$$

$\Omega$ is an arbitrary control volume and thus

$$\frac{\partial }{\partial t}(\rho 𝐯)+\mathrm{\nabla }\cdot (\rho 𝐯𝐯)+\mathrm{\nabla }p=\rho 𝐟$$

which is the momentum equation on partial differential form

Conservation of Energy

Gauss's divergence theorem applied to the surface integral term in the energy equation (Eqn. \ref{eq:governing:energy:int}) gives

$$\underset{\partial \Omega }{\iint }\rho h_o(𝐯\cdot 𝐧)dS=\underset{\Omega }{\iiint }\mathrm{\nabla }\cdot (\rho h_o𝐯)dV$$

Fixed control volume

$$\frac{d}{dt}\underset{\Omega }{\iiint }\rho e_{o}dV=\underset{\Omega }{\iiint }\frac{\partial }{\partial t}(\rho e_o)dV$$

The energy equation can now be written as

$$\underset{\Omega }{\iiint }[\frac{\partial }{\partial t}(\rho e_o)+\mathrm{\nabla }\cdot (\rho h_o𝐯)-\rho 𝐟\cdot 𝐯-\dot{q}\rho ]dV=0$$

$\Omega$ is an arbitrary control volume and thus

$$\frac{\partial }{\partial t}(\rho e_o)+\mathrm{\nabla }\cdot (\rho h_o𝐯)=\rho 𝐟\cdot 𝐯+\dot{q}\rho$$

which is the energy equation on partial differential form

Summary

The governing equations for compressible inviscid flow on partial differential form:

$$\frac{\partial \rho }{\partial t}+\mathrm{\nabla }\cdot (\rho 𝐯)=0$$

$$\frac{\partial }{\partial t}(\rho 𝐯)+\mathrm{\nabla }\cdot (\rho 𝐯𝐯)+\mathrm{\nabla }p=\rho 𝐟$$

$$\frac{\partial }{\partial t}(\rho e_o)+\mathrm{\nabla }\cdot (\rho h_o𝐯)=\rho 𝐟\cdot 𝐯+\dot{q}\rho$$

The Differential Equations on Non-Conservation Form

The Substantial Derivative

The substantial derivative operator is defined as

$$\frac{D}{Dt}=\frac{\partial }{\partial t}+𝐯\cdot \mathrm{\nabla }$$

where the first term of the right hand side is the local derivative and the second term is the convective derivative.

Conservation of Mass

If we apply the substantial derivative operator to density we get

$$\frac{D\rho }{Dt}=\frac{\partial \rho }{\partial t}+𝐯\cdot \mathrm{\nabla }\rho$$

From before we have the continuity equation on differential form as

$$\frac{\partial \rho }{\partial t}+\mathrm{\nabla }\cdot (\rho 𝐯)=0$$

which can be rewritten as

$$\frac{\partial \rho }{\partial t}+\rho (\mathrm{\nabla }\cdot 𝐯)+𝐯\cdot \mathrm{\nabla }\rho =0$$

and thus

$$\frac{D\rho }{Dt}+\rho (\mathrm{\nabla }\cdot 𝐯)=0$$

Eqn. \ref{eq:governing:cont:non} says that the mass of a fluid element with a fixed set of fluid particles is constant as the element moves in space.

Conservation of Momentum

We start from the momentum equation on differential form derived above

$$\frac{\partial }{\partial t}(\rho 𝐯)+\mathrm{\nabla }\cdot (\rho 𝐯𝐯)+\mathrm{\nabla }p=\rho 𝐟$$

Expanding the first and the second terms gives

$$\rho \frac{\partial 𝐯}{\partial t}+𝐯\frac{\partial \rho }{\partial t}+\rho 𝐯\cdot \mathrm{\nabla }𝐯+𝐯(\mathrm{\nabla }\cdot \rho 𝐯)+\mathrm{\nabla }p=\rho 𝐟$$

Collecting terms, we can identify the substantial derivative operator applied to the velocity vector and the continuity equation.

$$\rho \underset{=\frac{D𝐯}{Dt}}{\underbrace{[\frac{\partial 𝐯}{\partial t}+𝐯\cdot \mathrm{\nabla }𝐯]}}+𝐯\underset{=0}{\underbrace{[\frac{\partial \rho }{\partial t}+\mathrm{\nabla }\cdot \rho 𝐯]}}+\mathrm{\nabla }p=\rho 𝐟$$

which gives us the non-conservation form of the momentum equation

$$\frac{D𝐯}{Dt}+\frac{1}{\rho }\mathrm{\nabla }p=𝐟$$

Conservation of Energy

The last equation on non-conservation differential form is the energy equation. We start by rewriting the energy equation on differential form (Eqn. \ref{eq:governing:energy:pde}), repeated here for convenience

$$\frac{\partial }{\partial t}(\rho e_o)+\mathrm{\nabla }\cdot (\rho h_o𝐯)=\rho 𝐟\cdot 𝐯+\dot{q}\rho$$

Total enthalpy, $h_o$, is replaced with total energy, $e_o$

$$h_o=e_o+\frac{p}{\rho }$$

which gives

$$\frac{\partial }{\partial t}(\rho e_o)+\mathrm{\nabla }\cdot (\rho e_o𝐯)+\mathrm{\nabla }\cdot (p𝐯)=\rho 𝐟\cdot 𝐯+\dot{q}\rho$$

Expanding the two first terms as

$$\rho \frac{\partial e_o}{\partial t}+e_o\frac{\partial \rho }{\partial t}+\rho 𝐯\cdot \mathrm{\nabla }{e}_o+e_o\mathrm{\nabla }\cdot (\rho 𝐯)+\mathrm{\nabla }\cdot (p𝐯)=\rho 𝐟\cdot 𝐯+\dot{q}\rho$$

Collecting terms, we can identify the substantial derivative operator applied on total energy, $D{e}_o/Dt$ and the continuity equation

$$\rho \underset{=\frac{D{e}_o}{Dt}}{\underbrace{[\frac{\partial e_o}{\partial t}+𝐯\cdot \mathrm{\nabla }{e}_o]}}+e_o\underset{=0}{\underbrace{[\frac{\partial \rho }{\partial t}+\mathrm{\nabla }\cdot (\rho 𝐯)]}}+\mathrm{\nabla }\cdot (p𝐯)=\rho 𝐟\cdot 𝐯+\dot{q}\rho$$

and thus we end up with the energy equation on non-conservation differential form

\begin{equation} \rho\frac{De_o}{Dt} + \nabla\cdot(p\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho \label{eq:governing:energy:non} \end{equation}\\

%\section*{The Governing Equations on Differential Non-Conservation Form} % %\vspace*{1cm} % %\noindent Continuity: % %\begin{equation} %\frac{D\rho}{Dt}+\rho(\nabla\cdot\mathbf{v})=0 %\label{eq:governing:cont:non} %\end{equation}\\ % %\noindent Momentum: % %\begin{equation} %\frac{D\mathbf{v}}{Dt}+\frac{1}{\rho}\nabla p = \mathbf{f} %\label{eq:governing:mom:non} %\end{equation}\\ % %\noindent Energy: % %\begin{equation} %\rho\frac{De_o}{Dt} + \nabla\cdot(p\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho %\label{eq:governing:energy:non} %\end{equation}\\

Alternative Forms of the Energy Equation

Internal Energy Formulation

\noindent Total internal energy is defined as\\

\[e_o=e+\frac{1}{2}\mathbf{v}\cdot\mathbf{v}\]\\

\noindent Inserted in Eqn. \ref{eq:governing:energy:non}, this gives

\[\rho\frac{De}{Dt} + \rho\mathbf{v}\cdot\frac{D \mathbf{v}}{Dt} + \nabla\cdot(p\mathbf{v}) = \rho\mathbf{f}\cdot\mathbf{v} + \dot{q}\rho\]\\

\noindent Now, let's replace the substantial derivative $D\mathbf{v}/Dt$ using the momentum equation on non-conservation form (Eqn. \ref{eq:governing:mom:non}).\\

\[\rho\frac{De}{Dt} -\mathbf{v}\cdot\nabla p + \cancel{\rho\mathbf{f}\cdot\mathbf{v}} + \nabla\cdot(p\mathbf{v}) = \cancel{\rho\mathbf{f}\cdot\mathbf{v}} + \dot{q}\rho\]\\

\noindent Now, expand the term $\nabla\cdot(p\mathbf{v})$ gives\\

\[\rho\frac{De}{Dt} \cancel{-\mathbf{v}\cdot\nabla p} + \cancel{\mathbf{v}\cdot\nabla p} + p(\nabla\cdot\mathbf{v}) = \dot{q}\rho\Rightarrow \rho\frac{De}{Dt} + p(\nabla\cdot\mathbf{v}) = \dot{q}\rho\]\\

\noindent Divide by $\rho$\\

\begin{equation} \frac{De}{Dt} + \frac{p}{\rho}(\nabla\cdot\mathbf{v}) = \dot{q} \label{eq:governing:energy:non:b} \end{equation}\\

\noindent Conservation of mass gives\\

\[\frac{D\rho}{Dt}+\rho(\nabla\cdot\mathbf{v})=0\Rightarrow \nabla\cdot\mathbf{v} = -\frac{1}{\rho}\frac{D\rho}{Dt}\]

\noindent Insert in Eqn. \ref{eq:governing:energy:non:b}\\

\[\frac{De}{Dt} - \frac{p}{\rho^2}\frac{D\rho}{Dt} = \dot{q}\Rightarrow \frac{De}{Dt} + p\frac{D}{Dt} \left(\frac{1}{\rho}\right)= \dot{q}\]\\

\begin{equation} \frac{De}{Dt} + p\frac{D\nu}{Dt} = \dot{q} \label{eq:governing:energy:non:b} \end{equation}\\

\noindent Compare with the first law of thermodynamics: $de=\delta q-\delta w$\\

Enthalpy Formulation

\[h=e+\frac{p}{\rho}\Rightarrow \frac{Dh}{Dt}=\frac{De}{Dt}+\frac{1}{\rho}\frac{Dp}{Dt}+p\frac{D}{Dt}\left(\frac{1}{\rho}\right)\]\\

\noindent with $De/Dt$ from Eqn. \ref{eq:governing:energy:non:b}\\

\[\frac{Dh}{Dt}=\dot{q} - \cancel{p\frac{D}{Dt}\left(\frac{1}{\rho}\right)} +\frac{1}{\rho}\frac{Dp}{Dt}+\cancel{p\frac{D}{Dt}\left(\frac{1}{\rho}\right)}\]\\

\begin{equation} \frac{Dh}{Dt}=\dot{q} + \frac{1}{\rho}\frac{Dp}{Dt} \label{eq:governing:energy:non:c} \end{equation}\\

Total Enthalpy Formulation

\[h_o=h+\frac{1}{2}\mathbf{v}\mathbf{v}\Rightarrow\frac{Dh_o}{Dt}=\frac{Dh}{Dt}+\mathbf{v}\cdot\frac{D\mathbf{v}}{Dt}\]\\

\noindent From the momentum equation (Eqn. \ref{eq:governing:mom:non})\\

\[\frac{D\mathbf{v}}{Dt}=\mathbf{f}-\frac{1}{\rho}\nabla p\]\\

\noindent which gives\\

\[\frac{Dh_o}{Dt}=\frac{Dh}{Dt}+\mathbf{v}\cdot\mathbf{f} -\frac{1}{\rho}\mathbf{v}\cdot\nabla p \]\\

\noindent Inserting $Dh/Dt$ from Eqn. \ref{eq:governing:energy:non:c} gives\\

\[\frac{Dh_o}{Dt}=\dot{q} + \frac{1}{\rho}\frac{Dp}{Dt}+\mathbf{v}\cdot\mathbf{f} -\frac{1}{\rho}\mathbf{v}\cdot\nabla p = \frac{1}{\rho}\left[\frac{Dp}{Dt}-\mathbf{v}\cdot\nabla p\right] + \dot{q} + \mathbf{v}\cdot\mathbf{f}\]\\

\noindent The substantial derivative operator applied to pressure\\

\[\frac{Dp}{Dt}=\frac{\partial p}{\partial t}+\mathbf{v}\cdot\nabla p\]\\

\noindent and thus\\

\[\frac{Dp}{Dt}-\mathbf{v}\cdot\nabla p=\frac{\partial p}{\partial t}\]\\

\noindent which gives\\

\[\frac{Dh_o}{Dt}=\frac{1}{\rho}\frac{\partial p}{\partial t} + \dot{q} + \mathbf{v}\cdot\mathbf{f}\]\\

\noindent If we assume adiabatic flow without body forces\\

\[\frac{Dh_o}{Dt}=\frac{1}{\rho}\frac{\partial p}{\partial t}\]\\

\noindent If we further assume the flow to be steady state we get\\

\[\frac{Dh_o}{Dt}=0\]\\

\noindent This means that in a steady-state adiabatic flow without body forces, total enthalpy is constant along a streamline.\\