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.     

An explicit algebraic model of the planetary boundary layer turbulence: test computation of the neutrally stratified atmospheric boundary layer

An explicit algebraic model of Reynolds stresses and the turbulent heat flux vector for the planetary boundary layer in a neutrally stratified boundary layer of the atmosphere above a homogeneous rough surface is tested. The version of the algebraic model under consideration is constructed on the physical principles of the RANS (Reynolds-averaged Navier?Stokes) approximation for describing stratified turbulence, it employs three forecasting equations, and a correct reproduction of the main characteristics of a neutral atmospheric boundary layer — the components of the mean wind velocity, the wind turn angle, and the turbulent statistics is shown. Test computations show that the proposed model may be used for goal-oriented investigations of the atmospheric boundary layer.     

Employing Taylor and Heisenberg subfilter viscosities to simulate turbulent statistics in LES models

A turbulent subfilter viscosity for Large Eddy Simulation (LES) based on the Taylor statistical diffusion theory is proposed. This viscosity is described in terms of a velocity variance and a time scale, both associated to the inertial subrange. This new subfilter viscosity contains a cutoff wavenumber k_c, presenting an identical form (differing by a constant) to the Heisenberg subfilter viscosity. Therefore, both subfilter viscosities are described in terms of a sharp division between large and small wavenumbers of a turbulent flow and, henceforth, Taylor and Heisenberg subfilter viscosities are in agreement with the sharp Fourier filtering operation, frequently employed in LES models. Turbulent statistics of different orders, generated from atmospheric boundary layer simulations employing both Taylor and Heisenberg subfilter viscosities have been compared with observations and results provided by other simulations. The comparison shows that the LES model utilizing the approaches of Taylor and Heisenberg reproduces these turbulent statistics correctly in different vertical regions of a planetary convective boundary layer (CBL).     

Studies of the flow and turbulence fields in a turbulent pulsed jet flame using LES/PDF

A turbulent piloted jet flame subject to a rapid velocity pulse in its fuel jet inflow is proposed as a new benchmark case for the study of turbulent combustion models. In this work, we perform modelling studies of this turbulent pulsed jet flame and focus on the predictions of its flow and turbulence fields. An advanced modelling strategy combining the large eddy simulation (LES) and the probability density function (PDF) methods is employed to model the turbulent pulsed jet flame. Characteristics of the velocity measurements are analysed to produce a time-dependent inflow condition that can be fed into the simulations. The effect of the uncertainty in the inflow turbulence intensity is investigated and is found to be very small. A method of specifying the inflow turbulence boundary condition for the simulations of the pulsed jet flame is assessed. The strategies for validating LES of statistically transient flames are discussed, and a new framework is developed consisting of different averaging strategies and a bootstrap method for constructing confidence intervals. Parametric studies are performed to examine the sensitivity of the predictions of the flow and turbulence fields to model and numerical parameters. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed modelling methods.     

一种改进后的亚格子特征尺度表达式

针对大涡模拟涡粘性亚格子模式中使用网格尺度为特征尺度存在的不足,提出一种改进的亚格子特征尺度表达式,利用结合Sagaut混合尺度模式获得的改进模式与Smagorinsky模式、Germano模式和Sagaut模式等进行对比研究.时间发展混合层和中性大气边界层的大涡模拟结果表明,改进的有效亚格子特征尺度能合理反映亚格子脉动的时空分布特性,同时扩展了适用范围.在对耗散程度、流动演化过程、可解湍流强度和可解雷诺应力等方面的模拟,改进模式优于基准模式,表明亚格子脉动有效特征尺度在大涡模拟中具有一定的理论基础和广泛的应用价值.     

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.     

Investigating asymptotic suction boundary layers using a one-dimensional stochastic turbulence model

The turbulent asymptotic suction boundary layer is studied using a one-dimensional turbulence (ODT) model. ODT is a fully resolved, unsteady stochastic simulation technique. While flow properties reside on a one-dimensional domain, turbulent advection is represented using mapping events whose occurrences are governed by a random process. Due to its reduced spatial dimensionality, ODT achieves major cost reductions compared to three-dimensional (3D) simulations. A comparison to recent direct numerical simulation (DNS) data at moderate Reynolds number (Re = u_∞ / v₀ = 333, where u_∞ and v₀ are the free stream and suction velocity, respectively) suggests that the ODT model is capable of reproducing several velocity statistics, i.e. mean velocity and turbulent kinetic energy budgets, while peak turbulent stresses are under-estimated by ODT. Variation of the Reynolds number in the range Re ∈ [333,400,500,1000] shows that ODT can reproduce various trends observed as a result of increased suction in turbulent asymptotic suction boundary layers, i.e. the reduction of Reynolds stresses and enhanced skin friction. While up to Re = 500 our results can be directly compared to recent LES data, the simulation at Re = 1000 is currently not feasible through full 3D simulations, hence ODT may assist the design of future DNS or LES simulations at larger Reynolds numbers.     

Analysis of the turbulent boundary layer in the vicinity of a self-excited cylindrical Helmholtz resonator

This study investigates the changes in the structure of a turbulent boundary layer downstream of a flow-excited Helmholtz resonator. To this end, a fully developed turbulent boundary layer over a resonator mounted flush with a flat plate was simulated by implementing a large eddy simulation (LES). To assist in understanding the effect of the resonator on the flow structure, a sensitivity study was undertaken by changing the main geometrical parameters of the resonator. The results demonstrated that when the boundary layer thickness equals the orifice length, the cross-stream component of velocity fluctuations penetrates the boundary layer, resulting in a reduction of the turbulence intensity by up to 12%. Therefore, it is concluded that a Helmholtz resonator has the potential to reduce the instabilities within the boundary layer. These investigations also assist in identifying the optimal parameters to delay turbulence events within the grazing flow using Helmholtz resonators.     

Building resolving large-eddy simulations and comparison with wind tunnel experiments

We perform large-eddy simulations (LES) of the flow past a scale model of a complex building. Calculations are accomplished using two different methods to represent the edifice. The first method employs the standard Gal-Chen and Somerville terrain-following coordinate transformation, common in mesoscale atmospheric simulations. The second method uses an immersed boundary approach, in which fictitious body forces in the equations of motion are used to represent the building by attenuating the flow to stagnation within a time comparable to the time step of the model. Both methods are implemented in the same hydrodynamical code (EULAG) using the same nonoscillatory forward-in-time (NFT) incompressible flow solver based on the multidimensional positive definite advection transport algorithms (MPDATA). The two solution methods are compared to wind tunnel data collected for neutral stratification. Profiles of the first- and second-order moments at various locations around the model building show good agreement with the wind tunnel data. Although both methods appear to be viable tools for LES of urban flows, the immersed boundary approach is computationally more efficient. The results of these simulations demonstrate that, contrary to popular opinion, continuous mappings such as the Gal-Chen and Somerville transformation are not inherently limited to gentle slopes. Calculations for a strongly stratified case are also presented to point out the substantial differences from the neutral boundary layer flows.     

The compressible hybrid RANS/LES formulation using an additive operator

In this work, the compressible governing equations for the hybrid Reynolds-averaged/large-eddy simulations are formally derived by applying a hybrid filter to the Navier–Stokes equations. This filter is constructed by linearly combining the Reynolds-average (RANS) and large-eddy simulation (LES) operators with a continuous blending function. The derived hybrid equations include additional terms that represent the interactions between RANS and LES formulations. The relevance of these terms is investigated in flat-plate turbulent boundary layer simulations and indicate that these additional terms play a fundamental role in compensating for the turbulence that is neither modeled nor resolved in the transition region between RANS and LES. Results also show that when the additional terms are included, the calculations are not very sensitive to the blending function implemented in the hybrid filter. In the contrary, when these terms are neglected and a step function is implemented in the hybrid filter, nonphysical discontinuities are predicted in the flow statistics.     

超声速层流/湍流压缩拐角流动结构的实验研究

在Ma=3.0的超声速风洞中, 分别对上游边界层为超声速层流和湍流, 压缩角度为25°和28°的压缩拐角流动进行了实验研究. 采用纳米粒子示踪平面激光散射(NPLS)技术获得了流场整体和局部区域的精细结构, 边界层、剪切层、分离激波、回流区和再附激波等典型结构清晰可见, 测量了超声速层流压缩拐角壁面的压力系数. 从时间平均的流场结构中测量出分离激波、再附激波的角度和再附后重新发展的边界层的增长情况, 通过分析时间相关的流场NPLS图像, 可以发现流场结构随时间的演化特性. 实验结果表明: 在25°的压缩角度下, 超声速层流压缩拐角流动发生了典型的分离, 边界层迅速增长失稳转捩, 并引起一道诱导激波, 流场中出现了K-H涡、剪切层和微弱压缩波结构, 而超声速湍流压缩拐角流动没有出现分离, 湍流边界层始终表现为附着状态; 在28° 的压缩角度下, 超声速层流压缩拐角流动进一步分离, 回流区范围明显扩大, 诱导激波、分离激波向上游移动, 再附激波向下游移动, 分离区流动结构复杂, 相比之下, 超声速湍流压缩拐角流动的回流区范围明显较小, 边界层增长缓慢, 流场中没有出现诱导激波、K-H涡和压缩波, 流动分离区域的结构也相对简单, 但分离激波的强度则明显更强. 关键词： 压缩拐角 层流 湍流 流动结构     

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.     

The prediction of broadband noise from wind turbines

This paper describes a broadband noise prediction scheme for wind turbines. The source mechanisms included in the method are unsteady lift noise, unsteady thickness noise, trailing edge noise and the noise from separated flow. Special methods have been developed to model the inflow turbulence from the atmospheric boundary layer and acoustic radiation to the geometric near field of the rotor. Predictions are compared with measurements on 20 m and 80 m diameter wind turbines. The results show that the turbulence length scale in the atmospheric boundary layer is too large to give the measured noise levels. Very good agreement is obtained between predictions and measurements if the turbulence length scale is taken to be equal to the blade chord.     

Transition region where the large-scale and very large scale motions coexist in atmospheric surface layer: wind tunnel investigation

The turbulent structures in atmospheric surface layer (ASL) are investigated in wind tunnel with hot-wire anemometers in present study. The results show that there exist two typical turbulent structures characterised by their streamwise length scales, i.e. large-scale motions (LSMs) and very large scale motions (VLSMs) as reported recently in pipe flow, channel flow, zero-pressure-gradient turbulent boundary layers and near-neutral ASLs. Moreover, a transition region containing both LSMs and VLSMs is first reported in present investigation, and this region depends on the Reynolds numbers. In the transition region, as the wall-normal distance increases, the turbulent energy ratio of LSMs is gradually weakened but VSLMs strengthened.     

Correlation dimension of the flow and spatial development of dynamical chaos in a boundary layer

The relation between the spatial development of turbulence in a real hydrodynamical system (a boundary layer on a plate in a wind tunnel) and the increase of dimension downstream is established experimentally.     

Unsteady coupling of Navier-Stokes and radiative heat transfer solvers applied to an anisothermal multicomponent turbulent channel flow

Direct numerical simulations (DNS) of an anisothermal reacting turbulent channel flow with and without radiative source terms have been performed to study the influence of the radiative heat transfer on the optically non-homogeneous boundary layer structure. A methodology for the study of the emitting/absorbing turbulent boundary layer (TBL) is presented. Details on the coupling strategy and the parallelization techniques are exposed. An analysis of the first order statistics is then carried out. It is shown that, in the studied configuration, the global structure of the thermal boundary layer is not significantly modified by radiation. However, the radiative transfer mechanism is not negligible and contributes to the heat losses at the walls. The classical law-of-the-wall for temperature can thus be improved for RANS/LES simulations taking into account the radiative contribution.     

时变来流条件下大气边界层流动大涡模拟

采用大涡模拟方法,利用平坦地面模型和带楔形结构的粗糙地面模型研究近地面空间风场特征.结果表明,在时变来流条件下,平坦地面模型和粗糙地面模型均在阵风峰期产生下沉气流,在阵风谷期产生上升气流;垂直气流大小与下垫面粗糙度和水平来流风速变化幅度有关.大气边界层内风场的风速大小和方向总是随时间变化的,与稳定来流假设相比,随时间变化的来流条件下边界层的流动更有利于边界层内物质和能量的输运.     

Simulations of turbulent asymptotic suction boundary layers

A series of large-eddy simulations of a turbulent asymptotic suction boundary layer (TASBL) was performed in a periodic domain, on which uniform suction was applied over a flat plate. Three Reynolds numbers (defined as ratio of free-stream and suction velocity) of Re = 333, 400 and 500 and a variety of domain sizes were considered in temporal simulations in order to investigate the turbulence statistics, the importance of the computational domain size, the arising flow structures as well as temporal development length required to achieve the asymptotic state. The effect of these two important parameters was assessed in terms of their influence on integral quantities, mean velocity, Reynolds stresses, higher order statistics, amplitude modulation and spectral maps. While the near-wall region up to the buffer region appears to scale irrespective of Re and domain size, the parameters of the logarithmic law (i.e. von Kármán and additive coefficient) decrease with increasing Re, while the wake strength decreases with increasing spanwise domain size and vanishes entirely once the spanwise domain size exceeds approximately two boundary-layer thicknesses irrespective of Re. The wake strength also reduces with increasing simulation time. The asymptotic state of the TASBL is characterised by surprisingly large friction Reynolds numbers and inherits features of wall turbulence at numerically high Re. Compared to a turbulent boundary layer (TBL) or a channel flow without suction, the components of the Reynolds-stress tensor are overall reduced, but exhibit a logarithmic increase with decreasing suction rates, i.e. increasing Re. At the same time, the anisotropy is increased compared to canonical wall-bounded flows without suction. The reduced amplitudes in turbulence quantities are discussed in light of the amplitude modulation due to the weakened larger outer structures. The inner peak in the spectral maps is shifted to higher wavelength and the strength of the outer peak is much less than for TBLs. An additional spatial simulation was performed, in order to relate the simulation results to wind tunnel experiments, which – in accordance with the results from the temporal simulation – indicate that a truly TASBL is practically impossible to realise in a wind tunnel. Our unique data set agrees qualitatively with existing literature results for both numerical and experimental studies, and at the same time sheds light on the fact why the asymptotic state could not be established in a wind tunnel experiment, viz. because experimental studies resemble our simulation results from too small simulation boxes or insufficient development times.     

Generalised scale-by-scale energy-budget equations and large-eddy simulations of anisotropic scalar turbulence at various Schmidt numbers

A necessary condition for the accurate prediction of turbulent flows using large-eddy simulation (LES) is the correct representation of energy transfer between the different scales of turbulence in the LES. For scalar turbulence, transfer of energy between turbulent length scales is described by a transport equation for the second moment of the scalar increment. For homogeneous isotropic turbulence, the underlying equation is the well-known Yaglom equation. In the present work, we study the turbulent mixing of a passive scalar with an imposed mean gradient by homogeneous isotropic turbulence. Both direct numerical simulations (DNS) and LES are performed for this configuration at various Schmidt numbers, ranging from 0.11 to 5.56. As the assumptions made in the derivation of the Yaglom equation are violated for the case considered here, a generalised Yaglom equation accounting for anisotropic effects, induced by the mean gradient, is derived in this work. This equation can be interpreted as a scale-by-scale energy-budget equation, as it relates at a certain scale r terms representing the production, turbulent transport, diffusive transport and dissipation of scalar energy. The equation is evaluated for the conducted DNS, followed by a discussion of physical effects present at different scales for various Schmidt numbers. For an analysis of the energy transfer in LES, a generalised Yaglom equation for the second moment of the filtered scalar increment is derived. In this equation, new terms appear due to the interaction between resolved and unresolved scales. In an a-priori test, this filtered energy-budget equation is evaluated by means of explicitly filtered DNS data. In addition, LES calculations of the same configuration are performed, and the energy budget as well as the different terms are thereby analysed in an a-posteriori test. It is shown that LES using an eddy viscosity model is able to fulfil the generalised filtered Yaglom equation for the present configuration. Further, the dependence of the terms appearing in the filtered energy-budget equation on varying Schmidt numbers is discussed.     

A case study of dispersion in the lee of an obstacle

Summary The influence of hills, valleys and other obstacles on the atmospheric boundary layer is such that in their lee the mean wind velocity is reduced. Intermittent or fully separated flow conditions can arise, including regions of negligible mean wind. In these conditions dispersion is driven mainly by turbulent fluctuations, which prevail over the mean motion. In this work we review the general features of mean flow and turbulence intensity in the lee of a 2D obstacle. Using published wind tunnel data, we examine a case study of dispersion from a source placed in the lee and compare the measured concentrations with those of a random walk model simulation. Paper presented at the GNFAO/EURASAP Meeting, Turin, September 1989. To speed up publication, proofs were not sent to the authors and were supervised by the Scientific Committee.