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

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

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

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

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.     

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

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

Three-Step Difference Scheme for Solving Nonlinear Time-Evolution Partial Differential Equations

In this paper, a special three-step difference scheme is applied to the solution of nonlinear time-evolution equations, whose coefficients are determined according to accuracy constraints, necessary conditions of square conservation, and historical observation information under the linear supposition. As in the linear case, the schemes also have obvious superiority in overall performance in the nonlinear case compared with traditional finite difference schemes, e.g., the leapfrog（LF） scheme and the complete square conservation difference（CSCD） scheme that do not use historical observations in determining their coefficients, and the retrospective time integration（RTI） scheme that does not consider compatibility and square conservation. Ideal numerical experiments using the one-dimensional nonlinear advection equation with an exact solution show that this three-step scheme minimizes its root mean square error（RMSE） during the first 2500 integration steps when no shock waves occur in the exact solution, while the RTI scheme outperforms the LF scheme and CSCD scheme only in the first 1000 steps and then becomes the worst in terms of RMSE up to the 2500th step. It is concluded that reasonable consideration of accuracy, square conservation, and historical observations is also critical for good performance of a finite difference scheme for solving nonlinear equations.     

LONG VALID TIME ENERGY PERFECT CONSERVATIVE FIDELITY SPECTRAL SCHEMES OF BAROTROPIC PRIMITIVE EQUATIONS*

In accordance with a new compensation principle of discrete computations,the traditional meteorological global (pseudo-) spectral schemes of barotropic primitive equation (s) are transformed into perfect energy conservative fidelity schemes,thus resolving the problems of both nonlinear computational instability and incomplete energy conservation,and raising the computational efficiency of the traditional schemes.As the numerical tests of the new schemes demonstrate,in solving the problem of energy conservation in operational computations,the new schemes can eliminate the (nonlinear) computational instability and,to some extent even the (nonlinear) computational diverging as found in the traditional schemes,Further contrasts between new and traditional schemes also indicate that,in discrete operational computations,the new scheme in the case of nondivergence is capable of prolonging the valid in-tegral time of the corresponding traditional scheme,and eliminating certain kind of systematical computational "climate drift",meanwhile increasing its computational accuracy and reducing its amount of computation.The working principle of this paper is also applicable to the problem concerning baroclinic primitive equations.     

变步长显式完全平方守恒差分格式

综合隐式完全平方守恒差分格式和显式瞬时平方守恒差分格式的优点，针对一类非线性发展方程构造了一种通过自动调节时间步长来保持平方守恒性的显式差分格式。它基于加小耗散的思想，但又与小耗散法有所不同。本文取的耗散项不是一般的人工耗散，而是取能够弥补由于时间离散所产生的截断误差的所谓（时间）协调耗散。因此，该格式具有较高的时间精度。在数值试验中，该类格式可取得满意的效果。     

AN EXPLICIT COMPLETE SQUARE CONSERVATIVE SCHEME WITH ADJUSTABLE TIME INTERVALS*

In this paper,a kind of explicit difference scheme to solve nonlinear evolution equations,perfectly keeping the square conservation by adjusting the time step interval,is constructed,from the comprehensive maintenance of the advantages of the implicit complete square conservative scheme and the explicit instantaneous square conservative scheme.The new schemes are based on the thought of adding a small dissipation,but it is different from the small dissipation method.The dissipative term used in the new schemes is not a simple artificial dissipative term,but a so-called(time) harmonious dissipative term that can compensate for the truncation errors from the dissociation of the time differential term.Therefore,the new schemes may have a high time precision and may acquire a satisfactory effect in numerical tests.     

大气动力学方程组的定性理论及其应用

基于完整的湿大气动力学方程组,利用无穷维动力系统的新理论和新方法,系统讨论了强迫耗散的非线性大气系统的定性理论及其应用。将完整的强迫耗散非线性湿大气动力学方程组化为Hilbert空间中一个等价的算子方程,研究了算子的性质及其物理意义,在此基础上得到湿大气系统全局吸引子的存在性定理,揭示出系统向外源的非线性适应特征,并把结果推广到有地形动力作用和非定常外源强迫的情形。同时探讨了大气方程组惯性流形的存在,大气多平衡态产生的根源以及强迫、耗散和非线性对系统解的渐近行为的影响。在理论结果的基础上,提出强迫耗散的非线性动力系统中存在三类时间边界层、方程组简化准则、分解算法的算子约束原则以及支撑吸引子基底的少数自由度的构造方法,探讨了理论在非线性发展方程差分格式的设计和计算稳定性分析、多平衡态的数值分析、数值模式延伸预报的改进、短期气候预测以及一类中尺度系统分析与预测中的应用,指出描述长期过程动力学模式的必备条件,给出初值与模式相协调的合理解释。最后,对今后的研究方向作了展望。     

JFNK方法在求解全隐式一维非线性平流方程中的应用

JFNK(Jacobian-free Newton-Krylov)方法是由Newton迭代方法和Krylov子空间迭代方法构成的嵌套迭代方法。作者以全隐式差分的一维非线性平流方程（亦称无粘Burgers方程）探讨采用全隐式格式计算的必要性和JFNK方法的有效性。模拟结果表明, 隐式结果比显式和半隐式结果在稳定度和精度方面较大的优越性, 特别是模拟气流强的系统以及要素空间分布具有较大梯度的系统。     

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

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

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

针对性与非线性发展方程的几种差分格式，以一维线性和非线性平流方程为例，对线性与非线笥发展方程差分格式的计算稳定性进行了比较分析，揭示了差分格式结构和初值形式与计算稳定性的关系。理论分析和数值试验证明，线性与非线性发展方程差分格式的计算稳定性在本质上是完全不同的。     

THE FORMULATION OF FIDELITY SCHEMES OF PHYSICAL CONSERVATION LAWS AND IMPROVEMENTS ON A TRADITIONAL SPECTRAL MODEL OF BAROCLINIC PRIMITIVE EQUATIONS FOR NUMERICAL PREDICTION*

In this paper,two formulation theorems of time-difference fidelity schemes for general quadratic and cubic physical conservation laws are respectively constructed and proved,with earlier major conserving time-discretized schemes given as special cases.These two theorems can provide new mathematical basis for solving basic formulation problems of more types of conservative time-discrete fidelity schemes,and even for formulating conservative temporal-spatial discrete fidelity schemes by combining existing instantly conserving space-discretized schemes.Besides.the two theorems can also solve two large categories of problems about linear and nonlinear computational instability.The traditional global spectral-vertical finite-difference semi-implicit model for baroclinic primitive equations is currently used in many countries in the world for operational weather forecast and numerical simulations of general circulation.The present work,however,based on Theorem 2 formulated in this paper,develops and realizes a high-order total energy conserving semi-implicit time-difference fidelity scheme for global spectral-vertical finite-difference model of baroclinic primitive equations.Prior to this,such a basic formulation problem remains unsolved for long,whether in terms of theory or practice.The total energy conserving semi-implicit scheme formulated here is applicable to real data long-term numerical integration.The experiment of thirteen FGGE data 30-day numerical integration indicates that the new type of total energy conserving semi-implicit fidelity scheme can surely modify the systematic deviation of energy and mass conserving of the traditional scheme.It should be particularly noted that,under the experiment conditions of the present work,the systematic errors induced by the violation of physical laws of conservation in the time-discretized process regarding the traditional scheme designs(called type Z errors for short) can contribute up to one-third of the total systematic root-mean-square(RMS) error at the end of second week of the integration and exceed one half of the total amount four weeks afterwards.In contrast,by realizing a total energy conserving semi-implicit fidelity scheme and thereby eliminating corresponding type Z errors,roughly an average of one-fourth of the RMS errors in the traditional forecast cases can be reduced at the end of second week of the integration,and averagely more than one-third reduced at integral time of four weeks afterwards.In addition,experiment results also reveal that,in a sense,the effects of type Z errors are no less great than that of the real topographic forcing of the model.The prospects of the new type of total energy conserving fidelity schemes are very encouraging.     

THE FORMULATION OF FIDELITY SCHEMES OF PHYSICAL CONSERVATION LAWS AND IMPROVEMENTS ON A TRADITIONAL SPECTRAL MODEL OF BAROCLINIC PRIMITIVE EQUATIONS FOR NUMERICAL PREDICTION

In this paper,two formulation theorems of time-difference fidelity schemes for generalquadratic and cubic physical conservation laws are respectively constructed and proved,with earliermajor conserving time-discretized schemes given as special cases.These two theorems can providenew mathematical basis for solving basic formulation problems of more types of conservative time-discrete fidelity schemes,and even for formulating conservative temporal-spatial discrete fidelityschemes by combining existing instantly conserving space-discretized schemes.Besides.the twotheorems can also solve two large categories of problems about linear and nonlinear computationalinstability.The traditional global spectral-vertical finite-difference semi-implicit model for baroclinicprimitive equations is currently used in many countries in the world for operational weatherforecast and numerical simulations of general circulation.The present work,however,based onTheorem 2 formulated in this paper,develops and realizes a high-order total energy conservingsemi-implicit time-difference fidelity scheme for global spectral-vertical finite-difference model ofbaroclinic primitive equations.Prior to this,such a basic formulation problem remains unsolved forlong,whether in terms of theory or practice.The total energy conserving semi-implicit schemeformulated here is applicable to real data long-term numerical integration.The experiment of thirteen FGGE data 30-day numerical integration indicates that the newtype of total energy conserving semi-implicit fidelity scheme can surely modify the systematicdeviation of energy and mass conserving of the traditional scheme.It should be particularly notedthat,under the experiment conditions of the present work,the systematic errors induced by theviolation of physical laws of conservation in the time-discretized process regarding the traditionalscheme designs(called type Z errors for short)can contribute up to one-third of the totalsystematic root-mean-square(RMS)error at the end of second week of the integration and exceedone half of the total amount four weeks afterwards.In contrast,by realizing a total energyconserving semi-implicit fidelity scheme and thereby eliminating corresponding type Z errors,roughly an average of one-fourth of the RMS errors in the traditional forecast cases can be reducedat the end of second week of the integration,and averagely more than one-third reduced at integraltime of four weeks afterwards.In addition,experiment results also reveal that,in a sense,theeffects of type Z errors are no less great than that of the real topographic forcing of the model.The prospects of the new type of total energy conserving fidelity schemes are very encouraging.     

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

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

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

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

不同云微物理方案对一次飑线过程模拟的影响

本文利用WRF(Weather Research and Forecasting)模式的不同云微物理方案对2009年6月14日发生在华东地区的一次飑线过程进行1 km的高分辨率模拟,探讨不同云微物理方案对飑线模拟的影响。结果表明:双参数方案的模拟效果总体上优于单参数方案,其中WDM6方案模拟效果最佳,能够较好的模拟出强对流回波区、层云区的主要特征。在单参数方案中以WSM6方案最优。Kessler和LIN方案模拟飑线的回波范围偏小,强度偏弱。进一步对比热力和动力场发现,WDM6方案模拟的冷池的面积最大,强度最强,气压最高,飑线前部的入流处风速最大。不同云微物理方案对微物理场的影响较明显,相比单参数方案,双参数方案模拟的水凝物混合比更高,且能够模拟粒子数浓度,更准确地描述了云中的各类粒子特征。     

有效的正压原始方程拟能守恒保真（拟）谱模式

本工作遵循保真计算原理与方法,对正压原始方程气象传统全球（拟）谱模式方案进行改造,构造了正压原始方程拟能完全守恒（拟）谱模式新型保真计算方案,解决了正压原始方程的（非线性）计算稳定性问题和拟能守恒整体性质保持问题,改进了相应正压原始方程气象传统全球（拟）谱模式方案的计算效能。新型保真方案的数值实验表明,计算实践中,新方案在解决拟能守恒问题的同时,可解决（非线性）计算稳定问题,并在一定条件下可解决非线性计算收敛性问题。进一步的比较数值实验还表明,计算实践中,新型保真计算方案在提高相应气象传统方案的计算精度、     

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.     

长时效的正压原始方程能量完全守恒(拟)谱模式

遵循误差反演补偿新计算原理,对正压原始方程传统气象全球拟谱模式方案进行了改造,构造了正压原始方程能量完全守恒全球拟增模式新计算方案,解决了正压原始方程的（非线性）计算稳定性问题和能量守恒整体性质保持问题,改进了相应正压原始方程传统气象全球拟谱模式方案的计算效能。新方案的数值试验表明：在计算实践上,新方案在解决能量守恒问题的同时,可解决（非线性）计算稳定性问题,并在一定条件下可解决非线性计算收敛性问题。进一步的比较数值试验还表明：在计算实践上,新方案具有在提高相应传统气象方案的计算精度,减少其计算量的同时,延长其计算时效,解决其中一类特定“气候漂移”问题方面的效用。本工作原理也适用于斜压原始方程情形。