流体力学方程组
贡献者： 待更新
描述流体微元的物理量有 $\rho,p,T,e, \boldsymbol{\mathbf{u}} $ 等，而完整描述它们随时间的演化需要流体力学方程组。它们分别是质量守恒方程（式 3 ）、动量守恒方程（式 5 ）、动能守恒方程（式 13 ）、内能守恒方程（式 15 ）、本构方程（式 2 ）和状态方程 $\rho=\rho(p,T)$。
\begin{equation}
\begin{aligned}
& \frac{\partial \rho}{\partial t} +\nabla\cdot(\rho \boldsymbol{\mathbf{u}} )=0\\
&
\rho \frac{\mathrm{d}{u_j}}{\mathrm{d}{t}} =- \frac{\partial p}{\partial x_j} +\rho g_j+\mu \frac{\partial^{2}{u_j}}{\partial{x_i}^{2}} + \left(\mu_\nu+\frac{1}{3}\mu \right) \frac{\partial}{\partial{x_j}} \frac{\partial u_m}{\partial x_m} \\
&\rho \frac{\mathrm{d}{}}{\mathrm{d}{t}} \left(\frac{1}{2}| \boldsymbol{\mathbf{u}} |^2 \right) =\rho g_iu_i+u_j \frac{\partial T_{ij}}{\partial x_i} =\rho g_iu_i+u_j \left(- \frac{\partial p}{\partial x_j} + \frac{\partial \tau_{ij}}{\partial x_i} \right) \\
&\rho \frac{\mathrm{d}{e}}{\mathrm{d}{t}} =-p \frac{\partial u_m}{\partial x_m} +2\mu \left(S_{ij}-\frac{1}{3} \frac{\partial u_m}{\partial x_m} \delta_{ij} \right) ^2+\mu_\nu \left( \frac{\partial u_m}{\partial x_m} \right) ^2+ \frac{\partial}{\partial{x_i}} \left(k \frac{\partial T}{\partial x_i} \right)
\\
&T_{ij}=-p\delta_{ij}+\tau_{ij}=-p\delta_{ij}+2\mu \left(S_{ij}-\frac{2}{3}S_{mm}\delta_{ij} \right) +\mu_\nu S_{mm} \delta_{ij}\\
&\rho=\rho(p,T)~.
\end{aligned}
\end{equation}
上面的本构方程和状态方程可以代入随时间演化的方程，随时间演化的为前 $4$ 行（其中动量守恒方程是三分量的），共 $6$ 个方程；其中动能守恒方程可以由动量守恒方程推出，所以可以舍去，总共可以得到 $5$ 个分量的独立的方程。现在再来看变量个数，$\rho,p,T,e, \boldsymbol{\mathbf{u}} $ 一共 $7$ 个分量，其中 $e,\rho$ 都蕴含在状态方程之中：$\rho=\rho(p,T),e=e(p,T)$，所以实际上独立的变量个数是 $5$ 个。$5$ 个独立的方程和 $5$ 个独立的变量，构成了完备的流体力学方程组。
于是，这些方程可以用于作数值模拟。但是数值计算的计算能力和模拟精度是有限的。抛开数值计算，这些流体力学方程是相当复杂的，以至于即使人们研究相当简单的模型，也几乎无法得到解析解。因此我们会首先研究一些流体力学方程组能很好地解决和描述的问题，例如涡旋、伯努利公式、重力波等。