Transcript 微小振幅波の理論
微小振幅波の理論
酒井哲郎:海岸工学入門,森北出版
第11章(pp.128-133)
波峰(wave crest):波形の最高水位
波谷(wave trough):波形の最低水位
波高(wave height)H:波峰と波谷間の鉛直距離
波長(wave length)L:相続く波峰間(波谷間の)水平距離
周期(wave period)T:海面上のある固定点を相続く波峰(波谷)が通過
に要する時間
波速(celerity)C:波形の進行速度(位相速度)、C=L/T
水深(wave depth)h:水底面と平均水面間の鉛直距離
波峰
波長L
z
0
x
平均水面(静水面)
z=0
波高H
水深h
波谷
w
波の進行方向
u
水底面 z=-h
波の理論的取り扱い
波は水面が変形する.
波の振幅が小さい場合
微少振幅波(small amplitude wave)
微少振幅波理論
線形波(エアリー波)
波の振幅が大きい場合
有限振幅波(finite amplitude wave)
有限振幅波理論
ストークス波
クノイド波
ソリトン(孤立波)
流体運動の支配方程式
・x-z平面で考える.
・非圧縮性流体
・完全流体
連続の式
X方向運動方程式
(オイラーの式)
u w
0
x z
u
u
u
1 p
u w
t
x
z
x
流体運動の非線形性
z方向運動方程式
(オイラーの式)
w w
w
1 p
u w g
t
x
z
z
流体運動の非線形性
上記方程式系を適切な境界条件のもとで解けばよいが
波動場解析では,次のようにして解析を行う.
波動場は非回転(渦度がない)流れとして取り扱い可能
速度ポテンシャルFが存在する.
F
u
x
F
w
z
渦度に代入してみる.
w u F F 2F 2F
y
0
x z x z z x xz zx
連続の式に代入してみる.
u w F F 2F 2F
2 2 0
x z x x z z x
z
ラプラスの式
波動場ではラプラスの式を適切な境界条件のもとで解くことになる.
水表面の境界条件
水表面は形状が自由に変化する(自由水表面とよばれるゆえん).
z
h
x
静水面から水面までの距離をh(x,t)とする.hは場所xと時間tの関数である.
水表面の境界条件としては,水表面上で圧力が0であることを用いるが,水表
面の位置hが未知量である.
したがって,まず水表面がどこにあるのかを表す必要がある.
運動学的境界条件(kinematic boundary condition)
水面上の流体粒子はつねに水面上に存在する
h
h
uh
wh w( x,h, t )
t
x
t+dt
z
uhd t
h(x, t+dt)
h(x, t)
h
x
h
dt
t
whd t
uhd t
t
uh:水表面におけるx方向流速
wh:水表面におけるz方向流速
x
h
h
uh
wh
t
x
h
h F F h
wh uh
t
x z x x
: z h
初期に水表面に存在した流体粒子は時間が経過しても水表面に存在する.
x,y二次元平面では
h ( x, y, t )
h ( x, y, t )
h ( x, y, t )
uh ( x, y,h, t )
vh ( x, y,h, t )
wh ( x, y,h, t )
t
x
y
h F F h F h
t z x x y y
: z h
力学的境界条件(dynamic boundary condition)
流速は,ラプラスの式を解いて速度ポテンシャルが求まればそれを微分する
ことで求まる.それでは圧力はどうすれば求まるのか?
水理学Ⅰで学習した拡張されたベルヌイの定理(圧力方程式)を用いる.
F u2 v2 w2
P F (t )
t
2
p dp
P
:圧力関数
:質量力ポテンシャル
F u2 w2 p
gz F (t )
t
2
2
2
F 1 F F p
gz F (t )
t 2 x z
水表面(z=h)では圧力0(ゲージ圧)なので
2
2
F 1 F F
gh F (t )
t 2 x z
ここで F F
F (t)dt
t
0
: z h
と書き改めると上式は
2
2
F 1 F F
gh 0
t 2 x z
: z h
新たに書き直したFはラプラスの式を満足する.
x,y二次元平面では
2
2
2
F 1 F F F
gh 0
t 2 x y z
その他の境界条件
底面あるいは固体壁面:
底面あるいは固体壁面に対して垂直bな方向の流速成分は0である.
F
u
0
x
w
F
0
z
ラプラスの式
自由水面にお
ける境界条件
2 F 2F
2 0
2
x
z
h F F h
t z x x
: z h
流体運動の非線形性
2
2
F 1 F
F
gh 0
t 2 x y
流体運動の非線形性
圧力方程式
2
2
F 1 F F p
gz 0
t 2 x y
流体運動の非線形性
: z h
微少振幅波理論
前ページの水表面の境界条件はz=h(水表面)で与えるべきもので
あるが, hそのものが未知量なので取り扱いは複雑・困難である.
h ( x, t ) F ( x,h, t ) F ( x,h, t ) h ( x, t )
t
z
x
x
上記境界条件をz=0の回りでテーラー展開する.
h ( x, t ) F ( x,0, t ) F ( x,0, t ) h ( x, t )
t
z
x
x
F ( x,0, t )
F ( x,0, t ) h ( x, t )
h
h
z
z
z
x
x
この式はz=0(静水面)で成立する式である.
hを微少と仮定する.水面近傍の流速も微少となる.
h ( x, t ) F ( x,0, t ) F ( x,0, t ) h ( x, t )
t
z
x
x
F ( x,0, t )
F ( x,0, t ) h ( x, t )
h
h
z
z
z
x
x
運動学的条件は静水面(z=0)で以下のように与えられる.
h ( x, t ) F ( x,0, t )
t
z
:z 0
微少振幅波理論による運動学的境界条件
2
2
F ( x,h, t ) 1 F ( x,h, t ) F ( x,h, t )
gh ( x, t ) 0
t
2
x
y
上記境界条件をz=0の回りでテーラー展開する.
2
2
F ( x,0, t ) 1 F ( x,0, t ) F ( x,0, t )
gh ( x, t )
t
2
x
y
2
2
F
x
,0,
t
F
x
,0,
t
F
x
,0,
t
(
)
(
)
(
)
1
h
h
z
t
2 z
x
y
z
0
hを微少と仮定する.水面近傍の流速も微少となる.流速の積はさら
に微少になる.
2
2
F ( x,0, t ) 1 F ( x,0, t ) F ( x,0, t )
gh ( x, t )
t
2
x
y
2
2
F ( x,0, t )
1 F ( x,0, t )
F ( x,0, t )
h
h
h
z
t
2 z
x
z
y
0
力学的境界条件は静水面(z=0)で以下のように与えられる.
1 F ( x,0, t )
h ( x, t )
g
t
:z 0
微少振幅波理論による力学的境界条件
微少振幅波理論による波動場の支配方程式
ラプラスの式
自由水面にお
ける境界条件
2 F 2F
2 0
2
x
z
h F
t z
h
圧力方程式
1 F
g t
:z 0
:z 0
z=hからz=0になっ
たことに注意.
静水面に境界条件
を与えるので,取り
扱いが容易になる.
F p
gz 0
t
流体運動の非線形性はなくなり,すべて線形の微分方程式となった.
この方程式系から誘導される波は線形波と呼ばれる.
H
h cos ( kx t )
2
2
2
:波数
k
L
T
:角周波数
波動の基本解として上記の余弦波を考える.ラプラスの式および境界
条件は線形なので,変数分離法を用いて速度ポテンシャルを求めるこ
とができる.
F F ( z)sin ( kx t )
速度ポテンシャルをこのようにおいてラプラスの式に代入する.
h
H
w
sin ( kx t )
t
2
類似の形
F dF ( z)
w( z)
sin ( kx t )
z
dz
dF ( z)
w(0)
sin ( kx t ) : z 0
dz
ラプラスの式から次の常微分方程式が得られる.これを与えられた境
界条件の下で解く.
d 2 F ( z) 2
k F ( z) 0
2
dz
一般解
F ( z) Aexp(kz) B exp(kz)
境界条件を満足するように係数A,Bを決める.
最終的に次式が得られる.
Hg cosh k ( z h)
F
sin ( kx t )
2 cosh kh
Hg cosh k ( z h)
F
sin ( kx t )
2 cosh kh
H
h cos ( kx t )
2
1 F
h
g t
これらを使用して重要な関係式を以下のように得ることができる.
gk tanh kh
2
これは角周波数と波数(周期と波長)の関係を表すものである.
水面の波は周期と波長はそれぞれ好き勝手にはとれないことを意味する.
上式を分散関係式(dispersion relation equation)という.
速度ポテンシャルは分散関係式を用いて次のように書き換えることができる.
Hg cosh k ( z h)
F
sin ( kx t )
2 cosh kh
H cosh k ( z h)
F
sin ( kx t )
2k
sinh kh
波速
L
C
T k
分散関係式を変形する.
2 gk tanh kh
gT
2 h gT
2 h
C tanh kh
tanh
tanh
k
2
L
2
CT
kC
g
を代入して
g
gL
2 h
C tanh kh
tanh
k
2
L
2
gL
2 h
C
tanh
2
L
波長が決まれば波速も決まる.
浅海波(h/L<1/2)
gT 2
2 h
L
tanh
2
L
gT
2 h
gL
2 h
C
tanh
tanh
2
CT
2
L
深海波(h/L≧1/2)
2 h
2 h
深海波では tanh
tanh
1.0
CT
gT
C
2
極浅海波(長波)(h/L<1/25)
2 h
長波では tanh
L
C gh
となる.
L
gT 2
L
2
2 h となる.
L
L T gh
双曲線関数
e e
sinh x
2
x
x
e e
cosh x
2
x
x
e e
tanh x x x
e e
x
x