量纲分析：修订间差异

2026年9月20日 (日) 07:18的最新版本

量纲分析（dimensional analysis）通过长度、时间、质量等基本量的幂次关系，检查物理公式，找出不受单位选择影响的组合，并比较不同尺度的系统。它常能在完整求解之前确定某种比例关系。

摆长增加四倍，周期怎样变化

一个理想单摆的摆长为 $ℓ$ ，重力加速度为 $g$ ，摆球质量为 $m$ 。在相同振幅角下，把摆长从一米改成四米，周期是否也变为四倍？

先不求摆的运动方程，而是寻找能组成时间的量。长度的量纲记为 $L$ ，时间记为 $T$ ，质量记为 $M$ 。摆长的量纲是 $L$ ，重力加速度是 $L T^{- 2}$ ，质量是 $M$ 。假设周期 $P$ 可按这些量写成幂乘积 $P = C ℓ^{a} g^{b} m^{c} .$ $C$ 是无量纲系数， $a, b, c$ 是待求指数。左右量纲应相同，所以 $T = L^{a + b} T^{- 2 b} M^{c} .$ 逐项比较指数：左边没有质量因子，故 $c = 0$ ；时间指数给出 $- 2 b = 1$ ，故 $b = - 1 / 2$ ；长度指数为零，故 $a + b = 0$ ，得到 $a = 1 / 2$ 。于是 $P = C \sqrt{\frac{ℓ}{g}} .$ 在同一重力、相同振幅及这些模型假设下，摆长增加四倍，周期便增加 $\sqrt{4} = 2$ 倍。质量没有形成独立的时间尺度。

系数 $C$ 仍未确定。振幅角 $θ_{0}$ 没有量纲，但它可以变化，所以更一般的关系是 $P = \sqrt{ℓ / g} F (θ_{0})$ 。理想小振幅近似进一步给出方程 $θ^{″} + (g / ℓ) θ = 0$ ，其角频率为 $\sqrt{g / ℓ}$ ，这才得到 $P \approx 2 π \sqrt{ℓ / g}$ 。平方根尺度来自量纲， $2 π$ 和小振幅适用范围来自运动方程。

一米与一百厘米，哪些东西变了

摆长是一项物理量；长度 $L$ 是它的量纲；米和厘米是可选择的单位。换成厘米时，一米记为一百厘米， $g$ 的数值也由约 $9.8$ 米每二次方秒变为约 $980$ 厘米每二次方秒。比值 $ℓ / g$ 没变，所以算出的周期仍相同。

速度的量纲为 $L T^{- 1}$ ，加速度为 $L T^{- 2}$ ，力为 $M L T^{- 2}$ 。相乘时量纲指数相加，相除时相减；可以相加的物理项则必须量纲相同。例如位移公式 $s = v_{0} t + \frac{1}{2} a t^{2}$ 右边两项都是长度。如果漏写一个 $t$ ，成为 $a t / 2$ ，这一项是速度，便不能与位移相加。

但把正确公式中的系数改成十七和九，量纲仍然一致。量纲检查可排除不一致式子，不能单独确定数值系数。类似地，功和力矩虽然量纲相同，物理含义仍不同。

指数、对数与三角函数的输入也应采用无量纲组合。例如 $e^{- t / τ}$ 中 $τ$ 是时间尺度，比值与秒或分钟的选择无关；直接写“负三秒”的指数就缺少这种单位不变性。摄氏温度与开尔文读数还涉及零点平移，热力学公式应使用其要求的绝对温度，不能只改单位名称。

无量纲组合怎样由矩阵找到

将单摆变量 $P, ℓ, g, m$ 的量纲指数按列排好，按质量、长度、时间排列行，得到 $D = (\begin{matrix} 0 & 0 & 0 & 1 \\ 0 & 1 & 1 & 0 \\ 1 & 0 & - 2 & 0 \end{matrix}) .$ 例如第三列 $(0, 1, - 2)^{𝖳}$ 就表示 $g$ 的量纲。若希望乘积 $Π = P^{z_{1}} ℓ^{z_{2}} g^{z_{3}} m^{z_{4}}$ 无量纲，三种总指数都应为零，即 $D z = 0.$ 逐行求解： $z_{4} = 0$ ， $z_{2} + z_{3} = 0$ ， $z_{1} - 2 z_{3} = 0$ 。令 $z_{1} = 1$ ，得到 $z = (1, - 1 / 2, 1 / 2, 0)$ ，于是 $Π = P \sqrt{g / ℓ} .$ 它没有单位，且无论用米还是厘米算，数值相同。

矩阵有四列、秩为三，零空间是一维。加上无量纲振幅 $θ_{0}$ 后，多了一列全零，便有两组独立无量纲量： $P \sqrt{g / ℓ}$ 与 $θ_{0}$ 。前一节的未知函数 $F$ 正是这两个量之间的关系。

Buckingham Π 定理把这一过程推广到 $n$ 个变量：在量纲齐次、单位缩放适用及相应变量范围条件下，量纲矩阵秩为 $r$ ，就有 $n - r$ 组独立的无量纲幂乘积，可用来表达关系。对正变量取对数后，单位变化成为线性平移；零空间中的向量与这些单位变化方向正交，给出保持不变的组合。这是“变量数减秩”的线性代数来源。

组合不唯一，例如对正值的 $Π$ 使用 $Π^{2}$ 也可表达同一尺度信息。若能把输出单值地解出，才进一步写成一个无量纲量等于其他无量纲量的函数；定理本身并不决定这个函数。

从流动变量中找出雷诺数

考虑流体密度 $ρ$ 、速度尺度 $U$ 、长度尺度 $L_{0}$ 、动力黏度 $μ$ 。量纲分别是 $M L^{- 3}$ 、 $L T^{- 1}$ 、 $L$ 、 $M L^{- 1} T^{- 1}$ 。寻找无量纲乘积 $ρ^{a} U^{b} L_{0}^{c} μ^{d}$ ，逐项列出指数条件： $a + d = 0, - 3 a + b + c - d = 0, - b - d = 0.$ 先选 $a = 1$ ，第一式给出 $d = - 1$ ，第三式给出 $b = 1$ ，第二式再给出 $c = 1$ ，得到雷诺数 $Re = \frac{ρ U L_{0}}{μ} = \frac{U L_{0}}{ν}, ν = \frac{μ}{ρ} .$ $ν$ 称为运动黏度，量纲为 $L^{2} T^{- 1}$ ，与动力黏度不同。

取教学数值 $ρ = 1000 k g / m^{3}$ 、 $U = 0.2 m / s$ 、 $L_{0} = 0.01 m$ 、 $μ = 0.001 P a s$ ，代入得 $Re = 2000$ 。在流动方程中，惯性项的尺度为 $ρ U^{2} / L_{0}$ ，黏性项的尺度为 $μ U / L_{0}^{2}$ ；两者相除正是雷诺数。因此它比较的是惯性与黏性效应的相对大小。

具体流态还与几何和边界条件有关，不能只由二千这个数判断任意流动是否转捩。雷诺数的物理用途来自它在方程中的位置。

预测阻力时，输出也要列入变量

若目标是固定形状物体的阻力 $F_{D}$ ，还要把输出力列入清单，成为 $F_{D}, ρ, U, L_{0}, μ$ 五个量。对不可压缩、无相关自由表面作用的模型，量纲秩为三，可取两组独立组合： $Π_{1} = \frac{F_{D}}{ρ U^{2} L_{0}^{2}}, Π_{2} = Re .$ 分母 $ρ U^{2} L_{0}^{2}$ 的量纲为 $M L T^{- 2}$ ，恰好是力。于是可在单值范围写成 $F_{D} / (ρ U^{2} L_{0}^{2}) = Φ (Re)$ ，通过实验或流动计算确定 $Φ$ 。MIT：量纲分析与模型相似

重复变量法可以选 $ρ, U, L_{0}$ ，因为它们的三列量纲独立；分别与阻力、黏度组合即可消去量纲。若选择量纲相关的变量，便缺少消去某个基本量纲的自由度。

这里默认尺度 $ρ, U, L_{0}$ 为正，速度为零时应另外处理。若表面粗糙度 $k$ 也影响阻力，就需加入 $k / L_{0}$ ；改变物体形状时，还要保留长宽比等几何量。遗漏一个重要变量不会被矩阵运算自动发现，变量清单首先来自物理问题。

为什么缩小的模型未必能同时保持所有相似条件

设模型的长度为原型的 $λ$ 倍，速度比记为 $q = U_{m} / U_{p}$ 。同一种流体意味着运动黏度相同，于是 $\frac{{Re}_{m}}{{Re}_{p}} = q λ .$ 保持雷诺数，必须取 $q = 1 / λ$ 。长度缩为十分之一，就需把速度增为十倍。

若流动还有重力自由表面效应，需要比较弗劳德数 $Fr = U / \sqrt{g L_{0}}$ 。在同一重力下，它的比值为 $q / \sqrt{λ}$ ；保持它不变则要求 $q = \sqrt{λ}$ 。长度缩为十分之一时，速度反而应降到约 $0.316$ 倍。

长度缩为十分之一的模型在同一流体同一重力下，匹配雷诺数要求速度十倍，匹配弗劳德数要求速度约零点三一六倍 — 同一个长度缩尺，对两种无量纲数提出不同的速度要求。

图中两条箭头从同一缩尺条件出发。上方保持惯性与黏性比，下方保持速度与重力波动尺度的比，它们不能由同一速度同时实现。若优先匹配弗劳德数，雷诺数比就成为 $λ^{3 / 2} \approx 0.0316$ ；若还要匹配雷诺数，运动黏度也必须降到原型的约 $0.0316$ 倍。MIT：Marine Hydrodynamics，Lecture 19

量纲分析把这种冲突明确列了出来。实验中需根据主要机制决定匹配哪些量，并分析未匹配量的影响；只有当几何、边界条件和相关无量纲参数相应一致时，才可把模型结果换算到原型。

历史

Buckingham 于 1914 年发表关于物理相似系统的工作，系统表达了后来称为 Π 定理的无量纲化方法。NIST：The Life of (Buckingham) Pi 这一框架把已有的尺度与相似思想联系到单位变换，使多变量实验问题能够改写为更少的无量纲关系。

参考来源与延伸阅读

NIST：The Life of (Buckingham) Pi，历史与单位缩放说明。
MIT：Dimensional Analysis of Models and Data Sets，模型相似与尺度分析。
MIT：Marine Hydrodynamics，Lecture 19：雷诺数与弗劳德数同时匹配的限制。
先修：矩阵、向量空间；相关：数学建模、微分方程、一阶线性微分方程。

@@ 第1行： / 第1行： @@
-量纲分析研究物理量在基本量纲和单位变换下的结构，用于检查方程一致性、构造无量纲组合，以及分析不同尺度系统之间的相似关系。它能约束模型形式，却通常不能独自确定无量纲常数、具体函数或物理机制。
+'''量纲分析'''（dimensional analysis）通过长度、时间、质量等基本量的幂次关系，检查物理公式，找出不受单位选择影响的组合，并比较不同尺度的系统。它常能在完整求解之前确定某种比例关系。
-英文名称：Dimensional analysis。
+== 摆长增加四倍，周期怎样变化 ==
+一个理想单摆的摆长为 <math>\ell</math>，重力加速度为 <math>g</math>，摆球质量为 <math>m</math>。在相同振幅角下，把摆长从一米改成四米，周期是否也变为四倍？
-== English overview ==
+先不求摆的运动方程，而是寻找能组成时间的量。长度的量纲记为 <math>L</math>，时间记为 <math>T</math>，质量记为 <math>M</math>。摆长的量纲是 <math>L</math>，重力加速度是 <math>LT^{-2}</math>，质量是 <math>M</math>。假设周期 <math>P</math> 可按这些量写成幂乘积
-<div lang="en" class="math-english-summary">
+<math display="block">P=C\ell^ag^bm^c.</math>
-Dimensional analysis examines how quantities transform when fundamental units are changed. Dimensions describe types of quantities, such as length, mass, and time; units provide particular measuring scales. A valid physical equation must be dimensionally homogeneous, but homogeneity alone cannot establish the equation's physical correctness. Dimensionless numerical factors and functions still require theory, experiments, or additional information.
+<math>C</math> 是无量纲系数，<math>a,b,c</math> 是待求指数。左右量纲应相同，所以
+<math display="block">T=L^{a+b}T^{-2b}M^c.</math>
+逐项比较指数：左边没有质量因子，故 <math>c=0</math>；时间指数给出 <math>-2b=1</math>，故 <math>b=-1/2</math>；长度指数为零，故 <math>a+b=0</math>，得到 <math>a=1/2</math>。于是
+<math display="block">P=C\sqrt{\frac\ell g}.</math>
+在同一重力、相同振幅及这些模型假设下，摆长增加四倍，周期便增加 <math>\sqrt4=2</math> 倍。质量没有形成独立的时间尺度。
-This article begins with unit conversions and homogeneous equations, then derives scaling relationships for a pendulum and a fluid-flow example. The Buckingham Pi framework organizes dimensionless products using the null space of a dimension matrix. The number of independent products depends on the matrix rank, and different valid choices of products describe the same dimensional structure. We show why a mass variable drops out of an idealized pendulum scaling law only after the relevant-variable assumptions have been made, and why the factor of two pi requires a dynamical derivation. Reynolds number illustrates the ratio of competing physical effects rather than a universal threshold for every flow. Historical notes connect Buckingham's 1914 work with the broader practice of similarity analysis. The final discussion separates dimensional consistency, dynamic similarity, and empirical validation, with explicit attention to omitted variables and temperature-scale offsets.
+系数 <math>C</math> 仍未确定。振幅角 <math>\theta_0</math> 没有量纲，但它可以变化，所以更一般的关系是 <math>P=\sqrt{\ell/g}\,F(\theta_0)</math>。理想小振幅近似进一步给出方程 <math>\theta''+(g/\ell)\theta=0</math>，其角频率为 <math>\sqrt{g/\ell}</math>，这才得到 <math>P\approx2\pi\sqrt{\ell/g}</math>。平方根尺度来自量纲，<math>2\pi</math> 和小振幅适用范围来自运动方程。
-</div>
-== 量、量纲和单位是三件不同的事 ==
+== 一米与一百厘米，哪些东西变了 ==
-一段杆的长度是物理量，长度量纲通常记 L，米和厘米是该量可用的不同单位。把一米写成一百厘米，物理对象和量纲都不变，只是数值随单位改变。时间用 T、质量用 M 作为基本量纲符号时，速度量纲为 <math>LT^{-1}</math>，加速度为 <math>LT^{-2}</math>，力为 <math>MLT^{-2}</math>。
+摆长是一项物理量；长度 <math>L</math> 是它的量纲；米和厘米是可选择的单位。换成厘米时，一米记为一百厘米，<math>g</math> 的数值也由约 <math>9.8</math> 米每二次方秒变为约 <math>980</math> 厘米每二次方秒。比值 <math>\ell/g</math> 没变，所以算出的周期仍相同。
-量纲相同不代表物理意义相同。功与力矩都具有力乘长度的量纲，角度在常用 SI 量纲处理中是无量纲，能量与力矩仍是不同的物理概念。量纲分析是一种结构约束，不能把所有同量纲的量任意互换。反过来，数字看起来相同的两个量也未必可比较，例如三米与三秒。
+速度的量纲为 <math>LT^{-1}</math>，加速度为 <math>LT^{-2}</math>，力为 <math>MLT^{-2}</math>。相乘时量纲指数相加，相除时相减；可以相加的物理项则必须量纲相同。例如位移公式
+<math display="block">s=v_0t+\frac12at^2</math>
+右边两项都是长度。如果漏写一个 <math>t</math>，成为 <math>at/2</math>，这一项是速度，便不能与位移相加。
-基本量纲的选择取决于理论范围。本条示例只使用质量、长度、时间，涉及电磁、温度等问题时还需相应量纲或明确的单位体系。不能因为本例矩阵只有三行，就断言所有物理问题永远只有三种基本量。
+但把正确公式中的系数改成十七和九，量纲仍然一致。量纲检查可排除不一致式子，不能单独确定数值系数。类似地，功和力矩虽然量纲相同，物理含义仍不同。
-== 方程一致性检查 ==
+指数、对数与三角函数的输入也应采用无量纲组合。例如 <math>e^{-t/\tau}</math> 中 <math>\tau</math> 是时间尺度，比值与秒或分钟的选择无关；直接写“负三秒”的指数就缺少这种单位不变性。摄氏温度与开尔文读数还涉及零点平移，热力学公式应使用其要求的绝对温度，不能只改单位名称。
-加减号两侧参与相加的项必须具有相同量纲。例如 <math>s=v_0t+at^2/2</math> 中，每项量纲均为 L；若误写成 <math>s=v_0t+at/2</math>，最后一项量纲变成速度，立即发现不一致。乘除则对应量纲指数相加、相减，所以密度为 <math>ML^{-3}</math>，压力为 <math>ML^{-1}T^{-2}</math>。
-指数函数、对数函数的输入在物理公式中应是无量纲组合。例如衰减表达式 <math>N(t)=N_0e^{-t/\tau}</math> 合理，因为 t/τ 无量纲；若直接写 e 的“负三秒”次方，换用分钟后数值形式会变得不一致。三角函数中的相位通常以弧度表示，也应明确它来自一个无量纲角度或比值。
+== 无量纲组合怎样由矩阵找到 ==
+将单摆变量 <math>P,\ell,g,m</math> 的量纲指数按列排好，按质量、长度、时间排列行，得到
-量纲一致是必要检查，不是充分证明。式子 <math>s=17v_0t+9at^2</math> 仍然量纲一致，但不等于标准匀加速运动关系。分析不能从量纲本身判断系数应为一和二分之一，系数与函数形式还需要动力学推导、边界条件或实验资料。
-== 完整算例一：单摆周期的尺度 ==
-考虑理想单摆，在特定假设下，周期 P 可能依赖摆长 ℓ、重力加速度 g、质量 m 和振幅角 θ₀。先暂把角度固定，假设其余变量以幂乘积形式进入：<math>P=C\ell^ag^bm^c</math>。量纲给出
-<math display="block">T=L^a(LT^{-2})^bM^c=L^{a+b}T^{-2b}M^c.</math>
-比较指数：a+b=0，−2b=1，c=0，因此 a=1/2、b=−1/2，得到 <math>P=C\sqrt{\ell/g}</math>。这说明在给定变量假设下周期的尺度与摆长平方根成正比，质量不提供独立尺度。
-但 θ₀ 无量纲，量纲分析无法排除 C 依赖振幅，因此更一般地应写 <math>P=\sqrt{\ell/g}\,F(\theta_0)</math>。只有进一步使用理想小角度方程 <math>\theta''+(g/\ell)\theta=0</math>，解的角频率为 <math>\sqrt{g/\ell}</math>，才得到小振幅极限下 F≈2π。不能说量纲分析单独推导出了 2π，也不能把近似结果扩展到任意振幅。
-例如在同一 g 和相同小振幅条件下，摆长从一米变为四米，周期约变为原来的两倍；质量从一千克变为两千克，在该理想模型中不改变周期。若空气阻力、绳子弹性、摆球尺寸或其他因素不可忽略，就引入了新的变量和无量纲组合，原先的变量清单不再完整。
-== 从指数方程走向量纲矩阵 ==
-把每个变量在基本量纲中的指数写成一列，就得到量纲矩阵 D。对于 P、ℓ、g、m，按 M、L、T 排列行，有
 <math display="block">D=\begin{pmatrix}0&0&0&1\\0&1&1&0\\1&0&-2&0\end{pmatrix}.</math>
-一个乘积 <math>\Pi=P^{a_1}\ell^{a_2}g^{a_3}m^{a_4}</math> 无量纲，当且仅当 <math>Da=0</math>。这里求的是[[向量空间]]中的零空间，不是凭眼睛猜单位怎么抵消。该矩阵秩为三，四列减三得到一维零空间，可以选 <math>\Pi=P\sqrt{g/\ell}</math>。
+例如第三列 <math>(0,1,-2)^{\mathsf T}</math> 就表示 <math>g</math> 的量纲。若希望乘积 <math>\Pi=P^{z_1}\ell^{z_2}g^{z_3}m^{z_4}</math> 无量纲，三种总指数都应为零，即
+<math display="block">Dz=0.</math>
-把无量纲振幅 θ₀ 也列为变量，相当于再添一列全零，零空间维数变成二。两组独立组合可取 <math>\Pi_1=P\sqrt{g/\ell}</math> 和 <math>\Pi_2=\theta_0</math>，于是量纲一致的物理关系可以写成两个无量纲量之间的关系。在适用的单值情形下，才进一步写成 <math>\Pi_1=F(\Pi_2)</math>。
+逐行求解：<math>z_4=0</math>，<math>z_2+z_3=0</math>，<math>z_1-2z_3=0</math>。令 <math>z_1=1</math>，得到 <math>z=(1,-1/2,1/2,0)</math>，于是
+<math display="block">\Pi=P\sqrt{g/\ell}.</math>
-选出的组合并不唯一：<math>\Pi_1^2</math> 也无量纲，若在正值范围内使用，能表达相同的尺度关系。选择时通常追求容易解释和测量，而不是寻找唯一“正确名字”。需要区分独立组合的数量与某一组具体基的选择。
+它没有单位，且无论用米还是厘米算，数值相同。
-== Buckingham Π 定理讲了什么 ==
-在相应的单位变换和量纲齐次假设下，若一个关系涉及 n 个变量，量纲矩阵的秩为 r，就可以用 n−r 个独立无量纲组合来表达其量纲结构。这是 Buckingham Π 框架的核心。严格表述还需说明变量范围、正值幂乘积的可用性和关系的可表示条件；不能把一句口号当作任何任意函数都适用的无条件结论。
-其代数依据可从单位缩放理解。改变各基本单位的尺度时，每个变量数值按对应量纲指数缩放；乘积的总指数为零，就抵消所有这些缩放。取对数后，乘法缩放变为线性平移，而零空间向量给出不随这些平移改变的组合。秩—零化度定理解释了为何数目是 n−r，而不是永远用变量数减去三。
-定理的实用价值是压缩需要实验探索的维数。例如一个现象表面依赖许多有量纲参数，转换后可能只需研究两个或三个比值之间的函数关系。但未知函数仍然需要确定；量纲方法减少问题，不自动完成全部求解。若漏列表面粗糙度或温度，得到的压缩关系也可能无法覆盖实际现象。
-== 完整算例二：雷诺数从何而来 ==
-考虑流体密度 ρ、速度尺度 U、长度尺度 L₀、动力黏度 μ。它们量纲依次为 <math>ML^{-3}</math>、<math>LT^{-1}</math>、L 和 <math>ML^{-1}T^{-1}</math>。寻找乘积 <math>\rho^aU^bL_0^c\mu^d</math> 的无量纲条件，得到 a+d=0、−3a+b+c−d=0、−b−d=0。选 a=1，便有 d=−1、b=1、c=1，所以
-<math display="block">\operatorname{Re}=\frac{\rho UL_0}{\mu}.</math>
-同一数也可写 U L₀/ν，其中运动黏度 ν=μ/ρ。必须区分 μ 和 ν 的量纲，不能因为中文都含“黏度”就混用。
-从流动方程的量级估计，惯性项尺度为 <math>\rho U^2/L_0</math>，黏性项为 <math>\mu U/L_0^2</math>，两者比值就是 Re，因此它反映这两种效应的相对尺度。取教学数值 ρ=1000 千克每立方米、U=0.2 米每秒、L₀=0.01 米、μ=0.001 帕秒，得到 Re=2000。这只是计算演示，不凭这个数字就断言任意流动都处在相同状态。
-流态还与几何、入口扰动、边界条件和其他无量纲参数有关。某类管流的经验阈值不能直接移到绕球流、明渠流或微尺度流动。无量纲数有助于组织比较，但它不是跨所有场景通用的标签。
-== 补上输出变量：阻力关系的两组 Π ==
+矩阵有四列、秩为三，零空间是一维。加上无量纲振幅 <math>\theta_0</math> 后，多了一列全零，便有两组独立无量纲量：<math>P\sqrt{g/\ell}</math> 与 <math>\theta_0</math>。前一节的未知函数 <math>F</math> 正是这两个量之间的关系。
-前一例只找到了四个流动参数中的无量纲数，还没有说明希望预测什么。若要预测固定形状物体的阻力 <math>F_D</math>，并假设流体不可压缩、没有相关自由表面效应，变量清单应包含 <math>F_D,\rho,U,L_0,\mu</math>。现在是五个变量、量纲秩三，需两组独立组合：
-<math display="block">\Pi_1=\frac{F_D}{\rho U^2L_0^2},\qquad\Pi_2=\frac{\rho UL_0}{\mu}.</math>
-检查第一组的分母量纲：密度乘速度平方再乘长度平方，恰为 <math>MLT^{-2}</math>，与力相同。在可单值表达的范围，量纲结构允许写 <math>\Pi_1=\Phi(\Pi_2)</math>；它并没有提供未知函数 <math>\Phi</math> 的具体形状。
-采用重复变量法时，可以选密度、速度和长度，因为三列量纲独立；再分别用它们配合阻力与动力黏度消去量纲。若误选一组量纲相关的重复变量，指数方程可能不能给出所需独立组合。输出量必须进入总变量计数，不能只数输入，再照抄变量数减秩的公式。[https://ocw.mit.edu/courses/res-12-001-topics-in-fluid-dynamics-fall-2024/pages/essay-2-dimensional-analysis-of-models-and-data-sets-similarity-solutions-and-scaling-analysis/ MIT：变量清单与零空间方法]
+Buckingham Π 定理把这一过程推广到 <math>n</math> 个变量：在量纲齐次、单位缩放适用及相应变量范围条件下，量纲矩阵秩为 <math>r</math>，就有 <math>n-r</math> 组独立的无量纲幂乘积，可用来表达关系。对正变量取对数后，单位变化成为线性平移；零空间中的向量与这些单位变化方向正交，给出保持不变的组合。这是“变量数减秩”的线性代数来源。
-若物体表面粗糙度的长度尺度 <math>k</math> 会影响结果，就要增加 <math>k/L_0</math>；若形状还由长宽比决定，也要保留相应几何比值。“固定形状”意味着这些比值固定，不意味着它们在改变形状后仍可省略。量纲方法不能仅凭数据表告诉哪些机制重要，这一选择仍需要物理知识、实验或模型比较。
+组合不唯一，例如对正值的 <math>\Pi</math> 使用 <math>\Pi^2</math> 也可表达同一尺度信息。若能把输出单值地解出，才进一步写成一个无量纲量等于其他无量纲量的函数；定理本身并不决定这个函数。
-上述阻力无量纲化假定所选密度、长度与速度尺度为正。若直接取静止速度零，分母也为零，就不能沿用同一组合来表示该点；可研究趋近零的极限，或选择另一套非零特征尺度。一个公式在哪些参数域有效，与量纲是否抵消是两个独立问题，换标度时应明确说明。
+== 从流动变量中找出雷诺数 ==
+考虑流体密度 <math>\rho</math>、速度尺度 <math>U</math>、长度尺度 <math>L_0</math>、动力黏度 <math>\mu</math>。量纲分别是 <math>ML^{-3}</math>、<math>LT^{-1}</math>、<math>L</math>、<math>ML^{-1}T^{-1}</math>。寻找无量纲乘积 <math>\rho^aU^bL_0^c\mu^d</math>，逐项列出指数条件：
+<math display="block">a+d=0,\qquad -3a+b+c-d=0,\qquad-b-d=0.</math>
+先选 <math>a=1</math>，第一式给出 <math>d=-1</math>，第三式给出 <math>b=1</math>，第二式再给出 <math>c=1</math>，得到雷诺数
+<math display="block">\operatorname{Re}=\frac{\rho UL_0}{\mu}=\frac{UL_0}{\nu},\qquad\nu=\frac\mu\rho.</math>
+<math>\nu</math> 称为运动黏度，量纲为 <math>L^2T^{-1}</math>，与动力黏度不同。
-== 模型相似与缩尺实验 ==
+取教学数值 <math>\rho=1000\ \mathrm{kg/m^3}</math>、<math>U=0.2\ \mathrm{m/s}</math>、<math>L_0=0.01\ \mathrm m</math>、<math>\mu=0.001\ \mathrm{Pa\,s}</math>，代入得 <math>\operatorname{Re}=2000</math>。在流动方程中，惯性项的尺度为 <math>\rho U^2/L_0</math>，黏性项的尺度为 <math>\mu U/L_0^2</math>；两者相除正是雷诺数。因此它比较的是惯性与黏性效应的相对大小。
-若小模型和原型在几何、边界条件以及相关无量纲参数上满足相应相似要求，模型结果才有希望按比例转换到原型。例如同种流体下把长度缩成十分之一，若要保持雷诺数相同，就需把速度增为十倍，因为 Re 与 U L₀ 成正比。这个条件可能在设备上难以实现，或者同时破坏其他相似条件。
-若还涉及重力自由表面效应，通常还要考虑另一类无量纲比值；多个条件未必能在同一种流体与简单尺度变换下同时满足。工程实验需要明确优先保留哪些机制以及牺牲哪些条件，不能只说“做了一个缩小版，所以结果相同”。量纲分析在这里帮助指出相似条件之间的冲突，而不替代实验判断。
+具体流态还与几何和边界条件有关，不能只由二千这个数判断任意流动是否转捩。雷诺数的物理用途来自它在方程中的位置。
-=== 同时匹配雷诺数与弗劳德数的冲突 ===
+== 预测阻力时，输出也要列入变量 ==
-若自由表面的重力效应也重要，可取弗劳德数 <math>\operatorname{Fr}=U/\sqrt{gL_0}</math>。设模型长度为原型的 <math>\lambda</math> 倍，同一重力和同一种流体下，保持雷诺数要求速度比为 <math>1/\lambda</math>，保持弗劳德数却要求速度比为 <math>\sqrt\lambda</math>。除非比例为一，两种要求无法同时成立。
+若目标是固定形状物体的阻力 <math>F_D</math>，还要把输出力列入清单，成为 <math>F_D,\rho,U,L_0,\mu</math> 五个量。对不可压缩、无相关自由表面作用的模型，量纲秩为三，可取两组独立组合：
+<math display="block">\Pi_1=\frac{F_D}{\rho U^2L_0^2},\qquad\Pi_2=\operatorname{Re}.</math>
+分母 <math>\rho U^2L_0^2</math> 的量纲为 <math>MLT^{-2}</math>，恰好是力。于是可在单值范围写成 <math>F_D/(\rho U^2L_0^2)=\Phi(\operatorname{Re})</math>，通过实验或流动计算确定 <math>\Phi</math>。[https://ocw.mit.edu/courses/res-12-001-topics-in-fluid-dynamics-fall-2024/pages/essay-2-dimensional-analysis-of-models-and-data-sets-similarity-solutions-and-scaling-analysis/ MIT：量纲分析与模型相似]
-例如长度缩为十分之一，前者要求速度增至十倍，后者要求速度降到约零点三一六倍。若优先保持弗劳德数，则雷诺数比为 <math>\lambda^{3/2}\approx0.0316</math>。为了同时匹配两者且保持同一重力，模型流体的运动黏度也必须改成原型的约零点零三一六倍；这未必在设备与流体选择上可行。[https://ocw.mit.edu/courses/2-20-marine-hydrodynamics-13-021-spring-2005/b014d38899d92006e897830ac4c6e3e2_lecture19.pdf MIT：Marine Hydrodynamics，Lecture 19]
+重复变量法可以选 <math>\rho,U,L_0</math>，因为它们的三列量纲独立；分别与阻力、黏度组合即可消去量纲。若选择量纲相关的变量，便缺少消去某个基本量纲的自由度。
-这个矛盾不是量纲分析失败，而是它揭示缩尺实验不能自动保留全部动力学。实际试验可能分离估计不同阻力贡献或采用经验修正，但修正需要另外的理论与校准。即使雷诺数相同，入口条件、相对粗糙度或可压缩效应不同，也不能据一个数宣称完全动力相似；哪些无量纲参数可以忽略，应有适用范围的依据。
+这里默认尺度 <math>\rho,U,L_0</math> 为正，速度为零时应另外处理。若表面粗糙度 <math>k</math> 也影响阻力，就需加入 <math>k/L_0</math>；改变物体形状时，还要保留长宽比等几何量。遗漏一个重要变量不会被矩阵运算自动发现，变量清单首先来自物理问题。
-== 单位换算与常见陷阱 ==
+== 为什么缩小的模型未必能同时保持所有相似条件 ==
-长度改用厘米而时间仍用秒时，速度数值也必须相应换算，所有进入 Re 的量应使用自洽单位。一个无量纲结果如果随使用米还是厘米而变化，说明换算或公式有错误。百分数、角度制与弧度制虽然涉及无量纲量，数值约定仍不同，不能忽略转换因子。
+设模型的长度为原型的 <math>\lambda</math> 倍，速度比记为 <math>q=U_{\mathrm m}/U_{\mathrm p}</math>。同一种流体意味着运动黏度相同，于是
+<math display="block">\frac{\operatorname{Re}_{\mathrm m}}{\operatorname{Re}_{\mathrm p}}=q\lambda.</math>
+保持雷诺数，必须取 <math>q=1/\lambda</math>。长度缩为十分之一，就需把速度增为十倍。
-温度需要特别区分绝对温度和温差。摄氏与开尔文的绝对读数之间有平移，不能把 20 摄氏度直接当成热力学公式中二十开尔文；而一摄氏度的温差与一开尔文的温差数值相同。简单幂次缩放的量纲推理通常以适当的比例尺度变量为对象，不自动处理所有带原点平移的数值标度。
+若流动还有重力自由表面效应，需要比较弗劳德数 <math>\operatorname{Fr}=U/\sqrt{gL_0}</math>。在同一重力下，它的比值为 <math>q/\sqrt\lambda</math>；保持它不变则要求 <math>q=\sqrt\lambda</math>。长度缩为十分之一时，速度反而应降到约 <math>0.316</math> 倍。
-另一个陷阱是把“无量纲”理解成“必为常数”。振幅角、质量比、长度比都无量纲，仍可以连续变化，未知函数完全可能依赖它们。圆面积若只依赖半径 r，量纲可给出 <math>A=Cr^2</math>，但 C=π 来自几何，并不是量纲抵消自己产生的结果。
+[[File:Gezhi-teaching-similarity-scales.svg|frame|center|alt=长度缩为十分之一的模型在同一流体同一重力下，匹配雷诺数要求速度十倍，匹配弗劳德数要求速度约零点三一六倍|同一个长度缩尺，对两种无量纲数提出不同的速度要求。]]
-== 历史与学习自检 ==
+图中两条箭头从同一缩尺条件出发。上方保持惯性与黏性比，下方保持速度与重力波动尺度的比，它们不能由同一速度同时实现。若优先匹配弗劳德数，雷诺数比就成为 <math>\lambda^{3/2}\approx0.0316</math>；若还要匹配雷诺数，运动黏度也必须降到原型的约 <math>0.0316</math> 倍。[https://ocw.mit.edu/courses/2-20-marine-hydrodynamics-13-021-spring-2005/b014d38899d92006e897830ac4c6e3e2_lecture19.pdf MIT：Marine Hydrodynamics，Lecture 19]
-Buckingham 在 1914 年发表有关物理相似系统的工作，成为 Π 定理历史中的重要节点。[https://www.nist.gov/blogs/taking-measure/life-buckingham-pi NIST 的专题文章]介绍了这一工作及其单位缩放思想。相似、尺度与单位检查早已有多种实践，不能把整个量纲分析传统都视为 1914 年才突然出现。
-=== 定义提供的常数与量纲提供的形式 ===
+量纲分析把这种冲突明确列了出来。实验中需根据主要机制决定匹配哪些量，并分析未匹配量的影响；只有当几何、边界条件和相关无量纲参数相应一致时，才可把模型结果换算到原型。
-若只用密度 ρ、体积 V 来表达某均匀物体质量 m，量纲给出的形式是什么？设 <math>m=C\rho^aV^b</math>，比较 M 指数得 a=1，比较 L 指数得 −3a+3b=0，所以 b=1，得到 m=CρV。对于以质量除体积定义的平均密度，C=1 来自定义；若仅靠量纲，C 的值仍未确定。这个小例子区分了量纲结论与定义额外提供的信息。
-<div class="math-editorial-note">'''编者评注（AI 辅助）。''' 量纲分析的力量在于用很少计算暴露遗漏，而它的边界也应与这种力量一起学习。建议每次先列变量与单位，再构造量纲矩阵，最后明确写出“哪些常数或函数还不知道”。如果一个推导从量纲一致直接跳到某个精确系数，通常缺少了一步物理或数学论证。把这一步补清楚，比多背几个无量纲数名称更有价值。</div>
+== 历史 ==
+Buckingham 于 1914 年发表关于物理相似系统的工作，系统表达了后来称为 Π 定理的无量纲化方法。[https://www.nist.gov/blogs/taking-measure/life-buckingham-pi NIST：The Life of (Buckingham) Pi] 这一框架把已有的尺度与相似思想联系到单位变换，使多变量实验问题能够改写为更少的无量纲关系。
 == 参考来源与延伸阅读 ==