《数值天气预报》中非线性计算不稳定的举例.docVIP

《数值天气预报》中非线性计算不稳定的举例 　　摘要：数值天气预报是大气科学的一个重要分支，是一门实用性很强的应用基础学科。数值天气预报本质是用数值方法求解非线性的大气运动方程组。但当采用格点差分来表示微分方程中的非线性项时，易产生非线性计算不稳定现象。本文以一个简单的一维非线性平流方程的数值求解过程为例，给出隐式格式差分方程的基本计算方法并重新演示了非线性计算不稳定现象。由于举例更为简单，可加深学生对显式格式和隐式格式优缺点的理解。 　　关键词：大气科学；非线性计算不稳定；隐式格式 　　中图分类号：G642.0 文献标志码：A 文章编号：1674-9324（2013）42-0160-02 　　一、引言 　　根据物理定律，比如牛顿第二定律、质量守恒定律、能量守恒定律、气体试验定律等，可以得到支配大气运动的基本方程组。但由于大气运动基本方程组是一组高度非线性的偏微分方程组，很难求得其解析解，人们可通过数值的方法求其近似解（数值解），这就是数值天气预报。数值天气预报是一门实用性很强的应用基础学科[1]。通过差分方法求解大气运动基本方程组时，人们发现数值积分过程中会产生计算不稳定问题，这就需要采用恰当的差分格式或积分格式。在《数值天气预报》课程关于时间积分格式的讲解中[2]，重点介绍了“显示”格式，即差分方程的右端项全部为当前（或和过去）时刻的变量值，通过积数值分可求出方程左端的未来时刻变量值。但对于“隐式”及“半隐式”格式只是简单提及其概念和特点，比如隐式格式为用未来时刻变量值求出未来时刻变量值，它具有计算稳定、但计算复杂的特点。 　　在《数值天气预报》课程中，需讲解大量的公式推导和讲解，若不能配合简单而又形象的举例和图形，学生尤其是本科生作为授课对象，将很难理解和接受本课程中的相关内容，讲课的效果也将大打折扣。在非线性不稳定计算的举例中，教材中[2]虽给出了采用不同的初值和不同的差分方案对计算稳定性的影响，但并没有清楚地列出其求解过程，因此学生很难了解隐式格式差分的具体求解过程，对教材中列出的“显示”和“隐式”格式的各自优缺点更是难以理解。因此，需要对两种格式的计算过程进行相应的讲解，尤其是隐式格式。此外，教材中[2]对同一个微分方程构造的两个不同的差分方程中，除隐式格式和显式格式的差异外，还存在着对■采用了不同的差分格式，即显式格式采用中央差格式，隐式格式采用前差格式。本课程[2]已清楚的讲解到中央差格式虽具有较高的计算精度，但在时间差分计算时存在计算解的问题，若初值取得不当，则计算解会有较大的振幅。因此，从逻辑上讲，教材中给出的不同的差分方案的影响，实际上不仅仅来源于显式格式和隐式格式的差异，还来源于对时间微分采用不同差分格式的差异。这又加大了学生对显式格式和隐式格式特点的理解难度。 　　针对上述问题，本文将以简单的一维非线性平流方程为例，给出隐式格式差分方程的具体求解过程，重新探讨非线性计算不稳定现象，目的是使学生更好地了解显式格式和隐式格式差分方程的求解过程，深刻理解两种格式各自的优缺点。 　　二、非线性计算不稳定的计算实例 　　以大气科学中极具代表性的一维平流方程为例：■+u■=0，0≤x1。它可以写成 　　■+■（■）=0，0≤x1 （1） 　　或■+■（u■+■），0≤x1 （2） 　　以上两式与教材[2]基本一致，不同的是这里的x取值范围并不到1。在大气科学中，方程或模式的计算可在全球或某一纬圈上进行。在该情况下，没有纬向侧边界条件。对于上式而言，可认为u在x=1的取值等于u在x=0的取值，也即循环边界条件。在构造上述微分方程相应的差分方程过程中，对（1）式和（2）式分别采用显式格式和隐式格式： 　　uin+1=uin-■[（ui+1n+uin）2-（uin+ui-+1n）2] （3） 　　uin+1=uin-■[（■i+1+■i+■i-1）（■i+1-■i-1）] （4） 　　其中上标n为第n步，下标i为第i个格点，■i=（uin+1+uin）/2。可见，与教材中不同的是，（1）式和（2）式中■均取了前差格式，这样可避免由于三个时间层计算而出现的计算解问题，有利于问题的讨论更加集中。 　　同样给定两种不同的初值，两者仅相差一个常数： 　　ui0=sin2πiΔx （5） 　　ui0=1.5+sin2πiΔx （6） 　　计算中，Δx取=1/3，Δt=0.004，则|u■|≤umax=|u■|=2.5×0.004×3=0.0031，根据线性稳定性判据，其满足线性稳定性条件。将两种差分方案与两种初值两两组合，可得四种计算结果。在线性稳定性条件情况下，若结果仍出现不稳定现象，则一定是非线性计算不稳定。为检查计算稳定性，计算每一步中所有格点的总动能■∑ui2，若该值在某一有限区域内变