🙋
我在学校学过"量纲分析"，但这对实际工作有什么用呢？
🎓
用处大着呢。简单地说，量纲分析可以在有太多变量的复杂现象中理出头绪。比如说圆柱受流力取决于直径D、流速U、流体密度ρ、粘度μ这4个参数。5个变量的组合实验很麻烦吧？但用π定理就能压缩到只需关注2个无量纲数——雷诺数Re和阻力系数C_D——这样实验数据就能整理成一条曲线。
🙋
5个变量压缩到2个！这太神奇了…为什么能这样做呢？
🎓
因为有量纲的约束啊。F、D、U、ρ、μ都可以用质量M、长度L、时间T三个基本量纲来表示。物理方程左右两边的量纲必须相同，这就给了3个约束条件。根据π定理，$n - r = 5 - 3 = 2$，所以就得到2个π群。只要把这2个π群——Re和C_D——组织好，阻力问题就归结为$C_D = f(Re)$这一条曲线。
🙋
我看模拟器左边是原型，右边是模型。这就是"模型相似律"吧？
🎓
对。这是量纲分析最具实用价值的地方。既然$C_D = f(Re)$，只要模型和原型的Re相等，那C_D也就相等，这样从小模型的试验就能预测大原型的力。默认参数下原型D_p=1m、U_p=10m/s，模型D_m=0.10m，如果用同一种流体（比如水），为了Re相等需要把模型流速提高到100m/s——提高了10倍。这就是缩放的严酷现实。
🙋
模型变小了，作用在模型上的力应该会变小吧？但力比卡写的是1.000…
🎓
你眼力不错。同流体、Re相等时，流速提高的能量抵消了直径缩小的效应，力比F_m/F_p = (U_m/U_p)² × (D_m/D_p)²恰好等于1。10倍×10倍×0.1×0.1 = 1。模型体积虽然只有原型的千分之一，受力却一样大。这在传感器选型和夹具设计上很有实际意义。
🙋
那"运动粘度比"改成不是1会怎样？什么时候需要用不同流体？
🎓
原型在空气里、模型在水里试验的情况就是这样。空气的运动粘度大约是水的1/15，所以ν_m/ν_p≈0.07。拉一下这个滑块，你会看到模型流速U_m变小了。通过选用不同流体可以调整模型的试验流速到可行范围。这就是现实中风洞和水槽试验的常见技巧。
1. 设置原型大小和流速（D_p、U_p）：首先设定要预测的原型直径和流速。默认值为D_p = 1m、U_p = 10m/s，对应"水中直径1m的圆柱受到10m/s流速"的情形。
2. 设置模型大小（D_m）：输入模型试验中使用的模型直径。默认值为D_m = 0.10m（原型的1/10）。缩放比λ = D_m/D_p会在计算结果中显示。
3. 选择流体（ν_m/ν_p）：模型试验所用流体的运动粘度比。同流体时为1.0，空气→水时约为0.07等。
4. 查看结果："相似律所需模型流速U_m"是为了使Re与原型一致所需的模型试验流速。"力比F_m/F_p"表示模型受力是原型的多少倍。
步骤1：选定相关变量
影响圆柱周围阻力的变量有：阻力F、直径D、流速U、流体密度ρ、流体粘度μ。暂时忽略表面粗糙度和流动紊乱等因素。
步骤2：用基本量纲表示
F[MLT⁻²]、D[L]、U[LT⁻¹]、ρ[ML⁻³]、μ[ML⁻¹T⁻¹]。基本量纲有M、L、T三种。
步骤3：计算π群数目
根据白金汉π定理：$n - r = 5 - 3 = 2$个独立无量纲群。
步骤4：构成π群
选D、U、ρ为"重复变量"，对F和μ进行无量纲化：
$\pi_1 = \dfrac{\mu}{\rho U D} = \dfrac{1}{Re}$ → 习惯上取倒数得到 $Re = \dfrac{\rho U D}{\mu}$
$\pi_2 = \dfrac{F}{\rho U^2 D^2} \propto C_D$（标准的$C_D$分母为$\tfrac{1}{2}\rho U^2 A$，但形状系数剔除后本质相同）
步骤5：推导关系式
π定理的结论：阻力系数仅为Re的函数：$C_D = f(Re)$。通过实验求得这条曲线后，对任何D、U、ρ、μ的组合都能预测阻力。
步骤6：导出模型相似律
当$Re_m = Re_p$时$C_D$也相等。同流体（$\rho_m = \rho_p$）下有：
$\dfrac{U_m}{U_p} = \dfrac{D_p}{D_m}\cdot\dfrac{\nu_m}{\nu_p}, \qquad \dfrac{F_m}{F_p} = \left(\dfrac{U_m}{U_p}\right)^2 \left(\dfrac{D_m}{D_p}\right)^2$
默认参数的数值计算例：D_p=1.0m、U_p=10m/s、D_m=0.10m、ν_m/ν_p=1.0
缩放比$\lambda = D_m/D_p = 0.100$
$U_m = U_p \cdot (D_p/D_m) \cdot (\nu_m/\nu_p) = 10 \times 10 \times 1 = 100$ m/s
假定20°C水的ν_p = 1.0×10⁻⁶ m²/s：$Re = U_p D_p / \nu_p = 10 \times 1 / 10^{-6} = 1.0\times10^7$
模型侧：$Re_m = U_m D_m / \nu_m = 100 \times 0.1 / 10^{-6} = 1.0\times10^7$ 完全相等 ✓
力比$F_m/F_p = (100/10)^2 \times (0.1/1)^2 = 100 \times 0.01 = 1.000$
风洞试验中的航空器·汽车设计：直接用大飞机或汽车做试验成本太高、周期太长。根据白金汉π定理的相似律，可以用1/10～1/50的模型在风洞中试验，通过匹配Re数和Ma数来预测实机的升力、阻力和力矩。F1赛队利用风洞模型节省时间和成本的技术正是建立在这个原理上。
船舶模型试验（拖曳水槽）：船舶的造波阻力主要由弗劳德数Fr控制。通过调整模型试验的Fr使其与实船相同，可以用可操作尺寸的模型船评估巨大油轮或航母的阻力。粘性阻力（Re控制的部分）通过"ITTC1957摩擦线"等经验公式单独补正，这叫"Fr分割法"，是行业标准。
化工装置的放大：实验室烧杯规模开发出的反应工艺、混合过程要放大到工业规模（体积增加数百到数千倍）时，相似律至关重要。搅拌Re数、Nusselt数、阻滞数等各种π群需要优先级排序，选择是"几何相似"还是"动态相似"还是"热相似"，这决定了放大的可行性和成本。
河流·港湾水工模型：洪水流动、海啸爬升、港湾波浪——这些现象既用数值仿真检验，也用缩尺物理模型试验。根据Fr相似律缩放地形，调整流速和时间，就能在模型水槽中重现真实水工现象，验证堤防越流或防浪堤的效果。
最常见的误解是，"模型做小了，试验就容易了"这样的想法。实际上恰恰相反，Re相等条件下需要提高模型流速，$U_m = U_p \cdot (D_p/D_m)$，所以1/10模型需要10倍流速，1/100模型需要100倍流速。原型在水中10m/s，1/100模型就需要水中1000m/s——根本无法实现。这就是为什么要换流体（降低运动粘度）、提高压力（高压风洞）或优先满足某一种π群等实际技巧。
第二常见的误解是，"能把所有π群都同时满足"这样的期望。船舶自由液面流动中Re和Fr都很关键，但用同一流体同时满足两者需要ν_m/ν_p = (D_m/D_p)^(3/2)的流体——现实中不存在。航空器跨音速试验需要同时控制Re和Ma，也只能"挑一个优先，其他补正"。白金汉π定理教会的是"能把变量简化"，不是"能变魔法一步到位"。
最后一个常见错误是，"量纲分析就能决定物理"的假设。π定理确保了$C_D = f(Re)$这样的函数形式，但函数$f$的内容——C_D随Re如何变化、为什么在Re≈10⁵处有"阻力危机"突变、为什么超临界区近似稳定——这些必须通过实验或CFD来发现。模拟器显示的是"相似条件"和"力的缩放"，不能直接输出圆柱阻力的数值。π定理是整理工具，不能替代物理规律本身。
Re一致要求$U_m = U_p (D_p/D_m)$，Fr一致要求$U_m = U_p\sqrt{D_m/D_p}$。同时满足两者需$(D_p/D_m) = \sqrt{D_m/D_p} \cdot (\nu_p/\nu_m)$，即$\nu_m = \nu_p (D_m/D_p)^{3/2}$。1/100模型需要水的运动粘度的1/1000——这样的流体不存在。船舶模型试验因此采用"Fr优先"策略，Re影响的粘性阻力通过ITTC等经验公式另行补正，称为"Fr分割法"。
标准阻力系数定义为$C_D = F / (\tfrac{1}{2}\rho U^2 A)$，其中$A$是参考面积（圆柱为$D \times \ell$）。本工具为简化重点，只展示π群的本质组合$\pi_2 \propto F / (\rho U^2 D^2)$，省去1/2等系数和形状因子。量纲分析阶段不确定这些系数，后期可调整为$\pi_2 = (1/2) \times C_D \times $(形状因子)。动态相似和力比讨论用简化形式足够，实务中应用标准定义。
π群的个数(n-r)是固定的，但具体选法有自由度。比如可以用π_1'= Re²或π_2'= π_2/π_1，只要相互独立就行。习惯上选物理意义明确、工程中普遍认知的形式（Re、Fr、Ma、Nu等）。选择"重复变量"的方法也会影响最终π群的形式，初学者最安全的做法是按教科书步骤，最后对照已知的无量纲数验证。
本工具假定"用同一流体做模型试验"这样的标准情形，故$\rho_m/\rho_p = 1$固定。不同流体（如原型空气、模型水）的情况在运动粘度比ν_m/ν_p的滑块中表现，因为ν = μ/ρ。Re一致条件$\rho_m U_m D_m / \mu_m = \rho_p U_p D_p / \mu_p$可重写为$U_m D_m / \nu_m = U_p D_p / \nu_p$，密度差异被运动粘度吸收。力比的通用形式是$F_m/F_p = (\rho_m/\rho_p)(U_m/U_p)^2(D_m/D_p)^2$，密度比≠1时需另外乘以。