表面重力波推导

表面重力波推导
0 概念定义
$0.1$ 速度势
当流体作无旋运动时, 运动方程和连续方程将具有非常简单的形式. 首先, 速度旋度 为零 $(\nabla \times \mathbf{V}=0)$ 的条件保证在任何单连通域内存在单值速度势 $\phi (x,y,z,t)$; 使得
$$\mathbf{V}=\nabla \phi (x,y,z,t),$$
或
$$u=\frac{\partial \phi }{\partial x},v=\frac{\partial \phi }{\partial y},w=\frac{\partial \phi }{\partial z},$$
1 无限深度水域

1. 1 边界条件
   （1）水域底部的运动学边界条件
   $${\frac{\partial \phi }{\partial n}|}_{x=-d}=0$$
   （2）自由表面运动学边界条件
   $${\frac{\partial \phi }{\partial z}|}_{z=\zeta }=\frac{\partial \zeta }{\partial t}+{\frac{\partial \zeta }{\partial x}\frac{\partial \phi }{\partial x}|}_{z=\zeta }+{\frac{\partial \zeta }{\partial y}\frac{\partial \phi }{\partial y}|}_{z=\zeta },$$
   但在小振幅假定之下, 上式中的乘积项相对于其它两项可以忽略不计, 而且在同样的近似程 度下可认为
   $${{\frac{\partial \phi }{\partial z}|}_{z=\zeta }\approx \frac{\partial \phi }{\partial z}|}_{z=0},$$
   于是我们便得到了线性的自由表面运动学边界条件
   $$\frac{\partial \zeta }{\partial t}={\frac{\partial \phi }{\partial z}|}_{z=0},$$
   （3）自由表面动力学边界条件
   $${\frac{\partial \phi }{\partial t}|}_{z=\zeta }+g\zeta +{\frac{1}{2}(\nabla \phi \cdot \nabla \phi )|}_{z=\zeta }+\frac{P_s(x,y,t)}{\rho }=0$$
   可以略去乘积项 $(\nabla \phi \cdot \nabla \phi )$, 如果给定自由表面上的压强分布 $P_S(x,y,t)$, 我们便得到自 由表面上的动力学边界条件 为
   $$\frac{P_s(x,y,t)}{\rho }=-{\frac{\partial \phi }{\partial t}|}_{z=\zeta }-g\zeta$$
   考虑到小振幅的假定, 上式可改写为
   $$\frac{P_s(x,y,t)}{\rho }=-{\frac{\partial \phi }{\partial t}|}_{z=0}-g\zeta .$$
   考虑 $P_S(x,y,t)=$ const 的情形, 可引进新的势函数
   $${\phi }^{\prime }=\phi +\frac{P_0}{\rho }t,$$
   式中 $P_0$ 为自由表面党压强的数值, 自由表面动力学边界条件变为
   $${\frac{\partial {\phi }^{\prime }}{\partial t}|}_{s=0}+g\zeta =0.$$
   为简便计, 略去上两式中的搬号
   $${\frac{\partial \phi }{\partial t}|}_{z=0}+g\zeta =0$$
   在常表面压强的条件下,自由表面的运动学边界条件和动力学边界条件可以合并为
   $${\left(\frac{\partial \phi }{\partial z}+\frac{1}{g}\frac{{\partial }^2\phi }{\partial t^2}\right)|}_{z=0}=0$$
2. 2 求解方程
   因为任一点处的水面都以固定的圆频率 (例如, 记为 $\omega$ ) 作振动, 所以可以认为波动 场中的有关物理量也具有这一特性。特别, 任一点处的速度势对时间应该呈现简谐形式的 变化。因此波动势函数必具有如下的形式。
   $$\begin{array}{c}\phi =\cos \cot \cdot f(x,z)+\sin \omega t\cdot g(x,z),\\ \phi =A(x,z)\cos (\cot +\epsilon (x,z)),\end{array}$$
   上式中 ${\phi }_1(x,z)$ 为复值函数, $\phi$ 完全由下列方程和边界条件所确定
   $$\{\begin{array}{l}\Delta \phi =0,z<0,\\ {\left(\frac{\partial \phi }{\partial z}+\frac{1}{g}\frac{{\partial }^2\phi }{\partial t^2}\right)|}_{z=0}=0,\\ \phi \text{ 及其导数有界. }\end{array}$$
   对于无限深水域来说, 式 $\frac{\partial \phi }{\partial z}+\frac{1}{g}\frac{{\partial }^2\phi }{\partial t^2}$ 不仅当 $z=0$ 时等于零, 而且在整个水域它都为 零, 因此边界条件可以改写为方程的形式, 即
   $$\frac{\partial \phi }{\partial z}+\frac{1}{g}\frac{{\partial }^2\phi }{\partial t^2}=0,z<0$$
   则有方程和边界条件
   $$\{\begin{array}{l}\frac{{\partial }^2{\phi }_1}{\partial x^2}+\frac{{\partial }^2{\phi }_1}{\partial z^2}=0,z<0,\\ \frac{\partial {\phi }_1}{\partial z}-\frac{{\omega }^2}{g}{\phi }_1=0,z<0,\\ {\phi }_1\text{ 及 }{\phi }_1\text{ 的导数有界. }\end{array}$$
   可求得
   $$\begin{array}{c}\phi (x,z,t)=\cos \omega t\cdot e^{kz}\left(A_1\cos kx+B_1\sin kx\right)\\ +\sin \omega t\cdot e^{kz}\left(A_2\cos kx+B_2\sin kx\right)\\ {\omega }^2=gk\end{array}\}$$
   2 常值深度水域
3. 1 方程和边界条件
   此时 $\phi (x,z,t)$ 应满足如下的方程和边界条件:
   $$\{\begin{array}{l}\Delta \phi =0,z<0,\\ {\left(\frac{\partial \phi }{\partial z}+\frac{1}{g}\frac{{\partial }^2\phi }{\partial t^2}\right)|}_{z=0}=0,\\ {\frac{\partial \phi }{\partial z}|}_{z=-d}=0.\end{array}$$
4. 2 求解方程
   当任意点的水面都以固定频率 $\omega$ 随时间作简谐形式的起伏时, 波动势函数仍然具有 下式实部或虚部的形式
   $$\phi (x,z,t)=e^{i\omega t}{\phi }_1(x,z)$$
   代入上面的三个关系式后, 可得到决定 ${\phi }_1(x,z)$ 的相应方程和边界条件
   $$\begin{array}{rl}& \frac{{\partial }^2{\phi }_1}{\partial x^2}+\frac{{\partial }^2{\phi }_1}{\partial z^2}=0,z<0,\\ & {\left(\frac{\partial {\phi }_1}{\partial z}-\frac{{\omega }^2}{g}{\phi }_1\right)|}_{z=0}=0,\\ & {\frac{\partial {\phi }_1}{\partial z}|}_{z=-d}=0.\end{array}$$
   参照深水情形, 可取 ${\phi }_1(x,z)$ 的试解形式为
   $${\phi }_1(x,z)=Z(z)(A\cos kx+B\sin kx),$$
   式中 $A$ 和 $B$ 为任意的复常数, 于求得
   $$\begin{array}{l}{\phi }_1(x,z)=\mathrm{ch}k(z+d)(A\cos kx+B\sin kx)\\ {\omega }^2=gk\text{ th }kd.\end{array}\}$$
   取实部或虚部后，得到最普遍形式的势函数解
   $$\begin{array}{l}\phi (x,z,t)=\mathrm{ch}k(z+d)\cos \omega t\left(A_1\cos kx+B_1\sin kx\right)\\ +\mathrm{ch}k(z+d)\sin \cot \left(A_2\cos kx+B_2\sin kx\right)\\ {\omega }^2=gk\text{ th }kd\end{array}\}$$
   3 讨论
   从常深度水域中波的频散关系可得相速度关系
   $$c^2=\frac{g}{k}\text{ th }kd=\frac{gL}{2\pi }\text{ th }kd$$
   当水深 $d$ 无限时, th $kd=1$, 因此, 波速随波长得增加而增加, 有
   $$c^2=\frac{gL}{2\pi }$$
   当水深相对于波长很小时, th $kd\approx kd$, 因此, 波速随水深的增加而增加, 有
   $$c^2=gd$$