表面張力波

水面の表面張力波（さざ波）。

ノルウェー、オクスネスのリフィヨルド湾で見られたさざ波。

水と空気の界面に液滴が落ちた衝撃で生まれた表面張力波。

表面張力波（英: capillary wave）とは、流体の相境界上を伝播する波で、ダイナミクスと位相速度が表面張力の効果に支配されるもの。自然界に広く見られ、一般的にさざ波（英: ripple）と呼ばれる。水面の表面張力波の典型的な波長は数センチメートル以下で、位相速度は0.2〜0.3 m/sを超える。

流体界面の波の波長がそれよりも長くなると、表面張力のほか重力と慣性の効果を受ける表面張力重力波となる。一般的に見られる重力波はさらに波長が長くなったものである。

開けた水域で弱い風によって作られるさざ波は英語の海事用語で cat's paw wave と呼ばれ、その微風も cat's paw（猫足風）と呼ばれる。広い海原では、風によって引き起こされた小さいさざ波が成長してはるかに大きな海面波（風浪とうねり）が生じることがある。

分散関係

分散関係とは波の波長と周波数の関係をいう。表面張力の効果に完全に支配される純粋な表面張力波は、重力にも影響される表面張力重力波とは分散関係によって区別できる。

厳密な表面張力波

表面張力波の分散関係は以下となる。

$\omega ^{2}={\frac {\sigma }{\rho +\rho '}}\,|k|^{3}$

$\omega$ は角周波数、$\sigma$ は表面張力、$\rho$ は界面で接する流体のうち重い側の密度、$\rho '$ は軽い側の流体の密度、$k$ は波数を表す。波長は $\lambda ={\frac {2\pi }{k}}$ となる。流体と真空の界面（自由表面）の場合、分散関係は以下のように簡略化される。

$\omega ^{2}={\frac {\sigma }{\rho }}\,|k|^{3}$

表面張力重力波

深水表面で起きる表面張力重力波の分散関係。水面より上の領域は密度ゼロ ($\rho '=0$) としている。位相速度および群速度を $\scriptstyle {\sqrt[{4}]{g\sigma /\rho }}$ で割り、相対波長の逆数 $\scriptstyle {\frac {1}{\lambda }}{\sqrt {\sigma /(\rho g)}}$ の関数としてプロットしたもの。
青線 (A): 位相速度、赤線 (B): 群速度
実線: 表面張力重力波、破線: 重力波、一点鎖線: 表面張力波

一般には波は重力の影響も受けており、表面張力重力波と呼ばれる。無限の深さを持つ二流体の界面で起きる表面張力重力波の分散関係は次のようになる。

$\omega ^{2}=|k|\left({\frac {\rho -\rho '}{\rho +\rho '}}g+{\frac {\sigma }{\rho +\rho '}}k^{2}\right)$

ここで $g$ は重力加速度、 $\rho$ と $\rho '$ は二流体の密度である $(\rho >\rho ')$。第1項の係数 $(\rho -\rho ')/(\rho +\rho ')$ はアトウッド数である。

重力波領域

詳細は「エアリー波理論」を参照

波長が長い、すなわち波数 $k=2\pi /\lambda$ が小さい場合には、表面張力重力波の分散関係における第1項が支配的となり重力波に帰着する。この極限で波の群速度は位相速度の半分となる。このとき波束（群速度で伝播する）に含まれる波の山の一つ（位相速度で伝播する）に注目すると、その山は波束の背後から近づきつつ成長し、波束の腹を通り過ぎると減衰しながら前方に消えていく。

表面張力波領域

表面張力波の分散。

波長 $\lambda$ が短い、すなわち波数 $k$ が大きい波は表面張力波であり、前節と逆の振る舞いを示す。波の山は波束の前方で現れ、高さを増しながら波束の中心に近づき、波束の背後に消えていく。

最小位相速度

これら2つの極限の間には重力による分散が表面張力による分散を相殺する点がある。その特定の波長では群速度が位相速度と等しくなり、分散は生じない。それと正確に同じ波長において表面張力重力波の位相速度は最小値を取る。この臨界波長 $\lambda _{m}$ よりはるかに短い波長の波では表面張力が、はるかに長い波長の波では重力が支配的となる。$\lambda _{m}$ とそこから導かれる最小位相速度 $c_{m}$ は以下で与えられる。

${\begin{aligned}\lambda _{m}&=2\pi {\sqrt {\frac {\sigma }{(\rho -\rho ')g}}}\\c_{m}&={\sqrt {\frac {2{\sqrt {(\rho -\rho ')g\sigma }}}{\rho +\rho '}}}\end{aligned}}$

空気と水の界面では $\lambda _{m}=1.7$ cm、$c_{m}=0.23$ m/sとなる。

液体に小石か滴を落とすと様々な波長の波が同心円状に広がっていくが、それらが伝播するのはゆっくり広がる円の外側のみで、円の内側では流体は静止する。この円は最小群速度に対応する焦線である。

導出

リチャード・ファインマンの言によると「誰もが容易に目にすることができ、初等コースで波の例としてよく持ち出される[水波]は … 考えられる限り最悪の例であり、… 波が持ちうるあらゆる困難さを備えている」。実際、一般的な分散関係の導出は非常に複雑である。

系のエネルギーには重力、表面張力、流体運動の三つが寄与する。最初の二つはポテンシャルエネルギーであり、前掲の分散関係における括弧内の二項は（ $g$ と $\sigma$ を含むことから分かるように）これらに起因する。重力の効果をモデル化する際には、流体の密度が一定（すなわち非圧縮性）であり、$g$ も一定（波は重力が大きく変化するほどの高さに至らない）と仮定されている。表面張力に関しては、水平面を基準とした水面の鉛直変位（表面の導関数で表される）が小さいとされている。通常の波ではどちらも十分に良い近似となる。

三つ目の寄与は流体の運動エネルギーから来ている。三つのうちでは最も複雑であり、流体動力学的な枠組みが必要となる。ここでも非圧縮性（波の速度が流体中の音速よりはるかに小さいときにあてはまる）と、さらに渦なし流れが仮定される。それにより流れはポテンシャル流れとなる。これらも一般的な状況を概して良く近似する。そうして得られるポテンシャル方程式（ラプラス方程式となる）は適切な境界条件のもとで解くことができる。まず、水面から十分に遠方で流速は消失しなければならない（ここで想定している「深水」の状況がそれにあたる。そうでなければ結果は複雑になる。en:wind wave（風浪）を参照）。さらに流速の垂直成分は表面の運動と一致している必要がある。

最終的に、分散関係に対する運動エネルギーの寄与は括弧外の $|k|$ に現れる。この係数により、$k$ が低いときから高いときまですべての領域で分散性が生じる（例外は二つの分散性が相殺される $|k|$ の値とその近辺である）。

二つの半無限な流体領域の界面に発生する表面張力重力波の分散関係

二つの流体領域があり、それらの界面に表面張力が働くとする。界面は時間平均すると水平面をなす。二流体の密度は異なっており、下側と上側の密度をそれぞれ $\rho$ および $\rho '$ とする。流体は非粘性かつ非圧縮性であり、流れは渦なしだと仮定する。このような流れはポテンシャル流であり、下側と上側の流速はそれぞれ $\nabla \Phi$ および $\nabla \Phi '$ で与えられる。$\Phi (x,y,z,t)$ と $\Phi '(x,y,z,t)$ は速度ポテンシャル である。 エネルギーには重力のポテンシャル $V_{\mathrm {g} }$、表面張力のポテンシャル $V_{\mathrm {st} }$、運動エネルギー $T$ の三つの寄与がある。重力の項 $V_{\mathrm {g} }$ はもっとも単純であり、重力のポテンシャル密度 ($\rho gz,\rho 'gz$) を基準点から界面の鉛直座標 $z=\eta (x,y,t)$ まで積分することで ${\begin{aligned}V_{\mathrm {g} }&=\iint dx\,dy\;\int _{0}^{\eta }dz\;(\rho -\rho ')gz\\&={\frac {1}{2}}(\rho -\rho ')g\iint dx\,dy\;\eta ^{2}\end{aligned}}$ を得る。ただし界面の平均高さを $z=0$ とした。 変位 $\eta (x,y,t)$ によって界面の面積が増えると、表面張力エネルギーはそれに比例して増加する。 ${\begin{aligned}V_{\mathrm {st} }&=\sigma \iint dx\,dy\;\left[{\sqrt {1+\left({\frac {\partial \eta }{\partial x}}\right)^{2}+\left({\frac {\partial \eta }{\partial y}}\right)^{2}}}-1\right]\\&\approx {\frac {1}{2}}\sigma \iint dx\,dy\;\left[\left({\frac {\partial \eta }{\partial x}}\right)^{2}+\left({\frac {\partial \eta }{\partial y}}\right)^{2}\right]\end{aligned}}$ 上の最初の等式ではモンジュによる表現を用いた面積の計算が行われている。第二の等式は $\eta$ の導関数が小さいとき（界面があまり波打っていないとき）に成立する。 最後に流体の運動エネルギーからの寄与は以下で与えられる。 $T={\frac {1}{2}}\iint dx\,dy\;\left[\int _{-\infty }^{\eta }dz\;\rho \,\left|\mathbf {\nabla } \Phi \right|^{2}+\int _{\eta }^{+\infty }dz\;\rho '\,\left|\mathbf {\nabla } \Phi '\right|^{2}\right]$ ここで流体が非圧縮性であり、流れが渦なしであること（多くの場合、妥当な仮定である）を用いる。その結果 $\Phi (x,y,z,t)$ と $\Phi '(x,y,z,t)$ はいずれもラプラス方程式 $\nabla ^{2}\Phi =0$, $\nabla ^{2}\Phi '=0$ に従う。 これらを解くために適切な境界条件を与える。すなわち、界面から十分に遠方では $\Phi$ と $\Phi '$ はいずれも消失しなければならない（ここで想定されている「深水」の状況が当てはまる）。 グリーンの恒等式を用い、さらに界面の鉛直方向変位が小さい（そのため $z=\eta$ までの積分を $z=0$ までで近似することができる）と仮定すると、運動エネルギーは以下のように表せる。 $T\approx {\frac {1}{2}}\iint dx\,dy\;\left[\rho \,\Phi \,{\frac {\partial \Phi }{\partial z}}\;-\;\rho '\,\Phi '\,{\frac {\partial \Phi '}{\partial z}}\right]_{z=0}$ 分散関係を得るには、界面を $x$ 方向に伝播する正弦波 $\eta =a\,\cos \,(kx-\omega t)=a\,\cos \,\theta$ を考えれば十分である。振幅を $a$、波の位相を $\theta =kx-\omega t$ とした。速度ポテンシャルを界面の運動と結び付ける運動学的境界条件として、界面において両方の流体の鉛直速度成分は波の運動と一致しなければならない。 ${\frac {\partial \Phi }{\partial z}}={\frac {\partial \Phi '}{\partial z}}={\frac {\partial \eta }{\partial t}}$   ($z=0$) 各領域の速度ポテンシャルを求めるにあたって変数分離を試みると、それぞれのポテンシャル場は以下のように書かれる。 ${\begin{aligned}\Phi (x,y,z,t)&=+{\frac {1}{|k|}}{\text{e}}^{+|k|z}\,\omega a\,\sin \,\theta \\\Phi '(x,y,z,t)&=-{\frac {1}{|k|}}{\text{e}}^{-|k|z}\,\omega a\,\sin \,\theta \end{aligned}}$ 以上より、波のエネルギーに対する三つの寄与を水平面内で $x$ 方向に一波長分、$y$ 方向に単位幅にわたって積分すると以下のようになる。 ${\begin{aligned}V_{\mathrm {g} }&={\frac {1}{4}}(\rho -\rho ')ga^{2}\lambda \\V_{\mathrm {st} }&={\frac {1}{4}}\sigma k^{2}a^{2}\lambda \\T&={\frac {1}{4}}(\rho +\rho '){\frac {\omega ^{2}}{|k|}}a^{2}\lambda \end{aligned}}$ 分散関係は以下のラグランジアン $L=T-V$から求められる (ここで $V=V_{\mathrm {g} }+V_{\mathrm {st} }$) 。 $L={\frac {1}{4}}\left[(\rho +\rho '){\frac {\omega ^{2}}{|k|}}-(\rho -\rho ')g-\sigma k^{2}\right]a^{2}\lambda$ 線形波動理論のもとで正弦波の平均ラグランジアンは常に $L=D(\omega ,k)a^{2}$ の形を取る。したがって、唯一の自由なパラメータである $a$ についての変分条件から分散関係 $D(\omega ,k)=0$ が導かれる。ここで $D(\omega ,k)$ は上式の角かっこ内にあたり、分散関係は $\omega ^{2}=|k|\left({\frac {\rho -\rho '}{\rho +\rho '}}\,g+{\frac {\sigma }{\rho +\rho '}}\,k^{2}\right)$ となって前掲式と一致する。 結果として、水平面の単位面積当たり波の平均エネルギー $(T+V)/\lambda$ は ${\bar {E}}={\frac {1}{2}}\,\left[(\rho -\rho ')\,g+\sigma k^{2}\right]\,a^{2}$ である。また、線形波で一般的なようにポテンシャルと運動エネルギーは等しい（エネルギー等分配の法則は保たれている）。

ギャラリー

• 
  アメンボによって作られたさざ波。
• 
  そよ風によって湖水表面に作られたさざ波。

関連項目

• 毛細管現象
• 分散 (水波)
• 熱表面張力波
• 二層流
• ウェーブリップル

参考文献

• Longuet-Higgins,M. S. (1963). “The generation of capillary waves by steep gravity waves”. Journal of Fluid Mechanics 16 (1): 138–159. Bibcode: 1963JFM....16..138L. doi:10.1017/S0022112063000641. ISSN 1469-7645. 
• Lamb, H. (1994). Hydrodynamics (6th ed.). Cambridge University Press. ISBN 978-0-521-45868-9 
• Phillips, O. M. (1977). The dynamics of the upper ocean (2nd ed.). Cambridge University Press. ISBN 0-521-29801-6 
• Dingemans, M. W. (1997). Water wave propagation over uneven bottoms. Advanced Series on Ocean Engineering. 13. World Scientific, Singapore. pp. 2 Parts, 967 pages. ISBN 981-02-0427-2 
• Safran, Samuel (1994). Statistical thermodynamics of surfaces, interfaces, and membranes. Addison-Wesley 
• Tufillaro, N. B.; Ramshankar, R.; Gollub, J. P. (1989). “Order-disorder transition in capillary ripples”. Physical Review Letters 62 (4): 422–425. Bibcode: 1989PhRvL..62..422T. doi:10.1103/PhysRevLett.62.422. PMID 10040229. http://scholarship.haverford.edu/cgi/viewcontent.cgi?article=1062&context=physics_facpubs. 

外部リンク

• sklogwikiの表面張力波エントリ