研究报告：定义物理可实现的内流自然对流基准问题
一、 作者、机构及发表信息
本研究由 W. H. Leong, K. G. T. Hollands 和 A. P. Brunger 完成。其中，K. G. T. Hollands 为通讯作者。主要工作单位是加拿大滑铁卢大学机械工程系的太阳能热研究实验室 (Solar Thermal Research Laboratory, Department of Mechanical Engineering, University of Waterloo)。该研究发表于国际知名学术期刊《International Journal of Heat and Mass Transfer》的第30卷，出版时间为1987年。
二、 学术背景与研究目标
本研究属于计算流体力学 (Computational Fluid Dynamics, CFD) 验证与实验传热学领域。研究的核心动机源于CFD领域对“基准问题”的迫切需求。基准问题指的是一个得到广泛认可的、定义清晰的标准案例，用于验证和评估CFD代码的准确性和可靠性。一个理想的基准问题应易于陈述、能挑战复杂代码，且物理上可实现，即能在实验室中被精确复现和测量。
当时流行的自然对流基准问题——由 de Vahl Davis 提出的二维方腔问题（侧壁等温、上下壁面绝热）——存在物理不可实现的缺陷。首先，对于空气填充的空腔，完全绝热边界条件实际上无法实现。其次，假定的二维流动在现实中可能不稳定并发展为三维流动，这削弱了其作为基准问题的物理意义。尽管存在这些缺陷，该问题因其数学上的简洁性仍在CFD发展中扮演了重要角色。然而，随着CFD技术的成熟，开发一个物理上可实现的、可被实验精确检验的基准问题变得至关重要，尤其是在过渡流和湍流领域，湍流模型的准确性亟需可靠的实验数据进行验证。
因此，本研究旨在定义一个全新的、物理上可实现的内流自然对流基准问题，并设计建造相应的实验装置，以极高的精度（目标为1%）测量关键参数（努塞尔数，Nusselt Number），从而为CFD代码验证提供可靠的实验数据。最终目标是建立一个可以同时通过CFD模拟和物理实验进行严格验证的标准案例。
三、 详细工作流程
本研究的工作流程可分为四个主要阶段：基准问题的定义、实验装置的设计与构建、测量与数据采集、以及不确定度分析与CFD验证。
1. 基准问题的定义 研究团队决定在经典的二维方腔问题基础上进行改进，使其物理可实现。他们定义了以下核心特征： * 几何形状：将二维方形扩展为立方体腔，以容纳可能的三维流动，这更符合物理现实。 * 边界条件： * 一对相对的壁面保持为等温壁面，分别为热壁面 (Th) 和冷壁面 (Tc)。 * 其余四个壁面的温度设定为从Tc到Th的线性变化分布。这是关键改进，因为线性温度分布在实验室中比绝热条件更容易精确实现（通过精心设计的传导路径）。 * 该条件在CFD中也被称为“完美导热侧壁”条件。 * 流动工质：选择空气作为工作流体。原因包括：空气的热物性数据精确已知；通过改变空气压力即可在大范围内调节瑞利数 (Rayleigh Number, Ra)，而无需改变模型或温差；普朗特数 (Prandtl Number) 对温度不敏感。 * 研究工况：定义了三个倾斜角度作为不同的基准问题：0°（热面在下方，底部加热）、45°（倾斜加热）和90°（热面在侧面，侧向加热）。这三个角度覆盖了从稳定分层到复杂流动的不同物理机制。 * 测量参数：选择平均努塞尔数 (Nu) 作为核心测量参数。因为它反映了整体的换热性能，对全局流动响应敏感，且具有悠久的测量历史，有望达到高精度。
2. 实验装置的设计与构建 为实现上述定义的问题，研究人员设计并建造了一个精密的实验装置。 * 整体结构：立方体腔的边长 (L) 名义值为152.4 mm（6英寸）。装置核心包括：一个热流计、一个电加热板（同时作为冷板）以及两个匹配的半立方体（“半立方”）组件。两个半立方体沿腔体对角线平面拼接，形成完整的立方体。热板和背板通过焊接在背面的管道由独立的循环水系统加热或冷却。 * 线性温度分布的实现： * 侧壁设计与厚度：侧壁采用黄铜制造，其厚度 (t = 1.59 mm) 经过初步分析确定，以确保侧壁内部温度梯度极小，从而在其内表面（朝向空气）获得所需的线性温度分布。计算表明，最大偏离线性分布的温度差约为总温差 (ΔT) 的0.5%。 * 接触热阻的解决：初期测试发现，在侧壁与热板/背板的连接处存在显著的接触热阻，导致温度跳跃。为解决此问题： 1. 方案一：调换了冷热水流。并将冷板（电加热板）的温度调节至与侧壁-背板连接处的温度 (Tjc) 一致，消除了该处的温度跳跃。 2. 方案二：在侧壁与热板的连接处安装了一个辅助加热器（铜带嵌有镍铬丝），以补偿该处的热损失，成功将温度跳跃减小到总ΔT的0.1%以下。 * 瑞利数调节：为实现宽范围的瑞利数，整个模型被放置在一个压力容器内。通过调节容器内空气的压力（范围从0.1 kPa到100 kPa），可以在不改变几何尺寸或温差的情况下，使瑞利数在大约10^3到10^7的范围内变化。 * 热流量测量技术：采用了一种改进的“混合方法”测量通过腔体的总热流。核心公式为： q_conv = [(VI)_s - (VI)] - α (E_s - E) 其中，q_conv 是对流热流；(VI)_s 和 E_s 是模型处于“加热在上” (θ=180°) 静止状态（仅有导热和辐射）时测量的电功率和热流计电动势；(VI) 和 E 是目标倾角下测量的值；α 是热流计的标定常数。通过测量静止状态下的热流并从中减去，可以精确分离出纯对流贡献。 * 温度测量：使用铂电阻温度计和T型热电偶精确测量各壁面温度。热板和背板各嵌入两个铂电阻温度计，其温度取平均值。侧壁上布置了12个热电偶以监测其温度分布是否满足线性要求。
3. 测量与数据分析流程 * 实验步骤：设定好压力、倾角、壁面温度等参数后，首先进行热流计常数 α 的现场标定。然后，计算机程序自动调整冷板的电加热功率，直到冷板温度 (Tc) 与侧壁-背板连接处温度 (Tjc) 相等（控制在±0.001 K内）。达到稳态后，进行长达半小时的数据采集（共401个样本点），记录电压V、电流I、热流计输出E、所有温度及压力，并取平均值。 * 数据处理： * 努塞尔数 (Nu)：根据公式 (4) 计算，其中对流热流 q_conv 由上述混合方法获得，热导率k考虑了随温度和压力的变化。 * 瑞利数 (Ra)：根据公式 (5) 计算，其中包含了修正函数 f_k(T, p) 和 f_m(T, p) 来考虑空气物性随温度和压力的变化。 * 不确定度分析：采用Moffat的多样本法进行误差分析，并对该方法进行了适当修正以适应实验中Ra值围绕一个标称值（如3×10^5）变化的情况。总不确定度 U_y 是随机误差 e_r 和系统（偏差）误差 e_b 的均方根。表1详细列出了所有输入量（如V, I, ΔT, p, L等）的偏差极限。对于在特定Ra（记作Ra*）下报告的Nu（记作Nu*），还考虑了Nu随Ra变化的斜率 n 带来的误差。分析显示，实验测量Nu的总不确定度在0.5%至0.8%之间，达到了预设的1%精度目标。
4. CFD模拟验证 为了验证实验装置的准确性，研究团队使用商业CFD软件TASCFlow2D对Ra = 3×10^5 的三种倾角工况进行了数值模拟。 * 模拟细节：利用腔体关于中截面的对称性，只模拟了一半空间。采用了三种不同密度的网格（15×6×14, 21×11×21, 31×15×31），并应用重复理查森外推法 (Repeated Richardson Extrapolation) 来估计网格无限加密时的“精确”解，并以最密网格结果与外推结果差值的一半作为数值误差。 * 模拟设置：进行了两类模拟： * 常物性模拟：物性参数基于平均温度Tm取值。 * 变物性模拟：局部计算物性（热导率、粘度、比热容随温度变化，密度和体积膨胀系数遵循理想气体定律）。模拟中使用的Tm、ΔT和p值是根据多组实验的平均值设定的。
四、 主要研究结果
1. 侧壁温度线性分布的验证 图2展示了在Ra≈10^6、倾角90°条件下侧壁温度的测量结果。实验数据点与理论直线吻合极好，偏差在测量不确定度范围内（约±0.01 K）。在其他Ra和倾角下也观察到了类似的线性度，证明了实验装置成功实现了所定义的线性温度边界条件。
2. 临界瑞利数的验证 在0°倾角（底部加热）下，测量了流动开始失稳的临界瑞利数 Rac。对于线性温度侧壁边界条件的理论值是1708 (Catton)。实验在两种ΔT（2.4 K和4.3 K）下进行，测得的Rac分别为1708±33和1704±28。实验值与理论值高度吻合，这强有力地验证了侧壁边界条件满足设计要求，因为Rac对侧壁条件非常敏感（例如，绝热侧壁的Rac为1708/8=214）。
3. 误差分析与基准实验数据 在Ra* = 3×10^5下，对三个倾角（0°, 45°, 90°）进行了共计37组独立实验。结果汇总于表1。 * 测量得到的平均Nu*分别为：1.107 (0°)， 1.450 (45°)， 1.226 (90°)。 * 随机误差 e_r 在0.1%到0.6%之间，系统误差 e_b 在0.5%到0.6%之间。由此计算的总不确定度 U_{Nu} 在0.6%到0.8%之间，成功实现了低于1%的精度目标。
4. 与CFD模拟结果的对比 * 与常物性CFD模拟对比：模拟结果（表1第5项）与实验值接近，但除45°倾角外，其余两个倾角的差异超过了实验不确定度范围（比较比值 >1）。 * 与变物性CFD模拟对比：考虑到实验中ΔT/Tm比值（约1%）的影响，进行了变物性模拟。模拟结果（表1第7项）与实验值的差异显著缩小。差异百分比分别为：0°倾角 -0.34%， 45°倾角 -0.16%， 90°倾角 -0.02%。所有差异均小于实验不确定度（比较比值 ）。 * 这表明，为了在实验误差范围内与实验结果一致，CFD模拟必须考虑流体物性随温度的变化，即使温差相对于绝对温度很小（约1%）。
五、 研究结论与价值
本研究成功地定义并物理实现了一个用于验证CFD代码的内部自然对流基准问题。该问题基于一个具有线性温度分布侧壁的立方体空气腔，定义了三种倾角，使其成为物理上可实现的简单标准案例。
实验装置能够高精度地建立所需的边界条件，并测量平均努塞尔数。在Ra = 3×10^5下，测量不确定度达到约0.7%，测量结果与高精度CFD模拟（考虑变物性）的差异平均仅为0.2%，验证了实验装置的可靠性。临界瑞利数的测量结果也与理论值高度吻合，进一步证实了边界条件的准确性。
该研究的科学价值在于为CFD社区提供了一个严格、可复现的基准问题，弥补了传统基准问题物理不可实现的缺陷。其实验数据，特别是高精度的Nu测量值，为CFD代码，尤其是涉及变物性效应和复杂三维流动的代码，提供了关键的验证依据。其应用价值体现在推动CFD模型，特别是自然对流湍流模型，向更高精度发展。研究表明，即使温差很小，物性变化的影响也不容忽视，这对CFD实践具有重要指导意义。
六、 研究亮点
七、 其他有价值的内容
论文在讨论部分还指出了使用空气作为工质的其他优势：利用单一实验模型，仅通过改变压力即可覆盖极其宽广的瑞利数范围（从约10^3到10^7）。这为研究宽参数范围内的自然对流现象提供了极大的便利。文中提到，在更高瑞利数（10^4, 10^5, 10^6, 10^7）下的努塞尔数结果将在后续的姊妹论文中发表，预示着该基准问题将提供一套覆盖层流到湍流过渡区的完整实验数据库，对CFD代码的验证能力提出了更全面的挑战。此外，附录部分详细描述了CFD模拟的设置（网格策略、收敛准则、物性处理方法、外推技术），为其他研究者重复该数值基准提供了足够的信息。