半隐式半拉格朗日平方守恒计算格式的构造

在显式半拉格朗日完全平方守恒格式基础上, 构造出半隐式半拉格朗日完全平方守恒计算格式, 它继承了显式半拉格朗日完全平方守恒格式的优点, 并突破计算不稳定柯朗条件对时间步长的约束, 使时间步长大为增长.此外, 还给出这种新的计算格式在一维原始方程上的应用.     

一种依据初始场构造的参考大气廓线及对高原地区预报的影响

地形追随坐标系下的动力学方程组通常都要引入参考大气的概念来减少气压梯度力的计算误差.参考大气的一般取法是满足静止和静力平衡关系,实际大气被看成在参考大气上的偏差.实际大气与参考大气越接近,它们之间的偏差就越小,计算就越精确.参考大气的取法一般有等温大气、定常温度廓线、初始资料进行水平平均的垂直温度廓线等.本文在一个非静力的中尺度模式中选取经过水平平均的初始温度场的垂直廓线构造模式的参考大气廓线.但是初始资料在三维空间是离散的,这里根据三次样条函数的分段连续光滑的性质,构造了随起报时间不同而不同的,只随高度变化的解析形式的参考大气廓线.构造的参考大气廓线完全满足静力平衡关系,方程中的静力平衡部分可以得到最大限度地扣除.试验表明,与等温大气相比,这种根据初始状态选取的参考大气和实际资料更接近,在大地形附近气压梯度力的计算误差更小,500 hPa高原附近的24 h气压预报更精确.     

DESIGN AND TEST OF AN IMPROVED SCHEME FOR GLOBAL SPECTRAL MODEL WITH REDUCED TRUNCATION ERROR

An improved scheme of spectral model has been investigated in this paper. Through the introduction of a reference atmosphere, prognostic variables become smoother on the tilted sigma-surface over mountains or a frontal zone, and thus truncation errors are reduced. As a result, the problems about Gibbs waves and negative value of model topography on the sea, which are encountered in the current spectral model in the world, can be solved by the improved scheme. For comparison we have performed the numerical integrations of barotropic and baroclinic Rossby-Haurwitz waves, as well as the prediction of real cases with the improved scheme and the current scheme of spectral model. These results show that the scheme presented in the paper is better than the current scheme of spectral model.     

一个全球大气格点模式和超长波预报问题

本文是作者在“一个修改的σ坐标初始方程”北半球模式的基础上发展全球大气模式的第一个报告.文中给出了这个全球格点模式方程组的基本计算方案;并对一个全球大气形势预报的实例进行了分析,与半球模式预报的结果作了比较,讨论了有关超长波的预报问题.     

数值预报模式动力框架发展的若干问题综述

作者就大气数值模式动力框架发展的几个问题做了回顾和展望.关于模式地形的处理,讨论了"地形追随"(terrain-following)坐标和"台阶地形"(step-mountain)η坐标的优点、问题和对策.关于物理量守恒格式的构造,回顾了从"瞬时"守恒到隐式、显示和半隐式完全(包括时间离散)守恒格式的发展和近况,介绍了加速非线性全隐式问题迭代收敛的途径.还讨论了半拉格朗日守恒格式的构造问题.最后,对目前广泛使用的全球谱模式,讨论了其长处、局限和发展前景,并简单介绍了谱元方法.     

Sensitivity of Medium-Range Weather Forecasts to the Use of Reference Atmosphere

In this paper, the authors develop the earlier work of Chen Jiabin et al. (1986). In order to reduce spectral truncation errors, the reference atmosphere has been introduced in ECMWF model, and the spectrally-represented variables, temperature, geopotential height and orography, are replaced by their deviations from the reference atmosphere. Two modified semi- implicit schemes have been proposed to alleviate the computational instability due to the introduction of reference atmosphere. Concerning the deviation of surface geopotential height from reference atmosphere, an exact computational formulation has been used instead of the approximate one in the earlier work. To re duce aliasing errors in the computations of the deviation of the surface geopotential height, a spectral fit has been used slightly to modify the original Gaussian grid-point values of orography.A series of experiments has been performed in order to assess the impact of the reference atmosphere on ECMWF medium- range forecasts at the r     

一个包括地形和非绝热作用的原始方程数值模式

为了模拟大尺度动力过程并改进数值预报,开展了一个原始方程的半球数值模式的研究。 模式的主要特点是:将热力学方程进行改写,整个差分格式保持总能量守恒,但气压梯度项和静力关系的差分格式可以不受任何限制。模式包括地形、摩擦、水平扩散、降水、蒸发和各种非绝热等物理作用。幅射作用直接由导出的微分表达式计算;地面温度由地面热量平衡方程解出而不需要迭代计算。 最后给出这一模式的一次48小时预报个例的试验结果。     

半拉格朗日平方守恒计算格式的研究

作者研究了一种新的显式平方守恒的半拉格朗日计算格式.该格式是欧拉空间平方守恒格式在半拉格朗日空间的推广和发展,并采用了守恒插值方法.此外,还给出这种格式在一维原始方程上的应用.     

完全非内插半拉格朗日格式在一维 Burgers方程及二维浅水波方程上的应用

在作者过去提出的完全非内插半拉格朗日格式的基础上,针对半拉格朗日格式由于内插带来预报场人为的光滑性问题,进一步发展了这种计算格式,证明了此格式的计算稳定性。为检验这种新的计算格式的性能,在一维和二维问题上进行了应用。在一维问题中采用了一维无粘Burgers方程（方程中有突变点）;二维问题采用了浅水波方程,同时将这些计算结果与Ritchie方案及欧拉方案或一般半拉格朗日内插方案的计算结果进行了比较,发现新格式消除了内插和预报场的人为光滑,并且计算精度有一定程度的提高,这为以后将此格式推广到全球谱模式打下了基础。     

JFNK方法概述及其在大气全隐式非静力模式中的应用方案

首先介绍了近年来新发展的非线性方程全隐式数值求解的JFNK方法,及其在地球流体力学方面应用计算实例.可看到,无论在计算精度还是计算效率方面,全隐式数值求解远远超过常规的半隐式计算格式.其次,还讨论了JFNK方法在气象非静力模式中应用方案,并提出了用静力假定和半隐式差分格式来构造预条件处理器,变三维求解为二维求解,简化了方程组求解难度.该方案不仅可用于差分模式,也为用譜方法求解非静力模式提供可能.