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

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

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

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

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

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

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

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

对流边界层的大涡模拟研究

本文建立了一个均匀平坦地面上对流边界层的大涡模式,模式考虑了水汽,采用了考虑浮力和固壁影响订正的一阶闭合。并用所建模式进行了由热扰动发展的对流边界层的模拟及其对地表热状况变化响应的初步探讨性模拟工作。通过模拟认为,模式较好地反映了对流边界层的主要结构。     

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

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

黄河源区鄂陵湖湖面和湖边草地对流边界层湍流结构特征的大涡模拟研究

为研究黄河源区边界层湍流特征及其对物质和能量输送的影响,本文首次采用大涡模拟的方法,对比分析了黄河源区两种不同下垫面上(鄂陵湖和湖边草地)对流边界层(CBL)中精细的湍流结构特征.使用资料为2012年夏季黄河源区鄂陵湖流域野外观测实验的GPS探空资料、涡动相关观测资料.分析表明,模拟的黄河源区草地和湖上CBL的平均结构...     

一个对流边界层中的随机扩散模式

通过对对流边界层(CBL)湍流结构的分析,首次提出用两种不同尺度的湍流模拟CBL中的铅直扩散.在此基础上发展了一个随机扩散模式,并用它模拟了典型对流条件下两种高架连续点源的扩散.与Willis水槽模拟和Lamb等人数值模拟以及CONDORS计划外场试验的结果的比较表明,本模式能成功地模拟CBL中的横向积分浓度.与其它数值模式相比,还具有输入参数少、计算量极小和更加简单实用的优点.     

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

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

海面对流边界层垂直结构的数值模拟

冬季冷空气爆发进入暖海面时,海面边界层属与典型的对流边界层。此情况下,一般而言梯度输送理论(即一阶K闭合方法)是无能为力的。Wgnsaard(1984)提出了被动保守标量场的底直向上(bottom-up)和顶直向下(top-down)湍流交换模型。本文将它推广应用到风场相温度场,并用它模拟了冷气团入海期间海面边界层的结构。与实测资料比较,两者能吻合较好。     

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.     

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

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

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.     

On the turbulence structure in the stable boundary layer over the Greenland ice sheet

Turbulence measurements from a 30 m tower in the stably stratifiedboundary layer over the Greenland ice sheet are analyzed. The observationsinclude profile and eddy-correlation measurements at various levels. Atfirst, the analysis of the turbulence data from the lowest level (2 m aboveground) shows that the linear form of the non-dimensional wind profile(m) is in good agreement with the observations for z/L <0.4, whereL represents the Obukhov length. A linear regression yieldsm=1+5.8z/L. The non-dimensional temperature profile (h) at the2m level shows no tendency to increase with increasing stability. The datafrom the upper levels of the tower are analyzed in terms of both localscaling and surface-layer scaling. The m and the h values show atendency to level off at large stability (z/>0.4) where represents the local Obukhov length. Hence, the linear form of the functions is no longer appropriate under such conditions. The bestcorrespondence to the data can be achieved when using the expression ofBeljaars and Holtslag for m and h. The vertical profiles of theturbulent fluxes, the wind velocity variances and temperature variance arealso determined. The momentum flux profile and the profiles of wind speedvariances are in general agreement with other observations if a welldeveloped low-level wind maximum occurs, and the height of this maximum isused as a height scale.     

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＊）.     

对流边界层湍流通量及逆梯度输送参数化分析

"K"理论是众多气象预报模式中运用最广泛的湍流参数化方案之一,但无法解释"逆梯度"的输送,必须进行修正。最具代表性的修正方案有三种:方案Ⅰ(Deardroff方案)、方案Ⅱ(Holtslag和Moeng方案)和方案Ⅲ(刘烽方案)。本文利用香河的边界层观测资料对上述三种方案进行验证和比较,发现方案Ⅰ的结果在整个对流边界层(Convective Boundary Layer,CBL)呈系统性偏低,与观测不符;方案Ⅱ在CBL中上部能够再现逆梯度输送现象,基本能给出合理的湍流通量垂直分布,但在CBL的下部和上部与观测不符;方案Ⅲ的逆梯度项与高度有关,并在CBL中部达到最大,而其他两个方案中逆梯度项随高度不变。该方案不但在CBL中上部与方案Ⅱ的结果一致,并能合理表达整个CBL内的湍流通量分布,更接近观测结果。