On extraction of time-varying mean wind speed from wind record based on stationarity index

We have proposed in a previous study a non-stationary wind model to represent the typhoon record as a summation of a time-varying mean wind speed (TVM) and a stationary turbulence. This note further suggests a quantitative scheme, rather than the previous qualitative method, to find the best TVM for any given wind record. Trial TVMs are first extracted from the wind record by a data-processing technique named empirical mode decomposition. For each TVM, its corresponding turbulent component is computed by removing the TVM from the original wind record, and the degree of stationarity of the turbulence component is checked. The best TVM is taken as the one that leads to the maximum degree of stationarity. The degree of stationarity of turbulence is quantified by two indicators: ?? the ratio of horizontal wind variability and wind speed; and ?? the ratio of friction velocity at different Reynolds averaging periods. The applicability of the suggested scheme is validated with 550 typhoon and 3300 monsoon records of 10 minute duration and at different measurement heights. Threshold values for the two stationary indicators ?? and ?? are determined using field measurements and their sensitivities to the Reynolds averaging periods are discussed. Observations in this study demonstrate that the suggested scheme is proper for finding the TVM of a wind record. For a stationarity quantification of 10 minute duration record, the ?? indicator with 30 second Reynolds averaging period is recommended.     

Dynamical Implications of Block Averaging

We show that traditional Reynolds (block) averaging produces turbulence statistics whose time evolution is incompatible with the Navier–Stokes equation. Specifically, the zero integral scale that block averaging always produces leads to a trivial (zero-equals-zero) solution of the Navier–Stokes equation for autocovariances. We suggest alternative methods for analyzing turbulence time series that do not always generate a zero integral scale and, as a result, yield autocovariances whose time evolutions are compatible with the Navier–Stokes equation.     

The effect of idealized water waves on the turbulence structure and kinetic energy budgets in the overlying airflow

The influence of an idealized moving wavy surface on the overlying airflow is investigated using direct numerical simulations (DNS). In the present simulations, the bulk Reynolds number is Re = 8000 (; where U₀ is the forcing velocity of the flow, h the height of the domain and v the kinematic viscosity) and the phase speed of the imposed waves relative to the friction velocity, i.e., the wave age varies from very slow to fast waves. The wave signal is clearly present in the airflow up to at least 0.15λ (where λ is the wave length) and is present up to higher levels for faster waves. In the kinetic energy budgets, pressure transport is mainly of importance for slow waves. For fast waves, viscous transport and turbulent transport dominate near the surface. Kinetic energy budgets for the wave and turbulent perturbations show a non-negligible transport of turbulent kinetic energy directed from turbulence to the wave perturbation in the airflow. The wave-turbulent energy transport depends on the size, tilt, and phase of the wave-induced part of the turbulent Reynolds stresses.According to the DNS data, slow waves are more efficient in generating isotropic turbulence than fast waves.Despite the differences in wave-shape as well as in Reynolds number between the idealized direct numerical simulations and the atmosphere, there are intriguing similarities in the turbulence structure. Important information about the turbulence above waves in the atmosphere can be obtained from DNS—the data must, however, be interpreted with care.     

Comments about the drag laws used in hail growth simulations

The values of the drag coefficient (C_D), experimentally obtained by different authors for graupels and hailstones, are carefully analyzed, in order to investigate how the numerical simulations of hail trajectories, which markedly depend on the C_D values assigned to the growing particle, can be improved. It is found that new measurements of C_D should be performed in the whole range of the Reynolds number (Re) of the cloud ice particles. These measurements should be carried out applying the same criterion to define the characteristic variables (particle cross section and dimensions) used to calculate C_D and Re. Furthermore, it is found that the C_D fitting curves should be expressed as a function not only of Re but also of other parameters which represent the particle surface roughness and the departure from the spherical form. Based on this finding, a plausible definition of these variables and the way to introduce them in a hail trajectory calculation are proposed.     

Heat transport coefficients for constant energy flux models of broad leaves

Experimental determinations of the local heat transfer by forced convection from model leaves heated by a constant energy flux were made in the laboratory under laminar and turbulent flow conditions.The results are expressed in a logarithmic dimensionless plot of the local Nusselt number, Nu _d , against the local Reynolds number, Re _d . For the laminar case, Nu _d was only a linear function of Re _d ^1/2 downwind from the leading edge regions, although this relationship departed from that predicted theoretically due to the finite size and thickness of the model. For the turbulent case, a simple relationship between Nu _d and Re _d was found over a wide range of Reynolds numbers. The enhancement of heat transfer in the turbulent case depends primarily on the scale of turbulence rather than on the turbulent intensity.Past workers have discussed their results in relation to a factor , defined as the ratio between the heat transfer predicted by the Polhausen equation, and that measured. The results suggest that is not a unique parameter and may not be useful in describing the overall turbulent transfer process.     

Turbulent Flow at 190 m Height Above London During 2006–2008: A Climatology and the Applicability of Similarity Theory

Flow and turbulence above urban terrain is more complex than above rural terrain, due to the different momentum and heat transfer characteristics that are affected by the presence of buildings (e.g. pressure variations around buildings). The applicability of similarity theory (as developed over rural terrain) is tested using observations of flow from a sonic anemometer located at 190.3 m height in London, U.K. using about 6500 h of data. Turbulence statistics—dimensionless wind speed and temperature, standard deviations and correlation coefficients for momentum and heat transfer—were analysed in three ways. First, turbulence statistics were plotted as a function only of a local stability parameter z/Λ (where Λ is the local Obukhov length and z is the height above ground); the σ _i /u _* values (i = u, v, w) for neutral conditions are 2.3, 1.85 and 1.35 respectively, similar to canonical values. Second, analysis of urban mixed-layer formulations during daytime convective conditions over London was undertaken, showing that atmospheric turbulence at high altitude over large cities might not behave dissimilarly from that over rural terrain. Third, correlation coefficients for heat and momentum were analyzed with respect to local stability. The results give confidence in using the framework of local similarity for turbulence measured over London, and perhaps other cities. However, the following caveats for our data are worth noting: (i) the terrain is reasonably flat, (ii) building heights vary little over a large area, and (iii) the sensor height is above the mean roughness sublayer depth.     

Analysis of Turbulence Collapse in the Stably Stratified Surface Layer Using Direct Numerical Simulation

The nocturnal atmospheric boundary layer (ABL) poses several challenges to standard turbulence and dispersion models, since the stable stratification imposed by the radiative cooling of the ground modifies the flow turbulence in ways that are not yet completely understood. In the present work we perform direct numerical simulation of a turbulent open channel flow with a constant (cooling) heat flux imposed at the ground. This configuration provides a very simplified model for the surface layer at night. As a result of the ground cooling, the Reynolds stresses and the turbulent fluctuations near the ground re-adjust on times of the order of L/u _τ , where L is the Obukhov length scale and u _τ is the friction velocity. For relatively weak cooling turbulence survives, but when Re_L=Lu_t/n <~100{Re_L=Lu_\tau/\nu \lesssim 100} turbulence collapses, a situation that is also observed in the ABL. This criterion, which can be locally measured in the field, is justified in terms of the scale separation between the largest and smallest structures of the dynamic sublayer.     

On the structure of supercritical western boundary currents

The structure of supercritical western boundary currents is investigated using a quasi-geostrophic numerical model. The basic flow is of meridional Munk balance, and the input boundary is perturbed by the most unstable wave solution obtained from linear spatial instability calculations. Self-preserving (or equilibrium) solutions are obtained for the model runs at Re=30, 60, 90, and 120, and their energy and vorticity budgets are analyzed. In an analogy with the laboratory turbulence of wall boundary layers, the western boundary layer is divided into inner and outer layers. In the inner layer, the mean energy is dissipated via direct viscous dissipation, while in the outer layer it is converted to the eddy energy via turbulence production. The main scenario is that the mean energy is produced in the inner layer via ageostrophic pressure work divergence, and it is partly removed due to viscous action within a narrow region near the wall, defined here as viscous sub-layer. The remaining portion is converted to the eddy energy via turbulence production in the outer layer, which is in turn transported to the inner layer, then again to the viscous sub-layer where it is ultimately dissipated. In the near-wall side, the vorticity balance of the mean flow is maintained by viscous effect and Reynolds flux divergence, while in the offshore side it is maintained by beta effect and Reynolds flux divergence. The length scale of the supercritical boundary current is roughly , where L_M is the Munk length, as observed from a dimensional analysis.     

Resolved-scale turbulence in the atmospheric surface layer from a large eddy simulation

Large eddy simulation has encountered difficulties in handling turbulence in the atmospheric surface layer due to deficiencies in sub-grid scale models. This paper addresses the possibility of resolving the turbulence in the upper part of the surface layer by a low-aspect ratio of grid spacing. Results show that resolved-scale shear stresses dominate over the sub-grid scale components so that effects due to the sub-grid scale model can be ignored in this region. The effects of the lower boundary condition on the resolved-scale turbulence in the upper part of the surface layer are discussed. It is concluded that the normalized mean velocity shear and resolved turbulence in the upper part of the surface layer are not affected by the specification of the lower boundary condition. In addition, the present work proposes a new independent model parameter, the Smagorinsky Reynolds Number (Re_SM), and demonstrates that this number determines the resolved turbulence in the upper part of the surface layer.     

A comparison between two stochastic diffusion models in a complex three-dimensional flow

A Random Displacement Model (RDM) and a Langevin Equation Model (LEM) are used to simulate point releases in a complex flow around a building. The flow field is generated by a three-dimensional finite element model that uses the standardk- model to parameterize the turbulence. The RDM- and LEM-calculated concentration fields are compared, with particular emphasis on the structure in regions with high turbulence and/or recirculation. RDM and LEM results are similar qualitatively, but RDM tends to predict lower concentration levels. In part this is due to the higher early-time diffusion. However, the expected convergence at later times is prevented by the interaction of the diffusion with the strongly inhomogeneous mean flow.Notation a _i coefficient in the Langevin equation - b _ij coefficient in the Langevin equation - C ₀ the universal constant associated with the Lagrangian structure function - H building height (22.5 m) - K eddy viscosity - K _k eddy viscosity used in the definition of the off-diagonal Reynolds stresses - k turbulent kinetic energy - LEM Langevin Equation Model - p ₁ local unit vector in thexy-plane, orthogonal tos - p ₂ local unit vector, orthogonal to boths andp ₁ - RDM Random Displacement Model - s local unit vector in the streamline direction - T local decorrelation time (Lagrangian time scale) - U magnitude of the local Eulerian mean wind velocity - u _s total velocity in the streamline direction - u ₁ velocity component in thexy-plane, orthogonal to the streamline direction - u ₂ velocity component orthogonal to bothu _s andu ₁ - _i mean Eulerian wind velocity - W _i stochastic vector-valued Wiener process - x unit vector inx-direction - y unit vector iny-direction - z unit vector inz-direction - angle between thexy-plane and the mean wind streamline - angle between the projection in thexy-plane of the streamline and thex-axis - _ij the Kronecker delta function - rate of turbulence dissipation - _i/g_a the part ofa _i that contains mean wind and turbulence gradients - _ij inverse of a Reynolds stress tensor component - _ij shorthand for a quantity that defines a part of _i/g_a - _i shorthand for a quantity that defines a part of _i/g_a - _ij Reynolds stress tensor component     

Dispersion in sheared Gaussian homogeneous turbulence

By integrating the Fokker-Planck equation corresponding to a Lagrangian stochastic trajectory model, which is consitent with the selection criterion of Thomson (1987), an analytical solution is given for the joint probability density functionp(x_i, u_i, t) for the position (x _i) and velocity (u _i) at timet of a neutral particle released into linearly-sheared, homogeneous turbulence. The solution is compared with dispersion experiments conforming to the restrictions of the model and with a shortrange experiment performed in highly inhomogeneous turbulence within and above a model crop canopy. When the turbulence intensity, wind shear and covariance are strong, the present solution is better than simpler solutions (Taylor, 1921; Durbin, 1983) and as good as any numerical Lagrangian stochastic model yet reported.     

Studies of forced-convection heat and mass transfer of fluttering realistic leaf models

The forced-convection mass transfer - and by analogy, heat transfer - of various realistic leaf models at Reynolds numbers 2 x 10³<Re<4 x 104 was studied with an electrochemical method. The results are compared with similar measurements on plates and with transfer coefficients calculated from the laminar boundary-layer theory. In this way the validity of the commonly-used analytical expressions which represent the leaf by a rigid plate and neglect the effects of leaf curvature, fluttering, surface roughness and fluid turbulence, can be tested.The measurements show that for fluttering single leaves, the convective mass-transfer coefficients must be expected to be higher by a factor of 1.4 ± 0.1 than the ones calculated for rigid plates of equal size and shape. For a leaf in a crop, the increase might be as high as a factor of 2. The high transfer coefficients measured for elements of cedar foliage are also discussed.     

A numerical study of the convective boundary layer

Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The turbulent fields are obtained directly by solving the three-dimensional time-dependent Navier-Stokes equations (using the Boussinesq approximation to account for buoyancy effects), and no turbulence model is needed. Two levels of heating are considered, one quite vigorous, the other more moderate. Statistics for the vigorously heated case are found to agree reasonably well with laboratory, field, and large-eddy simulation results, when Deardorff's mixed-layer scaling is used. No indication of large-scale longitudinal roll cells is found in this convection-dominated flow, for which the inversion height to Obukhov length scale ratio –z _i /L _*=26. However, when heating is more moderate (so that –z _i /L _*=2), evidence of coherent rolls is present. About 10% of the total turbulent kinetic energy and turbulent heat flux, and 20% of the Reynolds shear stress, are estimated to be a direct consequence of the observed cells.     

Numerical Modelling of the Turbulent Flow Developing Within and Over a 3-D Building Array, Part I: A High-Resolution Reynolds-Averaged Navier—Stokes Approach

A study of the neutrally-stratified flow within and over an array of three-dimensional buildings (cubes) was undertaken using simple Reynolds-averaged Navier—Stokes (RANS) flow models. These models consist of a general solution of the ensemble-averaged, steady-state, three-dimensional Navier—Stokes equations, where the k-ε turbulence model (k is turbulence kinetic energy and ε is viscous dissipation rate) has been used to close the system of equations. Two turbulence closure models were tested, namely, the standard and Kato—Launder k-ε models. The latter model is a modified k-ε model designed specifically to overcome the stagnation point anomaly in flows past a bluff body where the standard k-ε model overpredicts the production of turbulence kinetic energy near the stagnation point. Results of a detailed comparison between a wind-tunnel experiment and the RANS flow model predictions are presented. More specifically, vertical profiles of the predicted mean streamwise velocity, mean vertical velocity, and turbulence kinetic energy at a number of streamwise locations that extend from the impingement zone upstream of the array, through the array interior, to the exit region downstream of the array are presented and compared to those measured in the wind-tunnel experiment. Generally, the numerical predictions show good agreement for the mean flow velocities. The turbulence kinetic energy was underestimated by the two different closure models. After validation, the results of the high-resolution RANS flow model predictions were used to diagnose the dispersive stress, within and above the building array. The importance of dispersive stresses, which arise from point-to-point variations in the mean flow field, relative to the spatially-averaged Reynolds stresses are assessed for the building array.     

Direct numerical simulation of a vigorously heated low Reynolds number convective boundary layer

Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The results are obtained by directly solving the three-dimensional time-dependen Navier-Stokes equations with the Boussinesq buoyancy approximation, resolving all relevant scales of motion (no turbulence closure is needed). The flow is capped by a stable temperature inversion and heated from below at a rate that produces an inversion-height to Obukhov-length ratio −z_i/L_* = 32. Temperature and velocity variance profiles are found to agree well with those from an earlier vigorously heated under-resolved computation at higher Reynolds number, and with experimental data of Deardorff and Willis (Boundary-Layer Meteorol., 32: 205–236, 1985). Significant helicity is found in the layer, and helical convection patterns of the scale of the inversion height are observed.     

Application of a Large-Eddy Simulation Database to Optimisation of First-Order Closures for Neutral and Stably Stratified Boundary Layers

Large-eddy simulation (LES) is a well-established numerical technique, resolving the most energetic turbulent fluctuations in the planetary boundary layer. By averaging these fluctuations, high-quality profiles of mean quantities and turbulence statistics can be obtained in experiments with well-defined initial and boundary conditions. Hence, LES data can be beneficial for assessment and optimisation of turbulence closure schemes. A database of 80 LES runs (DATABASE64) for neutral and stably stratified planetary boundary layers (PBLs) is applied in this study to optimize first-order turbulence closure (FOC). Approximations for the mixing length scale and stability correction functions have been made to minimise a relative root-mean-square error over the entire database. New stability functions have correct asymptotes describing regimes of strong and weak mixing found in theoretical approaches, atmospheric observations and LES. The correct asymptotes exclude the need for a critical Richardson number in the FOC formulation. Further, we analysed the FOC quality as functions of the integral PBL stability and the vertical model resolution. We show that the FOC is never perfect because the turbulence in the upper half of the PBL is not generated by the local vertical gradients. Accordingly, the parameterised and LES-based fluxes decorrelate in the upper PBL. With this imperfection in mind, we show that there is no systematic quality deterioration of the FOC in the strongly stable PBL provided that the vertical model resolution is better than 10 levels within the PBL. In agreement with previous studies, we found that the quality improves slowly with the vertical resolution refinement, though it is generally wise not to overstretch the mesh in the lowest 500 m of the atmosphere where the observed, simulated and theoretically predicted stably stratified PBL is mostly located. The submission to a special issue of the “Boundary-Layer Meteorology” devoted to the NATO advanced research workshop “Atmospheric Boundary Layers: Modelling and Applications for Environmental Security”.     

New Approaches in Two-equation Turbulence Modelling for Atmospheric Applications

The turbulence closure in atmospheric boundary-layer modelling utilizing Reynolds Averaged Navier–Stokes (RANS) equations at mesoscale as well as at local scale is lacking today a common approach. The standard kɛ model, although it has been successful for local scale problems especially in neutral conditions, is deficient for mesoscale flows without modifications. The k–ɛ model is re-examined and a new general approach in developing two-equation turbulence models is proposed with the aim of improving their reliability and consequently their range of applicability. This exercise has led to the replacement of the ɛ-transport equation by the transport equation for the turbulence inverse length scale (wavenumber). The present version of the model is restricted to neutrally stratified flows but applicable to both local scale and mesoscale flows. The model capabilities are demonstrated by application to a series of one-dimensional planetary boundary-layer problems and a two-dimensional flow over a square obstacle. For those applications, the present model gave considerably better results than the standard k–ɛ model.     

Stationarity of Turbulence in Light Winds during the Maritime Continent Thunderstorm Experiment

Non-stationary turbulence can invalidate eddy flux calculations. Two-hour longrecords of wind velocity, temperature and humidity are classified as stationaryor non-stationary based on the behaviour of the flux as a function of Reynoldsaveraging period; a number of indicators of stationarity are investigated. Thetwo-hour Maritime Continent Thunderstorm Experiment wind datasets are notcompletely stationary, as indicated by the lack of a spectral gap, but can beclassified as approximately stationary, when the mean wind speed is greaterthan 3.8 m s^-1, or the standard deviation of true wind direction is 10°, or the ratio of horizontal wind variance to wind speed is <0.25. In the stationary case the calculated friction velocity exhibits a 7% decrease on average when the Reynolds averaging period is doubled, while data classified as non-stationaryexhibit an increase of 32%. There is little non-stationary behaviour in the kinematicheat fluxes, and is independent of the non-stationarity of the friction velocity. Thekinematic heat fluxes show small decreases (around 3%) when the Reynolds averagingperiod is doubled.     

Direct Numerical Simulation of Stable Channel Flow at Large Stability

We consider a model for the stable atmospheric boundary at large stability, i.e. near the limit where turbulence is no longer able to survive. The model is a plane horizontally homogeneous channel flow, which is driven by a constant pressure gradient and which has a no-slip wall at the bottom and a free-slip wall at the top. At the lower wall a constant negative temperature flux is imposed. First, we consider a direct numerical simulation of the same channel flow. The simulation is computed with the neutral channel flow as initial condition and computed as a function of time for various values of the stability parameter h/L, where h is the channel height and L is related to the Obukhov length. We find that a turbulent solution is only possible for h/L < 1.25 and for larger values turbulence decays. Next, we consider a theoretical model for this channel flow based on a simple gradient transfer closure. The resulting equations allow an exact solution for the case of a stationary flow. The velocity profile for this solution is almost linear as a function of height in most of the channel. In the limit of infinite Reynolds number, the temperature profile has a logarithmic singularity at the upper wall of the channel. For the cases where a turbulent flow is maintained in the numerical simulation, we find that the velocity and temperature profiles are in good agreement with the results of the theoretical model when the effects of the surface layer on the exchange coefficients are taken into account. Frans Nieuwstadt, a recently retired member of the BLM Editorial Board and a well-known member of the boundary-layer/turbulence community, died unexpectedly on 18 May 2005. An obituary will appear in a later issue of BLM.     

Intermittency of Turbulence in the Atmospheric Boundary Layer: Scaling Exponents and Stratification Influence

We show the relationship between the intermittency of turbulence and the type of stratification for different atmospheric situations during the SABLES98 field campaign. With this objective, we first demonstrate the scaling behaviour of the velocity structure functions corresponding to these situations; next, we analyze the curvature of the scaling exponents of the velocity structure functions versus the order of these functions (ζ _p vs. p), where ζ _p are the exponents of the power relation for the velocity structure function with respect to the scale. It can be proved that this curve must be concave, under the assumption that the incompressible approximation does not break down at high Reynolds numbers. The physical significance of this kind of curvature is that the energy dissipation rate increases as the scale of the turbulent eddies diminishes (intermittency in the usual sense). However, the constraints imposed by stability, preventing full development of the turbulence, allow the function ζ _p versus p to show any type of curvature. In this case, waves of high frequency trapped by the stability, or bursts of turbulence caused by the breaking up of internal waves, may produce a redistribution of energy throughout the scaling range. Due to this redistribution, the variation with the scale of the energy dissipation rate may be smaller (decreasing the intermittency) and, even in more stable situations, this rate may diminish (instead of increasing) as the scale diminishes (convex form of the curve ζ _p vs. p).