关于CFD中的能量方程

weixin

总能量方程
能量守恒定义为绝热系统的总能量是一个常数。即总的能量不随时间变化，只能从一种形式转换为另一种形式且不能凭空消失。在CFD里面，通常只考虑动能（机械能）和分子内能（内能）。能量守恒可以表示为：

流体微团内能量的变化率=流入流体微团的净热流量+体积力和表面力对流体微团做功的功率

单位质量的动能定义为：

其中，U 为速度。单位质量的内能定义为：e。那么，流体微团动能和内能总和的时间变化率定义为：

上式即为"流体微团内能量的变化率"。

流入流体微团的净热流量来自于加热、辐射、热传导等。我们定义r 为产生的热源。那么，流体微团的净热源即为：

另外，由于热传导（热流具有方向性）对流体微团的加热为（热通量矢量定义为q）：

进一步的依据傅里叶定律有：

把下两式加和

因此，流入流体微团的净热流量为：

现在我们考虑重力矢量g，对流体微团做功的功率。我们知道，作用在运动物体上的力对物体做功的功率等于这个力乘以速度在这个力上的分量，即：

压力和剪切力对流体微团做功的功率可以表示为：

结合方程

我们有最终的总能量方程：

其中，∇⋅(τ⋅U)=(∇⋅τ)⋅U+τ:∇U∇⋅(τ⋅U)=(∇⋅τ)⋅U+τ:∇U。可以看出，上式采用的是非守恒形式。

内能方程
虽然上式为最终的能量方程。但是在CFD中通常的做法是抽离动能项来获得一个内能方程，内能方程可以从能量方程中减去动能（机械能）方程来获得。首先，我们有动量方程：

将动量方程中的每个速度分量方程乘以速度分量并加和有：

也即：

将总能量方程中减去上式，我们有最终的内能方程：

可以看出，内能方程中不含有体积力项。再一次的，方程(13)采用的是非守恒形式。

焓方程
对于可压缩流体，通常我们把总能量方程简化为焓方程。首先，我们有比焓h 以及总比焓h0 的定义：

回到非守恒形式的方程

其有守恒形式的总能量方程：

把比焓h 方程代入到上式中，我们有最终的守恒形式的比焓方程：

以及守恒形式的总比焓方程：

需要提及的是，在OpenFOAM 的新版本中，传热、可压缩、化学反应求解器求解的主要为守恒形式的总能量方程和比焓方程。并且，以下方程中的∇⋅(τ⋅U)∇⋅(τ⋅U)以及ρg⋅Uρg⋅U 被忽略。

但是在rhoCentralFoam中，则植入了此项。OpenFOAM 官方表示，能量方程的选择在某些情况下是至关紧要的。比如在进行震波捕获计算的时候，求解总能量方程要比求解内能方程结果精确的多。OpenFOAM 中还进一步的假设q=−αeff∇eq=−αeff∇e 或q=−αeff∇hq=−αeff∇h，这样以上两式则简化为（若无热源ρr=0）：

Boussinesq假定
在传热领域内，Boussinesq 假定认为在流动中温度的变化是非常小的，因此密度的变化也非常小。所以在流动的控制方程中，除了浮力项ρg，其他项的密度可以认为是常数。Boussinesq 假定使得方程的非线性特性降低，但由于Boussinesq 假定的限定条件，也使得调用Boussinesq 假定的求解器存在一定的限制。在工程中，Boussinesq 假定主要用于室温下的液体对流、建筑物对流、气相分散等。在温度变化比较大的情况下，不建议使用Boussinesq 假定。

Boussinesq 假定认为流体的密度可以这样计算：

其中，ρ 表示流体的密度，ρref 表示流体的参考密度，T 表示流体的问题，Tref 表示流体的参考问题，β 表示流体的体膨胀系数。考虑可压缩并附加重力的NS方程有：

依据Boussinesq 假定，我们认为浮力项外的密度为常数，因此提出ρ 有：

对于连续性方程，同样的，如果认为密度为常数，因此有连续性方程：

把连续性方程带入到动量方程中

有：

将流体的密度计算方程代入到上式中：

令，

有：

上式即为考虑Boussinesq 假定的动量方程，且可以看出调用Boussinesq 假定的动量方程为不可压缩的。

参考文献：
[1] 安德森, 吴颂平, 刘赵淼. 计算流体力学基础及其应用[M]. 机械工业出版社, 2007.
[2] Versteeg H.K, Malalasekera W. An introduction to computational fluid dynamics: the finite volume method[M]. Pearson Education, 2007.
[3] Energy Equation in OpenFOAM. http://cfd.direct/openfoam/energy-equation

来源：工程事微信公众号（ID：gctechshow）