A three dimensional implicit multiblock Navier-Stokes solver for hovering rotor vortical flow simulations has been developed. The governing equations used are cast in an attached blade rotating frame. Two formulations of the governing equations using the relative or absolute velocity as variables respectively are employed and investigated. The Osher's approximate Riemann solver is used for the convective fluxes evaluation. A modified MUSCL scheme is employed for improving the accuracy of the discretisation for the in viscid fluxes. A Block Incomplete Lower and Upper Decomposition (BILU) is adopted for solving the linear system resulted from the use of an implicit scheme. Special treatment for the terms, including extra flux terms and source terms, arising from the non-inertial reference system are implemented. A multiblock technique is used to obtain the exibility for quality grid generation. The suitability of different grid topologies for vortex wake capturing is demonstrated. Numerical tests show that significant improvement in computational efficiency is achieved by utilising the BILU implicit scheme in both fixed wing and hovering rotor calculations. Numerical simulations also demonstrate Navier-Stokes solutions give more accurate results than that from Euler solutions, especially in transonic tip speed cases. Computed results including surface pressure distributions and tip vortex trajectories are compared with the experimental data, which shows that the developed solver and the numerical scheme can simulate hovering rotor flows with good accuracy.

Magnetic resonance imaging (MRI) is an useful technique that provides a noninvasive method in clinical applications. For the patient of tetralogy of Fallot (TOF) after repaired, turbulence and regurgitation in blood flow may appear in pulmonary arteries. In this study, phase contrast MR imaging was applied and vortical flow patterns in the pulmonary arteries of patients after repair of TOF has been investigated. There are two major part of this study. Firstly we simulated vortical flow patterns of star, focus, and saddle which are most frequently appeared in blood flow. Quadrant index has been proposed for pattern analysis. In the second part we applied these parameters to in vivo data of repaired TOF patients, and compared with other parameters such as vorticity, coefficient of variance (CV), and regurgitant fraction (RF). Our result shows that the linear correlation between the mean of CV of velocity and mean of CV of vorticity in right pulmonary artery (RPA) as well as pulmonary trunk (PT) is larger than that in left pulmonary artery (LPA). This study shows that vorticity may provide some useful information of flow patterns and therefore helps doctors in clinical diagnosis

<p>With this contribution to the AVT-113/VFE-2 task group it was possible to prove the feasibility of high Reynolds number CFD computations to resolve and thus better understand the peculiar dual vortex system encountered on the VFE-2 blunt leading edge delta wing. Initial investigations into this phenomenon seemed to undermine the hypothesis, that the formation of the inner vortex system relies on the laminar state of the boundary layer at separation onset. As a result of this research, this initial hypothesis had to be expanded to account also for high Reynolds number cases, where a laminar boundary layer status at separation onset could be excluded. Furthermore, the data published in the same context shows evidence of secondary separation under the inner primary vortex. This further supports the supposition of a different generation mechanism of the inner vortical system other than a pure development out of a possibly laminar separation bubble. The unsteady computations performed on numerical grids with different levels of refinement led furthermore to the establishment of internal guidelines specific to the DES approach.</p>

With this contribution to the AVT-113/VFE-2 task group it was possible to prove the feasibility of high Reynolds number CFD computations to resolve and thus better understand the peculiar dual vortex system encountered on the VFE-2 blunt leading edge delta wing. Initial investigations into this phenomenon seemed to undermine the hypothesis, that the formation of the inner vortex system relies on the laminar state of the boundary layer at separation onset. As a result of this research, this initial hypothesis had to be expanded to account also for high Reynolds number cases, where a laminar boundary layer status at separation onset could be excluded. Furthermore, the data published in the same context shows evidence of secondary separation under the inner primary vortex. This further supports the supposition of a different generation mechanism of the inner vortical system other than a pure development out of a possibly laminar separation bubble. The unsteady computations performed on numerical grids with different levels of refinement led furthermore to the establishment of internal guidelines specific to the DES approach. / QC 20101111

The unsteady behavior of a tip leakage flow downstream of a simulated axial compressor rotor has been studied. The Virginia Tech low speed linear cascade wind tunnel was adapted to model the unsteady tip leakage flow produced by a rotor operating in the vortical wakes of a set of stator vanes. The cascade, consisting of 8 GE rotor B blades, has adjustable tip gap, inlet angle of 65.1 degrees, turning angle of 11.8 degrees and solidity of 1.076. The cascade Reynolds number, based on blade chord, was 393,000. A moving end wall was used to simulate the relative motion between rotor and casing, and vortex generators attached to the moving end wall were used to produce an idealized periodic unsteady vortical inflow similar to that shed by the junction of a row of inlet guide vanes. Measurements of the vortical inflow to the cascade produced by the generators and of the mean blade loading at the mid span are presented. The periodic and aperiodic behavior of the tip leakage flow downstream of the cascade, produced by this vortical disturbance, is also presented using phase and time averaged 3-component turbulence and pressure fluctuation measurements. These measurements are made for tip gap from 0.83% to 3.3% chord and streamwise locations from 0.772% to 1.117% blade spacing axially downstream of the cascade. The phase averaged inflow measurements reveal that the inflow produced by the vortex generators consists of a pair asymmetric counter-rotating vortices embedded in a thin (4.6% chord) endwall boundary layer. The vortices extend some 7.4% chord from the end wall. Their strength is about two orders smaller than the typical circulation of the tip leakage vortices produced by the cascade. Phase averaged single point three component hot-wire measurements downstream of the cascade reveal that the vortical inflow is, however, capable of producing significant large scale fluctuations in the size, strength, structure and position of the tip leakage vortex. These effects increase in magnitude with increase of tip gap. For small tip gaps these effects appear to be due to simple superposition between the inflow vortices and the tip leakage vortex. However for larger tip gaps these effects appear primarily a consequence of the inflow vortices interfering with the shedding of circulation from the blade tip. The fact that the circulation fluctuation is consistent with the inviscid unsteady loading prediction suggests that the inviscid response may be a major mechanism for generating the tip leakage unsteadiness. Although there is large periodic fluctuation in the tip leakage flow disturbed by the inflow, there is a larger aperiodic component. Two point correlation measurements and linear stochastic estimation are used to reveal the structure of this aperiodic part for a tip gap of 3.3% chord. The aperiodic fluctuation, containing most of the turbulence energy, is found appearing to be organized structures in large scale, and making the estimated instantaneous velocity field significantly different from the phase averaged periodic velocity field. Phase averaged pressure fluctuation measurements made using a microphone in the tip leakage vortex downstream of the cascade reveal that there are significant periodic fluctuating pressure waves and intense mean square fluctuation of the aperiodic fluctuating pressure. They are consistent with the measured periodic flow and aperiodic flow field respectively. These microphone measurements are validated using fluctuating pressure gradient estimates determined from the hot-wire measurements. / Ph. D.

<p>The overall goal of the research presented in this thesis is to extend the physical understanding of the unsteady external aerodynamics associated with highly maneuverable delta-wing aircraft by using and developing novel, more efficient computational fluid dynamics (CFD) tools. More specific, the main purpose is to simulate and better understand the basic fluid phenomena, such as vortex breakdown, that limit the performance of delta-wing aircraft. The problem is approached by going from the most simple aircraft configuration - a pure delta wing - to more complex configurations. As the flow computations of delta wings at high angle of attack have a variety of unusual aspects that make accurate predictions challenging, best practices for the CFD codes used are developed and documented so as to raise their technology readiness level when applied to this class of flows.</p><p>Initially, emphasis is put on subsonic steady-state CFD simulations of stand-alone delta wings to keep the phenomenon of vortex breakdown as clean as possible. For half-span models it is established that the essential characteristics of vortex breakdown are captured by a structured CFD code. The influence of viscosity on vortex breakdown is studied and numerical results for the aerodynamic coefficients, the surface pressure distribution and breakdown locations are compared to experimental data where possible.</p><p>In a second step, structured grid generation issues, numerical aspects of the simulation of this nonlinear type of flow and the interaction of a forebody with a delta wing are explored.</p><p>Then, on an increasing level of complexity, time-accurate numerical studies are performed to resolve the unsteady flow field over half and full-span, stationary delta wings at high angle of attack. Both Euler and Detached Eddy Simulations (DES) are performed to predict the streamwise oscillations of the vortex breakdown location about some mean position, asymmetry in the breakdown location due to the interaction between the left and right vortices, as well as the rotation of the spiral structure downstream of breakdown in a time-accurate manner. The computed flow-field solutions are visualized and analyzed in a virtual-reality environment.</p><p>Ultimately, steady-state and time-dependent simulations of a full-scale fighter-type aircraft configuration in steady flight are performed using the advanced turbulence models and the detached-eddy simulation capability of an edge-based, unstructured flow solver. The computed results are compared to flight-test data.</p><p>The thesis also addresses algorithmic efficiency and presents a novel implicit-explicit algorithm, the Recursive Projection Method (RPM), for computations of both steady and unsteady flows. It is demonstrated that RPM can accelerate such computations by up to 2.5 times.</p>

The overall goal of the research presented in this thesis is to extend the physical understanding of the unsteady external aerodynamics associated with highly maneuverable delta-wing aircraft by using and developing novel, more efficient computational fluid dynamics (CFD) tools. More specific, the main purpose is to simulate and better understand the basic fluid phenomena, such as vortex breakdown, that limit the performance of delta-wing aircraft. The problem is approached by going from the most simple aircraft configuration - a pure delta wing - to more complex configurations. As the flow computations of delta wings at high angle of attack have a variety of unusual aspects that make accurate predictions challenging, best practices for the CFD codes used are developed and documented so as to raise their technology readiness level when applied to this class of flows. Initially, emphasis is put on subsonic steady-state CFD simulations of stand-alone delta wings to keep the phenomenon of vortex breakdown as clean as possible. For half-span models it is established that the essential characteristics of vortex breakdown are captured by a structured CFD code. The influence of viscosity on vortex breakdown is studied and numerical results for the aerodynamic coefficients, the surface pressure distribution and breakdown locations are compared to experimental data where possible. In a second step, structured grid generation issues, numerical aspects of the simulation of this nonlinear type of flow and the interaction of a forebody with a delta wing are explored. Then, on an increasing level of complexity, time-accurate numerical studies are performed to resolve the unsteady flow field over half and full-span, stationary delta wings at high angle of attack. Both Euler and Detached Eddy Simulations (DES) are performed to predict the streamwise oscillations of the vortex breakdown location about some mean position, asymmetry in the breakdown location due to the interaction between the left and right vortices, as well as the rotation of the spiral structure downstream of breakdown in a time-accurate manner. The computed flow-field solutions are visualized and analyzed in a virtual-reality environment. Ultimately, steady-state and time-dependent simulations of a full-scale fighter-type aircraft configuration in steady flight are performed using the advanced turbulence models and the detached-eddy simulation capability of an edge-based, unstructured flow solver. The computed results are compared to flight-test data. The thesis also addresses algorithmic efficiency and presents a novel implicit-explicit algorithm, the Recursive Projection Method (RPM), for computations of both steady and unsteady flows. It is demonstrated that RPM can accelerate such computations by up to 2.5 times. / QC 20101019

The complex multiphase flow physics of spray-swirl interaction in both reacting and non-reacting environment is of fundamental and applied significance for a wide variety of applications ranging from gas turbine combustors to pharmaceutical drug nebulizers. Understanding the intricate dynamics between this two phase flow field is pivotal for enhancing mixing characteristics, reducing pollutant emissions and increasing the combustion efficiency in next generation combustors. The present work experimentally investigates the near and far-field break-up, dispersion and coalescence characteristics of a hollow cone spray in an unconfined, co¬annular isothermal swirling air jet environment. The experiments were conducted using an axial-flow hollow cone spray nozzle having a 0.5 mm orifice. Nozzle injection pressure (PN = 1 bar) corresponding to a Reynolds number at nozzle exit ReN = 7900 used as the test setting. At this setting, the operating Reynolds number of the co-annular swirling air stream number (Res) was varied in four distinct steps, i.e. Res = 1600, 3200, 4800 and 5600. Swirl was imparted to the co¬axial flow using a guided vane swirler with blade angle of Ф=45° (corresponding geometric swirl number SG = 0.8). Two types of laser diagnostic techniques were utilized: Particle / Droplet imaging analysis (PDIA) and shadowgraph to study the underlying physical mechanisms involved in the primary breakup, dispersion and coalescence dynamics of the spray. Measurements were made in the spray in both axial and radial directions and they indicate that Sauter Mean Diameter (SMD) in radial direction is highly reliant on the intensity of swirl imparted to the spray. The spray is subdivided into two zones as function of swirl in axial and radial direction: (1) near field of the nozzle (ligament regime) where variation in SMD arises predominantly due to primary breakup of liquid films (2) far-field of the nozzle where dispersion and collision induced coalescence of droplets is dominant. Each regime has been analyzed meticulously, by computing probability of primary break-up, probability of coalescence and spatio-temporal distribution of droplets which gives probabilistic estimate of aforementioned governing phenomena. In addition to this, spray global length scale parameters such as spray cone angle, break-up length, wavelength of liquid film has been characterized by varying Res while maintaining constant ReN.

This report documents the experimental study for four different three-dimensional turbulent flows. The investigation focuses on near wall measurements in these flows. Several experimental techniques are used in the studies; however, the bulk of the investigation focuses on a three-orthogonal-velocity-component fiber-optic laser Doppler anemometer (3D-LDA) system. The control volume of the 3D-LDA is on the order of 50 micro-meter in size, or a y⁺ distance of around 2.3 units (using average values of U_&#964 and &#957; from the experiment). An auxiliary small boundary layer wind tunnel (auxiliary tunnel) and a low speed linear compressor cascade wind tunnel (cascade tunnel) are utilized in this study. One of four flow experiments is done in the auxiliary tunnel the other three are in the cascade tunnel. The first three-dimensional turbulent flow is a vortical flow created by two half-delta wing vortex generators. Near wall secondary flow features are found. The second flow is an investigation of the first quarter chord tip gap flow in the cascade tunnel. Strong three-dimensional phenomena are found. The third flow investigated is the inflow to the compressor cascade with the moving wall. The experiment records shear layer interaction between the upstream flow and moving wall. Finally the fourth flow investigated is the inflow to the compressor cascade with the moving wall with half-delta wing vortex generators attached. Phase-averaged data reveal asymmetrical vortex structures just downstream of the vortex generators. This is the first time any near wall data has been taken on any of these flows. / Master of Science

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.