能量方程浅析

1.能量方程

根据热力学第一定律 ( 能量守恒定律 )，有：流体微团内能量的变化率=流人微团内的净热流量+体积力和表面力对微团做功的功率，可表述为：
$$E=B+C\text{\hspace{1em}(1-1)}$$

其中，E包含两部分：

• 由于分子随机运动而产生的(单位质量)内能 e；
• 流体微团平动时具有的动能。单位质量的动能为 $K={\mathbf{U}}^2/2$。

$$E=e+{\mathbf{U}}^2/2\text{\hspace{1em}(1-2)}$$

B部分的热流也包括两部分：

• 由于体积加热，如吸收或释放的辐射热（内部产生的能量）：$\rho r$，其中，r为单位质量的体积加热率；
• 由温度梯度导致的跨过表面的热输运，即热传导（流入流出的能量）：$\nabla \cdot \mathbf{q}=-k\cdot \nabla T$，k 为热导率
  $$B=\rho r+k\cdot \nabla T\text{\hspace{1em}(1-3)}$$

C部分为做功所消耗的能量，分为两部分：

• 体积力f做功：$\rho \mathbf{f}\cdot \mathbf{U}$；
• 表面力(压力p+应力$\tau$)做功：$-\nabla \cdot (p\mathbf{U})+\nabla \cdot (\tau \cdot \mathbf{U})$。

$$C=\rho \mathbf{f}\cdot \mathbf{U}-\nabla \cdot (p\mathbf{U})+\nabla \cdot (\tau \cdot \mathbf{U})\text{\hspace{1em}(1-4)}$$
其中涉及到的张量运算规则：

$$\tau \cdot \mathbf{U}=\left[\begin{array}{c}{\tau }_{xx}u+{\tau }_{yx}v+{\tau }_{zx}w\\ {\tau }_{xy}u+{\tau }_{yy}v+{\tau }_{zy}w\\ {\tau }_{xz}u+{\tau }_{yz}v+{\tau }_{zz}w\end{array}\right]\text{\hspace{1em}(1-5)}$$
则能量方程可表述为：
$$\frac{D[\rho (e+K)]}{Dt}=\rho r+k\cdot \nabla T+\rho \mathbf{f}\cdot \mathbf{U}-\nabla \cdot (p\mathbf{U})+\nabla \cdot (\tau \cdot \mathbf{U})\text{\hspace{1em}(1-6)}$$

2. 内能方程

上述为能量方程完整形式，CFD 中往往需要从中抽离一个动能项，从中得到一个动能方程。根据动量方程：

$$\frac{D(\rho \mathbf{U})}{Dt}=\nabla \cdot \tau -\nabla p+\rho \mathbf{f}\text{\hspace{1em}(2-1)}$$

将动量方程每个分量乘以速度，加和，可得：

$$\frac{D(\rho K)}{Dt}=(\nabla \cdot \tau )\cdot \mathbf{U}-\nabla p\cdot \mathbf{U}+\rho \mathbf{f}\cdot \mathbf{U}\text{\hspace{1em}(2-2)}$$

则内能方程为：
$$\frac{D(\rho e)}{Dt}=\rho r+k\cdot \nabla T-p\nabla \cdot \mathbf{U}+\tau :\nabla \mathbf{U}\text{\hspace{1em}(2-3)}$$

其守恒形式：

$$\frac{\partial (\rho e)}{\partial t}+\nabla \cdot (\rho \mathbf{U}e)=\rho r+k\cdot \nabla T-p\nabla \cdot \mathbf{U}+\tau :\nabla \mathbf{U}\text{\hspace{1em}(2-4)}$$

3. 温度方程

根据 Fourier 定律，流体的比焓和温度方程可表述为：

$$\frac{\partial (\rho h)}{\partial t}+\nabla \cdot (\rho \mathbf{U}h)=-p\nabla \cdot \mathbf{U}+\nabla \cdot (\lambda \nabla T)+\Phi +S_h\text{\hspace{1em}(3-1)}$$

其中，

• $p\nabla \cdot \mathbf{U}$ ：表面力对流体微元所做的功，一般可忽略；
• $\lambda$：流体的导热系数；
• $S_h$：流体的内热源；
• $\Phi$：耗散函数（由于粘性作用机械能转化为热能的部分）：

$$\begin{array}{rl}\Phi & =2\mu \left[(\frac{\partial u}{\partial x}{)}^2+[(\frac{\partial v}{\partial y}{)}^2+(\frac{\partial w}{\partial z}{)}^2\right]\\ & +\mu \left[(\frac{\partial v}{\partial x}+\frac{\partial u}{\partial y}{)}^2+(\frac{\partial w}{\partial y}+\frac{\partial v}{\partial z}{)}^2+(\frac{\partial u}{\partial z}+\frac{\partial w}{\partial x}{)}^2\right]\\ & +\lambda \nabla \cdot \mathbf{U}\end{array}\text{\hspace{1em}(3-2)}$$

对于理想流体，可取: $h=c_{p}T$，令 $c_p$ 为常数，将耗散函数纳入到源项 $S_T$ 中（$S_T=S_h+\Phi$），有：

$$\frac{\partial (\rho T)}{\partial t}+\nabla \cdot (\rho \mathbf{U}T)=\nabla \cdot (\frac{\lambda }{c_p}\nabla T)+S_T\text{\hspace{1em}(3-3)}$$

对于不可压缩流体，
$$\frac{\partial T}{\partial t}+\nabla \cdot (\mathbf{U}T)=\nabla \cdot (\frac{\lambda }{\rho c_p}\nabla T)+\frac{S_T}{\rho }\text{\hspace{1em}(3-4)}$$