運動媒質中的聲波
Sound Waves in Moving Media
南京大學聲學研究所 王新龍
運動改變声波的頻率,最常見之例當屬高速行駛列車鳴笛聲的音調(頻率)變化。此乃著名的多普勒效應(Doppler Effect)。其實,這是一種波動的運動學現象,涉及媒質、聲源和觀察者(或接收器)三者之間的相對運動。本文僅探討運動媒質中的聲傳播規律以及聲在運動媒質界面上的反射與透射等基本問題。
多普勒效應
試以兩端開口、長為l的聲管為例說明。眾所周知,管內聲波之所以共振,只因聲波的振盪週期恰等於管內聲波來回一個行程所需的時間。假如管與其內流體皆靜止不動,則管的最大共振週期為T0=2l/c0,最小共振頻率(基頻)為熟知的f0=1/T = c0/2l,其中c0為管內媒質的聲速。但是,若管中流體以速度V流動,聲波在管內一個來回耗時

所以,最小共振頻率 f 遂為

式中,f0是管內流體靜止時的共振頻率,M是媒質運動的馬赫數(Mach number),表徵流體的流動相對於聲速的快慢。馬赫數M大於1是超音速。從上式可知,管內空氣的流動使共振頻率f下降。
假設媒質攜聲源以速度V沿x的正向運動,聲源向x軸的正向輻射平面聲波。在相對媒質靜止的坐標系(稱為隨動坐標系)中,聲波的波前當以速度c0向前(x方向)傳播,聲波的頻率為f0=c0/l,其中l是聲波的波長。但是,對於靜止坐標系(或謂實驗室坐標系)中位於聲源正前方的觀察者而言,波前傳播速度是c=c0+V,但所觀察到的波長仍為l。因此,此靜止坐標系的觀察者所測得的聲波頻率是,

此種情形等價於媒質與聲源皆靜止,而觀察者以反向速度(-V)運動之情形。如果僅媒質以速度V沿x正向運動而聲源靜止(即在實驗室坐標系中),則波前在聲源的一個週期T內所行進的距離仍為l=cT=(c0+V)T,故其頻率為f=1/T=(c0+V)/l,仍由上式給出,與實驗室坐標系中的觀察者測量所得到的結果一致。此處須注意,波前傳播的距離(例如波長l)與觀察者運動與否有關。
伽利略坐標變換下的波矢與頻率
推而廣之,設媒質在三維空間以勻速V(三維矢量)運動,方向任意。在隨動坐標系中,空間坐標為r0=(x0, y0, z0),時間為t0。在實驗室(靜止)坐標系中,空間坐標為r=(x,y,z),時間為t。由於V是常矢量,兩個坐標體系通過伽利略變換建立聯繫:

(1)
在此變換下,坐標之間的空間梯度算符和時間微分算符存在以下關係
(2)設聲波的相位為φ,它既可以用隨動坐標系表示,也可以用靜止坐標系表示:
波矢定義為相位φ的負梯度。在實驗室和隨動坐標系下,聲波的波矢分別為,
但根據微商關係(2),恒有k=k0,即波矢具有伽利略變換的不變性。也因此,波長在兩坐標系下是不變的:λ =λ0。
另一方面,圓頻率定義為相位φ的時間變化率。設實驗室坐標系中的圓頻率為ω,隨動坐標系為ω0,
利用微分關係(2),得
利用馬赫數(Mach Number)

(3)
利用波數的不變性,此頻率關係可表為

(4a)
此即推廣了的多普勒公式,其中n是波矢k的方向矢量。設q是波矢k與媒質運動速度V之間的夾角,則上式可表為
(4b)
由此得到波傳播的相速度c = ω/k:
(5)
可見,在流動媒質中聲波的頻率及其傳播速度具有各向異性。當觀察者位於聲源的正前方(θ=0),聲速最大。而在聲源之後(θ=π),c=c0(1-M),聲速最小;對於超音速(M>1)運動的媒質,c<0,意謂在聲源之後的觀察者聽不到聲源發出的聲音。
運動媒質的聲波方程
數學上,聲速的各向異性乃聲波方程的空間對稱性破壞之故。在隨動坐標系r0=(x0, y0, z0)中,媒質靜止,聲壓場遵守聲速c0各向同性的聲波方程
利用公式(2)把隨動坐標變換到“靜止”坐標系r=(x, y, z),則上列聲波方程可在靜止坐標系中表示為
(6) 此即運動媒質中的聲波方程。可見,因媒質的運動(V≠0),波動方程不再保有空間對稱性。設其正向平面波解為代入波动方程(6),同樣可得到频率公式(4)。
誠然,運動媒質的聲波方程(6)也可直接從歐拉運動方程、連續性方程和狀態方程出發導出。惟須注意,因媒質的運動,靜止坐標系中的質點速度v'是質點振速 v(即質點圍繞平衡位置的振速,也即隨動坐標系中質點的速度)與媒質運動速度V之和:
振速v依然是擾動小量,但流速V則不一定。因此,在線性化近似下時間全導數不可僅僅用偏導數近似,即
結果,歐拉方程和連續性方程的线性化近似是
(7)
由此可立即導出聲壓的波動方程(6)。流動管道中的一維聲波傳播
設流體沿管軸(x方向)以速度V均勻流動。在此一維情形下,聲壓波動方程(6)簡化為
設某靜止坐標系中的聲源在管道某處輻射頻率為
ω、波數為
k的平面行波解:
p=
p0 exp(
jω
t-
jkx)。代入上列方程,得到波數與頻率的關係
其中±號表示聲波行進方向與流動同向或反向。此即開篇討論的多普勒效應。
在保持輻射頻率ω不變的前提下,如果媒質靜止(V=0),則波數為λ0=2πc0/ω。因此,運動媒質的波數k±對應的波數λ±為
可見,流動導致了波長的伸長或縮短。這一結論與前述的波長不變性並不矛盾,因為在當前的問題中觀察者和聲源都在同一的靜止(實驗室)坐標系中。
平面行波解
把平面行波解

(8)
代入方程(6),得到如下波數
k與頻率
ω的關係
(9)設隨動坐標系下的頻率為ω0,聲速為c0。由於在伽利略變換(1)下,波數(波長)是不變的,即k=k0=ω0/c0,所以得到實驗室坐標系的頻率ω和隨動坐標系的頻率ω0關係,與公式(4)完全一致。
流動媒質界面上的反射與透射
如下圖所示,媒質特性參數分別為ρ1c1和ρ2c2的兩媒質以z=0為分界面,兩者分別以速度V1和V2平行於界面(設為x方向)流動。
設入射、反射和透射聲壓依次為
注意,此處我們寧願採用與流速方向的夾角作為聲波傳播的方向角,而非常所採用的與法向的夾角作為方向角。根據公式(9),入射、反射和透射波數分別為

(10)
界面兩側聲壓連續的邊界條件可表為
其對任意的x成立,故必有

(11)
如此,聲壓連續方程簡化為

(12)
式中,
rp、
tp分別是聲壓的反射與透射係數。從等式(11)立即可得到以下重要結果

(13)
前者是通常反射定律的翻版——反射角等於入射角,後者是修正的Snell定律。
由於邊界兩側流動的存在,質點法向速度連續的邊界條件不再適用(請讀者自行思考),取而代之的是質點法向位移的連續(在靜止流體的情形下兩者等價)。根據前述線性化的歐拉方程,質點的z方向位移分量ζ與聲壓
p具有如下關係
於是,對應前述的入射、反射和透射聲壓的z方向位移分量依次為

(14)
ζ
1和ζ
2分別入射媒質和透射媒質總位移場的
z-分量。注意,導出過程中利用了公式(10)。根據邊界上ζ
1=ζ
2的邊界條件,得到如下方程

(15)
與方程(12)聯立求解,得到反射和透射係數的表達式

(16)
其中利用了恆等式(11)。
對於非流動的情形(V1=V2=0),公式(13)所給出的Snell定律化為尋常的形式,故而
因此,反射係數為
可見,法向入射的反射係數與流動與否無關。
法向聲阻抗率
在存在媒質流動的情形下,法向聲阻抗率的定義可以推廣為
其中ζ是位移矢量的法向分量。顯然,在無流動的情形下,此定義與通常的定義一致。 譬如,根據z-分量位移的表達式(14),可導出入射波和透射波的法向(z-方向)聲阻抗率分別為
可見,流動改變了行波的法向聲阻抗率。設
z1、
z2分別是界面入射側和透射側(總)的法向聲阻抗率。由於透射側媒質僅存在透射波,故作為入射聲波的負載阻抗,
z2=
zt。由聲壓和位移的連續性可知,法向聲阻抗率在界面兩側必須連續,即
把(12)和(14)的有關公式代入,即可得到如下重要關係
评论