非均匀网格的过渡湍流理论及其大气边界层数值模拟中的应用

本文将Ｓｕｌｌ提出的均匀网格下的过渡湍流理论推广到非均匀网格情形，推广的非均匀网格的过渡湍流理论满足Ｓｔｕｌｌ提出的对过渡矩阵系数的要求并具有清晰的物理意义。然后，将非均匀网络的过滤湍流理论应用于一维大气边界层数值模式中，对Ｗａｎｇａｒａ资料进行了模拟，并与均匀网格情形进行了对比。     

高阶矩模式模拟对流边界层的可行性研究

本文建立了二阶矩闭合的大气边界层数值模式,并模拟了Wangara资料第33天09:00—18:00时对流边界层的发展过程。计算结果与实测资料比较吻合,同时利用该模式研究了高阶矩模式模拟对流边界层的可行性,结果指出:二阶矩方程中的湍流扩散项对模拟对流边界层内的挟卷和反梯度现象起主要作用。而局地时间变化项仅影响对流边界层内的边界湍流扩散项。因此,假如处理好二阶矩方程中的湍流扩散项,高阶矩模式和大涡模式一样可以很好地模拟对流边界层的发生、发展。     

中尺度数值模拟中的边界层多尺度湍流参数化方案

该文在多尺度湍流理论的研究成果基础上, 将边界层湍流风谱与平均量的梯度相联系, 建立了边界层多尺度湍流参数化子模式, 之后放入MM5模式中进行了个例模拟研究, 并与MM5模式附带的M RF边界层参数化、Blackadar高分辨率边界层参数化的模拟结果进行比较和分析。结果表明, 多尺度湍流理论能够反映出实际大气边界层中热量垂直输送的规律, 将其用于中尺度数值预报模式的边界层物理过程参数化是可行的; 多尺度湍流参数化在地表层和边界层内各个层次上都着重考虑含能量最大的涡的作用以及水平热力不均匀性的影响, 因此在地形和下垫面比较复杂的区域, 对中尺度天气系统的模拟有进一步发展的前景。     

湍流动能闭合方法在中尺度模式中的应用

湍流动能闭合方法是近年来发展起来的用以模拟大气边界层的一种方法,本文对其做了简化,仅在边界层内使用这一方法,进一步减少计算量,使之更适用于中尺度数值模拟。对一维旺加拉（Wangara）资料进行的试验表明,这种方法保留了边界层计算的精度,又能节省机时。然后将这种方法应用于三维中尺度钟形山地形的模拟,描述了湍流动能的分布,以及湍流活动对山后回流区面积和强度的影响。     

一维大气边界层二阶闭合的有限元数值模式 I. 对流边界层模拟

利用有限元法用于二阶湍流闭合大气对流边界层模式,并模拟了Wangara试验33天的对流边界层的发展过程,计算结果表明,该模式能真实的模拟出对流边界层的发展过程,并且准确地反映对流边界层中湍流扩散和输送过程。     

尺度自适应大气边界层参数化改进及其对一次海雾的数值模拟研究

大气边界层湍流运动是地球大气运动最重要的能量输送过程之一。当数值模式分辨率接近活跃含能湍涡长度尺度时，湍流运动被部分解析，被称为“灰色区域”，传统的边界层方案不适合此时模式湍流问题的描述。为了提高模式边界层方案在包括“灰色区域”的不同网格尺度上的描述能力，适应不同分辨率模式的需要，在雷诺平均湍流理论基础上，修正Mellor-Yamada-Nakanishi-Niino (MYNN)方案湍流长度尺度参数和非局地湍流的参数表达，改进湿度和温度在“灰色区域”的湍流输送参数化及对网格尺度的自适应能力。利用改进的MYNN尺度自适应方案，分别采用3 km和1 km、1.5 km和0.5 km分辨率单向嵌套网格WRF中尺度模式，对2014年2月26日的一次黄海海雾过程进行模拟试验，检验不同分辨率下改进后的MYNN大气边界层参数化方案的合理性和对海雾的模拟效果。尺度自适应MYNN大气边界层参数化方案在千米级网格尺度上获得稳定、合理的湍流垂直输送计算结果。参照雾区卫星云图，不同分辨率模式低层云水混合比模拟结果具有稳定表现，模拟的雾区分布和温度、湿度等物理量结构都较好地再现了再分析“观测事实”，初步表明该参数化方案有较高的网格尺度自适应能力。      

利用尺度自适应大气边界层湍流参数化方案对一次陆地浓雾的数值模拟

为提高大气数值模式的模拟能力,改进大气边界层水汽、热量湍流输送计算和大雾天气的模拟效果,选用WRF三维非静力模式,采用具有局地和非局地垂直湍流尺度自适应计算能力的MYNN_SA参数化方案,对2017年12月28—29日我国华北—江淮地区的大范围浓雾过程进行了数值模拟研究,探讨从中尺度到灰区尺度分辨率范围,模式的尺度自适应大气边界层湍流参数化方案对稳定大气边界层发展、湍流输送和大雾天气模拟的影响。利用中国地面气象站观测资料和ERA5再分析数据,在接近灰区尺度的网格分辨率上,利用尺度自适应大气边界层湍流MYNN_SA参数化方案较之中尺度参数化MYNN方案,可明显改善次网格湍流输送计算,以及陆地浓雾的强度、空间分布和时间演变特征,可更精确地模拟云水混合比、逆温层和雾区的垂直结构。     

深凹地形边界层风场与湍流场结构及扩散规律的数值研究

建立了一个模拟小尺度复杂地形边界层平均流场和湍流场结构分布特征的三维非静力细网格高阶矩湍流闭合模式,并将此模式与一个拉格朗日粒子随机游动大气扩散数值模式相联结,以边界层气象模式的输出作为扩散模式的输入,构成了一个大气扩和模拟系统,成功地模拟了某地实际深凹露天矿区和理想凹坑地形两种情况下的平均流场、雷诺应力和湍流通量场,分析了其分布特征,并模拟了矿坑内污染物的散布规律     

局地环流与大气边界层湍流相互作用的数值研究

本文建立了一个二维定常大气中尺度数值模式，并用该模式讨论了由下垫面粗糙分布非均匀（局地动力强迫）和温度分布非均匀（局地热力强迫）而产生的局地大气环流与大气边界层湍流的相互作用。结果指出：下垫面存在温度分布非均匀时运动方程中的湍流交换项与水平气压梯项一样可促使局地热力环境的形成，由粗糙度分布非均匀强迫产生的局地环流与由下垫面温度分布非均匀强迫产和的局地热力环流间的非线笥相互作用是通过湍流交换实现的。     

青岛地区边界层结构的数值模拟

本文建立了一个非静力的三维细网格边界层模式，对青岛地区复杂下垫面条件下的边界层结构和湍流特征作了数值模拟，模式采用能量闭合方案，舍弃了静力近似，以提高空间分辨率和模式精度。模拟以实测资料为初始输入，对该地区风场和湍流场作了较细致的模拟分析，与观测资料相比，结果合理地一致。结果表明了陡峭地形和不规则海岸线对局地风场和湍流场影响很大。     

复杂下垫面模拟域大气边界层非局地闭合模拟研究

在非局地闭合大气边界层模拟研究基础上,选取黑河实验研究资料（HEIFE）,利用二维模式对绿洲和沙漠交界地区的内边界层结构进行了模拟研究,模式考虑了两种下垫面条件下动力、热力特性的不同,并加入了湿度及水汽通量的预报。模拟结果表明:二维非局地闭合模式能得到合理的TIBL结构的基本特征,包括平均场、湍流场等,并能模拟给出一些由下垫面特性形成的大气现象,如沙漠地区逆湿现象的生消过程、绿洲的冷岛效应、绿洲与沙漠因热力差异而形成的局部环流形势等细致结构,符合沙漠站白天观测事实的典型特征。     

Dynamical Initialization of Boundary-Layer Data

A procedure for the dynamic initialization of wind and temperature fields within the atmospheric boundary layer (ABL) is presented. The procedure uses transilient turbulence theory to assess the turbulent fluxes from observed wind and temperature data, which are used to estimate the eddy transfer coefficients, so enabling a local closure ABL model to be integrated forward in time to a new steady state. The method has been applied to initialize kytoon data taken at Kharagpur (22.3° N, 87.2° E) during the MONTBLEX field programme. Results of a case study for 17 June 1990 are discussed.     

Application of transilient turbulent theory to study interactions between the atmospheric boundary layer and forest canopies

The new Forest-Land-Atmosphere ModEl called FLAME is presented. The first-order, nonlocal turbulence closure called transilient turbulence theory (Stull, 1993) is applied to study the interactions between a forested land-surface and the atmospheric boundary layer (ABL). The transilient scheme is used for unequal vertical grid spacing and includes the effects of drag, wake turbulence, and interference to vertical mixing by plant elements. Radiation transfer within the vegetation and the equations for the energy balance at the leaf surface have been taken from Norman (1979). Among others, the model predicts profiles of air temperature, humidity and wind velocity within the ABL, sensible and latent heat fluxes from the soil and the vegetation, the stomata and aerodynamic resistances, as well as profiles of temperature and water content in the soil. Preliminary studies carried out for a cloud free day and idealized initial conditions are presented. The canopy height is 30 m within a vertical domain of 3 km. The model is able to capture some of the effects usually observed within and above forested areas, including the relative wind speed maximum in the trunk space and the counter gradient-fluxes in the lower part of the plant stand. Of special interest is the determination of the location and magnitude of the turbulent mixing between model layers, which permits one to identify the effects of large eddies transporting momentum and scalar quantities into the canopy. A comparison between model simulations and field measurements will be presented in a future paper.     

A comparison of local and non-local turbulence closure methods for the case of a cold air outbreak

Numerical experiments have shown that large-eddy-simulation models (LES) are able to reproduce the common features of convective boundary layers (CBL) quite well. Models which cannot resolve the convective motions due to their grid structure (1D-models or models with coarse horizontal and/or vertical resolution) have to take into account the effects of large eddies within their subgrid diffusion terms. Turbulent fluxes are frequently parameterized through first-order-closure methods (K-theory). Recently, non-local closure schemes have also been developed. In this paper we compare 1D-and 2D-models using different local and non-local first-order closure methods. The analysis is carried out for the case of an idealized cold air outbreak (CAO). One of the non-local closures is based on the so-called transilient turbulence theory. The reference states are given by a bulk-model and a 2D-model which resolves the large eddies explicitly. A comparison of the results is presented for characteristic quantities such as evolution of boundary-layer height and surface heat flux as well as mean wind and temperature profiles. It is found that simple local first-order closure does not give good agreement with the reference models. The results of the transilient turbulence model shows that a non-local closure is able to parameterize the effects of the large eddies. Comparable results are produced by a local closure where eddy diffusivities are parameterized by dimensionless gradient-functions.     

The Adaptation of the Atmospheric Boundary Layer to a Change in Surface Roughness

The adaptation of the atmospheric boundary layer to a change in the underlying surface roughness is an interesting problem and hence much research, theoretical, experimental, and numerical, has been undertaken. Within the atmospheric boundary layer an accurate numerical model for the turbulent properties of the atmospheric boundary layer needs to be implemented if physically realistic results are to be obtained. Here, the adaptation of the atmospheric boundary layer to a change in surface roughness is investigated using a first-order turbulence closure model, a one-and-a-half-order turbulence closure model and a second-order turbulence closure model. Perturbations to the geostrophic wind and the pressure gradients are included and it is shown that the second-order turbulence closure model, namely the standard k - model, is inferior to a lower-order closure model if a modification to limit the turbulent eddy size within the atmospheric boundary layer is not included within the model.     

Modelling the One-Dimensional Stable Boundary Layer with an E−ℓ Turbulence Closure Scheme

The atmospheric boundary layer (ABL) model of Weng and Taylor with E−ℓ turbulence closure is applied to simulate the one-dimensional stably stratified ABL. The model has been run for nine hours from specified initial wind, potential temperature and turbulent kinetic energy profiles, and with a specified cooling rate applied at the surface. Different runs are conducted for different cooling rates, geostrophic winds and surface roughnesses. The results are discussed and compared with other models, large-eddy simulations and published field data.     

湍流能量闭合法用于大气模式的初探

湍流能量（TKE）闭合方案是近年来大气边界层（PBL）模式中技推崇的方法,它具有11－2阶精度。大气环流模式的计算容量和垂直网络限制,不允许PBL方案过份精细。本文对TKE方案进行简化及网络稀疏化试验。并用Wangara资料验证。验证结果表明：此方案引入PBL内仅设置5层的大气环流模式是可行的。地面以上400米内模式结果与实测吻合。但中、上部位温计算值偏高。关于本方案在大尺度模式中的应用尚须进一步研究。     

Prediction of wind properties in urban environments using artificial neural network

Artificial neural network (ANN) modeling has been performed to predict turbulent boundary layer characteristics for rough terrain based on experimental tests conducted in a boundary-layer wind tunnel to simulate atmospheric boundary layer using passive roughness devices such as spires, barriers, roughness elements on the floor, and slots in the extended test section. Different configurations of passive devices assisted to simulate urban terrains. A part of the wind tunnel test results are used as training sets for the ANN, and the other part of the test results are used to compare the prediction results of the ANN. Two ANN models have been developed in this study. The first one has been used to predict mean velocity, turbulence intensity, and model length scale factor. Results show that ANN is an efficient, accurate, and robust modeling procedure to predict turbulent characteristics of wind. In particular, it was found that the ANN-predicted wind mean velocities are within 4.7%, turbulence intensities are within 6.2%, and model length scale factors are within 3.8% of the actual measured values. In addition, another ANN model has been developed to predict instantaneous velocities that enables calculating the power spectral density of longitudinal velocity fluctuations. Results show that the predicted power spectra are in a good agreement with the power spectra obtained from measured instantaneous velocities.