实在气体效果 — 故障排除指南
更详细的内容请查看 real-gas-effects.html。
常见故障和解决方法
实在气体计算中频繁出现哪些故障？
让我们看一下实在气体特有的问题模式。
1. 非物理的负压力·负温度
症状：PR-EOS计算中温度或压力变为负数并发散
这种情况经常发生吗？
非常常见。PR-EOS在低温·低密度区域可能返回非物理压力。特别是在范德瓦尔斯环内（气液共存区）压力可能变为负数。
对策：
- 切换到表格法（表格不覆盖非物理区域）
- 设置温度和压力的下限裁剪（Fluent的
temperature-limits） - 谨慎设置初始条件（避免低温区域）
- 先用理想气体得到初始解，再切换到实在气体
2. 拟临界转换处发散
症状：超临界流体跨越拟临界温度区域时残差增大
物性在这个区域急剧变化。
是的。$c_p$ 跳跃十倍，密度急剧变化，导致数值格式不稳定。
对策：
- 提高表格分辨率（拟临界线附近以0.1 K间隔）
- 减小伪时间步长（下弛缓系数设置为0.3以下）
- 使用基于压力的耦合求解器（比SIMPLE类更稳定）
- 添加人工扩散（数值粘性）以稳定化（虽会降低精度）
3. RGP表格范围超出错误
症状：Temperature/Pressure out of table range 错误导致停止
计算过程中条件超出了表格范围。
初期过渡或局部膨胀/压缩时温度·压力可能超出预期范围。
对策：
- 将表格范围设置比计算条件宽20-30%
- 启用表格外推设置（Fluent的RGP选项）
- 将初始条件设置在计算条件的中心值，以减小过渡振幅
4. 输运系数不匹配
症状：物性表中密度正确，但粘性或热传导率不连续
虽然EOS的密度正确，但输运系数有问题的情况。
在Fluent中使用PR-EOS时，有人直接对粘性使用Sutherland公式，但在超临界域Sutherland公式完全不准确。
对策：
- 将粘性·热传导率也包含在RGP表格中（REFPROP输出为4列表格）
- 使用Chung等人的相关式（高压域用）
- 对于超临界CO₂，推荐使用Fenghour相关式
5. 混合物的EOS参数
多成分混合物中二元相互作用参数 $k_{ij}$ 如何确定？
$k_{ij}$ 是进入PR-EOS混合法则的参数，通过拟合气液平衡数据确定。基本上使用文献值，但对于无数据体系只能设置 $k_{ij} = 0$。
应确认的内容：
- 在DDB（多特蒙德数据库）或NIST TDE中调查 $k_{ij}$ 的文献值
- 进行 $k_{ij}$ 敏感性分析（观察 $\pm$0.05变化如何影响结果）
- 如可能，与相平衡实验数据比对以确认EOS预测精度
实在气体故障排除需要热力学和数值方法两方面的知识。
是的。不仅需要理解CFD设置，还要理解状态方程的适用范围和物性数据的质量，否则容易出现看似合理但错误的结果。
实在气体计算中"音速变为负数"的故障为什么会发生
在使用实在气体模型的计算中，偶尔会看到"音速变为虚数（负数的平方根）导致计算失败"的错误。理想气体中音速 $a = \sqrt{\gamma RT}$ 总是正数，但实在气体从状态方程的偏微分求音速，因此数值上 $(\partial P/\partial \rho)_s$ 可能变为负数。很多情况下是表格补间范围超出，或初始条件进入状态方程"物理上无意义的区域"所致。对策是扩大状态表的范围，或确认初始条件在相图中避开气液二相区。
相关主题
价值
详细
错误