对流边界层中地面源的铅直扩散模拟

对流边界层(CBL)中高架源浓度轴线下倾而地面源浓度轴线上升。本文分析了形成这种相反的几何图象的机理，从而认为:只要能够恰当模拟CBL中铅直湍流结构的特征，一个粒子随机扩散模式应既能模拟高架源的扩散，也能模拟地面源的扩散。因此，用作者早先模拟高架源的模式模拟了地面源的铅直扩散，同样获得较好的结果，进一步证实了模式的有效性。本文还应用Thomson准则检验了模式，讨论了它在理论上的合理性和适用范围。     

对流边界层中的一个双高斯型PDF扩散模式

对流边界层（ＣＢＬ）中的污染扩散是非高斯型的。本文在下列三个假设下建立了双高斯型几率分布函数（ＰＤＦ）模式：１．对流边界层任一确定高度铅直速度Ｗ的几率分布函数ｐｗ由两个高斯分布迭加而成；２．从污染源释放的粒子具有源高的铅直速度几率分布，且其轨迹是线性的；３．粒子在地面的反射为全反射，在混合层高度Ｚｉ为全反射或有部分吸收。然后分析了三个高度上铅直速度Ｗ的一些统计特征量，比较了由ＰＤＦ模式计算的横风向积分浓度和Ｌａｍｂ的数值模拟［１－３］，Ｄｅａｒｄｏｒｆｆ的水槽模拟［４－６］结果，并用美国ＣＯＮＤＯＲＳ计划的外场试验资料［７］对ＰＤＦ模式进行了验证，结果均相当一致。     

一个直接考虑大涡、小涡作用的大气对流边界层模式

本文建立了一个大气对流边界层模式,构造此模式时,假设对流边界层中水平动量和热量在垂直方向上的湍流输送、扩散是由大涡与小涡组成的系统完成的。文中,用此模式作了一实例模拟,计算结果表明:模拟的结果与实测的结果相当一致。     

一个对流边界层大涡模式的建立与调试

本文介绍一个适合于对流边界层的大涡模式的建立及其调试结果。该大涡模式建立过程中注重于计算的节省，同时也强调原理与方法的简单和合理性。模式的调试表明，对于平坦均一地形的情况，模拟可以获得合理的结果。调试同时显示了模式对较低水平分辨率的适用条件，以及模式应用于模拟较大水平范围问题的可能性。     

对流边界层中粒子随机扩散模式和高斯模式的比较

用粒子随机扩散模式和高斯模式模拟的对应３种源高的横风向各分浓度的结果与水槽试验的结果进行比较，发现粒子随机扩散模式能较好地模拟出水槽试验结果，说明该模式能用来模拟对流边界层的扩散，而高斯模式不能完全反映对流边界层的扩散。比较两种模式的模拟结果，发现粒子随机扩散模式由于考虑了对流边界中由卷流和热泡引起的垂直方向上湍流的非均一性，比高斯模式能较好的模拟出污染物在垂直方向上充分混合较快、水平方向上扩散范     

用大涡模式研究对流边界层湍流

利用三维对流边界层的大涡模拟模式研究对流边界层湍流统计特征量,模拟结果与前人的同类模拟工作及实验观测结果相比较,得到了一致的结论.       

对流边界层中公路线源扩散的Monte-Carlo法模拟

开发了一个Ｍｏｎｔｅ－Ｃａｒｌｏ大气扩散模式且引入了一种检验随机数质量的填充二维点格的方法。通过与水槽试验资料对比,说明该模式具有较好的模拟效果。用该模式模拟公路线源的大气扩散,结果表明,高架路和一般公路的污染在地面浓度上的差异较大,而在垂直方向上差异较小。高架路引起的地面浓度值和扩散范围都较一般公路的大。垂直方向上尤其是２０米以上二者基本没有差异,浓度较大值都在５０～１００米以下。结合箱模式的思     

对流边界层中泡状结构的大涡模拟研究

本文提出“连续垂直运动区”的概念，对边界层中的对流泡状运动进行了重新定义。由此将对流边界层划分为上升泡区、下沉泡区和环境气流区，并利用大涡模拟提供的对流边界层数值模拟结果研究这种泡状运动的结构。研究结果较好地解释了一些有关泡状运动结构方面的似乎相互矛盾的观测事实，并揭示了对流泡状运动的一些新的特征。     

日间对流边界层中的非局地动量混合

日间对流边界层最显著的结构特征是在热力作用下所形成的组织化对流。与小尺度湍涡不同的是,组织化对流具有边界层尺度的垂直相干性,可实现垂直贯穿边界层的非局地物质和能量传输。本文针对对流边界层中的动量混合,探究组织化对流对动量输送的贡献。以高精度大涡模拟数据为研究资料,通过傅里叶变换、本征正交分解和经验模态分解3种滤波方法,分离组织化对流和背景湍涡,计算与两者相关的非局地和局地动量通量,发现与组织化对流相关的非局地动量通量是总通量的重要组成部分,并主导混合层中的垂直动量输送。而后,基于协谱和相位谱分析,探究组织化对流的空间结构对动量传输的影响,发现在热力主导的不稳定环境中,单体型环流结构对动量的传输效率较低。而在风切较强的近中性环境中,滚涡型组织化结构可使垂直和水平流向扰动速度的相位差减小,从而提升动量传输效率。研究结果表明,边界层方案需要包含非局地动量通量项,其参数化应考虑整体稳定度对传输效率的影响。     

论边界层中的大气扩散PDF模式

基于大气扩散Ｋ理论，用作为风速脉动均方差和拉氏时间尺度函数的湍流交换系数，得到了直接利用风速脉动几率密度而不用扩散参数的大气扩散ＰＤＦ模式。分别研究了对流边界层上升气流区与下降区垂直速度的统计特征，求得双正态ＰＤＦ模式。在给定ＣＢＬ自身参数如对流特征速度ｗ＊，顶高ｈｉ和源高度上的平均风速时，该模式计算出的无量纲浓度分布与室内外测试结果一致。     

基于大涡模拟评估GRAPES模式对对流边界层的模拟性能

为检验GRAPES半拉格朗日动力框架在大涡尺度上的模拟性能,为未来发展千米及其以下高分辨尺度的数值模式奠定基础,并构造GRAPES大涡模式以检验和发展边界层湍流参数化提供科学工具。通过在GRAPES模式中加入Smagorinsky-Lilly小尺度湍涡参数化,并将模式分辨率提高至50 m,构建GRAPES大涡模式(GRAPES_LES),以便分析GRAPES模式在大涡尺度上的适用性。同时利用广泛应用的已有大涡模式UCLA_LES作为参考,通过对干对流边界层湍流的模拟分析及与UCLA_LES模拟结果的对比,得出如下主要结论:GRAPES半拉格朗日动力框架能够模拟出与已有的大涡模式相似的边界层湍流特征;同时,通过分析也证明GRAPES存在由于采用半拉格朗日平流计算而带来过度耗散的问题:当使用相同的滤波尺度(Smagorinsky 常数)时,GRAPES_LES模拟出的速度场更为平滑,小尺度湍流结构过于光滑,通过对湍流能量的能谱分析更清楚地表明了这一点。进一步,对不同的Smagorinsky常数(对应不同的滤波尺度)进行了敏感性试验,表明可以通过改变滤波尺度,有效地缓解半拉格朗日框架隐含的耗散问题,得到更接近UCLA_LES所模拟的湍流特征。     

A numerical study of daily transitions in the convective boundary layer

We examine daily (morning–afternoon) transitions in the atmospheric boundary layer based on large-eddy simulations. Under consideration are the effects of the stratification at the top of the mixed layer and of the wind shear. The results describe the transitory behaviour of temperature and wind velocity, their second moments, the boundary-layer height Z _m (defined by the maximum of the potential temperature gradient) and its standard deviation σ _m , the mixed-layer height z _i (defined by the minimum of the potential temperature flux), entrainment velocity W _e, and the entrainment flux H _i . The entrainment flux and the entrainment velocity are found to lag slightly in time with respect to the surface temperature flux. The simulations imply that the atmospheric values of velocity variances, measured at various instants during the daytime, and normalized in terms of the actual convective scale w_*, are not expected to collapse to a single curve, but to produce a significant scatter of observational points. The measured values of the temperature variance, normalized in terms of the actual convective scale Θ_*, are expected to form a single curve in the mixed layer, and to exhibit a considerable scatter in the interfacial layer.     

Turbulent mixing of reactive gases in the convective boundary layer

We study the interactions of chemistry and turbulent mixing of tracersin the convective boundary layer with a second-order closure model,including higher order chemistry terms. In order to limit the number of predictive equations we prescribe the profiles for ¯w¯, ¯w¯ ¯ and the lengthscale l. However, for model validation we treat temperature and humidity asinert tracers, and compare the results with profiles observed during theAir Mass Transformation Experiment, and with similarity expressions for thesurface layer. We find good agreement of the mean profiles, but the (co-)variances are slightly underpredicted. Furthermore, the model usesdiagnostic equations expressing third moments of concentration in terms ofsecond moments and their vertical derivatives. They are compared withlarge-eddy model results, showing good agreement and, therefore, thesimplifications are justified. The model is applied to the transport of two gases subject to one bimolecular reaction. The importance of concentration correlations on themean transformation rate is studied. For two gases diffusing in oppositedirections we find for moderate and fast chemistry a 50% and90% decreased transformation rate due to the negatively correlatedconcentrations. These values are similar to large-eddy results of Schumannand Sykes et al. For two bottom-up tracers we find that the covariance ofboth reactive species is either positive or negative, increasing or reducingthe effective transformation rate depending on the Damköhler number (the ratio of the turbulent and the chemistry timescale). A significantdirect influence of chemistry on the flux divergence is found in bothcases. According to the model the effective transport to mid-levels of theboundary layer is increased when two reactive tracers diffuse in oppositedirections, and decreased in the case of two bottom-up tracers.     

Non-local mixing of momentum in the convective boundary layer

Recently Frech and Mahrt proposed a closure scheme which includes alarge-scale stress term to represent the effects of non-local momentummixing in the convective boundary layer. Here large-eddy simulation (LES)datasets are used to evaluate the performance of this scheme across a rangeof stabilities between neutral and highly convective conditions, and as afunction of baroclinity. Generally the inclusion of the non-local term inthe closure model leads to results in better agreement with LES, althoughsome modifications to the model formulation are suggested.     

城市建筑动力学效应对对流边界层影响的敏感性试验

本文将大涡模拟应用于城市对流边界层(CBL)湍流结构和流场特征的研究,在大涡模式中,拖曳系数取与建筑物高度及建筑物高度标准差有关的表达式以考虑次网格建筑物对风速和湍流动能(TKE)的面积平均影响.模拟结果表明,由于城市建筑物对气流的拖曳作用,使建筑物冠层及整个CBL内风速大幅度减小,城市冠层内部风速减小尤为明显,在夹卷层内,风速有一明显的跃变.在边界层中部对流运动已经发展成为较强的热泡,城市建筑物的动力学效应使热泡的水平尺度增大,CBL内平均上升气流速度和下沉气流速度减小,同时使CBL中上升气流所占比例比平坦地面增大.城市建筑物使CBL低层热通量、动量通量、速度方差和位温方差明显增大,但对近地层高度以上的湍流量影响不大.     

风切变对边界层对流影响的大涡模拟研究

利用"西北干旱区陆-气相互作用野外观测实验"加密观测期间在敦煌站的观测资料以及大涡模式,模拟了对流边界层的发展,以及示踪物从混合层向残留层传输的时空变化。模拟的对流边界层的结构及演变特征与实测结果基本一致。进一步通过有风切变和无风切变的敏感性数值试验,研究了风切变对垂直速度、位温和示踪物浓度的水平分布以及示踪物传输高度的影响。研究结果表明,在有风切变的试验中(甚至风切变仅存在于近地层中),对流边界层的增长加强,而且示踪物被传输的高度也较高。与浮力驱动的对流边界层相比,由浮力和风切变共同驱动的边界层中上升气流较弱而下沉气流较强,但前者的上升气流与下沉气流的分布在垂直方向上更为倾斜。由于夹卷作用的增强,浮力和风切变共同驱动的对流边界层较浮力驱动的对流边界层暖。在夹卷层,浮力和风切变共同驱动的边界层对流的上升气流和下沉气流都比浮力驱动的边界层对流中的强,而且垂直速度的概率密度函数分布也较对称,其位温和示踪物浓度的概率密度函数分布也比浮力驱动的边界层中的平直。对湍流动能收支的分析也表明风切变对湍流动能有重要影响,尤其对夹卷层中的湍流动能切变产生项影响较大。示踪物浓度的概率密度函数垂直分布显示,浮力驱动的边界层中示踪物浓度随高度变化较小,而浮力和风切变共同驱动的边界层中示踪物浓度随高度递减,但是示踪物传输的高度比较高。     

Negative water vapour skewness and dry tongues in the convective boundary layer: observations and large-eddy simulation budget analysis

This study focuses on the intrusion of dry air into the convective boundary layer (CBL) originating from the top of the CBL. Aircraft in-situ measurements from the IHOP_2002 field campaign indicate a prevalence of negative skewness of the water vapour distribution within the growing daytime CBL over land. This negative skewness is interpreted according to large-eddy simulations (LES) as the result of descending dry downdrafts originating from above the mixed layer. LES are used to determine the statistical properties of these intrusions: their size and thermodynamical characteristics. A conditional sampling analysis demonstrates their significance in the retrieval of moisture variances and fluxes. The rapid CBL growth explains why greater negative skewness is observed during the growing phase: the large amounts of dry air that are quickly incorporated into the CBL prevent a full homogenisation by turbulent mixing. The boundary-layer warming in this phase also plays a role in the acquisition of negative buoyancy for these dry tongues, and thus possibly explains their kinematics in the lower CBL. Budget analysis helps to identify the processes responsible for the negative skewness. This budget study underlines the main role of turbulent transport, which distributes the skewness produced at the top or the bottom of the CBL into the interior of the CBL. The dry tongues contribute significantly to this turbulent transport.     

A laboratory study of the turbulent velocity characteristics in the convective boundary layer

Based on the measurement of the velocity field in the convective boundary layer （CBL） in a convection water tank with the particle image velocimetry （PIV） technique, this paper studies the characteristics of the CBL turbulent velocity in a modified convection tank. The experiment results show that the velocity distribution in the mixed layer clearly possesses the characteristics of the CBL thermals, and the turbulent eddies can be seen obviously. The comparison of the vertical distribution of the turbulent velocity variables indicates that the modeling in the new tank is better than in the old one. The experiment data show that the thermal＇s motion in the entrainment zone sometimes fluctuates obviously due to the intermittence of turbulence. Analyses show that this fluctuation can influence the agreement of the measurement data with the parameterization scheme, in which the convective Richardson number is used to characterize the entrainment zone depth. The normalized square velocity wi^2/w＊^2. at the top of the mixed layer seems to be time-dependent, and has a decreasing trend during the experiments. This implies that the vertical turbulent velocity at the top of the mixed layer may not be proportional to the convective velocity （w＊）.