水平平板的自然对流
水平平板自然对流的理论基础
概要 — 上面与下面的决定性差异
水平板的自然对流，上向和下向完全不同吗？
完全不同。加热面向上时，板表面被加热的空气因浮力上升产生羽流（热上升流）。温热的空气不断离开，新鲜冷空气从侧面持续流入，所以热传递非常好。
相反加热面向下时，温热空气被板困住。热空气想上升但被板挡住，结果形成稳定分层结构，几乎不发生对流。热传递系数降至上面的一半以下。
具体在哪些场景会成为问题？
最常见的是PCB水平配置时的芯片冷却。芯片安装在上面时能通过自然对流有效冷却，但下面安装时极端冷却困难。服务器机架内的主板或嵌入式设备的基板天顶安装，这种差异会致命。
还有太阳热集热器（平板型）的上面热损失评估、厨房热板散热、工业厂房天花板保温设计等，水平板自然对流相关式都是必需的知识。
支配方程式与Rayleigh数
公式上怎么表述？首先是Rayleigh数吧？
对。自然对流强度的无量纲数是Rayleigh数 Ra。它表示浮力驱动与粘性·热扩散阻力的比值：
各符号的含义：
- $g$：重力加速度 [m/s²]
- $\beta$：体膨胀系数 [1/K]（理想气体则 $\beta = 1/T_f$，$T_f$ 为膜温度）
- $T_s$：板表面温度 [K]
- $T_\infty$：周围流体温度 [K]
- $L_c$：特性长度 [m]（后述）
- $\nu$：动粘性系数 [m²/s]
- $\alpha$：热扩散率 [m²/s]（$\alpha = k / \rho c_p$）
Ra数越大对流越剧烈？
完全正确。Ra数大=温度差大、板大、粘性小——浮力主导对流活跃。反之Ra数小则热传导主导，对流几乎不发生。水平板上面约 $Ra \approx 10^7$ 时从层流过渡到湍流。
Nu数相关式（McAdams·Lloyd-Moran）
这次讲相关式！上下面式子不同吧？
是的。McAdams（1954年）通过实验整理的著名相关式。首先是加热面向上（Hot Surface Up）：
以及加热面向下（Hot Surface Down）：
下面的系数0.27正好是上面0.54的一半！而且下面没有湍流式？
你观察得很仔细。下面因稳定分层对流很弱，几乎不发生湍流转变。因此1/4次幂规律能覆盖很宽的范围。系数恰好是一半意味着实际上热传递系数 $h$ 是上面的一半。相同温度差下冷却能力只有一半——设计上这是巨大的差异。
顺便说，从Nusselt数到热传递系数 $h$ 的转换是：
这里 $k$ 是流体（空气等）的热导率。
冷却面的情况是反过来吗？冷的板下面冷却效果好？
聪慧的观察。冷却面向下（Cold Surface Down）的机制与加热面向上相同——板下冷却的空气变重下沉产生羽流。所以用相同的相关式（系数0.54/0.15）。反之冷却面向上与加热面向下同理稳定分层，用系数0.27。
| 配置 | 物理现象 | 应用相关式 |
|---|---|---|
| 加热面向上（Hot Up） | 羽流上升 → 高Nu | 0.54 / 0.15 |
| 加热面向下（Hot Down） | 稳定分层 → 低Nu | 0.27 |
| 冷却面向下（Cold Down） | 下沉羽流 → 高Nu | 0.54 / 0.15 |
| 冷却面向上（Cold Up） | 稳定分层 → 低Nu | 0.27 |
特性长度 Lc = A/P 的含义
竖板用板高作特性长度。水平板用什么？
水平板用板面积 $A$ 除以周长 $P$ 的值作特性长度：
看几个典型形状的值：
| 板的形状 | 面积 $A$ | 周长 $P$ | $L_c = A/P$ |
|---|---|---|---|
| 正方形（边 $a$） | $a^2$ | $4a$ | $a/4$ |
| 矩形（$a \times b$） | $ab$ | $2(a+b)$ | $\frac{ab}{2(a+b)}$ |
| 圆板（直径 $D$） | $\pi D^2/4$ | $\pi D$ | $D/4$ |
为什么是面积÷周长？直观上不太能理解…
很好的问题。水平板自然对流中，流体通过板的边缘（边界）进出。板中心到边缘的"流体移动的代表距离"就是它要表示的。面积相对于周长大的板（如大正方形）——中心部分的流体替换缓慢。面积对周长比小的板（细长板）——各处都靠近边缘替换快。$L_c = A/P$ 标度化了这个"从边到中心的平均距离"。
上面·下面的流动结构差异
上面的羽流，具体什么形状产生？
加热面向上（层流区）时，板中央附近产生"蘑菇型"羽流准周期上升。冷空气从边缘流入，沿着板表向中心爬升并被加热，在中央上升形成"大规模环流"。
湍流区（$Ra > 10^7$）时，板全面随机产生羽流。出现类似Bénard对流的细胞结构，热传递大幅提升。1/3次幂律对应这种状态。
加热面向下则完全不同。热空气在板下形成薄的热边界层，但浮力指向板上方使流体无法脱离。只能通过极弱的横向扩散缓缓流出边缘。因此热传递极低。但当 $Ra > 10^5$ 时，"Taylor不稳定"类机制使板下也开始出现弱对流细胞。
原来板的朝向一个决定流动物理的全部…
McAdams 1954年的实验
William H. McAdams是MIT教授，在1954年的名著《Heat Transmission》第3版中系统整理了水平板自然对流相关式。他的实验很简单：在空气中水平放置加热铜板（10～30 cm见方），通过电加热器均匀输热。但"上下面分开整理"这个创意极具前瞻性，70多年后的今天仍被工程设计广泛采用。后来Lloyd & Moran（1974年）引入 $L_c = A/P$ 定义，使其适用于任意形状的板。
水平平板自然对流的数值计算方法
CFD分析的基本策略
用CFD求解水平板自然对流时什么思路？
基本是Navier-Stokes方程+能量方程的耦合求解。自然对流的速度场和温度场通过浮力项强烈耦合。要解的支配方程：
连续性方程：
动量方程（Boussinesq近似）：
能量方程：
关键是最后的浮力项 $-\rho_0 \beta (T - T_0) \mathbf{g}$。温度差驱动流动，所以必须同时求解温度和速度。
网格策略与边界层分辨
网格要多细？壁面附近很重要吧…
自然对流需要分辨速度边界层和温度边界层。水平板的温度边界层厚度约：
例如 $L_c = 0.05$ m、$Ra = 10^6$ 时，$\delta_T \approx 0.05 / 31.6 \approx 1.6$ mm。其中要配置最少10～15层网格。壁面第一层厚度约 0.05～0.1 mm。
超关键的一点是计算域的大小。域过小会阻碍诱引流，导致误差。目安是：
- 板周围水平方向：板尺寸的3～5倍以上
- 上面（羽流上游）：板尺寸的5～10倍
- 下面：板尺寸的3倍以上
域太小时Nu数被低估。实务中常用"把域扩大2倍后Nu变化不超1%"来确认域的独立性。
Boussinesq近似与完全可压缩性
Boussinesq近似什么时候用？
Boussinesq近似是"密度变化只影响浮力项，其他地方密度定常"的假设。适用条件：
空气（$T_\infty = 300$ K）时 $\Delta T < 30$ K 左右是安全。电子设备冷却（$\Delta T \sim 20$～$60$ K）接近极限。$\Delta T > 100$ K 的工业应用（炉、高功率电子）需用完全可压缩性（理想气体状态方程）。
湍流模型的选择
湍流模型选哪个？听说自然对流用k-ε会有问题？
自然对流湍流模型很讲究。标准 $k$-$\varepsilon$ 壁面行为不精确，自然对流Nu会高估30%以上。推荐顺序：
| 湍流模型 | 精度 | 计算成本 | 备注 |
|---|---|---|---|
| SST $k$-$\omega$ | 优秀 | 中 | 壁面分辨能力优越。最推荐 |
| $k$-$\varepsilon$（Low-Re型） | 优秀 | 中 | 不用壁面函数，一直解到壁面 |
| 标准 $k$-$\varepsilon$ +壁面函数 | 不良 | 低 | 无法准确捕捉自然对流的壁面流 |
| LES（大涡模拟） | 非常高 | 非常高 | 研究用。羽流动态详细分析 |
层流区（$Ra < 10^7$）用层流求解就够了，不需湍流模型。
水平平板自然对流的实务应用
PCB水平配置的芯片冷却设计
实际PCB冷却设计怎么用这个相关式？具体数字能讲讲吗？
好问题，来个具体例子。50 mm × 50 mm的BGA芯片装在水平基板上面，发热2 W，环温25℃。
第1步：特性长度
$L_c = A/P = (0.05 \times 0.05) / (4 \times 0.05) = 0.0125$ m
第2步：膜温度和物性值
表面温度假设65℃ → $T_f = (65+25)/2 = 45℃ = 318$ K
空气物性值（318 K）：$\nu = 1.75 \times 10^{-5}$ m²/s、$\alpha = 2.47 \times 10^{-5}$ m²/s、$k = 0.0274$ W/(m·K)、$\beta = 1/318 = 3.14 \times 10^{-3}$ 1/K
第3步：Ra数
第4步：Nu数和 $h$
$Ra = 5500$（<10^4略超范围但近似）：$Nu = 0.54 \times 5500^{0.25} = 4.65$
$h = 4.65 \times 0.0274 / 0.0125 = 10.2$ W/(m²·K)
第5步：表面温度确认
$T_s = T_\infty + Q/(hA) = 25 + 2 / (10.2 \times 0.0025) = 25 + 78.4 = 103.4℃$
与初始假设65℃不符 → 迭代收敛到 $T_s \approx 95℃$。
95℃很高啊！如果芯片在下面会更糟？
完全同意。下面的话 $Nu = 0.27 \times Ra^{0.25}$ 所以 $h$ 大约是一半的5 W/(m²·K)。表面温度会超过180℃，远超半导体结温上限（通常125～150℃）。自然对流无法冷却，必须加散热片、热导孔、铜嵌入、或设计筐体气流。
密集安装的情况怎样？芯片挨着靠…
好眼力。密集时热尾流（热后流）成问题。上游芯片的热空气流向下游芯片，下游冷却性能衰退。水平板各方都有诱引流所以比竖直排列没那么严重，但芯片间距小于芯片尺寸2倍时Nu下降10～20%有据可查。
太阳热集热器的热损失评估
太阳集热器也用？什么样的？
平板型太阳热集热器（1 m × 2 m）的例子。冬季条件：面板温60℃、外气5℃、$\Delta T = 55$ K。
$L_c = (1 \times 2) / (2 \times (1+2)) = 0.333$ m。膜温 $T_f = 305.5$ K 的空气物性值下：
$Ra \approx 2.8 \times 10^8$（湍流区）。上面（加热向上）的相关式：
$Nu = 0.15 \times (2.8 \times 10^8)^{1/3} = 98.2$
$h = 98.2 \times 0.0267 / 0.333 = 7.87$ W/(m²·K)
对流热损 $Q = h \cdot A \cdot \Delta T = 7.87 \times 2 \times 55 = 866$ W
这就是1000 W/m²日射的43%被损失！所以才要玻璃盖（温室效应）和选择吸收膜来降低对流损失。
验证与妥当性确认
CFD结果怎么验证对不对？
要分3个阶段验证：
- 网格收敛性：3个及以上网格密度，Nu数收敛到2%以内
- 域独立性：扩大计算域，结果无变化
- 与相关式对比：用McAdams式算Nu，CFD结果±15%以内为合理
标准的基准问题是de Vahl Davis（1983年）的正方形腔自然对流。$Ra = 10^3 \sim 10^6$ 的精确参考解公开发表，是验证求解器的绝佳工具。
FloTHERM与水平PCB热设计
Siemens EDA（前身Mentor Graphics）的FloTHERM 12有个功能：根据PCB放置角度自动切换自然对流Nu相关式。Panasonic在2019年用FloTHERM评估家用空调主板，水平配置部件的结温与测实值的误差±4℃以内。这证明了正确的相关式选择与网格设定下，3D CFD可成为实用的设计工具。
水平平板自然对流的软件对比
主流CFD工具的对应情况
水平板自然对流用哪个CFD软件比较好？
主要工具的对比表：
| 工具 | Boussinesq | 完全可压 | 辐射耦合 | 自然对流易用性 |
|---|---|---|---|---|
| Ansys Fluent | ○ | ○ | ○（DO/S2S） | ◎ 设置直观 |
| STAR-CCM+ | ○ | ○ | ○（S2S/DO） | ◎ 多面体网格 |
| COMSOL | ○ | ○ | ○ | ○ 多物理场耦合易 |
| OpenFOAM | ○ | ○ | △（需扩展） | ○ buoyantBoussinesqSimpleFoam |
| FloTHERM | ○ | — | ○ | ◎ 电子设备专用，最小设置 |
| Ansys Icepak | ○ | ○ | ○ | ◎ 电子筐体专用功能丰 |
电子设备冷却的话FloTHERM或Icepak比较省事？
完全同意。FloTHERM导入PCB的Gerber数据，输入芯片功耗就能运行，自然对流包含水平板相关式。但网格和湍流模型控制有限。严谨研究或高精度时Fluent/STAR-CCM+/OpenFOAM更灵活。
OpenFOAM的配置示例
OpenFOAM用哪个求解器？
定常用 buoyantBoussinesqSimpleFoam，非定常用 buoyantBoussinesqPimpleFoam。主要设置项：
transportProperties: 指定 $\beta$（体膨胀系数）、$\nu$（动粘性）、$Pr$、$Pr_t$g文件：设置重力向量如 (0, -9.81, 0)（板水平向）T边界条件：板面fixedValue（恒温）或fixedGradient（恒热流）- 外边界：
inletOutlet+totalPressure允许自然诱引流
湍流模型推荐 kOmegaSST。层流区（$Ra < 10^7$）可设 turbulenceModel 为 laminar。
水平平板自然对流的前沿研究
DNS/LES高精度分析
最新研究用什么方法？
水平加热板自然对流是DNS（直接数值模拟）的典型研究对象。Reese（2021年，J. Fluid Mech.）做了 $Ra = 10^{10}$ 的水平板DNS，详细可视化了湍流羽流产生机制。板端产生的大涡与板中央"细胞结构"相互作用的过程被清晰揭示。
LES（大涡模拟）能以较低成本捕捉羽流的非定常行为。水平板特别适用于：
- 羽流发生间隔和空间分布的统计
- 热尾流的动态行为（PCB密集实装评估）
- 遍历域（$10^6 \sim 10^8$）的精确Nu预测
RANS模型无法捉非定常效应。LES成本是RANS的50～100倍，但GPU求解器（Ansys Fluent GPU版、NekRS等）的发展使其趋向实用。
共轭热传递分析
还要考虑固体导热的情况？
PCB冷却高精度分析必须用共轭热传递（CHT）。板有限导热率，芯片下面和空隙间的温度差异无法用相关式捕捉。FR-4基板（$k \approx 0.3$ W/(m·K)）导热率低，芯片下方容易形成"热斑"。
CHT求解离散固体（基板+芯片）和流体（空气）两个区域，界面施加温度与热流连续条件。Fluent、STAR-CCM+、COMSOL都有CHT标配。
相关式手算和完整CFD+CHT什么时候各自用？
实务的使分：
- 概念设计·初期検讨：McAdams式手算。数分。±20%精度。
- 詳細设计·配置优化：FloTHERM/Icepak的3D CFD。数時間。±5～10%精度。
- 限界设计·規格適合：CHT分析+网格収束+実測比较。数日。±3～5%精度。
全部CFD時間。相関式大温度掴、問題案例CFD進的賢明。
水平平板自然对流的故障排除
常见失败与对策
初学者容易犯什么错？
见过很多典型失败：
| 失败 | 症状 | 对策 |
|---|---|---|
| 计算域太小 | Nu比相关式低20%以上 | 域扩大至板尺寸5倍以上 |
| Boussinesq近似滥用 | $\Delta T > 100$ K时精度恶化 | 改用完全可压缩（理想气体） |
| 用壁面函数 | $h$ 高估30%以上 | 改Low-Re模型或SST $k$-$\omega$ |
| 重力方向错 | 无对流或异常流向 | 再确认 $\mathbf{g}$ 向量 |
| 边界条件矛盾 | 质量守恒失败、发散 | 压力出口改 totalPressure |
| 初值设定不当 | 非物理过渡导致发散 | 流体初值$T_\infty$，壁面$T_s$ |
重力方向错真的常见？
听起来荒唐但太常见！CAD坐标系与CFD坐标系不同，或板倾斜时重力向量符号弄反。自然对流的唯一驱动力就是重力，方向错了整个物理无效。第1优先级检查项目。
收敛困难时的处方
自然对流分析收敛慢的印象…
事实。自然对流比强制对流收敛难。速度小，流体热耦合强，浮力项相对重。对策：
- 低松弛系数：速度0.3～0.5、压力0.2～0.3、温度0.7左右
- 初值精心设置：流体全$T_\infty$、壁面$T_s$。全体设$T_s$会初期浮力过大导致发散
- 非定常求解：定常收不了就用非定常推进。羽流本身是非定常，定常解可能不存在
- 收敛判定：只看残差不行。要监测温度·Nu的变化。残差在$10^{-4}$振荡说明羽流非定常性，不会更低。用时间平均值
"定常解不存在"——想想确实，羽流摇晃的…
完全对。特别$Ra > 10^7$湍流区，物理上不存在定常解。定常求解器的"收敛解"实际是时间平均状态。精确的Nu需非定常求解足够长时间再平均。这才是严谨。
更多