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

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

Impact of Turbulent Mixing in the Stratocumulus-Topped Boundary Layer on Numerical Weather Prediction

The impact of enhanced turbulent mixing induced by radiative cooling at the top of the stratocumulus-topped boundary layer (STBL) on numerical weather prediction is examined. An additional term involving top-down turbulent mixing via in-cloud radiative cooling is applied to the Yonsei University (YSU) planetary boundary layer (PBL) parameterization scheme using a top-down diffusivity profile and cloud-top entrainment. The modified scheme is evaluated in an advection fog case over the Yellow Sea of Korea using the Weather Research and Forecasting (WRF) model and in global medium-range forecasts using the Global/Regional Integrated Model system (GRIMs). In the fog case simulation, consideration of the additional top-down mixing parameterization in the YSU PBL simulates less formation and more rapid dispersion of the fog. As a result, the modified scheme simulates a drier and warmer boundary layer and a moister and cooler layer above the PBL. The modified algorithm also improves surface temperature prediction over the Yellow Sea accompanying early dissipation of the fog. In the global medium-range forecast experiment, the modified scheme simulates overall enhanced PBL mixing over the STBL in the tropics and subtropical ocean, showing drier and warmer regions near the surface and moister and cooler regions above the PBL, resulting in prediction of reduced low level cloud amount and increased downward shortwave radiation at the surface. The modified scheme appears to improve systematic bias in temperature and humidity in the lower troposphere compared to the control simulation.     

Turbulent Helicity in the Atmospheric Boundary Layer

We consider the assumption postulated by Deusebio and Lindborg (J Fluid Mech 755:654–671, 2014) that the helicity injected into the Ekman boundary layer undergoes a cascade, with preservation of its sign (right- or alternatively left-handedness), which is a signature of the system rotation, from large to small scales, down to the Kolmogorov microscale of turbulence. At the same time, recent direct field measurements of turbulent helicity in the steppe region of southern Russia near Tsimlyansk Reservoir show the opposite sign of helicity from that expected. A possible explanation for this phenomenon may be the joint action of different scales of atmospheric flows within the boundary layer, including the sea-breeze circulation over the test site. In this regard, we consider a superposition of the classic Ekman spiral solution and Prandtl’s jet-like slope-wind profile to describe the planetary boundary-layer wind structure. The latter solution mimics a hydrostatic shallow breeze circulation over a non-uniformly heated surface. A 180°-wide sector on the hodograph plane exists, within which the relative orientation of the Ekman and Prandtl velocity profiles favours the left rotation with height of the resulting wind velocity vector in the lowermost part of the boundary layer. This explains the negative (left-handed) helicity cascade toward small-scale turbulent motions, which agrees with the direct field measurements of turbulent helicity in Tsimlyansk. A simple turbulent relaxation model is proposed that explains the measured positive values of the relatively minor contribution to turbulent helicity from the vertical components of velocity and vorticity.     

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

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

城市大气边界层的数值研究

本文利用能量闭合的二维非线性、非定常模式,结合地面热量平衡方程,研究了城市热岛和城市热岛环流以及其它相关的城市边界层参量的发展演变规律;同时探讨了在城市污染面源的情况下,污染物的扩散、输送问题.模拟结果得到了天津市城市热岛观测资料的证实.     

Influence of a Two-scale Surface Roughness on a Neutral Turbulent Boundary Layer

Flow in the urban boundary layer is strongly influenced by the surface roughness, which is composed principally of isolated buildings or groups of buildings. Previous research has shown that the flow regime depends on the characteristic height of these obstacles (H), and the spacing between them (W). In reality, the urban boundary layer contains roughness elements with a wide range of length scales; in many practical situations these can be classified into large-scale roughness—buildings, or groups of buildings—and small-scale roughness, such as street furniture and elements on the façades and roofs. It is important to understand how the small-scale roughness might modify mass and momentum transfer in the urban boundary layer, but relatively little information is available concerning the potential interaction between large- and small-scale roughness elements in the different flow regimes. This problem has been studied using wind-tunnel experiments, by measuring vertical velocity profiles over a two-dimensional obstacle array, adding small-scale roughness elements to the top of larger parallel square bars. The experiments were performed for different cavity aspect ratios: the results show that the small-scale roughness increases the turbulence intensities and the momentum transfer when the large-scale obstacles are closely packed (H/W > 1) but it has very little effect for more widely-spaced obstacles (H/W < 1).     

天津城市边界层湍流统计特征

利用2008年2月1～28日在天津气象塔上观测的超声风、温资料和常规风、温及湿梯度资料,计算了天津城市边界层无量纲湍流速度方差、湍流温度方差、感热通量、动量通量和湍流动能。结果表明,在不稳定层结条件下,40m高度上无量纲湍流速度方差和湍流温度方差遵循莫宁-奥布克霍夫相似理论,在220m高度上只有垂直方向上的湍流速度方差...     

Large-eddy Simulation of the Dispersion of Solid Particles in a Turbulent Boundary Layer

A large-eddy simulation (LES) with the dynamic Smagorinsky-Germano subgrid-scale (SGS) model is used to study the dispersion of solid particles in a turbulent boundary layer. Solid particles are tracked in a Lagrangian way. The instantaneous velocity of the surrounding fluid is considered to have a large-scale part (directly computed by the LES) and a small-scale part. The SGS velocity of the surrounding fluid is given by a three-dimensional Langevin model written in terms of SGS statistics at a mesh level. An appropriate Lagrangian correlation time scale is considered in order to include the influences of gravity and inertia of the solid particle. Inter-particle collisions and the influence of particles on the mean flow are also taken into account. The results of the LES are compared with the wind-tunnel experiments of Nalpanis et al. (1993 J Fluid Mech 251: 661–685) and of Tanière et al. (1997 Exp in Fluids 23:463–471) on sand particles in saltation and in modified saltation, respectively.     

A Simple Physically-Based Model for Wind-Turbine Wake Growth in a Turbulent Boundary Layer

A systematic evaluation is conducted of the smoothed spectrum, which is a spectral estimate obtained by averaging over a window of contiguous frequencies. The technique is extended to the ogive, as well as to the cross-spectrum. It is shown that, combined with existing variance estimates for the periodogram, the variance—and therefore the random error—associated with these estimates can be calculated in a straightforward way. The smoothed spectra and ogives are biased estimates; with simple power-law analytical models, correction procedures are devised, as well as a global constraint that enforces Parseval’s identity. Several new results are thus obtained: (1) The analytical variance estimates compare well with the sample variance calculated for the Bartlett spectrum and the variance of the inertial subrange of the cospectrum is shown to be relatively much larger than that of the spectrum. (2) Ogives and spectra estimates with reduced bias are calculated. (3) The bias of the smoothed spectrum and ogive is shown to be negligible at the higher frequencies. (4) The ogives and spectra thus calculated have better frequency resolution than the Bartlett spectrum, with (5) gradually increasing variance and relative error towards the low frequencies. (6) Power-law identification and extraction of the rate of dissipation of turbulence kinetic energy are possible directly from the ogive. (7) The smoothed cross-spectrum is a valid inner product and therefore an acceptable candidate for coherence and spectral correlation coefficient estimation by means of the Cauchy–Schwarz inequality. The quadrature, phase function, coherence function and spectral correlation function obtained from the smoothed spectral estimates compare well with the classical ones derived from the Bartlett spectrum.     

Influence of Evaporating Droplets in the Turbulent Marine Atmospheric Boundary Layer

Semi-idealized model simulations are made of the nocturnal cold-air pool development in the approximately 1-km wide and 100–200-m deep Grünloch basin, Austria. The simulations show qualitatively good agreement with vertical temperature and wind profiles and surface measurements collected during a meteorological field expedition. A two-layer stable atmosphere forms in the basin, with a very strong inversion in the lowest part, below the approximate height of the lowest gap in the surrounding orography. The upper part of the stable layer is less strongly stratified and extends to the approximate height of the second-lowest gap. The basin atmosphere cools most strongly during the first few hours of the night, after which temperatures decrease only slowly. An outflow of air forms through the lowest gap in the surrounding orography. The outflow connects with a weak inflow of air through a gap on the opposite sidewall, forming a vertically and horizontally confined jet over the basin. Basin cooling shows strong sensitivity to surface-layer characteristics, highlighting the large impact of variations in vegetation and soil cover on cold-air pool development, as well as the importance of surface-layer parametrization in numerical simulations of cold-air-pool development.     

Experimental Investigation of Turbulent Momentum Transfer in a Neutral Boundary Layer over a Rough Surface

The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.     

Inertial Effects on the Vertical Transport of Suspended Particles in a Turbulent Boundary Layer

In many atmospheric flows, a dispersed phase is actively suspended by turbulence, whose competition with gravitational settling ultimately dictates its vertical distribution. Examples of dispersed phases include snow, sea-spray droplets, dust, or sand, where individual elements of much larger density than the surrounding air are carried by turbulent motions after emission from the surface. In cases where the particle is assumed to deviate from local fluid motions only by its gravitational settling (i.e., they are inertialess), traditional flux balances predict a power-law dependence of particle concentration with height. It is unclear, however, how particle inertia influences this relationship, and this question is the focus of this work. Direct numerical simulations are conducted of turbulent open-channel flow, laden with Lagrangian particles of specified inertia; in this way the study focuses on the turbulent transport which occurs in the lowest few meters of the planetary boundary layer, in regions critical for connecting emission fluxes to the fluxes felt by the full-scale boundary layer. Simulations over a wide range of particle Stokes number, while holding the dimensionless settling velocity constant, are performed to understand the role of particle inertia on vertical dispersion. It is found that particles deviate from their inertialess behaviour in ways that are not easily captured by traditional theory; concentrations are reduced with increasing Stokes number. Furthermore, a similarity-based eddy diffusivity for particle concentration fails as particles experience inertial acceleration, precluding a closed-form solution for particle concentration as in the case of inertialess particles. The primary consequence of this result is that typical flux parametrizations connecting surface emission models (e.g., saltation models or sea-spray generation functions) to elevated boundary conditions may overestimate particle concentrations due to the reduced vertical transport caused by inertia in between; likewise particle emission may be underestimated if inferred from concentration measurements aloft.     

Turbulent Pressure Statistics in the Atmospheric Boundary Layer from Large-Eddy Simulation

We use large-eddy simulation (LES) to study the turbulent pressure field in atmospheric boundary layers with free convection, forced convection, and stable stratification. We use the Poisson equation for pressure to represent the pressure field as the sum of mean-shear, turbulence–turbulence, subfilter-scale, Coriolis, and buoyancy contributions. We isolate these contributions and study them separately. We find that in the energy-containing range in the free-convection case the turbulence–turbulence pressure dominates over the entire boundary layer. That part dominates also up to midlayer in the forced-convection case; above that the mean-shear pressure dominates. In the stable case the mean-shear pressure dominates over the entire boundary layer.We find evidence of an inertial subrange in the pressure spectrum in the free and forced-convection cases; it is dominated by the turbulence–turbulence pressure and has a three-dimensional spectral constant of about 4.0. This agrees well with quasi-Gaussian predictions but is a factor of 2 less than recent results from direct numerical simulations at moderate Reynolds numbers. Measurements of the inertial subrange pressure spectral constant at high Reynolds numbers, which might now be possible, would be most useful.     

Disentangling Turbulent Gas Diffusion from Non-diffusive Transport in the Boundary Layer

Boundary-Layer Meteorology - An analysis based on the law of linear momentum conservation demonstrates unequivocally that the mass fraction is the scalar whose gradient determines gas diffusion,...     

沙漠-绿洲大气边界层结构的数值模拟

利用美国NCAR新版中尺度MM5V3．6非静力平衡模式,采用三重嵌套的模拟方法,模拟研究了沙漠绿洲的环流及边界层特征。结果表明：在没有西风气流背景的影响下,绿洲沙漠环流和大气边界层结构是对称的。沙漠绿洲改变了原有沙漠地区环流结构及温、湿场的分布,绿洲上空大气下沉,沙漠上空大气上升,从而产生了绿洲上空大气冷干,沙漠上空大气暖湿的边界层特征。绿洲边缘形成了由干到湿的强湿度梯度带围绕着绿洲,起到了保护绿洲的作用。在有西风背景气流的影响下,绿洲沙漠环流和大气边界层结构是非对称的。但是,西风背景气流的存在可以破坏绿洲系统对称的环流结构,不利于绿洲系统的稳定发展。     

Numerical Study of a Mesoscale Vortex in the Planetary Boundary Layer of the Meiyu Front

It was found that the heavy rainfall event along the Meiyu front in the lower reaches of the Yangtze River on 23 June 2009 was connected with a mesoscale disturbance vortex, which originated from the planetary boundary layer (PBL) and developed upward later and was discovered by using the Shuman-Shapiro filtering method. The mesoscale disturbance vortex in the PBL (PMDV) in this process corresponded well to the short-time rainstorm in the Doppler radar echo. Analysis of the high-resolution simulation results from the Advanced Weather Research and Forecasting Model (ARW) showed that there were several surface disturbances along the southern warm section of the Meiyu front prior to the generation of the PMDV. The PMDV interacted with the mesoscale convective system (MCS) and intensiˉed the local convective precipitation. The north and southwest flows in the PBL converged at the time of the PMDV formation. Meanwhile, a southwesterly jet on the top of the PBL to the south side of the vortex reinforced the ascending motion and convergence. Hence, it is concluded that the PMDV was generated when the strong cold air flows north of the shear line encountered the southwest flow south of the shear line. The convergence line in the PBL, the intensification of the southwest wind, and the southward aggression of the north wind were critical for the development of the PMDV. The release of latent heat was found crucial for the formation of the PMDV as it facilitated the convergence at low levels.     

A Three-Dimensional Correlation/Spectral Model for Turbulent Velocities in a Convective Boundary Layer

A three-dimensional model for correlation functions and spectra in theatmospheric, convective boundary layer (CBL) is presented. The modelincludes vertical inhomogeneities introduced by eddy-blocking at the ground.By assuming the disturbance to the turbulent flow resulting from the groundblocking is irrotational, an equation is developed which allows one to writethe inhomogeneous, two-dimensional (2D) cross spectra for the blocked flowin terms of the 2D cross spectra for a homogeneous flow. VonKármán's energy spectrum then is used to determine thehomogeneous, 2D cross spectra. Although there are only two adjustableparameters in the model, the variance and a length scale, the model is shownto agree quite well with a diversity of previous results for the CBL.     

夜间边界层的一维数值模拟

本文用一维数值模式,模拟了平坦地形情况下小风、晴空夜间边界层内风、温、湿的演变,计算结果与观测事实有较好的一致性。对有高云或少量中、低云存在时的计算结果分析表明,模式对风速的模拟仍然较好,对温度的模拟比实测值可偏低0.3—1℃,对相对湿度的模拟比实测值可偏高5—10％在符合模式使用范围的情况下,本模式可作为夜间边界层预报的一种方法参考使用。     

不同下垫面的大气边界层湍流扩散特征

利用海洋、戈壁、城郊和城市4种不同下垫面上的湍流观测资料,运用多尺度湍流概念和方法计算分析了风速谱和扩散参数等湍流特征量。结果表明:下垫面的粗糙度对大气边界层湍流特征影响很大,随着粗糙度的增大,速度谱的峰值频率向高频方向移动,扩散参数增大。