共查询到20条相似文献，搜索用时 187 毫秒
1.
大气数值模式离散垂直动量方程时,一般而言气压梯度力和重力不易维持严格的静力平衡关系。为精确平衡数值离散的垂直气压梯度力和重力,基于高精度多矩约束有限体积方法引入完全平衡数值公式,即以满足静力平衡关系的热力学参考态对重力源项进行数值离散构造,发展了适用于非静力大气的完全平衡多矩约束有限体积方法。一维标准数值试验表明,完全平衡多矩约束有限体积方法能在较粗糙的计算网格点上保持静力平衡参考态的数值计算误差在计算机的单精度（10-6）和双精度（10-14）水平,在具有小量级扰动的初始条件下,完全平衡多矩约束有限体积方法能较好地模拟扰动的传播,二维非静力热泡试验进一步验证了完全平衡多矩约束有限体积方法对非静力大气运动的模拟能力。数值试验结果验证了所发展方法的完全平衡属性和适用性,这为非静力大气模式发展提供了良好参考价值。 相似文献
2.
New approach to incorporating the impacts of non-hydrostatic perturbations in atmospheric models 
下载免费PDF全文
下载免费PDF全文 《大气和海洋科学快报》2017,(5)
通过对静力和非静力模式动力框架进行对比分析,研究发现:非静力扰动主要影响大气的垂直运动,对于水平运动的影响相对而言较小,且其影响可通过一个控制其打开和关闭状态的开关性变量来予以表示。基于上述发现,提出了一套新的非静力扰动整合方案。该方案将复杂耗时的偏微分方程数值求解问题转化为基于历史数据的局部函数建模问题,且能够比较方便地应用于一般的静力模式动力框架当中,大大降低了非静力大气环流模式的研发难度和工作量。 相似文献
3.
4.
全球(Z)双三次数值模式的设计与个例模拟(Ⅱ): 个例模拟结果分析 总被引：4，自引：3，他引：1
引入双三次数值模式:通过"双三次曲面拟合—时间步长积分—双三次曲面拟合—……",实现成为一种新算法的数值模式。在球面Z坐标系上,给出全球经纬网格双三次数值模式的一个设计。全球(Z)双三次数值模式采用非静力、全可压、干/湿绝热大气运动动力框架,并采用显式水平方向准拉格朗日/垂直方向欧拉时间积分方案,用NCEP再分析资料做初值场,初步模拟出从极地到热带的基本大气环流。模拟试验表明:因大气运动本质上是非线性的,当限定模式大气作线性"一次运动"时,个例模拟无须任何平滑;但限定模式大气作"三次运动"时,必须有适当平滑。本文认为,未来双三次数值模式可按场变量双三次曲面曲率,以判断采用符合物理诠释的局域(或单点)平滑;又若大气非线性运动不稳定区恰是天气系统发源地,则全球多重/时变套网格双三次数值模式是今后研究的方向。 相似文献
5.
一个描述小尺度深凹地形大气边界层的数值模式及数值模拟试验 总被引：1，自引：0，他引：1
文中针对小尺度深凹地形（深凹露天矿）的特点,建立了一个高分辨率的、三维非静力的大气边界层数值模式,模式中采用能量（E-ε）闭合方案。利用该模式对小尺度深凹地形边界层的三维结构进行了较充分的研究,并且有选择地对影响边界层结构的因子进行了数值试验。 相似文献
6.
引入双三次数值模式:双三次数值模式是通过作三次样条与双三次曲面拟合,实现各个大气要素量场的二阶可导,从而可对各个预报方程作时间积分。双三次数值模式适合采用大气运动原始方程组和采用水平方向准拉格朗日/垂直方向欧拉时间积分方案。且本文的全球(Z)双三次数值模式与个例模拟实际采用Navier-Stokes"浅薄大气"原始方程组,建立球面Z坐标系上的非静力、全可压、干/湿绝热大气运动动力框架。其离散化气压、气温预报方程与个例模拟揭示出大气运动中凝结降水,其天气学原因不仅是湿空气作Z坐标垂直上升运动,而且是湿空气被"减压/减温",即湿空气作P坐标"垂直上升运动",后者可由大气平流运动(如Rossby波)引起。 相似文献
7.
随着计算机和计算理论的发展,数值模式正在向全球化、精细化发展,以适应多尺度、多目的应用的要求,从而模糊了大气环流模式和中尺度数值模式的界限,其主要手段就是改进模式动力框架、离散化手段、计算方法、物理过程的普适性。在实际应用中如何提高模式的多尺度计算性能则是问题关键。该文从模式球面坐标系、网格构造、离散化方法的动力特性、垂直坐标和地形处理以及对物理过程的要求出发,探讨分析多尺度大气数值模式的特点:全球/区域可选、非静力近似、具有良好的频散关系和详细的物理过程,垂直高度坐标和"剪切"地形对多尺度通用模式的改良十分重要。除上述特点外,模式所采用的计算方法也应该最大限度地描述大气动力过程特性,采用高性能计算方案有利于多尺度预报。结合当前多尺度预报的国际研究热点和开发前沿,探讨我国新一代多尺度数值预报系统GRAPES的进一步发展及改进方向。 相似文献
8.
GRAPES模式中三维科氏力计算及其效果评估 总被引：1，自引：0，他引：1
作为一种连续可压缩流体,大气具有分层流体特性,其状态变化可由牛顿第二定律、热力学第一定律、连续方程和理想大气状态方程组成的偏微分方程组描述。为了更加精细地描述GRAPES全球模式的动力过程,使模式大气更接近真实大气,在全球非静力模式GRAPES中考虑三维科氏力作用,通过重新构建半隐式半拉格朗日求解大气动力方程组的亥姆霍兹方程系数,在不改变求解方案的前提下实现对GRAPES模式动力过程的更新。然后采用在静力平衡基础上建立的三维大气理想试验对新的动力过程进行数值试验,检验其计算效果和数值稳定性。结果显示,考虑三维科氏力的模式动力框架计算稳定,提高了三维标量和矢量场的计算精度,在水平1°×1°分辨率模式中,平衡流试验第15天计算结果标量场Π'的l1和l2误差分别为0.00023和0.0004,而三维矢量场 V 的l1和l2误差分别为0.002和0.003,均较原模式误差小一个数量级。在罗斯贝-豪威兹波、地形罗斯贝波和斜压波试验中,新框架均表现出很高的计算稳定性和良好的计算效果。 相似文献
9.
全球（Z）双三次数值模式的设计与个例模拟（Ⅰ）：模式动力框架设计 总被引：3，自引：2，他引：1
引入双三次数值模式：双三次数值模式是通过作三次样条与双三次曲面拟合,实现各个大气要素量场的二阶可导,从而可对各个预报方程作时间积分。双三次数值模式适合采用大气运动原始方程组和采用水平方向准拉格朗日／垂直方向欧拉时间积分方案。且本文的全球（Z）双三次数值模式与个例模拟实际采用Navier—Stokes“浅薄大气”原始方程组,建立球面Z坐标系上的非静力、全可压、干／湿绝热大气运动动力框架。其离散化气压、气温预报方程与个例模拟揭示出大气运动中凝结降水,其天气学原因不仅是湿空气作Z坐标垂直上升运动。而且是湿空气被“减压／减温”,即湿空气作P坐标“垂直上升运动”,后者可由大气平流运动（如Rossby波）引起。 相似文献
10.
11.
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. 相似文献
12.
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. 相似文献
13.
Xiaohua Yang 《Boundary-Layer Meteorology》1991,54(1-2):183-208
Nonhydrostatic effects in two-dimensional mesoscale sea breeze systems are investigated by numerical simulations. It is shown that nonhydrostatic effects are directly contributed by the vertical gradients of the vertical velocity variance as well as by the vertical accelerations. It is also shown that a K-type turbulence closure is not suitable in a nonhydrostatic primitive equation model, and a higher-order closure scheme should therefore be used. Results from hydrostatic and fully-nonhydrostatic models are compared for various surface and atmospheric background conditions, such as scale and strength of surface heating, geostrophic wind, stability, surface roughness contrast, Coriolis effect, etc. It is found that for strongly developed sea breeze cases, vertical gradients of vertical velocity variance contribute most to nonhydrostatic forcing in the lower layers, and that the resultant nonhydrostatic pressure gradient acts against the hydrostatic pressure gradient, so that nonhydrostatic simulations produce weaker systems than hydrostatic ones. For weak sea breeze systems, the difference between the two models tends to be small. 相似文献
14.
15.
Modeling nonhydrostatic atmospheric flow requires the solution of the vertical equation of motion and a prognostic or diagnostic equation for pressure. If the nonhydrostatic components of the flow are relatively small, they can be approximated and incorporated into a purely hydrostatic model, which usually is conceptually simpler and computationally more efficient. A method to do this for a linear model of local thermally-induced circulations is further developed and adapted to a non-linear numerical model of the neutral atmospheric boundary layer. A hydrostatic model and the quasi-nonhydrostatic version were used to simulate neutral flow over simple terrain features. One set of observations taken over a simple change in roughness and another set taken over a change in both roughness and terrain were simulated by both models to assess the capabilities of the quasi-nonhydrostatic technique.It is found that (as expected) the pressure deviation from the hydrostatic state is negligible for the roughness change, but it is an important aspect of neutral flow over terrain. Thus, for flow encountering a simple roughness change, the hydrostatic approximation is good, even for small horizontal scales. However, the quasi-nonhydrostatic model qualitatively produces the features in the observations for flow over a terrain change that the hydrostatic model cannot produce.Journal Paper No. J-12737 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa. Project No. 2779. 相似文献
16.
A parameterization of the nonhydrostatic pressure was modified and adapted to a nonlinear numerical model of the neutral atmospheric boundary layer. A hydrostatic model and the quasinonhydrostatic version were used to simulate neutral flow over a symmetrical hill of uniform roughness. Mean-flow quantities and some turbulence characteristics of the flow from both models are presented. These results were compared with observations, analytic theory, and other numerical models.The quasi-nonhydrostatic method produced qualitative features commonly observed in such flows that the hydrostatic model could not simulate. For instance, the observed velocity reduction at the hill base and the speedup at the summit both were simulated by the quasi-nonhydrostatic model. However, computation of vertical velocities from the incompressible continuity equation is inadequate above regions of recirculation and presents a limitation to the method.Journal Paper No. J-12741 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa. Project No. 2779. 相似文献
17.
深凹地形边界层风场与湍流场结构及扩散规律的数值研究 总被引：1，自引：0，他引：1
建立了一个模拟小尺度复杂地形边界层平均流场和湍流场结构分布特征的三维非静力细网格高阶矩湍流闭合模式,并将此模式与一个拉格朗日粒子随机游动大气扩散数值模式相联结,以边界层气象模式的输出作为扩散模式的输入,构成了一个大气扩和模拟系统,成功地模拟了某地实际深凹露天矿区和理想凹坑地形两种情况下的平均流场、雷诺应力和湍流通量场,分析了其分布特征,并模拟了矿坑内污染物的散布规律 相似文献
18.
J. Steppeler R. Hess U. Schättler L. Bonaventura 《Meteorology and Atmospheric Physics》2003,82(1-4):287-301
Summary ?Currently available computer power allows to run operational numerical weather prediction models at resolutions higher than
10 km. The aim of such high resolution modeling is the prediction of local weather, including orographically induced winds
and local precipitation patterns. In this range the hydrostatic approximation is no longer valid and nonhydrostatic models
have to be used instead. For several decades these models have been developed for research purposes only, but operational
application is now reality. In this paper, the numerical methods used in current nonhydrostatic forecast models will be reviewed
and some promising techniques in this field will be discussed. Special attention is given to aspects such as the choice of
the vertical coordinate, the efficiency of algebraic solvers for semi-implicit time discretizations, and accurate and non-oscillatory
advection schemes.
Received July 6, 2001; revision October 12, 2001 相似文献
19.
大气环境数值模拟研究新进展 总被引：14，自引：1，他引：13
近五年来,中国科学院大气物理研究所(简称大气所)在大气环境数值模拟方面取得了丰硕的成果,通过自主发展和引进,建立了完备的多尺度、多成分的大气环境数值模式,包括全球大气化学输送模式、区域和城市空气质量预报模式。大气所利用这些模式研究各种空间尺度上污染物浓度时空分布以及污染物的输送和演变,研究了多种污染过程的成因和污染变化规律,在污染物输送、低对流层臭氧高污染、区域及城市污染等方面取得了很多成果,并对区域或城市空气质量进行业务化实时预报。大气所还拓展了我国大气环境模拟研究的新领域:大气化学资料同化、污染模式集合预报、污染源反演新方法。初步建立了空气质量模式的资料同化系统(分别基于最优插值技术和集合卡曼滤波技术)和多模式集合预报体系,提高了模式预报水平;在污染源反演新方法方面进行了初步的探索。结合我国目前仍然面临着的大气环境问题,对今后大气环境数值模式的发展方向进行了展望。 相似文献
20.
A high resolution nonhydrostatic tropical atmospheric model is developed by using a ready-made regional atmospheric modeling system. The motivation is to investigate the convective activities associated with the tropical intraseasonal oscillation (ISO) through a cloud resolving calculation. Due to limitations in computing resources, a 2000 kmx2000 km region covering the forefront of an ISO-related westerly is selected as the model domain, in which a cloud-resolving integration with a 5-km horizontal resolution is conducted. The results indicate the importance of stratus-cumulus interactions in the organization of the cloud clusters embedded in the ISO. In addition, comparative integrations with 2-km and 5-km grid sizes are conducted, which suggest no distinctive differences between the two cases although some finer structures of convections are discernible in the 2-km case. The significance of this study resides in supplying a powerful tool for investigating tropical cloud activities without the controversy of cloud parameterizations. The parallel computing method applied in this model allows sufficient usage of computer memory, which is different from the usual method used when parallelizing regional model. Further simulation for the global tropics with a resolution around 5 km is being prepared. 相似文献