Parameterization of meandering phenomenon in a stable atmospheric boundary layer

Accounting for the current knowledge of the stable atmospheric boundary layer (ABL) turbulence structure and characteristics, a new formulation for the meandering parameters to be used in a Lagrangian stochastic particle turbulent diffusion model has been derived. That is, expressions for the parameters controlling the meandering oscillation frequency in low wind speed stable conditions are proposed. The classical expression for the meandering autocorrelation function, the turbulent statistical diffusion theory and ABL similarity theory are employed to estimate these parameters. In addition, this new parameterization was introduced into a particular Lagrangian stochastic particle model, which is called Iterative Langevin solution for low wind, validated with the data of Idaho National Laboratory experiments, and compared with others diffusion models. The results of this new approach are shown to agree with the measurements of Idaho experiments and also with those of the other atmospheric diffusion models. The major advance shown in this study is the formulation of the meandering parameters expressed in terms of the characteristic scales (velocity and length scales) describing the physical structure of a turbulent stable boundary layer. These similarity formulas can be used to simulate meandering enhanced diffusion of passive scalars in a low wind speed stable ABL.     

半干旱草原下垫面动量和感热总体输送系数参数化关系研究

本文利用锡林郭勒草原2008年春季近地层涡旋相关系统和铁塔的风、 温平均梯度观测资料, 分析了总体输送系数随梯度Richardson数的变化特征, 建立了动量总体输送系数随大气稳定度、近地层风速以及感热总体输送系数随大气稳定度和近地层气温的关系. 中性条件下, 半干旱草原植被下垫面动量总体输送系数与近地层大气动力状态之间存在明显的相互作用, 总体输送系数与近地面层风速之间满足二次曲线拟合关系; 风速较小时, 大气动力特征对地表粗糙度长度的改变不是很明显, 动量总体输送系数随气流增强而增大; 而当风速较大时, 强风速会使植被高度发生改变, 动量总体输送系数随气流增强而减小. 另外, 感热总体输送系数与近地层气温之间也存在二次曲线关系. 动量总体输送系数与近地层风速之间的关系、感热总体输送系数与近地层气温之间关系的建立为总体输送系数参数化提供了重要途径, 同时该方案避免了对动力学粗糙度长度和热力学粗糙度长度的求解. 关键词： 总体输送系数 参数化 湍流通量 相似性函数     

一个新的海边光学湍流外尺度和C_n~2的廓线模式

2016年12月13日至2017年1月2日期间,在茂名博贺海洋气象科学实验基地,采用自行研制的湍流气象探空仪,获取了30份海边温、湿、压、风速、风向和C_n~2等探空数据.基于HMNSP99外尺度模式,利用海边的探空数据拟合得到一个茂名大气光学湍流外尺度经验公式.同时对实验测得的高空湍流廓线数据进行统计平均,然后基于Hufnagel-Valley模式拟合得到符合海边湍流廓线规律的统计平均模式(C_n~2sea model).根据Tatarski高空湍流参数化方案,将用茂名外尺度公式估算的C_n~2分别与探空测量的C_n~2以及用其他外尺度模式估算的C_n~2进行了比较.对其进行统计性分析发现,利用新拟合的茂名外尺度公式、HMNSP99,Dewan以及(Coulman等外尺度模式计算的log10C_n~2)与实测值的整体相关系数分别为0.924,0.848,0.763和0.651,在变化趋势和量级上都表现出较好的一致性;以上四种外尺度模式估算结果的误差都很小,其整体平均绝对误差和平均相对误差分别为0.514和2.963%,0.627和3.612%,0.943和5.439%,0.766和4.417%,新拟合的外尺度模式的误差最小.进一步验证了新的海边外尺度和C_n~2廓线模式的可靠性和有效性,此外还发现高空大气光学湍流的发生与风切变和温度梯度具有十分密切的关系,为光电工程在海边场景应用所需的大气光学湍流廓线模式提供支持.     

Dynamics of momentumless turbulent wake in a shear flow of a linearly stratified medium

The numerical model of momentumless turbulent wake in a horizontally homogeneous shear flow of linearly stratified medium has been constructed. Based on this model, the investigation of the wake dynamics has been performed. The obtained data demonstrate the transformation of the zone of turbulent perturbations and internal waves generated by the wake under the action of shear flow, which leads to the deceleration of turbulence decay at large time values after the body passage.     

Characterisation of mixing rate due to high power ultrasound

For meaningful assessment of results from laboratory and pilot plant trials, it is often necessary to know the mixing characteristics within the ultrasonic reactors. Previous workers have used conductivity measurements in an attempt to characterize the residence time distribution in ultrasonic reactors, but these results do not provide direct data on the mixing within the high power region adjacent to the ultrasonic probe tip. We overcome this difficulty through direct visualization of the mixing process within the high energy region close to the tip of the ultrasonic probe. Our analysis proceeds by determining an approximate turbulent diffusivity in a batch reactor arrangement for different values of ultrasonic energy input. For input electrical power levels between 70 and 120 W and a processing volume of 30 ml, the effective turbulent diffusivity varied from about 0.2 x 10(-3) to 0.7 x 10(-3)m(2)/s. We demonstrate that such results can be coupled to a suitable dispersion model to estimate the actual residence time distribution in flow-though arrangements when the through-put adds little additional mixing energy. Therefore, coupling the effective turbulent diffusivity identified in a batch reactor with a suitable dispersion model for the reactor offers an alternative approach to the deduction of RTD when determining the actual RTD in the high intensity zone of steady flow sonochemical reactors is problematic.     

无人机视角下的林木爆燃火特性研究

森林火灾“爆燃”现象的特征是突然发生的高强度、高蔓延速度的燃烧。目前为止，关于“爆燃火”的原因还没有达成共识。以无人机视角下对林木爆燃火特性研究，以四川木里特大森林火灾为研究对象，通过分析凉山州某森林扑火部队3月31日木里森林火灾当天KWT（科卫泰）无人机航拍火场画面，结合无人机实时影像及实地调研数据，分析了峡谷地形林火蔓延时空特征，探讨了峡谷中风向风速变化时空分布规律，研究了地形变化条件下，不同海拔高度风速特征，建立了无人机倾角测量风速模型（其中为风速m·s-1，为无人机倾角°）。结果表明，高山峻岭特殊地形环境下每天4:00—12:00时间段为静风期，为峡谷林火扑救最佳时期；午后15:00—17:00和晚上20:00—22:00为山谷地形风速活跃期；仿真软件数据显示山顶、谷底与山腰不同海拔位置的风速风向不统一，谷底会产生乱流现象，且风速与海拔不存在正相关关系，小气候在复杂地形中占主导影响地位；在谷口至山谷深处的中间位置会出现气流速度的波峰状态，并易形成乱流，为爆燃火发生提供了客观必要条件。该研究可为复杂地形环境下，森林草原火灾扑救安全提供数据和技术支撑。     

Estimate for the doppler shift of a non-Gaussian signal upon coherent detection of scattered optical radiation in a turbulent atmosphere

A non-Gaussian model for estimating the radial velocity of turbulent flows in the atmosphere for coherent detection of scattered optical radiation is proposed. The model was obtained based on a theoretical approach that includes results of the statistical analysis of a pulse Doppler lidar signal in a turbulent medium, as well as on the perturbation-theory methods that have been developed in the theory of probability and mathematical statistics. It is shown that the estimate of the Doppler shift in the first-order perturbation theory is a sum of a regular component and two conditional fluctuation components—Gaussian and non-Gaussian ones. In the case of a homogeneous and isotropic turbulence, the estimate of the radial wind velocity is approximately equal to its true average value. The statistical uncertainty in measurements of the average radial wind velocity is determined by the behavior of conditional Gaussian and non-Gaussian components and significantly depends on the state of atmospheric turbulence. It is shown that basic equations of the non-Gaussian model in the limit case coincide with formulas of the local and nonlocal models, as well as with those of the Gaussian model.     

Jump at the onset of saltation

We reveal that the transition in the saturated flux for aeolian saltation is generically discontinuous by explicitly simulating particle motion in turbulent flow. This is the first time that a jump in the saturated flux has been observed. The discontinuity is followed by a coexistence interval with two metastable solutions. The modification of the wind profile due to momentum exchange exhibits a maximum at high shear strength.     

First and second order semi-Markov chains for wind speed modeling

The increasing interest in renewable energy, particularly in wind, has given rise to the necessity of accurate models for the generation of good synthetic wind speed data. Markov chains are often used for this purpose but better models are needed to reproduce the statistical properties of wind speed data. We downloaded a database, freely available from the web, in which are included wind speed data taken from L.S.I. -Lastem station (Italy) and sampled every 10 min. With the aim of reproducing the statistical properties of this data we propose the use of three semi-Markov models. We generate synthetic time series for wind speed by means of Monte Carlo simulations. The time lagged autocorrelation is then used to compare statistical properties of the proposed models with those of real data and also with a synthetic time series generated through a simple Markov chain.     

C波段机载合成孔径雷达海面风场反演新方法

针对基于散射计地球物理模型函数的机载合成孔径雷达(SAR)海面风场反演中存在的风向获取依赖于图像风条纹或数值预报、 散射计数据和浮标等背景场资料, 风向与SAR图像时空分辨率不匹配, 进而影响机载SAR海面风场反演精度等问题, 本文根据机载SAR对海探测特点, 研究一种适用于C波段机载SAR的海面风场反演新方法. 利用SAR图像距离向不同入射角的后向散射系数, 依据地球物理模型构造最小代价函数, 通过代价函数的求解直接从机载SAR数据同时反演出海面风速和风向. 利用论文提出的海面风场反演方法分别对仿真SAR数据和实测C波段机载SAR数据进行风向、 风速的反演误差分析及试验验证.研究结果表明, 该方法适用于机载SAR海面风场反演, 可不依赖背景风向直接反演出精度较高的风速和风向; 雷达后向散射系数误差是决定海面风速、风向反演精度的关键因素, 辐射定标精度越高则反演误差越小; 海面风速反演误差随着风速的提高而增大, 当海面风速大于18 m/s时, 风速反演误差显著增加, 而海面风向的反演误差与风速无明显关系. 关键词： 机载合成孔径雷达 海面风场 多入射角     

Quantitative wind tunnel studies using pressure- and temperature sensitive paints

The pressure sensitive paint (PSP) intensity and lifetime system is an optical measurement technique to investigate absolute pressure fields on model surfaces for basic research in laboratories, industrial wind tunnels or high speed rotating turbo machines. Detailed qualitative and quantitative information and understanding of flow phenomena can be obtained in speed ranges from U_∞=20 m/s up to Ma=5.0. A number of projects of industrial interest has been investigated in different wind tunnels covering low speed, transonic, trisonic and cryogenic facilities. The influence of the main error sources for the components of the PSP system have been checked. Comparison of experimental pressure fields obtained by means of PSP and the results of numerical calculations have been carried out. Different wind tunnel models ranging from basic configurations such as a cropped delta wing to a complex half model of a large propeller-driven transport aircraft with all flaps, rudders and shrouds, and rotating or oscillating models as well as Reynolds number effects on models have been investigated.     

Granular ratchets

We study the statistics of the horizontal component of atmospheric boundary layer wind speed and introduce a stochastic process which has similar properties. Motivated by the non-stationarity of wind velocity data, we describe statistical methods to verify the picture of natural atmospheric boundary layer turbulence to be composed of successively occurring close to ideal turbulence phases with different parameters. We focus on the fluctuation of wind speed around its mean behaviour and show that there is linear relationship between the standard deviation of the fluctuation and the current mean wind speed. Additionally, we analyse the increment statistics and investigate the time dependence of the parameters describing the increment distribution. Noting that the first order geometric auto regressive process has similar statistical properties, we investigate this similarity and analyse the extent to which this stochastic process is a suitable model for wind speed simulation.     

Density and velocity statistics in variable density turbulent mixing

We experimentally study variable–density mixing of miscible gases in an open-circuit wind tunnel using simultaneous particle image velocimetry and planar laser-induced fluorescence. Experiments of a high Atwood number (0.6) and low Atwood number (0.1) are performed to compare non-Boussinesq cases with the Boussinesq limit. The higher density gas is injected into the wind tunnel co-flow using a round jet configuration, and near-field and far-field measurements are performed to examine mixing in both momentum and buoyancy-dominated regimes. The effects of buoyancy are measurable and important in both large-scale mixing features and in turbulence quantities. The low Atwood number PDFs (probability density functions) show fast and uniform mixing. The high Atwood number PDFs of density have skewness towards the larger densities, indicating less mixing of the heavy fluid due to its inertia. The skewness in the density gradient PDFs at high Atwood number displays strong density local variations that can enhance mixing at molecular scales. Turbulent kinetic energy decreases with streamwise distance from the jet for low Atwood number but increases for high Atwood number due to larger buoyancy and density-driven shear. Over 3000 experimental realisations are used to calculate statistical characteristics of the mixing, including valuable and rarely given data such as Favre-averaged turbulent quantities: mass flux velocity, Reynolds stress, turbulent kinetic energy, and density-specific volume correlation. Buoyancy effects are observed in these quantities and the trends are compared qualitatively with direct numerical simulations.     

On hydromagnetic turbulent shear flow between moving parallel planes

Hydromagnetic turbulent shear flow of viscous, incompressible, electrically conducting fluid between two infinite uniformly porous moving parallel planes in the presence of axial and transverse magnetic field has been studied by the semi-empirical method. The expressions for the mean distributions for velocity and magnetic field for turbulent shear flow have been obtained for both the cases. As particular cases, the results have been obtained if the planes are fixed. The solutions obtained when the magnetic field is axial have been shown graphically for turbulent and laminar flows.     

Flame speed and self-similar propagation of expanding turbulent premixed flames

In this Letter we present turbulent flame speeds and their scaling from experimental measurements on constant-pressure, unity Lewis number expanding turbulent flames, propagating in nearly homogeneous isotropic turbulence in a dual-chamber, fan-stirred vessel. It is found that the normalized turbulent flame speed as a function of the average radius scales as a turbulent Reynolds number to the one-half power, where the average radius is the length scale and the thermal diffusivity is the transport property, thus showing self-similar propagation. Utilizing this dependence it is found that the turbulent flame speeds from the present expanding flames and those from the Bunsen geometry in the literature can be unified by a turbulent Reynolds number based on flame length scales using recent theoretical results obtained by spectral closure of the transformed G equation.     

Semi-analytical validation of a dynamic large-eddy simulation procedure for turbulent premixed flames via the G-equation

The performance of a dynamic subgrid model for the turbulent burning speed of a premixed flame is investigated for a series of idealized test cases where the flame front is wrinkled by a multiple-scale shear flow; a rigorous asymptotic subgrid model is also implemented for comparison. Explicit formulae for the flame wrinkled shape and turbulent speed are available to generate a reference database. The role of the subgrid wrinkling models is to achieve the same overall flame shape and propagation speed in a simulation where only the largest scales of the flow are explicitly accounted for. Very good results are obtained when the subgrid burning speed enhancement is estimated using the asymptotic subgrid model. On the other hand, the dynamic model attempts to exploit the scaling observable in the simulation to extrapolate the turbulent burning speed enhancement in the original system. The performance of this strategy is adequate for some regimes but poor for others; the source of the problem is traced back to the existence of a scaling transition that occurs as the flame propagating speed is adjusted during the large-eddy simulation. A modification to the scaling of the enhanced burning is implemented to account for the existence of the two distinct scaling ranges; it improves significantly the predictions of the dynamic model away from the transition, but results in the near-critical range remain predictably very poor compared with the rigorous asymptotic model results. These conclusions based on a priori performance for the reference steady data are confirmed by comparing unsteady large-eddy and direct simulations. Results based on rigorous mathematical tools are possible here because of the separation of length scales in the special class of idealized flow fields used in this study: their relevance to more realistic flows is also discussed.     

Thermal diffusivity measurements of spherical samples using active infrared thermography

This work presents a method for measuring the thermal diffusivity of spherical samples using active infrared thermography. The principal novelty of this method lies in the deduction of an analytical model to obtain the spatial and temporal distribution of temperature in spherical samples. The model is obtained from the classical theory of heat conduction or the 3D heat diffusion equation. In order to analyze the behavior of the model, an active infrared thermography is used in order to monitor the spatial and temporal temperature distribution. Three different materials are used as spherical samples and they are heated by radiation increasing this way its temperature. The recorded data is fitted to the model by adjusting the diffusivity parameter. The results of the diffusivity values obtained using this model are consistent with those obtained from a standard thermal properties analyzer.