发展方程计算稳定性的比较分析

对线性与非线性发展方程的几种差分格式,以一维线性和非线性平流方程为例,对线性与非线性发展方程差分格式的计算稳定性进行了比较分析,揭示了差分格式结构和初值形式与非线性计算稳定性的关系。     

强迫耗散非线性发展方程准完全平方守恒格式的构造

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散性发展方程的显式准完全平方守恒差分格式，理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的．值得推广应用。     

非线性Schrdinger方程差分格式的计算稳定性

运用判定非线性发展方程差分格式计算稳定性的Hirt启发性分析方法,对一类非线性Schrdinger方程差分格式的计算稳定性进行分析,得到了保证差分格式计算稳定的必要条件.数值试验结果进一步表明,得到的稳定性判据不仅是保证差分格式计算稳定的必要条件,而且在实际中也是非常有效的.     

发展方程的非线性计算不稳定问题

进一步讨论了有关非线性不稳定的一些问题,其主要内容有： １．考察了有代表性的三类发展方程,指出其对应的差分格式是否出现非线性计算不稳定,与原微分方程解的性质密切相关。 ２．进一步讨论了带周期边条件的守恒型差分格式的非线性计算稳定性问题,总结了克服非线性不稳定的有效措施。 ３．以非线性平流方程为例,着重分析了带非周期边条件的非守恒差分格式的非线性计算稳定性问题,给出了判别其计算稳定性的“综合分析判别法”。     

一类计算稳定性好的显式平流差分格式

通常的显式平流差分格式，如迎风格式，Lax-Wendroff格式等，均是有条件稳定的,其稳定条件与差分网格的时、空步长有关。本文对线性和拟线性平流方程分别构造了一种计算稳定性好的显式差分格式。对前一格式，本文严格证明了它的无条件稳定性及收敛性，并具有一阶精度；对后一格式，由于非线性方程的限制，本文用数值试验研究了它的计算稳定性。     

非线性计算稳定性的比较分析

非线性计算稳定性是数值天气预报中的一个重要问题。本文以一维非线性平流方程为模型,对六种常见的非线性差分格式进行定性和定量的稳定性比较分析,特别强调了非线性计算稳定性与格式结构和初值形式的依赖关系,从另一个角度揭示了非线计算不稳定的特性和机理。     

强迫耗散非线性发展方程显式差分格式的计算稳定性

基于计算准稳定的概念来分析强迫耗散非线性方程显式差分格式的计算稳定性，给出强迫耗散非线性大气方程组显示差分格式计算准稳定的判据，为设计强迫耗散非线性大气方程组计算稳定的显式差分格式提供了新的思路和理论依据。     

非线性发展方程的非守恒格式的计算稳定性问题

针对非线性发展方程的非守恒格式,以一维浅水波方程为例,对非守恒格式的计算稳定性进行了研究分析,探讨了非线性发展方程的非守恒格式与初值的关系.理论分析和数值试验表明,在格式结构已经确定的情况下,非守恒格式的计算稳定性主要由初值的形式所决定.     

耗散与计算稳定性(二)

本文以带耗散项的非线性方程——Burger方程为模型,用截谱形式的差分解代入相应的差分方程,分析了数值天气预报中常用的一类三时间层差分格式的计算稳定性问题,特别讨论了计算初值、时间步长和耗散因子对计算稳定性的影响.利用定性分析和数值试验的方法证实这类三时间层格式有时容易出现非线性计算不稳定,尤其是两步计算初值反号和耗散系数过小的时候是如此.但是耗散系数也不能取得过大,否则也会出现计算不稳定.     

大气海洋环境数值模拟中的若干计算问题

为了适应发展气候数值模拟和环境数值模拟的需要,综合介绍了与此密切有关的若干计算问题。首先,指出这类问题易于出现计算紊乱、非物理解和非线性计算不稳定,同时,简要介绍造成这些问题出现的三种机理：虚假频散、能量关系破坏和能谱的非线性转移。然后,把这类问题归结为一种“发展方程”,给出了与此有关的几个定理,阐述了计算稳定性、能量守恒性和算子非负性之间的密切关系。接着,分别介绍了多种完全平方守恒的差分格式的构造方法,其中包括隐式的和显式的完全平方守恒的差分格式,高时间精度和高阶紧致的完全平方守恒的差分格式,同时也给出若干具体算例说明这些格式的计算效果。最后,对全文作了小结,并指出今后应进一步研究的若干重要问题。     

The relationship between nonconservative schemes and initial values of nonlinear evolution equations

For the nonconservative schemes of the nonlinear evolution equations, taking the one-dimensionalshallow water wave equation as an example, the necessary conditions of computational stability are given.Based on numerical tests, the relationship between the nonlinear computational stability and the construc-tion of difference schemes, as well as the form of initial values, is further discussed. It is proved throughboth theoretical analysis and numerical tests that if the construction of difference schemes is definite, thecomputational stability of nonconservative schemes is decided by the form of initial values.     

A comparative study of conservative and nonconservative schemes

For the conservative and non-conservative schemes of nonlinear evolution equations, by taking the two-dimensional shallow water wave equations as an example, a comparative analysis on computational stability is carried out. The relationship between the nonlinear computational stability, the structure of the difference schemes, and the form of initial values is also discussed.     

The Impacts of Initial Perturbations on the Computational Stability of Nonlinear Evolution Equations

The impacts of initial perturbations on the computational stability of nonlinear evolution equations for non-conservative difference schemes and non-periodic boundary conditions are studied through theoretical analysis and numerical experiments for the case of onedimensional equations.The sensitivity of the difference scheme to initial values is further analyzed.The results show that the computational stability primarily depends on the form of the initial values if the difference scheme and boundary conditions are determined.Thus,the computational stability is sensitive to the initial perturbations.     

再论发展方程差分格式的构造和应用

【摘　要】本文把一大类大气、海洋方程归结为一种发展方程,具体构造了若干定时间步长的显式完全平方守恒差分格式。并证明在一定条件下,这类格式也具有能量守恒、“广义能量”守恒和“平均尺度”守恒的特性,它表明这类格式具有较好的计算稳定性和省时性。文中还探讨了显式平方守恒格式与隐式平方守恒格式之间的密切联系。最后给出了令人满意的用四波的R-H波作数值检验的结果。     

AN IMPROVED SEMI-IMPLICIT TIME DIFFERENCE SCHEME OF SPECTRAL MODEL AND NUMERICAL APPLICATIONS

In fact,the popular semi-implicit time difference scheme of spectral model still includes someimportant linear terms using time explicit difference scheme,and the major terms are directlyrelated to fast internal-and external-gravity waves in the atmospheric forecasting equation.Additionally,due to using time difference on two terms at different time.the popular schemeartificially introduces unbalance between pressure gradient force and Coriolis force terms whilenumerically computing their small difference between large quantities.According to thecomputational stability analysis conducted to the linear term time difference scheme in simpleharmonic motion equation,one improved semi-implicit time difference scheme is also designed inour study.By adopting a kind of revised time-explicit-difference scheme to these linear terms thatstill included in spectral model governing equations,the defect of spectral model which only partlyusing semi-implicit integrating scheme can be overcome effectively.Moreover,besides all spectralcoefficients of prognostic equations,especially of Helmholtz divergence equation,can be workedout without any numerical iteration,the time-step (computation stability) can also be enlarged(enhanced) by properly introducing an adjustable coefficient.     

AN IMPROVED SEMI-IMPLICIT TIME DIFFERENCE SCHEME OF SPECTRAL MODEL AND NUMERICAL APPLICATIONS*

In fact,the popular semi-implicit time difference scheme of spectral model still includes some important linear terms using time explicit difference scheme,and the major terms are directly related to fast internal-and external-gravity waves in the atmospheric forecasting equation.Additionally,due to using time difference on two terms at different time.the popular schemear tificially introduces unbalance between pressure gradient force and Coriolis force terms while numerically computing their small difference between large quantities.According to the computational stability analysis conducted to the linear term time difference scheme in simple harmonic motion equation,one improved semi-implicit time difference scheme is also designed in our study.By adopting a kind of revised time-explicit-difference scheme to these linear terms that still included in spectral model governing equations,the defect of spectral model which only partlyusing semi-implicit integrating scheme can be overcome effectively.Moreover,besides all spectral coefficients of prognostic equations,especially of Helmholtz divergence equation,can be worked out without any numerical iteration,the time-step (computation stability) can also be enlarged (enhanced) by properly introducing an adjustable coefficient.