Matching Items (287)
Filtering by
- All Subjects: Mechanical Engineering
- Genre: Masters Thesis
Description
Cavitation erosion is a significant cause of wear in marine components, such as impellers, propellers or rudders. While the erosion process has been widely studied on metals, the effect of cavitation on polymers is not well-understood. The stress response in metals differs greatly from that of polymers, e.g. rate and temperature effects are far more important, thus damage and wear mechanisms of polymers under cavitating flows are significantly different. In this work, heat-driven failure caused by viscous dissipation and void nucleation resulting from tensile stresses arising from stress wave reflections are investigated as two possible material failure mechanisms.
As a first step in developing a fundamental understanding of the cavitation erosion process on polymer surfaces, simulations are performed of the collapse of individual bubbles against a compliant surface e.g. metallic substrates with polyurea coatings. The surface response of collapse-driven impact loads is represented by a idealized, time-dependent, Gaussian pressure distribution on the surface. A two-dimensional distribution of load radii and durations is considered corresponding to characteristic of cavitating flows accelerated erosion experiments. Finite element simulations are performed to fit a response curve that relates the loading parameters to the energy dissipated in the coating and integrated with collapse statistics to generate an expected heat input into the coating.
The impulsive pressure, which is generated due to bubble collapse, impacts the material and generates intense shock waves. The stress waves within the material reflects by interaction with the substrate. A transient region of high tensile stress is produced by the interaction of these waves. Simulations suggests that maximum hydrostatic tension which cause failure of polyurea layer is observed in thick coating. Also, the dissipated viscous energy and corresponding temperature rise in a polyurea is calculated, and it is concluded that temperature has influence on deformation.
As a first step in developing a fundamental understanding of the cavitation erosion process on polymer surfaces, simulations are performed of the collapse of individual bubbles against a compliant surface e.g. metallic substrates with polyurea coatings. The surface response of collapse-driven impact loads is represented by a idealized, time-dependent, Gaussian pressure distribution on the surface. A two-dimensional distribution of load radii and durations is considered corresponding to characteristic of cavitating flows accelerated erosion experiments. Finite element simulations are performed to fit a response curve that relates the loading parameters to the energy dissipated in the coating and integrated with collapse statistics to generate an expected heat input into the coating.
The impulsive pressure, which is generated due to bubble collapse, impacts the material and generates intense shock waves. The stress waves within the material reflects by interaction with the substrate. A transient region of high tensile stress is produced by the interaction of these waves. Simulations suggests that maximum hydrostatic tension which cause failure of polyurea layer is observed in thick coating. Also, the dissipated viscous energy and corresponding temperature rise in a polyurea is calculated, and it is concluded that temperature has influence on deformation.
ContributorsPanwar, Ajay (Author) / Oswald, Jay (Thesis advisor) / Dooley, Kevin (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
Description
Recent studies of the occurrence of post-flutter limit cycle oscillations (LCO) of the F-16 have provided good support to the long-standing hypothesis that this phenomenon involves a nonlinear structural damping. A potential mechanism for the appearance of nonlinearity in the damping are the nonlinear geometric effects that arise when the deformations become large enough to exceed the linear regime. In this light, the focus of this investigation is first on extending nonlinear reduced order modeling (ROM) methods to include viscoelasticity which is introduced here through a linear Kelvin-Voigt model in the undeformed configuration. Proceeding with a Galerkin approach, the ROM governing equations of motion are obtained and are found to be of a generalized van der Pol-Duffing form with parameters depending on the structure and the chosen basis functions. An identification approach of the nonlinear damping parameters is next proposed which is applicable to structures modeled within commercial finite element software.
The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses.
The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses.
ContributorsSong, Pengchao (Author) / Mignolet, Marc P (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2015
Description
As the robotic industry becomes increasingly present in some of the more extreme environments such as the battle field, disaster sites or extraplanetary exploration, it will be necessary to provide locomotive niche strategies that are optimal to each terrain. The hopping gait has been well studied in robotics and proven to be a potential method to fit some of these niche areas. There have been some difficulties in producing terrain following controllers that maintain robust, steady state, which are disturbance resistant.
The following thesis will discuss a controller which has shown the ability to produce these desired properties. A phase angle oscillator controller is shown to work remarkably well, both in simulation and with a one degree of freedom robotic test stand.
Work was also done with an experimental quadruped with less successful results, but which did show potential for stability. Additional work is suggested for the quadruped.
The following thesis will discuss a controller which has shown the ability to produce these desired properties. A phase angle oscillator controller is shown to work remarkably well, both in simulation and with a one degree of freedom robotic test stand.
Work was also done with an experimental quadruped with less successful results, but which did show potential for stability. Additional work is suggested for the quadruped.
ContributorsNew, Philip Wesley (Author) / Sugar, Thomas G. (Thesis advisor) / Artemiadis, Panagiotis (Committee member) / Redkar, Sangram (Committee member) / Arizona State University (Publisher)
Created2015
Description
The football helmet is a device used to help mitigate the occurrence of impact-related traumatic (TBI) and minor traumatic brain injuries (mTBI) in the game of American football. The current design methodology of using a hard shell with an energy absorbing liner may be adequate for minimizing TBI, however it has had less effect in minimizing mTBI. The latest research in brain injury mechanisms has established that the current design methodology has produced a helmet to reduce linear acceleration of the head. However, angular accelerations also have an adverse effect on the brain response, and must be investigated as a contributor of brain injury.
To help better understand how the football helmet design features effect the brain response during impact, this research develops a validated football helmet model and couples it with a full LS-DYNA human body model developed by the Global Human Body Modeling Consortium (v4.1.1). The human body model is a conglomeration of several validated models of different sections of the body. Of particular interest for this research is the Wayne State University Head Injury Model for modeling the brain. These human body models were validated using a combination of cadaveric and animal studies. In this study, the football helmet was validated by laboratory testing using drop tests on the crown of the helmet. By coupling the two models into one finite element model, the brain response to impact loads caused by helmet design features can be investigated. In the present research, LS-DYNA is used to study a helmet crown impact with a rigid steel plate so as to obtain the strain-rate, strain, and stress experienced in the corpus callosum, midbrain, and brain stem as these anatomical regions are areas of concern with respect to mTBI.
To help better understand how the football helmet design features effect the brain response during impact, this research develops a validated football helmet model and couples it with a full LS-DYNA human body model developed by the Global Human Body Modeling Consortium (v4.1.1). The human body model is a conglomeration of several validated models of different sections of the body. Of particular interest for this research is the Wayne State University Head Injury Model for modeling the brain. These human body models were validated using a combination of cadaveric and animal studies. In this study, the football helmet was validated by laboratory testing using drop tests on the crown of the helmet. By coupling the two models into one finite element model, the brain response to impact loads caused by helmet design features can be investigated. In the present research, LS-DYNA is used to study a helmet crown impact with a rigid steel plate so as to obtain the strain-rate, strain, and stress experienced in the corpus callosum, midbrain, and brain stem as these anatomical regions are areas of concern with respect to mTBI.
ContributorsDarling, Timothy (Author) / Rajan, Subramaniam D. (Thesis advisor) / Muthuswamy, Jitendran (Thesis advisor) / Oswald, Jay (Committee member) / Mignolet, Marc (Committee member) / Arizona State University (Publisher)
Created2014
Description
Monte Carlo simulations are traditionally carried out for the determination of the amplification of forced vibration response of turbomachine/jet engine blades to mistuning. However, this effort can be computationally time consuming even when using the various reduced order modeling techniques. Accordingly, some investigations in the past have focused on obtaining simple approximate estimates for this amplification. In particular, two of these have proposed the use of harmonic patterns of the blade properties around the disk as an approximate alternative to the many random patterns of Monte Carlo analyses. These investigations, while quite encouraging, have relied solely on single degree of freedom per sector models of the rotor.
In this light, the overall focus of the present effort is a revisit of harmonic
mistuning of rotors focusing first the confirmation of the previously obtained findings with a more detailed model of the blisk in both conditions of an isolated blade-dominated resonance and of a veering between blade and disk dominated modes. The latter condition cannot be simulated by a single degree of freedom per sector model. Further, the analysis will consider the distinct cases of mistuning due to variations of material properties (Young's modulus) and geometric properties (geometric mistuning). In the single degree of freedom model, both mistuning types are equivalent but they are not, as demonstrated here, in more realistic models. The difference arises because changes in geometry induce not only changes in natural frequencies of the blades alone but of their modes and the importance of these two sources of variability is discussed with both Monte Carlo simulation and harmonic mistuning results.
The present investigation focuses also on the possible extension of the harmonic mistuning concept and of its quantitative information that can be derived from such analyses. From it, a novel measure of blade-disk coupling is introduced and assessed in comparison with the coupling index introduced in the past. In conclusions, the low cost of harmonic mistuning computations in comparison with full Monte Carlo simulations is
demonstrated to be worthwhile to elucidate the basic behavior of the mistuned rotor in a random setting.
In this light, the overall focus of the present effort is a revisit of harmonic
mistuning of rotors focusing first the confirmation of the previously obtained findings with a more detailed model of the blisk in both conditions of an isolated blade-dominated resonance and of a veering between blade and disk dominated modes. The latter condition cannot be simulated by a single degree of freedom per sector model. Further, the analysis will consider the distinct cases of mistuning due to variations of material properties (Young's modulus) and geometric properties (geometric mistuning). In the single degree of freedom model, both mistuning types are equivalent but they are not, as demonstrated here, in more realistic models. The difference arises because changes in geometry induce not only changes in natural frequencies of the blades alone but of their modes and the importance of these two sources of variability is discussed with both Monte Carlo simulation and harmonic mistuning results.
The present investigation focuses also on the possible extension of the harmonic mistuning concept and of its quantitative information that can be derived from such analyses. From it, a novel measure of blade-disk coupling is introduced and assessed in comparison with the coupling index introduced in the past. In conclusions, the low cost of harmonic mistuning computations in comparison with full Monte Carlo simulations is
demonstrated to be worthwhile to elucidate the basic behavior of the mistuned rotor in a random setting.
ContributorsSahoo, Saurav (Author) / Mignolet, Marc Paul (Thesis advisor) / Chattopadhyay, Aditi (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2014
Description
The Volume-of-Fluid method is a popular method for interface tracking in Multiphase applications within Computational Fluid Dynamics. To date there exists several algorithms for reconstruction of a geometric interface surface. Of these are the Finite Difference algorithm, Least Squares Volume-of-Fluid Interface Reconstruction Algorithm, LVIRA, and the Efficient Least Squares Volume-of-Fluid Interface Reconstruction Algorithm, ELVIRA. Along with these geometric interface reconstruction algorithms, there exist several volume-of-fluid transportation algorithms. This paper will discuss two operator-splitting advection algorithms and an unsplit advection algorithm. Using these three interface reconstruction algorithms, and three advection algorithms, a comparison will be drawn to see how different combinations of these algorithms perform with respect to accuracy as well as computational expense.
ContributorsKedelty, Dominic (Author) / Herrmann, Marcus (Thesis advisor) / Huang, Huei-Ping (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2015
Description
In this work, the hydrodynamics of Suction Stabilization is studied. Suction stabilization was found to stabilize floating platforms/floats in a much better way as compared to the conventional methods. This was achieved by an effective increment in the metacentric height due to the Inverse Slack Tank (IST) effect. The study involves the analysis of the existing designs and optimizing its performance. This research investigates the stability of such floats and the hydrodynamic forces acting on the same for offshore applications, such as wind turbines. A simple mathematical model for the condition of parametric resonance is developed and the results are verified, both analytically and experimentally.
ContributorsCherangara Subramanian, Susheelkumar (Author) / Redkar, Sangram (Thesis advisor) / Rajadas, John (Committee member) / Rogers, Bradley (Committee member) / Arizona State University (Publisher)
Created2014
Description
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 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.
ContributorsSowers, Timothy Wayne (Author) / Rajagopalan, Jagannathan (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2014
Description
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 (cw,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 (CO2) 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 cw,min is also proposed.
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 (CO2) 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 cw,min is also proposed.
ContributorsSingh, Prashant (Author) / Roy, Ramendra P (Thesis advisor) / Mignolet, Marc (Thesis advisor) / Lee, Taewoo (Committee member) / Arizona State University (Publisher)
Created2014
Optimization of the implementation of renewable resources in a municipal electric utility in Arizona
Description
A municipal electric utility in Mesa, Arizona with a peak load of approximately 85 megawatts (MW) was analyzed to determine how the implementation of renewable resources (both wind and solar) would affect the overall cost of energy purchased by the utility. The utility currently purchases all of its energy through long term energy supply contracts and does not own any generation assets and so optimization was achieved by minimizing the overall cost of energy while adhering to specific constraints on how much energy the utility could purchase from the short term energy market. Scenarios were analyzed for a five percent and a ten percent penetration of renewable energy in the years 2015 and 2025. Demand Side Management measures (through thermal storage in the City's district cooling system, electric vehicles, and customers' air conditioning improvements) were evaluated to determine if they would mitigate some of the cost increases that resulted from the addition of renewable resources.
In the 2015 simulation, wind energy was less expensive than solar to integrate to the supply mix. When five percent of the utility's energy requirements in 2015 are met by wind, this caused a 3.59% increase in the overall cost of energy. When that five percent is met by solar in 2015, it is estimated to cause a 3.62% increase in the overall cost of energy. A mix of wind and solar in 2015 caused a lower increase in the overall cost of energy of 3.57%. At the ten percent implementation level in 2015, solar, wind, and a mix of solar and wind caused increases of 7.28%, 7.51% and 7.27% respectively in the overall cost of energy.
In 2025, at the five percent implementation level, wind and solar caused increases in the overall cost of energy of 3.07% and 2.22% respectively. In 2025, at the ten percent implementation level, wind and solar caused increases in the overall cost of energy of 6.23% and 4.67% respectively.
Demand Side Management reduced the overall cost of energy by approximately 0.6%, mitigating some of the cost increase from adding renewable resources.
In the 2015 simulation, wind energy was less expensive than solar to integrate to the supply mix. When five percent of the utility's energy requirements in 2015 are met by wind, this caused a 3.59% increase in the overall cost of energy. When that five percent is met by solar in 2015, it is estimated to cause a 3.62% increase in the overall cost of energy. A mix of wind and solar in 2015 caused a lower increase in the overall cost of energy of 3.57%. At the ten percent implementation level in 2015, solar, wind, and a mix of solar and wind caused increases of 7.28%, 7.51% and 7.27% respectively in the overall cost of energy.
In 2025, at the five percent implementation level, wind and solar caused increases in the overall cost of energy of 3.07% and 2.22% respectively. In 2025, at the ten percent implementation level, wind and solar caused increases in the overall cost of energy of 6.23% and 4.67% respectively.
Demand Side Management reduced the overall cost of energy by approximately 0.6%, mitigating some of the cost increase from adding renewable resources.
ContributorsCadorin, Anthony (Author) / Phelan, Patrick (Thesis advisor) / Calhoun, Ronald (Committee member) / Trimble, Steve (Committee member) / Arizona State University (Publisher)
Created2014