1.2.1流动与传热的通用控制方程（OpenFOAM理论笔记系列）

最新推荐文章于 2026-07-16 09:19:53 发布

原创

最新推荐文章于 2026-07-16 09:19:53 发布 · 8k 阅读

1.2系统的控制方程

1.2.1流动与传热的通用控制方程

通用的流体流动与传热方程如下：

连续性方程：
${\partial \rho}\over{\partial t}}+\nabla\cdot(\rho \vec U)=0 \tag {1.4}$

动量方程：
${\partial (\rho\vec U)}\over{\partial t}}+\nabla\cdot(\rho \vec U\vec U)=\nabla\cdot\tau-\nabla p+\rho \vec g \tag {1.5}$

能量方程：
${\partial (\rho e)}\over{\partial t}}+\nabla\cdot(\rho e\vec U)=-\nabla\cdot\vec q-\nabla\cdot(p\vec U)+\nabla\cdot(\tau\cdot\vec U)+\rho\vec g\cdot\vec U+\rho Q \tag {1.6}$

状态方程：
$p=p(\rho,T)\space,\space Ｅ_i=E_i(rho,T) \tag {1.7}$

接下来我们对上述方程依次进行讨论：

对于连续性方程， $\rho$ 为流体的密度， $\vec U$ 为流体的速度。在不可压缩的情况下，密度不随时空坐标变化，其可化简为：
$\nabla \cdot \vec U =0 \tag{1.8}$
对于动量方程， $\tau$ 为粘性剪切应力， $\vec g$ 是重力加速度。当流体不可压缩时，密度可以提到各导数运算符外。为了能够求解动量方程，我们必须知道粘性剪切应力的表达形式，常见的流体为牛顿流体，其粘性剪切应力与形变遵循线性关系。考虑不可压缩流体，其形变率张量可以写为:
$\textbf{S}={1\over2}(\nabla\vec U+{\nabla\vec U}^T) \tag{1.9}$
粘性剪切应力与剪切率成线性关系，即：
$\tau=2\mu S=\mu\left[\nabla\vec U+(\nabla\vec U)^T\right] \tag{1.10}$
$\mu$ 即动力粘度。对于可压缩的牛顿流体， $\tau$ 的表达形式会更加复杂并且要把压力也归入粘性应力的作用中,我们用 $\sigma$ 表示完整的牛顿粘性应力则有：
$\sigma=-(p+{2\over3}\mu\nabla\cdot\vec U)\textbf{I}+\mu\left[\nabla\vec U+(\nabla\vec U)^T\right] \tag{1.11}$