The current demands for higher performance in gas turbine engines can be reached by raising combustion temperatures to increase thermal efficiency. Hot combustion temperatures create a harsh environment which leads to the consideration of the durability of the combustor and turbine sections. Improvements in durability can be achieved through understanding the interactions between the combustor and turbine. The flow field at a combustor exit shows non-uniformities in pressure, temperature, and velocity in the pitch and radial directions. This inlet profile to the turbine can have a considerable effect on the development of the secondary flows through the vane passage. This thesis presents a computational study of the flow field generated in a non-reacting gas turbine combustor and how that flow field convects through the downstream stator vane. Specifically, the effect that the combustor flow field had on the secondary flow pattern in the turbine was studied. Data from a modern gas turbine engine manufacturer was used to design a realistic, low speed, large scale combustor test section. This thesis presents the results of computational simulations done in parallel with experimental simulations of the combustor flow field. In comparisons of computational predictions with experimental data, reasonable agreement of the mean flow and general trends were found for the case without dilution jets. The computational predictions of the combustor flow with dilution jets indicated that the turbulence models under-predicted jet mixing. The combustor exit profiles showed non-uniformities both radially and circumferentially, which were strongly dependent on dilution and cooling slot injection. The development of the secondary flow field in the turbine was highly dependent on the incoming total pressure profile. For a case with a uniform inlet pressure in the near-wall region no leading edge vortex was formed. The endwall heat transfer was found to also depend strongly on the secondary flow field, and therefore on the incoming pressure profile from the combustor. / Master of Science

Detonation cycles are identified as an efficient alternative to the Brayton cycles used in power and propulsion applications. Rotating Detonation Engine (RDE) operating on a detonation cycle works by compressing the working fluid across a detonation wave, thereby reducing the number of compressor stages required in the thermodynamic cycle. Numerical analyses of RDEs are flexible in understanding the flow field within the RDE, however, three-dimensional analyses are expensive due to the differences in time-scale required to resolve the combustion process and flow-field. The alternate two-dimensional analyses are generally modeled with perfectly premixed fuel injection and do not capture the effects of improper mixing arising due to discrete injection of fuel and oxidizer into the chamber. To model realistic injection in a 2-D analysis, the current work uses an approach in which, a Probability Density Function (PDF) of the fuel mass fraction at the chamber inlet is extracted from a 3-D, cold-flow simulation and is used as an inlet boundary condition for fuel mass fraction in the 2-D analysis. The 2-D simulation requires only 0.4% of the CPU hours for one revolution of the detonation compared to an equivalent 3-D simulation. Using this method, a perfectly premixed RDE is comparing with a non-premixed case. The performance is found to vary between the two cases. The mean detonation velocities, time-averaged static pressure profiles are found to be similar between the two cases, while the local detonation velocities and peak pressure values vary in the non-premixed case due to local pockets fuel rich/lean mixtures. The mean detonation cell sizes are similar, but the distribution in the non-premixed case is closer due to stronger shock structures. An analytical method is used to check the effects of fuel-product stratification and heat loss from the RDE and these effects adversely affect the local detonation velocity. Overall, this method of modeling captures the complex physics in an RDE with the advantage of reduced computational cost and therefore can be used for design and diagnostic purposes. / Master of Science / The conventional Brayton cycle used in power and propulsion applications is highly optimized, at cycle and component levels. In pursuit of higher thermodynamic efficiency, detonation cycles are identified as an efficient alternative and gained increased attention in the scientific community. In a Rotating Detonation Engine (RDE), which is based on the detonation cycle, the compression of gases occurs across a shock wave. This method of achieving high compression ratios reduces the number of compressor stages required for operation. In an RDE (where combustion occurs between two coaxial cylinders), the fuel and oxidizer are injected axially into the combustion chamber where the detonation is initiated. The resultant detonation wave spins continuously in the azimuthal direction, consuming fresh fuel mixture. The combustion products expand and exhaust axially providing thrust/mechanical energy when coupled with a turbine. Numerical analyses of RDEs are flexible over experimental analysis, in terms of understanding the flow physics and the physical/chemical processes occurring within the engine. However, three-dimensional numerical analyses are computationally expansive, and therefore demanding an equivalent, efficient two-dimensional analysis. In most RDEs, fuel and oxidizer are injected from separate plenums into the chamber. This type of injection leads to inhomogeneity of the fuel-air mixture within the RDE which adversely affects the performance of the engine. The current study uses a novel method to effectively capture these physics in a 2-D numerical analysis. Furthermore, the performance of the combustor is compared between perfectly premixed injection and discrete, non-premixed injection. The method used in this work can be used for any injector design and is a powerful/efficient way to numerically analyze a Rotating Detonation Engine.

The detonation-based engine has the key advantage of increased thermodynamic efficiency over the traditional constant pressure combustor. These detonation-based engines are also known as Pressure Gain Combustion systems (PGC) and Rotating Detonation Combustor (RDC) is a form of PGC, in which the detonation wave propagates azimuthally around an annular combustor. Prior researchers have performed a high fidelity 3-D numerical simulation of a rotating detonation combustor (RDC) to understand the flow physics such as detonation wave velocity, pressure profile, wave structure; however, performing these 3-D simulations is computationally expensive. 2-D simulations are a potential alternative to reduce computational cost. In most RDCs, fuel and oxidizer are injected discretely from separate plenums, and this discrete fuel/air injection results in inhomogeneous mixing within the domain. Due to the discrete fuel injection locations, fuel/oxidizer will stratify to form localized pockets of rich and lean mixtures. The motivation of the present study is to investigate the impact of unmixedness and stratification length scales on the performance of an RDC using a 2-D numerical approach. Unmixedness, which is defined as the standard deviation of equivalence ratio normalized by the mean global equivalence ratio, is a measure of the degree of fuel-oxidizer inhomogeneity. To model the effect of unmixedness in a 2-D domain, a lognormal distribution of the fuel mass fraction is generated with a mean equivalence ratio of 1 and varying standard deviations at the inlet boundary as a numerical source term. Moreover, to model the effects of stratification length scales, fuel mass fraction at the inlet boundary cells is bundled for a given length scale, and the mass fractions for these bundles are updated based on the lognormal distribution after every three-time steps. Using this methodology, 2-D numerical analyses are carried out to investigate the performance of an RDC for an H2-air mixture with varying unmixedness and stratification length scales. Results show that mean detonation velocity decreases and wave speed variation increases with an increase in unmixedness. However, with an increase in stratification length scale mean velocity remain relatively unchanged but variation in local velocity increases. The detonation wave front corrugation also increases with an increase in mixture inhomogeneity. The mean detonation cell size increases with an increase in unmixedness. The cell shape becomes more distorted and irregular with an increase in stratification length scale and unmixedness. The combined effect of unmixedness and stratification length scale leads to a decrease in pressure gain. Overall, this concept is able to elucidate the effects of varying unmixedness and stratification length scales on the performance of an RDC. / Master of Science / Pressure Gain Combustion (PGC) system has gained significant focus in recent years due to its increased thermodynamic efficiency over a constant pressure Brayton Cycle. Rotating Detonation Combustor (RDC) is a type of PGC system, which is thermodynamically more efficient than the conventional gas turbine combustor. One of the main aspects of the detonation process is the rapid burning of the fuel-oxidizer mixture, which occurs so fast that there is not enough time for pressure to equilibrate. Therefore, the process is thermodynamically closer to a constant volume process rather than a constant pressure process. A constant volume cycle is thermodynamically more efficient than a constant pressure Brayton cycle. In an RDC, a mixture of fuel and air is injected axially, and a detonation wave propagates continuously through the circumferential section. Numerical simulation of an RDC provides additional flexibility over experiments in understanding the flow physics, detonation wave structure, and analyzing the physical and chemical processes involved in the detonation cycle. Prior researchers have utilized a full-scale 3-D numerical simulation for understanding the performance of an RDC. However, the major challenge with 3-D analyses is the computational expense. Thus, to overcome this, an inexpensive 2-D simulation is used to model the flow physics of an RDC. In most RDCs, the fuel and oxidizer are injected discretely from separate plenums. Due to the discrete fuel injection, the fuel/air mixture is never perfectly premixed and results in a stratified flow field. The objective of the current work is to develop a novel approach to independently investigate the effects of varying unmixedness and stratification length scales on RDC performance using a 2-D simulation.

The evaluation of how Computational Fluid Dynamics (CFD) package may be incorporated into a conceptual design method is performed. The repeatability of the CFD solution as well as the accuracy of the calculated aerodynamic coefficients and pressure distributions was also evaluated on two different wing-body models. The overall run times of three different mesh densities was also evaluated to investigate if the mesh density could be reduced enough so that the computational stage of the CFD cycle may become affordable to use in the conceptual design stage. A farfield method was derived and used in this analysis to calculate the lift and drag coefficients. The CFD solutions were also compared with two methods currently used in conceptual design - the vortex lattice based program Vorview and ACSYNT. The unstructured Euler based CFD package FELISA was used in this study. / Master of Science

Understanding the flow in stirred vessels can be useful for a wide number of industrial applications, like in mining, chemical and pharmaceutical processes. Remodeling and redesigning these processes may have a significant impact on the overall design characteristics, affecting directly product quality and maintenance costs. In most cases the flow around the rotating impeller blades interacting with stationary baffles can cause rapid changes of the flow characteristics, which lead to high levels of turbulence and higher shear rates. The flow is anisotropic and inhomogeneous over the entire volume. A better understanding and a detailed documentation of the turbulent flow field is needed in order to design stirred tanks that can meet the required operation conditions. This thesis describes efforts for accurate estimation of the velocity distribution and the turbulent characteristics (vorticity, turbulent kinetic energy, dissipation rate) in a cylindrical vessel agitated by a Rushton turbine (a disk with six flat blades) and in a tank typical of flotation cells. Results from simulations using FLUENT (a commercial CFD package) are compared with Time Resolved Digital Particle Image Velocimetry (DPIV) for baseline configurations in order to validate and verify the fidelity of the computations. Different turbulence models are used in this study in order to determine the most appropriate for the prediction of turbulent properties. Subsequently a parametric analysis of the flow characteristics as a function of the clearance height of the impeller from the vessel floor is performed for the Rushton tank as well as the flotation cell. Results are presented for both configurations along planes normal or parallel to the impeller axis, displaying velocity vector fields and contour plots of vorticity turbulent dissipation and others. Special attention is focused in the neighborhood of the impeller region and the radial jet generated there. This flow in this neighborhood involves even larger gradients and dissipation levels in tanks equipped with stators. The present results present useful information for the design of the stirring tanks and flotation cells, and provide some guidance on the use of the present tool in generating numerical solutions for such complex flow fields. / Master of Science

Inom den marina industrin har ett allt större fokus riktats mot att hitta lösningar för att minska fartygens energiförbrukning. Delvis till följd av globala trender såsom ökad miljömedvetenhet och högre bränslepriser, men framför allt på grund av nya internationella regelverk som begränsar de tillåtna utsläppen från fartyg. En åtgärd för att öka fartygs framdrivningseffektivitet är genom att installera energibesparande enheter (ESD). En Pre-Swirl Stator (PSS) är ett exempel på en sådan enhet, som består av ett antal statorblad som monteras framför propellern för att skapa en mer fördelaktig flödesregim och optimera propellerns arbetsmiljö. I denna studie genomförs en numerisk undersökning av en PSS som en potentiell lösning för att förbättra framdrivningseffektiviteten genom eftermontering på ett chemfartyg. Genom att analysera interaktionen mellan skrovets medströmsfält, statorbladen och propellern fås insikter kring hur olika designparametrar på PSS:n påverkar inflödet till propellern. Resultaten från CFD-analysen jämförs med och utan PSS i full skala för att avgöra om PSS:n har en positiv eller negativ effekt på propellerverkningsgraden. De designparametrar som undersökts är antal statorblad/designorientering, stigningsvinkel och NACA-profil. För designarbetet av PSS:n har CAD NX och CFD-programvaran STAR-CCM++ genom Kongsbergs egna HullProp Interface tillämpats. Resultaten visar att en PSS kan påverka chemfartygets framdrivningseffektivitet, där PSS:ns designparametrarna har en stor inverkan på om propellerverkningsgraden ökar eller minskar. Genom att installera en PSS kan framdrivningseffektiviteten förbättras och den största verkningsgradsökningen på 0,94 % erhölls för en asymmetrisk design. För att uppnå ökad framdrivningseffektivitet ska PSS:n kunna interagera med det inkommande flödet utan att skapa för stort motstånd. Samtidigt bör en jämn och stabil strömning av vatten genereras in mot propellern. Designparametrarna bör justeras för att undvika flödesseperation på statorbladen, eftersom detta leder till ökat motstånd och ojämn strömning in till propellern. Anfallsvinkeln mot propellerbladen får heller inte bli för stor till följd av statorbladen, då detta kan orsaka flödesseperation på propellerbladen, vilket resulterar i minskad tryckkraft och verkningsgrad. För vidare studier rekommenderas användning av mer avancerade CFD-metoder för att få en tydligare bild av den komplexa flödesdynamiken och verifiera resultaten. Detta kan leda till en bättre förståelse för flödet kring PSS och dess interaktion med propellern, innan en mer omfattande parameterstudie genomförs.

This thesis contributes an in-depth comparative study of u-tube and coaxial borehole heat exchangers. While it is widely accepted that the lower resistance of the coaxial heat exchanger should result in a performance advantage, the findings of several studies comparing the heat exchanger configurations did not definitively establish the mechanisms causing differences in performance. This study employs numerical modelling to consider heat exchangers over a broad range of time scales and under carefully controlled geometry and flow conditions, resulting in the identification of the key parameters influencing borehole heat exchanger performance. The first part of this study consists of a comparison of u-tube and coaxial heat exchangers under continuous loading. A detailed conjugate heat transfer numerical model was developed in OpenFOAM, designed to capture both short and long time scales of heat exchange, necessary to understand the nuanced differences between designs. A novel transient resistance analysis was employed to understand the dominant factors influencing performance. This study established that marginal differences exist between u-tube and coaxial borehole heat exchangers (BHEs) when operated continuously long term but that greater differences occur early in operation. The second phase of this investigation provided a framework for analysing borehole heat exchanger performance during intermittent operation, while also comparing u-tube and coaxial designs. During this study, it was found that reducing operating time, improving the the rate of the ground's recovery to its original temperature, and lowering the duty cycle improved BHE performance. Transit time was identified as a influential time scale, below which heating at the outlet was limited. Further, the benefits of operating below the transit time were mitigated by design-specific interaction between inlet and outlet flows. Finally, this study found that non-dimensionalizing operating time by transit time causes the differences between u-tube and coaxial performance to vanish, leading to the conclusion that differences in BHE performance are caused by variations in flow rather than thermal mass. / Thesis / Doctor of Philosophy (PhD) / This thesis provides an in-depth comparative study of two different designs of borehole heat exchanger, the u-tube and coaxial, which are used in geothermal applications to transfer heat to and from the ground. While many researchers anticipated that the coaxial design would perform better, several studies comparing the heat exchangers were not able to provide a clear answer about which heat exchanger performed best. This study addressed this gap by using detailed numerical simulations which showed that there was a marginal difference in performance between the two heat exchangers when operated for periods longer than a few hours, but that larger differences occurred early in operation (under 15 minutes). The results also showed that operating intermittently resulted in improvements in performance of the heat exchanger, particularly when operated for periods less than the time it takes fluid to travel the length of the piping.

The thermal efficiency of gas turbines increases with the temperature of the gas entering its turbine component. To enable high inlet temperatures, even those that far exceed the melting point of the turbine materials, the turbine must be cooled. One way is by internal cooling, where cooler air passes through U-ducts embedded inside turbine vanes and blades. Since the flow and heat transfer in these ducts are highly complicated, computational fluid dynamics (CFD) based on RANS have been used extensively to explore and assess design concepts. However, RANS have been found to be unreliable – giving accurate results for some designs but not for others. In this study, large-eddy simulations (LES) were performed for a U-duct with a trapezoidal cross section to assess four widely used RANS turbulence models: realizable k-ε (k-ε), shear-stress transport (SST), Reynolds stress model with linear pressure strain (RSM-LPS), and the seven-equation stress-omega full Reynolds stress model (RSM).<div><br></div><div>When examining the capability of steady RANS, two versions of the U-duct were examined, one with a staggered array of pin fins and one without pin fins. Results obtained for the heat-transfer coefficient (HTC) were compared with experimental measurements. The maximum relative error in the predicted “averaged” HTC was found to be 50% for k-ε and RSM-LPS, 20% for SST, and 30% for RSM-τω when there are no pin fins and 25% for k-ε, 12% for the SST and RSM-τω when there are pin fins. When there are no pin fins, all RANS models predicted a large separated flow region downstream of the turn, which the experiment does show to exist. Thus, all models predicted local distributions poorly. When there were pin fins, they behaved like guide vanes in turning the flow and confined the separation around the turn. For this configuration, all RANS models predicted reasonably well.<br></div><div><br></div><div>To understand why RANS cannot predict the HTC in the U-duct after the turn when there are no pin fins, LES were performed. To ensure that the LES is benchmark quality, verification and validation were performed via LES of a straight duct with square cross section where data from experiments and direct numerical simulation (DNS) are available. To ensure correct inflow boundary condition is provided for the U-duct, a concurrent LES is performed of a straight duct with the same trapezoidal cross section and flow conditions as the U-duct. Results obtained for the U-duct show RANS models to be inadequate in predicting the separation due to their inability to predict the unsteady separation about the tip of the turn. To investigate the limitations of the RANS models, LES results were generated for the turbulent kinetic energy, Reynolds-stresses, pressure-strain rate, turbulent diffusion, pressure diffusion, turbulent transport, and velocity-temperature correlations with focus on understanding their behavior induced by the turn region of the U-duct. As expected, the Boussinesq assumption was found to be incorrect, which led to incorrect predictions of Reynolds stresses. For RSM-τω, the modeling of the pressure-strain rate was found to match LES data well, but huge error was found on modeling the turbulent diffusion. This huge error indicates that the two terms in the turbulent diffusion – pressure diffusion and turbulent transport – should be modeled separately. Since the turbulent transport was found to be ignorable, the focus should be on modeling the pressure diffusion. On the velocity-temperature correlations, the existing eddy-diffusivity model was found to be over simplified if there is unsteady separation with shedding. The generated LES data could be used to provide the guidance for a better model.<br></div>

A simplified computational fluid dynamics and heat transfer (CFD-HT) model of an electronics enclosure was developed. The compact model was based on a server simulator, which dissipates a variable amount of heat at an adjustable air flow rate. Even though a server simulator does not accurately represent the geometry of an actual electronics enclosure, the modeling of such a unit deals with many of the same issues as the modeling of actual enclosures. Even at the server simulator level, a disparity in length scales prevents detailed modeling of intricate components most notably grilles, fins, and fans. Therefore, a compact model for each of these components was developed. Fan performance curves were determined experimentally for varying fan rotational speeds. In addition, component pressure drop characteristics were found experimentally for grilles and fin banks, and these empirical relationships were applied to the model as well. To determine the validity of the simplifications employed in the model, experimental outlet temperature and velocity measurements were taken to compare to those provided by the CFD-HT simulations.

For the past four decades, Mars missions have relied on Viking heritage technology for supersonic descent. Extending the use of propulsion, which is required for Mars subsonic deceleration, into the supersonic regime allows the ability to land larger payload masses. Wind tunnel and computational experiments on subscale supersonic retropropulsion models have shown a complex aerodynamic flow field characterized by the interaction of underexpanded jet plumes exhausting from nozzles on the vehicle with the supersonic freestream. Understanding the impact of vehicle and nozzle configuration on this interaction is critical for analyzing the performance of a supersonic retropropulsion system, as deceleration will have components provided by both the aerodynamic drag of the vehicle and thrust from the nozzles. This investigation focuses on the validity of steady state computational approaches to analyze supersonic retropropulsion flowfield structures and their effect on vehicle aerodynamics. Wind tunnel data for a single nozzle and a multiple nozzle configuration are used to validate a steady state, turbulent computational fluid dynamics approach to modeling supersonic retropropulsion. An analytic approximation to determine plume and bow shock structure in the flow field is also developed, enabling rapid assessment of flowfield structure for use in improved grid generation and as a configuration screening tool. Results for both the computational fluid dynamics and analytic approaches show good agreement with the experimental datasets. Potential limitations of the two methods are identified based on the comparisons with available data. Six additional geometries are defined to investigate the extensibility of the analytical model and determine the variation of supersonic retropropulsion performance with configuration. These validation geometries are split into two categories: three geometries with nozzles located on the vehicle forebody at varying nozzle cant angles, and three geometries with nozzles located on the vehicle aftbody at varying nozzle cant angles and number of nozzles. The forebody nozzle configurations show that nozzle cant angle is a significant driver in performance of a vehicle employing supersonic retropropulsion. Aerodynamic drag preservation for a given thrust level increases with increasing cant angle. However, increasing the cant angle reduces the contribution of thrust to deceleration. The tradeoff between these two contributions to the deceleration force is examined, noting that performance improvements are possible with modest nozzle cant angles. Static pitch stability characteristics are investigated for the lowest and highest cant angle configurations. The aftbody nozzle configuration results show that removing the plume flow from the region forward of the vehicle results in less interaction with the bow shock structure. This impacts aerodynamic performance, as the surface pressure remains relatively undisturbed for all thrust values examined. Static pitch stability characteristics for each of the aftbody nozzle configurations are investigated; noting that supersonic retropropulsion for these configurations exhibits a transition point from static stability to instability as a function of this center of mass location along the axis.