流体力学方程组

贡献者：待更新

预备知识　Navier-Stokes 方程

　　描述流体微元的物理量有 $ρ, p, T, e, u$ 等，而完整描述它们随时间的演化需要流体力学方程组。它们分别是质量守恒方程（式 3 ）、动量守恒方程（式 5 ）、动能守恒方程（式 13 ）、内能守恒方程（式 15 ）、本构方程（式 2 ）和状态方程 $ρ = ρ (p, T)$ 。

\begin{matrix} (1) & \begin{aligned} \frac{\partial ρ}{\partial t} + \nabla \cdot (ρ u) = 0 \\ ρ \frac{d u_{j}}{d t} = - \frac{\partial p}{\partial x_{j}} + ρ g_{j} + μ \frac{\partial^{2} u_{j}}{\partial {x_{i}}^{2}} + (μ_{ν} + \frac{1}{3} μ) \frac{\partial}{\partial x_{j}} \frac{\partial u_{m}}{\partial x_{m}} \\ ρ \frac{d}{d t} (\frac{1}{2} | u |^{2}) = ρ g_{i} u_{i} + u_{j} \frac{\partial T_{i j}}{\partial x_{i}} = ρ g_{i} u_{i} + u_{j} (- \frac{\partial p}{\partial x_{j}} + \frac{\partial τ_{i j}}{\partial x_{i}}) \\ ρ \frac{d e}{d t} = - p \frac{\partial u_{m}}{\partial x_{m}} + 2 μ {(S_{i j} - \frac{1}{3} \frac{\partial u_{m}}{\partial x_{m}} δ_{i j})}^{2} + μ_{ν} {(\frac{\partial u_{m}}{\partial x_{m}})}^{2} + \frac{\partial}{\partial x_{i}} (k \frac{\partial T}{\partial x_{i}}) \\ T_{i j} = - p δ_{i j} + τ_{i j} = - p δ_{i j} + 2 μ (S_{i j} - \frac{2}{3} S_{m m} δ_{i j}) + μ_{ν} S_{m m} δ_{i j} \\ ρ = ρ (p, T) . \end{aligned} \end{matrix}

上面的本构方程和状态方程可以代入随时间演化的方程，随时间演化的为前

4

行（其中动量守恒方程是三分量的），共

6

个方程；其中动能守恒方程可以由动量守恒方程推出，所以可以舍去，总共可以得到

5

个分量的独立的方程。现在再来看变量个数，

ρ, p, T, e, u

一共

7

个分量，其中

e, ρ

都蕴含在状态方程之中：

ρ = ρ (p, T), e = e (p, T)

，所以实际上独立的变量个数是

5

个。

5

个独立的方程和

5

个独立的变量，构成了完备的流体力学方程组。

　　于是，这些方程可以用于作数值模拟。但是数值计算的计算能力和模拟精度是有限的。抛开数值计算，这些流体力学方程是相当复杂的，以至于即使人们研究相当简单的模型，也几乎无法得到解析解。因此我们会首先研究一些流体力学方程组能很好地解决和描述的问题，例如涡旋、伯努利公式、重力波等。

致读者：小时百科一直以来坚持所有内容免费无广告，这导致我们处于严重的亏损状态。长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。因此，我们请求广大读者热心打赏 ，使网站得以健康发展。如果看到这条信息的每位读者能慷慨打赏 20 元，我们一周就能脱离亏损，并在接下来的一年里向所有读者继续免费提供优质内容。但遗憾的是只有不到 1% 的读者愿意捐款，他们的付出帮助了 99% 的读者免费获取知识，我们在此表示感谢。