An investigation of wind tunnel wall effects at high Mach numbers

This paper presents the results of an investigation of wind tunnel wall interference in a two-dimensional wind tunnel at high Mach numbers. The results are presented in the form of curves of lift coefficient versus the ratio of model chord to tunnel height, as functions of Mach number and angle of attack. The investigation was carried out by the authors at the Guggenheim Aeronautical Laboratory of the California Institute of Technology during the school year 1944-45.

Tests were carried out on the NACA low drag airfoil section 65,1-012 at Mach numbers from .60 to .80, and angles of attack of from 1 to 3 degrees. Models were 1", 2", 4" and 6" chord, giving values of the chord to tunnel height ration of .1 to .6. Schlieren photographs were made of shock waves where they occurred.

Study of the turbulence in the air-side and water-side boundary layers in experimental laboratory wind induced surface waves

This study detailed the structure of turbulence in the air-side and water-side boundary layers in wind-induced surface waves. Inside the air boundary layer, the kurtosis is always greater than 3 (the value for normal distribution) for both horizontal and vertical velocity fluctuations. The skewness for the horizontal velocity is negative, but the skewness for the vertical velocity is always positive. On the water side, the kurtosis is always greater than 3, and the skewness is slightly negative for the horizontal velocity and slightly positive for the vertical velocity. The statistics of the angle between the instantaneous vertical fluctuation and the instantaneous horizontal velocity in the air is similar to those obtained over solid walls. Measurements in water show a large variance, and the peak is biased towards negative angles. In the quadrant analysis, the contribution of quadrants Q2 and Q4 is dominant on both the air side and the water side. The non-dimensional relative contributions and the concentration match fairly well near the interface. Sweeps in the air side (belonging to quadrant Q4) act directly on the interface and exert pressure fluctuations, which, in addition to the tangential stress and form drag, lead to the growth of the waves. The water drops detached from the crest and accelerated by the wind can play a major role in transferring momentum and in enhancing the turbulence level in the water side.On the air side, the Reynolds stress tensor's principal axes are not collinear with the strain rate tensor, and show an angle α σ≈=-20°to-25°. On the water side, the angle is α σ≈=-40°to-45°. The ratio between the maximum and the minimum principal stresses is σ a/σ b=3to4 on the air side, and σ a/σ b=1.5to3 on the water side. In this respect, the air-side flow behaves like a classical boundary layer on a solid wall, while the water-side flow resembles a wake. The frequency of bursting on the water side increases significantly along the flow, which can be attributed to micro-breaking effects - expected to be more frequent at larger fetches. © 2012 Elsevier B.V.

Efficiency improvement study for small wind turbines through flow control

In this study, a constant suction technique for controlling boundary layer separation at low Reynolds numbers was designed and tested. This was later implemented on small wind turbines. Small wind turbines need to operate in low wind speeds, that is, in low Reynolds number regimes – typically in the range 104–105. Airfoils are prone to boundary layer separation in these conditions, leading to a substantial drop in aerodynamic performance of the blades. Under these conditions turbines will have reduced energy output. This paper presents experimental results of applying surface-suction over the suction-surface of airfoils for controlling boundary layer separation. The Reynolds numbers for the experiments are kept in the range 8×104–5×105. The air over the surface of the airfoil is drawn into the airfoil through a slit. It is found that the lift coefficient of the airfoils increases and the drag reduces. Based on the improved airfoil characteristics, an analysis of increase in Coefficient of Power (CP), versus input power for a small wind turbine blade with constant suction is presented.

Techniques for Enhancing Wind Energy Generation - A CFD Based Multibody Dynamics Approach in Horizontal Axis Wind Turbines

Wind energy has emerged as a major sustainable source of energy.The efficiency of wind power generation by wind mills has improved a lot during the last three decades.There is still further scope for maximising the conversion of wind energy into mechanical energy.In this context,the wind turbine rotor dynamics has great significance.The present work aims at a comprehensive study of the Horizontal Axis Wind Turbine (HAWT) aerodynamics by numerically solving the fluid dynamic equations with the help of a finite-volume Navier-Stokes CFD solver.As a more general goal,the study aims at providing the capabilities of modern numerical techniques for the complex fluid dynamic problems of HAWT.The main purpose is hence to maximize the physics of power extraction by wind turbines.This research demonstrates the potential of an incompressible Navier-Stokes CFD method for the aerodynamic power performance analysis of horizontal axis wind turbine.The National Renewable Energy Laboratory USA-NREL (Technical Report NREL/Cp-500-28589) had carried out an experimental work aimed at the real time performance prediction of horizontal axis wind turbine.In addition to a comparison between the results reported by NREL made and CFD simulations,comparisons are made for the local flow angle at several stations ahead of the wind turbine blades.The comparison has shown that fairly good predictions can be made for pressure distribution and torque.Subsequently, the wind-field effects on the blade aerodynamics,as well as the blade/tower interaction,were investigated.The selected case corresponded to a 12.5 m/s up-wind HAWT at zero degree of yaw angle and a rotational speed of 25 rpm.The results obtained suggest that the present can cope well with the flows encountered around wind turbines.The areodynamic performance of the turbine and the flow details near and off the turbine blades and tower can be analysed using theses results.The aerodynamic performance of airfoils differs from one another.The performance mainly depends on co-efficient of performnace,co-efficient of lift,co-efficient of drag, velocity of fluid and angle of attack.This study shows that the velocity is not constant for all angles of attack of different airfoils.The performance parameters are calculated analytically and are compared with the standardized performance tests.For different angles of ,the velocity stall is determined for the better performance of a system with respect to velocity.The research addresses the effect of surface roughness factor on the blade surface at various sections.The numerical results were found to be in agreement with the experimental data.A relative advantage of the theoretical aerofoil design method is that it allows many different concepts to be explored economically.Such efforts are generally impractical in wind tunnels because of time and money constraints.Thus, the need for a theoretical aerofoil design method is threefold:first for the design of aerofoil that fall outside the range of applicability of existing calalogs:second,for the design of aerofoil that more exactly match the requirements of the intended application:and third,for the economic exploration of many aerofoil concepts.From the results obtained for the different aerofoils,the velocity is not constant for all angles of attack.The results obtained for the aerofoil mainly depend on angle of attack and velocity.The vortex generator technique was meticulously studies with the formulation of the specification for the right angle shaped vortex generators-VG.The results were validated in accordance with the primary analysis phase.The results were found to be in good agreement with the power curve.The introduction of correct size VGs at appropriate locations over the blades of the selected HAWT was found to increase the power generation by about 4%

Wind-direction effects on urban-type flows

Practically all extant work on flows over obstacle arrays, whether laboratory experiments or numerical modelling, is for cases where the oncoming wind is normal to salient faces of the obstacles. In the field, however, this is rarely the case. Here, simulations of flows at various directions over arrays of cubes representing typical urban canopy regions are presented and discussed. The computations are of both direct numerical simulation and large-eddy simulation type. Attention is concentrated on the differences in the mean flow within the canopy region arising from the different wind directions and the consequent effects on global properties such as the total surface drag, which can change very significantly—by up to a factor of three in some circumstances. It is shown that for a given Reynolds number the typical viscous forces are generally a rather larger fraction of the pressure forces (principally the drag) for non-normal than for normal wind directions and that, dependent on the surface morphology, the average flow direction deep within the canopy can be largely independent of the oncoming wind direction. Even for regular arrays of regular obstacles, a wind direction not normal to the obstacle faces can in general generate a lateral lift force (in the direction normal to the oncoming flow). The results demonstrate this and it is shown how computations in a finite domain with the oncoming flow generated by an appropriate forcing term (e.g. a pressure gradient) then lead inevitably to an oncoming wind direction aloft that is not aligned with the forcing term vector.

Assessing wind profile effects on the global atmospheric torque

The impact of a new approach to the evaluation of surface gravity wave drag (GWD) is assessed. This approach uses linear theory, but incorporates the effects of wind profile shear and curvature, by means of a second-order WKB approximation. While the theory predicts the possibility of either drag enhancement or reduction, depending on the wind profile, results obtained with the ERA-40 reanalysis data clearly indicate the predominance of local drag enhancement. However, the global impact of shear on the atmospheric axial GWD torque comes mostly from regions with predominantly easterly flow, contributing to a slight reduction of the bias found in different studies of the global angular momentum budget. The relative correction due to shear on linear GWD is found not to depend too strongly on the levels chosen for the computation of the low-level wind derivatives.

Mountain waves in two-layer sheared flows: Critical-level effects, wave reflection, and drag enhancement

Internal gravity waves generated in two-layer stratified shear flows over mountains are investigated here using linear theory and numerical simulations. The impact on the gravity wave drag of wind profiles with constant unidirectional or directional shear up to a certain height and zero shear above, with and without critical levels, is evaluated. This kind of wind profile, which is more realistic than the constant shear extending indefinitely assumed in many analytical studies, leads to important modifications in the drag behavior due to wave reflection at the shear discontinuity and wave filtering by critical levels. In inviscid, nonrotating, and hydrostatic conditions, linear theory predicts that the drag behaves asymmetrically for backward and forward shear flows. These differences primarily depend on the fraction of wavenumbers that pass through their critical level before they are reflected by the shear discontinuity. If this fraction is large, the drag variation is not too different from that predicted for an unbounded shear layer, while if it is small the differences are marked, with the drag being enhanced by a considerable factor at low Richardson numbers (Ri). The drag may be further enhanced by nonlinear processes, but its qualitative variation for relatively low Ri is essentially unchanged. However, nonlinear processes seem to interact constructively with shear, so that the drag for a noninfinite but relatively high Ri is considerably larger than the drag without any shear at all.

Asymptotic gravity wave drag expressions for non-hydrostatic rotating flow over a ridge

Asymptotic expressions are derived for the mountain wave drag in flow with constant wind and static stability over a ridge when both rotation and non-hydrostatic effects are important. These expressions, which are much more manageable than the corresponding exact drag expressions (when these do exist) are found to provide accurate approximations to the drag, even when non-hydrostatic and rotation effects are strong, despite having been developed for cases where these effects are weak. The derived expressions are compared with approximations to the drag found previously, and their asymptotic behaviour in various limits is studied.

Internal wave drag in stratified flow over mountains on a beta plane

The impact of the variation of the Coriolis parameter f on the drag exerted by internal Rossby-gravity waves on elliptical mountains is evaluated using linear theory, assuming constant wind and static stability and a beta-plane approximation. Previous calculations of inertia-gravity wave drag are thus extended in an attempt to establish a connection with existing studies on planetary wave drag, developed primarily for fluids topped by a rigid lid. It is found that the internal wave drag for zonal westerly flow strongly increases relative to that given by the calculation where f is assumed to be a constant, particularly at high latitudes and for mountains aligned meridionally. Drag increases with mountain width for sufficiently wide mountains, reaching values much larger than those valid in the non-rotating limit. This occurs because the drag receives contributions from a low wavenumber range, controlled by the beta effect, which accounts for the drag amplification found here. This drag amplification is shown to be considerable for idealized analogues of real mountain ranges, such as the Himalayas and the Rocky mountains, and comparable to the barotropic Rossby wave drag addressed in previous studies.

A linear model of gravity wave drag for hydrostatic sheared flow over elliptical mountains

An analytical model of orographic gravity wave drag due to sheared flow past elliptical mountains is developed. The model extends the domain of applicability of the well-known Phillips model to wind profiles that vary relatively slowly in the vertical, so that they may be treated using a WKB approximation. The model illustrates how linear processes associated with wind profile shear and curvature affect the drag force exerted by the airflow on mountains, and how it is crucial to extend the WKB approximation to second order in the small perturbation parameter for these effects to be taken into account. For the simplest wind profiles, the normalized drag depends only on the Richardson number, Ri, of the flow at the surface and on the aspect ratio, γ, of the mountain. For a linear wind profile, the drag decreases as Ri decreases, and this variation is faster when the wind is across the mountain than when it is along the mountain. For a wind that rotates with height maintaining its magnitude, the drag generally increases as Ri decreases, by an amount depending on γ and on the incidence angle. The results from WKB theory are compared with exact linear results and also with results from a non-hydrostatic nonlinear numerical model, showing in general encouraging agreement, down to values of Ri of order one.

Resonant gravity-wave drag enhancement in linear stratified flow over mountains

High-drag states produced in stratified flow over a 2D ridge and an axisymmetric mountain are investigated using a linear, hydrostatic, analytical model. A wind profile is assumed where the background velocity is constant up to a height z1 and then decreases linearly, and the internal gravity-wave solutions are calculated exactly. In flow over a 2D ridge, the normalized surface drag is given by a closed-form analytical expression, while in flow over an axisymmetric mountain it is given by an expression involving a simple 1D integral. The drag is found to depend on two dimensionless parameters: a dimensionless height formed with z_1, and the Richardson number, Ri, in the shear layer. The drag oscillates as z_1 increases, with a period of half the hydrostatic vertical wavelength of the gravity waves. The amplitude of this modulation increases as Ri decreases. This behaviour is due to wave reflection at z_1. Drag maxima correspond to constructive interference of the upward- and downward-propagating waves in the region z < z_1, while drag minima correspond to destructive interference. The reflection coefficient at the interface z = z_1 increases as Ri decreases. The critical level, z_c, plays no role in the drag amplification. A preliminary numerical treatment of nonlinear effects is presented, where z_c appears to become more relevant, and flow over a 2D ridge qualitatively changes its character. But these effects, and their connection with linear theory, still need to be better understood.

Drag produced by trapped lee waves and propagating mountain waves in a two-layer atmosphere

The surface drag force produced by trapped lee waves and upward propagating waves in non-hydrostatic stratified flow over a mountain ridge is explicitly calculated using linear theory for a two-layer atmosphere with piecewise-constant static stability and wind speed profiles. The behaviour of the drag normalized by its hydrostatic single-layer reference value is investigated as a function of the ratio of the Scorer parameters in the two layers l_2/l_1 and of the corresponding dimensionless interface height l_1 H, for selected values of the dimensionless ridge width l_1 a and ratio of wind speeds in the two layers. When l_2/l_1 → 1, the propagating wave drag approaches 1 in approximately hydrostatic conditions, and the trapped lee wave drag vanishes. As l_2/l_1 decreases, the propagating wave drag progressively displays an oscillatory behaviour with l_1 H, with maxima of increasing magnitude due to constructive interference of reflected waves in the lower layer. The trapped lee wave drag shows localized maxima associated with each resonant trapped lee wave mode, occurring for small l_2/l_1 and slightly higher values of l_1 H than the propagating wave drag maxima. As l1a decreases, i.e. the flow becomes more non-hydrostatic, the propagating wave drag decreases and the regions of non-zero trapped lee wave drag extend to higher l_2/l_1. These results are confirmed by numerical simulations for l_2/l_1 = 0.2. In parameter ranges of meteorological relevance, the trapped lee wave drag may have a magnitude comparable to that of propagating wave drag, and be larger than the reference single-layer drag. This may have implications for drag parametrization in global climate and weather-prediction models.

Constraints on wave drag parameterization schemes for simulating the quasi-biennial oscillation. Part I: gravity wave forcing

Parameterization schemes for the drag due to atmospheric gravity waves are discussed and compared in the context of a simple one-dimensional model of the quasi-biennial oscillation (QBO). A number of fundamental issues are examined in detail, with the goal of providing a better understanding of the mechanism by which gravity wave drag can produce an equatorial zonal wind oscillation. The gravity wave–driven QBOs are compared with those obtained from a parameterization of equatorial planetary waves. In all gravity wave cases, it is seen that the inclusion of vertical diffusion is crucial for the descent of the shear zones and the development of the QBO. An important difference between the schemes for the two types of waves is that in the case of equatorial planetary waves, vertical diffusion is needed only at the lowest levels, while for the gravity wave drag schemes it must be included at all levels. The question of whether there is downward propagation of influence in the simulated QBOs is addressed. In the gravity wave drag schemes, the evolution of the wind at a given level depends on the wind above, as well as on the wind below. This is in contrast to the parameterization for the equatorial planetary waves in which there is downward propagation of phase only. The stability of a zero-wind initial state is examined, and it is determined that a small perturbation to such a state will amplify with time to the extent that a zonal wind oscillation is permitted.

Constraints on wave drag parameterization schemes for simulating the quasi-biennial oscillation. Part II: combined effects of gravity waves and equatorial planetary waves

This study examines the effect of combining equatorial planetary wave drag and gravity wave drag in a one-dimensional zonal mean model of the quasi-biennial oscillation (QBO). Several different combinations of planetary wave and gravity wave drag schemes are considered in the investigations, with the aim being to assess which aspects of the different schemes affect the nature of the modeled QBO. Results show that it is possible to generate a realistic-looking QBO with various combinations of drag from the two types of waves, but there are some constraints on the wave input spectra and amplitudes. For example, if the phase speeds of the gravity waves in the input spectrum are large relative to those of the equatorial planetary waves, critical level absorption of the equatorial planetary waves may occur. The resulting mean-wind oscillation, in that case, is driven almost exclusively by the gravity wave drag, with only a small contribution from the planetary waves at low levels. With an appropriate choice of wave input parameters, it is possible to obtain a QBO with a realistic period and to which both types of waves contribute. This is the regime in which the terrestrial QBO appears to reside. There may also be constraints on the initial strength of the wind shear, and these are similar to the constraints that apply when gravity wave drag is used without any planetary wave drag. In recent years, it has been observed that, in order to simulate the QBO accurately, general circulation models require parameterized gravity wave drag, in addition to the drag from resolved planetary-scale waves, and that even if the planetary wave amplitudes are incorrect, the gravity wave drag can be adjusted to compensate. This study provides a basis for knowing that such a compensation is possible.

The physics of orographic gravity wave drag

The drag and momentum fluxes produced by gravity waves generated in flow over orography are reviewed, focusing on adiabatic conditions without phase transitions or radiation effects, and steady mean incoming flow. The orographic gravity wave drag is first introduced in its simplest possible form, for inviscid, linearized, non-rotating flow with the Boussinesq and hydrostatic approximations, and constant wind and static stability. Subsequently, the contributions made by previous authors (primarily using theory and numerical simulations) to elucidate how the drag is affected by additional physical processes are surveyed. These include the effect of orography anisotropy, vertical wind shear, total and partial critical levels, vertical wave reflection and resonance, non-hydrostatic effects and trapped lee waves, rotation and nonlinearity. Frictional and boundary layer effects are also briefly mentioned. A better understanding of all of these aspects is important for guiding the improvement of drag parametrization schemes.