量纲分析与π定理
量纲分析与π定理的理论基础
量纲分析的基本原理
老师，我对量纲分析的认识仅限于"检查单位"，在CFD中这么重要吗？
不仅如此。量纲分析是最小化支配物理问题的独立参数数量的强大方法。Buckingham的π定理是其核心。
π定理的论述如下。如果一个物理现象由 $n$ 个物理量描述，这些物理量包含 $k$ 个基本量纲（质量M、长度L、时间T等），则该现象可以由 $n - k$ 个独立无量纲数完全描述。
具体是什么意思？
比如考虑圆管内的压力损失 $\Delta p$。相关的物理量有 $\Delta p, \rho, U, D, \mu, L, \varepsilon$ 这7个。基本量纲是M、L、T这3个，所以 $7 - 3 = 4$ 个无量纲数可以表示。
这里 $\text{Re} = \rho U D / \mu$ 是雷诺数，$L/D$ 是长径比，$\varepsilon/D$ 是相对粗糙度。
从7个变量减到4个，确实很大！
主要无量纲数及其物理意义
我们来整理CFD中常见的无量纲数。每一个都表示"两种力或效应的比"。
| 无量纲数 | 定义 | 物理意义 |
|---|---|---|
| 雷诺数 Re | $\rho UL/\mu$ | 惯性力/粘性力 |
| 马赫数 Ma | $U/c$ | 流速/音速（可压缩性指标） |
| 弗劳德数 Fr | $U/\sqrt{gL}$ | 惯性力/重力 |
| 韦伯数 We | $\rho U^2 L/\sigma$ | 惯性力/表面张力 |
| 斯特劳哈尔数 St | $fL/U$ | 非定常振荡的特性 |
| 普朗特数 Pr | $\nu/\alpha = c_p\mu/k$ | 动量扩散/热扩散 |
| 努塞尔特数 Nu | $hL/k$ | 对流换热/导热 |
| 格拉晓夫数 Gr | $g\beta\Delta T L^3/\nu^2$ | 浮力/粘性力 |
我知道雷诺数，但没想到有这么多。
根据问题选择合适的无量纲数是量纲分析的关键。例如自然对流中，瑞利数 $\text{Ra} = \text{Gr} \cdot \text{Pr}$ 是支配参数。
Navier-Stokes方程的无量纲化
把量纲分析应用到Navier-Stokes方程会怎样？
用代表长度 $L$、代表速度 $U$、代表时间 $L/U$、代表压力 $\rho U^2$ 进行无量纲化，非压缩性N-S方程变成这样：
从这个式子可以看出，Re越大粘性项的贡献越小，Fr越大重力项的影响越小。CFD的解最终由这些无量纲数决定。
所以在风洞试验中匹配Re数很重要！
正是这样。几何形状相似且Re相同时，流场的无量纲解是相同的。这就是动态相似律。
弗劳德的模型船实验——尺度律的发现
19世纪，造船工程师威廉·弗劳德解决了"能否从模型船的阻力测量预测实船性能"这一问题。他发现的规律是"弗劳德数（Fr = V/√(gL)）相等时波形在几何上相似"，这是π定理原型之一。现代船舶设计仍然以模型水池实验和CFD的结合为基础，弗劳德数的一致是必要条件。没有量纲分析，就没有从1/100缩尺模型实验结果推外推到实船的科学依据。
量纲分析与π定理的数值计算方法
CFD中无量纲化的实现
CFD中怎样使用量纲分析？求解器的输入还是有量纲的吧？
好问题。这里有两种方法。
方法1：有量纲计算，后处理时无量纲化
- 物理量直接用SI单位输入
- 结果通过代表量进行无量纲化（例：$C_p = \Delta p / (\frac{1}{2}\rho U_\infty^2)$）
- Ansys Fluent和STAR-CCM+主要采用这种方式
方法2：从一开始就求解无量纲化的方程
- 用代表量正规化的变量输入
- 研究用代码或OpenFOAM的自定义求解器常用这种方式
- 优势是只需设定Re数作为参数
商业软件一般用有量纲的方式啊。
利用相似律的标度化方法
量纲分析最大的应用就是通过相似律进行标度化。确定实机与模型间该匹配的无量纲数，由此决定实验条件。
比如汽车风洞试验用1/5的缩尺模型。
用同样的空气（相同的 $\rho, \mu$），需要 $U_m = 5 \times U_f$。实车100 km/h对应模型500 km/h，会引起可压缩性效应不可忽视。
完全相似也很难实现啊。
对。不能实现完全相似时用部分相似。只匹配支配的无量纲数，对影响较小的无量纲数做妥协。船舶中匹配弗劳德数（造波阻力），不匹配的Re数用补正公式处理是典型例子。
CFD中无量纲数的计算与监测
从CFD结果计算无量纲数时有什么注意的？
代表量的选择会影响结果的解释。下面这些要点很重要。
| 无量纲量 | 代表量的选择 | 注意事项 |
|---|---|---|
| 阻力系数 $C_D$ | 前向投影面积 $A$、$U_\infty$ | 面积定义要统一 |
| 压力系数 $C_p$ | 自由流速 $U_\infty$、静压 $p_\infty$ | 基准压力的选择 |
| 摩擦系数 $C_f$ | 壁面剪切应力 $\tau_w$ | 局部值 vs. 平均值的区别 |
| $y^+$ | $u_\tau = \sqrt{\tau_w/\rho}$ | 壁面第一层网格的质量指标 |
特别是 $y^+$ 是判断湍流模型适用性的最重要无量纲数。壁函数模型时 $30 < y^+ < 300$，壁面分辨模型时 $y^+ \approx 1$ 是目标。
$y^+$ 在这里与量纲分析关联，我没想到。
无量纲数是"大幅减少计算次数"的工具
量纲分析的实际威力在于"参数数量的削减"。比如圆管内压力损失取决于"流速V、密度ρ、粘度μ、管径D、长度L、粗糙度ε"这6个变量，但用π定理可以压缩成"Re数和ε/D（相对粗糙度）"这2个变量的函数，这就是穆迪图。6个变量各独立变化需6³=216次试验，无量纲整理后只需数十次。CFD参数化研究中也是同样道理，先用无量纲数整理就能大幅削减计算情况数。
量纲分析与π定理的实际应用
π定理的应用步骤
怎样用π定理导出无量纲数，教我步骤！
分步骤讲。以热交换器的换热问题为例。
Step 1：列出相关物理量
换热系数 $h$、管径 $D$、流速 $U$、密度 $\rho$、粘度 $\mu$、热导率 $k$、比热 $c_p$，共7个。
Step 2：确认基本量纲
M、L、T、$\Theta$（温度）4个。
Step 3：确定无量纲数个数
$7 - 4 = 3$ 个π群。
Step 4：选择重复变量
选 $\rho, U, D, k$ 4个（与基本量纲个数相同）。
Step 5：构造各π群
Dittus-Boelter关联式 $\text{Nu} = 0.023\,\text{Re}^{0.8}\,\text{Pr}^{0.4}$ 就是从这儿来的！
对。π定理确定无量纲数的"形式"，实验或数值计算确定系数和指数。
CFD分析中的单位系管理
单位系弄错过导致计算结果不对……
常见的失误。注意这些问题。
| 问题模式 | 具体例子 | 影响 |
|---|---|---|
| SI/CGS混用 | 粘度 Pa·s 和 poise 混淆 | Re数偏差100倍 |
| mm单位形状+SI物性 | 长度mm、密度kg/m³ | 压力偏差10⁶倍 |
| 角度 deg/rad | 湍流强度输入 | 边界条件完全不同 |
| 温度 K/C | 浮力计算中的ΔT | Gr数不正确 |
Ansys Fluent用mm单位导入形状的话，要么把密度改成 tonne/mm³、粘度改成 MPa·s，要么把形状转成m单位再导入。
Fluent用mm形状时物性值要转换啊。
验证用基准问题
量纲分析理解的确认用基准问题。
1. 圆柱周围流动（Re依赖性确认）
- Re = 40：稳态孪生涡，$C_D \approx 1.5$
- Re = 100：卡门涡街形成，$\text{St} \approx 0.16$
- Re = 1000：湍流转变开始
- 不同尺度（$D$ 和 $U$ 组合）下匹配Re数，验证得到相同 $C_D$
2. 管内层流（Hagen-Poiseuille流动）
- 理论解：$f = 64/\text{Re}$（Darcy摩擦系数）
- 助流区长度：$L_e/D \approx 0.06\,\text{Re}$
- 验证CFD结果是否符合这个关系
用有理论解的问题验证，确认量纲分析预测的依赖性是否再现就行了。
单位检验是最强的调试方法
遇到"计算结果不对"时，先检查单位是工程师的铁则。从量纲分析的角度讲，方程式两边量纲不一致就一定有计算错误。反过来说"量纲对了不代表计算对，但量纲不对就肯定错了"。实务中常见是SI和英制单位混用。1999年火星气候轨道器撞上火星大气烧毁事件原因就是牛顿·秒和磅力·秒的单位混淆。量纲分析涉及生命。
量纲分析与π定理的软件比较
主要CFD工具中无量纲数的处理
各个CFD软件怎样处理无量纲数？
商业软件的无量纲数参考值设定和报告功能不同。比较一下。
ANSYS Fluent
Fluent用 Reference Values 面板设定代表量（速度、长度、面积、密度、温度等）。这个设定是 $C_D$、$C_L$、$C_p$、Nu等无量纲量计算的基准。
- Report > Forces 输出阻力·升力系数
- Surface Integrals 确认壁面 $y^+$ 分布
- Custom Field Function 可定义任意无量纲数
- TUI命令：
/report/reference-values
ANSYS CFX
CFX用CEL（CFX Expression Language）直接定义无量纲数。
Re = density * velocity * length / viscosity这样的式子可作变量注册- 监测点能追踪无量纲数的收敛历史
- 耦合求解器压力基准选择对结果影响小
Simcenter STAR-CCM+
STAR-CCM+用Field Function定义无量纲数。
- 内置 $C_D$、$C_L$ 报告功能
- Dimensionless Wall Distance（$y^+$）默认计算
- 参数化研究中改变Re数的设计试验（DOE）容易
- Java宏可实现自动化
OpenFOAM
OpenFOAM用 functionObjects 计算无量纲数。
forceCoeffs：计算阻力·升力系数。在controlDict中指定lRef、Aref、magUInfyPlus：壁面 $y^+$ 计算（postProcess -func yPlus）fieldAverage：时间平均的无量纲量codedfunctionObject 能任意定义自定义无量纲数
单位系设置比较
各软件的单位系处理有什么不同？
| 工具 | 单位系 | 注意事项 |
|---|---|---|
| Fluent | SI/CGS/自定义可选 | 内部全部SI。输入时转换 |
| CFX | SI固定 | 输入时用CEL指定有单位的数值 |
| STAR-CCM+ | SI/自定义可选 | Units Manager全体管理 |
| OpenFOAM | 隐含SI | 有单位检验功能（dimensionSet） |
OpenFOAM的 dimensionSet 独特，用 [kg m s K mol A cd] 的7个指数管理量纲。比如粘度是 [0 2 -1 0 0 0 0]（动粘度 $\nu$ 情况）。量纲不匹配的演算编译时出错，这是防虫机制。
OpenFOAM的量纲检验严格，防虫效果好啊。
利用相似律进行参数化研究
量纲分析的实务价值在于参数化研究的效率化。
比如要用多种流量和流体评价配管的压力损失，有量纲参数的组合会非常多，但整理成Re数和相对粗糙度 $\varepsilon/D$ 这2个参数，穆迪图这样1页图就搞定了。
穆迪图本身就是量纲分析的产物啊！
对。CFD也是同样想法，在无量纲数空间中设计参数化研究的话，计算情况数能大幅减少。STAR-CCM+的Design Manager或modeller连携特别适合这个用途。
商用代码的"无量纲化"设定要当心
Fluent和StarCCM+内部不是用无量纲化形式计算，而是用SI单位计算。而研究代码或部分专门求解器是无量纲化后计算，输出也是无量纲值。这个区别不清楚的话，比较结果时会发现数值完全对不上。另外，CFD验证论文中常用的"基准Re数"在商用代码中重现，需要明确指定物性值·速度·形状寸法来匹配Re数。"我要Re=1000" 这样直接输入的地方没有。
量纲分析与π定理的前沿研究
不完全相似与补正方法
不能实现完全动态相似的情况怎么办？
实务中完全相似不能取得的情况更多。这些应对法有几种。
1. 部分相似 + 经验补正
船舶分野的弗劳德法典型。模型试验中匹配弗劳德数，Re数不匹配造成的摩擦阻力差用ITTC 1957的摩擦线补正。
2. 变形模型（Distorted model）
河川模型这样，水平·竖直缩尺比变化的方法。弗劳德数保存但几何相似破裂。
3. 代替流体的使用
航空领域用极低温氮气（NTF：National Transonic Facility）提高密度来接近Re数到实机。
物理制约用创意克服啊。
机器学习与量纲分析的融合
近年注目的是量纲分析和机器学习的结合。
无量纲学习（Dimensionless Learning）这个方法。机器学习模型的输入输出用无量纲数的话，汎化性能飞跃向上。
Buckingham的π定理作为约束条件组入神经网络，实现与物理法则矛盾不了的预测。Bakarji & Tartakovsky (2021) 的研究提出了从数据自动发现支配的无量纲数的算法。
量纲分析和AI的结合有意思！
多尺度问题中的量纲分析
尺度不同的物理共存问题中，多个无量纲数的相对大小很重要。
比如有微气泡的湍流边界层中：
- 边界层全体：$\text{Re}_\delta \sim 10^4$（惯性支配）
- 气泡近处：$\text{We} \sim O(1)$（表面张力和惯性拮抗）
- 气泡内部：$\text{Re}_b \sim 10^0$（粘性支配）
这样的情况需要各尺度上基于合适的无量纲数的建模，和AMR（Adaptive Mesh Refinement）组合用局部分辨率变化方法。
量纲分析与π定理的故障排除
单位系和标度化引起的故障
量纲分析相关的故障实际有什么样的？
单位系不一致和标度化错误是CFD中最多的初歩错误之一。具体的例子看一下。
1. 压力值完全不对
症状：大气压条件却压力是 $10^9$ Pa 量级
原因：CAD形状是mm单位，物性值用SI（m基准）输入。长度成为1000倍，$\Delta p \propto \rho U^2$ 所以压力成为 $10^6$ 倍。
对策：Fluent的情况，Mesh > Scale把形状换成m单位，或用整合的单位系（密度：tonne/mm³、粘度：MPa·s）。
2. Re数与预期不符
症状：应该层流却有湍流振荡，或反过来应该湍流却是层流解
原因：代表长度或代表速度的定义不适当。比如应该用管径 $D$ 却用半径 $r$（Re偏离2倍）。
对策：
- 手工计算 $\text{Re}$ 确认：$\text{Re} = \rho U D / \mu$
- 管内转变：$\text{Re}_D \approx 2300$
- 外部流（平板）：$\text{Re}_x \approx 5 \times 10^5$
Re数要自己计算验算很重要啊。
3. 阻力系数 $C_D$ 与文献值不符
症状：球的 $C_D$ 与文献值有大偏差
可能的原因：
- Reference Values的面积错了（投影面积 $\pi D^2/4$ vs. 表面积 $\pi D^2$）
- 速度的基准不是自由流速
- 计算域太小，阻挡效应出现（推奨：上游5D以上、下游15D以上、侧方5D以上）
面积的定义是否和文献一致是关键啊。
$y^+$ 相关故障
$y^+$ 是壁面近处网格品质的最重要指标。
| 湍流模型 | 需要的 $y^+$ | 网格战略 |
|---|---|---|
| 标准壁函数 (Standard Wall Function) | 30 - 300 | 壁离开配置第一层 |
| 增强壁处理 (Enhanced Wall Treatment) | $\approx 1$ | 壁贴紧棱晶层 |
| $k$-$\omega$ SST（壁面分辨） | $\approx 1$ | 10层以上棱晶层推奖 |
| SA (Spalart-Allmaras) | $\approx 1$ | 边界层内充足层数 |
$y^+$ 的事前推定式：
平板的经验式 $C_f \approx 0.058\,\text{Re}_L^{-0.2}$ 用的话，分析前能概算第一层厚度。
计算前预测$y^+$设计网格啊。
标度化错误的检测清单
计算结果评价时的"量纲分析的"检查项目。
- 手工计算Re数，确认流动规制（层流/转变/湍流）是否预期通
- $C_D$、$C_L$ 是否数量级适妥当（鈍头物体 $C_D \sim O(1)$、流线形 $O(10^{-2})$）
- 压力损失是否 $\Delta p \sim \frac{1}{2}\rho U^2$ 量级内
- 用非压缩性假定的场合，Ma数0.3以下确认
- 有温度变化，Gr/Re² 大小判断强制/自然对流支配
量纲分析的知识有，结果的妥当性能迅速判断啊。
对。不盲目信计算結果，用無量綱数的数量级推定"奇怪"意識到，CFD工程師的基本素养。
Re数一致也"同样流动"不成立有时候
"模型试验中Re数和实机一致即相似律成立"——多数情况对但陷阱有。比如可压缩性效果（马赫数）、表面张力（韦伯数）、重力（弗劳德数）同时重要问题，全无量纲数同時一致是物理不可能。風洞实验高Re数出的高圧空気使、Re数合数——起。次元分析的素的問題「的無量纲数優先」物理的判断。
错误