临界热流密度 — CAE术语解说

正确。这是加热面向液体的热传递达到最大的临界点，超过这个点后，加热面被蒸汽膜覆盖，发生"膜沸腾"，热传递系数急剧下降。例如，大气压下水在水平向上加热的情况下，CHF约为1.1 MW/m²。这种转变通常用Zuber理论模型很好地解释。

这是针对无限水平平板的经典模型。基于气泡的Helmholtz不稳定性，CHF

随着压力升高，蒸发潜热下降，气相密度增加。在Zuber公式中，

直接预测非常困难。CHF强烈依赖于微观尺度的气泡动力学和加热面的湿润性，用DNS（直接数值仿真）重现需要微米级以下的网格，实际计算成本几乎不可行。因此，在实务中采用相关式和亚网格模型的方法。

例如，欧洲原子核学会（ERECOM）开发的"NURESIM"项目中的模型。这基于CFD单元内的局部蒸汽体积分数和热流密度，判断是否达到CHF。判定式中包含查表法（LUT）或

典型工作流程是这样的：首先，用RANS（雷诺平均纳维-斯托克斯）模型和欧洲多相流流动样式图（EMM）预测流动样式。接着，在每个计算单元中比较局部热流密度和从相关式得出的局部CHF值。当局部热流密度超过局部CHF值的单元达到一定比例时，判定整个系统达到CHF。可以用Ansys Fluent的用户定义函数（UDF）或STAR-CCM+的字段函数实现这个逻辑。

通常采用"最小临界热流密度比（MDNBR）"作为设计基准。这是将预期的最大热流密度除以从实验相关式求得的临界热流密度，取其最小值。该值不应低于1.0。轻水反应堆通常设定MDNBR > 1.3的余度。

根据流路形状和适用范围选择。圆管内强制对流沸腾常用"W-3相关式"（西屋电气公司）。沸水型轻水反应堆（BWR）燃料棒束使用"通用电气（GE）相关式"或"美国电力研究所（EPRI）相关式"，加压水型轻水反应堆（PWR）使用"WRB-1相关式"等。近年来，国际上更多参考基于更广泛数据库的"Groeneveld查表"。

至少要检查以下3点：1) 使用的相关式的适用范围（压力、质量流量、干度、流路形状）是否包含设计条件。2) 热流密度和流量分布在物理上是否合理（考虑入口效应、旁通流等）。3) 进行网格敏感性分析，确认CHF预测值不依赖于网格尺寸。特别是壁面附近的网格分辨率很重要。

没有专门的"CHF求解器"标准功能。但两者都具有高级多相流沸腾仿真功能。Ansys Fluent可以将"欧拉多相模型"与"RPI（伦斯勒理工学院）沸腾模型"结合，计算详细的气泡行为。STAR-CCM+除了"VOF（流体体积）法"和"欧拉多相"外，还提供专门的"沸腾"物理模型，可以追踪接近膜沸腾转变的状态。

系统代码和亚通道分析代码中，CHF评估是核心功能。例如，美国广泛使用的"TRACE"（美国核管理委员会）、"RELAP5"，法国的"CATHARE"，日本的"SPIRAL"等。这些是一维流管模型，但内置了大量基于实验的高精度CHF相关式（Katto、Biasi、Groeneveld LUT等），对整个堆芯的安全分析至关重要。

确实有限制。COMSOL的"二相流、移流界面"模块可以追踪气液界面，但沸腾引起的激烈相变以及气泡合并和分裂导致的CHF现象，用现有计算资源直接仿真非常困难。使用COMSOL的实际做法是，计算流场和温度场后，用户另外基于相关式计算局部CHF值，在后处理中比较。

主要有两个原因。1) 源项过陡峭：蒸发、凝结源项太大，数值上不稳定。对策是线性化源项或使用隐式求解器。2) 网格和时间步长不当：要解析气液界面，需要足够细的网格（通常界面宽度的3～5个单元以上）和满足CFL数≤1的足够小的时间步长。

首先检查相关式的输入变量单位是否正确。特别是压力（Pa vs. MPa）和热流密度（W/m² vs. kW/m²）容易出错。其次，检查是否正确计算了相关式所需的"局部条件"。例如，相关式要求"平衡干度"，但输入的是CFD结果的"非平衡（实际）蒸汽率"。热流密度分布如果有极端峰值，可能会大幅低于局部CHF。

Groeneveld查表（2006年版）以压力、质量流速、干度和内径为独立变量。重要的是不能在表的范围外进行外推。例如，表的压力范围是0.1～21 MPa。另外，查表是针对圆管的，应用于矩形流路或棒束时，需要用等效直径进行适当的修正（束修正系数）。忽略这一点会导致过高或过低估计。