Matching Items (102)
Filtering by
- All Subjects: Mechanical Engineering
- Creators: Herrmann, Marcus
- Creators: Lee, Hyunglae
- Creators: Mignolet, Marc
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
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
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
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
Description
A new method of adaptive mesh generation for the computation of fluid flows is investigated. The method utilizes gradients of the flow solution to adapt the size and stretching of elements or volumes in the computational mesh as is commonly done in the conventional Hessian approach. However, in the new method, higher-order gradients are used in place of the Hessian. The method is applied to the finite element solution of the incompressible Navier-Stokes equations on model problems. Results indicate that a significant efficiency benefit is realized.
ContributorsShortridge, Randall (Author) / Chen, Kang Ping (Thesis advisor) / Herrmann, Marcus (Thesis advisor) / Wells, Valana (Committee member) / Huang, Huei-Ping (Committee member) / Mittelmann, Hans (Committee member) / Arizona State University (Publisher)
Created2011
Description
The flow around a golf ball is studied using direct numerical simulation (DNS). An immersed boundary approach is adopted in which the incompressible Navier-Stokes equations are solved using a fractional step method on a structured, staggered grid in cylindrical coordinates. The boundary conditions on the surface are imposed using momentum forcing in the vicinity of the boundary. The flow solver is parallelized using a domain decomposition strategy and message passing interface (MPI), and exhibits linear scaling on as many as 500 processors. A laminar flow case is presented to verify the formal accuracy of the method. The immersed boundary approach is validated by comparison with computations of the flow over a smooth sphere. Simulations are performed at Reynolds numbers of 2.5 × 104 and 1.1 × 105 based on the diameter of the ball and the freestream speed and using grids comprised of more than 1.14 × 109 points. Flow visualizations reveal the location of separation, as well as the delay of complete detachment. Predictions of the aerodynamic forces at both Reynolds numbers are in reasonable agreement with measurements. Energy spectra of the velocity quantify the dominant frequencies of the flow near separation and in the wake. Time-averaged statistics reveal characteristic physical patterns in the flow as well as local trends within dimples. A mechanism of drag reduction due to the dimples is confirmed, and metrics for dimple optimization are proposed.
ContributorsSmith, Clinton E (Author) / Squires, Kyle D (Thesis advisor) / Balaras, Elias (Committee member) / Herrmann, Marcus (Committee member) / Adrian, Ronald (Committee member) / Stanzione, Daniel C (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2011
Description
A numerical study of incremental spin-up and spin-up from rest of a thermally- stratified fluid enclosed within a right circular cylinder with rigid bottom and side walls and stress-free upper surface is presented. Thermally stratified spin-up is a typical example of baroclinity, which is initiated by a sudden increase in rotation rate and the tilting of isotherms gives rise to baroclinic source of vorticity. Research by (Smirnov et al. [2010a]) showed the differences in evolution of instabilities when Dirichlet and Neumann thermal boundary conditions were applied at top and bottom walls. Study of parametric variations carried out in this dissertation confirmed the instability patterns observed by them for given aspect ratio and Rossby number values greater than 0.5. Also results reveal that flow maintained axisymmetry and stability for short aspect ratio containers independent of amount of rotational increment imparted. Investigation on vorticity components provides framework for baroclinic vorticity feedback mechanism which plays important role in delayed rise of instabilities when Dirichlet thermal Boundary Conditions are applied.
ContributorsKher, Aditya Deepak (Author) / Chen, Kangping (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2011
Description
Multiphase flows are an important part of many natural and technological phe- nomena such as ocean-air coupling (which is important for climate modeling) and the atomization of liquid fuel jets in combustion engines. The unique challenges of multiphase flow often make analytical solutions to the governing equations impos- sible and experimental investigations very difficult. Thus, high-fidelity numerical simulations can play a pivotal role in understanding these systems. This disserta- tion describes numerical methods developed for complex multiphase flows and the simulations performed using these methods. First, the issue of multiphase code verification is addressed. Code verification answers the question "Is this code solving the equations correctly?" The method of manufactured solutions (MMS) is a procedure for generating exact benchmark solutions which can test the most general capabilities of a code. The chief obstacle to applying MMS to multiphase flow lies in the discontinuous nature of the material properties at the interface. An extension of the MMS procedure to multiphase flow is presented, using an adaptive marching tetrahedron style algorithm to compute the source terms near the interface. Guidelines for the use of the MMS to help locate coding mistakes are also detailed. Three multiphase systems are then investigated: (1) the thermocapillary motion of three-dimensional and axisymmetric drops in a confined apparatus, (2) the flow of two immiscible fluids completely filling an enclosed cylinder and driven by the rotation of the bottom endwall, and (3) the atomization of a single drop subjected to a high shear turbulent flow. The systems are simulated numerically by solving the full multiphase Navier- Stokes equations coupled to the various equations of state and a level set interface tracking scheme based on the refined level set grid method. The codes have been parallelized using MPI in order to take advantage of today's very large parallel computational architectures. In the first system, the code's ability to handle surface tension and large tem- perature gradients is established. In the second system, the code's ability to sim- ulate simple interface geometries with strong shear is demonstrated. In the third system, the ability to handle extremely complex geometries and topology changes with strong shear is shown.
ContributorsBrady, Peter, Ph.D (Author) / Herrmann, Marcus (Thesis advisor) / Lopez, Juan (Thesis advisor) / Adrian, Ronald (Committee member) / Calhoun, Ronald (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2011