基于混合热格子玻尔兹曼模型的液滴蒸发数值模拟

杨 帆，王亚隆，郭雪岩

（上海理工大学 能源与动力工程学院/上海市动力工程多相流动与传热重点实验室，上海 200093）

相变传热广泛存在于自然界及工业生产过程，其中液滴蒸发是相变传热研究的重要分支。近年来，随着液滴蒸发在农业、医疗、冷却、印刷、镀膜等方面的广泛应用［1-2］，液滴蒸发引起了越来越多学者的关注。国内外学者通过实验观察、理论推导和数值模拟对液滴蒸发进行了一系列的研究，得到了很多有意义的研究结果。

在过去的几十年里，格子玻尔兹曼方法（LBM）已经发展成为流动和传热领域一种非常有效的数值模拟方法［3］。LBM由于具有边界处理简单、并行计算等优点，被广泛应用于模拟单相流、多相流、相变［4-5］等流动和传热现象。目前存在多种多相流LBM模型，包括由Shan等［6］提出的伪势模型，由 Gunstensen 等［7］提出的颜色模型，由Swift等［8］提出的自由能模型，以及由He等［9］提出的相场模型。近年来，越来越多的学者致力于采用格子玻尔兹曼方法模拟蒸发、沸腾、冷凝等气液相变现象。Sun等［10］采用一种三维混合格子玻尔兹曼模型模拟了水平加热壁面上单个气泡和多个气泡的生长等相变现象。Hazi等［11］基于熵的能量方程和非理想气体的热力学关系导出了目标温度方程，提出了一种改进的伪势模型，给出了相应的热格子玻尔兹曼方程，并采用该模型模拟了缓慢流动的流体在水平平板上的沸腾等相变现象。Hatani等［12］采用Lee提出的基于Cahn-Hilliard扩散界面理论的格子玻尔兹曼多相模型模拟了气体的冷凝等相变现象。Gong等［13］基于单弛豫时间的伪势模型提出了一种新的模型用于模拟恒温壁面上和恒热流加热壁面上气泡的产生、长大和脱离壁面的相变现象，而后又模拟了在不同浸润性壁面上气泡的产生等相变现象。Li等［14］采用多弛豫时间的伪势模型和有限差分法的耦合模型模拟了核态沸腾以及液滴的蒸发等相变现象。

在前人研究基础上，本文采用单组分伪势模型与有限差分耦合的混合热格子玻尔兹曼模型(TLBM)，采用真实流体P-R状态方程，对液滴的蒸发过程进行了研究。首先，通过模拟方腔内液滴蒸发过程，验证了所采用方法及程序的有效性。模拟了液滴撞击并附着在高温壁面后的蒸发相变过程，分析了壁面温度、液滴邦德数和雷诺数对液滴撞击壁面后蒸发相变过程的影响。

1 混合热格子玻尔兹曼模型

采用Gong等［13］提出的伪势模型模拟流体的流动。该模型中流体粒子分布函数的演化方程为

式中，F为流体受到的合力。

式（1）中 Δfi(x,t)是通过精确差分方法［15］进行求解，即

本文采用的真实流体P-R状态方程［13］为

流体粒子受到的重力表示为

另外，根据 Ba 等［16］和 Liu 等［17］所采用的方法，计算壁面一层的密度，并给定壁面浸润性。

除了流体流动的模型，本文参考Márkus等［18］使用的能量方程模拟相变过程，即

时间离散化采用四阶龙格库塔法［14］，即

式中，Tt、Tt+δt分别为当前时刻、下一时刻的温度。

2 模拟结果与分析

2.1 液滴蒸发的D2定律

为了验证本文所采用方法及程序的有效性，首先对方形区域内液滴蒸发过程进行模拟，考察液滴直径随时间的变化是否符合D2定律［20］。

本算例中重力为0，导热系数和定压比热容为常数。计算域为方形区域，网格数151 × 151，初始化计算域内液滴直径D0为40个格子单位，液滴温度设置为饱和温度Ts=0.86Tc，周围汽相温度和边界温度均设置为高温Tv=0.90Tc。

图1给出了(D/D0)2随时间的变化，其中D为蒸发过程中液滴直径，D0为液滴初始直径。从图中可以看到，(D/D0)2随时间几乎呈线性变化，这与Nishiwaki等［21］的实验结果相吻合，也印证了本文所采用模型的正确性。

2.2 液滴撞击高温壁面并蒸发

计算域为方形区域，网格数为151 × 151，左右边界设置为周期性边界条件，上边界设置为压力出口，下边界设置为高温壁面。

为了更好地描述撞击高温壁面液滴的物性，引入雷诺数Re和邦德数Bo，其中Re和液滴黏性有关，Re越大，液滴黏性越小。

邦德数Bo和液滴所受重力有关，液滴所受到的重力越大，Bo越大。

图1 液滴撞击高温壁面并蒸发的物理模型Fig.1 Physical model of the droplet collision on the high-temperature wall and its evaporation

液滴在蒸发过程中，若液滴初始体积相同，蒸发完全时间越短，蒸发速率越大。

2.2.1 壁面温度对液滴撞击高温壁面后蒸发过程的影响

图3 液滴撞击不同温度的高温壁面后液滴体积随时间的变化Fig.3 Changes of the droplet volume with time when the droplet collides on the wall with different temperatures

2.2.2 重力对液滴撞击高温壁面后蒸发过程的影响

图4 铺展直径Fig.4 Spreading diameter

图5 不同邦德数的液滴撞击高温壁面后铺展直径、体积随时间的变化Fig.5 Changes of spreading diameter and the droplet volume with time when the droplet with different Bond numbers collides on the high-temperature wall

2.2.3 黏性对液滴撞击高温壁面后蒸发过程的影响

考虑Bo=15.15、温度为 0.86Tc、初始直径相同，Re分别为10.08、12.60、16.80的液滴撞击温度为 0.896Tc的高温壁面，雷诺数越大，液滴的黏性就越小。分析了不同雷诺数的液滴对液滴撞击高温壁面后蒸发过程的影响。图6给出了不同雷诺数的液滴撞击同一高温壁面后、 液滴体积随时间的变化。从图6（a）中可以看到，由于重力相同，所以液滴撞击到壁面所需时间几乎相同。另外，液滴雷诺数越大，液滴黏性越小，液滴撞击到壁面上的最大铺展直径越大。从图6（b）中可以看到，随着雷诺数的增大，液滴撞击到壁面后蒸发速率也增大。这同样是由于雷诺数越大，液滴的黏性越小，液滴撞击到壁面后与高温壁面的接触直径越大，即换热面积增大，换热增强，液滴蒸发得越快。

3 结 论

本文采用混合热格子玻尔兹曼模型模拟了液滴在方形区域内的蒸发过程以及液滴撞击高温壁面并蒸发的过程，得出以下结论：

(1)在方形区域内的蒸发过程中，液滴直径随时间呈线性变化，这与实验结果相吻合。

(2)在液滴撞击高温壁面并蒸发的过程中，壁面温度对液滴撞击壁面后的蒸发速率有显著的影响。壁面温度越高，液滴在撞击壁面后与壁面的换热就越剧烈，液滴在壁面上的蒸发速率就越大。

(3)在液滴撞击高温壁面并蒸发的过程中，重力越大或液滴黏性越小，液滴撞击壁面后的最大铺展直径越大，即液滴与壁面的换热面积就越大，换热就越强，液滴在壁面上的蒸发速率就越大。

图6 不同雷诺数的液滴撞击高温壁面后铺展直径、体积随时间的变化Fig.6 Changes of spreading diameter and the droplet volume with time when the droplets with different Reynolds numbers collides on the high-temperature wall