波浪理论及其线性化

最新推荐文章于 2023-06-07 10:09:36 发布

Tianan1234

最新推荐文章于 2023-06-07 10:09:36 发布

阅读量3.4k

点赞数 2

分类专栏：海岸动力学

本文链接：https://blog.csdn.net/dengtianan/article/details/115560841

版权

海岸动力学专栏收录该内容

3 篇文章 13 订阅

订阅专栏

文章目录

一、波浪运动的控制方程及其边界条件
二、波浪运动方程的线性化及微幅波理论

一、波浪运动的控制方程及其边界条件

$\begin{cases} G.D.E: &\nabla^2\Phi=\frac {\partial ^2 \Phi} {\partial x^2} + \frac {\partial ^2 \Phi} {\partial z^2}=0, \ \ (-h\le z\le\eta,\ -\infty\le x\le+\infty)\\\\ B.B.C: &\frac {\partial\Phi} {\partial z}=0,\ (z=-h)\\\\ D.F.S.B.C:&\frac {\partial\Phi} {\partial t}+\frac 1 2[(\frac{\partial\Phi} {\partial x})^2+(\frac{\partial\Phi} {\partial z})^2]+gz=0,\ \ (z=\eta)\\\\ K.F.S.B.C:&\frac {\partial\eta} {\partial t}+\frac {\partial\eta} {\partial x}\frac {\partial\Phi} {\partial x}-\frac {\partial\Phi} {\partial z}=0,\ \ (z=\eta)\\\\ L.B.C:&\Phi(x,z,t)=\Phi(x-ct,z) \end{cases}$
【说明】
G.D.E: 波浪运动微分方程（控制方程）
B.B.C: 底部边界条件
D.F.S.B.C: 自由表面动力边界条件
K.F.S.B.C: 自由表面运动边界条件
L.B.C: 侧边界条件

上述原始的波浪运动控制方程在求解的过程中有两个难点：

(1). 两个自由表面边界条件D.F.S.B.C和K.F.S.B.C均有非线性项（分别为 $+\frac 1 2[(\frac{\partial\Phi} {\partial x})^2+(\frac{\partial\Phi} {\partial z})^2]$ 和 $+\frac {\partial\eta} {\partial x}\frac {\partial\Phi} {\partial x}$ ，难以求解；
(2). 自由水面位移 $\eta$ 在边界条件上的值是未知的，即自由表面边界条件不确定。

二、波浪运动方程的线性化及微幅波理论

为了解决上述两个难点，我们将基本方程线性化，实际上是对两个自由表面边界条件D.F.S.B.C和K.F.S.B.C线性化。

【假设条件】：假设波浪运动是缓慢，振幅 $\eta$ 远比波长 $L$ 和水深 $h$ 小，的考虑 $H < < L, H < < h$ ；根据波浪微小振幅假设，只保留一阶量，二阶量可忽略，且任意两个波浪要素之间的乘积也可忽略，于是有：

D.F.S.B.C: $\frac {\partial\Phi} {\partial t}+\cancel{\frac 1 2[(\frac{\partial\Phi} {\partial x})^2 +(\frac{\partial\Phi} {\partial z})^2]}+gz=0,\ \ (z=\eta)\tag{1}$
K.F.S.B.C: $\frac {\partial\eta} {\partial t}+\cancel{\frac {\partial\eta} {\partial x}\frac {\partial\Phi} {\partial x}}-\frac {\partial\Phi} {\partial z}=0,\ \ (z=\eta)\tag{2}$
【数学工具】：
泰勒展开式： $f(x_0+\Delta x,y)=f(x_0,y)+\frac {\partial f(x_0,y)} {\partial x}\Delta x+\frac {\partial^2 f(x_0,y)} {\partial x^2} \frac {(\Delta x)^2} {2!}+...+\frac {\partial^n f(x_0,y)} {\partial x^n} \frac {(\Delta x)^n} {n!}$
只保留一阶： $f(x_0+\Delta x,y)=f(x_0,y)+\frac {\partial f(x_0,y)} {\partial x}\Delta x \tag{3}$
由于原波浪方程中 $\eta$ 的高度是未知不确定的，而在微幅波中我们认为 $\eta$ 足够小，于是可以用泰勒展开式将自由水面（ $z=\eta$ )关联到静水面上（ $z = 0$ )

D.F.S.B.C: $[\frac {\partial\Phi} {\partial t}+gz]_{z=\eta}=0$ 泰勒展开： $[\frac {\partial\Phi} {\partial t}+gz]_{z=0}+\frac {\partial \cancel{[\frac {\partial \Phi} {\partial t}}+gz]} {\partial z}\eta=0$ $[\frac {\partial\Phi} {\partial t}]_{z=0}+\frac {\partial gz} {\partial z}\eta=0$ $[\frac {\partial\Phi} {\partial t}]_{z=0}+g\eta=0$ ${\eta=-\frac 1 g \frac {\partial\Phi} {\partial t}, (z=0)}\tag{4}$
【结论】：在线性波中，自由表面的位移和速度势函数的时间变化率成线性关系

K.F.S.B.C: $[\frac {\partial\eta} {\partial t}-\frac {\partial\Phi} {\partial z}]_{z=\eta}=0$ 泰勒展开： $[\frac {\partial\eta} {\partial t}-\frac {\partial\Phi} {\partial z}]_{z=0}+\frac {\partial (\cancel{\frac {\partial \eta} {\partial t}}-\cancel{\frac {\partial \Phi} {\partial z}})} {\partial z}\eta =0$ $\frac {\partial\eta} {\partial t}-\frac {\partial\Phi} {\partial z}=0, (z=0)$ $\frac {\partial\eta} {\partial t}=\frac {\partial\Phi} {\partial z}, (z=0)\tag{5}$