非共沸混合物相变方向的模拟求助（可有偿）

wwyycc

我用人家论文中改进的LEE模型方法写了UDF模拟相变过程时，得到的气液相的组分含量还有换热系数和实验结果差很多。想找人交流一下模型改进和UDF编写这方面的问题。
对于非共沸混合物饱和温度不断变化，气液相间组分含量也是不断变化的这一情况，大佬是如何改进相变模型的呢？如何确定饱和温度和组分的传质速率的呢？

lwh8826

非共沸混合物相变模拟：模型改进与 UDF 完整解决方案
问题：论文中改进的 LEE 模型写 UDF 模拟相变，气液相组分含量和换热系数与实验结果差很多。
根因：标准 LEE 假设单一固定饱和温度 + 纯组分，对非共沸混合物系统性失效。
---
一、为什么 LEE 模型差很多（根因）
标准 LEE 模型（Fluent 默认）：
蒸发: m_vap = r·α_l·ρ_l·(T − T_sat)/T_sat     （T > T_sat）
冷凝: m_cond = r·α_v·ρ_v·(T_sat &#8722; T)/T_sat   （T < T_sat）
对非共沸混合物（氨-水、R32/R134a、烃类混合物等）：
1. 饱和温度不是定值：非共沸有温度滑移，泡点 T_bub ≠ 露点 T_dew，T_sat 随组分/压力连续变化
2. 组分在相变中分离：轻组分优先蒸发、重组分优先冷凝，气液相组分含量动态变化
3. LEE 用固定 T_sat + 不传组分 → 相变驱动力错、组分不输运 → 组分含量和换热系数全偏
结论：必须抛弃"固定 T_sat 的纯 LEE"，改用"动态泡露点 + 组分平衡传质"。
---
二、总体框架（推荐架构）
每个网格单元（UDF 源项）：
  P, T, 组分 x_i（液相摩尔分数）/ y_i（气相摩尔分数）
      ↓
  ① 动态饱和温度 T_sat(x_i, P)   ← 闪蒸/泡露点
      ↓
  ② 相变方向判断：T > T_sat → 蒸发；T < T_sat → 冷凝
      ↓
  ③ 总相变速率 m（LEE 扩展，或能量平衡）
      ↓
  ④ 组分平衡分配 m_i（K_i = y_i/x_i）
      ↓
  ⑤ 三个源项：质量 + 能量(潜热+显热) + 组分(输运)
---
三、核心①：动态饱和温度（泡点/露点）计算
泡点温度（液相给定 x_i，求 T）：
∑ x_i·γ_i·P_sat_i(T) = P
露点温度（气相给定 y_i，求 T）：
∑ y_i / [γ_i·P_sat_i(T)] = 1/P
其中 P_sat_i(T) 用 Antoine 方程：
log10(P_sat_i) = A_i &#8722; B_i/(C_i + T)
UDF 泡点迭代（Raoult 定律版，γ=1）：
/* 输入: 液相摩尔分数 x[N], 压力 P; 输出: 泡点 T_bub */
real bubble_temperature(real x[], real P) {
    real T = P_sat_average_guess(x, P);   /* 各组分饱和温度的质量平均初值 */
    for (iter = 0; iter < 50; iter++) {
        real sum = 0.0;
        for (i = 0; i < N; i++)
            sum += x[i] * Antoine_p_sat(i, T);
        real f = sum - P;                  /* 泡点方程残差 */
        /* 数值求导 df/dT（中心差分） */
        real T2 = T + 1e-3, sum2 = 0.0;
        for (i = 0; i < N; i++)
            sum2 += x[i] * Antoine_p_sat(i, T2);
        real dfdt = (sum2 - sum) / 1e-3;
        real dT = -f / (dfdt + 1e-12);
        if (fabs(dT) < 1e-6) break;
        T += dT;
    }
    return T;
}
非共沸两相区：泡点 < 界面温度 < 露点。相变驱动用泡点（液相蒸发）和露点（气相冷凝）可分开处理，或用一个界面平衡温度。
---
四、核心②：组分传质速率（平衡分配）
气液界面假设局部平衡，用平衡常数：
K_i = y_i / x_i = γ_i·P_sat_i(T) / P
各组分传质速率（改进 LEE，组分加权）：
液相组分 i 蒸发:  m_i = r·α_l·ρ_l·(T &#8722; T_sat)/T_sat · x_i·w_i
气相组分 i 冷凝:  m_i = r·α_v·ρ_v·(T_sat &#8722; T)/T_sat · y_i·w_i
w_i 是组分权重（体现轻组分优先蒸发：轻组分 w 大、重组分 w 小），
或用平衡分配：液相组分转气相时按 y_i（平衡气相组成）分配。
严格做法（界面平衡，能量主导）：
总相变速率 m 由能量平衡决定:
  m·h_fg(T, x) = Q_传热          (潜热 = 界面能量守恒)
组分按平衡分配:
  m_i = m · x_i_interface 或 m·K_i 归一化
---
五、核心③：三个源项（UDF 必须全加）
质量源项（连续性方程）：
液相:  S_m_l = &#8722;m    （蒸发损失）
气相:  S_m_v = +m    （冷凝增加）
能量源项（潜热 + 显热）：
S_energy = &#8722;m·h_fg(T, x)
&#9888;&#65039; 非共沸混合物 h_fg 不是常数，随组分/温度变化：
h_fg(T,x) = ∑ x_i·h_fg_i(T)     （各组分潜热按组分加权）
组分输运源项（关键！不加速率组分不传质）：
液相组分 i:  S_Y_li = &#8722;m_i
气相组分 i:  S_Y_vi = +m_i
UDF 源项框架（DEFINE_SOURCE）：
/* 液相质量源项 */
DEFINE_SOURCE(liq_mass_source, c, t, dS, eqn) {
    real P = C_P(c,t), T = C_T(c,t);
    real x[N];  get_liq_molefractions(c, t, x);      /* 取液相组分 */
    real Tsat = bubble_temperature(x, P);            /* 动态泡点 */
    real alpha_l = C_VOF(c, t, 0);                   /* 液相体积分数 */
    real rho_l = C_R(c, t, 0);
    real r = relaxation_factor;                      /* 松弛因子 1e3~1e6 */
    real m;
    if (T > Tsat) m = r * alpha_l * rho_l * (T - Tsat) / Tsat;  /* 蒸发 */
    else          m = 0.0;
    dS[eqn] = ...;   /* 源项对变量的导数，保证收敛 */
    return -m;       /* 液相损失 */
}
/* 能量源项、组分源项同理，返回 ±m·h_fg 和 ±m_i */
---
六、换热系数处理
局部换热系数：
h = q'' / (T_wall &#8722; T_sat(x, P))
&#9888;&#65039; 换热系数偏差大的 3 个关键：
1. T_sat 必须动态算（用泡露点），固定 T_sat → 温差错 → h 错
2. 潜热 h_fg 随组分加权（非共沸不是常数）
3. 松弛因子 r 标定：r 影响相变速率，太大震荡、太小偏差，对网格/工况标定
---
七、调参策略
参数        建议        作用
松弛因子 r        10&#179;~10&#8310;，网格/工况标定        相变速率收敛性
平衡常数模型        Raoult → 活度系数(γ) → 状态方程        组分分配精度
组分输运        必须加源项        组分含量正确性
时间步长        满足相变源项刚性        避免换热系数虚高
潜热        组分加权        换热系数正确性
---
八、验证方法（三步）
1. 纯组分验证：先用纯工质跑 LEE，确认基础（T_sat 定值、单组分）吻合
2. 平衡验证：无相变（冻结流）验证泡露点计算和组分平衡
3. 相变验证：非共沸工况，对比实验的组分含量 + 换热系数 + 温度分布
---
一句话总结
非共沸相变 = 动态泡露点（闪蒸迭代）+ 组分平衡传质（K_i 分配）+ 三源项全加（质量/能量/组分）+ 组分加权潜热。抛弃固定 T_sat 的纯 LEE，换热系数和组分含量自然对上。
---
*如需针对具体工质（氨-水/R32-R134a 等）定制闪蒸函数或组分分配，告知工质和工况即可深化。*