[Ocean Modelling for Begineers] Ch5. 2D Shallow-Water Modelling

本章利用二维浅水模型研究表面重力波的不同物理过程，如湖水中风驱动流体，正压不稳定机制（？the barotropic instability mechanism）。本章将为读者介绍使用不同的对流格式模拟物体运动轨迹和动量方程中非线性项。

5.1 浅水中长波

5.1.1 二维浅水波动方程

做以下假设：

1. 水体中密度为常数
2. 忽略底部摩擦
3. 忽略科氏力（周期与惯性周期相比非常小）
4. 忽略非线性项（波速远大于流体微团速度）

控制方程简化为

\[\begin{eqnarray} \begin{aligned} & \frac{\partial u}{\partial t} = -g \frac{\partial \eta}{\partial x} \cr & \frac{\partial v}{\partial t} = -g \frac{\partial \eta}{\partial y} \cr & \frac{\partial \eta}{\partial t} = -\frac{\partial (uh)}{\partial x} - \frac{\partial (vh)}{\partial y} \end{aligned} \end{eqnarray}\]

5.1.2 Arakawa C网格

Arakawa C网格（Arakawa and Lamb, 1977）是一种交错网格形式，其中流速点分布在相邻水位点周围（如图5.1）。这种网格形式经常被海洋学模型所采用，也被后面模型网格形式。注意流速节点u，v并不在同一节点位置。

5.1.3 有限差分方程

二维浅水方程离散格式，其中j为y方向单元序号，k为x方向单元序号。

\[\begin{equation} \begin{aligned} & u^{n+1}\_{j,k} = u^{n}\_{j,k} - \Delta t g(\eta^n_{j, k+1} - \eta^n_{j, k})/\Delta x \cr & v^{n+1}\_{j,k} = v^{n}\_{j,k} - \Delta t g(\eta^n_{j+1, k} - \eta^n_{j, k})/\Delta y \cr & \eta^*\_{j, k} = \eta^n_{j, k} - \Delta t \\{ (u^{n+1}\_{j,k} h_e - u^{n+1}\_{j,k-1} h\_w)/\Delta x - (v^{n+1}\_{j,k} h_n - v^{n+1}\_{j-1,k} h\_s)/\Delta y \\} \end{aligned} \end{equation}\]

二维格式中一阶Shapiro过滤算子

\[\begin{equation} \eta_{j,k}^{n+1} = (1-\epsilon)\eta_{j,k}^\* + 0.25 \epsilon(\eta_{j,k-1}^\* + \eta_{j,k+1}^\* + \eta_{j-1,k}^\* + \eta_{j+1,k}^\*) \end{equation}\]

5.1.4 陆地与海岸线

陆地定义为水深为0或者为负值，并且在模拟过程中流速分量恒为0。海岸线处隐式的设定法线方向流量为0。由于Arakawa-C网格的优势，这种方法可以按照以下两种情形设置，如果陆地在相邻节点东侧，那么该点u分量为0；同样的，如果陆地在相邻节点北侧，那么该点v分量为0。图5.2给出了在Arakawa-C网格中陆地和海岸线分布形式。干湿算法也可以类似按照一维情形处理。

5.1.5 稳定性

\[\begin{equation} \Delta t \le \frac{min(\Delta t, \Delta y)}{ \sqrt{2gh_{max}} } \end{equation}\]

5.2 练习8：浅水中长波

5.2.1 目标

本练习目的是模拟环形表面重力波在二维平面上传播过程。

5.2.2 任务描述

计算区域为 \(500m \times 500m\) 矩形网格，空间步长为10m，全场水深为10m，四周为陆边界。

初始时刻在中心点有1m的初始水位，随着时间递进开始向四周传播形成环形波。时间步长为0.1s，共模拟100s

5.2.3 代码范例