Cumulus clouds and convective boundary layers: a tropical perspective on two turbulent shear flows

This paper is a review prepared for the second Marseille Colloquium on the mechanics of turbulence, held in 2011, 50 years after the first. The review covers recent developments in our understanding of the large-scale dynamics of cumulus cloud flows and of the atmospheric boundary layer in the low-wind convective regime that is often encountered in the tropics. It has recently been shown that a variety of cumulus cloud forms and life cycles can be experimentally realized in the laboratory, with the transient diabatic plume taken as the flow model for a cumulus cloud. The plume is subjected to diabatic heating scaled to be dynamically similar to heat release from phase changes in clouds. The experiments are complemented by exact numerical solutions of the Navier-Stokes-Boussinesq equations for plumes with scaled off-source heating. The results show that the Taylor entrainment coefficient first increases with heating, reaches a positive maximum and then drops rapidly to zero or even negative values. This reduction in entrainment is a consequence of structural changes in the flow, smoothing out the convoluted boundaries in the non-diabatic plume, including the tongues engulfing the ambient flow. This is accompanied by a greater degree of mixedness in the core flow because of lower dilution by the ambient fluid. The cloud forms generated depend strongly on the history of the diabatic heating profile in the vertical direction. The striking effects of heating on the flow are attributable to the operation of the baroclinic torque due to the temperature field. The mean baroclinic torque is shown to peak around a quasi-cylindrical sheet situated midway between the axis of the flow and the edges. This torque is shear-enhancing and folds down the engulfment tongues. The increase in mixedness can be traced to an explosive growth in the enstrophy, triggered by a strong fluctuating baroclinic torque that acts as a source, especially at the higher wave numbers, thus enhancing the mixedness. In convective boundary layers field measurements show that, under conditions prevailing in the tropics, the eddy fluxes of momentum and energy do not follow the Monin-Obukhov similarity. Instead, the eddy momentum flux is found to be linear in the wind speed at low winds; and the eddy heat flux is, to a first approximation, governed by free convection laws, with wind acting as a small perturbation on a regime of free convection. A new boundary layer code, based on heat flux scaling rather than wall-stress scaling, shows promising improvements in predictive skills of a general circulation model.

On the behavior of pliable plate dynamics in wind : application to vertical axis wind turbines

Numerous studies have shown that flexible materials improve resilience and durability of a structure. Several studies have investigated the behavior of elastic plates under the influence of a free stream, such as studies of the fluttering flag and others of shape reconfiguration, due to a free stream.

The principle engineering contribution of this thesis is the design and development of a vertical axis wind turbine that features pliable blades which undergo various modes of behavior, ultimately leading to rotational propulsion of the turbine. The wind turbine design was tested in a wind tunnel and at the Caltech Laboratory for Optimized Wind Energy. Ultimately, the flexible blade vertical axis wind turbine proved to be an effective way of harnessing the power of the wind.

In addition, this body of work builds on the current knowledge of elastic cantilever plates in a free stream flow by investigating the inverted flag. While previous studies have focused on the fluid structure interaction of a free stream on elastic cantilever plates, none had studied the plate configuration where the trailing edge was clamped, leaving the leading edge free to move. Furthermore, the studies presented in this thesis establish the geometric boundaries of where the large-amplitude flapping occurs.

Dynamic Stall on Vertical Axis Wind Turbines

In this study the dynamics of flow over the blades of vertical axis wind turbines was investigated using a simplified periodic motion to uncover the fundamental flow physics and provide insight into the design of more efficient turbines. Time-resolved, two-dimensional velocity measurements were made with particle image velocimetry on a wing undergoing pitching and surging motion to mimic the flow on a turbine blade in a non-rotating frame. Dynamic stall prior to maximum angle of attack and a leading edge vortex development were identified in the phase-averaged flow field and captured by a simple model with five modes, including the first two harmonics of the pitch/surge frequency identified using the dynamic mode decomposition. Analysis of these modes identified vortical structures corresponding to both frequencies that led the separation and reattachment processes, while their phase relationship determined the evolution of the flow.

Detailed analysis of the leading edge vortex found multiple regimes of vortex development coupled to the time-varying flow field on the airfoil. The vortex was shown to grow on the airfoil for four convection times, before shedding and causing dynamic stall in agreement with 'optimal' vortex formation theory. Vortex shedding from the trailing edge was identified from instantaneous velocity fields prior to separation. This shedding was found to be in agreement with classical Strouhal frequency scaling and was removed by phase averaging, which indicates that it is not exactly coupled to the phase of the airfoil motion.

The flow field over an airfoil undergoing solely pitch motion was shown to develop similarly to the pitch/surge motion; however, flow separation took place earlier, corresponding to the earlier formation of the leading edge vortex. A similar reduced-order model to the pitch/surge case was developed, with similar vortical structures leading separation and reattachment; however, the relative phase lead of the separation mode, corresponding to earlier separation, necessitated that a third frequency to be incorporated into the reattachment mode to provide a relative lag in reattachment.

Finally, the results are returned to the rotating frame and the effects of each flow phenomena on the turbine are estimated, suggesting kinematic criteria for the design of improved turbines.

Optimizing the energy output of vertical axis wind turbines for fluctuating wind conditions

The control of a wind turbine to the mean wind speed in a gusty wind results in very poor performance. Fluctuations in wind speed with time constants shorter than the response time of a wind turbine results in operation away from optimum design conditions. The effectiveness of a turbine operating in a gusty wind is shown though the use of an unsteady performance coefficient, C e. This performance coefficient is similar in form to a power coefficient. However in order to accommodate unsteady effects, Ce is defined as a ratio of energy extracted to the total wind energy available over a set time period. The turbine's response to real wind data is modelled, in the first instance, by assuming a constant rotational speed operation. It is shown that a significant increase in energy production can be realized by demanding a Tip Speed Ratio above the steady state optimum. The constant speed model is then further extended to incorporate inertial and controller effects. Parameters dictating how well a turbine can track a demand in Tip Speed Ratio have been identified and combined, to form a non-dimensional turbine response parameter. This parameter characterizes a turbine's ability to track a demand in Tip Speed Ratio dependent on an effective gust frequency. A significant increase in energy output of 42% and 245% is illustrated through the application of this over-speed control. This is for the constant rotational speed and Tip Speed Ratio feedback models respectively. The affect of airfoil choice on energy extraction within a gusty wind has been considered. The adaptive control logic developed enables the application of airfoils demonstrating high maximum L/D values but sharp stalling characteristics to be successfully used in a VAWT design.

Turbulent flow structure in experimental laboratory wind-generated gravity waves

This paper is the third part of a report on systematic measurements and analyses of wind-generated water waves in a laboratory environment. The results of the measurements of the turbulent flow on the water side are presented here, the details of which include the turbulence structure, the correlation functions, and the length and velocity scales. It shows that the mean turbulent velocity profiles are logarithmic, and the flows are hydraulically rough. The friction velocity in the water boundary layer is an order of magnitude smaller than that in the wind boundary layer. The level of turbulence is enhanced immediately beneath the water surface due to micro-breaking, which reflects that the Reynolds shear stress is of the order u *w 2. The vertical velocities of the turbulence are related to the relevant velocity scale at the still-water level. The autocorrelation function in the vertical direction shows features of typical anisotropic turbulence comprising a large range of wavelengths. The ratio between the microscale and macroscale can be expressed as λ/Λ=a Re Λ n, with the exponent n slightly different from -1/2, which is the value when turbulence production and dissipation are in balance. On the basis of the wavelength and turbulent velocity, the free-surface flows in the present experiments fall into the wavy free-surface flow regime. The integral turbulent scale on the water side alone underestimates the degree of disturbance at the free surface. © 2012 Elsevier B.V.

Placement of Small Vertical Axis Wind Turbine to Maximize Power Generation Influenced by Architectural and Geographic Interfaces in Urban Areas

Current methods for large-scale wind collection are unviable in urban areas. In order to investigate the feasibility of generating power from winds in these environments, we sought to optimize placements of small vertical-axis wind turbines in areas of artificially-generated winds. We explored both vehicular transportation and architecture as sources of artificial wind, using a combination of anemometer arrays, global positioning system (GPS), and weather report data. We determined that transportation-generated winds were not significant enough for turbine implementation. In addition, safety and administrative concerns restricted the implementation of said wind turbines along roadways for transportation-generated wind collection. Wind measurements from our architecture collection were applied in models that can help predict other similar areas with artificial wind, as well as the optimal placement of a wind turbine in those areas.

Horizontal and vertical distribution of cirripede cyprid larvae in an upwelling system off the Portuguese coast

The distribution of cirripede cyprids in relation to associated oceanographic conditions was obtained from a grid survey and intensive vertical sampling at a fixed station located 21 km off the northwest Portuguese coast in May 2002. Analysis of cyprid length composition allowed separation of 3 species groups. Chthamalus montagui, Pollicipes pollicipes and Balanus perforatus were largely restricted to the neuston layer and showed only low-amplitude vertical migration. Most C. stellatus cyprids only appeared in the upper 20 m at night, a migration which did not appear to be affected by physical conditions in the water column, but some differences in the vertical migration pattern between days were probably related to varying light penetration. C. montagui is the most abundant adult species found along the Portuguese coast, but C. stellatus cyprids, at densities of up to 8.7 ind. m–3, were the most common sampled in all depth strata at the fixed station. Cyprid horizontal distribution was mainly restricted to an offshore band along the inner shelf, where highest densities were 11 to 15 ind. m–3. This distribution pattern was considered to result from upwelling-favourable wind conditions, creating fronts along the shelf in which the cyprids become concentrated. Cyprid vertical migration, in association with current vertical shear and onshore movement of fronts during upwelling-relaxation periods, may be the mechanisms returning cyprids to the coast to settle. The regularity of these events in the region falls within the period of cyprid viability.

On the generation mechanisms of short-scale unbalanced modes in rotating two-layer flows with vertical shear

We report on the results of a laboratory investigation using a rotating two-layer annulus experiment, which exhibits both large-scale vortical modes and short-scale divergent modes. A sophisticated visualization method allows us to observe the flow at very high spatial and temporal resolution. The balanced long-wavelength modes appear only when the Froude number is supercritical (i.e. $F\,{>}\,F_\mathrm{critical}\,{\equiv}\, \upi^2/2$), and are therefore consistent with generation by a baroclinic instability. The unbalanced short-wavelength modes appear locally in every single baroclinically unstable flow, providing perhaps the first direct experimental evidence that all evolving vortical flows will tend to emit freely propagating inertia–gravity waves. The short-wavelength modes also appear in certain baroclinically stable flows. We infer the generation mechanisms of the short-scale waves, both for the baro-clinically unstable case in which they co-exist with a large-scale wave, and for the baroclinically stable case in which they exist alone. The two possible mechanisms considered are spontaneous adjustment of the large-scale flow, and Kelvin–Helmholtz shear instability. Short modes in the baroclinically stable regime are generated only when the Richardson number is subcritical (i.e. $\hbox{\it Ri}\,{<}\,\hbox{\it Ri}_\mathrm{critical}\,{\equiv}\, 1$), and are therefore consistent with generation by a Kelvin–Helmholtz instability. We calculate five indicators of short-wave generation in the baroclinically unstable regime, using data from a quasi-geostrophic numerical model of the annulus. There is excellent agreement between the spatial locations of short-wave emission observed in the laboratory, and regions in which the model Lighthill/Ford inertia–gravity wave source term is large. We infer that the short waves in the baroclinically unstable fluid are freely propagating inertia–gravity waves generated by spontaneous adjustment of the large-scale flow.

An analytical model of mountain wave drag for wind profiles with shear and curvature

An analytical model is developed to predict the surface drag exerted by internal gravity waves on an isolated axisymmetric mountain over which there is a stratified flow with a velocity profile that varies relatively slowly with height. The model is linear with respect to the perturbations induced by the mountain, and solves the Taylor–Goldstein equation with variable coefficients using a Wentzel–Kramers–Brillouin (WKB) approximation, formally valid for high Richardson numbers, Ri. The WKB solution is extended to a higher order than in previous studies, enabling a rigorous treatment of the effects of shear and curvature of the wind profile on the surface drag. In the hydrostatic approximation, closed formulas for the drag are derived for generic wind profiles, where the relative magnitude of the corrections to the leading-order drag (valid for a constant wind profile) does not depend on the detailed shape of the orography. The drag is found to vary proportionally to Ri21, decreasing as Ri decreases for a wind that varies linearly with height, and increasing as Ri decreases for a wind that rotates with height maintaining its magnitude. In these two cases the surface drag is predicted to be aligned with the surface wind. When one of the wind components varies linearly with height and the other is constant, the surface drag is misaligned with the surface wind, especially for relatively small Ri. All these results are shown to be in fairly good agreement with numerical simulations of mesoscale nonhydrostatic models, for high and even moderate values of Ri.