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.     

物理守恒律保真格式构造与数值预报斜压原始方程传统谱模式改进研究

文中构造并证明了一般二次和三次物理守恒律时间差分保真格式两个构造定理,以往一些主要时间离散守恒格式构造方案可作为两个定理特例给出。它们不仅可为解决更加广泛类别的时间离散保真格式构造基本问题提供适用数学基础,而且也为结合已有瞬时空间离散守恒格式,解决更加广泛类别的时-空离散意义下保真格式构造基本问题提供适用的数学基础。此外,文中两个定理还可解决两大类问题的线性和非线性计算不稳定性问题。斜压原始方程传统半隐式全球谱-垂直有限差分模式目前是世界上许多国家的业务预报和大气环流模式。本工作利用文中新构定理,构造并且实现了斜压原始方程全球谱-垂直有限差分模式半隐式高阶全能量守恒方案。以往该项基本问题无论在理论还是实践上长期以来一直都未能得到解决。     

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.