浮力驱动流动的数学描述
在热对流的研究中，浮力驱动流动（Buoyancy-driven flow），通常被称为自然对流（Natural Convection），其核心动力源于流体内部由于温度梯度引起的密度差异。与机械装置（如风扇或泵）驱动的强制对流不同，浮力驱动流动的运动完全由流体自身的物理性质变化所驱动。为了精确描述这一复杂的物理过程，我们需要建立一套基于流体力学基本定律的数学模型。
浮力驱动流动的数学描述建立在质量守恒、动量守恒和能量守恒三大定律之上。在连续介质假设下，这些定律表现为偏微分方程组。
1. 质量守恒方程（连续性方程）
对于不可压缩流体（或在密度变化极小的情况下），质量守恒要求流体在空间中任何一点的流入量必须等于流出量。其数学表达式为：
$$\nabla \cdot \mathbf{u} = 0$$
其中，$\mathbf{u}$ 是流体的速度矢量。
2. 动量守恒方程（Navier-Stokes 方程）
动量方程描述了流体受力与运动之间的关系。在浮力驱动流动中，除了压力梯度力和粘性力外，最关键的是由于密度不均匀产生的重力项（浮力项）：
$$\rho \left( \frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} \right) = -\nabla p + \mu \nabla^2 \mathbf{u} + \rho \mathbf{g}$$
其中：
- $\rho$ 是流体密度；
- $p$ 是压力；
- $\mu$ 是动力粘度；
- $\mathbf{g}$ 是重力加速度矢量。
3. 能量守恒方程
由于浮力是由温度变化引起的，因此必须引入能量方程来描述温度场的变化：
$$\rho c_p \left( \frac{\partial T}{\partial t} + \mathbf{u} \cdot \nabla T \right) = k \nabla^2 T$$
其中：
- $T$ 是温度；
- $c_p$ 是定压比热容；
- $k$ 是热导率。
Boussinesq 近似：数学描述的核心简化
在实际的数学建模中，直接求解上述方程组非常困难，因为密度 $\rho$ 是温度 $T$ 的函数，这使得方程组呈现高度的非线性。为了简化计算，工程和物理学中广泛采用 Boussinesq 近似。
Boussinesq 近似的核心思想是：在流体的连续性方程和动量方程的大部分项中，将密度视为常数 $\rho_0$；仅在与重力相关的浮力项中，保留由于温度变化引起的密度微小波动。
根据热胀冷缩原理，密度与温度的关系可以线性化表示为：
$$\rho = \rho_0 [1 - \beta(T - T_0)]$$
其中：
- $\rho_0$ 是参考温度 $T_0$ 下的密度；
- $\beta$ 是流体的体胀系数。
将此关系代入动量方程的重力项中，得到简化后的浮力项：
$$\rho \mathbf{g} \approx \rho_0 \mathbf{g} - \rho_0 \beta (T - T_0) \mathbf{g}$$
由于 $\rho_0 \mathbf{g}$ 项可以合并到压力梯度项 $-\nabla p$ 中，最终的简化动量方程变为：
$$\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} = -\frac{1}{\rho_0} \nabla p' + \nu \nabla^2 \mathbf{u} + \beta (T - T_0) \mathbf{g}$$
其中 $p'$ 是修正后的压力。这种近似在温差较小（$\Delta T \ll T_0$）且流体接近不可压缩的场景下（如空气、水）具有极高的准确性。
无量纲化与特征参数
为了深入理解浮力驱动流动的特性，数学描述通常会通过无量纲化处理，引入几个关键的无量纲数。这些参数决定了流动的状态（是层流还是湍流）以及热传递的效率。
1. Grashof 数 ($Gr$)
Grashof 数反映了浮力驱动力与粘性阻力的比值，类似于强制对流中的雷诺数：
$$Gr = \frac{g \beta (T_s - T_\infty) L^3}{\nu^2}$$
其中 $L$ 是特征长度，$\nu$ 是运动粘度，$T_s$ 是壁面温度，$T_\infty$ 是环境温度。
2. Prandtl 数 ($Pr$)
Prandtl 数描述了动量扩散能力（粘性）与热扩散能力（热传导）的比值：
$$Pr = \frac{\nu}{\alpha} = \frac{\mu c_p}{k}$$
其中 $\alpha$ 是热扩散率。
3. Rayleigh 数 ($Ra$)
Rayleigh 数是自然对流研究中最核心的参数，它定义为 $Ra = Gr \cdot Pr$。它综合考虑了浮力、粘性和热扩散的影响：
$$Ra = \frac{g \beta (T_s - T_\infty) L^3}{\nu \alpha}$$
- 低 $Ra$ 数：浮力不足以克服粘性和热扩散，流动表现为稳定的层流。
- 高 $Ra$ 数：浮力占据主导，流动会发生失稳，转变为复杂的湍流状态。
典型应用示例：垂直加热平板的流动
为了直观理解数学描述的应用，我们考虑一个经典的物理模型：一个垂直放置的、温度为 $T_s$ 的无限大平板，周围环绕着温度为 $T_\infty$ 的流体。
- 物理过程：靠近平板的流体受热，密度减小，产生向上的浮力。
- 边界层形成：由于粘性作用，流体在靠近平板处速度较慢，形成速度边界层；由于热传导作用，形成温度边界层。
- 数学求解方向：
- 在边界层内，由于 $L$ 较大，我们可以简化 Navier-Stokes 方程，忽略曲率项。
- 通过引入无量纲坐标 $\eta = y / \delta$（$\delta$ 为边界层厚度），将偏微分方程转化为常微分方程。
- 最终可以推导出 Nusselt 数 ($Nu$) 与 Rayleigh 数 ($Ra$) 之间的经验关系式，例如：$Nu = C \cdot Ra^n$。
通过这种数学描述，工程师可以无需进行大规模实验，仅通过计算 $Ra$ 数的大小，就能预判散热设备的散热性能，并指导散热片的几何设计。
总结
浮力驱动流动的数学描述是一个从物理现象到控制方程，再到简化近似，最后进行无量纲化分析的过程。通过 Boussinesq 近似，我们将复杂的密度变化问题转化为温度场与速度场的耦合问题；而通过 Rayleigh 数，我们将复杂的流动特性浓缩为一个关键的数值指标。掌握这些数学工具，是进行热对流数值模拟（CFD）和工程设计的基石。