船隻尾波的鳥瞰圖
從前方觀測的船隻尾波
尾波 (英語 :wake )是固體 在划過流體 (特別是液體 )表面時在尾部產生的V形傳播的波 ,例如水鳥 或船舶 勻速游過水體 時在水面激起的後方波紋。因為由英國 的開爾文男爵 ——物理學家 威廉·湯姆森(William Thomson,1824~1907)最先對船波進行數學研究,因此也稱為開爾文船波 (Kelvin wake或Kelvin ship wave )。
船形物體的尾波形狀和福祿數
F
r
{\displaystyle Fr}
有密切關係。
F
r
=
V
g
l
{\displaystyle Fr={\frac {V}{\sqrt {gl}}}}
其中g 為重力常數,V 是船速,l 是船的長度。
令船的長度
l
=
k
⋅
V
2
g
{\displaystyle l=k\cdot {\frac {V^{2}}{g}}}
則
F
r
=
1
k
{\displaystyle Fr={\frac {1}{\sqrt {k}}}}
.
對於長度大而速度低的輪船,Fr數小,開爾文船波主要是長波,其波前與速度矢量的夾角比較小。
而小快艇,長度小,速度高,Fr 數大,開爾文船波則以短波長的水波為主,而波前則與速度矢量成較大的夾角。[ 1]
開爾文船波動研究,對於船舶的設計有重要意義,因為船舶的馬力,有一部分消耗在激起船波。利用Fr數與速度成正比,與長度的平方根成反比的規律,可以利用小的模型,縮小船長
M
2
{\displaystyle M^{2}}
倍,同時縮小速度M倍,可以在實驗室中模擬海上舟。[ 2]
Integrand of Kelvin Wake Integral
Kelvin Ship Wake Integrand contour Maple plot
當船隻以速度V駛過深水湖面,波形的幅度在相對於船隻為靜止的極坐標(
ρ
,
ϕ
{\displaystyle \rho ,\phi }
中在船隻的速度矢量方向,
ϕ
=
0
{\displaystyle \phi =0}
),由下列公式表示[ 3]
K
(
ϕ
,
ρ
)
=
∫
−
π
/
2
π
/
2
cos
ρ
cos
(
θ
+
ϕ
)
cos
2
θ
d
θ
{\displaystyle K(\phi ,\rho )=\int _{-\pi /2}^{\pi /2}\cos \rho {\frac {\cos(\theta +\phi )}{\cos ^{2}\theta }}d\theta }
其中
ρ
=
g
r
/
V
2
{\displaystyle \rho =gr/V^{2}}
1
ρ
=
V
2
g
r
{\displaystyle {\frac {1}{\rho }}={\frac {V^{2}}{gr}}}
是福祿數 的平方
F
r
2
{\displaystyle Fr^{2}}
g
{\displaystyle g}
為重力常數
l
{\displaystyle l}
為船的長度。
上列K函數是下列多鞍點積分的正數部分:
K
(
ϕ
,
ρ
)
=
ℜ
(
∫
−
∞
∞
exp
(
i
ρ
f
(
θ
,
ρ
)
d
θ
)
{\displaystyle K(\phi ,\rho )=\Re (\int _{-\infty }^{\infty }\exp(i\rho f(\theta ,\rho )d\theta )}
其中,多鞍點積分的核函數為
f
(
θ
,
ϕ
)
=
−
cos
(
θ
+
ϕ
)
cos
2
θ
{\displaystyle f(\theta ,\phi )=-{\frac {\cos(\theta +\phi )}{\cos ^{2}\theta }}}
此核函數是一個多鞍點函數,振盪劇烈如圖
求其極點,
d
f
(
θ
,
ϕ
)
d
θ
=
sin
(
θ
+
ϕ
)
cos
(
θ
)
2
−
2
cos
(
θ
+
ϕ
)
sin
(
θ
)
cos
(
θ
)
3
=
0
{\displaystyle {\frac {df(\theta ,\phi )}{d\theta }}={\frac {\sin(\theta +\phi )}{\cos(\theta )^{2}}}-{\frac {2\cos(\theta +\phi )\sin(\theta )}{\cos(\theta )^{3}}}=0}
解之,得
θ
1
=
arctan
(
(
1
/
4
)
(
1
+
(
1
−
8
tan
(
ϕ
)
2
)
)
tan
(
ϕ
)
)
=
−
arctan
(
(
1
/
4
)
(
−
1
+
(
1
−
8
tan
(
ϕ
)
2
)
)
tan
(
ϕ
)
)
{\displaystyle \theta _{1}=\arctan({\frac {(1/4)(1+{\sqrt {(1-8\tan(\phi )^{2}))}}}{\tan(\phi )}})=-\arctan({\frac {(1/4)(-1+{\sqrt {(}}1-8\tan(\phi )^{2}))}{\tan(\phi )}})}
由此
ϕ
1
=
19.47
{\displaystyle \phi _{1}=19.47}
度,
ϕ
2
=
−
19.47
{\displaystyle \phi _{2}=-19.47}
度
這就是凱爾文船波的V型波包線 的夾角,最早由凱爾文男爵發現,而且角度與船速無關.[ 4] [ 5] 至於波紋本身則與船速矢量的夾角為
θ
=
π
−
19.47
=
35.3
{\displaystyle \theta =\pi -19.47=35.3}
°[ 1]
Kelvin Wake (Maple density plot)
開爾文船波波形
開爾文船波積分
K
(
ϕ
,
ρ
)
{\displaystyle K(\phi ,\rho )}
必須通過數值積分計算。開爾文男爵根據被積分函數在積分區間內劇烈震盪的特點,提出了駐相法(Method of Stationary Phase)。
原理:當被積分函數劇烈震盪時,除了在極點外,震盪的被積分函數正負相抵消,因此可以將此被積分函數在極點的值作為整個積分的近似,駐相法乃是拉普拉斯方法 的推廣。[ 6]
被積分函數
f
(
θ
,
ϕ
)
=
−
c
o
s
(
θ
+
ϕ
)
c
o
s
2
θ
{\displaystyle f(\theta ,\phi )=-{\frac {cos(\theta +\phi )}{cos^{2}\theta }}}
的兩個極點是:
θ
p
=
a
r
c
t
a
n
(
(
1
/
4
)
∗
(
1
+
(
1
−
8
∗
t
a
n
(
ϕ
)
2
)
)
t
a
n
(
ϕ
)
)
{\displaystyle \theta _{p}=arctan({\frac {(1/4)*(1+{\sqrt {(1-8*tan(\phi )^{2}))}}}{tan(\phi )}})}
θ
m
=
−
a
r
c
t
a
n
(
(
1
/
4
)
∗
(
−
1
+
(
1
−
8
∗
t
a
n
(
ϕ
)
2
)
)
t
a
n
(
ϕ
)
)
{\displaystyle \theta _{m}=-arctan({\frac {(1/4)*(-1+{\sqrt {(}}1-8*tan(\phi )^{2}))}{tan(\phi )}})}
令
f
m
=
f
(
θ
m
,
ϕ
)
=
s
i
n
(
(
1
/
2
)
∗
ϕ
−
(
1
/
2
)
∗
a
r
c
s
i
n
(
3
∗
s
i
n
(
ϕ
)
)
)
s
i
n
(
(
1
/
2
)
∗
ϕ
+
(
1
/
2
)
∗
a
r
c
s
i
n
(
3
∗
s
i
n
(
ϕ
)
)
)
{\displaystyle f_{m}=f(\theta _{m},\phi )={\frac {sin((1/2)*\phi -(1/2)*arcsin(3*sin(\phi )))}{sin((1/2)*\phi +(1/2)*arcsin(3*sin(\phi )))}}}
f
p
=
f
(
θ
p
,
ϕ
)
=
c
o
s
(
(
1
/
2
)
∗
ϕ
+
(
1
/
2
)
∗
a
r
c
s
i
n
(
3
∗
s
i
n
(
ϕ
)
)
)
c
o
s
(
−
(
1
/
2
)
∗
ϕ
+
(
1
/
2
)
∗
a
r
c
s
i
n
(
3
∗
s
i
n
(
ϕ
)
)
)
{\displaystyle f_{p}=f(\theta _{p},\phi )={\frac {cos((1/2)*\phi +(1/2)*arcsin(3*sin(\phi )))}{cos(-(1/2)*\phi +(1/2)*arcsin(3*sin(\phi )))}}}
f
b
a
r
:=
1
/
2
∗
(
f
p
+
f
m
)
{\displaystyle fbar:=1/2*(f_{p}+f_{m})}
D
2
F
=
d
2
F
(
θ
,
ϕ
)
d
θ
2
{\displaystyle D2F={\frac {d^{2}F(\theta ,\phi )}{d\theta ^{2}}}}
D
2
F
p
=
D
2
F
(
θ
p
,
ϕ
)
{\displaystyle D2F_{p}=D2F(\theta _{p},\phi )}
D
2
F
m
=
D
2
F
(
θ
m
,
ϕ
)
{\displaystyle D2F_{m}=D2F(\theta _{m},\phi )}
Δ
:=
(
3
/
4
∗
(
f
m
−
f
p
)
)
(
2
/
3
)
{\displaystyle \Delta :=(3/4*(f_{m}-f_{p}))^{(}2/3)}
u
=
Δ
1
/
2
2
∗
(
1
D
2
F
p
+
1
−
D
2
F
m
)
{\displaystyle u={\sqrt {\frac {\Delta ^{1/2}}{2}}}*({\frac {1}{\sqrt {D2F_{p}}}}+{\frac {1}{\sqrt {-D2F_{m}}}})}
v
=
2
Δ
1
/
2
∗
(
1
D
2
F
p
−
1
−
D
2
F
m
)
{\displaystyle v={\sqrt {\frac {2}{\Delta ^{1/2}}}}*({\frac {1}{\sqrt {D2F_{p}}}}-{\frac {1}{\sqrt {-D2F_{m}}}})}
K
(
ϕ
,
ρ
)
≈
2
∗
π
∗
(
u
∗
c
o
s
(
ρ
∗
f
b
a
r
)
∗
A
i
r
y
A
i
(
−
ρ
(
2
/
3
)
∗
Δ
)
/
ρ
(
1
/
3
)
+
v
∗
s
i
n
(
ρ
∗
f
b
a
r
)
∗
A
i
r
y
A
i
(
1
,
−
ρ
(
2
/
3
)
∗
Δ
)
/
ρ
(
2
/
3
)
)
{\displaystyle K(\phi ,\rho )\approx 2*\pi *(u*cos(\rho *fbar)*AiryAi(-\rho ^{(}2/3)*\Delta )/\rho ^{(}1/3)+v*sin(\rho *fbar)*AiryAi(1,-\rho ^{(}2/3)*\Delta )/\rho ^{(}2/3))}
開爾文船波的波峰,由下列兩個參數方程式描述[ 7]
x
:=
X
∗
s
i
n
(
β
)
∗
(
1
−
(
1
/
2
)
∗
s
i
n
(
β
)
2
)
{\displaystyle x:=X*sin(\beta )*(1-(1/2)*sin(\beta )^{2})}
y
:=
X
∗
s
i
n
(
β
)
2
∗
c
o
s
(
β
)
/
(
2
∗
M
)
{\displaystyle y:=X*sin(\beta )^{2}*cos(\beta )/(2*M)}
Frank J. Oliver, NIST Handbook of Mathematical Functions, 2010, Cambridge University Press
Jame Lighthill Waves in Fluids, Cambridge University Press 1979