数值海洋与大气模式(一):控制方程和差分离散化

[复制链接]
1.控制方程
海洋与大气的控制方程是由一组偏微分方程表示的。以大气的基本控制方程为例,其中(1)(2)(3)式为运动方程,(4)式为热力学方程,(5)式为水汽平流方程,(6)式为连续方程,可由质量守恒推导而来。具体项的物理含义和数学推导不再赘述。倘若能精确求解如下方程组,即可对天气进行预测或进行数值实验。但该方程组的复杂程度目前难以求得其解析解,因此需要借助数值方法求解下列方程组。

+ H% [( Y! B- O; B# B" ]" v                               
登录/注册后可看大图
上述方程组在进行求解时,考虑到计算开支和应用场景,通常会对模型进行简化。如果在水平尺度远大于垂直尺度时,例如对海岸模拟时,可将控制方程近似得到浅水方程组(Shallow Water Equations)。该方程组的形式相对简单,其解可表示惯性重力波、平流和罗斯贝波等现象,因此该方程组也被称为浅水波方程。因其形式简单,并在河道、河口和近海有较多应用空间,该方程组十分适合用于测试和检验数值方法。在不考虑粘性、风应力等,最简单的浅水方程组的形式如下所示。

1 f9 Q5 b0 r7 d- ^2 w0 Q                               
登录/注册后可看大图
2.离散化与差分离散
介绍完了方程,下面就该到了求解。学过编程语言的同学们肯定知道,计算机只能处理离散过程,我们可以很简单的通过任何一种编程语言去计算
' {- r2 \2 d1 k8 Q1 V
                               
登录/注册后可看大图
,然而我们却无法让计算机直接计算
$ j3 b4 Q! s3 }
                               
登录/注册后可看大图
。因此,数值方法的关键就是如何计算导数。在学习微积分的时候,导数可以这样表示

& w% K$ ?* o$ b+ F- i                               
登录/注册后可看大图
而有限差分法的核心就是使用差分近似微分,当差分的间距取得足够小时,就可以近似计算出导数的值。此时,我们再去计算上式

' i4 \/ L( `* s: c7 ~% v                               
登录/注册后可看大图
时,就可以近似计算

3 `6 H# F" u" T" `3 {% a) l9 y: C                               
登录/注册后可看大图
。这样就完成了一个简单的离散化,并用差分近似微分进行了计算求解。而应用在偏微分方程的求解时,还需要探讨几个问题。
以平流方程为例,该方程形式如下,是一种双曲型偏微分方程。和我们要求解的流体力学方程形式相似。

6 y! X5 r+ c4 K0 Q, A0 u                               
登录/注册后可看大图
不妨取
5 M0 I1 i9 {: _0 f8 n
                               
登录/注册后可看大图
这样的离散点,则差分方程的解
# E2 m* I+ f) K% t: H
                               
登录/注册后可看大图
表示为

) ]/ y: s% Y. G: g( u& O4 N5 P                               
登录/注册后可看大图
,即上标代表时间各点,下标代表空间各点。则平流方程的差分方程可表示成如下形式。
1 A8 u( ~. V; ?( }
                               
登录/注册后可看大图
在对方程求解前,必须探讨两个问题,当时间步长和空间步长趋近于无穷小时,差分方程是否能收敛于偏微分方程?差分方程的解是否能收敛于偏微分方程的解么?前者是在探讨差分方程的一致性,后者是在探讨差分方程的收敛性。
$ i- D( V: l7 y. m
                               
登录/注册后可看大图
不妨将这些项进行泰勒展开,然后代回差分方程,可得下式
  V  I- u1 o% K$ u" V5 F/ ]0 I) i
                               
登录/注册后可看大图
不难发现,差分方程比微分方程多了一些项,而这些项都是随着步长趋于0而趋于0的,因此可以很好的证明一致性。而多出的这些项被称之为截断误差,是差分近似微分时造成的误差。因为我们知道,做泰勒展开时,只有保留所有项才能完全相等,然而展开项是无穷多的。在实际做运算时,针对特定的精度,只能取到相应阶数。而被舍弃掉的高阶小项被称之为截断误差。

# `! x! G  Q/ ~' ^) {  \/ h                               
登录/注册后可看大图

# j( B' v! C% g, V( f/ o; g  N2.2 收敛性
上式差分方程可以化成如下形式,其中
- ?# `6 w8 |. S
                               
登录/注册后可看大图
。不难看出,化成这种形式后,可以发现这种形式的差分,本质上是使用n时刻的值去插值出n+1时刻的值。当

+ x$ s! W7 E) ~9 l$ `1 R# R                               
登录/注册后可看大图
时,就等价于上述我们所说的差分格式。

( K: `- _: u: w4 J. i! A* a                               
登录/注册后可看大图
放缩得之后可推出,只有在
; y) E* a2 w% P- A* n
                               
登录/注册后可看大图
才能保证

% h- Z* J) |* J                               
登录/注册后可看大图
有界。

: H' A4 }- q1 U" Q' j: O                               
登录/注册后可看大图

+ E0 R; p7 W) M                               
登录/注册后可看大图
也被称为Courant数,其绝对值小于1的条件被称为CFL判据。从插值的角度分析,只有使用n时刻的值做内插得到n+1时刻的值是稳定的,而外插是不稳定的。稳定时,解的收敛性是有保障的。
2.3 差分格式
不同差分格式的稳定性和收敛性是有显著差别的。选择合适的差分格式,不仅能改善计算精度,也会影响到计算效率。因此,在进入数值模式的求解之前,有必要理解差分格式。差分可分为时间差分和空间差分,时间差分是用于近似对t的导数。空间差分则是针对于对空间变量的导数。时间差分可分为显式差分、隐式差分和半隐式差分。
显式差分的格式如下,右端均为n时刻的值。这样做会比较便于计算,因为在计算n+1时刻的值时,n时刻的值是已知的,只需要代入计算即可。从初始条件开始(t=0时刻),一轮一轮的迭代,即可求得所有值。但收敛条件苛刻,在很多复杂问题上,原始的显式差分不具备计算稳定性,并且会破坏方程的种种守恒性。
: ]8 B1 n+ o- a1 K$ G& s! v- N$ G
                               
登录/注册后可看大图
隐式差分的格式如下,和上述显式差分格式的区别为右端的时间项是n+1时刻的值。隐式差分在数学上,有很好的收敛性和守恒性。但求解会比显式复杂很多,需要通过矩阵迭代求解。

# @6 E" G: E: I7 C: f- A                               
登录/注册后可看大图
半隐式差分结合了上述两种格式,可以做到无条件稳定,可以设置时间步长时不需要考虑Courant数,在真实场景中应用广泛。
8 ?4 U# Z( T0 \5 J4 s9 n6 x
                               
登录/注册后可看大图
如果时间采用中央差分,则可以得到蛙跳格式。蛙跳格式在海洋与大气模式中使用较多,等后文结合具体模式再详细讨论。由之前的泰勒展开可知,中央差分是二阶精度,因此该格式的时间和空间是二阶精度。每一次时间迭代时,都会跳过一层,这也是蛙跳格式名字的由来。该方法每计算一层,需要前两层的值已知,是一种显式格式。
) M, R" b; }5 U; }; f* Y2 b. Z( x
                               
登录/注册后可看大图
针对具体的海洋大气方程进行离散化时,变量比上述例子更多。因此我们还需要考虑网格问题,不同的变量放在不同的位置对计算精度和效率也会有较大的影响。交错网格和非交错网格对计算有着显著的差别。这些会在后文进行一一说明。
版权声明
本文创作的初衷是用于帮助数值模式的学习者。欢迎转载,转载请私信并注明作者和出处,请勿用于任何商业用途。
参考书目
大气模式、资料同化和可预报性. Kalnay E, 蒲朝霞译. 气象出版社, 2005. Numerical Weather and Climate Prediction. Thomas Tomkins Warner. 2011. Ocean Modelling for Beginners. Jochen Kämpf. Springer Berlin Heidelberg, 2009.

7 @  |! J; Q  k# Y& E: M7 a$ k
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
jnsbmj
活跃在2022-11-6
快速回复 返回顶部 返回列表