A Two-Step Shape-Preserving Advection Scheme

This paper proposes a new two-step non-oscillatory shape-preserving positive definite finite difference advection transport scheme, which merges the advantages of small dispersion error in the simple first-order upstream scheme and small dissipation error in the simple second-order Lax-Wendroff scheme and is completely different from most of present positive definite advection schemes which are based on revising the upstream scheme results. The proposed scheme is much less time consuming than present shape-preserving or non-oscillatory advection transport schemes and produces results which are comparable to the results obtained from the present more complicated schemes. Elementary tests are also presented to examine the behavior of the scheme.     

Numerical simulation of plume-advection over a hill with coupled models. Part I: Numerical solution of the transport-equation

Summary Four numerical methods suitable for the calculation of transport are compared on a linear two-dimensional advection equation. It is found that a mass-compensation algorithm which removes negative values has also a beneficial effect on the accuracy of schemes which are of higher order but are not positive. It is shown by numerical experiments, that the upstream spline advection scheme combined with a mass compensation is prefereable for some applications rather than other high accurate positive schemes like the flux-corrected transport scheme.With 3 Figures     

水汽输送的算法精度与模式暴雨预报能力的对比分析

利用中尺度模式ＭＭ４，研究了中尺度模式的暴雨模拟对于不同不汽输送方案的敏感性。采用５种平流差分格式来计算水汽的水平输送，用该模式对３次暴雨过程作了数值试验。并把不同方案的降水模拟结果与实况进行了对比。结果表明，中尺度模式的暴雨模拟对水汽输送方案十分敏感，正定高精度的Ｐｒａｔｈｅｒ格式的模拟结果最接近实况。     

高精度迎风偏斜格式的比较与分析

利用一种具有任意阶精度的一般显式有限差分公式构造出高精度迎风偏斜格式，并利用 Fourier分析法评估了这些迎风偏斜格式的耗散误差与频散误差。结果表明，偶阶精度格式的数值相速度快于实际相速度，而奇阶精度格式的数值相速度慢于实际相速度。并且，偶阶精度格式的耗散误差与频散误差低于相邻的奇阶精度格式。为了检验这些格式的计算性能，在一维问题上进行了应用。首先，考虑恒定风场条件下的一维平流试验。主要选择两种不同的初始条件来评价数值格式的精度，这两种试验问题分别是高斯函数、方波函数。试验结果表明，随着数值格式精度的提高，数值格式的误差逐渐减小。而对于高于六阶精度的格式来说，改进的程度并不是很大。其次，应用各阶格式到具有两种不同初始条件的无粘Burgers方程。数值结果表明，随着数值格式阶数的增加，数值结果也得到了明显改进。而对于高于六阶精度的格式来说，进一步的变化并不明显。总之，在兼顾效率与精度条件下六阶迎风偏斜格式是最好的。     

一种新的正定平流方案及其在水汽预报方程中的应用

本文在Smolarkiewicz(1983)方案的基础上,提出了一种新的正定平流方案。它是无条件计算稳定的,形式简单、物理意义清楚且阻尼小,能够很容易地应用到多维模式中。文中通过一简单的平流试验,与其它几种有名的正定平流方案进行了比较;最后将这种方案应用到国家气象中心细网格有限区域业务预报模式水汽预报方程中,避免了原水汽计算方案产生严重负水汽的问题,并对24h降水预报场进行分析。     

Numerical analysis and case experiment for forecasting capability by using high accuracy moisture advectional algorithm

Summary The paper concerns the sensitivity of the regional mesoscale model (MM4) to a moisture transport algorithm with different accuracies, and the effective improvement of the capability to simulate meso-scale synoptical event by using the revised version of the algorithm. The two advective schemes of second-order moment conservational Prather scheme (PRS hereinafter) and the polynomial-fitted fourth-order Bott scheme and the upstream scheme (UPS) are introduced into the MM4 to calculate moisture transport in order to provide other types of the algorithm with different accuracies for choice in addition to the original B grid second-order central scheme. Of the three schemes, the first two ensure the positive definite character with high accuracy (almost no numerical disspation and computational dispersion); and the last one also keeps the positive definite feature, but possesses stronger capability of numerical dissipation. Case studies have been undertaken of several rainstorm and typhoon events by virtue of the different-accuracy moisture transport algorithm, other things being the same. Results reveal great differences between them. An significant improvement is found in determining meso-scale synoptic system structure, rain cluster activity, rainstorm belt and its center, severity and its persistence, and the false alarm of or failure to report subrainstorm center when the PRS scheme is employed. Obviously, these differences and improvement are only due to moisture re-distribution for different schemes of moisture transport and the related variation of latent heat. The paper investigates the influences of different numerical schemes. The high accuracy and positive solution PRS application which produces a significant improvement in the model result, demonstrating the sensitivity of the model to moisture distribution in the model represents an important measure to improve the simulation by means of revising the numerically calculating method.With 11 Figures     

PRM标量平流方案在GRAPES全球预报系统中的应用

如何更好地模拟水物质的分布,对于数值天气预报效果的改进,特别是对于更好地模拟降水过程,具有重要的意义。半拉格朗日模式中的标量平流计算要求做到高精度、守恒、正定和保形,但GRAPES_GFS (Global-Regional Assimilation and PrEdiction System, Global Forecast System) 中采用的QMSL(Quasi-Monotone Semi-Lagrangian)平流方案在水汽的强梯度、不连续区域计算精度较低,且不能做到严格守恒。本研究借鉴计算流体力学领域的研究进展,将一个基于分段有理函数的物质平流方案PRM(Piecewise Rational Method)引入GRAPES_GFS中,按照通量形式求解水汽方程,并对极区进行了混合等技术处理。通过一系列理想试验对两种平流方案进行了对比,证明了PRM方案精度较高,特别是在水汽梯度大的区域优势明显,频散、耗散误差较小,守恒、保形性也要好于QMSL方案。通过对GRAPES_GFS中批量预报试验效果的检验,验证了PRM方案可以有效地改进模式对水物质分布的模拟,提高了降水的预报效果,对模式综合预报性能的提升也有明显作用。     

相变修正方案在GRAPES模式标量平流中的应用

如何更好地模拟水物质的空间分布和小尺度变化,对于数值天气预报效果的改进,特别是对于更好地模拟降水过程,具有重要的意义.计算机的飞速发展使数值模式的分辨率不断提高,云的显式计算成为可能,这样就要求水物质在平流的过程中必须要做到高精度、守恒、保形.水物质场是正定标量的场,具有空间和时间变化幅度大、存在强梯度甚至不连续的特点,水物质场的合理模拟一直是数值预报中的一个难题.GRAPES模式中的标量平流方案采用PRM分段有理函数方法,比较好地解决了该半拉格朗日模式中水物质平流的高精度、守恒、保形问题,但是当有凝结潜热发生时,由于半拉格朗日平流方案求解上游点时的插值,在云边缘区域会造成虚假的云水,进而导致不合理的相变过程.为了解决以上问题,本研究在GRAPES模式中PRM平流方案的案础上,加入了非线性半拉格朗日相变潜热的修正方案,旨在改进GRAPES模式对水物质平流问题的模拟,提高降水的预报效果.该研究通过理想试验,验证了非线性半拉格朗口相变修正方案可以有效地限制云边缘由于半拉格朗日平流方案插值产生的虚假柑变;然后将该方案加入GRAPES模式的PRM水物质平流方案中,通过实际个例模拟验证了加入非线性半拉格朗日方案以后,模式可以更好地模拟水物质的平流过程,且对云中热力场及水物质分布地模拟更加合理,同时预报出的雨带中心区与实况更加符合.     

Application of a Shape-Preserving Advection Scheme to the Moisture Equation in an E-grid Regional Forecast Model

This paper presents a methodology which is very useful to design shape-preserving advection finite difference scheme on general E-grid horizontal arrangement of variables through introducing a two-step shape-preserving positive definite advection scheme in the moisture equation of the LASG-REM (LASG regional E-grid eta-coordinate forecast model). By trial-forecasting six local heavy raincases, the efficiency of the shape-preserving advection scheme in practical application has been examined. The LASG-REM with the shape-preserving advection scheme has a good forecasting ability for local precipitation.     

The influence of numerical advection schemes on the results of ocean general circulation models

The dependence of results from coarse-resolution models of the North Atlantic circulation on the numerical advection algorithm is studied. In particular, the sensitivity of parameters relevant for climate simulations as e.g., meridional transport of mass and heat and main thermocline thickness is investigated. Three algorithms were considered: (a) a central difference scheme with constant values for horizontal and vertical diffusion, (b) an upstream scheme with no explicit diffusion, and (c) a flux-corrected transport (FCT) scheme with constant and strictly isopycnal diffusion. The temporal evolution of the three models on time scales of centuries is markedly different, the upstream scheme resulting in much shorter adjustment time whereas the central difference scheme is slower and controlled by vertical diffusion rather than advection. In the steady state, the main thermocline structure is much less diffusive in the FCT calculation which also has much lower heat transport. Both horizontal circulation and overturning in the meridional-vertical plane are strongest in the upstream-model. The results are discussed in terms of the effective vertical (diapycnal) mixing in the different models. A significant increase in vertical resolution would be required to eliminate the high sensitivity due to the numerical algorithms, and allow physically motivated mixing formulations to become effective.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil