格拉晓夫数(Grashof Number)——自然对流的驱动力标尺

格拉晓夫数(Grashof number,简称 Gr)表征浮升力与粘性力的相对强弱,是自然对流分析中最重要的无量纲参数。没有外加驱动的流动问题,比如暖气片周围的空气爬升、芯片在无风扇条件下的散热,流动由密度差驱动,此时判断”流动会不会发生、处于什么流态”,首先就要看 Gr 数。它之于自然对流,大致相当于雷诺数之于强制对流。

定义与表达式

$$
Gr = \frac{\text{浮升力}}{\text{粘性力}} = \frac{g \beta \Delta T L^3}{\nu^2}
$$

其中:

$g$:重力加速度(m/s²);
$\beta$:流体热膨胀系数(K⁻¹),反映温度变化引起的密度变化;
$\Delta T$:特征温差(K),如壁面与主流温差;
$L$:特征长度(m),如平板高度或管道直径;
$\nu$:运动粘度(m²/s)。

注:

对于理想气体,$\beta = 1/T$($T$ 为绝对温度);

若流动由浓度差驱动(如溶质扩散),$\Delta T$ 需替换为浓度差 $\Delta C$,此时称溶质 Grashof 数($Gr_c = \frac{g \beta_c \Delta C L^3}{\nu^2}$)。

物理意义:自然对流的控制参数

Gr 数决定自然对流是否发生及其流态:

Gr < 10⁸:粘性力主导,流动为层流,热传递以导热为主;例如微电子散热中的低速气流。
10⁸ ≤ Gr ≤ 10¹⁰:浮升力与粘性力相当,可能发生层流向湍流过渡;典型场景如暖气片附近的空气流动。
Gr > 10¹⁰:浮升力主导,流动发展为湍流,热传递效率显著提升;常见于大型换热设备或大气环流。

与雷诺数(Re)的关系:在混合对流中,Gr/Re² 表征自然对流与强制对流的相对重要性:

Gr/Re² ≫ 1 → 自然对流主导(如无风时的热源散热);
Gr/Re² ≪ 1 → 强制对流主导(如风扇冷却);
Gr/Re² ≈ 1 → 两者均需考虑。

在 CFD 中的应用

在计算流体动力学(CFD)中,Gr 数直接影响数值模型的选择与边界条件设置:

湍流模型修正:高 Gr 数需启用浮升力修正(如 OpenFOAM 中的 buoyantKEpsilon 模型);忽略修正可能导致热边界层预测偏差。
网格设计要求:近壁区域需加密网格以解析薄热边界层(尤其高 Pr 流体);网格尺寸应满足 $y^+ \approx 1$(粘性子层分辨率)。
边界条件设置:壁面采用 buoyantPressure 或 fixedFluxPressure 处理浮升力引起的压力梯度;温度边界需匹配密度变化(如 temperatureDependent 粘度)。

工程应用场景

领域	典型问题	Gr 数的作用
建筑环境	室内通风设计	优化窗户位置促进自然对流(Gr > 10⁹ 时效率最佳)
电子散热	芯片无风扇冷却	评估散热片间距(Gr 决定最优肋片高度)
能源工程	太阳能集热器效率	控制流道尺寸避免层流抑制传热(Gr/Re² > 0.1)
化工设备	反应釜温度均匀性	预测自然对流混合强度(Gr ∝ 流动动能)

典型案例:

竖直平板散热:Nu(努塞尔数)经验公式为 $Nu = C (Gr \cdot Pr)^n$($C, n$ 为常数),直接关联换热效率;

核反应堆安全:事故工况下冷却剂自然循环能力由 Gr 数评估。

与其他无量纲数的关联

Gr 数常与以下参数联用:

Rayleigh 数(Ra):$Ra = Gr \cdot Pr$,综合表征浮升力与热扩散的竞争,更直接决定自然对流稳定性;
普朗特数(Pr):影响速度与热边界层相对厚度(高 Pr 流体热边界层更薄);
Archimedes 数(Ar):$Ar = Gr / Re^2$,专用于混合对流分析。

小结

Grashof 数是解析自然对流问题的核心标尺:

机理层面:量化浮升力/粘性力平衡,预测流态与换热模式;
设计层面:指导换热设备几何优化(如散热片间距、反应器高径比);
数值层面:约束 CFD 模型选择与网格划分策略。

掌握 Gr 数,就掌握了电子散热、建筑通风、能源设备等自然对流工程问题的分析入口。

本文整理自个人 CFD 学习笔记系列