共查询到19条相似文献，搜索用时 171 毫秒
1.
A Two-Step Shape-Preserving Advection Scheme 总被引：29，自引：3，他引：29
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. 相似文献
2.
Improvement of the semi-Lagrangian advection scheme in the GRAPES model: Theoretical analysis and idealized tests 总被引：2，自引：0，他引：2
ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided. 相似文献
3.
A moist thermodynamic advection parameter, defined as an absolute value of the dot product of hori- zontal gradients of three-dimensional potential temperature advection and general potential temperature, is introduced to diagnose frontal heavy rainfall events in the north of China. It is shown that the parameter is closely related to observed 6-h accumulative surface rainfall and simulated cloud hydrometeors. Since the parameter is capable of describing the typical vertical structural characteristics of dynamic, thermodynamic and water vapor fields above a strong precipitation region near the front surface, it may serve as a physical tracker to detect precipitable weather systems near to a front. A tendency equation of the parameter was derived in Cartesian coordinates and calculated with the simulation output data of a heavy rainfall event. Results revealed that the advection of the parameter by the three-dimensional velocity vector, the covariance of potential temperature advection by local change of the velocity vector and general potential temperature, and the interaction between potential temperature advection and the source or sink of general potential temperature, accounted for local change in the parameter. This indicated that the parameter was determined by a combination of dynamic processes and cloud microphysical processes. 相似文献
4.
Generalized Square Conservative Multistep Finite Difference Scheme Incorporating Historical Observations 
下载免费PDF全文
下载免费PDF全文 In this paper,a multistep finite difference scheme has been proposed,whose coefficients are determined taking into consideration compatibility and generalized quadratic conservation,as well as incorporating historical observation data.The schemes have three advantages:high-order accuracy in time,generalized square conservation,and smart use of historical observations.Numerical tests based on the one-dimensional linear advection equations suggest that reasonable consideration of accuracy,square conservation,and inclusion of historical observations is critical for good performance of a finite difference scheme. 相似文献
5.
Implementation of the Semi-Lagrangian Advection Scheme on a Quasi-Uniform Overset Grid on a Sphere 总被引：4，自引：0，他引：4
The semi-Lagrangian advection scheme is implemented on a new quasi-uniform overset （Yin-Yang） grid on the sphere. The Yin-Yang grid is a newly developed grid system in spherical geometry with two perpendicularly-oriented latitude-longitude grid components （called Yin and Yang respectively） that overlapp each other, and this effectively avoids the coordinate singularity and the grid convergence near the poles. In this overset grid, the way of transferring data between the Yin and Yang components is the key to maintaining the accuracy and robustness in numerical solutions. A numerical interpolation for boundary data exchange, which maintains the accuracy of the original advection scheme and is computationally efficient, is given in this paper. A standard test of the solid-body advection proposed by Williamson is carried out on the Yin-Yang grid. Numerical results show that the quasi-uniform Yin-Yang grid can get around the problems near the poles, and the numerical accuracy in the original semi-Lagrangian scheme is effectively maintained in the Yin-Yang grid. 相似文献
6.
The breeding method has been widely used to generate ensemble perturbations in ensemble forecasting due to its simple concept and low computational cost. This method produces the fastest growing perturbation modes to catch the growing components in analysis errors. However, the bred vectors(BVs) are evolved on the same dynamical flow, which may increase the dependence of perturbations. In contrast, the nonlinear local Lyapunov vector(NLLV) scheme generates flow-dependent perturbations as in the breeding method, but regularly conducts the Gram–Schmidt reorthonormalization processes on the perturbations. The resulting NLLVs span the fast-growing perturbation subspace efficiently, and thus may grasp more components in analysis errors than the BVs.In this paper, the NLLVs are employed to generate initial ensemble perturbations in a barotropic quasi-geostrophic model.The performances of the ensemble forecasts of the NLLV method are systematically compared to those of the random perturbation(RP) technique, and the BV method, as well as its improved version—the ensemble transform Kalman filter(ETKF)method. The results demonstrate that the RP technique has the worst performance in ensemble forecasts, which indicates the importance of a flow-dependent initialization scheme. The ensemble perturbation subspaces of the NLLV and ETKF methods are preliminarily shown to catch similar components of analysis errors, which exceed that of the BVs. However, the NLLV scheme demonstrates slightly higher ensemble forecast skill than the ETKF scheme. In addition, the NLLV scheme involves a significantly simpler algorithm and less computation time than the ETKF method, and both demonstrate better ensemble forecast skill than the BV scheme. 相似文献
7.
The radon transport test, which is a widely used test case for atmospheric transport models, is carried out to evaluate the tracer advection schemes in the Grid-Point Atmospheric Model of IAP-LASG （GAMIL）. Two of the three available schemes in the model are found to be associated with significant biases in the polar regions and in the upper part of the atmosphere, which implies potentially large errors in the simulation of ozone-like tracers. Theoretical analyses show that inconsistency exists between the advection schemes and the discrete continuity equation in the dynamical core of GAMIL and consequently leads to spurious sources and sinks in the tracer transport equation. The impact of this type of inconsistency is demonstrated by idealized tests and identified as the cause of the aforementioned biases. Other potential effects of this inconsistency are also discussed. Results of this study provide some hints for choosing suitable advection schemes in the GAMIL model. At least for the polax-region-concentrated atmospheric components and the closely correlated chemical species, the Flux-Form Semi-Lagrangian advection scheme produces more reasonable simulations of the large-scale transport processes without significantly increasing the computational expense. 相似文献
8.
DESIGN AND TEST OF AN IMPROVED SCHEME FOR GLOBAL SPECTRAL MODEL WITH REDUCED TRUNCATION ERROR 总被引：3，自引：0，他引：3
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. 相似文献
9.
Tingting MA Guoxiong WU Yimin LIU Zhihong JIANG Jiahui YU 《Journal of Meteorological Research》2019,(3):400-415
The external source/sink of potential vorticity (PV) is the original driving force for the atmospheric circulation. The relationship between surface PV generation and surface PV density forcing is discussed in detail in this paper. Moreover, a case study of the extreme winter freezing rain/snow storm over South China in January 2008 is performed, and the surface PV density forcing over the eastern flank of the Tibetan Plateau (TP) has been found to significantly affect the precipitation over South China in this case. The TP generated PV propagated eastward in the middle troposphere. The associated zonal advection of positive absolute vorticity resulted in the increasing of cyclo-nic relative vorticity in the downstream region of the TP. Ascending air and convergence in the lower troposphere developed, which gave rise to the development of the southerly wind. This favored the increasing of negative meridio-nal absolute vorticity advection in the lower troposphere, which provided a large-scale circulation background conducive to ascending motion such that the absolute vorticity advection increased with height. Consequently, the ascending air further strengthened the southerly wind and the vertical gradient of absolute vorticity advection between the lower and middle troposphere in turn. Under such a situation, the enhanced ascending, together with the moist air transported by the southerly wind, formed the extreme winter precipitation in January 2008 over South China. 相似文献
10.
A QUANTITATIVE STUDY FOR ABNORMAL FREEZING RAINS AND SNOWSTORMS IN SOUTHERN CHINA IN EARLY 2008 总被引：2，自引：1，他引：1
A quantitative diagnosis is carried out for the upward branch of a local meridional circulation over southern China(SC) during the abnormal snowstorms with severe freezing rain from 10 January to 3 February 2008.The diagnostic study shows that the upward branch is mainly associated with the zonal advection of westerly momentum and meridional temperature advection instead of the latent heating(which is commonly the dominant factor in many other storm cases).The corresponding weather analyses indicate that(1) the zonal advection of westerly momentum represents the effect of the upper-level divergence on the anticyclone-shear side in the entrance of a 200 hPa westerly jet with a westward deviation from its climatological location over southwestern Japan;(2) the meridional temperature advection represents the interaction between the mid-lower layer(850 to 400 hPa) warm advection over SC(ahead of temperature and pressure troughs with the latter trough deeper than the former in the Bay of Bengal) and cold advection over north China(steered by an underlying flow at 500 hPa);(3) the relatively weak vapor transport(compared to that of spring,summer and autumn) from the Bay of Bengal and the South China Sea to SC and the existence of a temperature inversion layer in the lower troposphere over SC diminish the effect of latent heating.With the significant increase of vapor transport after 24 January,the role of latent heating is upgraded to become the third positive contributor to the upward branch over SC. 相似文献
11.
Ｅ网格变量分布下差分格式的性质 总被引：8，自引：5，他引：8
本文用线性化的浅水波方程系统地分析、讨论了Ｅ网格变量分布下差分计算的性质；指出了目前有关讨论Ｅ网格差分性质在方法上存在的问题；给出了一种解决Ｅ网格中两个Ｃ网格分离的方法；提出了一种Ｅ网格用于海洋计算时，海洋侧边界条件的处理方法。初步的数值试验证明了这些方法对抑制Ｅ网格中的计算躁声是有效的。 相似文献
12.
通过多种理想试验对正定、保形守恒型有理函数插值半拉格朗日平流方案分别在平面直角坐标以及阴阳网格球面坐标中进行了计算性能分析,并采用多种误差模对守恒型有理函数插值半拉格朗日平流方案的网格收敛性进行评估。结果表明,采用分段有理函数插值的守恒型半拉格朗日平流方案可以有效消除不连续分布处的数值振荡、保证正定性,物理场平滑分布时维持1-2阶收敛速度;而在不连续点或大梯度区域以及应用分维技术的多维算法都会通过有理函数的降阶特性,影响平流计算的收敛阶数,并且,在球面坐标中受球面曲率的影响,守恒型有理函数插值半拉格朗日平流算法的网格收敛速度有所降低。 相似文献
13.
对LASG-REM模式的改进及云雨数值试验 总被引：1，自引：0，他引：1
改进中尺度动力预报模式中较为简单的云雨过程处理方案,用预报方程对大尺度(即网格尺度)的水汽、云水和雨水作显式求解,对次网格尺度的积云对流进行参数化处理,并以1998年7月21日00时(UTC)的欧亚探空资料和地面资料为初始场,对武汉地区的暴雨过程进行数值试验。结果表明,改进后的LASG—REM模式对降水有很好的模拟能力,能够模拟出空中云水、雨水的分布特征。 相似文献
14.
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. 相似文献
15.
两步保形平流方案在高分辨率球面经纬网格下的跳点差分试验 总被引：1，自引：0，他引：1
尝试用一个基于有限差分的两步保形平流方案（TSPAS）替代美国国家大气研究中心第5版公用大气模式（NCAR-CAM5）欧拉-谱动力框架中的半拉格朗日传输（SLT）平流方案,并针对不同分辨率版本,计算水汽等物质的传输过程。通过设计保证守恒性的跳点差分法,解决了基于经纬网格的模式在高分辨率下由于极区纬向格距过小而造成的差分方案积分时间步长过短的问题,使得TSPAS在高分辨率下仍然可以使用较大的时间步长。理想试验表明：（1）跳点TSPAS在守恒性、数值精度（量级上）、保形性等方面均较好地保持了原TSPAS的特点,而积分步长可增大到与CAM5-SLT相同,比原TSPAS提高超过一个数量级;（2）较半拉格朗日传输方案,尽管TSPAS的耗散更大,但是在相同时间步长下更省时,且改进了半拉格朗日传输求解平流方程造成的不守恒问题。CAM5的比较试验表明：（1）采用跳点TSPAS的模式结果与非跳点计算结果相当,在高纬度跳点区域也并未出现模拟异常;（2）在高分辨率下,采用跳点TSPAS方案的CAM5模拟结果与原模式结果相近,并初步显示对东亚青藏高原南侧的虚假降水有所抑制。该工作确保了基于欧拉通量差分型的TSPAS在高分辨率模式下的大时间步长稳定积分,为解决有限差分方法在高分辨率模式发展中的相关计算问题提供了启示,为后续针对东亚地区的模式性能改进奠定了基础。 相似文献
16.
阴阳网格上质量守恒计算性能分析 总被引：3，自引：1，他引：2
质量守恒数值计算是球面准均匀阴阳网格构造全球大气环流模式的重要条件,也是提高阴阳网格应用质量的重要技术手段。本文针对通量形式平流方程,在球面坐标上采用多种理想数值试验对阴阳网格上的三种守恒计算方案和边界插值非守恒计算方案进行了比较检验。发现,质量守恒方案不仅对全球数值积分重要,还影响数值计算精度,满足局地守恒条件的全球强迫守恒方法可以获得较高的精度;网格内质量均匀分布的阴阳网格边界通量一致性守恒强迫计算方案,实现了在不增加计算误差条件下保证局地和全球守恒的目的,且具有很小的计算负担,可以作为阴阳网格上全球质量强迫守恒的有效计算方案;而网格质量的线性分布可以有效提高阴阳网格的数值积分计算精度,但在一定程度上会增加计算负担。 相似文献
17.
AREM及其对2003年汛期降水的模拟 总被引：35，自引：7，他引：35
文中回顾了一个能考虑陡峭复杂地形的有限区域 η坐标数值预报模式 (RegionalEta coordinateModel,简称REM )十几年来的发展和应用过程 ,列举了该模式系统在气象、水文、环境和军事保障等科研和业务单位的主要应用。同时介绍了在此基础上发展升级的暴雨数值预报模式 (AdvancedREM ,简称AREM )。由于 (A)REM抓住了中国区域的地形和天气特点 ,已成为模拟和预报中国天气和环境灾害的重要工具 ,在环境治理、抗洪抢险、防灾减灾以及军事气象保障中发挥了重要作用。通过分析新版本AREM对 2 0 0 3年夏季中国东部区域降水过程的模拟和预报结果 ,反映了AREM对中国东部区域降水的模拟和预报能力 ,再现了 2 0 0 3年夏季中国东部主要降水过程的雨量分布和演变特征。ARM对主要降水区域平均降水量的逐日预报与观测基本一致。它不仅对发生在淮河流域的较大范围的降水预报如此 ,而且对有些受地形强迫影响的局部强降水中心 ,预报结果也与观测基本一致。 相似文献
18.
包含正定水汽高精度算法的MM5模式的暴雨实例试验 总被引：2，自引：1，他引：1
本文首先对非静力中尺度模式ＭＭ５的水汽输送方案作了半理想试验。结果显示,ＭＭ５原有的水汽平 流输送算法在水汽梯度较大的地方产生水汽负值,影响模式的模拟结果和模式的稳定性,对此,本文把一个正定的水汽输送算法Ｐｒａｔｈｅｒ算法引入ＭＭ５,并用实例作了数值试验,检验了模式的模拟能力,结果表明,对ＭＭ５的改进是敏感的有效的,Ｐｒａｔｈｅｒ格式增强了模式暴雨落区,暴雨中心的模拟能力。 相似文献
19.
相变修正方案在GRAPES模式标量平流中的应用 总被引：2，自引：1，他引：1
如何更好地模拟水物质的空间分布和小尺度变化,对于数值天气预报效果的改进,特别是对于更好地模拟降水过程,具有重要的意义.计算机的飞速发展使数值模式的分辨率不断提高,云的显式计算成为可能,这样就要求水物质在平流的过程中必须要做到高精度、守恒、保形.水物质场是正定标量的场,具有空间和时间变化幅度大、存在强梯度甚至不连续的特点,水物质场的合理模拟一直是数值预报中的一个难题.GRAPES模式中的标量平流方案采用PRM分段有理函数方法,比较好地解决了该半拉格朗日模式中水物质平流的高精度、守恒、保形问题,但是当有凝结潜热发生时,由于半拉格朗日平流方案求解上游点时的插值,在云边缘区域会造成虚假的云水,进而导致不合理的相变过程.为了解决以上问题,本研究在GRAPES模式中PRM平流方案的案础上,加入了非线性半拉格朗日相变潜热的修正方案,旨在改进GRAPES模式对水物质平流问题的模拟,提高降水的预报效果.该研究通过理想试验,验证了非线性半拉格朗口相变修正方案可以有效地限制云边缘由于半拉格朗日平流方案插值产生的虚假柑变;然后将该方案加入GRAPES模式的PRM水物质平流方案中,通过实际个例模拟验证了加入非线性半拉格朗日方案以后,模式可以更好地模拟水物质的平流过程,且对云中热力场及水物质分布地模拟更加合理,同时预报出的雨带中心区与实况更加符合. 相似文献