查看原文
其他

量纲分析和量纲制

郑伟谋 中国物理学会期刊网 2022-07-04

|作者:郑伟谋(中国科学院理论物理研究所

本文选自《物理》2021年第12期


摘要  物理定律不依赖于测量单位的选择。量纲分析探讨这种不变性及其后果和应用。无纲量为单位变换下的不变量,物理规律最终必然只能用无纲量表达。从一个问题中的物理变量可构造出的无纲量数要少于原始变量数,带来简化,构造的无纲量可更深刻反映物理量间的内在关系。量纲概念足够深刻,但方法足够简单,应该是物理训练的重要内容。文章阐述量纲分析的基本概念、原理及其应用,大部分内容来自文献,着重讨论量纲制及其与单位制的关系,企图厘正文献中的一些混乱。特别指出,仅就量纲分析操作而言,可以只用MLT量纲制。


关键词  量纲,量纲分析,无纲量,量纲制



物理定律的对称性意味着物理定律在各种变换下的不变性。一个简单原理是物理定律不依赖于测量单位的选择。量纲是物理量不受单位变换而改变的品性,量纲分析探讨这种不变性及其后果和应用。一类特殊的物理量是无纲量,为单位变换下的不变量,物理规律最终必然只能用无纲量表达。从一个问题中的物理变量可构造出的无纲量数要少于原始变量数,带来简化,而且,无纲量可更深刻反映物理量间的内在关系。


量纲分析方法是探讨科学规律、解决科学和工程问题的一个普适工具,非常值得学习和掌握。量纲分析既可以用于实验设计和数据整理,也可以在求解问题前就对问题有个定量和定性的把握,且有助于加深对物理规律的认识。面对复杂问题时,建立数学模型可能非常困难,或者方程非常复杂难以求解,或者难以理解所得解的意义。有时需要做试验,而实际尺寸很难在试验条件中实现,必须缩小尺寸做模型试验,必须满足一定的相似条件,这种条件必须建立在量纲分析和相似论的基础上。


量纲分析很难说是从何时开始的。Dimension一词,1833年泊松首先使用,在此之前用齐次性homogeneity。1822年傅里叶明确表述,物理定律应与单位无关,写在其名著《热的解析理论》(Analytical Theory of Heat)中[1]。这导致一个重要结论:任何有意义的定律,对于每一个计量单位都必须是齐次方程式。很多科学大师如牛顿、欧拉和麦克斯韦等用量纲的概念处理问题。1871年瑞利关于天空蓝色的解释和后来对风中绳弦发声的研究[2,3],1883年雷诺关于雷诺数的工作,都是量纲分析的早期范例。量纲分析的正规形式为Π-定理:如果一个物理关系含有n个独立变数与m个基本量纲,则它可用n-m个无量纲量表示。利用物理定律的量纲平衡(齐次)原理可确定各物理量之间的关系,从而从定性到半定量至定量地解决问题。未查到中文“量纲”的起源,“量纲分析”在日文中称为“次元解析”。关于量纲分析的文献很多,无法一一列举,如中文书[4,5,6],最近又出了一本[7],布里奇曼的英文书[8]是经典,稍晚的可举[9],书中有丰富的文献和习题,还有[10]。


01

从一个故事和一个例子谈起    


传统炸弹,通过在有限的空间里短时释放大量高温高压气体获得机械效能。原子弹是没有伴随气体的强爆炸点源,没人了解其机械效能,当时就有爆炸专家猜测原子弹威力没有期望的那么大。大家认为全英国只有一个人可以解决这个问题,他就是剑桥大学的泰勒(G. I. Taylor)教授。Taylor认为点源强爆炸瞬间释放巨大的能量E,将急剧地压缩和加温其周围的空气,使之以超声速的球形冲击波向外急速膨胀。他列出了问题的流体力学偏微分方程组,但当时无法求解这个非线性方程组。这时Taylor想到了量纲分析,设空气绝热系数为γ,在时刻t爆炸的球形冲击波的波阵面半径为R,波阵面为球形边界面,内部是超热火球,外部是密度为ρ0的正常大气,其压强比球内压强小几个数量级,可以忽略不计。最终这个复杂问题只由5个参量描述,可写为

此处代表某种函数关系。


取质量M、长度L和时间T为基本量纲,上述5个参量的量纲为[R]=L,[E]=ML2T-2,[γ]=1,[ρ0]=ML-3,[t]=T。Taylor知道,5个参量中γ≡Π2已经是无量纲,由其余4个仅可构造出唯一无量纲量Π1=RE-1/5ρ01/5t-2/5,物理过程最终只由这两个量描述,可写Π1=S2),此处S为任意函数,即波阵面半径R=S(γ)E1/5ρ0-1/5t2/5。由之还可得波阵面的速度v(t)=(2/5)S(γ)E1/5ρ0-1/5t3/5,随时间的增大而衰减并趋于零,波阵面的加速度a(t)=-(6/25)S(γ)E1/5ρ0-1/5t-8/5,波阵面处的压强p(t)=C(γ)E2/5ρ03/5t-6/5,等等。1941年6月27日,Taylor给英国有关机构提交了保密报告。1947年美国公开发表第一颗原子弹爆炸从0.10 ms到1.93 ms等时间间隔的14帧火球照片。Taylor根据火球照片在1950年公开发表原子弹的爆炸当量估计:E=7.19×1013 J≈1.7万吨TNT当量,美国军方很是恼火,因为这个爆炸当量当时还是高度绝密的参数。这类推理大概也可用于分析地下和水下核试验。在简化处理时当然忽略了许多细节,如初始爆炸的机械和化学过程、火球并非理想球形、火球受地面限制等。


量纲分析的目的是组织变量成为无纲量并减少变量数,它带来许多好处,首先是节约财力和时间。设想已知一浸没在液流中的物体所受阻力F依赖于其线度L、流体速度v、密度ρ和粘度μF=f(Lv, ρ, μ)。设物体形状复杂,准确求解不可能,只能做试验定F。如果每个自变量取10个值,试验点达104,但量纲分析后的关系约化为F/ρv2L2=g(Re),Re≡ρvL/μ,仅一个自变量即雷诺数。函数g当然不同于f,但包含后者的所有信息。适当的10个试验点可以得到类似精度的g,而且,原先104试验点可以标在同一条g曲线上。第二个好处是帮助试验设计和理论分析,识别重要物理过程,剔除次要变量,例如,雷诺数衡量惯性力和粘滞力的相对权重。第三个好处是提供标度律和相似性,设计模型替代原型,例如,机翼举力试验可保持同雷诺数在小风洞中用小模型做,且知道Fρv2L2。此例稍后还会详细讨论。


02

量纲分析的基本概念和方法    

2.1 量纲制、量纲表示矢量、量纲表示矩阵和Π-定理


物理学是一门高度定量化的学科。物理量离不开测量,测量给出物理量的值。各物理量间由定义或物理定律联系,可选少数几个物理量作为基本量,并为之规定基本量度单位,构成单位制。约定基本量不可能依基本物理定律彼此导出。此处特别将基本量集合称为量纲制,不问其单位;单位制决定量纲制,但反之不然。其他物理量通过定义或基本物理定律作为导出量。量纲“意会”容易“言传”难,基本量不依赖于其单位选择的属性是其量纲,如量纲之一长度,地球半径的量纲是长度。选定基本量和单位制后,导出量的单位可用基本单位表出,这种表达式称为导出量的量纲式,导出量与此相关的不依赖于单位选择的属性是其量纲,描述导出量与基本量间与量值无关的简约关系。例如,在MLT量纲制有三个基本量:质量、长度和时间,对应的三个基本量纲分别记作质量M、长度L和时间T,如L不可能依基本物理定律由MT导出,而速度v=dx/dt是导出量,是长度除以时间的结果,其量纲记作[v]=LT-1,单位改变时v的测量值改变,但量纲不变。物理量的量纲只关注物理量在基本单位改变时不变的属性,不同质的物理量可以有相同的量纲,如功和力矩。


在MLT量纲制中,一般的导出量z的量纲可写作基本量纲的幂次形式[z]=LαMβTγ,此处幂次αβγ也称作因次,早期用于表示量纲。设导出量z在单位000下的值为z0,在单位111下的值为z1,则z=z00α0β0γ=z11α1β1γ,即


X1X2,⋯,Xm是完全的基本量纲,一般量z的量纲可写作

形式上写

可将量纲表示为m维矢量空间中的矢量(a1a2,⋯,am)T,而lnXj为正交基矢,此处上标T记将行矢量转化为列矢量的转置操作。无纲量 ζ 的量纲为[ζ]=X10X20Xm0,约定[ζ]=1,对应于矢量空间的原点。显然,无纲量的值与单位制的选择无关,是不变的常数。


因为物理定律不依赖于测量单位的选择,物理关系式中出现的求和如u+v,求和项的量纲必须一致,即[u]=[v],亦即量纲一致性或齐次性原理,数学上与齐次函数有关。如果函数f(x)对于非零的λf(λx)=λkf(x),则称fk-阶齐次函数。显然,齐次函数的单项式为f(x)=xk。齐次函数的基本定理即欧拉定理:x·∇f(x)=kf(x)。由量纲的齐次原理可引伸出另一条原理即量纲的幂次原理:任何导出量的量纲是基本量纲幂次的单项式,因为导出量的物理定律定义式中各相加项的量纲必须一致。乘积量如z=uv的量纲,如果,则自然。超越函数如ez,sinz有幂级数展开如ez=1+z+z2/2+⋯,要求宗量z一般应为无纲量。以关系式A=Bsinα-(D+D1)/E+F为例,应有[α]=1,[D]=[D1],[A]=[B]=[D/E]=[F]。如果违反量纲一致性,关系式必定是错的。


在量纲的矢量空间表示下,可讨论几个重要推广。首先,一个量的量纲对应该矢量空间中的一个矢量。所谓几个物理量的量纲彼此独立,是指它们的量纲表示矢量线性无关。一组N个物理量的量纲表示矢量,可写成N×m的量纲表示矩阵的形式。其独立量纲数,是指N个量纲表示矢量中线性无关矢量数的最大值,即量纲表示矩阵的秩,不妨将之记作m。进一步还可以讨论基矢变换:m个基矢一般可用m个量纲彼此独立的一组物理量替换,为区分原先的基本量纲量,可称后者为参考量纲量。参考量纲量的量纲表示矩阵必定是满秩的。设参考量纲量为y(1)y(2),⋯,y(m),则参考量纲量的量纲表示矩阵为A=(aij)。以标记分量为lnXi的列矢量,即,相应地有,则。设参考量纲量y(j)在基本单位{i}下的值为y0(j),在基本单位{i}下的值为y1(j),记单位缩放因子i/iγi,则参照(1),有

标记分量为lny0(i)的列矢量,如果选择满足,则lny1(i)=0,即可以通过选择适当的基本单位缩放因子使得所有参考量纲量y1(i)=1,i=1,2,⋯,m。因为矩阵A是满秩的,总可以求解:


如果对于基本量纲量,,而对于参考量纲量,,则由,有,同样,βi可由{αj}通过逆变换得到。设xi是量纲为基本量纲Xi的量,即[xi]=Xi,则为无纲量,即[πz]=1,且可写为

如果z的量纲表示矢量与变量uv,⋯,w的量纲表示矢量线性相关,则z的量纲[z]可用量纲[u],[v],⋯,[w]的有理幂次的乘积表示,量纲式(2)是一特例。设一组n个物理量的量纲表示矩阵的秩为mn,则由此n个量可构造出l=n-m个无纲量。


量纲分析的基础是1914年由Buckingham提出的一条定理[11]。1922年布里奇曼在其专著Dimensional Analysis中将之命名为Π-定理[8]。Vaschy在1892年及Riabouchinsky在1911年曾独立发表了类似工作。布里奇曼由于高压物理的研究获得1946年诺贝尔物理学奖。设想一个物理问题由n个独立变量z(1)z(2),⋯,z(n)完全描述:f(z(1)z(2),⋯,z(n))=0。此处“独立”的意义,一般指n个变量中的任意n-1个取确定值时,余下的一个仍可在一个连续范围内变化。此处的独立性,应与前面的量纲表示矢量的线性无关性明确区分,不可混淆。设以上n个独立变量的量纲表示矩阵的秩为m=n-l,不妨令z(l+1)z(l+1),⋯,z(n)m个参考量纲量,并记作y(1)y(2),⋯,y(m)。参照(3)式,将前lz(k)表作,则

前面已经证明,适当选择单位制可使得所有y(i)的值为1,此时

表明物理问题的函数关系一定可以写成无纲量的函数的形式,此处不同于f,但包含后者的所有信息。这就是Π-定理:含n个独立变量的物理关系式,如果其量纲表示矩阵的秩为m,则物理关系一定可以用l=n-m个无纲量的函数表示。这是物理规律不依赖于单位选择的自然结果,反映了物理定律一个简单却非常重要的不变性原理。此处l个无纲量分别由l个独立变量得到,彼此自然不等价,即无纲量有独立性,或者说,无纲量中l-1个固定时,余下的一个仍可在一个连续范围内变化。


2.2 量纲分析的步骤


量纲分析的步骤如下:

第1步,选择量纲制,列出问题所有的独立关键参量,设共有n个;

第2步,确定所有n个参量的量纲;

第3步,确定量纲表示矩阵的秩m。往往从n个变量中适当选取m个量纲独立的参考量纲量;

第4步,构造出l=n-m个不等价的无纲量πj。往往用m个参考量纲量,对余下的l个参量逐一构造无纲量;

第5步,再次验证所有无纲量的确无量纲,适当整理取某些无纲量的简化组合,使之是科学界通用的已命名无纲量,通常让每个待考察的原始变量只出现在一个无纲量中,或方便于考察极限情形;

第6步,写出问题最终的无纲量关系。在理想情形下,可表现为简单的幂次形式,通常也称之为标度律(scaling law)。


上述步骤给出的问题自由度一般为l=n-m。如果实际问题的n个独立变量中有n1个取固定值,如问题中出现的有纲物理常数,设此n1变量的量纲表示矩阵的秩为m1,则由之可构造出l1=n1-m1个取固定值的无纲量,问题的有效自由度降为n-m-l1


量纲分析的确有些技巧,量纲分析解决复杂问题往往给人一种神奇的美的享受,量纲分析的范例不少是科学大作。如何确定问题的关键物理量,涉及如何简化问题,物理因素有哪些,孰轻孰重,需要有高深的科学修养、丰富的物理知识以及对问题的深刻理解,特别是对于新出现的科学问题。加进关系不大的多余物理量,增加分析的复杂性,但丢失关键物理量导致失败,显然不能解决问题。基本量纲量和参考量纲量的选择有随意性,一种选择会比另一种方便,无纲量的具体形式也可随之不同,但不同的无纲量组彼此等价,最终的无纲量数即物理问题的自由度是不变的。



03

应用例    


定性思考和半定性实验,力求对问题的性质和解的概貌有所估计,这种能力靠物理直觉和洞察力,靠经验和功力,需要长期思考和培养。“学而不思则罔,思而不学则殆。”


3.1 证明勾股弦定理


一个直角三角形的面积A由斜边c和一个锐角α完全决定,各量量纲为[c]=L,[A]=L2,[α]=1。量纲表示矩阵的秩为1,可构造两个无纲量:π0=A/c2π1=α,可写π0=Φ(π1),即A=c2Φ(α),函数Φ的形式可以不知道。将原先的大三角形,如图分割成与原三角形相似的两个三角形,则应有A=c2Φ(α)=a2Φ(α)+b2Φ(α),得c2=a2+b2

勾股弦定理
3.2 管道流动的阻力


考虑圆管道内流动流体所受的管壁摩擦切应力τ,主要的相关物理量为管直径d、管壁线粗糙度ϵ、流体密度ρ、粘度μ和中心流速v,在MLT量纲制中,[d]=[ϵ]=L,[ρ]=ML-3,[v]=LT-1,[μ]=ML-1T-1,[τ]=ML-1T-2,量纲表示矩阵的秩为3,无纲量数为6-3=3。容易证实,可取dρv为参考量纲量。为求对应于τ的无纲量π0,写π0=τdαρβvγ,分别计算M、L、T的总幂次,得量纲平衡方程:

1+β=0,-1+α-3β+γ=0,-2-γ=0,

α=0,β=-1,γ=-2,即π0=τ/ρv2。当然可以写下参考量纲量的量纲表示矩阵,求逆后再系统地构造无纲量。不过,通常不难手调而直接写出无纲量,其余两个为π1=ρvd/μ≡Re,π2=ϵ/d,其中Re为著名的雷诺数。根据Π-定理,可写π0=f(π1,π2),即τ=ρv2f(Re,ϵ/d)。管道流动有层流和湍流不同形式,τ 的行为很不同。


3.3 单摆的周期


考虑重力场中的摆运动,讨论摆周期,除周期T外,关键变量还应有摆长l、摆球质量m、重力加速度g和总能量E。在MLT基本量纲量下,mglTE的量纲分别为[m]=M,[g]=LT-2,[l]=L,[T]=T,[E]=ML2T-2。选mgl为参考量纲量,可验证量纲表示矩阵的确满秩,且M=[m],L=[l],T=[l/g]-1/2。运用Π-定理,由TE构造出无纲量共2个,为Tπ1=T(g/l)1/2Eπ2=E/(mlg),应可写π1=Φ(π2),即T=(l/g)1/2Φ(E/mgl)。引入ω(g/l)1/2,如果π2≪1,可写Φ(π2)=C0+C1π2+⋯ ≈C0,准至最低阶,得T=C0(l/g)1/2=C0/ω,表明周期TmE无关,并不显然。实际的周期由第一类完全椭圆积分给出,为T=4K(k2)/ω,此处参数k取决于摆的能量。


3.4 液滴的振动

忽略重力场,考虑球形小液滴在自身表面张力γ作用下的振动周期T,此处的振动指液滴形状改变。振动周期T依赖于γ及液滴的球半径和密度ρT=f(γrρ)。在MLT量纲制中,各量的量纲为[T]=T,[γ]=MT-2,[r]=L,[ρ]=ML-3,由之可构造唯一的无纲量,得。在此问题中也许需要顾及如粘滞性和可压缩性,液滴足够小时它们也许不重要,当然需要经验和实验检验。


3.5 液体表面张力测量


测液体表面张力γ的一个方法,是测从半径为的管口脱落下来的液滴质量m。基本理论基于平衡状态,粘滞不影响,可略。简单理论认为mg=2π,如果也考虑液滴上方液柱,设管口处液体压强为p,有mg+pπr2=2π。如果紧邻管口处的液滴表面为圆柱形,则根据弯曲液面内外压差公式,p=γ/r,代入上式得mg=πrγ。显然在液滴脱落的瞬间表面不可能为圆柱形,但很难准确计算,不过仍可用量纲分析。此处五个关键物理量为ρgγmr,可构造两个无量纲量γr/mgγ/ρgr2≡(a/r)2,得一般关系γr/mg=ϕ(a/r),此处ϕ为某个函数,实验发现它在a/r的一个范围内为常数,γ=0.263mg/r,相当于mg=1.12π。对于很小和很大的a/rγ的修正因子略小。


3.6 悬臂梁形变


求悬臂梁在尖端加负载P时的形变位移δ。设梁长为L,截面惯性矩为I,弹性模量为Y,且δ=f(PILY)。在MLT量纲制中,[δ]=L,[L]=L,[P]=MLT-2,[I]=L4,[Y]=ML-1T-2,量纲表示矩阵的转置如下表所示,其第一行乘以2与第三行相加得全零行,秩小于3,实为2,无纲量数为5-2=3,无纲量为π0=δ/Lπ1=P/YL2π2=L4/I,根据Π-定理,可写δ=LF(P/YL2L4/I)。考虑梁的小形变正比于负载,则应有,如果进而要求δ∝1/I,得δL(P/YL2)(L4/I)=PL3/YI



悬臂梁的刚度σ由引起梁端发生单位位移形变的力衡量。设长、宽、高分别为lbh,杨氏模量为Y。如果允许剪切,则还需考虑剪切模量μ。在MLT量纲制中各量的量纲为[σ]=MT-2,[l]=[b]=[h]=L,[Y]=[μ]=ML-1T-2,量纲仅出现Lα和组合MT-2,量纲矢量矩阵的秩仅为2,由6个量可构造4个无纲量:π0=σ/Ylπ1=b/lπ2=h/lπ3=Y/μ,得容易看到,梁刚度与宽度成正比,应可分出因子b/l。刚度对高度的依赖不像如此简单,借助其他知识进一步可写为


3.7 LC回路的振荡周期


LC回路的电压变化LC回路的振荡周期τ的有关变量可取电压V、电流I、电感L和电容c在电荷-能量-时间的QET量纲制中,各量的量纲为[τ]=T,[V]=EQ-1,[I]=QT-1,[L]=ET2Q-2,[c]=Q2E-1量纲矢量矩阵的秩为3,可构造无纲量,前者仅是τ和电路常数的组合。于是,可写周期,此处已取无纲量组合以使函数自变量不再含τ。对于稳恒回路,函数Φ取常数,得熟知的结果。在MLT量纲制中,电荷、能量和时间的量纲为[Q]=M1/2L3/2T-1,[E]=ML2T-2,[T]=T,量纲矢量矩阵满秩,所以两个量纲制完全等价。


3.8 星体公转周期


考虑自由空间质量为m1m2的两个星体彼此公转的周期T。记引力常数为G,星体距离为R,设T=f(m1m2GR)。在MLT量纲制中,各量量纲为[m1]=[m2]=M,[G]=M-1L3T-2,[R]=L,可构造无纲量π0=TGm1/R3π1=m2/m1,所以考虑行星绕日运动,如果m1m2,则T2/R3为常数,即行星轨道的开普勒第三定律。此处引力常数进入分析,涉及引力质量和惯性质量转换的问题。


3.9 水面波


决定水面波的因素有二:表面张力使弯曲的表面展平,重力使倾斜的表面恢复水平。设波速为v,波长为λ,水密度为ρ,表面张力系数为γ,重力加速度为g。在MLT量纲制中,[v]=LT-1,[λ]=L,[ρ]=ML-3,[γ]=MT-2,[g]=LT-2。表面张力起作用时,略g,唯一的无纲量对应于。重力起作用时,略γ,得,即正比于自由落体通过一个波长距离所获得的速度,且不出现密度;值得注意,此时3个量vgλ的量纲表示矩阵的秩仅为2。此处未计及水深h的影响,相当于取hλ的极限。考虑无纲量λ2ρg/γ≫1的浅水波情形,量纲法给出。海啸是由地震引起的海洋巨浪,其波长达数百千米,与平均海深比,属浅水波,海啸波长曾是估计平均海深的最早手段。海浪随着海滩水深减小而波速变小,后浪逐前浪,逼高的浪头可达数十米。由量纲法可写,记无纲量ηλ/h,则η→0时,,而η→∞时,,与理论结果一致。


3.10 弹性球的撞击形变


考虑表面涂墨的弹性球沿垂直方向撞击平板弹跳在平板上的印迹。设球有理想弹性,平板完全平滑,质量很大,不移动也不发生形变。关键量除待求的印迹半径d外,有球直径D、撞击速度v、密度ρ及描述球弹性的模量Y和泊松比,质量可由体积和密度算出而不必考虑,最终写d=d(vρDYγ)。可以设想,重力和空气阻力也有影响,但为可略的次要因素。在MLT量纲制中,各量量纲为[d]=[D]=L,[v]=LT-1,[ρ]=ML-3,[Y]=[P]=ML-1T-2,[γ]=1。量纲表示矩阵的秩为3,无纲量数为6-3=3,无纲量为π0=d/Dπ1=Y/ρv2π2=γ,根据Π-定理,可写π0=f(π1π2),即d=Df(Y/ρv2γ),印迹的5个依赖变量骤减为2。量纲分析给不出函数f的具体形式,需要通过理论或实验定f。原始的变量空间5维,工作量大,降为2维后工作量大减,并且,如果实物球的直径很大,现在可用小球做模拟试验,而且可用不同材质的球如铁球或橡胶球做。模拟试验表明,dγ的依赖很弱,即实验点π0π1几乎落在一条曲线上,可写d=DF(Y/ρv2)。


假若关键变量中遗漏了Y,将错误地得出d=Df(γ),后果显然是致命的。不过,如果详细考察实验数据,有可能发现有关键变量遗漏,虽然不能确定缺失的具体变量。相反,如果关键变量中添加了重力加速度g,将导致新增无纲量即无纲重力π3=gD/v2,关系式变为d=(Y/ρv2γgD/v2),问题的复杂性不必要地增加。但是,如果v不足够大使得无纲重力π3不小,也许重力会有可观察的效应,需要检验。总之,简化问题的叙述,突出本质性因素,是解决物理问题的第一步。参考量纲量的选取及构造出的无纲量不唯一,但数目唯一,彼此等价,例如,两个无纲量的乘积仍为无纲量。通常让单个无纲量只代表单个关注的自变量,该无纲量由无量纲化相应的单个原始变量得到,即一无纲量对一有纲变量。


3.11 平行板边界层问题


在无重力场的沿x方向速度为U的均匀流中,放置宽为l(沿z方向)无限长度沿x方向的薄板,流体运动由纳维—斯托克斯(NS)方程和连续性方程描述(ρ(∂tu+u∙∇u)=-∇p+μ2u,∇∙u=0),问题看似简单,实则很难,原因在于其非线性的本质,来自流体与其自身相互作用的反馈机制。1904年普朗特(Prandtl)提出革命性的边界层概念,将流动分为靠近表面的粘滞流边界层和远处的无粘流,雷诺数足够高则边界层很薄。边界层压强项可略,方程为

边界条件对于x>0,y>0为u(0,y)=Uu(x,∞)=Uu(x,0)=v(x,0)=0。相关量为lxyuvUν=μ/ρ,各有量纲[l]=[x]=[y]=L,[u]=[v]=[U]=LT-1,[ν]=L2T-1,可构造无纲量Re≡ρUl/μ=Ul/νπx=Ux/νπy=Uy/νπu=u/Uπv=v/U。无量纲化相当于取U=ν=1,方程和边界条件变为uux+vuy=uyyux+vy=0;u(0,y)=u(x,∞)=1,u(x,0)=v(x,0)=0,此处用下标标记偏导,如ux=∂u/∂x。雷诺数Re作为量纲分析的直接结果出现,衡量惯性力和粘滞力之比。量纲分析指导将原始方程写成无量纲化的形式。


无量纲化后的方程和边界条件仍存在对称性,它们在变换xa2xyayuuvv/a下不变,表明x/x0=(y/y0)2。引入,令u=ψy=f′,则由连续性方程,NS方程变为2f′′′+ff′′=0,边界条件成为f(0)=f′(0)=0,,整个问题化为单参量。



04

量纲制的基本量个数    


在CGS单位制或高斯单位制中基本单位只有3个,在MKSA单位制中有4个,而在国际单位制SI中有7个。基本单位的个数,就是基本量纲量的个数。增加基本量纲量如引入温度或电流,便于与力学过程相区分,分析热学或电学过程所起的作用,但增加基本量纲量的个数,根据Π-定理,似乎将减少无纲量的数目。物理问题的自由度显然不会因单位制的选择而改变,增加基本量纲量一般也增加相关的物理量,如引入基本量温度后需要相应引入有纲量玻尔兹曼常数,问题的自由度最终不变。仅就量纲分析操作而言,可以只用MLT量纲制。例如,对于SI单位制的公式,在MLT量纲制中电磁学常数量纲可取[ϵ0]=1,[μ0]=[c-2]=T2L-2。物理量如空间矢量,量纲制或单位制如坐标系,物理量不依赖于单位制,但物理量的表示随单位制而变。物理量的量纲不依赖于量纲制中基本量单位的选择,但依赖于量纲制如MLT或MLTΘ(Θ表示温度)的选择。


4.1 流场中物体的传热


量纲分析的一位先驱瑞利在1915年讨论流体传热问题:在匀速运动的不可压无粘性流体中放置一物体,其与远处流体的温差为θ,求单位时间内由物体传出的热量H。瑞利认为关键物理量有6个,为H及物体线度l、远处流速v、温差θ、流体单位体积热容c和热导率κ,选五元量纲制:除MLT外另加温度Θ和热量Q,有[l]=L,[v]=LT-1,[θ]=Θ,[c]=QL-3Θ-1,[H]=QT-1,[κ]=QL-1T-1Θ-1。量纲表示矩阵的秩为4,无纲量数为2,得无纲量π1=H/lκθπ2=lvc/κ,所以H=lκθΦ(lvc/κ)。瑞利的推理引起争议。其实随着新增的基本量纲量温度Θ和热量Q,引入了新关键量即热功当量J和玻尔兹曼常数kB,有[J]=ML2T-2Q-1,[kB]=ML2T-2Θ-1,量纲表示矩阵的秩变为5,无纲量数为8-5=3,新增无纲量π3=Jcl3/kB。如果只用MLT量纲制,以等价的kBθ取代θ,则[l]=L,[v]=LT-1,[kBθ]=ML2T-2,[c]=L-3,[κ]=L-1T-1,[H]=ML2T-3。量纲表示矩阵的秩为3,无纲量数为6-3=3,得无纲量π1′=H/kBθπ2′=v/l2κπ3′=l3c,不难看出π2=π2π3′,现在H=kB(lκθ)Ψ(v/l2κcl3),与瑞利分析修正后的结果一致。两种量纲制下量纲表示矩阵的转置如以下二表所示。



4.2 管道流动的传热


考虑圆管道内稳恒流动流体的每单位温差下单位时间内通过单位面积的传热量即传热率h问题,主要的相关物理量为管直径d、管长l、流体密度ρ、粘度μ、热导率κ、平均流速v及温差θ 和比热c,在MLT量纲制中,[h]=L-2T-1,[d]=[l]=L,[ρ]=ML-3,[μ]=ML-1T-1,[v]=LT-1,[θ]=ML2T-2,[κ]=L-1T-1,[c]=M-1,量纲表示矩阵的秩为3,无纲量数为9-3=6,对应h的无纲量π0=hd/κ≡Nu,称为Nusselt数,π1=ρvd/μ≡Re,是雷诺数,π2=/κ≡Pr,是Prandtl数,描述粘滞和传热的相互影响,π4=d/l,描述几何形比,π5=/v2对应θ,π6=κl2/v。根据Π-定理,可写Nu=f(Re,Pr,d/l/v2κl2/v),即τ=ρv2f(Re,ϵ/d)。量纲分析直接给出了重要的无量纲数。在MLTΘ量纲制中,π6对应于新增的有纲玻尔兹曼常数kB


4.3 饮料加冰冷却


热饮料加小方冰块冷却比加同质量的大方冰块快。通常的解释是:小冰块的总表面积大,传热多,所以冷却快。但是,线度为l的冰块面积∝l2,总体积V的冰块数~V/l3,冷却速率应该∝1/l,得冷却时间Tl。不过通常的解释并不准确,量纲分析的结果是Tl2。问题的有关量除lT外,为温差θ、冰的热导率κ和饮料体积比热c。在MLT量纲制中,[l]=L,[T]=T,[θ]=ML2T-2,[c]=L-3,[κ]=L-1T-1,量纲表示矩阵的秩为3,无纲量数为5-3=2,无纲量为π0=Tlκπ1=l3c,根据Π-定理,可写T=(1/lκ)F(cl3)。考虑cl3≪1,对F展开至线性项,得T∝(1/)cl3=cl2/κ。如果采用MLTΘ量纲制,则引入新有纲量玻尔兹曼常数kB,各量量纲为[l]=L,[T]=T,[θ]=Θ,[c]=ML-1T-2Θ-1,[κ]=MLT-3Θ-1,[kB]=ML2T-2Θ-1,量纲表示矩阵的秩为4,无纲量数为6-4=2,无纲量为π0′=/l2cπ1′=l3c/kB,根据Π-定理,可写T=(l2c/κ)Ψ(l3c/kB)。不难看出,π0′~π0/π1π1′~π1,不同量纲制的结果等价。两种量纲制下量纲表示矩阵的转置如以下二表所示。代谢率与动物体重的3/4次幂呈正比,即克莱伯的实验总结定律(Kleiber's law),而德国生理学家鲁伯纳根据表面散热速率得出的却是2/3次幂,不符合实际。加冰冷却问题中的幂次估计的不一致性,也许提供有益的启示。


05

理论物理应用    

5.1 黑体辐射


自由电磁场的辐射压p与辐射能密度u间有关系。玻尔兹曼运用热力学由内能U=3pVu=(∂U/∂V)T=T(∂S/∂V)T-p=T(∂p/∂T)V-p=,推出斯特藩定律u=σT4,且给出熵。1893年维恩用当时不很流行的热力学方法处理黑体辐射。维恩设想提高辐射谱密度的两个独立过程,一为升温,另一为绝热压缩,认为只要最终温度相等,则二者谱密度分布相同。维恩考虑半径为R体积为V的球腔以速度v=dR/dt作准静态绝热压缩,由等熵=常数,得VT3RT为常数。运用多普勒频移和速度的关系,他得到dλ/λ=dR/R。记绝热压缩dR前后的温度分别为T0T,则λT=λ0T0。将能量密度u=σT4分配到各波长,应可写u(λT)=T5f(λT),或者(νT)=T3F(T/ν)。由量纲分析看,函数F的宗量应取无纲量,需要引入普朗克常数h和玻尔兹曼常数kB,将T/ν改造成无纲量/kBT。另外,。取νchkBT为独立变量,可构造另一无纲量/3c-3,最后,(νT)=h(ν/c)3F(/kBT),与维恩公式一致,由之对ν积分得斯特藩的T4定律,不同于维恩,此处直接由量纲分析得到。
Jeans曾详细考察过黑体辐射问题,他设想黑体腔内充稀薄电子气,保持温度θ,电子间彼此作用,也与辐射场作用,电子吸收辐射又发射辐射,并达平衡。辐射场能量密度u,应与电子质量m、温度θ、电荷e和光速c有关,而与电子数密度无关,因为电子平均动能仅取决于温度。在MLT量纲制中,各量的量纲为[u]=ML-1T-2,[m]=M,[e]=M1/2L3/2T-1,[θ]=ML2T-2,[c]=LT-1,可构造两个无纲量π0=ue6θ-4π1=θ/mc2。于是,u=e-6θ4F(θ/mc2),含任意函数F。函数F的自变量π2~v2/c2,相当于电子热速度与光速之比,通常温度下很小,有理由写ue-6θ4,与斯特藩定律一致。黑体辐射问题离不开量子常数h,但以上推导不出现h。斯特藩系数隐含h,以上给出的斯特藩系数含e,显然需要小心解读。


5.2 原子半径


考虑氢原子,关键量应有电子质量me、电荷e及真空介电常数ϵ0,实际上要定出原子半径还需普朗克常数h。为方便,将SI的MLTI量纲置换成MLTQ,即电流换成电荷,此时[e]=Q,[me]=M,[ϵ0]=M-1L-3T2Q2,[h]=ML2T-1。由此四量可构造长度量r=ϵ0h2/mee2,但仅me,e,ϵ0不足以构造出长度量。如果添加光速c,因为[c]=LT-1且由前四量可构造时间量t=ϵ0h3/mee4,得无纲量e2/ϵ0hc,对应精细结构常数αe2/4πϵ0hc,长度量的一般形式取r=ϵ0h2/mee2Φ(α),对于非相对论情形,可取Φ(α)=1,得前面的结果。当原子特征速度与光速相比不是太小时,相对论修正结果含光速c。在MLT量纲制中不必考虑ϵ0,[e]=M1/2L3/2T-1,[me]=M,[h]=ML2T-1,得长度量r=h2/mee2,与上述长度量一致。


5.3 电子经典半径


如果电子的静能量mec2完全来自电磁能,则电子有多大?这是电子经典半径问题,结果是re=e2/4πϵ0mec2≈2.8×10-15 m,但实验值小于10-22 m,即量子半径远小于经典半径。前面知道,仅meeϵ0不足以构造出长度量r,或者说仅reϵ0不足以构造质量me。电磁相互作用由光子传递,加上光速c很自然,由量纲平衡得估计me=De2/ϵ0rc2,与上述电子经典半径式一致。电子半径实验值比经典估计小7个量级,指示量子处理的必要性。


5.4 静电场能量


自由静电场能量的有关量包括电荷e、电场强度E和待求的能量密度u。在MLT量纲制中,各量的量纲为[e]=M1/2L3/2T-1,[E]=[e/r2]=M1/2L-1/2T-1,[u]=ML-1T-2。量纲矢量矩阵的秩为2,可构造唯一的无纲量π0=u/E2,所以uE2,实际的比例常数为1/2。如果量纲制增加基本变量电荷Q,则引入新变量即真空介电常数ϵ0,在MLTQ量纲制中,各量的量纲为[e]=Q,[E]=QL-2,[u]=ML-1T-2,[ϵ0]=Q2T2M-1L-3。现在量纲矢量矩阵的秩变为3,但仍可构造唯一的无纲量π0′=u/ϵ0E2,与π0等价。如果考虑在介质中,则需计及介质介电常数ϵ,因为[ϵ]=1,结果含ϵ 的任意函数,须另外确定。


5.5 理想气体的压强


气体分子运动论认为压强p来自分子对器壁的碰撞,理想气体分子无大小,设分子质量为m,数密度为n,温度θ衡量分子的运动能量。在MLT量纲制中,各量的量纲为[p]=ML-1T-2,[m]=M,[n]=L-3,[θ]=ML2T-2,量纲矢量矩阵的秩为3,可构造唯一的无纲量π0=p/,得p。如果用MLTΘ量纲制,各量的量纲为[p]=ML-1T-2,[m]=M,[n]=L-3,[θ]=Θ,量纲矢量矩阵满秩,不存在无纲量。但实际上Θ引入有纲量玻尔兹曼常数kB,且[kB]=ML2T-2Θ-1,量纲矢量矩阵的秩为4,可构造唯一的无纲量π0=p/nkBθ,得一致结果pnkBθ


5.6 宏观性质和原子间力


现在分析物体宏观性质可能如何依赖于原子间力。设原子间力取f=Br-n-Ar-2,此处n>2固定,AB为正但对不同材料可变,A项描写远程吸引,B项描写近程排斥。假定关键量有压强p、温度θ、原子质量m及力系数AB。在MLT量纲制中,各量的量纲为[p]=ML-1T-2,[θ]=ML2T-2,[m]=M,[A]=ML3T-2,[B]=MLn+1T-2,可构造两个无纲量。于是,可写πp=F(πθ),压强p和温度θ不再独立,此处的任意函数F对于不同的AB仍然同一,πp和πθ相当于范德瓦耳斯的对应态,是约化压强和温度。对于任何其他物理量q,如热导率,可用mAB无纲量化成πq,并可写πq=Φ(πpπθ)。
考虑固体比热的一个问题:决定固体弹性的原子间力是否同样决定红外光学频率。为分析此问题,考虑特征频率ν、压缩率κ、原子的质量m和体积数密度n,在MLT量纲制中,各自的量纲为[ν]=T-1,[κ]=[1/p]=M-1LT2,[n]=L-3,[m]=M,可构造唯一的无纲量2κn1/3,得2κn1/3=常数。以铜的实际值估计,常数值为0.18,仍处在1上下,可视作合理。
如果企图建立引力的电磁理论,则尝试考虑引力常数G、电子电荷e、质量m和光速c,在MLT量纲制中,各自的量纲为[G]=M-1L3T-2,[e]=M1/2L3/2T-1,[m]=M,[c]=LT-1,由唯一的无纲量得G∝(e/m)2。以实际值估计的比例常数为2.4×10-43,极不合理,表明应该放弃将引力归于电磁的想法。虽然量纲法不能确定比例因子,但实际估算值可指示分析的合理性。任何物理关系应该反映量纲关系,量纲分析可考察物理关系的真实性。


5.7 范德瓦耳斯方程和对应态


范德瓦耳斯方程p=NT/(V-Nb)-a(N/V)2是描述临界现象的一个模型。在三相点附近,密度很高,吸引势相对于陡峭的排斥势而言很弱且几乎彼此相消,每个分子可看作处于均匀的负势背景中。此背景势应与密度成正比,对内能有贡献-aN2/V,相应有压强修正-a(N/V)2,得范德瓦耳斯方程。虽然范德瓦耳斯方程最初目的在于说明临界现象,但却更适于描述三相点附近。范德瓦耳斯方程可粗略地描述气液凝聚。低温的等温线有一个最大和一个最小。上升温度至一个特定值,等温线的最大和最小重合为一点,此点满足拐点条件:

第一个方程表明压缩率无限。解此二方程可得

此处Tc即为临界温度。引入无量纲约化变量

可得对应态定律:

即范德瓦耳斯方程在临界温度附近有“普适”形式,与体系的物质构成细节无关,其意义深刻:即使范德瓦耳斯方程不可用,用对应态的无量纲量处理实验数据如将不同实验的数据绘制在一张图上,可揭示普适类。对应态强调物理问题中出现的各种特征尺度。



06

有纲物理常数的压缩与恢复    


物理常数如重力加速度,常常以关键物理参量的形式进入物理问题。例如,考虑重力场中垂直悬挂弹簧末端质量为m的质点的运动,设弹簧的自然长度为l,劲度系数为k,求弹簧长度x(t)随时间t的变化。牛顿方程为md2x/dt2=-k(x-l)+mg,初条件为x(0)=l,dx(0)/dt=0。物理方程中常数的压缩,相当于经常看到的“可以选择单位使得m=g=k=1”的说法,此时牛顿方程变为d2x/dt2=-x+l+1,初条件仍为x(0)=l,dx(0)/dt=0。量mgk的量纲分别为MLT-2,MT-2,容易验证其量纲表示矩阵满秩,且M=[m],L=[gm/k],T=[m/k]1/2。常数压缩过的方程解为x=l+1-cost。由之不难恢复物理常数。常数压缩实质上是作无量纲化。因为相对于常数压缩的单位选择,总有mαgβkγ=1,不管αβγ取何值,所以任何量可乘以这样的因子,只是需保证量纲正确。例如,由解给出的频率ω=1,因[ω]=T-1ω=1意味着ω=(k/m)1/2。考虑到[x]=L,常数恢复后,,满足要求:量纲正确,且m=g=k=1时退化为约化解。
现在再就单摆运动讨论物理常数的压缩与恢复。适当选择单位使得m=g=l=1,则体系的拉氏量退化为,运动方程退化为,动量为。至于物理常数的恢复,如由[ω]=t-1,得
谐振子的拉氏量,适当选择单位使得m=k=1,则体系的拉氏量退化为,运动方程退化为,解~e±it,频率ω=1。考虑到[m]=M,[x]=L,[k]=MT-2,[ω]=T-1,常数恢复后,满足量纲要求,且m=k=1时退化为约化解ω=1。在氢原子的量子问题中,适当选择单位使得m=ℏ=kee2=1,则体系能级为
物理常数压缩的范例,是基于三大物理常数ℏ,Gc的普朗克自然单位:ℏ=G=c=1。在MLT单位系中,[ℏ]=ML2T-1,[G]=M-1L3t-2,[c]=LT-1,所以,,分别有长度、时间和质量的量纲,为普朗克基本长度、时间和质量,即所谓的绝对单位制。出现在公式中的物理量,实际上是相对于这三者的无纲量。有纲常数如引力常数通过自然定律引入,值随单位的选择而变。将与普适自然定律联系的基本常数取作1,可以反过来将度量单位决定于自然定律。选择哪些常数或定律,未必唯一;普朗克选择常数G=c=h=1。有时人们要求太多,如追问如此定出的长度单位是否有本质性的意义。普朗克的方案可以扩充,如另选玻尔兹曼常数kB=1,引入温度单位。


07

相似性理论和模型试验    


量纲分析可用于复杂问题,通常描述问题的方程未知,更不消说解,但量纲分析将问题归结为少数无纲量的任意函数。工程技术领域中常用模拟试验代替真实试验如小模型风洞试验,经济又方便,且可超越实际限制。模型试验包括用模型动物进行的医学试验。1868年,弗劳德(W. Froude)彻底改革船舶设计,提出弗劳德数(v2/gl)用以指导小模型试验。他曾被嘲笑说,小模型也许带来无穷的乐趣,但没有任何实际意义,不过工程界最终转变了态度。模型试验的理论基础是基于量纲分析的相似性理论。虽然量纲分析中出现的任意函数的具体形式未知,只要保持无纲量不变,则模型的行为与原型等价。如果适当选取无纲量的组合,让它们在原型中固定,则所涉及函数的自变量数可以减少,模型试验可简化。当然,此处的相似性未必是几何相似性。只有同量纲量可以比较,体系间的物理相似性由无纲量刻画。另外,相对于描述问题的方程,定解条件如初条件和边条件的相似性应该一并考虑。


7.1 流体中运动物体的阻力


以在流体中运动的物体受阻问题为例讨论,特例包括空中飞行的炮弹或飞机、下落的雨滴和水下潜艇,物体的形状复杂,但可明确给定。假定重力影响可忽略,除待分析的阻力τ外,至少应考虑物体的速度v、代表线度l及流体密度ρ、粘度μ、声速vph,还有描述形状的各种形比参数,仅以单个符号γ记。形比为无纲量,在MLT量纲制中,其余6个量的量纲为[τ]=MLT-2,[v]=[vph]=LT-1,[l]=L,[ρ]=ML-3,[μ]=ML-1T-1,可构造出无纲量π0=τ/v2l2ρπ1=ρvl/μ≡Re,π2=v/vph,一般关系可表作π0=F(π1π2γ),即τ=ρl2v2F(π1π2γ)。
π2≪1的低速情形下,π2影响可略,要求阻力正比于速度,得π0π1=f(γ),即τ=vlμf(γ),具有斯托克斯阻力公式的形式,相当于1/π0的小Re展开。此时保持完全几何形比相似的条件下,模型试验只有单参数。再增大速度,虽然π2的影响仍可略,但密度的效应已需考虑,此时应取π0=f2(π1γ),即τ=v2l2ρf2(Re,γ)。保持完全几何形比相似的条件下,模型试验还需保证雷诺数一致。再继续增大速度,物体速度和流体介质声速可比,此时介质难以逃脱物体,粘滞效应次要,可设想π0=f3(π2γ),即τ=v2l2ρf3(v/vphγ)。此时模型试验受到严重限制。
如果考虑斯托克斯低速小球的阻力与重力平衡达到稳态速度,则τ∝(ρ0-ρ)d3g,相当于变量τ g替代,此处ρ0为物体密度,此时π2π2=ρ0/ρπ0π0=(ρ0/ρ-1)dg/v2。流体力学公式为,相当于,与Π-定理不矛盾。
7.2 管道泊肃叶流动


考虑水平圆管道内流体的流动,主要物理量为管直径d、流体密度ρ、粘度μ、平均流速和压强梯度γ≡Δpl。在MLT量纲制中,[d]=L,[ρ]=ML-3,[μ]=ML-1T-1,[]=LT-1,[γ]=ML-2T-2,可构造两个无纲量,前者是流体力学中非常重要的无纲参数雷诺数。记管道截面Sd2/4,则根据Π-定理,可写π2=Φ(π1),进而管道截面总压力梯度ΓγS=函数Φ(Re)称为管道阻力系数。流量,由Reπ2=γd2S/μQv,得Qv=πγd4/4μReΦ(Re)≡(γd4/μ)Ψ(Re)。
管道流动问题已归结为求函数Φ(Re),实验测定Φ时,可用不同的流体、流速、密度、管径和压差,实验数据按π2π1=Re应落在同一条曲线上。管道流动有层流和湍流不同形式,函数Φ(Re)相应有两支,小Re支对应层流。
泊肃叶流考虑的问题很类似,主要物理量为管直径d、流体密度ρ、粘度η、平均流速、管长l和压降ΔP物理量数仍为6,现在[ΔP]=ML-1T-2,[l]=L,量纲表示矩阵的秩仍为3。无纲量除Re=ρdv/η外,另有ΔP/ρv2d/l最终ΔP=ρv2Φ(Re,d/l)。


08

标度律    


伽利略曾明确地表述:“体形越小,其相对强度越大。因此,一只小狗能够背负两三只与自身同等大小的狗,但我相信,一匹马连一匹与自身同等大小的马都驮不了。”这些动物的机体构成相似,密度ρ及其他材料性能彼此相似。衡量材料机械强度的量是体模量B,量纲如压强即如单位面积受力。动物所受的重力与体积即线度l的立方成正比,承重的截面积Al2,应该将ρgl3/l2=ρgl与体模量比较,强度要求ρgl<B,小动物胜出。


鸵鸟是世界上最大的鸟,不会飞是因为翅膀退化吗?飞翔的必要条件是举力克服体重W,按高雷诺数估计的举力f=CSρv2,式中S为翅膀面积,ρ为空气密度,v是起飞速度,C为无纲比例因子,得起飞条件。从简单的几何相似性出发,因ml3Sl2,有。燕子的最小飞行速度约20 km/h,鸵鸟的l约大25倍,临界起飞速度应大5倍为100 km/h,鸵鸟奔跑不慢,但最高速度不过40 km/h,不足以起飞。


狭义的标度律,通常表现为单项式。标度律,从严谨的量纲分析来看,是必要而未必充分,相当于简化版的量纲分析。如果问题足够简单且有能力把握问题的关键量,使得无纲量只有一个,此时标度律是确定的。从完整的量纲分析出发,再区分退化情形,得到标度律,才是稳妥的做法,可以参考前面关于饮料加冰冷却的例子及流体中运动物体的阻力的例子讨论。标度律,几何上对应于分形。标度律被用于描述各种复杂现象,如前面提到的代谢率与动物体重的克莱伯3/4次幂定律。


8.1 承重圆柱的标度律


建筑物的承重圆柱,其重量随着体积增长而增长,其构成材质的机械性能变化不大。设柱长为l,截面半径为r,密度为ρ,弹性模量为Y。重压下发生侧弯的弯弧半径为R,张角为θ,则θ~l/R,柱高相对变化为Δl/l=[-2Rsin(θ/2)]/~θ2。将承重侧弯的圆柱看作是垂直放置的弹簧,则弹性能~θ2。弹性曲线的能量正比于平方曲率沿曲线的积分R-2l,而Y与压强同量纲,由量纲分析给出的圆柱形变弹性能~Yr4θ2/l。仅考虑自重,形变引起的重力势能改变为mgΔl~r2lρg·2=r2l2ρgθ2。稳定条件要求形变弹性能大于重力势能改变,γYr4θ2/l>r2l2ρgθ2,此处γ为无纲比例因子,临界柱长为lc=(γYr2/ρg)1/3~r2/3,表明柱长增大比半径慢。各量量纲[l]=[r]=L,[θ]=1,[ρ]=ML-3,[g]=LT-2,[Y]=ML-1T-2,无纲量π0=Y/ρglπ1=θπ2=r/l,为得到上面的标度律,需要一些额外的知识。梁柱的负载由其横截面积支撑,在构成材质不变的条件下,其强度取决于其横截面积,随着线度增长强度变弱,需增大横截面以保证强度。人们往往习惯于线性思考,忘记标度律。树干的半径和高度,或者动物四肢的半径和长度,包括心率与代谢率等,会随体形的增长而发生不同比例的变化,1/4-次幂标度律是生物学的普遍特征[12]
8.2 大气声波:等温或绝热?


声波中压强和密度在变,意味着温度也在变。如果运动足够快,来不及传热,过程绝热;如果运动足够慢,可充分传热,过程等温,前者和后者似乎分别对应高、低频。然而,实际并不如此简单,高频波的波长短,传热距离也短,有利于等温,低频波有相反趋势。量纲分析可解决直觉上的矛盾。傅里叶定律描述热传导:q=-κT,热量守恒:,联合得热扩散的方程为∂T/∂t=D2T,此处扩散系数D=κ/ρCv有量纲[D]=L2T-1。声波中波长λ和频率ν是关键量,分别有量纲[λ]=L,[ν]=T-1,可构造无纲量πD/λ2ν。如果D很大,表明热扩散很快,声波等温,否则绝热。但是D为有纲量,无法简单说大小,必须以无纲量π衡量大小,决定过程是否等温。大气声速约340 m/s,且色散不显著。因为声速vph=λν,等温判据可写作/vph2>1,或ν>vph2/Dνc。大气的νc~1000 Hz,频率在此以下的声波绝热。
8.3 瑞利的蓝天解释


最早解释蓝天的是Tyndall,他做红外实验需要极清洁的空气,发现不洁空气中的尘埃会散射光,打进一束强光可检测空气是否清洁。他还发现散射光偏蓝,透射光偏红,并大胆联系到蓝天解释。Tyndall的解释不完全对,日光非白光,接近黑体谱,散射体主要不是尘埃,看到的不是波长更短的紫色,这与人眼的光敏感性有关。瑞利进一步审视了Tyndall的解释。光如何被尘埃散射?此处不涉及光是波或是粒子的量子问题,即可将之看作经典电磁波。受电磁波驱动,原子如受激偶极子,即原子吸收电磁场能量而振动,又在各个方向再发射辐射。这是瑞利散射,是原子或小粒子对光的一种弹性散射,粒子尺度可只比波长小一量级。设幅度Ei均匀的电磁波投射到体积V~L3上,如果线度Lλ,远小于波长,整个体积将同相位受激振荡,有散射幅Es~L3。散射强度Is=Es2,在各方向随距离r以1/r2衰减:Is~L6/r2。考虑到EsEi及量纲平衡,应有,得瑞利的ν4-律。白云,看起来灰白,其水滴尺寸比分子大得多,散射以米氏散射为主,不大依赖于波长。
8.4 科尔莫戈罗夫湍流能谱


1941年科尔莫戈罗夫(Kolmogorov)发表了各向同性稳态湍流能谱即湍流能量依湍涡尺度分布的公式。虽然未必与实验一致,却是湍流理论的奠基性工作。记相应于湍涡尺度l的波数k=2π/l,湍涡k1k2k3的能量转移满足选择定则k3=k1+k2,有两类极端情形:k3≈2k1k3=k1+δ≈k1,Kolmogorov考虑前者。在所谓的波数惯性区内耗散不重要,Kolmogorov公式适用于惯性区。设湍涡尺度为d,特征速度为u,则特征时间τd/u。稳态下能量在不同尺度间的转移率ε应对所有尺度等同,提供给dmax的能量等于在dmin耗散的能量,且[ε]=[E/mt]=L2T-3。湍涡速度u应只依赖于dε,由量纲法得u~()1/3,特征时间τ~d2/3ε-1/3。最小湍涡取决于粘性ν,因为湍涡越小粘性越重要。由[ν]=[μ/ρ]=L2T-1,得dmin~ν3/4ε-1/4。体系尺度L~dmax,相应的速度U~(εL)1/3,得长度标度比dmax/εdmin=L/3/4ε-1/4=(UL/ν)3/4=Re3/4,此处推导用到ε~U3/L,进而速度标度比umax/umin~(εdmax/εdmin)1/3=Re1/4。定义E(k)为能量密度,即E(k)dk有能量/质量的量纲,所以[E(k)]=L3T-2E(k)~k-3/5ε2/3


09

结 语    


物理定律不依赖于测量单位的选择,是最基本也最简单的对称性。物理规律最终必然只能用无纲量表达,由此带来独立有效物理量数的约化即降低问题的自由度,且更深刻反映物理量间的内在关系。通过简单分析一个问题中各物理量的量纲,可以给出问题的有效变量数目和可能形式,限制这些物理量间的某种可能关系。虽然量纲法不能预言明确的函数关系,但其结论也表现出与任何具体的函数形式无关的普适性。量纲分析方法是探讨科学规律、解决科学和工程问题的一个普适工具,非常值得学习和掌握。但十分奇怪的是,虽然量纲分析的原理极其简单,如此重要的量纲分析方法,却一般不在物理课程中系统讲授。量纲分析的范例不少是科学大作。量纲法不能事先确定关键物理量是哪些,决定量纲分析成败的,是如何确定问题的关键物理量,涉及如何简化问题,物理因素有哪些,孰轻孰重,需要有高深的科学修养、丰富的物理知识以及对问题的深刻理解,特别是对于新出现的科学问题,如处理涉及生命的现象。一旦单位制确定,量纲制也确定,但反之不然。仅就量纲分析操作而言,可以只用MLT量纲制。一个问题的有效变量数目和形式,不会因量纲制不同而有本质改变。量纲法属于定性和半定量方法,许多物理人理解定性方法似乎是不得以而为之,其实数学上的定性方法回答有无,往往并非定量方法可比。量纲分析在各个领域均有应用,流体力学是传统领域,生物医学则是新兴领域。本文的大部分内容只是一些文献的整理和总结,也企图厘正文献中的一些混乱。


参考文献

[1] Fourier J B J. Analytical Theory of Heat. New YorkDover Pub.1955

[2] 谈庆明. 量纲分析. 合肥:中国科技大学出版社,2007

[3] 孙博华. 量纲分析与Lie. 北京:高等教育出版社,2016

[4] 赵凯华. 定性与半定量物理学. 北京:高等教育出版社,2007

[5] 梁灿彬,曹周健. 量纲理论与应用. 北京:科学出版社,2020

[6] Rayleigh J S WStrutt B J W. The Theory of Sound. MacMillan1877

[7] Rayleigh J S W. Nature19159566

[8] Bridgman P W. Dimensional Analysis2nd ed. New HavenYale University Press1931

[9] White F M. Fluid Mechanics4th edCh. 5 Dimensional Analysisand Similarity. McGraw-Hill College1998

[10] Sonin A A. The Physical Basis of Dimensional Analysis2nd Ed. CambridgeMIT2001. http://web.mit.edu/2.25/www/pdf/DA unified.pdf

[11] Buckingham E. Phys. Rev.19144345

[12] West S G 著,张培 译. 规模:复杂世界的简单法则.中信出版集团,2018





《物理》在淘宝店和微店上线,扫码即可购买过刊和现刊。


  

▼往期精彩回顾▼


1.黑体辐射公式的多种推导及其在近代物理构建中的意义(Ⅰ)

2.怀念鲍家善先生 | 都有为
3.多拓扑荷“磁束子”的发现
4.无愁河上的浪荡汉子们|量子多体中的呐喊与彷徨之十一
5.二维材料高次谐波及光波电子学
6.强场亚周期光脉冲研究
7.爱因斯坦:我们为什么研究物理学?
8.细推物理须行乐—— 祝贺李政道95华诞

9.何处不关联?

10.30年前预言的奇怪量子效应,终于被证实

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存