蒸发模型
蒸发的理论基础
概述
老师,CFD的蒸发模型计算什么?
液滴或液膜的蒸发速度预测模型。燃料喷射后的液滴蒸发、喷雾干燥、冷却塔、涂装干燥等,处理伴随相变化的热物质传递。拉格朗日粒子追踪(DPM)中各液滴蒸发计算是主流。
液滴的蒸发是什么物理现象?
液滴表面的蒸汽向周围气体扩散的过程。液滴表面处于饱和浓度,对应蒸气压,远处浓度较低,浓度差是驱动力。同时蒸发潜热吸收导致液滴温度下降,热传递与物质传递耦合的问题。
支配方程
请告诉我蒸发速度的公式。
古典的Abramzon-Sirignano(1989)模型中,液滴的质量变化率表示为。
这里 $d_p$ 是液滴直径,$D_{AB}$ 是蒸汽的二成分扩散系数,$Sh^*$ 是修正Sherwood数,$B_M$ 是Spalding质量传递数。
$Y_s$ 是液滴表面的蒸汽质量分率吧?
完全正确。$Y_s$ 由液滴温度的饱和蒸气压通过Clausius-Clapeyron式求得,由摩尔分率转换为质量分率。$Y_\infty$ 是远处(CFD单元平均)的蒸汽质量分率。
液滴温度变化由热平衡求得。
$B_T$ 是Spalding热传递数,$Nu^*$ 是修正Nusselt数。右边第1项是对流加热,第2项是蒸发冷却(潜热吸收)。
d-二乘法则
我听说过d-二乘法则。
稳态蒸发时液滴直径的平方随时间线性减小。
$K$ 是蒸发速率常数,由液体类型和环境条件决定。这种线性减小就是d-二乘法则,是蒸发模型验证最常用的指标。
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定律求得。
$x_i$ 是液相中的摩尔分率,$\gamma_i$ 是活度系数,$p_i^{sat}$ 是纯成分的饱和蒸气压。轻质成分先蒸发,液滴组成随时间变化。
考虑液滴内部的温度分布吗?
简易模型(Uniform Temperature)假设液滴内部温度均匀。精密模型(Diffusion Limit)求解液滴内部温度和组成的径向分布。Fluent可选"Infinite Diffusion"或"Diffusion-Limited"。
工具别的实现
| 工具 | 蒸发模型 | 多成分 | 特点 |
|---|---|---|---|
| Ansys Fluent | Convection/Diffusion Controlled | Raoult法则 | 物种输送耦合 |
| STAR-CCM+ | Abramzon-Sirignano | 多成分支持 | 拉格朗日框架 |
| OpenFOAM (sprayFoam) | 多种蒸发模型 | 支持 | 可定制 |
| CONVERGE | 多成分蒸发 | 详细化学耦合 | AMR支持 |
燃烧分析时蒸发模型和化学反应模型的耦合很重要。液滴蒸发产生的蒸汽着火、燃烧过程的精确复现,需要物种输送的分辨率和时间步管理。
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 A | n-十二烷高温环境喷射 | 蒸汽渗透长、液相到达距离 |
| 喷雾干燥 | 实机数据 | 出口温度、残留水分 |
先在单液滴确认d-二乘法则,再在喷雾系中验证全体行为是标准方法。
衣料干燥的CFD——干衣机制造商的节能设计革命
家用衣料干衣机开发中CFD的本格活用是2010年代开始的。模拟滚筒内温湿度分布、空气流动、布料含水率变化的耦合模型,通过优化滚筒形状和通气孔位置,实现消费电力15%削减。难点是"布多孔质内毛细管蒸发"和"滚筒强制对流蒸发"两阶段的模型化。欧洲能源标签规制严格化,这种精细模拟成为家电制造商竞争力源泉。
蒸发的软件比较
商用工具比较
请比较实现蒸发模型的工具。
| 工具 | 单成分 | 多成分 | Flash蒸发 | 特点 |
|---|---|---|---|---|
| Ansys Fluent | Convection/Diffusion | Raoult, UNIFAC | Thermal breakup | VOF-to-DPM耦合 |
| STAR-CCM+ | Abramzon-Sirignano | 多成分 | Flash boiling | 拉格朗日/欧拉交换 |
| OpenFOAM | sprayFoam内置 | 支持 | 有限 | 可自由添加模型 |
| CONVERGE | 详细蒸发 | 多成分 | 支持 | AMR + 燃烧耦合 |
用途别推荐
根据用途选择吗?
Flash boiling是什么?
高温高压燃料从喷口喷射时,急速减压使液温超过饱和温度,液滴内部急速沸腾。GDI引擎高温条件下发生,促进喷雾细化。与常规表面蒸发机理不同,需要专用Flash boiling模型。
Fluent 2021以后装备Thermal Breakup模型,STAR-CCM+装备Flash Boiling Spray Model。
工艺模拟器与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液滴蒸发的直接计算
机器学习的活用
用AI改进蒸发模型的研究也有吗?
以DNS数据为教师数据,用神经网络构建液滴蒸发速度的代理模型研究有。多成分蒸发中组成-温度-压力高维参数空间的表查询,用学习模型效率化的尝试也有。
超临界蒸发——液体火箭引擎的"另一个世界"
液氢、液氧为推进剂的液体火箭引擎中,燃烧室压力超过推进剂临界压力的超临界状态下发生喷射、混合、燃烧。此条件液滴蒸发概念消失,液体到气体的"连续过渡"支配,古典蒸发模型根本失效。采用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启用 |
| OpenFOAM | sprayFoam的evaporationModel中AntoineCoeffs单位系(Pa或bar)注意 |
| CONVERGE | 确认自适应网格细化等级在蒸发区域充分 |
液滴不蒸发——边界条件设置错误的诊断
蒸发CFD中"液滴温度达到沸点但质量不减"的现象,多数情况是液滴周围气体边界条件问题。液滴蒸发速度与周围气体水蒸气分压(湿度)和平衡蒸气压的差成正比,入口边界湿度100%时蒸发驱动力为零。Ranz-Marshall等热质传递相关式在液滴Reynolds数0.1~1000范围有效,超微小液滴(< 5 µm)或高速液滴(> 200 m/s)需要补正系数。发现蒸发不发生时,推荐首先确认各单元蒸发源项的等值线图。
价值
详细
错误