正压涡度方程中考虑粘性项与平滑过程的一些意见

本文首先讨论了如何依据空间步长来选取平滑系数。其次考虑在涡度方程中加上粘性项后对微分方程解的影响,讨论如何选取适当的粘性系数。第三讨论涡度方程加粘性项后对差分格式计算稳定性的影响,证明了在常用时间中心差分格式中应把粘性项取在第(n-1)时间层,这样计算才是稳定的。反之,如把粘性项取在第n时间层,则差分格式计算不稳定。并且时间向前差的格式在加上粘性项后也可以稳定,但要求的粘性系数过大。最后讨论了用涡度方程作数值预告时加上平滑过程与加上粘性项之间的关系,从计算上看,二者并不等价。但可以把平滑过程与扩散方程类比,以确定适当的平滑系数。     

数值模式中的回溯时间积分格式

在对微分方程初值问题的数值计算中 ,当采用前向或后向差分格式时有两个时间层次 ,而当用中央差分格式时则有 3个时间层次 .由于把预报问题提为初值问题 ,因此只有一个初始场的信息被利用 .本文从扩展信息源的角度出发 ,对微分方程的时间积分提出了一种新格式——回溯时间积分格式 ,并对在大气、海洋、水文等常用的平流方程进行了实例计算 .制约动力系统的微分方程可以写为 : u/ t=F ( u,t) ( 1 )其中 ,u为状态变量 ,F为除局地变化项外的所有项 ,称为源函数 ,设时间间隔为△ t,式 ( 1 )的时间中央差格式为 :u( t △ t) =u( t-△ t) 2△ t…     

两种变时间步长的显式差分算法

针对一类大气、海洋问题，从算子方程出发，在显式完全平方守恒差分格式的基础上提出两种变时间步长的差分算法，使数值计算更加有效。在人工调节时间步长法中，通过人为控制方程中算子能量变化的办法调整时间步长；在自动调节时间步长法中，通过数值试验自动比较出所构造差分格式稳定性的好坏，将它与泰勒展开法及分解算法进行混合应用，验明了自动调节时间步长法具有综合判断差分格式好坏的功能，并得出一些有意义的结论。     

正压原始方程套网格试验

本文用“单向”(one—way)方法进行正压原始方程套网格试验。粗、细网格均采用空间分解准拉格朗日格式和由廖洞贤等提出的与此差分格式一致的内边界条件(即方案二)。试验结果表明,用方案二作出的24小时细网格预报,与按半动量格式和粗网格直接给细网格提供内边界条件(即方案一)所得的结果相比,去除了在流出处出现的不稳定,系统的流入、流出也较好。这两种方案再加上Davies“松弛区域”的方法(分别为方案四,方案三),其结果与上述类似,且方案四比方案二还要好些。     

正压原始方程套网格试验

本文用“单向”(one—way)方法进行正压原始方程套网格试验。粗、细网格均采用空间分解准拉格朗日格式和由廖洞贤等提出的与此差分格式一致的内边界条件(即方案二)。试验结果表明,用方案二作出的24小时细网格预报,与按半动量格式和粗网格直接给细网格提供内边界条件(即方案一)所得的结果相比,去除了在流出处出现的不稳定,系统的流入、流出也较好。这两种方案再加上Davies“松弛区域”的方法(分别为方案四,方案三),其结果与上述类似,且方案四比方案二还要好些。     

耗散与计算稳定性(二)

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

显式自由表面模型时间分裂格式的一个积分方案

讨论了有关显式自由表面模型的时间分裂格式中的分裂误差。在传统的积分方案中，由于采用蛙跳格式，斜压动量方程具有三层的时间离散结构，但是，正压动量方程在相应的时间层上只有二层的离散结构。因此，正压方程与斜压方程的时间离散结构是不相容的。分析结果表明，这是导致分裂误差的原因。本提出一种能保持时间离散结构相容的正压和斜压方程的时间积分方案。该积分方案的分裂误差很小因此特别适合于需要长期积分的气候研究。该方案在北太平洋的一个自由表面模型中得到了应用。     

数值预报讲座 第二讲 北半球五层原始方程模式

1982年2月,北京气象中心开始向全国发布一组新的数值预报图,其中北半球的形势预报采用了修正σ坐标的原始方程五层模式。该模式是对北京气象中心原先使用的三层原始方程模式(简称Λ模式)的一个改进。其特点是在垂直分解上从三层扩充为五层;在水平分解上从540公里改为381公里;在静力关系差分上进一步包括了“时间微分”方案;在时间积分上将原来的欧拉后差     

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

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

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

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

The shallow water equations on the sphere and their Lagrange-Galerkin-solution

Summary ?The shallow water equations are formulated on the sphere in a three-dimensional coordinate system with the aid of tangential velocity components and differential operators. We introduce a modified semi-Lagrangian scheme for the discretization in time. The discretization in space is solved by linear finite elements. The grids we use are regular refinements of a macro triangulation which itself is derived from a highly symmetric polyeder also known as a bucky or soccer ball. The good numerical results show that this combination is a promising approach. The numerical algorithm is stable and its strength is the conservation of mass and energy. Received April 13, 2001; Revised December 18, 2001     

水平气压梯度力差分格式误差的定量分析

通过分析发现地形坐标系中气压梯度力差分格式的计算误差应分为二类，其中，第二类误差在水平坐标面倾斜时出现。这种误差可理解为将水平坐标面上的气压(或位势高度)插回到等高面(或等压面)上的插值误差。 构造气压梯度力差分格式的目的正是为了减小这种第二类误差，而不是其他误差。 误差分析表明静力扣除法和Corby 格式的第二类误差都很小，比一般方法的误差小约一个量级。分析还表明第二类误差产生的主要原因是由于气压随高度非线性变化。依据这一分析思路，本文改进了#As(z)#a坐标系中的静力扣除法。 改进后的方法比改进前的误差小约一个量级。     

Implementation of a Conservative Two-step Shape-Preserving Advection Scheme on a Spherical Icosahedral Hexagonal Geodesic Grid

An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme(TSPAS), was generalized and implemented on a spherical icosahedral hexagonal grid(also referred to as a geodesic grid) to solve the transport equation. The C grid discretization was used for the spatial discretization. To implement TSPAS on an unstructured grid,the original finite-difference scheme was further generalized. The two-step integration utilizes a combination of two separate schemes(a low-order monotone scheme and a high-order scheme that typically cannot ensure monotonicity) to calculate the fluxes at the cell walls(one scheme corresponds to one cell wall). The choice between these two schemes for each edge depends on a pre-updated scalar value using slightly increased fluxes. After the determination of an appropriate scheme, the final integration at a target cell is achieved by summing the fluxes that are computed by the different schemes. The conservative and shape-preserving properties of the generalized scheme are demonstrated. Numerical experiments are conducted at several horizontal resolutions. TSPAS is compared with the Flux Corrected Transport(FCT) approach to demonstrate the differences between the two methods, and several transport tests are performed to examine the accuracy, efficiency and robustness of the two schemes.     

面向四维变分资料同化的NSAS积云深对流参数化方案的简化及线性化研究

GRAPES全球四维变分资料同化系统需要积云深对流参数化方案的线性化与伴随方案,直接采用原始复杂参数化方案进行线性化并不可行,需要发展简化光滑方案来减缓非线性与非连续性特征。GRAPES全球模式采用NSAS积云对流参数化方案,积云深对流对环境的反馈主要通过补偿下沉来实现,研究突出补偿下沉作用,忽略降水蒸发、动量反馈等贡献,形成简化方案。采用输入温、湿度廓线加入不同幅度小扰动方法,评估参数化方案计算的温度、比湿时间倾向对输入扰动的敏感性,检验非线性与非连续特征。提出避免或减缓非连续“开关”的方法,在简化方案的基础上发展了简化光滑方案。简化光滑方案与原始积云深对流方案相比,在对流触发上一致,在对流的位温与比湿倾向、降水的时序模拟等方面相似,而在减缓非线性、避免非连续性方面显著优于原始方案。基于简化光滑方案发展的线性化方案表明,对小于2倍分析增量幅度的扰动,线性化方案可以较好地模拟非线性方案的扰动发展。发展的简化光滑方案具有合理性和实用性。      

数值模式中长波辐射参数化方案的对比试验

应用一个简化的长波辐射参数方案，与郭晓岚－钱永甫模式原有长波辐射参数化方案进行了对比，分析了它们对１００ｈＰａ，３００ｈＰａ，５００ｈＰａ，７００ｈＰａ及地面的高度场、温度场、流场及气压场模拟结果的影响，发现不仅两者结果相当一致，而且与实测结果也符合很好。     

On the use of the Boussinesq equations,the reduced system,and the primitive equations for the computation of geophysical flows

Motivated in part by the mathematical problems associated with the application of open boundary conditions to the hydrostatic primitive equations (PE), Browning et al. (1990, Dyn. Atmos. Oceans, 14: 303–332) proposed the use of the reduced system (RS) of equations to replace PE for oceanographic problems. The RS are essentially the Boussinesq equations (BO) with the non-hydrostatic terms in the vertical momentum equation multiplied by a constant δ² ? 1. This artificially alters the physics (e.g. changing the intemal-inertial wave properties) to facilitate numerical integration, but the changes are assumed to have negligible effects on the dynamics of interest. We assess the accuracy and utility of the RS (following the guidelines for the choice of δ) by comparing numerical finite difference solutions of RS, PE and BO for initial-value problems involving three-dimensional instability of an ocean front and atmospheric frontal development in a two-dimensional Eady wave. Both explicit (BO, PE) and semi-implicit (BOSI, PESI) time-difference schemes are used for the Boussinesq and primitive equations. For RS, the same explicit scheme as for BO is used where δ ? 1 allows larger time steps than with the other explicit models. It is found that relative to BO solutions, the errors for RS are small but increase rapidly and monotonically with increasing δ (over a range consistent with the guidelines) and are greater than the errors for the other models. The use of BOSI allows time steps at least as large as those for RS and results in smaller errors than RS. For these problems, BOSI is the preferable model to replace PE.     

Instability in Lagrangian stochastic trajectory models,and a method for its cure

A number of authors have reported the problem of unrealistic velocities (“rogue trajectories”) when computing the paths of particles in a turbulent flow using modern Lagrangian stochastic (LS) models, and have resorted to ad hoc interventions. We suggest that this problem stems from two causes: (1) unstable modes that are intrinsic to the dynamical system constituted by the generalized Langevin equations, and whose actual triggering (expression) is conditional on the fields of the mean velocity and Reynolds stress tensor and is liable to occur in complex, disturbed flows (which, if computational, will also be imperfect and discontinuous); and, (2) the “stiffness” of the generalized Langevin equations, which implies that the simple stochastic generalization of the Euler scheme usually used to integrate these equations is not sufficient to keep round-off errors under control. These two causes are connected, with the first cause (dynamical instability) exacerbating the second (numerical instability); removing the first cause does not necessarily correct the second, and vice versa. To overcome this problem, we introduce a fractional-step integration scheme that splits the velocity increment into contributions that are linear (U _i ) and nonlinear (U _i U _j ) in the Lagrangian velocity fluctuation vector U, the nonlinear contribution being further split into its diagonal and off-diagonal parts. The linear contribution and the diagonal part of the nonlinear contribution to the solution are computed exactly (analytically) over a finite timestep Δt, allowing any dynamical instabilities in the system to be diagnosed and removed, and circumventing the numerical instability that can potentially result in integrating stiff equations using the commonly applied explicit Euler scheme. We contrast results using this and the primitive Euler integration scheme for computed trajectories in a drastically inhomogeneous urban canopy flow.     

A more accurate scheme for calculating Earth’s skin temperature

The theoretical framework of the vertical discretization of a ground column for calculating Earth’s skin temperature is presented. The suggested discretization is derived from the evenly heat-content discretization with the optimal effective thickness for layer-temperature simulation. For the same level number, the suggested discretization is more accurate in skin temperature as well as surface ground heat flux simulations than those used in some state-of-the-art models. A proposed scheme (“op(3,2,0)”) can reduce the normalized root–mean–square error (or RMSE/STD ratio) of the calculated surface ground heat flux of a cropland site significantly to 2% (or 0.9 W m^?2), from 11% (or 5 W m^?2) by a 5-layer scheme used in ECMWF, from 19% (or 8 W m^?2) by a 5-layer scheme used in ECHAM, and from 74% (or 32 W m^?2) by a single-layer scheme used in the UCLA GCM. Better accuracy can be achieved by including more layers to the vertical discretization. Similar improvements are expected for other locations with different land types since the numerical error is inherited into the models for all the land types. The proposed scheme can be easily implemented into state-of-the-art climate models for the temperature simulation of snow, ice and soil.     

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.