共查询到20条相似文献，搜索用时 15 毫秒
1.
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. 相似文献
2.
非跳点网格在模式动力—物理过程的耦合方面具有独特的优势,但是由于二阶精度差分方案下非跳点网格频散误差较大而很少被使用于数值天气预报模式。随着近年来数值模式计算精度的不断提高,非跳点网格在频散关系方面的计算误差是否会发生变化还有待研究。本文在高阶精度差分格式下通过浅水波方程对跳点网格和非跳点网格的频散关系进行理论分析和数值试验,主要得到以下结论:（1）在低波数区跳点网格的频散关系基本不随计算精度的提高而变化,但是非跳点网格下的频散关系则随着计算精度的提高而更加接近真实解。在四阶精度下,非跳点网格的频散关系已经非常接近跳点网格。（2）差分精度提高以后,在高波数区非跳点网格仍然存在频率极大值,而且极值中心随着计算精度的提高而逐渐向更高波数区移动。跳点网格在计算精度提高以后高波数区的频率仍然随波数单调增加,且更接近真实解。（3）在高阶精度非跳点网格模拟试验的基础上,结合高阶扩散项对高频短波进行滤除,可以得到与二阶精度跳点网格相接近的模拟结果。总之,在高阶精度有限差分方案下利用非跳点网格构造模式动力框架是一种比较可行的做法。 相似文献
3.
An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy. 相似文献
4.
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. 相似文献
5.
相变修正方案在GRAPES模式标量平流中的应用 总被引：2，自引：1，他引：1
如何更好地模拟水物质的空间分布和小尺度变化,对于数值天气预报效果的改进,特别是对于更好地模拟降水过程,具有重要的意义.计算机的飞速发展使数值模式的分辨率不断提高,云的显式计算成为可能,这样就要求水物质在平流的过程中必须要做到高精度、守恒、保形.水物质场是正定标量的场,具有空间和时间变化幅度大、存在强梯度甚至不连续的特点,水物质场的合理模拟一直是数值预报中的一个难题.GRAPES模式中的标量平流方案采用PRM分段有理函数方法,比较好地解决了该半拉格朗日模式中水物质平流的高精度、守恒、保形问题,但是当有凝结潜热发生时,由于半拉格朗日平流方案求解上游点时的插值,在云边缘区域会造成虚假的云水,进而导致不合理的相变过程.为了解决以上问题,本研究在GRAPES模式中PRM平流方案的案础上,加入了非线性半拉格朗日相变潜热的修正方案,旨在改进GRAPES模式对水物质平流问题的模拟,提高降水的预报效果.该研究通过理想试验,验证了非线性半拉格朗口相变修正方案可以有效地限制云边缘由于半拉格朗日平流方案插值产生的虚假柑变;然后将该方案加入GRAPES模式的PRM水物质平流方案中,通过实际个例模拟验证了加入非线性半拉格朗日方案以后,模式可以更好地模拟水物质的平流过程,且对云中热力场及水物质分布地模拟更加合理,同时预报出的雨带中心区与实况更加符合. 相似文献
6.
GRAPES模式中三维科氏力计算及其效果评估 总被引：1，自引：0，他引：1
作为一种连续可压缩流体,大气具有分层流体特性,其状态变化可由牛顿第二定律、热力学第一定律、连续方程和理想大气状态方程组成的偏微分方程组描述。为了更加精细地描述GRAPES全球模式的动力过程,使模式大气更接近真实大气,在全球非静力模式GRAPES中考虑三维科氏力作用,通过重新构建半隐式半拉格朗日求解大气动力方程组的亥姆霍兹方程系数,在不改变求解方案的前提下实现对GRAPES模式动力过程的更新。然后采用在静力平衡基础上建立的三维大气理想试验对新的动力过程进行数值试验,检验其计算效果和数值稳定性。结果显示,考虑三维科氏力的模式动力框架计算稳定,提高了三维标量和矢量场的计算精度,在水平1°×1°分辨率模式中,平衡流试验第15天计算结果标量场Π'的l1和l2误差分别为0.00023和0.0004,而三维矢量场 V 的l1和l2误差分别为0.002和0.003,均较原模式误差小一个数量级。在罗斯贝-豪威兹波、地形罗斯贝波和斜压波试验中,新框架均表现出很高的计算稳定性和良好的计算效果。 相似文献
7.
伴随技术是四维变分同化（4DVar）系统中计算代价函数梯度的最佳办法,切线性和伴随模式的效果和效率直接影响着4DVar系统的发展。基于GRAPES（Global and Regional Assimilation PrEdiction System）全球切线性和伴随模式1.0版本,利用GRAPES全球模式2.0版本在并行框架和性能等方面的改善,重新优化和设计了GRAPES全球切线性伴随模式2.0版本,提高了GRAPES全球切线性和伴随模式的效果和效率,优化了切线性模式程序结构,使其计算时间最优可控制在非线性模式的1.2倍以内;采用在切线性模式中保存基态的方法,重构了伴随模式的程序结构,使其计算时间最优控制在非线性模式的1.5倍以内;在GRAPES全球切线性物理过程的设计中,将线性物理过程的轨迹基态计算和切线性扰动计算解耦,提高了GRAPES全球切线性和伴随模式的计算效果和效率。 相似文献
8.
基于半拉格朗日 (semi-Lagrangian) 方案的数值天气预报模式, 求解半拉格朗日轨迹上游点变量, 通常采用传统直线逐点拉格朗日多项式插值, 由已知模式格点 (欧拉网格点) 的数值插值获得。对于三维空间上游点的插值, N阶精度需要O(N3) 运算量。N增大, 运算量将大幅增加, 特别耗费计算机机时, 而采用Cascade插值法 (降阶插值法) 则只需要O(N) 运算量。它的显著特点是:用曲线代替直线, 通过一系列中间过渡网格点, 在曲线上用一维拉格朗日插值, 使得相邻拉格朗日格点或中间过渡点的插值不再是孤立的, 而且可以重复使用某些中间结果, 达到减少运算量的目的。将这种方法合理应用于GRAPES模式, 并根据模式的特点, 对Cascade插值过程中独立变量的距离分段计算, 从而有利于实现并行计算。计算结果表明Cascade插值法与传统直线逐点插值法相比, 计算效率平均提高约30%, 同时不降低精度。 相似文献
9.
In this paper,a numerical simulation of the geostrophic adjustment process with C-grid network is illustrated.A difference scheme which has the energy and potential vorticity conserving relation consistent with the differential equations is given,and the effect of some time difference schemes on dispersion of the gravity-inertia wave is discussed.An improved forward-backward time integration scheme is proposed for keeping the computational stability.The effect of various boundary conditions for a finite region model On the gravity-inertia wave is shown by some calculated results. 相似文献
10.
阴阳网格上质量守恒计算性能分析 总被引：3，自引：1，他引：2
质量守恒数值计算是球面准均匀阴阳网格构造全球大气环流模式的重要条件,也是提高阴阳网格应用质量的重要技术手段。本文针对通量形式平流方程,在球面坐标上采用多种理想数值试验对阴阳网格上的三种守恒计算方案和边界插值非守恒计算方案进行了比较检验。发现,质量守恒方案不仅对全球数值积分重要,还影响数值计算精度,满足局地守恒条件的全球强迫守恒方法可以获得较高的精度;网格内质量均匀分布的阴阳网格边界通量一致性守恒强迫计算方案,实现了在不增加计算误差条件下保证局地和全球守恒的目的,且具有很小的计算负担,可以作为阴阳网格上全球质量强迫守恒的有效计算方案;而网格质量的线性分布可以有效提高阴阳网格的数值积分计算精度,但在一定程度上会增加计算负担。 相似文献
11.
The calculation scheme of the smoothed-level and hybrid (SLEVE-hybrid for short) coordinates in numerical forecasting model is not limited to one. It is divided into the semi-analytical scheme and the finite differential scheme in
terms of the various differential methods of the coordinate deformation variables. Comparing the dynamic equation and the long-time batch simulation results of the two schemes, the present study draws the following conclusions. The first-
order finite difference accuracy of the coordinate deformation variables in the finite differential scheme is theoretically lower than that in the semi-analytical scheme. The larger the vertical gradient of the layer thickness is, the larger the
relative errors of the finite differential scheme are. The long-time batch simulation test in the GRAPES model dynamic core demonstrates that the bias of the temperature and the geopotential height in the semi-analytical scheme is smaller
under the default layering, while the simulation difference of the two schemes is greatly reduced when the layering is more uniform. 相似文献
12.
辐射方案对GRAPES全球模式的影响 总被引：10，自引：4，他引：6
辐射过程是中期数值天气预报模式和气候模式中一个重要的物理过程。本文针对GRAPES (Global/Regional Assimilation and PrEdiction System) 全球模式非静力半隐式半拉格朗日的特点, 进行了不同辐射方案的敏感性比较试验, 以期找到一个适合于GRAPES全球模式的辐射方案。用到的方案包括WRF模式中的RRTM长波辐射方案、 NASA的Goddard短波辐射方案以及ECMWF的长短波辐射方案。试验结果表明, ECMWF的辐射方案更适合于GRAPES模式, 30天平均的北半球500 hPa高度场距平相关系数较WRF的辐射方案提高了近2%。另外, 针对实际运行中存在的辐射计算不稳定问题, 对辐射方案进行了气压订正和温度订正, 并对2005年7月影响我国的海棠台风进行了模拟试验, 气压订正后的120小时模式预报在对副热带高压和台风的描述上更接近于分析, 可用预报时效也有了明显提高。 相似文献
13.
本文提出了一种用于双向套网格模式的变格距差分计算方案。该方案在不同格距的网格区采用不同精度的差分格式，它自然地连接粗细网格，避免了一般套网格方案在粗细网格相重合点上进行的重复计算。用解析法和数值试验证明了:它与其它一些变格距差分格式相比，对短波的穿透能力有明显改进，虚假的反射也较小。应用该方案建立了正压原始方程双向套网格模式，并采用空间分解和时间分解计算方法。这不仅使二维问题转化为二个一维问题，而且二维套网格也可简化为一维均匀网格和一维套网格两部分，从而使计算和程序简化。用理想场为初值所做的一系列数值试验表明，该模式中的波可以自由进出粗细网格区，计算稳定。最后，还用该模式做了台风路径预报试验，给出了一些试验结果。 相似文献
14.
为改善球面经纬度网格在高分辨率应用的苛刻限制,提高全球大气动力模式的时间积分效率,选取以阴阳网格为基础构建通量型非静力大气模式动力框架,采用有限体积法和通量型平流显式算法积分方案,保证模式的守恒计算性能。该动力框架在标准三维大气理想试验中进行了中期积分试验,对动力框架的计算效果、性能进行检验。在三维平衡流试验,罗斯贝-豪威茨波试验和山脉罗斯贝波试验中均表现出很好的稳定性和三维计算效果,其中水平2.5°分辨率模式的平衡流垂直速度误差为10-5量级,而经向速度误差在10-2量级,罗斯贝-豪威茨波保持基本波形稳定传播,而地形罗斯贝波试验则给出背风坡激发低槽在发展过程中不断向下游和南半球传播。 相似文献
15.
《热带气象学报(英文版)》2021,(2)
The calculation scheme of the smoothed-level and hybrid(SLEVE-hybrid for short) coordinates in numerical forecasting model is not limited in number. It is divided into the semi-analytical scheme and the finite differential scheme in terms of the various differential methods of the coordinate deformation variables. Having compared the dynamic equation and the long-time batch simulation results of the two schemes, the present study draws the following conclusions. The first-order finite difference accuracy of the coordinate deformation variables in the finite differential scheme is theoretically lower than that in the semi-analytical scheme. The larger the vertical gradient of the layer thickness is, the larger the relative errors of the finite differential scheme are. The long-time batch simulation test in the GRAPES model dynamic core demonstrates that the bias of the temperature and the geopotential height in the semianalytical scheme is smaller under the default layering, while the simulation difference of the two schemes is greatly reduced when the layering is more uniform. 相似文献
16.
针对GRAPES四维变分同化系统的升级, 研究了GRAPES模式垂直扩散方案线性化问题。通过2005年8月7—27日21个个例的批量试验, 发现在GRAPES模式垂直扩散方案源代码的基础上逐句线性化得到的切线性垂直扩散方案即使能通过正确性测试试验, 在少数情况下也会存在很大误差。切线性模式计算的扰动气压场和扰动风场可能出现明显异常, 这种异常与垂直扩散方案中地表动量通量的强非线性有关。如果在切线性垂直扩散方案中忽略地表动量通量扰动, 既可以避免异常的出现, 又不影响其他正常时刻的计算精度。修改后的切线性垂直扩散方案能够在所有变量上一致地提高切线性模式的计算精度。 相似文献
17.
阴阳网格上的质量守恒算法对于阴阳网格在全球模式构建和应用具有重要意义,是模式长期稳定积分和保证计算效果的重要性能指标。本研究在已有的质量均匀分布假定下阴阳网格守恒强迫算法的基础上,构建网格内质量的双线性分布和边界通量线性分布的质量守恒强迫算法,以提高阴阳网格平流计算的精度和模式积分的稳定性。运用CIP-CSLR平流方案对通量形式平流方程数值求解,分别通过"余弦钟"平流试验、正弦波试验和变形流试验对质量双线性分布、边界通量线性分布的新方案与质量和通量均匀分布的原方案进行了对比,标准化误差和标量场分布均表明新方案可有效提高阴阳网格守恒算法的计算效果,且计算负担没有明显增加,具有较好的实用价值。 相似文献
18.
19.
设计了适用于四维变分同化系统的扰动预报模式GRAPES_PF。根据GRAPES的地形追随坐标非静力原始方程组,采用小扰动分离方法推导微分形式的线性扰动预报方程组,并利用与GRAPES非线性模式相似的数值求解方案求解线性扰动微分方程组。在设计扰动预报模式时采用了两个时间层半隐式半拉格朗日方案对动量方程、热力学方程、水汽方程和连续方程进行时间差分,空间差分方案的变量分布水平方向采用Arakawa C跳点网格,垂直方向采用Charney/Phillips跳层。利用代数消元法进一步推导得到只包含未来时刻扰动Exner气压的亥姆霍兹方程,进而通过广义共轭余差法（GCR）求解,在此基础上得到未来时刻扰动量的预报值。基于所开发的扰动模式开展了数值试验。首先在非线性模式中施加一个中尺度初始扰动高压,得到初始扰动在非线性模式中的后续演变,然后将相同的初始扰动作为扰动模式的初值进行时间积分,将扰动模式预报的结果与非线性模式的结果做了对比。结果表明,所开发的扰动模式GRAPES_PF较好地模拟了惯性重力内波的传播过程:初始高压扰动激发了一个迅速向外传播的惯性重力内波,在气压场向风场适应的过程中,水平风场、垂直运动、位温和湿度等变量均出现了扰动增量,与非线性模式得到的结果相当接近。GRAPES_PF作为GRAPES非线性模式的合理线性模式为建立基于线性扰动预报的区域四维变分同化系统奠定了科学基础。 相似文献
20.
Sensitivities of parameterization schemes were conducted based on the Global/Regional Assimilation and Prediction System (GRAPES) model. Surface observations were used to evaluate the simulations and to improve the model’s ability to simulate the extreme precipitation over southern China on 20 July 2016. The results showed that GRAPES captured the large-scale precipitation over southern China but failed to predict the extreme precipitation over Xinyi. The model showed a systematic cold biases by adopting different parameterization schemes. In particular, the ECMWF analyses data showed a strong cold bias over Guangdong province and Guangxi Region. Observational nudging results showed that the surface temperature could largely help to alleviate the cold bias. The alleviation in the warm sector accounted for main improvement by the nudging scheme, and the RMSE was reduced by 1.56 degree from 3.25 degree to 1.69 degree by 1-h simulation and with 1.3 degree alleviation by 2-h simulation. Sensitivities using different parameterizations and the nudging scheme showed that the model’s underestimation of the precipitation was still present despite improvements in the predictions of surface temperature. 相似文献