平板强制对流
平板强制对流的理论基础
边界层这个主角
当均匀来流 \( U_\infty \) 沿平板流动时，壁面附近会发展出速度从0变化到 \( U_\infty \) 的速度边界层，以及温度从壁温 \( T_w \) 变化到主流温度 \( T_\infty \) 的温度边界层。换热的阻力几乎全部集中在这层薄薄的区域内，而决定两层厚度之比的是普朗特数 \( Pr = \nu/\alpha \)（空气 \( Pr \approx 0.7 \)：两层厚度相当；油类 \( Pr \gg 1 \)：温度边界层更薄；液态金属 \( Pr \ll 1 \)：温度边界层更厚）。平板是一切对流换热问题的"基准问题"，也是少数能把CFD、实验与关联式三者放在一起对照的构型，因此它同样被用作V&V的验证基准算例。
层流区——Blasius解与局部Nu数
以距前缘距离 \( x \) 构造的雷诺数 \( Re_x = U_\infty x / \nu \) 低于临界值时（平板的经验值 \( Re_{x,c} \approx 5\times10^5 \)，但主流湍流度会使其大幅变动），边界层为层流，由相似解（Blasius解）可导出局部努塞尔数（等温壁、\( Pr \gtrsim 0.6 \)）。
$$ Nu_x = 0.332\, Re_x^{1/2} Pr^{1/3}, \qquad \overline{Nu}_L = 0.664\, Re_L^{1/2} Pr^{1/3} $$
平均值正好是局部值的2倍——换热系数在前缘处最大（边界层最薄），向下游按 \( x^{-1/2} \) 衰减。对等热流壁面，系数0.332要换成0.453。"壁面条件（等温还是等热流）不同，系数就不同"，这是用关联式核对CFD结果时首先要确认的前提。
湍流区与混合平均
转捩之后的湍流边界层换热能力大幅增强，具有代表性的局部式，以及贯通层流区与湍流区求平均的混合平均式如下（前提为 \( Re_{x,c} = 5\times10^5 \)、\( 0.6 < Pr < 60 \)）。
$$ Nu_x = 0.0296\, Re_x^{4/5} Pr^{1/3}, \qquad \overline{Nu}_L = (0.037\, Re_L^{4/5} - 871)\, Pr^{1/3} $$
按惯例，所有物性均取膜温度 \( T_f = (T_w + T_\infty)/2 \) 来评估。即便在理想条件下，关联式的误差也有±10～20%左右，所以不能把"关联式当作真值"，而应视为"有实验支撑的基准带"。
手算与CFD的计算方法
手算实例——空气、0.5m平板
公式我看懂了，可是代进真实数字后h到底是多大，我完全没有概念……
亲手算一遍，这份感觉能受用一辈子。取空气（膜温度按50℃计：\( \nu \approx 1.8\times10^{-5} \) m²/s、\( k \approx 0.028 \) W/m·K、\( Pr = 0.70 \)）、\( U_\infty = 10 \) m/s、板长 \( L = 0.5 \) m。先算 \( Re_L = 10 \times 0.5 / 1.8\times10^{-5} \approx 2.8\times10^5 \)——低于临界值，所以整块板都是层流。平均努塞尔数 \( \overline{Nu}_L = 0.664 \times \sqrt{2.8\times10^5} \times 0.70^{1/3} \approx 0.664 \times 529 \times 0.888 \approx 312 \)，平均换热系数 \( \bar{h} = \overline{Nu}_L\, k / L \approx 312 \times 0.028 / 0.5 \approx 17 \) W/m²K。空气强制对流的h是几十W/m²K的量级——心里有这个尺度，CFD一算出h=500，你立刻就能察觉"哪里不对劲"。
CFD求解时的网格与壁面处理
- 低雷诺数型（壁面积分）——壁面第一层单元取 \( y^+ \le 1 \)，边界层内布15～25层。这是精确求出换热的基本做法
- 壁面函数型——把第一层单元放在 \( 30 \lesssim y^+ \lesssim 300 \)。虽然粗糙，但适合大规模模型。中间的 \( y^+ \approx 5\sim30 \) 最糟糕，两种假设都不成立，换热会严重偏离
- 入口条件——主流湍流度会强烈改变转捩位置。若要与关联式对比，要么把湍流度压低（<1%）以形成层流区，要么使用转捩模型（γ-Reθ等）
- 前缘的分辨率——边界层最薄的前缘附近，在流向上也要加密。前缘处的数值钝化会影响整块板的平均h
作为验证基准算例的用法
平板是"CFD设置健全性测试"中性价比最高的问题。把局部 \( Nu_x \) 分布与关联式叠在一起，确认三点：①层流区是否重现出 \( x^{-1/2} \) 的斜率，②转捩位置是否与预期一致，③湍流区的量值是否落在±10%左右以内。如果带着在这里都对不上的设置（网格、壁面处理、湍流模型）直接进入实际几何，那么实际几何出现偏差时就再也无法切分原因了。
实务应用指南
关联式的三种用途
- 概算设计——散热面积、冷却风量的初步估算。在做CFD之前先把答案的量级抓住
- CFD的验算——把分析得到的平均h与关联式对比，相差2倍以上就要怀疑设置
- 提供边界条件——把 \( h \) 作为不做精细求解的面上的换热系数，交给热传导分析或一维热网络
了解与实物产生偏差的主要原因
| 因素 | 影响 | 对策 |
|---|---|---|
| 主流湍流度偏高（实机为百分之几） | 转捩提前，平均h增大 | 降低混合平均式中的临界Re数重新评估，或保守地采用全湍流式 |
| 存在未加热起始段 | 温度边界层发展滞后，h上升 | 使用带起始段修正的公式 |
| 表面粗糙度 | 使湍流区的h增大（光滑壁公式为下限） | 改用基于摩擦系数的粗糙壁关联式 |
| 有限宽度与端部效应 | 三维流动使h出现局部波动 | 宽长比较小时用CFD确认 |
| 物性的评估温度 | 温差大时物性误差会直接传递 | 严格按膜温度评估，必要时加物性修正系数 |
面对"既非等温也非等热流"的实物
实际的散热面处于等温壁（0.332系）与等热流壁（0.453系）之间。两式之差约为36%——也就是说，仅壁面条件的理想化本身就带来这个量级的不确定幅度。需要精度时，正道是进一步做与固体侧热传导耦合的共轭传热分析，而把关联式保留下来作为耦合结果的合理性检查。
各工具中的处理
CFD工具的设置要点
| 工具 | 要点 |
|---|---|
| Ansys Fluent | Enhanced Wall Treatment（y+自适应壁面处理）。转捩用Transition SST（γ-Reθ）模型。通过壁面热流报告提取局部Nu |
| OpenFOAM | 壁面函数用nutkWallFunction系＋alphatWallFunction。低雷诺数型则用y+<1的网格并去掉壁面函数。后处理用wallHeatFlux工具 |
| STAR-CCM+ | All y+壁面处理。内置转捩模型（Gamma ReTheta）。用派生部件功能提取局部h分布 |
| COMSOL | 用共轭传热接口做流固一体分析。层流区可直接求解 |
| 一维／热网络工具 | 内置平板关联式（自动计算h）。使用前要确认其内嵌假设的适用范围（Re、Pr、壁面条件） |
工具间对比前先统一"h的定义"
换热系数由 \( h = q_w / (T_{ref} - T_w) \) 定义，但参考温度 \( T_{ref} \) 的取法（主流温度、绝热壁温、入口温度）在工具与关联式之间可能并不一致。特别是高速流（黏性加热起作用时）或内外流混合的模型中，定义之差会造成几十个百分点的表观偏差。工具之间以及与关联式的对比，用不依赖定义的壁面热流 \( q_w \) 本身来做才安全。
前沿研究动向
转捩预测的高级化
从层流到湍流的转捩是平板换热最大的不确定因素，仅用局部变量判定转捩的γ-Reθ型模型仍在持续改进，并向敏感环境（湍流度、压力梯度、粗糙度）扩展。基于LES与DNS的转捩机理（旁路转捩、Tollmien-Schlichting波）数据库已经成型，被用作RANS转捩模型的标定基础。从实务角度看，"转捩位置的预测误差＝平均h误差的主因"这一格局并未改变，因此凡涉及转捩的项目，把转捩位置的敏感度分析纳入标准流程是稳妥的做法。
向非理想条件的扩展
在更接近实机的条件下——高湍流度主流、粗糙表面、脉动流、与冲击射流的复合——平板换热数据的积累与关联式的扩展仍在推进。此外，从一开始就包含共轭效应（板内热传导形成壁温分布，进而影响边界层）的共轭Blasius型问题，其解析解与数值解也在完善，摆脱"等温／等热流二选一"的趋势在理论层面同样在推进。
数据驱动的关联式构建
用符号回归或神经网络从实验与CFD数据中构建换热关联式的研究日益增多。它有可能把适用范围拓宽到超越经典关联式，但外推的风险与物理一致性（量纲一致、极限行为）的保证仍是课题，因此"保留无量纲数的框架，只把系数与指数交给数据驱动去优化"这种折中路线，是目前最有希望的实用落点。
故障排查
按症状的原因与对策
| 症状 | 可能原因 | 对策 |
|---|---|---|
| CFD的Nu比层流关联式高出几十个百分点 | 意外发生转捩（湍流度过大）、y+不匹配 | 降低入口湍流度／使用转捩模型。检查y+分布 |
| CFD的Nu远低于关联式 | 边界层分辨不足（层数、前缘）、壁面函数超出适用范围 | 重新设计边界层网格。改用低雷诺数型 |
| 与实验对不上（但CFD与关联式一致） | 实验侧的湍流度、粗糙度、端部效应 | 重新核对实验条件的无量纲数。评估偏差是否落在不确定度范围内 |
| 不同工具之间h相差二三成 | 参考温度的定义差异 | 改用热流对比。查文档确认定义 |
| 局部Nu分布在前缘处发散式增大 | 前缘奇异性（边界层厚度→0）的离散化 | 把前缘附近的数值排除在评价之外（理论上也有Nu_x→∞） |
| 手算与结果的量级对不上 | 物性的评估温度有误、特征长度取错（x与L混淆） | 按膜温度重新评估物性。确认是否把局部式与平均式混用了 |
拿不准时，先怀疑转捩
当CFD、关联式和实验给出三个不同的值时，该从哪里下手呢？
平板对不上时的头号嫌疑就是转捩位置。层流与湍流的h相差数倍，如果三方对"转捩发生在哪里"的认识都不统一，数值对不上是理所当然的。诊断要靠局部分布——CFD看壁面热流沿 \( x \) 方向的分布，实验看壁温分布，找出h急剧上升的位置（即转捩）。把三方的转捩位置对齐之后再比一次，残差通常就收敛到±10～20%这种"正常的不确定度"范围内了。在比较数值之前，先比较分布的形状——这是对流换热验证的铁律。
帮助
更多
错误