Matching Items (482)
Filtering by
- All Subjects: Mechanical Engineering
- Member of: Theses and Dissertations
- Resource Type: Text
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 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
Description
Identification of early damage in polymer composite materials is of significant importance so that preventative measures can be taken before the materials reach catastrophic failure. Scientists have been developing damage detection technologies over many years and recently, mechanophore-based polymers, in which mechanical energy is translated to activate a chemical transformation, have received increasing attention. More specifically, the damage can be made detectable by mechanochromic polymers, which provide a visible color change upon the scission of covalent bonds under stress. This dissertation focuses on the study of a novel self-sensing framework for identifying early and in-situ damage by employing unique stress-sensing mechanophores. Two types of mechanophores, cyclobutane and cyclooctane, were utilized, and the former formed from cinnamoyl moeities and the latter formed from anthracene upon photodimerization. The effects on the thermal and mechanical properties with the addition of the cyclobutane-based polymers into epoxy matrices were investigated. The emergence of cracks was detected by fluorescent signals at a strain level right after the yield point of the polymer blends, and the fluorescence intensified with the accumulation of strain. Similar to the mechanism of fluorescence emission from the cleavage of cyclobutane, the cyclooctane moiety generated fluorescent emission with a higher quantum yield upon cleavage. The experimental results also demonstrated the success of employing the cyclooctane type mechanophore as a potential force sensor, as the fluorescence intensification was correlated with the strain increase.
ContributorsZou, Jin (Author) / Dai, Lenore L (Thesis advisor) / Chattopadhyay, Aditi (Thesis advisor) / Lind, Mary L (Committee member) / Mu, Bin (Committee member) / Yu, Hongyu (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
A level set approach for denoising and adaptively smoothing complex geometry stereolithography files
Description
Stereolithography files (STL) are widely used in diverse fields as a means of describing complex geometries through surface triangulations. The resulting stereolithography output is a result of either experimental measurements, or computer-aided design. Often times stereolithography outputs from experimental means are prone to noise, surface irregularities and holes in an otherwise closed surface.
A general method for denoising and adaptively smoothing these dirty stereolithography files is proposed. Unlike existing means, this approach aims to smoothen the dirty surface representation by utilizing the well established levelset method. The level of smoothing and denoising can be set depending on a per-requirement basis by means of input parameters. Once the surface representation is smoothened as desired, it can be extracted as a standard levelset scalar isosurface.
The approach presented in this thesis is also coupled to a fully unstructured Cartesian mesh generation library with built-in localized adaptive mesh refinement (AMR) capabilities, thereby ensuring lower computational cost while also providing sufficient resolution. Future work will focus on implementing tetrahedral cuts to the base hexahedral mesh structure in order to extract a fully unstructured hexahedra-dominant mesh describing the STL geometry, which can be used for fluid flow simulations.
A general method for denoising and adaptively smoothing these dirty stereolithography files is proposed. Unlike existing means, this approach aims to smoothen the dirty surface representation by utilizing the well established levelset method. The level of smoothing and denoising can be set depending on a per-requirement basis by means of input parameters. Once the surface representation is smoothened as desired, it can be extracted as a standard levelset scalar isosurface.
The approach presented in this thesis is also coupled to a fully unstructured Cartesian mesh generation library with built-in localized adaptive mesh refinement (AMR) capabilities, thereby ensuring lower computational cost while also providing sufficient resolution. Future work will focus on implementing tetrahedral cuts to the base hexahedral mesh structure in order to extract a fully unstructured hexahedra-dominant mesh describing the STL geometry, which can be used for fluid flow simulations.
ContributorsKannan, Karthik (Author) / Herrmann, Marcus (Thesis advisor) / Peet, Yulia (Committee member) / Frakes, David (Committee member) / Arizona State University (Publisher)
Created2014
Description
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.
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.
ContributorsPandit, Harshad Rajendra (Author) / Chen, Kang P (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2014
Description
Dimensional Metrology is the branch of science that determines length, angular, and geometric relationships within manufactured parts and compares them with required tolerances. The measurements can be made using either manual methods or sampled coordinate metrology (Coordinate measuring machines). Manual measurement methods have been in practice for a long time and are well accepted in the industry, but are slow for the present day manufacturing. On the other hand CMMs are relatively fast, but these methods are not well established yet. The major problem that needs to be addressed is the type of feature fitting algorithm used for evaluating tolerances. In a CMM the use of different feature fitting algorithms on a feature gives different values, and there is no standard that describes the type of feature fitting algorithm to be used for a specific tolerance. Our research is focused on identifying the feature fitting algorithm that is best used for each type of tolerance. Each algorithm is identified as the one to best represent the interpretation of geometric control as defined by the ASME Y14.5 standard and on the manual methods used for the measurement of a specific tolerance type. Using these algorithms normative procedures for CMMs are proposed for verifying tolerances. The proposed normative procedures are implemented as software. Then the procedures are verified by comparing the results from software with that of manual measurements.
To aid this research a library of feature fitting algorithms is developed in parallel. The library consists of least squares, Chebyshev and one sided fits applied on the features of line, plane, circle and cylinder. The proposed normative procedures are useful for evaluating tolerances in CMMs. The results evaluated will be in accordance to the standard. The ambiguity in choosing the algorithms is prevented. The software developed can be used in quality control for inspection purposes.
To aid this research a library of feature fitting algorithms is developed in parallel. The library consists of least squares, Chebyshev and one sided fits applied on the features of line, plane, circle and cylinder. The proposed normative procedures are useful for evaluating tolerances in CMMs. The results evaluated will be in accordance to the standard. The ambiguity in choosing the algorithms is prevented. The software developed can be used in quality control for inspection purposes.
ContributorsVemulapalli, Prabath (Author) / Shah, Jami J. (Thesis advisor) / Davidson, Joseph K. (Committee member) / Takahashi, Timothy (Committee member) / Arizona State University (Publisher)
Created2014
Description
The focus of this investigation is on the formulation and a validation of reduced order models (ROMs) for the prediction of the response of structures with embedded piezoelectric actuators. The ROMs considered here are those constructed in a nonintrusive manner from a commercial finite element software, NASTRAN is adopted here. Notwithstanding the popularity of piezoelectric materials in structural dynamics related applications such as structural health monitoring and energy harvesting, not all commercial finite element software allow directly their modeling. In such cases, e.g., with NASTRAN, one can proceed with an analogy and replace the electric actuation in the piezoelectric material by a fictitious thermal effect producing the same strain. This process recasts the determination of a ROM for a structure with embedded piezoelectric actuator into a similar ROM but for a heated structure, the framework of which has recently been developed. Yet, the temperature field resulting from the analogy would be quite different from the one considered in past effort and would excite a broad array of structural modes. Accordingly, as a preamble to considering a beam with a piezoelectric layer, a simpler plate model is considered that is subjected to a uniform temperature but a complex pressure loading that excites the entire set of modes of the plate in the broad frequency band considered. The very good match of the predictions obtained by this ROM in comparison to their full finite element counterparts provides the necessary confidence to next address a beam with embedded piezoelectic actuator. The test model considered for this validation is a built-up nano beam analyzed recently in nonlinear geometric conditions by full finite elements and by a non-intrusive ROM procedure under harmonic variations of the piezoelectic voltage. This structural model and its loading conditions are very different from those considered in past applications of nonintrusive ROMs, thus the excellent results obtained here provide further support of the broad generality of the nonintrusive ROM methodology, including of the appropriateness of the "dual modes" basis functions.
ContributorsVyas, Varun (Author) / Mignolet, Marc (Thesis advisor) / Hollkamp, Joseph (Committee member) / Liu, Yongming (Committee member) / Arizona State University (Publisher)
Created2014
Description
With the ever-increasing demand for high-end services, technological companies have been forced to operate on high performance servers. In addition to the customer services, the company's internal need to store and manage huge amounts of data has also increased their need to invest in High Density Data Centers. As a result, the performance to size of the data center has increased tremendously. Most of the consumed power by the servers is emitted as heat. In a High Density Data Center, the power per floor space area is higher compared to the regular data center. Hence the thermal management of this type of data center is relatively complicated.
Because of the very high power emission in a smaller containment, improper maintenance can result in failure of the data center operation in a shorter period. Hence the response time of the cooler to the temperature rise of the servers is very critical. Any delay in response will constantly lead to increased temperature and hence the server's failure.
In this paper, the significance of this delay time is understood by performing CFD simulation on different variants of High Density Modules using ANSYS Fluent. It was found out that the delay was becoming longer as the size of the data center increases. But the overload temperature, ie. the temperature rise beyond the set-point became lower with the increase in data center size. The results were common for both the single-row and the double-row model. The causes of the increased delay are accounted and explained in detail manner in this paper.
Because of the very high power emission in a smaller containment, improper maintenance can result in failure of the data center operation in a shorter period. Hence the response time of the cooler to the temperature rise of the servers is very critical. Any delay in response will constantly lead to increased temperature and hence the server's failure.
In this paper, the significance of this delay time is understood by performing CFD simulation on different variants of High Density Modules using ANSYS Fluent. It was found out that the delay was becoming longer as the size of the data center increases. But the overload temperature, ie. the temperature rise beyond the set-point became lower with the increase in data center size. The results were common for both the single-row and the double-row model. The causes of the increased delay are accounted and explained in detail manner in this paper.
ContributorsRamaraj, Dinesh Balaji (Author) / Gupta, Sandeep (Thesis advisor) / Hermann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2015
Description
Gallium-based liquid metals are of interest for a variety of applications including flexible electronics, soft robotics, and biomedical devices. Still, nano- to microscale device fabrication with these materials is challenging because of their strong adhesion to a majority of substrates. This unusual high adhesion is attributed to the formation of a thin oxide shell; however, its role in the adhesion process has not yet been established. In the first part of the thesis, we described a multiscale study aiming at understanding the fundamental mechanisms governing wetting and adhesion of gallium-based liquid metals. In particular, macroscale dynamic contact angle measurements were coupled with Scanning Electron Microscope (SEM) imaging to relate macroscopic drop adhesion to morphology of the liquid metal-surface interface. In addition, room temperature liquid-metal microfluidic devices are also attractive systems for hyperelastic strain sensing. Currently two types of liquid metal-based strain sensors exist for inplane measurements: single-microchannel resistive and two-microchannel capacitive devices. However, with a winding serpentine channel geometry, these sensors typically have a footprint of about a square centimeter, limiting the number of sensors that can be embedded into. In the second part of the thesis, firstly, simulations and an experimental setup consisting of two GaInSn filled tubes submerged within a dielectric liquid bath are used to quantify the effects of the cylindrical electrode geometry including diameter, spacing, and meniscus shape as well as dielectric constant of the insulating liquid and the presence of tubing on the overall system's capacitance. Furthermore, a procedure for fabricating the two-liquid capacitor within a single straight polydiemethylsiloxane channel is developed. Lastly, capacitance and response of this compact device to strain and operational issues arising from complex hydrodynamics near liquid-liquid and liquid-elastomer interfaces are described.
ContributorsLiu, Shanliangzi (Author) / Rykaczewski, Konrad (Thesis advisor) / Alford, Terry (Committee member) / Herrmann, Marcus (Committee member) / Hildreth, Owen (Committee member) / Arizona State University (Publisher)
Created2015