A major concern in the operation of present-day gas turbine engines is the ingestion of hot mainstream gas into rotor-stator disk cavities of the high-pressure turbine stages. Although the engines require high gas temperature at turbine entry for good performance efficiency, the ingested gas shortens the lives of the cavity internals, particularly that of the rotor disks. Steps such as installing seals at the disk rims and injecting purge (secondary) air bled from the compressor discharge into the cavities are implemented to reduce the gas ingestion. Although there are advantages to the above-mentioned steps, the performance of a gas turbine engine is diminished by the purge air bleed-off. This then requires that the cavity sealing function be achieved with as low a purge air supply rate as possible. This, in turn, renders imperative an in-depth understanding of the pressure and velocity fields in the main gas path and within the disk cavities. In this work, experiments were carried out in a model 1.5-stage (stator-rotor-stator) axial air turbine to study the ingestion of main air into the aft, rotor-stator, disk cavity. The cavity featured rotor and stator rim seals with radial clearance and axial overlap and an inner labyrinth seal. First, time-average static pressure distribution was measured in the main gas path upstream and downstream of the rotor as well as in the cavity to ensure that a nominally steady run condition had been achieved. Main gas ingestion was determined by measuring the concentration distribution of tracer gas (CO2) in the cavity. To map the cavity fluid velocity field, particle image velocimetry was employed. Results are reported for two main air flow rates, two rotor speeds, and four purge air flow rates.

Passive flow control achieved by surface dimpling can be an effective strategy for reducing drag around bluff bodies - an example of substantial popular interest being the flow around a golf ball. While the general effect of dimples causing a delay of boundary layer separation is well known, the mechanisms contributing to this phenomena are subtle and not thoroughly understood. Numerical models offer a powerful approach for studying drag reduction, however simulation strategies are challenged by complex geometries, and in applications the introduction of ad hoc turbulence models which introduce additional uncertainty. These and other factors provide much of the motivation for the current study, which focused on the numerical simulations of the flow over a simplified configuration consisting of a dimpled flat plate. The principal goals of the work are to understand the performance of the numerical methodology, and gain insight into the underlying physics of the flow. Direct numerical simulation of the incompressible Navier-Stokes equations using a fractional step method was employed, with the dimpled flat plate represented using an immersed boundary method. The dimple geometry utilizes a fixed dimple aspect ratio, with dimples arranged in a single spanwise row. The grid sizes considered ranged from approximately 3 to 99 million grid points. Reynolds numbers of 3000 and 4000 based on the inlet laminar boundary layer thickness were simulated. A turbulent boundary layer was induced downstream of the dimples for Reynolds numbers which did not transition for the flow over an undimpled flat plate. First and second order statistics of the boundary layer that develops agree reasonably well with those for turbulent channel flow and flat plate boundary layers in the sublayer and buffer layers, but differ in the outer layer. Inspection of flow visualizations suggest that early transition is promoted by thinning of the boundary layer, initiation of shear layer instabilities over the dimples, flow separation and reattachment, and tripping of the boundary layer at the trailing edge of the dimples.

Efficient performance of gas turbines depends, among several parameters, on the mainstream gas entry temperature. At the same time, transport of this high temperature gas into the rotor-stator cavities of turbine stages affects the durability of rotor disks. This transport is usually countered by installing seals on the rotor and stator disk rims and by pressurizing the cavities by injecting air (purge gas) bled from the compressor discharge. The configuration of the rim seals influences the magnitude of main gas ingestion as well as the interaction of the purge gas with the main gas. The latter has aerodynamic and hub endwall heat transfer implications in the main gas path. In the present work, experiments were performed on model single-stage and 1.5-stage axial-flow turbines. The turbines featured vanes, blades, and rim seals on both the rotor and stator disks. Three different rim seal geometries, viz., axially overlapping radial clearance rim seals for the single-stage turbine cavity and the 1.5-stage turbine aft cavity, and a rim seal with angular clearance for the single-stage turbine cavity were studied. In the single-stage turbine, an inner seal radially inboard in the cavity was also provided; this effectively divided the disk cavity into a rim cavity and an inner cavity. For the aft rotor-stator cavity of the 1.5-stage turbine, a labyrinth seal was provided radially inboard, again creating a rim cavity and an inner cavity. Measurement results of time-average main gas ingestion into the cavities using tracer gas (CO2), and ensemble-averaged trajectories of the purge gas flowing out through the rim seal gap into the main gas path using particle image velocimetry are presented. For both turbines, significant ingestion occurred only in the rim cavity. The inner cavity was almost completely sealed by the inner seal, at all purge gas flow rates for the single-stage turbine and at the higher purge gas flow rates for 1.5-stage turbine. Purge gas egress trajectory was found to depend on main gas and purge gas flow rates, the rim seal configuration, and the azimuthal location of the trajectory mapping plane with respect to the vanes.

Wind measurements are fundamental inputs for the evaluation of potential energy yield and performance of wind farms. Three-dimensional scanning coherent Doppler lidar (CDL) may provide a new basis for wind farm site selection, design, and control. In this research, CDL measurements obtained from multiple wind energy developments are analyzed and a novel wind farm control approach has been modeled. The possibility of using lidar measurements to more fully characterize the wind field is discussed, specifically, terrain effects, spatial variation of winds, power density, and the effect of shear at different layers within the rotor swept area. Various vector retrieval methods have been applied to the lidar data, and results are presented on an elevated terrain-following surface at hub height. The vector retrieval estimates are compared with tower measurements, after interpolation to the appropriate level. CDL data is used to estimate the spatial power density at hub height. Since CDL can measure winds at different vertical levels, an approach for estimating wind power density over the wind turbine rotor-swept area is explored. Sample optimized layouts of wind farm using lidar data and global optimization algorithms, accounting for wake interaction effects, have been explored. An approach to evaluate spatial wind speed and direction estimates from a standard nested Coupled Ocean and Atmosphere Mesoscale Prediction System (COAMPS) model and CDL is presented. The magnitude of spatial difference between observations and simulation for wind energy assessment is researched. Diurnal effects and ramp events as estimated by CDL and COAMPS were inter-compared. Novel wind farm control based on incoming winds and direction input from CDL's is developed. Both yaw and pitch control using scanning CDL for efficient wind farm control is analyzed. The wind farm control optimizes power production and reduces loads on wind turbines for various lidar wind speed and direction inputs, accounting for wind farm wake losses and wind speed evolution. Several wind farm control configurations were developed, for enhanced integrability into the electrical grid. Finally, the value proposition of CDL for a wind farm development, based on uncertainty reduction and return of investment is analyzed.

Production from a high pressure gas well at a high production-rate encounters the risk of operating near the choking condition for a compressible flow in porous media. The unbounded gas pressure gradient near the point of choking, which is located near the wellbore, generates an effective tensile stress on the porous rock frame. This tensile stress almost always exceeds the tensile strength of the rock and it causes a tensile failure of the rock, leading to wellbore instability. In a porous rock, not all pores are choked at the same flow rate, and when just one pore is choked, the flow through the entire porous medium should be considered choked as the gas pressure gradient at the point of choking becomes singular. This thesis investigates the choking condition for compressible gas flow in a single microscopic pore. Quasi-one-dimensional analysis and axisymmetric numerical simulations of compressible gas flow in a pore scale varicose tube with a number of bumps are carried out, and the local Mach number and pressure along the tube are computed for the flow near choking condition. The effects of tube length, inlet-to-outlet pressure ratio, the number of bumps and the amplitude of the bumps on the choking condition are obtained. These critical values provide guidance for avoiding the choking condition in practice.

ABSTRACT The heat recovery steam generator (HRSG) is a key component of Combined Cycle Power Plants (CCPP). The exhaust (flue gas) from the CCPP gas turbine flows through the HRSG − this gas typically contains a high concentration of NO and cannot be discharged directly to the atmosphere because of environmental restrictions. In the HRSG, one method of reducing the flue gas NO concentration is to inject ammonia into the gas at a plane upstream of the Selective Catalytic Reduction (SCR) unit through an injection grid (AIG); the SCR is where the NO is reduced to N2 and H2O. The amount and spatial distribution of the injected ammonia are key considerations for NO reduction while using the minimum possible amount of ammonia. This work had three objectives. First, a flow network model of the Ammonia Flow Control Unit (AFCU) was to be developed to calculate the quantity of ammonia released into the flue gas from each AIG perforation. Second, CFD simulation of the flue gas flow was to be performed to obtain the velocity, temperature, and species concentration fields in the gas upstream and downstream of the SCR. Finally, performance characteristics of the ammonia injection system were to be evaluated. All three objectives were reached. The AFCU was modeled using JAVA - with a graphical user interface provided for the user. The commercial software Fluent was used for CFD simulation. To evaluate the efficacy of the ammonia injection system in reducing the flue gas NO concentration, the twelve butterfly valves in the AFCU ammonia delivery piping (risers) were throttled by various degrees in the model and the NO concentration distribution computed for each operational scenario. When the valves were kept fully open, it was found that it led to a more uniform reduction in NO concentration compared to throttling the valves such that the riser flows were equal. Additionally, the SCR catalyst was consumed somewhat more uniformly, and ammonia slip (ammonia not consumed in reaction) was found lower. The ammonia use could be decreased by 10 percent while maintaining the NO concentration limit in the flue gas exhausting into the atmosphere.

Multi-touch tablets and smart phones are now widely used in both workplace and consumer settings. Interacting with these devices requires hand and arm movements that are potentially complex and poorly understood. Experimental studies have revealed differences in performance that could potentially be associated with injury risk. However, underlying causes for performance differences are often difficult to identify. For example, many patterns of muscle activity can potentially result in similar behavioral output. Muscle activity is one factor contributing to forces in tissues that could contribute to injury. However, experimental measurements of muscle activity and force for humans are extremely challenging. Models of the musculoskeletal system can be used to make specific estimates of neuromuscular coordination and musculoskeletal forces. However, existing models cannot easily be used to describe complex, multi-finger gestures such as those used for multi-touch human computer interaction (HCI) tasks. We therefore seek to develop a dynamic musculoskeletal simulation capable of estimating internal musculoskeletal loading during multi-touch tasks involving multi digits of the hand, and use the simulation to better understand complex multi-touch and gestural movements, and potentially guide the design of technologies the reduce injury risk. To accomplish these, we focused on three specific tasks. First, we aimed at determining the optimal index finger muscle attachment points within the context of the established, validated OpenSim arm model using measured moment arm data taken from the literature. Second, we aimed at deriving moment arm values from experimentally-measured muscle attachments and using these values to determine muscle-tendon paths for both extrinsic and intrinsic muscles of middle, ring and little fingers. Finally, we aimed at exploring differences in hand muscle activation patterns during zooming and rotating tasks on the tablet computer in twelve subjects. Towards this end, our musculoskeletal hand model will help better address the neuromuscular coordination, safe gesture performance and internal loadings for multi-touch applications.

The flow of liquid PDMS (10:1 v/v base to cross-linker ratio) in open, rectangular silicon micro channels, with and without a hexa-methyl-di-silazane (HMDS) or poly-tetra-fluoro-ethylene (PTFE) (120 nm) coat, was studied. Photolithographic patterning and etching of silicon wafers was used to create micro channels with a range of widths (5-50 μm) and depths (5-20 μm). The experimental PDMS flow rates were compared to an analytical model based on the work of Lucas and Washburn. The experimental flow rates closely matched the predicted flow rates for channels with an aspect ratio (width to depth), p, between one and two. Flow rates in channels with p less than one were higher than predicted whereas the opposite was true for channels with p greater than two. The divergence between the experimental and predicted flow rates steadily increased with increasing p. These findings are rationalized in terms of the effect of channel dimensions on the front and top meniscus morphology and the possible deviation from the no-slip condition at the channel walls at high shear rates.

In addition, a preliminary experimental setup for calibration tests on ultrasensitive PDMS cantilever beams is reported. One loading and unloading cycle is completed on a microcantilever PDMS beam (theoretical stiffness 0.5 pN/ µm). Beam deflections are actuated by adjusting the buoyancy force on the beam, which is submerged in water, by the addition of heat. The expected loading and unloading curve is produced, albeit with significant noise. The experimental results indicate that the beam stiffness is a factor of six larger than predicted theoretically. One probable explanation is that the beam geometry may change when it is removed from the channel after curing, making assumptions about the beam geometry used in the theoretical analysis inaccurate. This theory is bolstered by experimental data discussed in the report. Other sources of error which could partially contribute to the divergent results are discussed. Improvements to the experimental setup for future work are suggested.

In order to achieve higher gas turbine efficiency, the main gas temperature at turbine inlet has been steadily increased from approximately 900°C to about 1500°C over the last few decades. This temperature is higher than the maximum acceptable temperature for turbine internals. The hot main gas may get ingested into the space between rotor and stator, the rotor-stator disk cavity in a stage because of the pressure differential between main gas annulus and the disk cavity. To reduce this ingestion, the disk cavity is equipped with a rim seal; additionally, secondary (purge) air is supplied to the cavity. Since the purge air is typically bled off the compressor discharge, this reducing the overall gas turbine efficiency, much research has been carried out to estimate the minimum purge flow necessary (c_w,min) for complete sealing of disk cavities.

In this work, experiments have been performed in a subscale single-stage axial turbine featuring vanes, blades and an axially-overlapping radial-clearance seal at the disk cavity rim. The turbine stage is also equipped with a labyrinth seal radially inboard. The stage geometry and the experimental conditions were such that the ingestion into the disk cavity was driven by the pressure asymmetry in the main gas annulus. In the experiments, time-averaged static pressure was measured at several locations in the main annulus and in the disk cavity; the pressure differential between a location on the vane platform close to lip (this being the rim seal part on the stator) and a location in the 'seal region' in the cavity is considered to be the driving potential for both ingestion and egress. Time-averaged volumetric concentration of the tracer gas (CO₂) in the purge air supplied was measured at multiple radial locations on the stator surface. The pressure and ingestion data were then used to calculate the ingestion and egress discharge coefficients for a range of purge flow rates, employing a simple orifice model of the rim seal. For the experiments performed, the egress discharge coefficient increased and the ingestion discharge coefficient decreased with the purge air flow rate. A method for estimation of c_w,min is also proposed.

Hydraulic fracturing is an effective technique used in well stimulation to increase petroleum well production. A combination of multi-stage hydraulic fracturing and horizontal drilling has led to the recent boom in shale gas production which has changed the energy landscape of North America.

During the fracking process, highly pressurized mixture of water and proppants (sand and chemicals) is injected into to a crack, which fractures the surrounding rock structure and proppants help in keeping the fracture open. Over a longer period, however, these fractures tend to close due to the difference between the compressive stress exerted by the reservoir on the fracture and the fluid pressure inside the fracture. During production, fluid pressure inside the fracture is reduced further which can accelerate the closure of a fracture.

In this thesis, we study the stress distribution around a hydraulic fracture caused by fluid production. It is shown that fluid flow can induce a very high hoop stress near the fracture tip. As the pressure gradient increases stress concentration increases. If a fracture is very thin, the flow induced stress along the fracture decreases, but the stress concentration at the fracture tip increases and become unbounded for an infinitely thin fracture.

The result from the present study can be used for studying the fracture closure problem, and ultimately this in turn can lead to the development of better proppants so that prolific well production can be sustained for a long period of time.