首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
包含正定水汽高精度算法的MM5模式的暴雨实例试验   总被引:2,自引:1,他引:1  
叶志宁  葛孝贞 《气象科学》2000,20(2):136-142
本文首先对非静力中尺度模式MM5的水汽输送方案作了半理想试验。结果显示,MM5原有的水汽平 流输送算法在水汽梯度较大的地方产生水汽负值,影响模式的模拟结果和模式的稳定性,对此,本文把一个正定的水汽输送算法Prather算法引入MM5,并用实例作了数值试验,检验了模式的模拟能力,结果表明,对MM5的改进是敏感的有效的,Prather格式增强了模式暴雨落区,暴雨中心的模拟能力。  相似文献   

2.
采用隐式立方样条计算平流过程的数值模式及理想实验   总被引:1,自引:0,他引:1  
肖锋  程麟生 《大气科学》1992,16(5):538-547
本文介绍了一种用隐式立方样条求解平流方程的数值方法,并从理论上对其无条件稳定性进行了证明,在此基础上建立了一个在地形坐标系下的两维原始方程模式,模式在行星边界层参数化中引入了湍流动能方程,在模式顶部引入了吸收层.数值实验表明:模式有较好的计算稳定性,对较高的模式水平分辨率和复杂地形均有较强的适应能力;对复杂地形和下垫面非均匀热源条件下中尺度系统的模拟能获得合理的结果,并具有较高的精确度.  相似文献   

3.
4.
Previous studies concerning the interaction of dual vortices have been made generally in the determin-istic framework. In this paper, by using an advection equation model, eight numerical experiments whose integration times are 30 h are performed in order to analyze the interaction of dual vortices and the vortex self-organization in a coexisting system of deterministic and stochastic components. The stochastic compo-nents are introduced into the model by the way that the Iwayama scheme is used to produce the randomly distributed small-scale vortices which are then added into the initial field. The different intensity of the small-scale vortices is described by parameter K being 0.0, 0.4, 0.6, 0.8, and 1.0, respectively. When there is no small-scale vortex (K=0.0), two initially separated meso-beta vortices rotate counterclockwise mutu-ally, and their quasi-final flow pattern is still two separated vortices; after initially incorporating small-scale vortices (K=0.8, 1.0), the two separated meso-beta vortices of initially same intensity gradually evolve into a major and a secondary vortex in time integration. The major vortex pulls the secondary one, which gradually evolves into the spiral band of the major vortex. The quasi-final flow pattern is a self-organized vortex with typhoon-like circulation, and the relative vorticity at its center increases with increasing in K value, suggesting that small-scale vortices feed the self-organized vortex with vorticity. This may be a pos-sible mechanism responsible for changes in the strength of the self-organized vortex. Results also show that the quasi-final pattern not only relates with the initial intensity of the small-scale vortices, but also with their initial distribution. In addition, three experiments are also performed in the case of various boundary conditions. Firstly, the periodic condition is used on the E-W boundary, but the fixed condition on the S-N boundary; secondly, the fixed condition is set on all the boundaries; and thirdly, the periodic condition is chosen on all the boundaries. Their quasi-final flow patterns in the three experiments are the same with each other, exhibiting a larger scale typhoon-like circulation. Based on these results mentioned above, authors think that the transition of vortex self-organization study from the deterministic system to the coexisting system of deterministic and stochastic components is worth exploring.  相似文献   

5.
罗哲贤  李春虎 《气象学报》2007,65(6):856-863
以往双涡相互作用的动力学一般都在决定性的框架内研究。文中用一个平流方程模式,实施积分时间为30 h的8组试验,分析决定性和随机性共存系统中双涡相互作用和涡旋自组织的问题。随机性通过以下方式引入模式:先用Iwayama方案生成随机分布的小尺度涡,再将这些小尺度涡加入初始场。试验中,初始随机分布小尺度涡的强度参数K分别取0.0、0.4、0.6、0.8和1.0。结果表明,没有小尺度涡的条件下(K=0.0),初始分离的两个β中尺度涡逆时针互旋,其准终态流型是两个分离的涡;引进小尺度涡后,K取0.8、1.0时,初始分离强度相同的两个β中尺度涡,逐渐形成主次之分。主涡将次涡拉伸成为螺旋带,其准终态流型是一个自组织起来的类似于台风环流的涡旋。准终态涡中心的相对涡度值随K值的加大而加大。结果还表明,准终态流型不仅与初始小尺度涡的强度参数有关,而且与初始小尺度涡的分布有关。此外,在相同初始场的情况下,还实施了3类不同边畀条件的试验:第1类,在东西边界取周期条件,在南北边界取固定条件;第2类,在所有边界均取固定条件;第3类,在所有边界均取周期条件。这3类试验的准终态流型相同,都显示出一个类似于台风涡旋的环流。根据这些结果可以初步认为,涡旋自组织的研究从决定性动力学向随机动力学的过渡是值得探索的。  相似文献   

6.
We revisit the analytical model for atmospheric boundary-layer flow over a hill covered with a canopy of Finnigan and Belcher (Q J R Meteorol Soc 130:1–29, 2004). Remaining within the overall scope of that analysis we extend in two ways. First we include the impacts of the advection terms within the upper canopy in a simple, but approximate, manner. Second we establish a modification for the associated pressure perturbation. Both extensions allow us to extend the parameter range wherein the analytical framework can be expected to reasonably hold. The within-canopy advection terms act to provide a downstream shift, and decreased magnitude, to the flow perturbations within the canopy as compared to the predictions from the original analysis. Through continuity, similar, but smaller, impacts are seen above the canopy. Together these act to reduce the differences in the streamwise positions of the topographic speed-up seen above and within the canopy. The modified pressure perturbation also acts to decrease the magnitude of the flow perturbations. The predicted topographic influence on the flow is reduced from that given in the original analysis but, importantly, the positions where the topographic influences most strongly affect the flow, and by extension the scalar concentration fields, are also changed. Predictions from the revised analysis are shown to be in good agreement with wind-tunnel data for flow over an isolated narrow ridge covered by a partially dense canopy.  相似文献   

7.
Numerical simulations of scalar transport in neutral flow over forested ridges are performed using both a 1.5-order mixing-length closure scheme and a large-eddy simulation. Such scalar transport (particularly of CO2) has been a significant motivation for dynamical studies of forest canopy–atmosphere interactions. Results from the 1.5-order mixing-length simulations show that hills for which there is significant mean flow into and out of the canopy are more efficient at transporting scalars from the canopy to the boundary layer above. For the case with a source in the canopy this leads to lower mean concentrations of tracer within the canopy, although they can be very large horizontal variations over the hill. These variations are closed linked to flow separation and recirculation in the canopy and can lead to maximum concentrations near the separation point that exceed those over flat ground. Simple scaling arguments building on the analytical model of Finnigan and Belcher (Q J Roy Meteorol Soc 130:1–29, 2004) successfully predict the variations in scalar concentration near the canopy top over a range of hills. Interestingly this analysis suggests that variations in the components of the turbulent transport term, rather than advection, give rise to the leading order variations in scalar concentration. The scaling arguments provide a quantitative measure of the role of advection, and suggest that for smaller/steeper hills and deeper/sparser canopies advection will be more important. This agrees well with results from the numerical simulations. A large-eddy simulation is used to support the results from the mixing-length closure model and to allow more detailed investigation of the turbulent transport of scalars within and above the canopy. Scalar concentration profiles are very similar in both models, despite the fact that there are significant differences in the turbulent transport, highlighted by the strong variations in the turbulent Schmidt number both in the vertical and across the hill in the large-eddy simulation that are not represented in the mixing-length model.  相似文献   

8.
A non-hydrostatic numerical model, the Regional Atmospheric Modeling System (RAMS), has been used to investigate the development of katabatic jumps in Coats Land, Antarctica. In the control run with a 5 m s-1downslope directed initial wind, a katabatic jump develops near the foot of the idealized slope. The jump is manifested as a rapid deceleration of the downslope flow and a change from supercritical to subcritical flow, in a hydraulic sense, i.e., the Froude number (Fr) of the flow changes from Fr > 1 to Fr> 1. Results from sensitivity experiments show that an increase in the upstream flow rate strengthens the jump, while an increase in the downstream inversion-layer depth results in a retreat of the jump. Hydraulic theory and Bernoulli's theorem have been used to explain the surface pressure change across the jump. It is found that hydraulic theory always underestimates the surface pressure change, while Bernoulli's theorem provides a satisfactory estimation. An analysis of the downs balance for the katabatic jump indicates that the important forces are those related to the pressure gradient, advection and, to a lesser extent, the turbulent momentum divergence. The development of katabatic jumps can be divided into two phases. In phase I, the t gradient force is nearly balanced by advection, while in phase II, the pressure gradient force is counterbalanced by turbulent momentum divergence. The upslope pressure gradient force associated with a pool of cold air over the ice shelf facilitates the formation of the katabatic jump.  相似文献   

9.
复杂山区上空垂直速度场和热力对流活动的理想数值模拟   总被引:2,自引:1,他引:2  
利用英国气象局高分辨率的边界层数值模式BLASIUS,针对中国西北一个复杂山区进行了一系列的理想数值模拟,分析了在不同天气条件下山区上空的垂直速度场分布和对流特征以及地形对热力对流活动的影响,同时讨论了与地形有关的对流触发机制。模式结果表明,复杂山区的垂直运动在稳定层结和风速较大的情况下较易预测,而在中性层结下,山区上空的垂直运动分布随机性强。在Froude数小于0.5的条件下,气流往往被山峰阻塞而在迎风坡造成地形强迫和辐合性抬升,从而易在迎风坡触发深对流活动;在背风坡则由于迎风坡的绕流重新辐合而造成垂直运动。绕流的辐合是触发深对流活动的另一重要因子。在大风或Froude数较大的条件下,地形重力波容易在山地下游被激发。地形重力波与对流活动的相互影响在模式中清楚可见。在适当的条件下,重力波除了可以与对流活动相耦合从而使气团上升到更高的高度外,重力波的走向很可能会影响到深对流系统的传播路径。研究还发现稳定度对相邻两条对流线之间的距离长短也有影响。稳定度较小时,相邻两条对流线之间的平均宽度趋向变大而单个对流线的强度也相应变大。定量化的结论和理论升华值得进一步的数值模拟研究。  相似文献   

10.
Summary The development of a cold front influenced by orography and large scale forcing is examined with a two-dimensional meso-scale model. The model is based on the primitive equations and uses the hydrostatic and anelastic approximations. Gradients of the basic flow and temperature field in the third dimension are taken into account during the simulations. Low diffusive numerical schemes and radiation boundary conditions reduce the numerical errors to an acceptable minimum for a two day simulation and avoid reflections at the upper and lateral boundaries. Frontogenetical forcing is included in the simulations by specifying either a vertically sheared or horizontally convergent basic zonal flow field. Model runs with an idealized cold front were carried out over flat terrain and in the presence of a bell shaped mountain ridge.The simulations show a weakening of the cold front on the windward side of the mountain ridge and a strong reintensification on the leeward side relative to the control runs without topography. Analysis of frontogenesis terms demonstrates the importance of convergence in the ageostrophic circulation and of along-front temperature advection for the development of the cold front. The strong intensification of the cold front on the leeward side of the mountain ridge can only partly be explained by superposition with the mountain induced wave. It is mainly caused by ageostrophic deformation forcing in the strong downward flow of this wave.The results also show that the cold front passage over the mountain ridge is not a continuous process. The formation of a new frontal structure on the leeward side of the mountain ridge, well separated from the primary one, is observed while the initial cold front still exists in the upslope region. Generally nonlinear interactions between the mountain wave and the cold front are the important mechanisms to explain these phenomena.With 18 Figures  相似文献   

11.
Previous studies concerning the interaction of dual vortices have been made generally in the deterministic framework. In this paper, by using an advection equation model, eight numerical experiments whose integration times are 30 h are performed in order to analyze the interaction of dual vortices and the vortex self-organization in a coexisting system of deterministic and stochastic components. The stochastic components are introduced into the model by the way that the Iwayama scheme is used to produce the randomly distributed small-scale vortices which are then added into the initial field. The different intensity of the small-scale vortices is described by parameter K being 0.0, 0.4, 0.6, 0.8, and 1.0, respectively. When there is no small-scale vortex (K=0.0), two initially separated meso-beta vortices rotate counterclockwise mutually, and their quasi-final flow pattern is still two separated vortices; after initially incorporating small-scale vortices (K=0.8, 1.0), the two separated meso-beta vortices of initially same intensity gradually evolve into a major and a secondary vortex in time integration. The major vortex pulls the secondary one, which gradually evolves into the spiral band of the major vortex. The quasi-final flow pattern is a self-organized vortex with typhoon-like circulation, and the relative vorticity at its center increases with increasing in K value, suggesting that small-scale vortices feed the self-organized vortex with vorticity. This may be a possible mechanism responsible for changes in the strength of the self-organized vortex. Results also show that the quasi-final pattern not only relates with the initial intensity of the small-scale vortices, but also with their initial distribution. In addition, three experiments are also performed in the case of various boundary conditions. Firstly, the periodic condition is used on the E-W boundary, but the fixed condition on the S-N boundary; secondly, the fixed condition is set on all the boundaries; and thirdly, the periodic conditio  相似文献   

12.
Interdecadal variability in the North Pacific region is investigated in a 500-y control integration of the Hamburg ECHAM+LSG coupled ocean-atmosphere general circulation model. The spectrum is predominantly red, but a significant peak with a period of about 18 y is detected in the spectrum of sea surface temperature (SST). This peak is shown to be associated with an irregular oscillation that involves both the model ocean and atmosphere. The SST, sea-level pressure, and geopotential height at 500 hPa all undergo a primarily standing oscillation with an extensive monopole structure centered near the date line. The surface anticyclone is situated to the northeast of the warm SST anomaly, and there is a small westward tilt with height; temporal changes are approximately in phase. The anomalous surface heat flux accompanying the warm phase of SST is primarily out of the ocean, but is compensated by anomalous warm advection by surface currents, allowing the SST anomaly to persist. Oceanic thermocline anomalies propagate northward in the western Pacific, and lag the atmosphere indicating a disequilibrium with the atmosphere; sub-surface thermal advection appears to play an important role. A comparison is made between the model's 18-y oscillation and oscillatory components identified in an analysis of the GISST observational SST dataset, which have periods of approximately 6 and roughly 30 y.  相似文献   

13.
Summary Melting snow is often patchy, and advection of heat from warm patches of snow-free ground can accelerate the melting. For large-scale atmospheric modelling applications, snowmelt rates and fluxes of heat and moisture over surfaces with snow cover distributed on scales too small to be resolved by the model grid have to be parametrized. In this paper, a boundary-layer model is used to model advection over surfaces with short vegetation and varying fractions of snow cover. Boundary-layer model results are used to assess the performance of a tile model, which calculates separate fluxes for snow-covered and snow-free fractions of the surface given area-average temperatures, humidities and windspeeds at a reference height in the atmosphere. It is found that the tile model can give good results, but its performance depends strongly on the choice of reference height. Received August 21, 1997 Revised February 23, 1998  相似文献   

14.
T. Zhou  Z. Li 《Climate Dynamics》2002,19(2):167-180
The East Asia summer monsoon (EASM) is simulated with a variable resolution global atmospheric general circulation model (GCM) developed at the Laboratoire de Météorologie Dynamique, France. The version used has a local zoom centered on China. This study validates the model's capability in reproducing the fundamental features of the EASM. The monsoon behaviors over East Asia revealed by the ECMWF reanalysis data are also addressed systematically, providing as observational evidence. The mean state of the EASM is generally portrayed well in the model, including the large-scale monsoon airflows, the monsoonal meridional circulation, the cross-equatorial low-level jets, the monsoon trough in the South China Sea, the surface cold high in Australia, and the upper-level northeasterly return flow. While the performance of simulating large-scale monsoonal climate is encouraging, the model's main deficiency lies in the rainfall. The marked rainbelt observed along the Yangtze River Valley is missed in the simulation. This is due to the weakly reproduced monsoonal components in essence and is directly related to the weak western Pacific subtropical high, which leads to a fragile subtropical southwest monsoon on its western flank and results in a weaker convergence of the southwest monsoon flow with the midlatitude westerlies. The excessively westward extension of the high, together with the distorted Indian low, also makes the contribution of the tropical southwest monsoon to the moisture convergence over the Yangtze River Valley too weak in the model. The insufficient plateau heating and the resulting weak land-sea thermal contrast are responsible for the weakly reproduced monsoon. It is the deficiency of the model in handling the low-level cloud cover over the plateau rather than the horizontal resolution and the associated depiction of plateau topography that results in the insufficient plateau heating. Comparison with the simulation employing regular coarser mesh model reveals that the local zoom technique improves, in a general manner, the EASM simulation.  相似文献   

15.
北半球阻塞高压的维持Ⅰ:准地转和Ertel位涡分析   总被引:2,自引:3,他引:2  
刘辉  吴国雄  曾庆存 《气象学报》1995,53(2):177-185
从北半球不同地区的4个阻塞高压个例,研究了阻塞高压维持机制及其地域性差异。对阻塞高压区域300hPa时间平均准地转让涡低值区的维持机制的分析表明,该位涡低值区的维持机制存在显著的地域性差异。其中,在2个太平洋阻高个例中,时间平均流的位涡平流使低位涡向阻高西北部输送,从而有利于位涡低值区经向度的维持;时变扰动的位涡输送则有抵消平均流作用的倾向,从而对位涡低值区的维持起耗散作用。在大西洋和东亚阻高个例中,平均流的位涡平流使低位涡区向东移动,而时变扰动的位涡输送则有抵消平均流位涡平流的作用,从而有利于位涡低值区在原地稳定维持。以上主要是相对涡度输送的差异所致。等熵面Ertel位涡分析表明,阻高区域330K时间平均等熵位涡低值区的维持机制与300hPa时间平均准地转位涡低值区的维持机制十分相似,从而表明以上等压面准地转位涡分析可以近似用来代表等熵Ertel位涡分析。  相似文献   

16.
The viscous semigeostrophic solutions obtained for the baroclinic Eady wave fronts are analyzed for the generation of the cross-frontal temperature gradient in the boundary layer. In the case of free-slip boundaries, the cross-frontal gradient is maximally generated at the surface by meridional temperature advection. In the case of no-slip boundaries, surface friction reduces the meridional temperature advection in the boundary layer: The maximum generation occurs above the surface layer and the temperature gradient at the surface is maintained by vertical diffusion. The no-slip solution is compared with the Ekman-layer model solution. Errors are quantified for the use of the Ekman-layer model in the mature state of frontogenesis.The surface frontogenesis is found to be affected by diffusivity both directly and indirectly. The direct effect of diffusivity is represented explicitly by the diffusion term in the potential temperature equation. The indirect effect of diffusivity is related implicitly to the temperature advection caused by the viscous part of the ageostrophic motion whose horizontal velocity component is defined by the frictional wind deflection (away from the geostrophy). The direct effect of diffusivity is frontolytical, whilst theindirect effect of diffusivity is frontogenetic in the mesoscale vicinity of the front. The indirect effect of diffusivity contributes dominantly to the mesoscale surface frontogenesis for the free-slip case, but it is offset by the divergence of the dynamic part of the ageostrophic motion at the surface level for the non-slip case.  相似文献   

17.
We study turbulent flow over two-dimensional hills. The Reynolds stresses are represented by a second-order closure model, where advection, diffusion, production and dissipation processes are all accounted for. We solve a full set of primitive non-hydrostatic dynamic equations for mean flow quantities using a finite-difference numerical method. The model predictions for the mean velocity and Reynolds stresses are compared with the measured data from a wind-tunnel experiment that simulates the atmospheric boundary layer. The agreement is good. The performance of the second-order closure model is also compared withthat of lower level turbulence models, including the eddy-viscositymodel and algebraic Reynolds stress models. It is concluded that thepresent closure is a considerable improvement over the other modelsin representing various physical effects in flow over hills. Thefeasibility of running a finite-difference numerical simulationincorporating a full second-order closure model on an IBM workstationis also demonstrated.  相似文献   

18.
Measurements of mean wind flow and turbulence parameters have been made over Cooper's Ridge, a 115 m high elongated ridge with low surface roughness. This paper describes measurements of the streamwise and vertical variations in the mean field for a variety of atmospheric stability conditions. In near-neutral conditions, the normalised speedup over the ridge compares well with measurements from Askervein (Mickleet al., 1988). The near-neutral results are also compared to an analytical flow model based on that of Huntet al. (1988a). Measured streamwise variations show less deceleration at the foot of the hill and slightly more acceleration at the crest of the hill than does the model. In non-neutral conditions, the speedup over the ridge reduces slightly in unstable conditions and increases by up to a factor of two in stable conditions. The model is modified to allow boundary-layer stability to change the upwind wind profile and the depths of the inner and middle layers. Such a modification is shown to describe the observations of speedup well in unstable and weakly stable conditions but to overestimate the speedup in moderate to strongly stable conditions. This disagreement can be traced to the model's overestimation of the upstream scaling velocity at the height of the middle layer through its use of a stable wind profile form which has greater shear than that of the observed profiles, in possible combination with the three-dimensionality of the ridge which would allow enhanced flow around, rather than over, the feature in more stable conditions.  相似文献   

19.
The Ekman boundary layer over orography: An analysis of vertical motion   总被引:2,自引:0,他引:2  
A model of the planetary boundary layer is used to determine the field of vertical motion over large-scale orography. This model represents Ekman boundary-layer dynamics modified by the inclusion of accelerations of the geostrophic wind under the geostrophic momentum approximation. The orography is represented by a circular mountain. The inviscid solution is provided by the sum of a constant translation and a steady, uniform potential vorticity, anticyclonic vortex. The boundary-layer solution vanishes on the mountain, but is matched to the inviscid solution as the top of the boundary layer is approached. The vertical velocity field at the top of the boundary layer is determined by integration of the continuity equation. The field of motion is largely determined by descent from above into the anticyclonic circulation, as in the classical Ekman model. Contributions that arise from the inclusion of accelerations are associated with boundary-layer advection and ageostrophic divergence that produce vorticity tendencies. Finally, the boundary-layer vertical motion is shown to be comparable in magnitude to the vertical motion forced by inviscid flow over the orography, although the distributions of each are significantly different. Effects of mountain asymmetry and a changing pressure field, that can be treated more fully by numerical model simulations, are not considered in the present study.On leave at the University of Colorado, 1990.  相似文献   

20.
Summary Two-dimensional simulations of plume advection over a hill are presented. The calculations are carried out in two steps. In a first step the flow is calculated by a mesoscale-model. This model solves the equations of motions in a cartesian grid. In a second step, the transport equation is solved using the calculated flow-parameters of the mesoscale model. The numerical formulation of the transport-equation was presented in Part I. This calculation is performed in a coordinatetransformed grid. The merit of this combination is discussed in the paper.With 11 Figures  相似文献   

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

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