瑞利数 — CAE术语解说

很好的问题。它表示浮力与粘性力、热扩散之间的竞争关系。具体来说，是浮力产生的动能生成率与粘性耗散率和热扩散引起的温度梯度缓和率之比。用公式表达是：

观察很敏锐。格拉斯霍夫数

考虑一个从下方加热的水平液体层。当瑞利数较小（约1708以下）时，热只能通过传导传递，流体保持静止，处于稳定状态。但当超过临界瑞利数（约1708）后，浮力战胜粘性和热扩散的阻力，自发地产生有规律的六边形对流单元（Bénard对流）。这就是"对流不稳定性"的开始。在电子产品散热中，基板和散热器之间的空气层可能出现这种现象。

完全正确。通常我们分别输入物性参数（密度、粘度、热膨胀系数、热导率、比热）和边界条件（温度差、代表长度）。求解器在内部计算瑞利数，从而判断流动的物理状态。比如，在Ansys Fluent的"Boussinesq近似"模型中，直接输入热膨胀系数β，结合操作温度和重力加速度的设置，浮力项就能被正确计算，瑞利数也就随之确定了。

对数值解的精度影响很大。高瑞利数（例如10^6以上）的流动非常复杂，会形成很薄的速度和温度边界层。要解析这些边界层，需要在墙面附近使用极细的网格（y+<1）。此外，求解器的压力-速度耦合算法（如Fluent的Coupled算法）和对流项的离散化格式（高阶格式）的选择影响收敛速度。设置不当会导致非物理振荡或发散。

取决于温度差ΔT的大小。Boussinesq近似适用于密度变化较小的情况（比如空气中ΔT为几十K），计算较为稳定。而对于大的温度差（几百K）或高度方向上下降明显的模型，密度变化不容忽视，需要用非压缩理想气体法则或完整的状态方程。当瑞利数公式中的βΔT项变得过大时，Boussinesq近似的假设（仅在浮力项中考虑密度变化）就会失效。

可以用于设计初期的筛选和CFD结果的验证。首先，根据预期的温度差和代表长度（如散热器高度）手工计算大致的瑞利数。若Ra超过10^8，流动进入湍流区，CFD需要选择湍流模型（k-ε或SST k-ω）。分析完成后，从仿真结果计算努塞尔数（无量纲热传递系数），与实验相关式比对（如竖直平板层流：

要根据具体情况。如果评估翅片间的自然对流，就取翅片间距作为L。如果想用竖直平板相关式评估整个散热器的散热，就用翅片高度。关键是与对比对象（文献实验式或以前的仿真案例）保持一致的定义。有些环境试验规范如JIS C 60068或MIL-STD-810会规定特定条件下的自然对流热传递系数，可以从这些条件反推代表长度。

温度差较大时必须考虑。空气的动力粘度随温度上升而增加，热膨胀系数与绝对温度成反比。若忽略这些，用常温物性，瑞利数可能被严重低估，可达数倍之差。在CFD设置中，应将物性定义为温度的多项式函数或理想气体法则。Ansys的材料库中有基于NIST数据的温度相关物性，可以利用。

物理基础相同，但设置哲学和初值生成有差异。Fluent采用传统的"操作温度"概念，用户设定它，然后用相对于它的温度差计算浮力。Star-CCM+的"密度变化模型"选择更明确（如定密度、Boussinesq、理想气体等），参考密度和参考温度可以分别设置。COMSOL选择"非等温流"物理接口后，重力项和温度场自动耦合，所有物性从材料库继承。COMSOL在处理弱非线性初值化时有优势。

OpenFOAM选定求解器（如`buoyantBoussinesqSimpleFoam`）时就启用了Boussinesq近似。物性在`transportProperties`文件中定义。最大区别是商用软件隐藏的"参考压力点"设置，用户常需显式指定。参考点设置错误时，整体压力场偏移，影响浮力驱动流的初期加速，导致收敛缓慢。这与瑞利数的大小无关，是重要的配置。

这些软件的特色是用"经验相关式"处理自然对流热传递。FloTHERM可在器件表面设置"自然对流"热阻，基于瑞利数相关式（如

首要原因是数值瑞利数接近物理临界值（1708），或远远超过它，而网格和求解器跟不上。当Ra接近临界值时，解对扰动敏感，微小的非对称性或舍入误差都会被放大，造成振荡或收敛缓慢。一个对策是有意在初值中引入实际存在的微小扰动（如0.1K的温度不均一），以激发对流单元。

分阶段方法有效。先将重力设为零，用纯热传导条件计算到收敛，再以此结果为初值逐步增加重力（比如从1/10倍增加到实值），这称为"参数扫描"。Fluent的"混合初值化"功能可自动完成这个过程。也可以降低求解器的松弛系数（从默认0.3降到0.1），或先用一阶迎风格式离散对流项，收敛后再改为高阶格式。

可以作为粗略指标。Ra超过10^9的强湍流自然对流本质上是非定常的，会产生涡旋。定常分析残差持续振荡时，通常表示物理上存在非定常现象。若强行用定常分析收敛，会得到非物理的平均化结果。时间步应小于流体流过代表长度L的时间（L/U，U为特征速度）。瑞利数越高，流速越大，需要的时间步越小。

第一步，用自己的CFD条件重新计算瑞利数，确认与对比文献的Ra一致。特别注意代表长度L和温度差ΔT的定义是否相同。第二步，检查物性参数。文献实验式通常在特定温度（如膜温）下评估物性，确保你的CFD也是在那个温度评估的。第三步，核对边界条件：墙面是"等温"还是"恒热流"，相关式本身可能不同。