普朗特数(Prandtl Number)——动量与热扩散的竞速

普朗特数(Prandtl number,简称 Pr)表征动量扩散与热扩散的相对速率,直接影响流动与传热的耦合行为。它有一个独特之处:Pr 纯由流体物性组成,与流动尺度、速度无关,同一状态下空气永远约 0.7、机油上千、液态金属约 0.01。这个”物性指纹”决定了速度边界层与热边界层谁厚谁薄,进而在边界层相似性、换热经验关联式和 CFD 数值格式选择中反复出现。

定义与表达式

$$
Pr = \frac{\text{动量扩散率}}{\text{热扩散率}} = \frac{\nu}{\kappa} = \frac{\mu c_p}{k}
$$

其中:

$\nu$:运动粘度(m²/s),反映动量传递能力;
$\kappa$:热扩散率(m²/s),$\kappa = \frac{k}{\rho c_p}$,反映热量传递能力;
$\mu$:动态粘度(Pa·s);
$c_p$:定压比热容(J/(kg·K));
$k$:热导率(W/(m·K))。

物理上,Pr 描述流体中速度边界层与热边界层的相对厚度。例如:

Pr ≫ 1(如机油 Pr≈1000):速度边界层发展快于热边界层(热量传递慢);
Pr ≪ 1(如液态金属 Pr≈0.01):热边界层远厚于速度边界层(热量传递快);
Pr ≈ 1(如空气 Pr≈0.7):两者发展速率相近。

量纲一致性:Pr 是无量纲数,公式中各物理量单位可相互约去,确保普适性。

物理意义:边界层行为的决定因子

流动与传热的耦合机制
- Pr 数直接关联 Navier-Stokes 方程(动量守恒)与能量方程(热量守恒);
- 在边界层内,速度分布 $u(y)$ 与温度分布 $T(y)$ 的相似性由 Pr 决定:

Pr 范围	边界层特征	典型流体
Pr ≪ 1	热边界层较厚,温度分布平坦	液态金属(汞、钠)
Pr ≈ 1	速度与热边界层厚度相当	空气(Pr=0.7)
Pr ≫ 1	热边界层极薄,温度梯度大	机油、甘油

对换热效率的影响
- 高 Pr 流体(如油类)传热效率低,需强化换热(如增大流速或扰动);
- 低 Pr 流体(如液态金属)传热效率高,常用于核反应堆散热。

在 CFD 中的关键作用

数值格式稳定性
- 网格 Peclet 数(Peᵢ)需满足 Peᵢ = Reᵢ · Pr ≤ 2,否则中心差分格式失稳,需改用迎风格式;
- OpenFOAM 中 limitedLinear 等格式会根据局部 Peᵢ 自动切换限制器。
壁面函数选择
- 高 Pr 流体需精细的近壁网格(y⁺≈1),否则热边界层分辨率不足;
- 低 Pr 流体可使用标准壁面函数(如 nutkWallFunction)。
多场耦合问题
- 空化模拟(如 interPhaseChangeFoam):Pr 影响气-液界面的相变速率;
- 热毛细对流(Marangoni 效应):$Ma = Re \cdot Pr$,控制表面张力驱动的流动。

工程应用场景

领域	应用案例	Pr 的作用
能源工程	换热器设计	优化管径与流速,匹配 Pr 以提升传热效率
化工反应器	流体混合与反应速率控制	低 Pr 流体(如熔盐)加速反应物扩散
电子散热	液态金属冷却芯片	利用低 Pr 特性实现高效热管理
环境科学	大气或水体污染物扩散模拟	Pr 修正湍流模型中的热-动量扩散比

典型案例:

液体桥热毛细对流(Marangoni 流动):$Ma = |\sigma_t| h \Delta T / (\rho \nu \kappa) = Re \cdot Pr$,其中 $\sigma_t$ 为表面张力温度系数。高 Ma 数导致非稳态涡结构;

微流控芯片:通过控制 Pr 实现快速温度均一化(需 Pr<1)。

与其他无量纲数的关联

Pe = Re · Pr:Peclet 数,表征对流与热扩散的相对强弱;
Ma = Re · Pr:Marangoni 数,表征表面张力驱动的流动强度;
Pr 与 Sc(施密特数):$Sc = \nu / D$($D$ 为质量扩散率),二者分别关联热扩散与质量扩散。

小结

普朗特数 Pr 是解析流动-传热耦合问题的核心参数:

机理层面:揭示速度场与温度场的协同演化规律;
数值层面:约束离散格式稳定性与网格设计;
工程层面:指导换热设备优化及新材料(如液态金属)应用。

在空化、燃烧、电子散热这类流动与传热深度耦合的问题里,Pr 都是绕不开的分析起点。

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