Acoustic sounding of diurnal variations and gravity waves in the planetary boundary layer

In the introductory part the paper presents some SODAR results on diurnal variations of different boundary layer structures from measurements during the KOPEX-experiment at a city station and a suburban station. In Sect. 1 the investigations about wave processes in the planetary boundary layer are discussed. The SODAR measurements were carried out in different international cooperations at the Kopisty Atmospheric Observatory of the Institute of Physics of the Atmosphere Prague (northern Czechoslovakia) in the year 1986 (KOPEX-86), at the Observatory Juliusruh of the Heinrich-Hertz-Institut Berlin at the island of Rügen (JULEX-87) and at the Meteorological Observatory of the nuclear power plant Jaslovske Bohunice in Slovakia, about 60 km north-east of Bratislava (JABEX-88).Formerly: Heinrich-Hertz-Institute, Berlin-Adlershof, Germany     

Aspect sensitivity of sound backscattering in the atmospheric boundary layer

The dependence of backscatter intensity on the inclination of a SODAR beam was measured in the nocturnal boundary layer. The aspect sensitivity was revealed, but the angular dependence was not so strong as that of radiowave backscattering in the upper atmosphere. The mechanism of the phenomenon is under discussion as well as its influence on the accuracy of SODAR measurements of C _T ² values which are very useful for air-pollution meteorology.     

Variances of the vertical and horizontal wind measured by tower instruments and SODAR

At the Karlsruhe Nuclear Research Center, a monostatic Doppler SODAR has been operated near a 200 m high meteorological tower. Vertical profiles of the components of the wind vector and of the variances of the vertical wind and of the horizontal wind direction have been sampled continuously from the SODAR and the tower instruments as 30-min mean values. During seven episodes lasting more than 2 h, the instantaneous vertical wind speed was measured simultaneously by a sonic anemometer and the SODAR, and spectra were calculated. An intercomparison of the spectra and wind data measured directly by the SODAR and the tower instruments has been performed by a linear regression and correlation analysis referring to different height levels of measurement and stability classes.Systematic and statistically distributed differences among the data measured by the SODAR and the in situ instruments are discussed. Our investigation shows that the variance of only the vertical wind speed can be measured reliably by a SODAR. The variances of the horizontal wind direction and the elevation angle of the wind vector should not be used in a routine manner from a SODAR.     

Similarity methods to derive turbulence quantities and mixed-layer depth from SODAR measurements in the convective boundary layer: A review

A brief review is made of similarity methods which have been developed to calculate turbulence parameters and mixed-layer depth from SODAR measurements. Emphasis is placed on surface fluxes and mixed-layer depth. A variety of parametric relations valid in the convective boundary layer are surveyed and relevant similarity methods are outlined. The methods are primarily intended for flat and relatively homogeneous terrain. It is shown that, in most cases, horizontal homogeneity provides a less stringent constraint.Finally, some limitations of the existing methods are presented and discussed. It is proposed that the most serious errors are mainly associated with the poor accuracy of some parameters measured by the SODAR and also with the limitations of the theoretical basis employed.     

SODAR study of the boundary layer over moscow for air-pollution application

The 24-hour SODAR measurements of the atmospheric boundary-layer parameters over Moscow were carried out during one month in the summer of 1990 and one month in the winter of 1991 at three points within the city and suburbs. The simultaneous measurements of the concentrations of several pollutants were also carried out at one of these points. The relations between local concentrations of some pollutants and the parameters measured by SODAR (mixing height and wind speed) are revealed. This fact confirms the usefulness of SODAR observations for air-pollution meteorology. The statistics of the boundary-layer parameters in summer and winter seasons as well as the data on their spatial variations over the city were obtained.     

确定大气边界层顶高度的新方法及数值实验

提出了一种确定大气边界层顶高度的数值微分新方法,该方法使用了正则化技术,把对弯角廓线求导数的数值微分问题转化为求目标泛函极小值的问题,采用双参数模型函数方法来选择正则化参数,最后利用最大梯度法确定边界层顶高度.首先通过两个数值实验验证了新方法的有效性,实验结果显示,随着掩星资料噪音的增多,由差分法和结合L曲线方案的数值微分方法得到的边界层顶高度波动增大,而通过双参数模型函数方法得到的高度很稳定,这说明新方法能够很好地过滤噪音,从而保留廓线的主要信息.随后基于2007-2011年1,4,7,10月的COSMIC弯角数据,利用新方法分析了全球海洋大气边界层顶高度的季节特征,并与用掩星资料自带的大气边界层顶高度数据zbalmax得到的季节分布进行对比.结果表明,两者的季节分布特征十分一致:海温相对周围海域高的区域,边界层顶高度较高,反之,边界层顶高度较低;在暖流经过的海域,边界层顶高度较高,在寒流经过的海域,边界层顶的高度相对较低.     

超声速平板圆台突起物绕流实验和数值模拟研究

高速飞行器表面不可避免的存在突起物并形成复杂流场, 从而引起飞行器气动特性和热载荷的变化; 同时, 突起物是流动控制的重要方法之一, 合适的突起物形状及安装位置对于改善冲压发动机进气道性能有重要意义. 本文采用基于纳米粒子的平面激光散射技术(NPLS)研究了马赫3.0来流边界层为层流的平板上三个不同高度圆台突起物绕流流场, 主要关注了突起物后方的尾迹边界层, 并采用高精度的显式五阶精度加权紧致非线性格式(WCNS-E-5)离散求解Navier-Stokes方程模拟了该流场. 获得了超声速圆台绕流精细流场结构, 观察到突起物后方尾迹区域边界层发展的过程. 结合实验和数值模拟结果可以发现, 当圆台高度接近或者小于当地边界层厚度时, 突起物对边界层的扰动非常弱, 圆台后方尾迹边界层能够维持较长距离的层流状态, 在边界层转捩阶段也有清晰的发卡涡结构出现; 反之, 边界层受到的扰动明显增大, 在突起物后方很快发展为湍流; 风洞噪声对本文研究圆台引起的边界层扰动有一定影响, 实验获得的边界层转捩位置要比数值结果靠前. 基于NPLS流场图像, 采用间歇性方法分析了圆台突起物后方边界层的特性, 对于高度大于边界层厚度的圆台其间歇性曲线较为接近并且更加饱满, 边界层的脉动也更为强烈.     

On the use of spires for generating inflow conditions with energetic coherent structures in large eddy simulation

While it has long been a practice to place spires near the inlet of a wind tunnel to quickly develop a turbulent boundary layer with similarities to an atmospheric boundary layer, this has not been the case for creating turbulent boundary layer inflow in large eddy simulations (LESs) of turbulent flows. We carry out LES with the curvilinear immersed boundary method to simulate the flow in a wind tunnel with a series of spires in order to investigate the feasibility of numerically developing inflow conditions from a precursory spire LES and assessing the similarities of the turbulence statistics to those of an atmospheric boundary layer. The simulated mean velocity field demonstrates that a turbulent boundary layer with height equal to the spire height develops very quickly, within five spire heights downstream. The major attribute of using spires for precursory simulations is the spatially evolving coherent structures that form downstream of the spires offering a range of length scales at both the vertical and streamwise directions allowing multiple turbulent inflow conditions to be extracted from a single simulation. While the distribution of length scales far from the spires resembles an atmospheric boundary layer, some turbulence statistics have some significant differences.     

Formation of the "superconducting" core in turbulent thermal convection

A sinusoidal temperature perturbation is superimposed on the bottom plate of a cylindrical convection cell, and its decay is measured at the cell midheight. Rayleigh numbers up to Ra=10(13) and aspect ratios 1 and 4 are considered. The technique allows a dynamic measurement of the height of the layer interposed between the superconducting core and the boundary. This deduced height is in good agreement with results from recent numerical simulations.     

Heat and mass transfer characteristics during rapid solidification of Fe-Cu peritectic alloys

The viscose flow and microstructure formation of Fe-Cu peritectic alloy melts are investigated by analyzing the velocity and temperature fields during rapid solidification, which is verified by rapid quenching experiments. It is found that a large temperature gradient exists along the vertical direction of melt puddle, whereas there is no obvious temperature variation in the tangent direction of roller surface. After being sprayed from a nozzle, the alloy melt changes the magnitude and direction of its flow and velocity rapidly at a height of about 180 μm. The horizontal flow velocity increases rapidly, but the vertical flow velocity decreases sharply. A thermal boundary layer with 160–300 μm in height and a momentum boundary layer with 160–240 μm in thickness are formed at the bottom of melt puddle, and the Reynolds number Re is in the range of 870 to 1070 in the boundary layer. With the increase of Re number, the cooling rate increases linearly and the thickness of thermal boundary layer increases monotonically. The thickness of momentum boundary layer decreases slowly at first, then rises slightly and decreases sharply. If Re < 1024, the liquid flow has remarkable effects on the microstructure formation due to dominant momentum transfer. The separated liquid phase is likely to form a fiber-like microstructure. If Re>1024, the heat transfer becomes dominating and the liquid phase flow is suppressed, which results in the formation of fine and uniform equiaxed microstructures. Supported by the National Natural Science Foundation of China (Grant Nos. 50121101 and 50395105)     

MAX-DOAS仪器观测大气Ring效应及对气溶胶参数的敏感性研究

Ring效应是指大气中O₂和N₂分子对太阳光的转动拉曼散射致使太阳夫琅禾费结构变浅(被填充)的现象。大气气溶胶能够改变光子在大气中的光程和大气散射性质,进而影响到光子发生转动拉曼散射的几率(RSP),最终影响填充效应。通过观测RSP在不同气溶胶状态下的变化,可以反演得到气溶胶参量信息。采用地基多轴差分吸收光谱(multi-axis differential optical absorption spectroscopy, MAX-DOAS)方法在晴朗无云天气下对Ring效应进行了观测,并把测量值和模型值进行了对比,两者一致性较好;选取大气辐射传输模型McArtim研究了在不同大气条件下Ring效应对气溶胶参数等的灵敏度,结果表明在大多数测量情况下,气溶胶光学厚度和边界层高度对RSP影响较大,在90°仰角时,AOD从0.1增加到1,RSP减少了24.6%,边界层高度从1 km增加到3 km,RSP增加了4.4%。研究表明,Ring效应对气溶胶光学厚度和边界层高度较为敏感,这为反演气溶胶的垂直分布提供了一种新方法。     

Generation of Turbulent Inflow Data for Spatially-Developing Boundary Layer Simulations

A method for generating three-dimensional, time-dependent turbulent inflow data for simulations of complex spatially developing boundary layers is described. The approach is to extract instantaneous planes of velocity data from an auxiliary simulation of a zero pressure gradient boundary layer. The auxiliary simulation is also spatially developing, but generates its own inflow conditions through a sequence of operations where the velocity field at a downstream station is rescaled and re-introduced at the inlet. This procedure is essentially a variant of the Spalart method, optimized so that an existing inflow–outflow code can be converted to an inflow-generation device through the addition of one simple subroutine. The proposed method is shown to produce a realistic turbulent boundary layer which yields statistics that are in good agreement with both experimental data and results from direct simulations. The method is used to provide inflow conditions for a large eddy simulation (LES) of a spatially evolving boundary layer spanning a momentum thickness Reynolds number interval of 1530–2150. The results from the LES calculation are compared with those from other simulations that make use of more approximate inflow conditions. When compared with the approximate inflow generation techniques, the proposed method is shown to be highly accurate, with little or no adjustment of the solution near the inlet boundary. In contrast, the other methods surveyed produce a transient near the inlet that persists several boundary layer thicknesses downstream. Lack of a transient when using the proposed method is significant since the adverse effects of inflow errors are typically minimized through a costly upstream elongation of the mesh. Extension of the method for non-zero pressure gradients is also discussed.     

机体表面湍流边界层噪声特性及预测方法研究

湍流边界层噪声是飞机巡航过程中的主要外部噪声源,对舱内噪声水平的影响尤为重要。因此,对飞机机体表面湍流边界层噪声的研究具有重要意义。本文通过试验获得了某型民机巡航过程中的湍流边界层噪声,试飞工况为3500ft/0.78、3500ft/0.7、2500ft/0.67、1500ft/0.66。对实测数据进行分析,发现湍流边界层噪声与动压、边界层厚度等参数有关。同时,利用计算流体力学的方法得到了飞机机体表面的压力分布,并分析了压力梯度对湍流边界层噪声的影响。最后,基于工程预测方法对湍流边界层噪声进行了预测,对于不存在逆压梯度的区域,预测结果与试验结果吻合较好,仅部分频段存在一定偏差。通过对模型的参数进行优化,改善了预测结果。     

Mie散射激光雷达研究大气边界层特性

基于Mie散射大气激光雷达,应用Klett算法对回波信号进行了反演,得到大气消光系数和后向散射系数。然后,通过对大气后向散射系数曲线进行拟合得到大气边界层混合层高度以及卷夹层厚度等特征参数。对比雷达实测数据与探空气球所测数据,证明二者存在较好相关性。基于实测数据讨论了四川盆地大气边界层特性对气候的影响。     

Vertical profile of aerosol extinction based on the measurement of O_4 of multi-elevation angles with MAX-DOAS

A method for aerosol extinction profile retrieval using ground-based multi-axis differential optical absorption spectroscopy(MAX-DOAS) is studied, which is based on a look-up table algorithm. The algorithm uses parametric method to represent aerosol extinction profiles and simulate different atmospheric aerosol states through atmospheric radiation transfer model. Based on the method, aerosol extinction profile was obtained during six cloud-free days. The O_4 differential air mass factor(dAMF) measured by MAX-DOAS is compared with the corresponding model results under different atmospheric conditions(R~2= 0.78). The aerosol optical thickness, aerosol weight factor in boundary layer, and the height of the boundary layer are obtained after the process of minimization and look-up table method. The retrieved aerosol extinction in boundary layer is compared with PM2.5 data measured by ground point instrument. The diurnal variation trends of the two methods are in good agreement. The correlation coefficients of the two methods are 0.71 when the aerosol optical thickness is smaller than 0.5. The results show that the look-up table method can obtain the aerosol state of the troposphere and provide validation for other instrument data.     

激光雷达探测北京城区夏季大气边界层

为了研究北京城区夏季大气边界层结构变化特征及大气边界层内气溶胶消光特性,2004年8月利用便携式米散射激光雷达对北京城区夏季大气边界层进行了系统观测。反演了观测站上空大气气溶胶的消光特性垂直分布以及大气边界层的高度。分析了气象条件和人类活动对大气边界层结构的影响。观测数据表明： 北京城区夏季大气边界层有明显的日变化特征,早晚比较低,日间有一个从低变高再变低的过程,中午前后达到最高。结合气象参数对测量数据进行的统计分析表明：北京城区夏季大气边界层高度相对稳定,多分布在1.8km以下,平均值为0.68km;大气边界层内存在浓度较高的气溶胶粒子,平均光学厚度（3km以内）在0.30左右。     

Sound from boundary layer flow over steps and gaps

This study is concerned with the radiated sound from boundary layer flows over small forward and backward steps and gap configurations of similar dimension. These measurements were performed in the Virginia Tech Anechoic Wall Jet Facility for step heights that ranged from approximately 10 percent to 100 percent of the incoming boundary layer height. The results show the influence of step height and boundary layer edge velocity on the far-field sound from forward and backward steps. Neither source shows clear dipole directivity and at least the larger step heights considered in this study are shown to not be acoustically compact. A new mixed scaling normalization is proposed for the far-field spectra from both types of step. Backward steps are shown to be much weaker producers of far-field sound than similarly sized forward steps. The implications of this behavior are discussed with respect to the far-field sound measured from various gap flows.     

Turbulent plane wall jets over smooth and rough surfaces

Large-eddy simulations were carried out to study the effects of surface roughness on a plane wall-jet using the Lagrangian dynamic eddy-viscosity subgrid-scale model, at Re = 7500 (based on the jet bulk velocity and height). Results over both smooth and rough surfaces were validated by experimental data at the same Reynolds number. As the jet is injected into the still environment, large-scale rollers are generated in the shear layer between the high-momentum fluid of the jet and the surrounding and are convected downstream with the flow. To understand the extent to which the outer-layer structures modify the flow in the inner layer and the extent to which the effect of roughness spreads away from the wall, both instantaneous and mean flow fields were investigated. The results revealed that, for the Reynolds number and roughness height considered in this study, the effect of roughness is mostly confined to the near-wall region of the wall jet. There is no structural difference between the outer layer of the wall jet over the smooth and rough surfaces. Roughness does not affect the size of the outer-layer structures or the scaling of the profiles of Reynolds stresses in the outer layer. However, in the inner layer, roughness redistributes stresses from streamwise to wall-normal and spanwise directions toward isotropy. Contours of joint probability-density function of the streamwise and wall-normal velocity fluctuations at the bottom of the logarithmic region match those of the turbulent boundary layer at the same height; while the traces of the outer-layer structure were detected at the top of the logarithmic region, indicating that they do not affect the flow very close to the wall, but still modify a major portion of the inner layer. This modification must be taken into consideration when the inner layer of a wall jet is compared with the conventional turbulent boundary layer.     

Revisiting the boundary layer structure used in Craig and Gordon's model of isotope fractionation in evaporation

Craig and Gordon's (CG) model of isotope fractionation in evaporation was derived more than 40 years ago and was based on the turbulent boundary layer structure model acceptable at that time. That view suggested that turbulent flows consist of eddies with a wide range of length scales moving randomly in the flow domain. There is evidence that some parameters in CG model do not fully correspond to data in the literature. Owing to advances in fluid dynamics research techniques, it has been shown in recent decades that the apparent chaotic flow in turbulent boundary layers is in fact governed by few well-organised structures. This article reviews the major characteristics of these coherent structures based on available results from experimental, numerical and theoretical studies of turbulent and laminar boundary layers. The review points on some differences between past and present views of the boundary layer structure and on their relation and possible influence on power laws in CG model.