超声波作用对液滴冷却冻结过程的影响

高蓬辉, 丛 健, 张东海, 周晋鹏, 张正函

(1. 中国矿业大学 力学与土木工程学院, 江苏 徐州 221116;2. 中国矿业大学 深部岩土力学与地下工程国家重点实验室, 江苏 徐州 221116;3. 中国矿业大学 国际学院, 江苏 徐州 221116)

1 前 言

超声波辅助冷却冻结在食品冻结保鲜[1]、医用治疗[2]、生物化学[3-4]等领域得到广泛的应用，具有较好的应用前景。如何利用超声波实现高效冻结成为近年来广泛关注的问题之一。Wang 等[5]通过基于熵方法的纯物质冻结相场模型追踪冻结过程中固-液交界面的位置移动。赵凯璇等[6]建立了真空环境中单液滴闪蒸冻结模型，分析了环境压力和温度对于液滴温度和半径变化的影响。耿县如等[7]讨论了液滴大小、气体流速、环境温度对液滴温度和凝固时间的影响。吕明等[8]结合Rayleigh-Plesset 方程对液滴内空化气泡生长机理和过程进行了模拟分析。Virone 等[9]设计了一种超声波空穴反应器，给出了空化效应所产生气泡的膨胀、破裂方程以及溶液结核速度的计算公式，分析了溶液温度和结核速度的变化规律。Mathieu等[10]研究了在超声波作用下水冻结产生的冰核数量与液体温度和超声波压力间的关系，建立了超声波空化作用下单个气泡的Navier-Stokes equations 方程(简称N-S 方程)，获得了实际液滴空化效应产生的冰核数。余德洋[11]利用显微成像技术分析了纯水、蔗糖溶液在超声波作用下的结晶以及超声波对冰晶二次成核和冰晶生长速度的影响。

已有文献探讨了超声波强化冻结成核和冰晶的生长机理，但对于超声波作用下冻结过程中热质传递机理研究较少，尚未阐明超声波空化效应和热效应在冷却冻结过程中的作用机理。本文将在声场理论的基础上，结合液滴冷却冻结热质传递规律，建立超声波作用下液滴冻结的热质传递模型，研究超声波作用下液滴冻结过程中液滴温度、液固界面、直径大小的变化规律，确定超声波空化效应导致传质引起的热传递量与超声波热效应产热量的变化关系，为超声波辅助冻结的高效利用提供参考和指导。

2 液滴冷却冻结模型

2.1 物理模型及实验装置

超声波分为2 大类：一类为高频率、低能量诊断超声波；另一类为强化液体结晶的低频率、高能量超声波[12]。超声波空化效应[13]是液体中微气泡在声场作用下快速膨胀、闭合、崩溃，强化了液滴表面与周围环境间的传质，有助于降低界面能[14]和提高液滴的成核速度[15]。超声波作用下液滴冻结模型如图1 所示。

图1 超声波作用下液滴冷却冻结示意图Fig.1 Schematic diagram of droplet cooling and freezing under ultrasonic

超声波作用下液滴冻结实验装置见图2，将液滴悬挂于低温箱中，通过容器底部布置的超声波发生器输入超声波，利用显微镜对液滴冻结进行观测。实验中采用的设备为金相显微镜SGO-BH200M (物镜放大倍数为5～50 倍)、高速摄像机 (拍摄帧数15 000 fps)、热电偶 (测量范围-40～150 )℃、超声波发生器 (超声波频率20～100 kHz，强度0～2 000 W·m-2)以及恒温恒湿箱 (初始温度设置为-15 ℃)，液滴的初始直径2.0 mm。

图2 超声波作用下液滴冻结实验系统Fig.2 Experiment setup of droplet freezing under ultrasonic

2.2 数学模型

对于相变冻结过程，Feuillebois 等[16]、张跃凯[17]根据冻结条件的差异将冻结过程进行了不同的划分。结合超声波作用下液滴的实际冷却冻结过程，将液滴的冷却冻结分为2 个阶段：1) 蒸发冷却过程，液滴蒸发冷却至凝固点，液滴表面状态从液态转变为固态；2) 液滴深度冻结过程，从表面向内部进行冻结，液滴由液固状态逐渐转变为固态。

液滴表面蒸发过程中的传质为

将水蒸气视为理想气体，根据理想气体状态方程，由式(1)得

超声波作用下液滴在蒸发过程中的热量交换包括超声波热效应的产热量Qheat、液滴表面对流换热量Qconv、液滴表面传质引起的传热量Qevap、液滴内部的能量变化为Qin，其能量方程为

液滴内部的能量变化：

液滴表面对流换热量：

式(5)中液滴与周围环境的对流传热系数[18]：

液滴表面传质引起的传热量：

式(7)中液滴表面的对流传质系数[19]：

当液滴表面冻结成固体时，其表面与周围气流之间的有效扩散系数为

当液滴处于固液共存状态时，液滴中液体所占的比例为κ 时，液滴的总质扩散系数DGLS为

超声波热效应产热量：

超声波作用在液滴上总的衰减系数由黏滞衰减系数αη、热传导衰减系数αξ和散射衰减系数αs组成，各部分的计算式如式(12)～(14)：

考虑液滴孔隙率的影响，液滴对超声波总的吸收系数为[20]

根据能量守恒，能量变化量等于进入液滴的能量和离开液滴的能量的差值，由式(4)、(5)、(7)和(11)得到液滴冷却能量方程：

对于液滴的深度冻结阶段，考虑液滴表面已全部由液态转变为固态。液滴的凝固释放潜热Qsolid，冻结过程应满足能量平衡关系式

液滴的凝固放热量Qsolid为

同理得到液滴冻结过程中能量方程：

3 结果与讨论

图3 超声波作用和无超声波作用液滴冻结过程中直径随时间的变化规律Fig.3 Profiles of droplet diameter in freezing processes with and without ultrasonic

利用超声波作用下液滴冷却冻结模型分析不同的初始液滴直径、超声波强度及频率下，液滴的直径、温度、固液交界面位置随时间的变化规律，探讨超声波热效应产热量和超声波作用下液滴表面传质引起传热量的变化规律，确定超声波作用时间对冻结的影响。模拟计算过程中主要参数见表1。

表1 主要参数Table 1 Parameters for model simulation

图4 不同大小液滴直径随冻结时间的变化规律Fig.4 Profiles of droplet diameter as a function of freezing time

3.1 液滴物理状态

3.1.1 液滴直径超声波和无超声波作用下液滴冷却冻结过程中直径的变化规律见图3。超声波强化了液滴与周围环境的质传递，因此超声波作用下液滴直径随时间变化剧烈，减小速度较快。超声波声强周期性变化，导致液滴和环境之间的压力差会出现正负交替，从而使液滴直径出现一定的往复变化。当液滴进入凝固升华阶段时，水的体积会膨胀，所以液滴直径先增大后减小，最终趋于稳定。可见，超声波可有效强化液滴的冷却冻结，未加超声波时液滴的冻结最终直径为0.993 8 mm，液滴的直径减少了0.62%，加载超声波的液滴最终直径为0.977 3 mm，减少了2.27%。

在不同的液滴初始直径2.0、1.6 和1.0 mm 下，液滴直径随冻结时间变化规律如图4 所示。从图4可以看出，当液滴完全冻结时液滴直径分别减少0.12、0.04 和0.022 7 mm。液滴初始直径越大，最终直径减少越多。原因在于液滴的初始直径越大，液滴与环境相接触的表面积越大，液滴表面传质越强。

3.1.2 液滴温度

超声波和无超声波作用下液滴冷却冻结过程中温度变化规律见图5。从图5 中可以看出在超声波作用下，液滴的温度下降较快（超声波作用下液滴的蒸发阶段时间比无超声波作用的时间缩短1/2 左右）。冷却冻结相同时间，无超声波作用时液滴的最低温度可以达270.9 K；在超声波作用下液滴的最低温度可以达到264.9 K。因为超声波空化效应增强了液滴与周围环境间的传质，导致传质引起的传热量增加，从而促进了液滴的冷却冻结，这一结果与文献[21-22]的研究结论一致。

图5 超声波作用和无超声作用波液滴冻结过程中温度随时间的变化规律Fig.5 Profiles of droplet temperature in freezing processes with and without ultrasonic

图6 超声波作用下液滴冻结过程中液滴直径D和温度Td 随时间的变化规律Fig.6 Profiles of droplet diameter D and temperature Td in freezing processes with ultrasonic

液滴冻结过程中液滴直径与液滴温度随时间的变化规律如图6 所示。分析过程中液滴过冷度设为5 K，在到达液滴凝固点前，液滴快速蒸发，温度不断降低，直径也逐渐减小。从图6 中可以看出液滴在冻结过程中温度变化与液滴直径的变化趋势基本一致，当液滴经过凝固点后，液滴膨胀直径增大，同时液滴凝固释放的热量会导致液滴的温度出现短暂上升随后下降。

3.2 液固界面

当超声波的频率为40 000 Hz，强度分别为1 000、700、400 W·m-2，液滴初始半径r=1.0 mm 时，冻结过程中液固界面位置的变化见图7。图7 表明超声波强度较大时液滴冷却冻结过程中热质传递较强，液固交界面的移动速度较快，液滴由外到内完全冻结所需要时间较短，曲线的斜率越大。当超声波强度分别为1 000、700 和400 W·m-2时，液滴液固交界面移动到液滴中心的时间分别为46、48 和54 s。

图7 不同超声波强度下液滴冻结过程中液固交界面位置随时间的变化规律Fig.7 Profiles of droplet solid-liquid interface position in freezing processes under different ultrasonic intensities

图8 不同大小液滴冻结过程中液固交界面随时间的变化规律Fig.8 Profiles of droplet solid-liquid interface position in freezing processes for different droplet radius

不同大小液滴在冷却冻结过程中液固交界面位置随时间的变化规律见图8，根据曲线的斜率可以得到液滴的固-液交界面位置随时间的变化速度。液滴的直径越小，其液固交界面的移动速度越快，液滴冻结所需的时间越短。因为在相同的超声波强度、频率和外界条件下，液滴的直径越小，液滴内部单位体积所产生气泡的数量越多，液滴表面与周围环境的传质过程越剧烈，使液滴冷却冻结得到强化。

3.3 空化效应引起的传热量和热效应的产热量

超声波辅助冷却冻结过程中，超声波空化效应增强了液滴与周围环境之间的传质和传热，而超声波热效应会产生热量（相当于增加了内热源）。蒸发冻结前期，超声波空化效应导致传质引起的传热量较大，超声波有利于液滴的快速蒸发冻结。当液滴表面开始冻结，液态水所占比例下降，超声波空化效应导致传质引起的传热逐渐减弱，在液滴冻结的某个时刻，超声波热效应的产热量将大于超声波空化效应传质引起的传热量，此后加载超声波不能强化液滴的冷却冻结。

不同超声波频率下，超声波空化效应导致传质引起的传热量和超声波热效应产热量的变化规律见图9。图中虚线表示不同超声波频率下超声波热效应产热量，实线表示超声波空化效应导致传质引起的传热量。在初始阶段，空化效应有助于液滴表面的传质，强化了液滴表面的传热。随着冻结的进行，液滴表面逐渐由液态转变为固态，液滴外表面的传质减弱，导致传质引起传热量增加的幅度逐渐减小，之后出现下降；而超声波热效应的产热量随时间逐渐增加。当到达临界点C1、C2、C3、C4、C5和C6即超声波热效应产热量与空化效应传质引起的传热量相等的位置时，超声波热效应的产热量与空化效应传质引起的传热量相等，此后加载超声波不能起到强化冷却冻结的效果。

图9 液滴在不同超声波频率下传质的传热量Qevap和超声波热效应产热量Qheat 的变化规律Fig.9 Profiles of heat by mass transfer Qevap and heat by ultrasonic thermal effect Qheat under different ultrasonic frequencies

图10 液滴在不同超声波强度下的传质传热量Qevap和超声波热效应产热量Qheat 的变化规律Fig.10 Profiles of heat by mass transfer Qevap and heat by ultrasonic thermal effect Qheat under different ultrasonic intensities

不同超声波强度下超声波空化效应导致传质引起的传热量和超声波热效应产热量的变化见图10。图中虚线表示超声波热效应的产热量，实线表示超声波空化效应导致传质引起的传热量。液滴与周围环境间传质导致的传热量曲线与超声波热效应产热量曲线的交点(临界点)C1、C2、C3和C4，随着超声波强度的降低向后移动。同样临界点之后液滴内超声波热效应的产热量将大于超声波空化效应传质引起的传热量，此后加载超声波不利于液滴的冻结。

3.4 模型验证

图11 液滴冻结过程中固液交界面随时间的变化规律Fig.11 Profiles of droplet solid-liquid interface as a function of time

在液滴冷却冻结实验的基础上，将超声波和无超声波作用下液滴冻结过程中固液界面的实验结果与模拟结果进行了对比，见图11。实线表示超声波作用下液滴液固界面位置随时间变化的模拟曲线。实验中在-15 ℃的环境中，超声波作用下液滴完全冻结需要52 s，无超声波作用需要65 s，可以发现超声波作用下液滴更容易冻结。对比实验与模拟结果，两者最大误差10.6%，最小误差0.3%，且两者变化趋势相同，表明该理论模型能够较好地模拟预测超声波作用下液滴的冷却冻结。

4 结 论

在声场理论和热质守恒的基础上，建立了超声波作用下液滴冷却冻结数学模型，分析了不同条件下超声波对冻结过程的影响，给出了在液滴冷却冻结过程中超声波空化效应导致传质引起传热量与超声波热效应产热量的变化规律。主要结论为：

(1) 获得了超声波和无超声波作用下液滴冷却冻结过程中直径、温度随时间的变化规律，明确了超声波作用有助于液滴的蒸发冷却。

(2) 超声波空化效应强化了液滴与周围环境的传质，有助于液滴的冷却冻结。当超声波的强度为1 000、700 和400 W·m-2时，完全冻结所需时间分别为46、48、54 s。

(3) 揭示了超声波热效应产热量与超声波空化效应导致传质引起传热量的变化规律，确定了超声波辅助冷却冻结的作用时间临界点，为超声波冷却冻结的高效利用提供了依据和参考。

符号说明：

c — 超声波在液滴中的传播速度，m·s-1

cp— 液滴的比定压热容，J·kg-1·K-1

cV— 液滴的比定容热容，J·kg-1·K-1

D — 液滴的直径，m

Dv— 溶液与空气间的质扩散系数，m2·s-1

DGA— 冰体与空气间的质扩散系数，m2·s-1

DGS— 冰体与空气的有效质扩散系数，m2·s-1

DGLS— 液滴表面总的质扩散系数，m2·s-1

f — 超声波频率，Hz

h — 对流传热系数，W·m-2·K-1

Hm— 传质换热系数，m·s-1

I — 超声波强度，W·m-2

K — 孔隙率

Levap— 液滴的蒸发潜热，J·kg-1

m˙ — 质扩散速率，kg·m-2

M — 水蒸气的摩尔质量，kg·mol-1

Lsolid— 液滴的凝固潜热，J·kg-1

N — 液体界面处溢出的气泡数

pd— 液滴的压力，Pa

p∞— 环境压力，Pa

Pr — 普朗特数

r — 液滴固液界面的位置，mm

u — 气流速度，m·s-1

R — 通用气体常数，J·mol-1·K-1

R0— 液滴的初始直径，m

Rd— 液滴的半径，m

Rg— 气体常数，J·mol-1·K-1

Re — 雷诺数

Sn— 总的液滴表面更新面积，m2

Sc — 施密特数

t — 液滴的冻结时间，s

Td— 液滴的温度，K

T∞— 环境温度，K

Vd— 液滴的体积，m3

αs— 散射衰减系数，m-1

αξ— 热传导衰减系数，m-1

αη— 黏滞衰减系数，m-1

β — 水对超声波的吸收系数，m-1

γ — 比热容比

ε — 液滴的空隙率

φ — 液滴表面更新比例

η — 液滴的动力粘度，Pa·s

κ — 液滴的液体部分所占比例

λ — 热传导系数，W·m-1·K-1

ρ — 液滴密度，kg·m-3

ρvapour— 液滴表面的水蒸气质量浓度，kg·m-3

τ — 曲折因数