The gasification of carbonaceous fuel materials offers significant potential for the production of both energy and chemical products. Advancement of gasification technologies may be expedited through the use of computational fluid dynamics, as virtual reactor design offers a low cost method for system prototyping. To that end, a series of numerical studies were conducted to identify a computational modeling strategy for the simulation of coal gasification in fluidized bed reactors. The efforts set forth by this work first involved the development of a validatable hydrodynamic modeling strategy for the simulation of sand and coal fluidization. Those fluidization models were then applied to systems at elevated temperatures and polydisperse systems that featured a complex material injection geometry, for which no experimental data exists. A method for establishing similitude between 2-D and 3-D multiphase systems that feature non-symmetric material injection were then delineated and numerically tested. Following the development of the hydrodynamic modeling strategy, simulations of coal gasification were conducted using three different chemistry models. Simulated results were compared to experimental outcomes in an effort to assess the validity of each gasification chemistry model. The chemistry model that exhibited the highest degree of agreement with the experimental findings was then further analyzed identify areas of potential improvement. / Ph. D. / Efficient utilization of coal is critical to ensuring stable domestic energy supplies while mitigating human impact on climate change. This idea may be realized through the use of gasification systems technologies. The design and planning of next-generation coal gasification reactors can benefit from the use of computational simulations to reduce both development time and cost. This treatise presents several studies where computational fluid dynamics was applied to the problem of coal gasification in a bubbling fluidized bed reactor with focuses on accurate tracking of solid material locations and modeling of chemical reactions.

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 research examines different techniques and proposes a Computational Fluid Dynamics (CFD) model as a robust tool for flow characterization of load leveling valves. The load leveling valve is a critical component of an air suspension system since it manages air spring pressure, a key function that directly impacts vehicle dynamic performance in addition to maintaining a static ride height. Efficiency of operation of a load leveling valve is established by its flow characteristics, a metric useful in determining suitability of the valve for application in a truck-suspension configuration and for comparison among similar products. The disk-slot type load leveling valve was chosen as the subject of this study due to its popularity in the heavy truck industry. Three distinct methods are presented to model and evaluate flow characteristics of a disk-slot valve. First is a theoretical formulation based on gas dynamic behavior through an orifice; second is an experimental technique in which a full pneumatic apparatus is used to collect instantaneous pressure data to estimate air discharge; and third is a CFD approach. Significant discrepancies observed between theoretically estimated results and experimental data suggest that the theoretical model is incapable of accurately capturing losses that occur during air flow. These variations diminish as the magnitude of discharge coefficient is altered. A detailed CFD model is submitted as an effective tool for load leveling valve flow characterization/analysis. This model overcomes the deficiencies of the theoretical model and improves the accuracy of simulations. A 2-D axisymmetric approximation of the real fluid domain is analyzed for flow characteristics using a Realizable k-ϵ turbulence model, scalable wall functions, and a pressure-based coupled algorithm with a second order discretization function. The CFD-generated results were observed to be in agreement with the experimental findings. CFD is found to be advantageous in the evaluation of flow characteristics as it furnishes precise data without the need to experimentally evaluate a physical model/prototype of the valve, thereby benefitting suspension engineers involved in the development and testing of load leveling valve designs. This document concludes with a sample case study which uses CFD to characterize flow in a modified disk-slot load leveling valve, and discusses the results in light of application on a heavy truck. / MS / A majority of heavy trucks in North America equipped with air suspensions use a device known as a load leveling valve. This is a mechanical control system which manages pressure in air springs to maintain a preset/constant static ride height irrespective of the payload, doing so by sensing the distance between the truck frame and the axle. The rate of airflow to/from air springs in response to a road disturbance or load shift is critical to the stability of the truck when on the road. This rate of airflow for a given set of conditions constitutes flow characteristics of a load leveling valve. Accurate measurement of flow characteristics is necessary to understand the actual effect of the use of a particular valve on a truck-suspension configuration. This research addresses that requirement by presenting three distinct methods to model and evaluate flow characteristics of a load leveling valve, conducted on the disk-slot valve for its popularity in the heavy truck industry. First is a theoretical formulation based on flow of gas through an orifice; second is an experimental technique in which a full pneumatic apparatus is used to collect instantaneous pressure data to estimate air discharge; and third is a Computational Fluid Dynamics (CFD) approach. Significant discrepancies observed between theoretically estimated results and experimental data suggest that the theoretical model is incapable of accurately capturing losses that occur during air flow. The disparities also justify the adoption of CFD as an alternate method. A comprehensive CFD model is proposed as a capable tool for load leveling valve flow analysis/characterization. This model overcomes the deficiencies of the theoretical model and improves the accuracy of simulations. CFD-generated results are found to be in agreement with the experimental findings, highlighting its effectiveness at flow characterization. The ability of a CFD model to furnish precise data without the need to experimentally evaluate a physical model/prototype of the valve promises to benefit suspension engineers involved in the development and testing of load leveling valve designs. This document concludes with a sample case study which uses CFD to characterize flow in a modified disk-slot valve, and discusses the results in light of application on a heavy truck.

The rapid increase in computational power over the last decade has unlocked the possibility of providing patient-specific healthcare via simulation and data assimilation. In the past 2 decades, computational approaches to simulating cardiovascular flows have advanced significantly due to intense research and adoption of methods in medical device companies. A significant source of friction in porting these tools to the hospital and getting in the hands of surgeons is due to the expertise required to handle the geometry pre-processing and meshing of models. Meshless meth- ods reduce the amount of corner cases which makes it easier to develop robust tools surgeons need. To accurately simulate modifications to a region of vasculature as in surgical planning, the entire system must be modeled. Unfortunately, this is computationally prohibitive even on to- day’s machines. To circumvent this issue, the Radial-Basis Function Finite Difference (RBF-FD) method for solution of the higher-dimensional (2D/3D) region of interest is tightly-coupled to a 0D Lumped-Parameter Model (LPM) for solution of the peripheral circulation. The incompress- ible flow equations are updated by an explicit time-marching scheme based on a pressure-velocity correction algorithm. The inlets and outlets of the domain are tightly coupled with the LPM which contains elements that draw from a fluid-electrical analogy such as resistors, capacitors, and in- ductors that represent the viscous resistance, vessel compliance, and flow inertia, respectively. The localized RBF meshless approach is well-suited for modeling complicated non-Newtonian hemo- dynamics due to ease of spatial discretization, ease of addition of multi-physics interactions such as fluid-structure interaction of the vessel wall, and ease of parallelization for fast computations. This work introduces the tight coupling of meshless methods and LPMs for fast and accurate hemody- namic simulations. The results show the efficacy of the method to be used in building robust tools to inform surgical decisions and further development is motivated.

The performance of a vapor compression system is known to be affected by the ingestion of liquid droplets in the compressor. In these multiphase flows, the liquid and the vapor phase are tightly coupled. Therefore the interphase heat, mass and momentum transfer as well as droplet dynamics including droplet breakup and droplet-wall interactions play a vital role in governing these flows. Only thermodynamic analyses or two-dimensional mean-line calculations are not sufficient to gain an in-depth understanding of the complex multiphase flow field within the compressor. The objective of this research was to extend the current understanding of the operation of a multistage centrifugal compressor under two-phase flow conditions, by performing three-dimensional computational analysis. In this work, two-phase flow of a single constituent (refrigerant R134a) through a two-stage, in-line centrifugal compressor was analyzed using CFD. The CFD model accounted for real gas behavior of the vapor phase. Novel user defined routines were implemented to ensure accurate calculations of interphase heat, mass and momentum transfer terms and to model droplet impact on the compressor surfaces. An erosion model was developed and implemented to locate the erosion "hot spots" and to estimate the amount of material eroded. To understand the effects of increasing liquid carryover, the mass flow rate of the liquid phase was increased from 1% to 5% of the vapor mass flow rate. The influence of droplet size on the compressor performance was assessed by varying the droplet diameter at the inlet from 100 microns to 400 microns. The results of the two-phase flow simulations were compared with the simulation involving only the vapor phase. Liquid carryover altered the flow field within the compressor, and as a result, both impellers were observed to operate at off-design conditions. This effect was more pronounced for the second impeller. The overall effects of liquid carryover were detrimental to the compressor performance. The erosion calculations showed maximum erosion potential on the blade and shroud of the first impeller. The results from this investigation provided new and useful information that can be used to support improved design solutions. / Ph. D. / The performance of a compressor is known to be affected by the ingestion of liquid droplets, and thus, it is a research topic of interest for both academia as well as industry. This work extends the current understanding of the operation of a multistage centrifugal compressor under two-phase flow conditions, by employing high-fidelity computational fluid dynamics (CFD). In this research, the two-phase flow of refrigerant R134a through a two-stage, in-line centrifugal compressor was analyzed. The CFD model used in this research incorporated real gas behavior of the vapor phase, as well as the interphase heat, mass and momentum transfer processes. An erosion model was also developed and implemented to locate the erosion "hot spots" on the compressor surfaces, and to estimate the amount of material eroded. The effects of increasing the liquid carryover, as well as the influence of droplet size on the compressor performance were assessed. Liquid carryover altered the flow field within the compressor. As a result, the compressor operated at off-design conditions. The overall effects of liquid carryover were detrimental to the compressor performance. The erosion calculations showed maximum erosion potential on the blade and shroud of the first impeller. The results from this investigation provided new and useful information that can be used to support improved design solutions.

The labyrinth seal is widely used in turbo machines to reduce leakage flow. The stability of the rotor is influenced by the labyrinth seal because of the driving forces generated in the seal. The working fluid usually has a circumferential velocity component before entering the seal; the ratio of circumferential velocity and shaft synchronous surface velocity is defined as pre-swirl rate. It has been observed that pre-swirl rate is an important factor affecting driving forces in the labyrinth seal thus affecting the stability of the rotor. Besides the pre-swirl, the eccentricity, the clearance, and the configuration of tooth locations are all factors affecting the rotordynamic properties of the labyrinth seal. So it is of interest to investigate the exact relationships between those factors and the seal's rotordynamic properties. In this research, three types of labyrinth seals have been modeled: the straight eye seal, the stepped eye seal, and the balance drum seal. For the straight eye seal, a series of models were built to study the influence of eccentricity and clearance. The other two seals each have only one model. All models were built with Solid Works and meshed with ANSYS-ICEM. Flows in those models were simulated by numerically solving the Reynolds-Averaged Navier-Stokes (RANS) equations in the ANSYS-CFX and then rotordynamic coefficients for each seal were calculated based on the numerical results. It had previously been very difficult to generate a pre-swirl rate higher than 60% in a numerical simulation. So three ways to create pre-swirl in ANSYS-CFX were studied and finally the method by specifying the inlet velocity ratio was employed. Numerical methods used in this research were introduced including the frame transfer, the k-ε turbulence model with curvature correction, and the scalable wall function. To obtain the optimal mesh and minimize the discretization error, a systematical grid study was conducted including grid independence studies and discretization error estimations. Some of the results were compared with previous bulk-flow or experimental results to validate the numerical model and method. The fluid field in the labyrinth seal must be analyzed before conducting rotordynamic analysis. The predicted pressure distributions and leakages were compared with bulk-flow results. A second small vortex at the downstream edge of each tooth was found in the straight eye seal. This has never been reported before and the discovery of this small vortex will help to improve seal designs in the future. The detailed flows in discharged region and in chambers were also discussed. Radial and tangential forces on the rotor were solved based on the fluid field results. It is shown that the traditional first-order rotordynamic model works well for low pre-swirl cases but does not accurately reflect the characteristics for high pre-swirl cases. For example compressor eye seals usually have pre-swirl rates bigger than 70% and the second order model is required. Thus a second-order model including inertia terms was built and applied to the rotordynamic analysis in this research. The influence of pre-swirl, eccentricity and clearance were studied using the straight eye seal model. The rotordynamic characteristics of the stepped eye seal and the balance drum seal were studied considering high pre-swirl rates. Some relationships between influencing factors and the four rotordynamic coefficients were concluded. The results also showed that for all the three seals higher pre-swirl leads to higher cross-coupled stiffness which is one of the main factors causing rotor instability. The rotor stability analysis was conducted to study the influence of drum balance seal on the stability. The rotor was designed with typical dimensions and natural frequencies for a centrifugal compressor rotor. The parameters for bearing and aerodynamic force were also set according to general case in compressors to minimize the effects from them. The result shows that the high pre-swirl rate in balance drum seal leads to rotor instability, which confirmed the significant effect of pre-swirl on the seal and the rotor system. / Ph. D.

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.