一次冷锋过境后的海风三维结构数值模拟

为研究大尺度系统风对海风的影响以及海风三维结构特征,利用山东省123个地面自动站资料、青岛地区三十多个内陆及沿海、海岛观测站以及奥帆赛场3个浮标站资料,对2006年8月21日青岛一次海风个例进行了分析,并利用美国俄克拉荷马大学风暴分析预测中心开发的ARPS(the Advanced Regional Prediction System)模式,对海风过程进行了数值模拟研究。结果发现:在较强的离岸风背景下,当内陆气温高于海面气温2℃左右时,海风也可以发生。海风首先在海岸线附近的海上开始,发展的同时向内陆及远海地区推进。海风低层环流很浅,主要位于500 m以下。在较强的偏北离岸风下,海风向内陆推进的距离很短。偏北的大尺度系统风由于渤海冷下垫面的影响,不利于青岛海风的维持。海风开始时,在1500～2500 m高度处同时有反环流出现,但直到傍晚前后,海风的垂直环流圈才发展得比较清晰,其高度也更接近地面。海风消亡后,高层的垂直环流圈及反环流维持3 h左右才逐渐消亡。     

地表通量输送对飑线过程影响的数值模拟研究

使用水平分辨率20km的RSM-RANAL再分析资料和WRF模式,对2004年7月12日影响上海地区的一次飑线过程进行分析和数值模拟,结果表明:(1)雷暴发生区域的西北侧约300km处有一冷锋,地面风场辐合;中低层为一致西南急流,雷暴发生区域上空300hPa位于西风急流右后侧辐散区中;地面存在明显的温度梯度。(2)雷暴发生前浙皖交界山区处的CAPE增至1000J/kg,CIN减至30J/kg以下;杭州湾附近的CAPE值大于2000J/kg。这些地方之后均有强对流单体生成,并成为飑线的一部分。(3)WRF模式对这次飑线过程的主要特征模拟较好,包括地面风场和飑线的结构。(4)敏感性试验指出陆面过程对本次飑线过程的形成发展起了重要作用。一方面通过潜热通量输送增大边界层湿度;另一方面通过感热通量输送改变了边界层的层结结构,使得低层辐合(或抬升)比较容易释放不稳定。在强对流天气的预报中地表条件不可忽略。     

Application of Dynamic Subgrid-scale Models for Large-eddy Simulation of the Daytime Convective Boundary Layer over Heterogeneous Surfaces

The sensitivity of large-eddy simulation (LES) to the representation of subgrid-scale (SGS) processes is explored for the case of the convective boundary layer (CBL) developing over surfaces with varying degrees of spatial heterogeneity. Three representations of SGS processes are explored: the traditional constant Smagorinsky–Lilly model and two other dynamic models with Lagrangian averaging approaches to calculate the Smagorinsky coefficient (C _S ) and SGS Prandtl number (Pr). With initial data based roughly on the observed meteorology, simulations of daytime CBL growth are performed over surfaces with characteristics (i.e. fluxes and roughness) ranging from homogeneous, to striped heterogeneity, to a realistic representation of heterogeneity as derived from a recent field study. In both idealized tests and the realistic case, SGS sensitivities are mostly manifest near the surface and entrainment zone. However, unlike simulations over complex domains or under neutral or stable conditions, these differences for the CBL simulation, where large eddies dominate, are not significant enough to distinguish the performance of the different SGS models, irrespective of surface heterogeneity.     

Subgrid-Scale Dynamics of Water Vapour, Heat, and Momentum over a Lake

We examine the dynamics of turbulence subgrid (or sub-filter) scales over a lake surface and the implications for large-eddy simulations (LES) of the atmospheric boundary layer. The analysis is based on measurements obtained during the Lake-Atmosphere Turbulent EXchange (LATEX) field campaign (August–October, 2006) over Lake Geneva, Switzerland. Wind velocity, temperature and humidity profiles were measured at 20 Hz using a vertical array of four sonic anemometers and open-path gas analyzers. The results indicate that the observed subgrid-scale statistics are very similar to those observed over land surfaces, suggesting that the effect of the lake waves on surface-layer turbulence during LATEX is small. The measurements allowed, for the first time, the study of subgrid-scale turbulent transport of water vapour, which is found to be well correlated with the transport of heat, suggesting that the subgrid-scale modelling of the two scalars may be coupled to save computational resources during LES.     

Footprints in Homogeneously and Heterogeneously Driven Boundary Layers Derived from a Lagrangian Stochastic Particle Model Embedded into Large-Eddy Simulation

A Lagrangian stochastic (LS) model, which is embedded into a parallelised large-eddy simulation (LES) model, is used for dispersion and footprint evaluations. For the first time an online coupling between LES and LS models is applied. The new model reproduces concentration patterns, which were obtained in prior studies, provided that subgrid-scale turbulence is included in the LS model. Comparisons with prior studies show that the model evaluates footprints successfully. Streamwise dispersion leads to footprint maxima that are situated less far upstream than previously reported. Negative flux footprints are detected in the convective boundary layer (CBL). The wide range of applicability of the model is shown by applying it under neutral and stable stratification. It is pointed out that the turning of the wind direction with height leads to a considerable dependency of source areas on height. First results of an application to a heterogeneously heated CBL are presented, which emphasize that footprints are severely affected by the inhomogeneity.     

Numerical Study of the Response of an Atmospheric Surface Layer to a Spatially Nonuniform Plant Canopy

High-accuracy large-eddy simulations of neutral atmospheric surface-layer flow over a gapped plant canopy strip have been performed. Subgrid-scale (SGS) motions are parameterized by the Sagaut mixed length SGS model, with a modification to compute the SGS characteristic length self-adaptively. Shaw’s plant canopy model, taking the vertical variation of leaf area density into account, is applied to study the response of the atmospheric surface layer to the gapped dense forest strip. Differences in the region far away from the gap and in the middle of the gap are investigated, according to the instantaneous velocity magnitude, the zero-plane displacement, the potential temperature and the streamlines. The large-scale vortex structure, in the form of a roll vortex, is revealed in the region far away from the gap. The nonuniform spatial distribution of plants appears to cause the formation of the coherent structure. The roll vortex starts in the wake of the canopy, and results in strong fluctuations throughout the entire canopy region. Wind sweeps and ejections in the plant canopy are also attributed to the large vortex structure.     

A Hybrid Spectral/Finite-Volume Algorithm for Large-Eddy Simulation of Scalars in the Atmospheric Boundary Layer

Pseudospectral methods are frequently used in the horizontal directions in large-eddy simulation of atmospheric flows. However, the same approach often creates unphysical oscillations for scalar fields if there are horizontal heterogeneities in the sources and/or sinks, as is usual in air pollution problems. A hybrid approach is developed to combine the use of pseudospectral representation of the velocity field and bounded finite-volumes for the scalar concentration. An interpolation scheme that yields a divergence-free interpolated velocity field is derived and implemented, and its importance is illustrated by two sample applications.     

Large-Eddy Simulation of Flows over Random Urban-like Obstacles

Further to our previous large-eddy simulation (LES) of flow over a staggered array of uniform cubes, a simulation of flow over random urban-like obstacles is presented. To gain a deeper insight into the effects of randomness in the obstacle topology, the current results, e.g. spatially-averaged mean velocity, Reynolds stresses, turbulence kinetic energy and dispersive stresses, are compared with our previous LES data and direct numerical simulation data of flow over uniform cubes. Significantly different features in the turbulence statistics are observed within and immediately above the canopy, although there are some similarities in the spatially-averaged statistics. It is also found that the relatively high pressures on the tallest buildings generate contributions to the total surface drag that are far in excess of their proportionate frontal area within the array. Details of the turbulence characteristics (like the stress anisotropy) are compared with those in regular roughness arrays and attempts to find some generality in the turbulence statistics within the canopy region are discussed.     

蒙特卡罗不确定性分析在O3模拟中的初步应用

利用蒙特卡罗不确定性分析方法,分析了嵌套网格空气质量模式系统(NAQPMS/IAP)中154个模式输入变量不确定性对臭氧模拟的影响,量化了模式在北京奥运会期间臭氧模拟的不确定性,并确定出了主要不确定性因子。结果表明:(1)在奥运会期间(2008年8月8日～2008年8月24日),北京城区臭氧模拟的平均不确定性为19ppb,不确定性存在明显的日变化特征,白天不确定性大,夜间不确定性小。(2)在白天,北京城区近地层臭氧模拟最重要的不确定性来源是局地前体物排放,其次是NO2光解系数、风向、北京周边前体物排放和垂直扩散系数。另外,地面约150m以上,对臭氧模拟影响最大的不确定性因子是风向和北京周边前体物排放;夜间,北京城区近地层臭氧模拟的最大不确定性来源是局地NOx排放和垂直扩散系数。     

2008年北京奥运会开幕式化学危害预测数值模拟

针对2008年北京奥运会安保,建立了一个复杂的系统,包括两个气象预报模式MM5和RAMS6.0以及一个复杂地形上有毒云团的扩散模式CDM。开发了MM5和RAMS6.0接口模块,发展了基于RAMS6.0气象场预报结果的CDM,实现了对未来36h内的精细气象场和扩散场进行模拟,气象场预报的分辨率为500～1000m,扩散场的分辨率可达到50～200m。基于T213数据,应用MM5进行了水平分辨率为3km的气象场初步预报,将该预报结果转化为RAMS6.0识别的数据格式,应用RAMS6.0对气象场进行分辨率为500～1000m的精细预报,得出北京奥运会期间的风场与湍流量的预报结果。基于该气象场,根据想定,应用CDM对北京奥运会国家体育场"鸟巢"附近发生化学事件时进行了模拟研究,得出危害范围、危害等级、危害开始时间、危害持续时间等时空分布,为奥运会开幕式化学危害应急提供技术支持。