船隻尾波的鳥瞰圖 從前方觀測的船隻尾波 尾波 （英語 ：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}}
[ 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 )}
原理：當被積分函數劇烈震盪時，除了在極點外，震盪的被積分函數正負相抵消，因此可以將此被積分函數在極點的值作為整個積分的近似，駐相法乃是拉普拉斯方法 的推廣。[ 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 
. doi:10.1023/B:INAM.0000041400.70961.1b.
