Wind tunnel analysis of the flow topology over complex topographies

El principal objetivo de este trabajo es aportar conocimiento para contestar la pregunta: ¿hasta que punto los ensayos en túnel aerodinámico pueden contribuir a determinar las características que afectan la respuesta dinámica de los aerogeneradores operando en terreno complejo?. Esta pregunta no es nueva, de hecho, el debate en la comunidad científica comenzó en el primer tercio del siglo pasado y aún está intensamente vivo. El método generalmente aceptado para enfrentar el mencionado problema consiste en analizar un caso de estudio determinado en el cual se aplican tanto ensayos a escala real como análisis computacionales y ensayos en túnel aerodinámico. Esto no es ni fácil ni barato. Esta es la razón por la cual desde el experimento de Askervein en 1988, los modelizadores del flujo atmosférico tuvieron que esperar hasta 2007 a que el experimento de Bolund fuese puesto en marcha con un despliegue de medios técnicos equivalentes (teniendo en cuenta la evolución de las tecnologías de sensores y computación). El problema contempla tantos aspectos que ambas experiencias fueron restringidas a condiciones de atmósfera neutra con efectos de Coriolis despreciables con objeto de reducir la complejidad. Este es el contexto en el que se ha desarrollado la presente tesis doctoral. La topología del flujo sobre la isla de Bolund ha sido estudiada mediante la reproducción del experimento de Bolund en los túneles aerodinámicos A9 y ACLA16 del IDR. Dos modelos de la isla de Bolund fueron fabricados a dos escalas, 1:230 y 1:115. El flujo de entrada en el túnel aerodinámico simulando la capa límite sin perturbar correspondía a régimen de transición (transitionally rough regime) y fue usado como situación de referencia. El modelo a escala 1:230 fue ensayado en el túnel A9 para determinar la presión sobre su superficie. La distribución del coeficiente de presión sobre la isla proporcionó una visualización y estimación de una región de desprendimiento sobre el pequeño acantilado situado al frente de la misma. Las medidas de presión instantánea con suficiente grado de resolución temporal pusieron de manifiesto la no estacionariedad en la región de desprendimiento. El modelo a escala 1:115 fue ensayado utilizando hilo caliente de tres componentes y un sistema de velocimetría por imágenes de partículas de dos componentes. El flujo fue caracterizado por el ratio de aceleración, el incremento normalizado de energía cinética turbulenta y los ángulos de inclinación y desviación horizontal. Los resultados a lo largo de la dirección 270°y alturas de 2 m y 5 m presentaron una gran similitud con los resultados a escala real del experimento de Bolund. Los perfiles verticales en las localizaciones de las torres meteorológicas mostraron un acuerdo significativo con los resultados a escala real. El análisis de los esfuerzos de Reynolds y el análisis espectral en las localizaciones de los mástiles meteorológicos presentaron niveles de acuerdo variados en ciertas posiciones, mientras que en otras presentaron claras diferencias. El mapeo horizontal del flujo, para una dirección de viento de 270°, permitió caracterizar el comportamiento de la burbuja intermitente de recirculación sobre el pequeño acantilado existente al frente de la isla así como de la región de relajación y de la capa de cortadura en la región corriente abajo de Bolund. Se realizaron medidas de velocidad con alta resolución espacial en planos perpendiculares a la dirección del flujo sin perturbar. Estas medidas permitieron detectar y caracterizar una estructura de flujo similar a un torbellino longitudinal con regiones con altos gradientes de velocidad y alta intensidad de turbulencia. Esta estructura de flujo es, sin duda, un reto para los modelos computacionales y puede considerarse un factor de riesgo para la operación de los aerogeneradores. Se obtuvieron y analizaron distribuciones espaciales de los esfuerzos de Reynolds mediante 3CHW y PIV. Este tipo de parámetros no constituyen parte de los resultados habituales en los ensayos en túnel sobre topografías y son muy útiles para los modelizadores que utilizan simulación de grades torbellinos (LES). Se proporciona una interpretación de los resultados obtenidos en el túnel aerodinámico en términos de utilidad para los diseñadores de parques eólicos. La evolución y variación de los parámetros del flujo a lo largo de líneas, planos y superficies han permitido identificar como estas propiedades del flujo podrían afectar la localización de los aerogeneradores y a la clasificación de emplazamientos. Los resultados presentados sugieren, bajo ciertas condiciones, la robustez de los ensayos en túnel para estudiar la topología sobre terreno complejo y su comparabilidad con otras técnicas de simulación, especialmente considerando el nivel de acuerdo del conjunto de resultados presentados con los resultados a escala real. De forma adicional, algunos de los parámetros del flujo obtenidos de las medidas en túnel son difícilmente determinables en ensayos a escala real o por medios computacionales, considerado el estado del arte. Este trabajo fue realizado como parte de las actividades subvencionadas por la Comisión Europea como dentro del proyecto FP7-PEOPLE-ITN-2008WAUDIT (Wind Resource Assessment Audit and Standardization) dentro de la FP7 Marie-Curie Initial Training Network y por el Ministerio Español de Economía y Competitividad dentro del proyecto ENE2012-36473, TURCO (Determinación en túnel aerodinámico de la distribución espacial de parámetros estadísticos de la turbulencia atmosférica sobre topografías complejas) del Plan Nacional de Investigación (Subprograma de investigación fundamental no orientada 2012). El informe se ha organizado en siete capítulos y un conjunto de anexos. En el primer capítulo se introduce el problema. En el capítulo dos se describen los medios experimentales utilizados. Seguidamente, en el capítulo tres, se analizan en detalle las condiciones de referencia del principal túnel aerodinámico utilizado en esta investigación. En el capítulo tres se presentan resultados de ensayos de presión superficial sobre un modelo de la isla. Los principales resultados del experimento de Bolund se reproducen en el capítulo cinco. En el capítulo seis se identifican diferentes estructuras del flujo sobre la isla y, finalmente, en el capitulo siete, se recogen las conclusiones y una propuesta de lineas de trabajo futuras. ABSTRACT The main objective of this work is to contribute to answer the question: to which extend can the wind tunnel testing contribute to determine the flow characteristics that affect the dynamic response of wind turbines operating in highly complex terrains?. This question is not new, indeed, the debate in the scientific community was opened in the first third of the past century and it is still intensely alive. The accepted approach to face this problem consists in analysing a given case study where full-scale tests, computational modelling and wind tunnel testing are applied to the same topography. This is neither easy nor cheap. This is is the reason why since the Askervein experience in 1988, the atmospheric flow modellers community had to wait till 2007 when the Bolund experiment was setup with a deployment of technical means equivalent (considering the evolution of the sensor and computing techniques). The problem is so manifold that both experiences were restricted to neutral conditions without Coriolis effects in order to reduce the complexity. This is the framework in which this PhD has been carried out. The flow topology over the Bolund Island has been studied by replicating the Bolund experiment in the IDR A9 and ACLA16 wind tunnels. Two mock-ups of the Bolund island were manufactured at two scales of 1:230 and 1:115. The in-flow in the empty wind tunnel simulating the incoming atmospheric boundary layer was in the transitionally rough regime and used as a reference case. The 1:230 model was tested in the A9 wind tunnel to measure surface pressure. The mapping of the pressure coefficient across the island gave a visualisation and estimation of a detachment region on the top of the escarpment in front of the island. Time resolved instantaneous pressure measurements illustrated the non-steadiness in the detachment region. The 1:115 model was tested using 3C hot-wires(HW) and 2C Particle Image Velocimetry(PIV). Measurements at met masts M3, M6, M7 and M8 and along Line 270°were taken to replicate the result of the Bolund experiment. The flow was characterised by the speed-up ratio, normalised increment of the turbulent kinetic energy, inclination angle and turning angle. Results along line 270°at heights of 2 m and 5 m compared very well with the full-scale results of the Bolund experiment. Vertical profiles at the met masts showed a significant agreement with the full-scale results. The analysis of the Reynolds stresses and the spectral analysis at the met mast locations gave a varied level of agreement at some locations while clear mismatch at others. The horizontal mapping of the flow field, for a 270°wind direction, allowed to characterise the behaviour of the intermittent recirculation bubble on top of the front escarpment followed by a relaxation region and the presence of a shear layer in the lee side of the island. Further detailed velocity measurements were taken at cross-flow planes over the island to study the flow structures on the island. A longitudinal vortex-like structure with high mean velocity gradients and high turbulent kinetic energy was characterised on the escarpment and evolving downstream. This flow structure is a challenge to the numerical models while posing a threat to wind farm designers when siting wind turbines. Spatial distribution of Reynold stresses were presented from 3C HW and PIV measurements. These values are not common results from usual wind tunnel measurements and very useful for modellers using large eddy simulation (LES). An interpretation of the wind tunnel results in terms of usefulness to wind farm designers is given. Evolution and variation of the flow parameters along measurement lines, planes and surfaces indicated how the flow field could affect wind turbine siting. Different flow properties were presented so compare the level of agreement to full-scale results and how this affected when characterising the site wind classes. The results presented suggest, under certain conditions, the robustness of the wind tunnel testing for studying flow topology over complex terrain and its capability to compare to other modelling techniques especially from the level of agreement between the different data sets presented. Additionally, some flow parameters obtained from wind tunnel measurements would have been quite difficult to be measured at full-scale or by computational means considering the state of the art. This work was carried out as a part of the activities supported by the EC as part of the FP7- PEOPLE-ITN-2008 WAUDIT project (Wind Resource Assessment Audit and Standardization) within the FP7 Marie-Curie Initial Training Network and by the Spanish Ministerio de Economía y Competitividad, within the framework of the ENE2012-36473, TURCO project (Determination of the Spatial Distribution of Statistic Parameters of Flow Turbulence over Complex Topographies in Wind Tunnel) belonging to the Spanish National Program of Research (Subprograma de investigación fundamental no orientada 2012). The report is organised in seven chapters and a collection of annexes. In chapter one, the problem is introduced. In chapter two the experimental setup is described. Following, in chapter three, the inflow conditions of the main wind tunnel used in this piece of research are analysed in detail. In chapter three, preliminary pressure tests results on a model of the island are presented. The main results from the Bolund experiment are replicated in chapter five. In chapter six, an identification of specific flow strutures over the island is presented and, finally, in chapter seven, conclusions and lines for future works related to the presented one are included.

Heat flow and heat conductivity coefficient measured in bottom sediments collected in Cruise 13 of R/V Akademik Sergey Vavilov in the Pechora Sea

In Cruise 13 of R/V Akademik Sergey Vavilov in the Pechora Sea, six heat flow varied from 50 to 75 mW/m**2. Deep heat flow in the Pechora Sea was calculated equal to 45 mW/m**2, which is confirmed by results of geological and geophysical studies and corresponds to Middle Baikal age of the basement. A model of structure of the lithosphere in the Pechora Sea is suggested. Total thickness of the lithosphere in the basin (190 km) determined from geothermal data agrees well with that in transition zones from the continent to the ocean. According to estimates of deep heat flow in the region obtained, thickness of the mantle (160 km), of the basaltic (15 km), and of the granitic (15 km) layers of the lithosphere were also evaluated. Temperature values at boundaries of the sedimentary layers were calculated over a geological and geophysical profile crossing the Pechora Sea basin. Temperatures obtained agree with the temperature interval of hydrocarbon generation and correspond to Permian-Triassic sedimentary sequences, which are the most productive ones in the Pechora Sea region from the point of view of oil and gas potential.

Effects of corrugation frequency and aspect ratio on natural convection within an enclosure having sinusoidal corrugation over a heated top surface

Natural convection of a two-dimensional laminar steady-state incompressible fluid flow in a modified rectangular enclosure with sinusoidal corrugated top surface has been investigated numerically. The present study has been carried out for different corrugation frequencies on the top surface as well as aspect ratios of the enclosure in order to observe the change in hydrodynamic and thermal behavior with constant corrugation amplitude. A constant flux heat source is flush mounted on the top sinusoidal wall, modeling a wavy sheet shaded room exposed to sunlight. The flat bottom surface is considered as adiabatic, while the both vertical side walls are maintained at the constant ambient temperature. The fluid considered inside the enclosure is air having Prandtl number of 0.71. The numerical scheme is based on the finite element method adapted to triangular non-uniform mesh element by a non-linear parametric solution algorithm. The results in terms of isotherms, streamlines and average Nusselt numbers are obtained for the Rayleigh number ranging from 10^3 to 10^6 with constant physical properties for the fluid medium considered. It is found that the convective phenomena are greatly influenced by the presence of the corrugation and variation of aspect ratios.

Conduction-radiation effect on natural convection flow in a fluid-saturated non-Darcy porous medium enclosed by non-isothermal walls

The effect of conduction-convection-radiation on natural convection flow of Newtonian optically thick gray fluid, confined in a non-Darcian porous media square cavity is numerically studied. For the gray fluid consideration is given to Rosseland diffusion approximation. Further assuming that (i) the temperature of the left vertical wall is varying linearly with height, (ii) cooled right vertical and top walls and (iii) the bottom wall is uniformly-heated. The governing equations are solved using the Alternate Direct Implicit method together with the Successive Over Relaxation technique. The investigation of the effect of governing parameters namely the Forschheimer resistance (Γ), the Planck constant (Rd), and the temperature difference (Δ), on flow pattern and heat transfer characteristics has been carried out. It was seen that the reduction of flow and heat transfer occurs as the Forschheimer resistance is increased. On the other hand both the strength of flow and heat transfer increases as the temperature ratio, Δ, is increased.

Similarity solutions for flow and heat transfer of non-Newtonian fluid over a stretching surface

Similarity solutions are carried out for flow of power law non-Newtonian fluid film on unsteady stretching surface subjected to constant heat flux. Free convection heat transfer induces thermal boundary layer within a semi-infinite layer of Boussinesq fluid. The nonlinear coupled partial differential equations (PDE) governing the flow and the boundary conditions are converted to a system of ordinary differential equations (ODE) using two-parameter groups. This technique reduces the number of independent variables by two, and finally the obtained ordinary differential equations are solved numerically for the temperature and velocity using the shooting method. The thermal and velocity boundary layers are studied by the means of Prandtl number and non-Newtonian power index plotted in curves.

Time dependent rotational flow of a viscous fluid over an infinite porous disk with a magnetic field

Both the semi-similar and self-similar flows due to a viscous fluid rotating with time dependent angular velocity over a porous disk of large radius at rest with or without a magnetic field are investigated. For the self-similar case the resulting equations for the suction and no mass transfer cases are solved numerically by quasilinearization method whereas for the semi-similar case and injection in the self-similar case an implicit finite difference method with Newton's linearization is employed. For rapid deceleration of fluid and for moderate suction in the case of self-similar flow there exists a layer of fluid, close to the disk surface where the sense of rotation is opposite to that of the fluid rotating far away. The velocity profiles in the absence of magnetic field are found to be oscillatory except for suction. For the accelerating freestream, (semi-similar flow) the effect of time is to reduce the amplitude of the oscillations of the velocity components. On the other hand the effect of time for the oscillating case is just the opposite.

An Analysis of a Boundary Layer Diffusion Flame Over a Porous Flat Plate in a Confined Flow

A numerical analysis of the gas dynamic structure of a two-dimensional laminar boundary layer diffusion flame over a porous flat plate in a confined flow is made on the basis of the familiar boundary layer and flame sheet approximations neglecting buoyancy effects. The governing equations of aerothermochemistry with the appropriate boundary conditions are solved using the Patankar-Spalding method. The analysis predicts the flame shape, profiles of temperature, concentrations of variousspecies, and the density of the mixture across the boundary layer. In addition, it also predicts the pressure gradient in the flow direction arising from the confinement ofthe flow and the consequent velocity overshoot near the flame surface. The results of thecomputation performed for an n-pentane-air system are compared with experimental data andthe agreement is found to be satisfactory.

Three-Dimensional Unsteady Flow With Heat and Mass Transfer Over a Continuous Stretching Surface

A numerical solution of the unsteady boundary layer equations under similarity assumptions is obtained. The solution represents the three-dimensional unsteady fluid motion caused by the time-dependent stretching of a flat boundary. It has been shown that a self-similar solution exists when either the rate of stretching is decreasing with time or it is constant. Three different numerical techniques are applied and a comparison is made among them as well as with earlier results. Analysis is made for various situations like deceleration in stretching of the boundary, mass transfer at the surface, saddle and nodal point flows, and the effect of a magnetic field. Both the constant temperature and constant heat flux conditions at the wall have been studied.

Flow And Heat-Transfer Over An Upstream Moving Wall With A Magnetic-Field And A Parallel Free Stream

The flow and heat transfer over an upstream moving non-isothermal wall with a parallel free stream have been considered. The magnetic field has been applied in the free stream parallel to the wall and the effect of induced magnetic field has been included in the analysis. The boundary layer equations governing the steady incompressible electrically conducting fluid flow have been solved numerically using a shooting method. This problem is interesting because a solution exists only when the ratio of the wall velocity does not exceed a certain critical value and this critical value depends on the magnetic field and magnetic Prandtl number. Also dual solutions exist for a certain range of wall velocity.

A study of over-reflection in Hartmann flow

It is observed that Hartmann flow sustains wave propagation in its centre region for waves whose phase speed is less than the maximum flow speed. Similar to the previous observations it is found that viscous boundary layers around the critical level and at the wall replace the exponential regions and wave sinks required for over-reflection in the inviscid flow. The uniform magnetic field stabilizes the flow for small-wave-number disturbances along thez-direction. Over-reflection is confined to a few ranges of phase speeds for which the two boundary layers are close together rather than widely separated. These ranges correspond exactly to those for which unstable eigenmodes exist. Over-reflection is associated with a wave phase tilt opposite in direction to the shear.

Mixed convection heat transfer from the bottom tip of a cylinder spinning about a vertical axis in a saturated porous medium

The present work is a numerical study of heat transfer characteristics from the bottom tip of a cylinder spinning about a vertical axis in an infinitely saturated porous medium. The problem is axisymmetric. The non-dimensionalized governing equations are solved using the SIMPLER algorithm on a staggered grid. The influence of rotational Reynolds numbers and Darcy numbers on the heat transfer for a Grashof number of 104 and Prandtl number of 7.0 is studied. It is found that for very high Darcy numbers, over a wide range of rotational Reynolds numbers, the heat transfer takes place mainly due to conduction. The convective heat transfer takes place for lower Darcy numbers and for higher rotational Reynolds numbers. Moreover, there is a rapid increase in the overall Nusselt number below a certain Darcy number with increase in the rotational Reynolds numbers. The effect of the Darcy number and the rotational Reynolds number on the heat transfer and fluid flow in the porous medium is depicted in the form of streamline and isotherm plots. The variation of the overall Nusselt number with respect to the Darcy number for various rotational Reynolds numbers is plotted. The variation of the local Nusselt number with respect to the radial coordinate at the heated tip of the vertical cylinder is plotted for various Darcy and rotational Reynolds numbers.

Flow induced by an asymmetrically placed disk rotating coaxially inside a cylindrical casing

In this paper, the flow due to a rotating disk non-symmetrically placed with respect to the height of the enclosing stationary cylinder is analyzed numerically. The full Navier-Stokes equations expressed in terms of stream function and vorticity are solved by successive over-relaxation for different disk radii, its distance from the bottom casing and rotational Reynolds numbers. It is observed that the flow pattern is strongly influenced by the size and the position of the disk. When the disk is very close to the top casing and small in radius, there are two regions of different scales and the vortices in the region of small scale are trapped between the disk and the top casing. Further, the variation of the moment coefficient is determined for different positions and sizes of the rotating disk. The calculations shows that the frictional torque increases rapidly, when the disk approaches the top casing. This finding is of importance for the design of vertical rotating disk reactors applied in chemical vapor deposition.

Levitation of a drop over a film flow

A vertical jet of water impinging on a horizontal surface produces a radial film flow followed by a circular hydraulic jump. We report a phenomenon where fairly large (1 mi) drops of liquid levitate just upstream of the jump on a thin air layer between the drop and the film flow. We explain the phenomenon using lubrication theory. Bearing action both in the air film and the water film seems to be necessary to support large drops. Horizontal support is given to the drop by the hydraulic jump. A variety of drop shapes is observed depending on the volume of the drop and liquid properties. We show that interaction of the forces due to gravity, surface tension, viscosity and inertia produces these various shapes.

Computational and experimental studies of fluid flow and heat transfer in a Calandria based reactor

CFD investigations are carried out to study the heat flux and temperature distribution in the calandria using a 3–Dimensional RANS code. Internal flow computations and experimental studies are carried out for a calandria embedded with a matrix of tubes working together as a reactor. Numerical investigations are carried on the Calandria reactor vessel with horizontal inlets and outlets located on top and the bottom to study the flow pattern and the associated temperature distribution. The computations have been carried out to simulate fluid flow and convective heat transfer for assigned near–to working conditions with different moderator injection rates and reacting heat fluxes. The results of computations provide an estimate of the tolerance bands for safe working limits for the heat dissipation at different working conditions by virtue of prediction of the hot spots in the calandria. The isothermal CFD results are validated by a set of experiments on a specially designed scaled model conducted over a range of flows and simulation parameters. The comparison of CFD results with experiments show good agreement.

Surface pressure and heat transfer measurements in the unsteady separated hypersonic flow field over double cones

Reliable bench mark experimental database in the separated hypersonic flow regime is necessary to validate high resolution CFD codes. In this paper we report the surface pressure and heat transfer measurements carried out on double cones (first cone semi-apex angle = 15, 25 deg.; second cone semi-apex angle= 35, 68 deg.) at hypersonic speeds that will be useful for CFD code validation studies. The surface pressure measurements are carried out at nominal Mach number of 8.35 in the IISc hypersonic wind tunnel. On the other hand the surface heat transfer measurements are carried out at a nominal Mach number of 5.75 in the IISc hypersonic shock tunnel. The flow separation point on the first cone, flow reattachment on the second cone and the wild fluctuation of the transmitted shock on the second cone surface (25/68 deg. double cone) in the presence of severe adverse pressure gradient are some of the flow features captured in the measurements. The results from the CFD studies indicate good agreement with experiments in the attached flow regime while considerable differences are noticeable in the separated flow regime.