蒸发模型

分类:流体分析(CFD) | 综合版 2026-04-06
CAE visualization for evaporation model theory - technical simulation diagram
蒸发模型

蒸发的理论基础

概述

🧑‍🎓

老师,CFD的蒸发模型计算什么?


🎓

液滴或液膜的蒸发速度预测模型。燃料喷射后的液滴蒸发、喷雾干燥、冷却塔、涂装干燥等,处理伴随相变化的热物质传递。拉格朗日粒子追踪(DPM)中各液滴蒸发计算是主流。


🧑‍🎓

液滴的蒸发是什么物理现象?


🎓

液滴表面的蒸汽向周围气体扩散的过程。液滴表面处于饱和浓度,对应蒸气压,远处浓度较低,浓度差是驱动力。同时蒸发潜热吸收导致液滴温度下降,热传递与物质传递耦合的问题。


支配方程

🧑‍🎓

请告诉我蒸发速度的公式。


🎓

古典的Abramzon-Sirignano(1989)模型中,液滴的质量变化率表示为。


$$ \dot{m} = \pi d_p \rho_g D_{AB} Sh^* \ln(1 + B_M) $$

🎓

这里 $d_p$ 是液滴直径,$D_{AB}$ 是蒸汽的二成分扩散系数,$Sh^*$ 是修正Sherwood数,$B_M$ 是Spalding质量传递数。


$$ B_M = \frac{Y_s - Y_\infty}{1 - Y_s} $$

🧑‍🎓

$Y_s$ 是液滴表面的蒸汽质量分率吧?


🎓

完全正确。$Y_s$ 由液滴温度的饱和蒸气压通过Clausius-Clapeyron式求得,由摩尔分率转换为质量分率。$Y_\infty$ 是远处(CFD单元平均)的蒸汽质量分率。


🎓

液滴温度变化由热平衡求得。


$$ m_p c_p \frac{dT_p}{dt} = \pi d_p k_g Nu^* \ln(1 + B_T) \frac{(T_\infty - T_p)}{B_T} - \dot{m} h_{fg} $$

🎓

$B_T$ 是Spalding热传递数,$Nu^*$ 是修正Nusselt数。右边第1项是对流加热,第2项是蒸发冷却(潜热吸收)。


d-二乘法则

🧑‍🎓

我听说过d-二乘法则。


🎓

稳态蒸发时液滴直径的平方随时间线性减小。


$$ d^2(t) = d_0^2 - K t $$

🎓

$K$ 是蒸发速率常数,由液体类型和环境条件决定。这种线性减小就是d-二乘法则,是蒸发模型验证最常用的指标。


Coffee Break 杂谈

D2法则——液滴蒸发的"黄金律"与现实的偏差

液滴蒸发速度由"D2法则"控制,液滴直径的平方随时间线性减小。D^2 = D0^2 - K*t(K:蒸发常数)由1953年Godsave和Spalding几乎同时发表的燃料液滴燃烧解析而来。然而现实液滴受内部循环、Marangoni对流、气相热辐射等影响,D2法则可能产生最大30%的误差。多液滴密集的喷雾环境中,相邻液滴的蒸汽作用导致Group evaporation现象,偏差更大。

蒸发的数值计算手法

数值求解的细节

🧑‍🎓

请告诉我蒸发模型的数值计算要点。


🎓

在拉格朗日液滴追踪中,每个时间步对每个液滴计算蒸发量。蒸发失去的质量作为源项反映到气相的物种输送方程(two-way coupling)。


🎓

计算流程如下。

1. 在液滴位置补间气体温度和蒸汽浓度

2. 计算液滴表面的饱和蒸气压(Antoine式或Clausius-Clapeyron式)

3. 计算Spalding传递数 $B_M$, $B_T$

4. 计算蒸发速率 $\dot{m}$ 和液滴温度变化率

5. 更新液滴质量和直径

6. 将气相的质量、动量、能量源反映到气相


多成分液滴的蒸发

🧑‍🎓

汽油这样混合多种成分的液滴怎么处理?


🎓

多成分蒸发模型中,各成分的蒸气压由Raoult定律求得。


$$ p_i = x_i \gamma_i p_i^{sat}(T) $$

🎓

$x_i$ 是液相中的摩尔分率,$\gamma_i$ 是活度系数,$p_i^{sat}$ 是纯成分的饱和蒸气压。轻质成分先蒸发,液滴组成随时间变化。


🧑‍🎓

考虑液滴内部的温度分布吗?


🎓

简易模型(Uniform Temperature)假设液滴内部温度均匀。精密模型(Diffusion Limit)求解液滴内部温度和组成的径向分布。Fluent可选"Infinite Diffusion"或"Diffusion-Limited"。


工具别的实现

工具蒸发模型多成分特点
Ansys FluentConvection/Diffusion ControlledRaoult法则物种输送耦合
STAR-CCM+Abramzon-Sirignano多成分支持拉格朗日框架
OpenFOAM (sprayFoam)多种蒸发模型支持可定制
CONVERGE多成分蒸发详细化学耦合AMR支持
🎓

燃烧分析时蒸发模型和化学反应模型的耦合很重要。液滴蒸发产生的蒸汽着火、燃烧过程的精确复现,需要物种输送的分辨率和时间步管理。


Coffee Break 杂谈

Langmuir-Knudsen模型——薄膜蒸发的分子论基础

液体表面分子的蒸发和凝聚速度由Langmuir-Knudsen模型从分子运动论描述,涉及Knudsen层的物理,其中气液界面的平均自由程与膜厚相当。在微设备(MEMS)薄膜蒸发器和纳米液滴蒸发中,连续体假设破裂,需要该模型。然而CFD实现时"蒸发系数σ_e(0.01~1范围内约4位数的不确定性)"的不确定性是最大的课题,相同水蒸发在文献值下计算结果差异巨大。

蒸发的实务应用

实践指南

🧑‍🎓

请说明带蒸发的喷雾分析步骤。


🎓

以喷雾干燥机分析为例。


🎓

1. 建模:干燥室、喷霾口位置、排气口

2. 网格生成:喷口附近密集(1~2 mm),远处粗糙(5~10 mm)

3. 气相设置:热风温度、湿度、流速作为入口条件

4. 液滴投入:从喷口投入Rosin-Rammler分布液滴

5. 蒸发模型:启用Convection/Diffusion Controlled蒸发

6. 物种输送:添加水蒸汽的输送方程

7. 后处理:液滴轨迹、蒸发完成位置、出口温湿度


物性值的重要性

🧑‍🎓

蒸发分析中特别重要的物性是什么?


🎓

以下物性对蒸发速率有直接影响。


物性重要度说明
饱和蒸气压 $p_{sat}(T)$最高Antoine式参数精度是关键
蒸发潜热 $h_{fg}$考虑温度依赖性
二成分扩散系数 $D_{AB}$Chapman-Enskog理论或实测
液体比热 $c_{p,l}$影响液滴温度变化
液体密度 $\rho_l$影响液滴直径计算
🧑‍🎓

Antoine式的参数在哪里能获得?


🎓

NIST Chemistry WebBook是最可信的数据源。DIPPR(Design Institute for Physical Properties)数据库也是业界标准。


常见验证方法

🧑‍🎓

如何验证蒸发模型?


🎓

常用以下基准实验。


实验条件测量量
单液滴蒸发悬滴法d-二乘法则的蒸发速率常数 $K$
Ranz-Marshall对流中单液滴Sherwood数、Nusselt数
ECN Spray An-十二烷高温环境喷射蒸汽渗透长、液相到达距离
喷雾干燥实机数据出口温度、残留水分
🎓

先在单液滴确认d-二乘法则,再在喷雾系中验证全体行为是标准方法。


Coffee Break 杂谈

衣料干燥的CFD——干衣机制造商的节能设计革命

家用衣料干衣机开发中CFD的本格活用是2010年代开始的。模拟滚筒内温湿度分布、空气流动、布料含水率变化的耦合模型,通过优化滚筒形状和通气孔位置,实现消费电力15%削减。难点是"布多孔质内毛细管蒸发"和"滚筒强制对流蒸发"两阶段的模型化。欧洲能源标签规制严格化,这种精细模拟成为家电制造商竞争力源泉。

蒸发的软件比较

商用工具比较

🧑‍🎓

请比较实现蒸发模型的工具。


🎓
工具单成分多成分Flash蒸发特点
Ansys FluentConvection/DiffusionRaoult, UNIFACThermal breakupVOF-to-DPM耦合
STAR-CCM+Abramzon-Sirignano多成分Flash boiling拉格朗日/欧拉交换
OpenFOAMsprayFoam内置支持有限可自由添加模型
CONVERGE详细蒸发多成分支持AMR + 燃烧耦合

用途别推荐

🧑‍🎓

根据用途选择吗?


🎓
用途推荐工具理由
柴油喷射+蒸发+燃烧CONVERGE, FluentAMR支持、详细化学耦合
GDI(汽油直喷)Flash boilingSTAR-CCM+, FluentFlash蒸发模型装备
喷雾干燥Fluent, STAR-CCM+DPM + 蒸发 + 壁面碰撞
冷却塔FluentWet-Bulb模型集成
学术研究OpenFOAM新蒸发模型实现
🧑‍🎓

Flash boiling是什么?


🎓

高温高压燃料从喷口喷射时,急速减压使液温超过饱和温度,液滴内部急速沸腾。GDI引擎高温条件下发生,促进喷雾细化。与常规表面蒸发机理不同,需要专用Flash boiling模型。


🎓

Fluent 2021以后装备Thermal Breakup模型,STAR-CCM+装备Flash Boiling Spray Model。


Coffee Break 杂谈

工艺模拟器与CFD——精馏塔设计的双刀流

化工厂精馏塔设计中,Aspen Plus等工艺模拟器计算全体物质能量平衡,CFD分析托盘或填料的详细行为,"分层方法"已成定式。大型化学制造商(BASF、Dow等)用CFD分析托盘液膜蒸发与实际分离效率的对应关系,建立数据库,成为下一代精馏塔设计的差异化来源。Siemens Star-CCM+的欧拉多相+蒸发模型在精馏托盘基准事例丰富,评价实绩直接成为客户的信任基础。

蒸发的先端研究

先端技术和研究动向

🧑‍🎓

蒸发模型的最新研究有什么?


🎓

来看几个方向。


超临界蒸发

🎓

火箭引擎或柴油引擎高温高压条件下,液滴周围气体超过燃料的临界点。此条件下液气界面消失,常规蒸发模型(d-二乘法则)不适用。采用Peng-Robinson或Soave-Redlich-Kwong状态方程的Real-Fluid蒸发模型在研究中。


🧑‍🎓

常规蒸发模型不能用了。


🎓

Matagni & Bellan(MIT)或Okong'o & Bellan的DNS研究是先驱,超临界液滴的"蒸发"表述为界面无密度过渡。CONVERGE或OpenFOAM该类模型实现研究进行中。


多成分蒸发的详细化学

🎓

汽油等数百成分混合燃料定义代理燃料(5~20成分替代混合物),追踪各成分蒸发。与NIST/LLNL化学反应机制的一致性很重要。


DNS液滴蒸发的直接计算

🎓

DNS直接计算液滴周围速度场、温度场、浓度场,用于现有模型验证和闭合改进的研究增加。Schlottke & Weigand(Stuttgart)的VOF-DNS代表性。


机器学习的活用

🧑‍🎓

用AI改进蒸发模型的研究也有吗?


🎓

以DNS数据为教师数据,用神经网络构建液滴蒸发速度的代理模型研究有。多成分蒸发中组成-温度-压力高维参数空间的表查询,用学习模型效率化的尝试也有。


Coffee Break 杂谈

超临界蒸发——液体火箭引擎的"另一个世界"

液氢、液氧为推进剂的液体火箭引擎中,燃烧室压力超过推进剂临界压力的超临界状态下发生喷射、混合、燃烧。此条件液滴蒸发概念消失,液体到气体的"连续过渡"支配,古典蒸发模型根本失效。采用Peng-Robinson等实在气体状态方程与CFD组合的超临界蒸发分析是JAXA、NASA从2000年代精力投入的领域。超过Widom线时物性急变(密度变化数十%)使CFD收敛困难,适应性状态方程自动切换是实装关键。

蒸发的故障处理

故障排除

🧑‍🎓

请告诉蒸发模型常见故障。


🎓

逐一看。


1. 液滴不蒸发

🎓

症状:液滴不蒸发,直接飞过。


🎓

对策

  • 确认物种输送启用,蒸汽成分定义正确
  • 确认液滴的沸点和潜热设置正确
  • 确认环境温度高于液滴沸点(低温下蒸发极慢)
  • 确认DPM的蒸发模型启用

2. 液滴温度非物理性低

🧑‍🎓

液滴温度变成负数了…


🎓

对策

  • 确认蒸发潜热值正确(注意单位系)
  • 确认液滴尺寸不过小(微小液滴瞬时蒸发温度急降)
  • 确认气相温度场收敛
  • 确认时间步不过大

3. 蒸发速率与实验不符

🎓

对策

  • 确认Antoine式参数(与NIST数据对比)
  • 确认二成分扩散系数 $D_{AB}$ 值和温度依赖性
  • 确认对流条件正确(Reynolds数对应Nu, Sh补正)
  • 多成分液滴情况确认Raoult定律活度系数

4. 气相蒸汽浓度非物理性

🎓

对策

  • 确认2-way coupling正确设置
  • 确认包裹数充分多(少了统计误差大)
  • 确认网格不过粗(源项集中一个单元导致局部过饱和)

5. 工具特有注意

工具注意点
Fluent液滴物种分率初值设置正确(多成分情况和为1)
STAR-CCM+确认Lagrangian phase interaction的mass/heat transfer启用
OpenFOAMsprayFoam的evaporationModel中AntoineCoeffs单位系(Pa或bar)注意
CONVERGE确认自适应网格细化等级在蒸发区域充分
Coffee Break 杂谈

液滴不蒸发——边界条件设置错误的诊断

蒸发CFD中"液滴温度达到沸点但质量不减"的现象,多数情况是液滴周围气体边界条件问题。液滴蒸发速度与周围气体水蒸气分压(湿度)和平衡蒸气压的差成正比,入口边界湿度100%时蒸发驱动力为零。Ranz-Marshall等热质传递相关式在液滴Reynolds数0.1~1000范围有效,超微小液滴(< 5 µm)或高速液滴(> 200 m/s)需要补正系数。发现蒸发不发生时,推荐首先确认各单元蒸发源项的等值线图。

关联模拟器

用这个领域的交互模拟器体会理论

模拟器列表

关联领域

热分析V&V·品质保证结构分析
本文的评价
感谢您的回答!
有参考
价值
想要更
详细
报告
错误
有参考价值
0
想要更详细
0
报告错误
0
Written by NovaSolver Contributors
Anonymous Engineers & AI — 网站地图
查看个人资料