J3-04 · 仿真分析 CAE / 换热器/部件级 CFD
待认领listed
→
众筹中recruiting
→
已开课open · Course 售卖
一课题可并存多讲师版本
课程大纲
学完能做什么
- 能按流型/尺度选出 VOF、Mixture 或欧拉-欧拉两相方法
- 能按所选两相框架配置壁面沸腾（欧拉两流体配 RPI 三段热流分配、VOF 配 Rohsenow 型或界面解析型）与冷凝（Lee 型）相变模型
- 能控制 Courant 数、时间步与界面锐化保证收敛
- 能判断表面张力/润湿性对结果的敏感度
- 能把两相 CFD 结果与部件试验对标并识别不可信结果
内容大纲
- 两相方法选型
- 界面追踪 VOF（分层/塞状/大界面场景）
- Mixture 混合模型（弥散相、低算力）
- 欧拉-欧拉双流体（泡状/雾状、含相间作用）
- 按流型（A2-03）与网格能力匹配方法，别硬套
- 框架决定可用的相变模型：RPI 三段热流分配只在欧拉两流体下提供（Fluent 需 Eulerian 或 Multi-Fluid VOF，STAR-CCM+ 需 EMP Wall Boiling）；经典 VOF 能配的是界面解析型相变（Lee/Tanasawa/Schrage）或 Rohsenow 型亚网格核态沸腾（发动机水套过冷沸腾的主流做法）——先定框架再选相变模型，顺序反了必然卡在前处理
- 相变模型
- 蒸发/冷凝质量源（Lee 模型）与经验系数标定
- 壁面沸腾 RPI（只在欧拉两流体框架下提供）：对流加淬灭加蒸发三段热流分配
- 冷凝：膜状 vs 珠状、过冷液膜的处理
- 相变潜热与能量源的闭合
- 数值稳定性与时间步控制
- Courant 数与显式/隐式 VOF 的时间步限制
- 界面锐化/压缩与伪流（spurious currents）
- 表面张力 CSF 模型与网格分辨率要求
- 亚松弛、瞬态推进与收敛判据设置
- 整车两相场景
- 制冷剂直冷板/蒸发器内的沸腾与干涸（dryout）
- 冷凝器过热段-两相-过冷段的相变分布
- 两相浸没冷却的沸腾换热与气泡管理
- CO₂ 跨临界近临界物性剧变的挑战
- 对标、动态失稳与结果可信度
- 与可视化/焓差部件试验对标
- 相变经验系数按工况标定，而非用默认值
- 网格/时间步敏感性与不可信结果的识别
- 与 1D 系统（J1-02 制冷循环一维）的边界与工况交接；与相邻方法课的边界：本课是相变两相（难在相变源项与界面锐化），无相变的自由液面 + 旋转域（难在交界面界面穿越与液膜破碎）见 J3-06，两课互为前置/关联
- 稳态两相 CFD 在结构上就预测不了动态失稳：要复现并联通道的 Ledinegg 流量漂移或密度波振荡，必须做瞬态，并把上游可压缩容积（储液腔、闪蒸段）与真实边界建进模型，且边界条件类型（定压差 vs 定流量）本身就决定算不算得出失稳（机理见 A2-03）
- 试验上看到支路流量摆动、而仿真是稳态收敛的，不构成『仿真对上了』的证据——稳态解天然不含失稳模式，应先换瞬态 + 真实边界与上游容积重算，再谈误差对标
关键公式
∂α/∂t + ∇·(αu) = ṁ/ρ
VOF 相体积分数输运方程，右端为相变质量源，α 从 0 到 1 描述液气界面
Lee 模型 ṁ = C·α·ρ·(T−T_sat)/T_sat
蒸发(T>T_sat)/冷凝(T<T_sat)质量传输；系数 C 必须按工况标定，非通用常数
q_wall = q_conv + q_quench + q_evap
RPI 壁面沸腾把壁面热流分成对流、淬灭、蒸发三段，是沸腾换热建模主线
F_σ = σ·κ·∇α（CSF）
连续表面张力体积力，κ 为界面曲率；网格不足会诱发伪流淹没真实流动
Co = u·Δt/Δx
Courant 数，界面追踪需 Co<1（常需更严），是控时间步的第一红线
能量源 = ṁ·h_fg
相变潜热闭合：蒸发吸热、冷凝放热，与质量源同步施加
关键概念
VOFMixture 模型欧拉-欧拉双流体Lee 相变模型RPI 壁面沸腾表面张力 CSFCourant 数界面锐化伪流 spurious currents干涸 dryout膜状/珠状冷凝相变潜热
推荐工具与标准
Fluent（VOF/Eulerian/Wall Boiling） STAR-CCM+（两相与沸腾） OpenFOAM（interFoam 等两相求解器） CoolProp/REFPROP（工质物性） 本站冷媒工况计算器
ASME V&V 20（CFD 验证与确认） AIAA CFD 验证与确认指南（G-077） 说明：两相/沸腾 CFD 无专用标准，经验系数须以部件试验标定
工程案例
某纯电平台电池制冷剂直冷板，末端电芯偏热，怀疑制冷剂提前干涸。用 VOF 加界面解析型相变（Lee/Tanasawa）或 Rohsenow 型壁面沸腾建直冷板内两相流（若改走欧拉两流体 + RPI，须再叠 CHF/流型转换扩展才算得到环状流液膜耗尽型干涸），复现进口过冷液—中段沸腾—末端干涸的相分布，量化干涸位置随冷媒流量与热流的移动，据此调分液与流道把干涸点推到芯体外，末端温差收敛。（直冷走冷媒无水路，勿与液冷混模）
动手做
交付物 · 取一段制冷剂蒸发流道：① 选两相方法并说明依据（流型/尺度）；② 配置 Lee 或壁面沸腾模型，用一组试验数据标定相变系数；③ 做时间步/网格敏感性，给出干涸位置与换热分布，交付一份含对标误差与“结果可信度”结论的报告。
常见误区
- 用 Lee 模型默认系数不标定，相变速率随手一填，结果全靠系数“凑”
- 时间步过大 Co≫1，界面弥散成一团、假收敛还自我感觉良好
- 网格分辨率不足，表面张力伪流淹没真实流动
- 大界面分层流硬套弥散 Mixture 模型（或反之），流型与方法不匹配
- 只看漂亮的两相云图，不与试验对标就下结论
- 在经典 VOF 下到处找 RPI 壁面沸腾选项，找不到就以为软件缺功能——错的是框架与相变模型没配对，该换的是框架不是模型
相关课题