<p>Highly loaded airfoils in low-pressure turbines (LPTs) suffer from laminar flow separation from the suction side of the airfoils aft of the throat of the passages. This separation harms the performance of the engine by reducing the power extraction from the turning air and ultimately reduces the overall turbine efficiency. Flow control techniques have been investigated to eliminate flow separation in aerodynamic surfaces to abate the losses associated with it. This Master of Science Thesis investigates the design, implementation and testing of pulsated injection actuation in a low-Reynolds flow over a wall-mounted hump.</p> <p>Furthermore, this Thesis expands on the existing expertise in the infrared (IR) thermography measurement technique at the Purdue Experimental Turbine Aerothermal Lab. This is done through an investigation of the factors affecting the IR measurement technique and the development of an optical instrument (borescope) to implement in an annular cascade wind tunnel. IR thermography is used on the wall-mounted hump blowdown tests to detect the separation point in the boundary layer using two techniques: by an investigation of the surface temperature distribution and an investigation of the heat transfer behavior at the surface. Finally, the borescope is commissioned through the first testing campaign of the LPT airfoils, and are processed to thermally investigate the passage.</p> <p>This thesis succeeds in expanding the IR capabilities within PETAL, and at demonstrating pulsated injection as an effective method to eliminate flow separation. Furthermore, IR successfully detects flow separation on the wall-mounted hump through the two methods presented, as well as detecting the boundary layer reattachment caused by the flow control technique. The limitations of the thermal methodology, as well as those of the optical probe are addressed, and the uncertainties in the measurements are quantified. Finally, steps to continue the studies are suggested at the end of each methodology chapter, including the potential redesign of the IR borescope to improve the quality of measurements. </p>

Many flow cases in fluid dynamics face undesirable flow separation due to ad-verse pressure gradients on wall boundaries. This occurs, for example, due togeometrical reasons as in a highly curved turbine-inlet duct or on flow-controlsurfaces such as wing trailing-edge flaps within a certain angle-of-attack range.Here, flow-control devices are often used in order to enhance the flow and delayor even totally eliminate flow separation. Flow control can e.g. be achieved byusing passive or active vortex generators (VGs) for momentum mixing in theboundary layer of such flows. This thesis focusses on such passive and activeVGs and their modelling for computational fluid dynamics investigations. First, a statistical VG model approach for passive vane vortex genera-tors (VVGs), developed at the Royal Institute of Technology Stockholm andthe Swedish Defence Research Agency, was evaluated and further improvedby means of experimental data and three-dimensional fully-resolved computa-tions. This statistical VVG model approach models those statistical vortexstresses that are generated at the VG by the detaching streamwise vortices.This is established by means of the Lamb-Oseen vortex model and the Prandtllifting-line theory for the determination of the vortex strength. Moreover, thisansatz adds the additional vortex stresses to the turbulence of a Reynolds-stresstransport model. Therefore, it removes the need to build fully-resolved three-dimensional geometries of VVGs in a computational fluid dynamics mesh. Usu-ally, the generation of these fully-resolved geometries is rather costly in termsof preprocessing and computations. By applying VVG models, the costs arereduced to that of computations without VVGs. The original and an improvedcalibrated passive VVG model show sensitivity for parameter variations suchas the modelled VVG geometry and the VVG model location on a flat plate inzero- and adverse-pressure-gradient flows, in a diffuser, and on an airfoil withits high-lift system extracted. It could be shown that the passive VG modelqualitatively and partly quantitatively describes correct trends and tendenciesfor these different applications. In a second step, active vortex-generator jets (VGJs) are considered. They were experimentally investigated in a zero-pressure-gradient flat-plate flow atTechnische Universitä̈t Braunschweig, Germany, and have been re-evaluated for our purposes and a parameterization of the generated vortices was conducted. Dependencies of the generated vortices and their characteristics on the VGJsetup parameters could be identified and quantified. These dependencies wereused as a basis for the development of a new statistical VGJ model. This modeluses the ansatz of the passive VVG model in terms of the vortex model, theadditional vortex-stress tensor, and its summation to the Reynolds stress ten-sor. Yet, it does not use the Prandtl lifting-line theory for the determinationof the circulation but an ansatz for the balance of the momentum impact thatthe VGJ has on the mean flow. This model is currently under developmentand first results have been evaluated against experimental and fully-resolvedcomputational results of a flat plate without pressure gradient. / <p>QC 20120511</p>

This study presents a numerical investigation of ducted tidal turbines, employing three-dimensional Reynolds-averaged Navier-Stokes simulations. Bidirectional ducted turbines are modelled with and without aperture, referred to as ducted and open- centre turbines respectively. The work consists of two investigations. In the first, the turbine rotors are represented by actuator discs, a simplification which captures changes in linear momentum and thus the primary interaction of the turbine with the flow through and around the duct, while greatly reducing computational complexity. In the second investigation, the turbine rotors are represented through a CFD-integrated blade element momentum model, employing realistic rotor data, capturing swirl and blade drag in addition to the extraction of linear momentum. Both modelling techniques were employed to investigate the performances of bare, ducted, and open-centre turbines, relating these to the flow fields exhibited. For axial flow, substantial decreases in power generated by the ducted and open-centre turbines were found, relative to a bare turbine of equal total device diameter. For open-centre turbines, an increase in aperture size leads to a further reduction in power generated. Increased blockage was shown to positively affect the performance of all devices. Two further measures of performance were employed: power density, normalising the power by the rotor area, and basin efficiency, relating the power generated to the overall power removed from the flow. Moderate increases in power density can be achieved for the ducted and open-centre devices, while their basin efficiencies are of similar value to that of the bare turbine. For yawed inflow, the performance of the bare turbine decreases, whilst that of the ducted and open-centre turbines increases. This is due to an increased flow velocity following flow acceleration around the inlet lip of the duct and also an increase in effective blockage as ducts present greater projected frontal area when approached non-axially.

No description available.

Global renewable energy requirements rapidly increase with the transition to a fossil-free society. As a result, intermittent energy resources, such as wind- and solar power, have become increasingly popular. However, their energy production varies over time, both in the short- and long term. Hydropower plants are therefore utilized as a regulating resource more frequently to maintain a balance between production and consumption on the electrical grid. This means that they must be operated away from the design point, also known as the best-efficiency-point (BEP), and often are operated at part-load (PL) with a lower power output. Moreover, some plants are expected to provide a spinning reserve, also referred to as speed-no-load (SNL), to respond rapidly to power shortages. During this operating condition, the turbine rotates without producing any power. During the above mentioned off-design operating conditions, the flow rate is restricted by the closure of the guide vanes. This changes the absolute velocity of the flow and increases the swirl, which is unfavorable. The flow field can be described as chaotic, with separated regions and recirculating fluid. Shear layer formation between stagnant- and rotating flow regions can be an origin for rotating flow structures. Examples are the rotating-vortex-rope (RVR) found during PL operation and the vortical flow structures in the vaneless space during SNL operation, which can cause the flow between the runner blades to stall, also referred to as rotating stall. The flow structures are associated with pressure pulsations throughout the turbine, which puts high stress on the runner and other critical parts and shortens the turbine's lifetime. Numerical models of hydraulic turbines are highly coveted to investigate the detrimental flow inside the hydraulic turbines' different sections at off-design operating conditions. They enable the detailed study of the flow and the origin of the instabilities. This knowledge eases the design and assessment of mitigation techniques that expand the turbines' operating range, ultimately enabling a wider implementation of intermittent energy resources on the electrical grid and a smoother transition to a fossil-free society. This thesis presents the numerical study of the Porjus U9 model, a scaled-down version of the 10 MW prototype Kaplan turbine located along the Luleå river in northern Sweden. The distributor contains 20 guide vanes, 18 stay vanes and the runner is 6-bladed. The numerical model is a geometrical representation of the model turbine located at Vattenfall Research and Development in Älvkarleby, Sweden. The commercial software ANSYS CFX 2020 R2 is used to perform the numerical simulations. Firstly, the draft tube cone section of the U9 model is numerically studied to investigate the sensitivity of a swirling flow to the GEKO (generalized kω) turbulence model. The GEKO model aims to consolidate different eddy viscosity turbulence models. Six free coefficients are changeable to tune the model to flow conditions and obtain results closer to an experimental reference without affecting the calibration of the turbulence model to basic flow test cases. Especially, the coefficients affecting wall-bounded flows are of interest. This study aims to analyze if the GEKO model can be used to obtain results closer to experimental measurements and better predict the swirling flow at PL operation compared to other eddy viscosity turbulence models. Results show that the near-wall- and separation coefficients predict a higher swirl and give results closer to experimentally obtained ones. Secondly, a simplified version of the U9 model is investigated at BEP and PL operating conditions and includes one distributor passage with periodic boundary conditions, the runner and the draft tube. The flow is assumed axisymmetric upstream of the runner, hence the single distributor passage. Previous studies of hydraulic turbines operating at PL show difficulties predicting the flow's tangential velocity component as it is often under predicted. Therefore, a parametric analysis is performed to investigate which parameters affect the prediction of the tangential velocity in the runner domain. Results show that the model predicts the flow relatively well at BEP but has problems at PL; the axial velocity is overpredicted while the tangential is underpredicted. Moreover, the torque is overpredicted. The root cause for the deviation is an underestimation of the head losses. Another contributing reason is that the runner extracts too much swirl from the flow, hence the low tangential velocity and the high torque. Sensitive parameters are the blade clearance, blade angle and mass flow. Finally, the full version of the U9 model is analyzed at SNL operation, including the spiral casing, full distributor, runner and draft tube. During this operating condition, the flow is not axisymmetric; vortical flow structures extend from the vaneless space to the draft tube and the flow stalls between the runner blades. A mitigation technique with independent control of each guide vane is presented and implemented in the model. The idea is to open some of the guidevanes to BEP angle while keeping the remaining ones closed. The aim is to reduce the swirl and prevent the vortical flow structures from developing. Results show that the flow structures are broken down upstream the runner and the rotating stall between the runner blades is reduced, which decreases the pressure- and velocity fluctuations. The flow down stream the runner remains mainly unchanged.

Flow separation and reattachment happens in many fluid mechanical situations occurring in engineering applications as well as in nature. The flow over a backward-facing step represents a geometrically simple flow situation exhibiting both flow separation and reattachment. Broadly speaking there are only two important parameters in the problem, the Reynolds number(Re) based on the step height(h),and a geometrical parameter, referred to as the Expansion ratio(ER), defined as the downstream channel height to the upstream channel height. In spite of the relative simplicity of this geometry, the flow downstream is quite complex. The main focus of the present work is to elucidate the unsteady three-dimensional coherent structures present in this flow at large Re, Re>36,000,based on the step height(h). For this, we use velocity field measurements from Particle Image Velocimetry (PIV)in conjunction with hotwire anemometry measurements. The time-averaged structure of this flow is first studied in detail, including the effect of Reynolds number(Re) and Expansion Ratio(ER), on it. These studies show that at sufficiently large Re (Re>20,000), the reattachment length becomes independent of Re. The detailed internal structure of the separation bubble is also found to be independent of Re, but for Revalues that are relatively larger(Re>36,000). At large Re, the main effect of ER ,is found to be on the reattachment length, which increases with ER and saturates for ER values greater than about 1.8. The detailed internal structure of the separation bubble has been mapped at high Re and is found to be nearly the same for all ER, when the streamwise length is normalized by the reattachment length. In order to elucidate the unsteady coherent vortical structures, PIV measurements are done in two orthogonal planes downstream of the backward-facing step. These measurements are done for ER= 1.50 at large Re(Re=36,000) and in a large aspect ratio facility(AR= span length/step height= 24); the latter being important to avoid any effects due to span-wise confinement. In the spanwise plane parallel to the lower wall(x-z plane),instantaneous velocity fields show counter rotating vortex pairs, which is a signature of the three-dimensional vortical structures in this plane. Using conditional averaging, this counter-rotating vortex pair signature is captured right from upstream of the step, to well after reattachment. Spatial correlations are used to get the length scale of these coherent vortical structures, which varies substantially from the attached boundary layer before separation to the region after reattachment. The variation of these structures in the cross-stream (vertical) direction at reattachment and beyond gives an idea about their three dimensional shape. The circulation of these counter-rotating pairs is measured from the conditionally aver-aged fields, and is found to increase with streamwise distance reaching normalized circulation values (Γ/Uoh) of about 0.5 around reattachment. Velocity spectra downstream of the step show peaks corresponding to both the shear layer frequency(Stsl)and a relatively lower frequency that corresponds to large-scale shedding from the separation bubble (Stb); the latter in particular being quasi-periodic. Small amplitude sinusoidal forcing at the shedding frequency(Stb) is applied close to the step, by blowing and suction, to make the quasi-periodic shedding more regular. Measurements show that this has a very small effect on both the mean separation bubble and on the counter-rotating structures in the x-z plane. This mild forcing however enables phase locked PIV measurements to be made which shows the bubble shedding phenomenon in the cross-stream plane(side view or x-y plane). The phase-averaged velocity fields show significant variations from phase to phase. Although there is some hint of structures being shed, from these phase-averaged fields, it is not very clear. One of the primary reasons is the fact that the flow is effectively spanwise averaged, as the three-dimensional structures are not locked in the spanwise direction. To get a three dimensional view of the sheddin gphenomenon, it is necessary to lock the spanwise location with respect to the three-dimensional vortical structures before averaging across the different phases. We use the condition, u’<- urms, to locate the central plane between the counter-rotating structures, which in effect are the “legs” of the three-dimensional structure. With this condition, we effectively get a slice of the shedding cycle cutting through the “head” of the three-dimensional structure. Apart from this cut, we also get a cut between adjacent structures from the weak sweep events, with the condition u’<- urms. Using these conditions, on the phase-locked velocity fields, we effectively lock the structures in time, as well as in the spanwise direction. With this ,a clearer picture of the shedding process emerges. The flow is highly three-dimensional near reattachment and the shedding of the separation bubble is modulated in the spanwise direction owing to the three-dimensional hairpin like vortical structures in the flow. The separation bubble is seen bulged out and lifted high at locations where the head of the hairpin vortex passes, owing to the strong ejection of fluid caused by the vortical structure. On the other hand, outside the hairpin vortices, weak sweep events push the flow towards the wall and make it shallow and less prominent, with the shedding being very weak in this plane. From these observations, a three-dimensional picture of the flow is proposed.

Draft tube is very important part of hydraulic turbines. Only optimum work together with turbine can bring highest performance of this machine set. Hence it is necessary to deal with character of flow in the draft tube for different operating conditions. Efficiency of the draft tube depends on many phenomena of flow. Observing these phenomena and finding their relation with energetic transformation in the draft tube is a suitable tool to judge quality of draft tube performance. Incorrect design of the draft tube can sometimes cause lower efficiency of whole machine set. The goal of this thesis is finding the main reasons causing draft tube efficiency drop for given operating conditions.

Peripheral vascular grafts are used for the treatment of peripheral arterial disease and arteriovenous grafts for vascular access in end stage renal disease. The development of neo-intimal hyperplasia and thrombosis in the distal anastomosis remains the main reason for occlusion in that region. The local haemodynamics produced by a graft in the host vessel is believed to significantly affect endothelial function. Single spiral flow is a normal feature in medium and large sized vessels and it is induced by the anatomical structure and physiological function of the cardiovascular system. Grafts designed to generate a single spiral flow in the distal anastomosis have been introduced in clinical practice and are known as spiral grafts. In this work, spiral peripheral vascular and arteriovenous grafts were compared with conventional grafts using ultrasound and computational methods to identify their haemodynamic differences. Vascular-graft flow phantoms were developed to house the grafts in different surgical configurations. Mimicking components, with appropriate acoustic properties, were chosen to minimise ultrasound beam refraction and distortion. A dual-beam two-dimensional vector Doppler technique was developed to visualise and quantify vortical structures downstream of each graft outflow in the cross-flow direction. Vorticity mapping and measurements of circulation were acquired based on the vector Doppler data. The flow within the vascular-graft models was simulated with computed tomography based image-guided modelling for further understanding of secondary flow motions and comparison with the experimental results. The computational assessments provided a three-dimensional velocity field in the lumen of the models allowing a range of fluid dynamic parameters to be predicted. Single- or double-spiral flow patterns consisting of a dominant and a smaller vortex were detected in the outflow of the spiral grafts. A double- triple- or tetra-spiral flow pattern was found in the outflow of the conventional graft, depending on model configuration and Reynolds number. These multiple-spiral patterns were associated with increased flow stagnation, separation and instability, which are known to be detrimental for endothelial behaviour. Increased in-plane mixing and wall shear stress, which are considered atheroprotective in normal vessels, were found in the outflow of the spiral devices. The results from the experimental approach were in agreement with those from the computational approach. This study applied ultrasound and computational methods to vascular-graft phantoms in order to characterise the flow field induced by spiral and conventional peripheral vascular and arteriovenous grafts. The results suggest that spiral grafts are associated with advanced local haemodynamics that may protect endothelial function and thereby may prevent their outflow anastomosis from neo-intimal hyperplasia and thrombosis. Consequently this work supports the hypothesis that spiral grafts may decrease outflow stenosis and hence improve patency rates in patients.