水平壁面自然对流（加热/冷却）分析
水平壁面自然对流是热对流学中的基础模型，其驱动力源于流体密度随温度变化产生的浮力效应。当水平放置的固体壁面与周围流体存在温差时，近壁面流体受热或冷却，导致密度改变。若壁面温度高于流体温度（加热情况），近壁面流体密度降低，在重力作用下向上运动，形成上升流；若壁面温度低于流体温度（冷却情况），近壁面流体密度增加，向下沉降，形成下降流。这种由温度梯度引发的流体宏观运动，构成了自然对流的核心特征。
控制方程与无量纲参数
分析水平壁面自然对流，通常基于纳维-斯托克斯方程（Navier-Stokes equations）和能量方程，并引入 Boussinesq 近似。该近似假设流体密度仅在重力项中随温度线性变化，而在其他项中视为常数，从而简化了数学处理。
描述该流动状态的关键无量纲参数包括：
- 格拉晓夫数 ($Gr$)：表征浮力与粘性力之比。对于水平壁面，$Gr_L = \frac{g \beta \Delta T L^3}{\nu^2}$，其中 $g$ 为重力加速度，$\beta$ 为体积膨胀系数，$\Delta T$ 为壁面与流体的温差，$L$ 为特征长度，$\nu$ 为运动粘度。
- 普朗特数 ($Pr$)：表征动量扩散与热扩散之比，反映了流体物性对流动和传热的影响。
- 努塞尔数 ($Nu$)：表征对流传热与导热之比，是计算对流换热系数的核心参数。
流动边界层结构分析
在水平壁面自然对流中，边界层的发展具有显著的方向性特征。
加热壁面情形
当壁面被加热时，流体受热上升。在壁面入口端（前缘），边界层较薄，流速较低。随着流体沿壁面向下游移动，浮力持续加速流体，边界层厚度逐渐增加，而局部流速先增大后趋于稳定。由于浮力作用，速度剖面在远离壁面处出现峰值，随后因粘性耗散而减小。这种流动结构导致局部努塞尔数沿流动方向逐渐减小，因为边界层增厚降低了温度梯度。
冷却壁面情形
当壁面被冷却时，流体受冷下沉。与加热情形类似，边界层从入口端开始发展并逐渐增厚。然而，由于流体向下运动，流动稳定性受到重力方向的影响。在大多数气体和液体中，冷却壁面的流动结构与加热壁面具有相似的对数律特征，但具体的速度分布和温度分布会因 $Pr$ 数不同而有所差异。值得注意的是，对于高 $Pr$ 数流体，动量边界层远薄于热边界层；而对于低 $Pr$ 数流体（如液态金属），热边界层则远薄于动量边界层。
传热关联式与工程应用
在工程实践中，直接求解偏微分方程成本较高，因此常采用经验或半经验关联式来预测平均努塞尔数。对于水平平板，常用的关联式取决于流动状态（层流或湍流）以及 $Ra$ 数（瑞利数，$Ra = Gr \cdot Pr$）的范围。
例如，对于层流状态（$10^4 < Ra_L < 10^9$），平均努塞尔数可近似表示为：
$$ \overline{Nu}_L = 0.665 + \frac{0.357 Ra_L^{1/5}}{[1 + (0.492/Pr)^{9/16}]^{8/15}} $$
该公式适用于 $0.6 < Pr < 10^4$ 的范围。在计算时，需根据壁面平均温度确定流体物性。
影响因素与数值模拟建议
影响水平壁面自然对流换热的因素主要包括：
- 温差大小：温差直接决定浮力强度，进而影响 $Gr$ 数和流动速度。
- 流体物性：$Pr$ 数决定了热边界层与速度边界层的相对厚度，对传热效率有显著影响。
- 几何尺寸：特征长度 $L$ 影响边界层的发展程度，较长的壁面会导致更厚的边界层和更低的平均换热系数。
在进行数值模拟（CFD）时，建议采用 $k-\epsilon$ 或 $k-\omega$ 湍流模型处理湍流区域，并启用浮力项。网格划分需在壁面附近足够精细，以准确捕捉边界层内的梯度和流动细节。此外，验证模拟结果时，应确保网格无关性，并将计算得到的局部和平均努塞尔数与上述经典关联式进行对比，以评估模型的准确性。
总结
水平壁面自然对流是理解复杂热对流现象的基础。通过深入分析其物理机制、控制方程及边界层结构，工程师可以更准确地预测换热性能。在实际应用中，结合无量纲参数和经典关联式，能够高效地解决大多数工程问题，而数值模拟则为非标准几何或复杂工况提供了强大的分析工具。