首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
阴阳网格上质量守恒计算性能分析   总被引:3,自引:1,他引:2  
李江浩  彭新东 《大气科学》2013,37(4):852-862
质量守恒数值计算是球面准均匀阴阳网格构造全球大气环流模式的重要条件,也是提高阴阳网格应用质量的重要技术手段。本文针对通量形式平流方程,在球面坐标上采用多种理想数值试验对阴阳网格上的三种守恒计算方案和边界插值非守恒计算方案进行了比较检验。发现,质量守恒方案不仅对全球数值积分重要,还影响数值计算精度,满足局地守恒条件的全球强迫守恒方法可以获得较高的精度;网格内质量均匀分布的阴阳网格边界通量一致性守恒强迫计算方案,实现了在不增加计算误差条件下保证局地和全球守恒的目的,且具有很小的计算负担,可以作为阴阳网格上全球质量强迫守恒的有效计算方案;而网格质量的线性分布可以有效提高阴阳网格的数值积分计算精度,但在一定程度上会增加计算负担。  相似文献   

2.
A numerical model is used to study the properties of pollutant dispersion over a large uniformly-sloped surface in the stable atmospheric boundary layer. By simulating the structure of boundary layer flow to improve the advective wind field and vertical eddy diffusivities included in the advection-diffusion equa-tion, this numerical model permits an estimation of the distribution of pollutant concentration for more real-istic atmospheric diffusion conditions.  相似文献   

3.
A model was developed for pollutant dispersion from a point source simulating the Hadera (Israel) power plant stack. The model is based on the NCAR mesoscale meteorological MM4 model that provides the wind fields and coefficients of turbulent diffusion. The model was implemented using an implicit numerical scheme with changing directions. A comparison between the model calculations and an analytical solution for the advection-diffusion equation shows good agreement. Relatively low numerical diffusion of the adopted advection scheme was noted. Results for the hilly region of central Israel are presented for a summer case.  相似文献   

4.
A global baroclinic primitive equation model using the spectral technique has been constructed for short‐ and medium‐range numerical weather prediction. The spectral technique, which is a special case of the Galerkin method, employs spherical harmonic basis functions in the evaluation of all horizontal derivatives. The use of a transform technique allows all the horizontal operations to be performed efficiently and allows physical processes to be evaluated in real space. The model employs a semi‐implicit algorithm for time integration and finite differencing in the vertical. Physical processes include orography, moist convection, large scale precipitation and boundary layer processes.  相似文献   

5.
The influence of atmospheric stability on the behaviour of the third moment of flow velocities observed inside a deciduous forest canopy is examined. Results suggest that buoyancy plays a dominant role in dictating the magnitude of gusts observed inside tall vegetation. Furthermore, an examination of the turbulence recorded throughout leaf fall inside the same forest indicates that larger velocity skewnesses are observed inside a canopy in full leaf than inside a sparse canopy. The behaviour of the measured terms in the non-dimensionalized rate equation of the third moment of canopy flow velocities is also examined. Turbulent diffusion and turbulence gradient interaction terms are largest in stable conditions in the upper canopy layer while these are most important in unstable conditions in the lower canopy layer. In all stability regimes, the turbulent diffusion term is the main source of skewness. The turbulence gradient interaction term, the residual and buoyant production terms all contribute to destroy skewness in stable conditions.  相似文献   

6.
1. IntroductionThe Southern Oscillation and El Ni?o/La Ni?aare anomalous events which occur in the tropical at-mosphere and ocean respectively, and their physicalmechanisms have been successively studied in the re-cent half century. Bjerknes (1969) pointed out thatthe close link between Walk circulation and SST inthe equatorial Pacific is the manifestation of the inter-action of tropical atmospheric and oceanic motions.Therefore, Southern Oscillations and El Ni?o cycleswere described as …  相似文献   

7.
Summary In the present paper, an attempt is made for generalized the atmospheric diffusion operator. This can be accomplished by employing the realizability procedure, to identify a surface operator, that ensures self-adjointness’ of the atmospheric diffusion operator. The dispersion modeling in low wind speeds assumes importance because of the high frequency of occurrence and episodic nature of these poor diffusion conditions. A steady-state mathematical model for hermitized model has been calculated for the dispersion of air pollutants in low winds by taking into account the diffusion in the three coordinate directions and advection along the mean wind. The eddy diffusivities have been parameterized in terms of downwind distance for near source dispersion (Arya, 1995). The constants involved in this parameterization are the squares of intensities of turbulence. An analytical solution for resulting advection-diffusion equation with the physically relevant boundary conditions has been obtained. The solution has been used to simulate the field tracer data collected at IIT Delhi in low wind convective conditions.  相似文献   

8.
张祎  宇如聪  李建  陈昊明 《气象学报》2013,71(6):1089-1102
尝试用一个基于有限差分的两步保形平流方案(TSPAS)替代美国国家大气研究中心第5版公用大气模式(NCAR-CAM5)欧拉-谱动力框架中的半拉格朗日传输(SLT)平流方案,并针对不同分辨率版本,计算水汽等物质的传输过程。通过设计保证守恒性的跳点差分法,解决了基于经纬网格的模式在高分辨率下由于极区纬向格距过小而造成的差分方案积分时间步长过短的问题,使得TSPAS在高分辨率下仍然可以使用较大的时间步长。理想试验表明:(1)跳点TSPAS在守恒性、数值精度(量级上)、保形性等方面均较好地保持了原TSPAS的特点,而积分步长可增大到与CAM5-SLT相同,比原TSPAS提高超过一个数量级;(2)较半拉格朗日传输方案,尽管TSPAS的耗散更大,但是在相同时间步长下更省时,且改进了半拉格朗日传输求解平流方程造成的不守恒问题。CAM5的比较试验表明:(1)采用跳点TSPAS的模式结果与非跳点计算结果相当,在高纬度跳点区域也并未出现模拟异常;(2)在高分辨率下,采用跳点TSPAS方案的CAM5模拟结果与原模式结果相近,并初步显示对东亚青藏高原南侧的虚假降水有所抑制。该工作确保了基于欧拉通量差分型的TSPAS在高分辨率模式下的大时间步长稳定积分,为解决有限差分方法在高分辨率模式发展中的相关计算问题提供了启示,为后续针对东亚地区的模式性能改进奠定了基础。  相似文献   

9.
张耀科 《气象学报》1965,37(1):34-43
本文首先讨论了如何依据空间步长来选取平滑系数。其次考虑在涡度方程中加上粘性项后对微分方程解的影响,讨论如何选取适当的粘性系数。第三讨论涡度方程加粘性项后对差分格式计算稳定性的影响,证明了在常用时间中心差分格式中应把粘性项取在第(n-1)时间层,这样计算才是稳定的。反之,如把粘性项取在第n时间层,则差分格式计算不稳定。并且时间向前差的格式在加上粘性项后也可以稳定,但要求的粘性系数过大。最后讨论了用涡度方程作数值预告时加上平滑过程与加上粘性项之间的关系,从计算上看,二者并不等价。但可以把平滑过程与扩散方程类比,以确定适当的平滑系数。  相似文献   

10.
采用地形坐标系下适合于复杂地形且满足地-气耦合的静力学三维中-β尺度气象预报模式,使用08时的实测流场经变分法初始化后得到的初始流场,模拟研究了上海地区的局地流场时空演变特征。模式较细致地考虑了湍流、长短波辐射、凝结、蒸发、地表热力平衡等因子,在此基础上,使用有限元法数值求解一阶闭合的三维欧拉平流扩散方程,模拟了上海地区的SO2浓度分布。结果表明:上海的局地流场具有明显的日变化特征,受海陆温差的影响明显,日间风速较夜间风速大,风向呈顺时针旋转;在不稳定层结情况下,高架点源排放的高空污染物对地面的污染物浓度有贡献,会造成地面污染物浓度升高。  相似文献   

11.
JFNK方法在求解全隐式一维非线性平流方程中的应用   总被引:1,自引:0,他引:1       下载免费PDF全文
JFNK(Jacobian-free Newton-Krylov)方法是由Newton迭代方法和Krylov子空间迭代方法构成的嵌套迭代方法。作者以全隐式差分的一维非线性平流方程(亦称无粘Burgers方程)探讨采用全隐式格式计算的必要性和JFNK方法的有效性。模拟结果表明, 隐式结果比显式和半隐式结果在稳定度和精度方面较大的优越性, 特别是模拟气流强的系统以及要素空间分布具有较大梯度的系统。  相似文献   

12.
We examine the performance of two steady-state models, a numerical solution of the advection-diffusion equation and the Gaussian plume-model-based AERMOD (the American Meteorological Society/Environmental Protection Agency Regulatory Model), to predict dispersion for surface releases under low wind-speed conditions. A comparison of model estimates with observations from two tracer studies, the Prairie Grass experiment and the Idaho Falls experiment indicates that about 50% of the concentration estimates are within a factor of two of the observations, but the scatter is large: the 95% confidence interval of the ratio of the observed to estimated concentrations is about 4. The model based on the numerical solution of the diffusion equation in combination with the model of Eckman (1994, Atmos Environ 28:265–272) for horizontal spread performs better than AERMOD in explaining the observations. Accounting for meandering of the wind reduces some of the overestimation of concentrations at low wind speeds. The results deteriorate when routine one-level observations are used to construct model inputs. An empirical modification to the similarity estimate of the surface friction velocity reduces the underestimation at low wind speeds.  相似文献   

13.
显式自由表面模型时间分裂格式的一个积分方案   总被引:5,自引:0,他引:5  
讨论了有关显式自由表面模型的时间分裂格式中的分裂误差。在传统的积分方案中,由于采用蛙跳格式,斜压动量方程具有三层的时间离散结构,但是,正压动量方程在相应的时间层上只有二层的离散结构。因此,正压方程与斜压方程的时间离散结构是不相容的。分析结果表明,这是导致分裂误差的原因。本提出一种能保持时间离散结构相容的正压和斜压方程的时间积分方案。该积分方案的分裂误差很小因此特别适合于需要长期积分的气候研究。该方案在北太平洋的一个自由表面模型中得到了应用。  相似文献   

14.
基于2005年盘锦芦苇湿地近地面层湍流通量和微气象梯度的连续观测,研究了芦苇湿地近地面层湍流通量参数化方案。结果表明,盘锦芦苇湿地近地面层经常维持弱稳定和弱不稳定层结。在不稳定层结(-0.4  相似文献   

15.
大气涡旋的螺旋结构   总被引:1,自引:1,他引:1  
利用速度场的分解,说明大气中常见的斑图.速度场可分为变形场和旋转场.若大气动力学方程中只有气压梯度力和科里奥利力平衡,此时速度场只有旋转场,地面天气图上气旋反气旋斑图是闭合的.若再加上摩擦力,则斑图为实际的气旋反气旋螺旋.对三维有水平辐散辐合的涡旋,其斑图为三维螺旋型式,常见是漏斗状.台风小口朝上,龙卷风大口朝上.  相似文献   

16.
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.  相似文献   

17.
Experimental studies have shown that tracers injected over short time intervals can often be measured downwind near the ground for a much longer time than would be expected considering the duration of the release and the wind speed. A two-dimensional advection-diffusion model was used in a numerical study to investigate this effect. The results from the numerical study showed that the rate of decrease in concentrations downwind at a receptor after a source has stopped emitting was more sensitive to magnitude of the wind speed and alongwind mixing coefficient than the rate of increase in concentrations when the source first starts emitting. The lower the wind speed or the greater the mixing, the more slowly the concentration decreased as the region upwind of the receptor continued to contribute tracer until it was depleted by advection and vertical diffusion. These atmospheric processes are able, qualitatively, to account for the slowly decreasing air concentrations.  相似文献   

18.
Eddy-covariance data have been analyzed to investigate the influence of local stability on heat transfer within open canopies. The flux–gradient relationship for heat is derived from the temperature variance equation, and the stability dependence of the flux–gradient relationship is examined and discussed. The results indicate that the strong stability dependence of the nondimensional standard deviation of temperature, and the small contributions of turbulent transport to the temperature variance, lead to a strong stability dependence of the nondimensional temperature gradient within open canopies. Quadrant analysis and hole size analysis were performed for momentum and heat fluxes in the subcanopy, and the results indicate that the contribution of each quadrant to the total flux depends on both the local stability and canopy depth. The intermittency of the turbulent flux does not show a clear dependence on local stability. As the contribution of ejections to the heat flux increases, the vertical flux of the temperature variance changes sign from negative to positive, leading to small temperature variance transport in unstable conditions. Multi-resolution analysis indicates that heat and momentum are transported with different dominant time scales in very unstable conditions, suggesting a different role of local buoyancy in heat and momentum transfer.  相似文献   

19.
A study has been made of the evolution of the planetary boundary layer height (PBLH), the heat flux, and momentum flux using Clarke's Wangara data for a period of two days and two nights, 33/34–34/35. The observed Wangara data are considered as being an output of the first two internal levels of a general circulation model, at 1000 and 2000 m height. A time-dependent equation has been used to forecast, explicitly, the PBLH for both convectively unstable and stable periods. A comparison is made between observed and computed values whenever possible. In the unstable case, the Deardorff model (1974) has been used for the prediction of the unstable PBLH. Part of the stable case study involved a formulation of a time-dependent model for the prediction of the stable PBLH. The solution obtained from the model compared favourably with the results of a model suggested by Khakimov. Contribution number 382.  相似文献   

20.
A study has been made of the evolution of the planetary boundary layer height (PBLH), the heat flux, and momentum flux using Clarke’s Wangara data for a period of two days and two nights, 33/34-34/35. The observed Wangara data are considered as being an output of the first two internal levels of a general circulation model, at 1000 and 2000 m height. A time-dependent equation has been used to forecast, explicitly, the PBLH for both convectively unstable and stable periods. A comparison is made between observed and computed values whenever possible. In the unstable case, the Deardorff model (1974) has been used for the prediction of the unstable PBLH. Part of the stable case study involved a formulation of a time-dependent model for the prediction of the stable PBLH. The solution obtained from the model compared favourably with the results of a model suggested by Khakimov.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号