共查询到20条相似文献，搜索用时 31 毫秒
1.
Qingchang QIN Xueshun SHEN Chungang CHEN Feng XIAO Yongjiu DAI Xingliang LI 《大气科学进展》2019,36(10):1129-1142
A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrain-following grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future. 相似文献
2.
A global transport model is proposed in which a multimoment constrained finite volume (MCV) scheme is applied to a Yin-Yang overset grid. The MCV scheme defines 16 degrees of freedom (DOFs) within each element to build a 2D cubic reconstruction polynomial. The time evolution equations for DOFs are derived from constraint conditions on moments of line-integrated averages (LIA), point values (PV), and values of first-order derivatives (DV). The Yin-Yang grid eliminates polar singularities and results in a quasi-uniform mesh. A limiting projection is designed to remove nonphysical oscillations around discontinuities. Our model was tested against widely used benchmarks; the competitive results reveal that the model is accurate and promising for developing general circulation models. 相似文献
3.
自适应网格在大气海洋问题中的初步应用 总被引：15，自引：4，他引：15
自适应网格法是80年代兴起的通过求解椭圆型方程的边值问题来数值生成网格的一种新方法。它是在任意形状的区域上求偏微分方程的数值解的一种非常有效的工具。该方法抛弃了等距均匀的差分网格，代之以能够自动地适应所研究问题中解的特征的疏密程度不均的曲线网格。如在边界上计算网格与实际边界相重合，在区域内部可任意调节网格点的疏密程度等。本文扼要地介绍了自适应网格的原理及其构造方法。并将其应用于生成南海区域的计算网格以及数值预报台风路径的自适应网格。 相似文献
4.
A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid. Two kinds of moments(i.e., point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell. The PV moment is updated by a conventional semi-Lagrangian method, while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation. Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme, a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment. To achieve exactly positive-definite preserving, two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme. The resulting scheme is inherently conservative, non-negative, and allows a Courant number larger than one.Moreover, the spatial reconstruction can be performed within a single cell, which is very efficient and economical for practical implementation. In addition, a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry, which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry. Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity, and it has the potential to transport the tracers accurately in a real atmospheric model. 相似文献
5.
Forecasts and simulations are varied owing to different allocation of 3-dimensional variables in mesoscale models. No attempts have been made to address the issue of optimizing the simulation with a 3-dimensional variables distribution that should come … 相似文献
6.
A 3D dynamic core of the non-hydrostatic model GRAPES(Global/Regional Assimilation and Prediction System) is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. Three-dimensional Coriolis forcing is introduced to the new core, and full representation of the Coriolis forcing makes it straightforward to share code between the Yin and Yang subdomains. Similar to that in the original GRAPES model, a semi-implicit semi-Lagrangian scheme is adopted for temporal integration and advection with additional arrangement for cross-boundary transport. Under a non-centered second-order temporal and spatial discretization, the dry nonhydrostatic frame is summarized as the solution of an elliptical problem. The resulting Helmholtz equation is solved with the Generalized Conjugate Residual solver in cooperation with the classic Schwarz method. Even though the coefficients of the equation are quite different from those in the original model, the computational procedure of the new core is just the same. The bi-cubic Lagrangian interpolation serves to provide Dirichlet-type boundary conditions with data transfer between the subdomains. The dry core is evaluated with several benchmark test cases, and all the tests display reasonable numerical stability and computing performance. Persistency of the balanced flow and development of both the mountain-induced Rossby wave and Rossby–Haurwitz wave confirms the appropriate installation of the 3D Coriolis terms in the semi-implicit semi-Lagrangian dynamic core on the Yin-Yang grid. 相似文献
7.
阴阳网格上的质量守恒算法对于阴阳网格在全球模式构建和应用具有重要意义,是模式长期稳定积分和保证计算效果的重要性能指标。本研究在已有的质量均匀分布假定下阴阳网格守恒强迫算法的基础上,构建网格内质量的双线性分布和边界通量线性分布的质量守恒强迫算法,以提高阴阳网格平流计算的精度和模式积分的稳定性。运用CIP-CSLR平流方案对通量形式平流方程数值求解,分别通过"余弦钟"平流试验、正弦波试验和变形流试验对质量双线性分布、边界通量线性分布的新方案与质量和通量均匀分布的原方案进行了对比,标准化误差和标量场分布均表明新方案可有效提高阴阳网格守恒算法的计算效果,且计算负担没有明显增加,具有较好的实用价值。 相似文献
8.
Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh. 相似文献
9.
Hann-Ming H. Juang 《Meteorology and Atmospheric Physics》1992,50(1-3):75-88
Summary The set of fully compressible nonhydrostatic equations governing a broad spectrum of atmospheric motion was transformed fromz coordinates to sigma coordinates under a hydrostatic base state. The hydrostatic base state may be either time-independent, such as a hydrostatic balance with-out motion or with motion such as a thermal wind balance, or time-dependent such as might be obtained from the result of integrating a hydrostatic model. The transformed set of equations can be used to predict and study all scales of at mospheric phenomena.The set of perturbation equations was also derived under the same condition. The computational sensitivity in computation of pressure gradient force in sigma coordinate can be improved by computing the pressure gradient interms of perturbations under a certain hydrostatic state.The hydrostatic regional spectral model developed by Juang and Kanamitsu (1991) was modified to be a nonhydrostatic spectral model based on the nonhydrostatic equations in sigma coordinates with time-dependent hydrostatic base states. A semi-implicit time integration scheme was used. Two experiments were performed to test this nonhydrostatic spectral model with acceptable results.With 7 Figures 相似文献
10.
阴阳网格上质量守恒计算性能分析 总被引：3，自引：1，他引：2
质量守恒数值计算是球面准均匀阴阳网格构造全球大气环流模式的重要条件,也是提高阴阳网格应用质量的重要技术手段。本文针对通量形式平流方程,在球面坐标上采用多种理想数值试验对阴阳网格上的三种守恒计算方案和边界插值非守恒计算方案进行了比较检验。发现,质量守恒方案不仅对全球数值积分重要,还影响数值计算精度,满足局地守恒条件的全球强迫守恒方法可以获得较高的精度;网格内质量均匀分布的阴阳网格边界通量一致性守恒强迫计算方案,实现了在不增加计算误差条件下保证局地和全球守恒的目的,且具有很小的计算负担,可以作为阴阳网格上全球质量强迫守恒的有效计算方案;而网格质量的线性分布可以有效提高阴阳网格的数值积分计算精度,但在一定程度上会增加计算负担。 相似文献
11.
本文基于自适应网格的原理，发展了一套网格随时间变的自适应网格模式。该模式保持着 IAP 模式原有的整体积分性质，并具有在任何时刻可随意调整网格点的疏密程度及网格点与边界相重合的特点。可用于台风路径的数值预报等方面。 相似文献
12.
With the development of numerical weather prediction technology,the traditional global hydrostatic models used in many countries of the world for operational weather forecasting and numerical simulations of general circulation have become more and more unfit for high-impact weather prediction.To address this,it is important to invest in the development of global nonhydrostatic models.Few existing nonhydrostatic global models use consistently the grid finite difference scheme for the primitive equations of dynamical cores,which can subsequently degrade the accuracy of the calculations.A new nonhydrostatic global spectral model,which utilizes the Eulerian spectral method,is developed here from NCAR Community Atmosphere Model 3.0(CAM3.0).Using Janjic's hydrostatic/nonhydrostatic method,a global nonhydrostatic spectral method for the primitive equations has been formulated and developed.In order to retain the integrity of the nonhydrostatic equations,the atmospheric curvature correction and eccentricity correction are considered. In this paper,the Held-Suarez idealized test and an idealized baroclinic wave test are first carried out,which shows that the nonhydrostatic global spectral model has similar climate states to the results of many other global models for long-term idealized integration,as well as better simulation ability for short-term idealized integration.Then,a real case experiment is conducted using the new dynamical core with the full physical parameterizations of subgrid-scale physical processes.The 10-day numerical integration indicates a decrease in systematic error and a better simulation of zonal wind,temperature,and 500-hPa height. 相似文献
13.
基于动态自适应网格的开源软件Gerris受到越来越多海洋和水文研究者的关注.概述了Gerris开发背景、研究现状和特点,详细阐述了Gerris数值方案,包括动态自适应网格、动态负载平衡技术原理、广义正交曲线坐标系、内嵌复杂固体边界和地形数据的处理方法,并探讨了Gerris在海洋数值模拟中的初步应用.结果表明,Gerris动态自适应网格在多尺度问题模拟中的优势独特,在海洋数值模拟应用中可通过自适应网格提高地理特征的精度,通过GTS(或KDT)格式的数据来处理地形和网格,达到同时兼顾精确性和易用性的目的,使得Gerris与其他海洋模式进行有机结合成为重要发展方向. 相似文献
14.
Xiaohua Yang 《Boundary-Layer Meteorology》1993,65(4):401-424
This paper describes a nonhydrostatic and incompressible mesoscale model formulation using a terrain-following coordinate system. A tensor transformation procedure is used to derive a diagnostic equation for the nonhydrostatic pressure field. The model features a simplified second-order turbulence closure scheme. The two-dimensional version of the nonhydrostatic model, as well as the corresponding hydrostatic model, are applied to simulate stably stratified airflow over mesoscale bell-shaped mountain ridges. The results show that the nonhydrostatic model is capable of simulating nonhydrostatic dynamics of mesoscale lee wave systems such as the trapped wave phenomenon. 相似文献
15.
Summary The fourth generation of numerical weather prediction (NWP) models is currently under development at the Deutscher Wetterdienst
(DWD) consisting of a global grid point model (GME) and limited-area Lokal-Modell (LM). The nonhydrostatic fully compressible
LM has been designed to meet high-resol ution regional fore-cast requirements at meso-β and meso-γ scales. The initial LM
implementation is based on the NCAR/Penn State MM5 with the addition of a novel generalized terrain-following coordinate and
rotated lat-lon grid. A fully 3D semi-implicit time-stepping scheme has been implemented by retaining the full buoyancy term
instead of the approximate form found in MM5. In contrast with earlier schemes, mass-lumping is not applied to simplify the
elliptic operator on an Arakawa-C/Lorenz grid. The resulting variable-coefficient elliptic problem is solved using a minimal
residual Krylov iterative method with line relaxation preconditione rs. The new semi-implicit scheme is compared with a variant
of the Klemp–Wilhelmson split-explicit scheme (horizontal explicit, vertical implicit) on the basis of computational efficiency
and accuracy at resolutions ranging from 7 km to 400 m. Both idealized 3D mountain wave flows and naturally occuring flows
are analyzed. Below the tropopause, the 3D semi-implicit scheme can be more efficient for low Mach number M ≪ 1 flows when the number of small time steps Δt
s
of the split-explicit approach increases with the sound-speed Courant number.
Revised October 7, 1999 相似文献
16.
17.
非常定自适应网格模式在台风路径数值预报中的应用 总被引：7，自引：0，他引：7
本文将作者在前文中介绍的非常定自适应网格模式用于台风路径的数值预报。由于自适应网格在台风中心附近安排了较密的网格，有效地提高了网格的分辨率，使得台风环流的结构在预报48小时以后仍能保持，其预报的路径与实况相比令人鼓舞。 相似文献
18.
对流层中层中尺度对流涡旋在台风榴莲(2001)生成中的作用——模拟诊断分析 总被引：2，自引：1，他引：1
本文基于PSU/NCAR MM5中尺度模式对台风榴莲 (2001) 生成过程成功的数值模拟, 利用模式输出的较高时空分辨率资料, 对台风榴莲生成过程中对流层中层中尺度对流涡旋 (MCV) 的作用进行了诊断分析。结果表明, 中层MCV在台风榴莲生成中的作用有三个重要方面: 第一, 中尺度组织化作用: 伴随中层MCV的垂直次级环流圈, 使得区域内的积云对流热塔趋向于逐步在中心区域集中, 热塔相互之间容易发生相互作用, 通过合并过程有些热塔得到加强, 而有些趋于消亡。同时, 热塔聚集后的群体效应反馈作用又使得中层MCV加强或维持, 进一步促进热塔的合并以及向轴对称化发展; 第二, 存贮效应: 因为中层MCV的生命史比积云对流热塔长, 能够将消亡对流热塔所携带的热量、 水汽、 涡度加以存贮和保留, 使得中层MCV区域向有利于TC生成的方向发展, 最终成为TC环流的 “胚胎”; 第三, 中层MCV与对流层低层的槽 (涡旋) 以及对流热塔之间通过相互作用, 共同实现中低层系统的垂直耦合。 相似文献
19.
复杂下垫面地域边界层结构的三维细网格数值模拟 总被引：8，自引：0，他引：8
建立了一个可供复杂下垫面地域使用的非静力的三维细网格边界层模式，就复杂下垫面条件下的边界层结构和湍流特征作了以实例为对照的数值模拟试验，模式采用能量闭合方案，舍弃了静力近似。以实测资料为初台输入，同时还做了一些数值试验，分别获得了采用静力与非静力模式和不同闭合方案对ＰＢＬ模型的结果。 相似文献
20.
网格嵌套技术对一次中尺度对流系统降水过程模拟的影响 总被引：6，自引：5，他引：1
利用非静力中尺度模式MM5对2002年7月22日12:00~23日12:00(世界时,下同)长江流域的一次梅雨锋暴雨过程进行数值模拟试验,主要讨论了网格嵌套技术对降水和中尺度对流系统的影响。结果表明:三重嵌套在D1,D2域选用积云参数化方案后,模拟的雨区收缩,虚假降水中心相对减少,降水强度及分布更接近观测值。在模式非线性动力、热力及湿物理过程共同驱动下,通过嵌套网格的双向相互作用,使可分辨云尺度的细网格域D3将其信息通过嵌套边界向次细网格域D2传递,然后再通过D2域边界向粗网格域D1域传递。同样,动力、热力反馈也会反向进行。结果将有助于改进各网格域的预报效果。但对D1网格域系统位置及其发展演变过程的影响相对小些;另外,通过双向多重嵌套,可提高模式预报区域的分辨率,特别是提高模式关键预报区域的分辨率,这也就有可能改进预报的中尺度物理场,使其能够较真实地描写大气实况。 相似文献