连续介质热力学的能量守恒定律

力学有四大基本定律（“基本公理”，Eringen, 1980, P89-92）：质量守恒原理、动量平衡原理、动量矩平衡原理、能量守恒原理。本节讨论连续介质的能量守恒定律，从全局形式出发，通过对动能、内能、热能和外力功的表述，得到局部形式的能量守恒定律。

1. 全局能量守恒定律

动能与内能之和的变化率等于外力的功率与单位时间内进入或流出系统的所有其他能量之和：
$$\frac{\mathrm{D}}{\mathrm{D}t}(\mathcal{K}+\mathcal{E})=\mathcal{W}+\sum _i{\mathcal{U}}_i(\mathrm{4.2.1})$$

式中$\mathcal{K}$、$\mathcal{E}$、$\mathcal{W}$分别为动能、内能和功率。${\mathcal{U}}_i$表示进入或流出系统的第$i$类能量（例如电能、化学能、热能）。这就是能量守恒的全局形式。

式(4.2.1)还可写为：

$$\dot{\mathcal{K}}+\dot{\mathcal{E}}=\mathcal{W}+Q(\mathrm{4.2.2})$$

2. 热力学第一定律

为简单起见，我们只考虑具有热能$Q$的系统，并考虑热静力学（$\dot{\mathcal{K}}=0$），于是式(4.2.1)写为：
$$\dot{\mathcal{E}}=\mathcal{W}+Q$$

令$\mathrm{d}\mathcal{E}=\dot{\mathcal{E}}\mathrm{d}t$、$\delta \mathcal{W}=\mathcal{W}\mathrm{d}t$、$\delta \mathcal{Q}=Q\mathrm{d}t$，得到热力学第一定律：
$$\mathrm{d}\mathcal{E}=\delta \mathcal{W}+\delta \mathcal{Q}$$

3. 连续介质的局部能量守恒定律

动能为：
$$\mathcal{K}=\frac{1}{2}{\int }_V\rho v_k{v}_k\mathrm{d}v(\mathrm{4.2.3})$$

对于连续介质我们假设存在一个单位质量的内能密度$\epsilon$：
$$\mathcal{E}={\int }_V\rho \epsilon \mathrm{d}v(\mathrm{4.2.4})$$

计算动能和内能的物质导数：
$$\dot{\mathcal{K}}={\oint }_V\left[\rho v_k{a}_k\mathrm{d}v+\frac{1}{2}{v}_k{v}_k\frac{\mathrm{D}}{\mathrm{D}t}(\rho \mathrm{d}v)\right](\mathrm{4.2.}{3}^{\prime })$$

$$\dot{\mathcal{E}}={\int }_V\left[\rho \dot{\epsilon }\mathrm{d}v+\epsilon \frac{\mathrm{D}}{\mathrm{D}t}(\rho \mathrm{d}v)\right](\mathrm{4.2.}{4}^{\prime })$$

在连续介质$V+S$中，热量$Q$可通过系统表面$\mathcal{S}$进入系统，也可以通过系统内部的单位质量的分布热源 $h$ 产生。令$\mathbf{q}$为系统表面处的热量矢量并指向系统外部，在该点$\mathbf{x}$处的外法线方向为$\mathbf{n}$，于是可写出$Q$的表达式：

$$Q={\oint }_S\mathbf{q}\cdot \mathbf{n}\mathrm{d}a+{\int }_V\rho h\mathrm{d}v(\mathrm{4.2.3})$$

另一方面，外力的功率$\mathcal{W}$ 由表面力和体积力的功率之和组成，考虑非极性介质（Eringen, 1980, 中译本P121）：

$$\mathcal{W}={\oint }_V\rho \mathbf{f}\cdot \mathbf{v}\mathrm{d}v+{\int }_S{\mathbf{t}}_{(\mathrm{n})}\cdot \mathbf{n}\mathrm{d}a={\oint }_V\rho f_k{v}_k\mathrm{d}v+{\int }_S{t}_{lk}{v}_k{n}_l\mathrm{d}a(\mathrm{4.2.5})$$

式中表面力的计算用到了关系式${\mathbf{t}}_{(\mathrm{n})}=t_{kl}{n}_k{\mathbf{i}}_l$（Eringen，1980，P105）。

可以通过Gauss定理将外力功$\mathcal{W}$和热量$Q$的面积分转换为体积分：

$$\mathcal{W}={\int }_V\left(t_{lk,l}{v}_k+t_{lk}{v}_{k,l}+\rho f_k{v}_k\right)\mathrm{d}v(\mathrm{4.2.6})$$

$$Q={\int }_V\left(q_{k,k}+\rho h\right)\mathrm{d}v(\mathrm{4.2.}{6}^{\prime })$$

将式(4.2.3’)(4.2.4’)(4.2.6)(4.2.6’)代入(4.2.2)，并使用局部动量平衡原理和局部能量守恒原理，整理得：

$$\rho \dot{\epsilon }=t_{kl}{v}_{l,k}+q_{k,k}+\rho h(\mathrm{4.2.8})$$

根据相关讨论（Eringen，1980，P84-84），可知速度梯度可分解为变形率张量和自旋张量，即$v_{k,l}=d_{kl}+{\omega }_{kl}$，可证明$t_{kl}{v}_{l,k}=t_{kl}{d}_{l,k}$，于是(4.2.8)还可写为：

$$\rho \dot{\epsilon }=t_{kl}{d}_{l,k}+q_{k,k}+\rho h(\mathrm{4.2.9})$$

机械能$t_{kl}{d}_{l,k}$也称为应力功率。(4.2.9)表明单位时间内的内能变化是由应力功率和热量（包括$q$和$h$）产生的。这就是连续介质的局部能量平衡方程。