|
1.控制方程 海洋与大气的控制方程是由一组偏微分方程表示的。以大气的基本控制方程为例，其中(1)(2)(3)式为运动方程，(4)式为热力学方程，(5)式为水汽平流方程，(6)式为连续方程，可由质量守恒推导而来。具体项的物理含义和数学推导不再赘述。倘若能精确求解如下方程组，即可对天气进行预测或进行数值实验。但该方程组的复杂程度目前难以求得其解析解，因此需要借助数值方法求解下列方程组。 上述方程组在进行求解时，考虑到计算开支和应用场景，通常会对模型进行简化。如果在水平尺度远大于垂直尺度时，例如对海岸模拟时，可将控制方程近似得到浅水方程组(Shallow Water Equations)。该方程组的形式相对简单，其解可表示惯性重力波、平流和罗斯贝波等现象，因此该方程组也被称为浅水波方程。因其形式简单，并在河道、河口和近海有较多应用空间，该方程组十分适合用于测试和检验数值方法。在不考虑粘性、风应力等，最简单的浅水方程组的形式如下所示。 2.离散化与差分离散介绍完了方程，下面就该到了求解。学过编程语言的同学们肯定知道，计算机只能处理离散过程，我们可以很简单的通过任何一种编程语言去计算 ，然而我们却无法让计算机直接计算 。因此，数值方法的关键就是如何计算导数。在学习微积分的时候，导数可以这样表示 而有限差分法的核心就是使用差分近似微分，当差分的间距取得足够小时，就可以近似计算出导数的值。此时，我们再去计算上式 时，就可以近似计算 。这样就完成了一个简单的离散化，并用差分近似微分进行了计算求解。而应用在偏微分方程的求解时，还需要探讨几个问题。 以平流方程为例，该方程形式如下，是一种双曲型偏微分方程。和我们要求解的流体力学方程形式相似。 不妨取 这样的离散点，则差分方程的解 表示为 ，即上标代表时间各点，下标代表空间各点。则平流方程的差分方程可表示成如下形式。 在对方程求解前，必须探讨两个问题，当时间步长和空间步长趋近于无穷小时，差分方程是否能收敛于偏微分方程？差分方程的解是否能收敛于偏微分方程的解么？前者是在探讨差分方程的一致性，后者是在探讨差分方程的收敛性。 不妨将这些项进行泰勒展开，然后代回差分方程，可得下式 不难发现，差分方程比微分方程多了一些项，而这些项都是随着步长趋于0而趋于0的，因此可以很好的证明一致性。而多出的这些项被称之为截断误差，是差分近似微分时造成的误差。因为我们知道，做泰勒展开时，只有保留所有项才能完全相等，然而展开项是无穷多的。在实际做运算时，针对特定的精度，只能取到相应阶数。而被舍弃掉的高阶小项被称之为截断误差。
, P* H3 A6 i" t' }2.2 收敛性上式差分方程可以化成如下形式，其中 。不难看出，化成这种形式后，可以发现这种形式的差分，本质上是使用n时刻的值去插值出n+1时刻的值。当 时，就等价于上述我们所说的差分格式。  0 k8 T3 s/ B% ^7 M9 k# T8 `9 R7 F8 {$ [
登录/注册后可看大图 也被称为Courant数，其绝对值小于1的条件被称为CFL判据。从插值的角度分析，只有使用n时刻的值做内插得到n+1时刻的值是稳定的，而外插是不稳定的。稳定时，解的收敛性是有保障的。 2.3 差分格式不同差分格式的稳定性和收敛性是有显著差别的。选择合适的差分格式，不仅能改善计算精度，也会影响到计算效率。因此，在进入数值模式的求解之前，有必要理解差分格式。差分可分为时间差分和空间差分，时间差分是用于近似对t的导数。空间差分则是针对于对空间变量的导数。时间差分可分为显式差分、隐式差分和半隐式差分。 显式差分的格式如下，右端均为n时刻的值。这样做会比较便于计算，因为在计算n+1时刻的值时，n时刻的值是已知的，只需要代入计算即可。从初始条件开始(t=0时刻)，一轮一轮的迭代，即可求得所有值。但收敛条件苛刻，在很多复杂问题上，原始的显式差分不具备计算稳定性，并且会破坏方程的种种守恒性。 隐式差分的格式如下，和上述显式差分格式的区别为右端的时间项是n+1时刻的值。隐式差分在数学上，有很好的收敛性和守恒性。但求解会比显式复杂很多，需要通过矩阵迭代求解。 半隐式差分结合了上述两种格式，可以做到无条件稳定，可以设置时间步长时不需要考虑Courant数，在真实场景中应用广泛。 如果时间采用中央差分，则可以得到蛙跳格式。蛙跳格式在海洋与大气模式中使用较多，等后文结合具体模式再详细讨论。由之前的泰勒展开可知，中央差分是二阶精度，因此该格式的时间和空间是二阶精度。每一次时间迭代时，都会跳过一层，这也是蛙跳格式名字的由来。该方法每计算一层，需要前两层的值已知，是一种显式格式。 针对具体的海洋大气方程进行离散化时，变量比上述例子更多。因此我们还需要考虑网格问题，不同的变量放在不同的位置对计算精度和效率也会有较大的影响。交错网格和非交错网格对计算有着显著的差别。这些会在后文进行一一说明。 版权声明本文创作的初衷是用于帮助数值模式的学习者。欢迎转载，转载请私信并注明作者和出处，请勿用于任何商业用途。 参考书目大气模式、资料同化和可预报性. Kalnay E, 蒲朝霞译. 气象出版社, 2005. Numerical Weather and Climate Prediction. Thomas Tomkins Warner. 2011. Ocean Modelling for Beginners. Jochen Kämpf. Springer Berlin Heidelberg, 2009.
7 ]8 |* I! E( o% E | 
- k# V5 U) X5 ]8 T
/ ^4 J# i5 `9 M" M( h
, n W8 B2 T0 H/ v$ o2 T, }
, F# ~( `5 J- q" E* I* v4 J- Q/ H
8 J4 y& f* V& E0 P
' Y. b9 y" [& V5 z1 s7 m. z7 [+ A
