Theses and Dissertations
This collection includes both ASU Theses and Dissertations, submitted by graduate students, and the Barrett, Honors College theses submitted by undergraduate students.
Displaying 1 - 10 of 13
Filtering by
- All Subjects: Aerospace Engineering
Description
The Magnetoplasmadynamic (MPD) thruster is an electromagnetic thruster that produces a higher specific impulse than conventional chemical rockets and greater thrust densities than electrostatic thrusters, but the well-known operational limit---referred to as ``onset"---imposes a severe limitation efficiency and lifetime. This phenomenon is associated with large fluctuations in operating voltage, high rates of electrode erosion, and three-dimensional instabilities in the plasma flow-field which cannot be adequately represented by two-dimensional, axisymmetric models. Simulations of the Princeton Benchmark Thruster (PBT) were conducted using the three-dimensional version of the magnetohydrodynamic (MHD) code, MACH. Validation of the numerical model is partially achieved by comparison to equivalent simulations conducted using the well-established two-dimensional, axisymmetric version of MACH. Comparisons with available experimental data was subsequently performed to further validate the model and gain insights into the physical processes of MPD acceleration. Thrust, plasma voltage, and plasma flow-field predictions were calculated for the PBT operating with applied currents in the range $6.5kA < J < 23.25kA$ and mass-flow rates of $1g/s$, $3g/s$, and $6g/s$. Comparisons of performance characteristics between the two versions of the code show excellent agreement, indicating that MACH3 can be expected to be as predictive as MACH2 has demonstrated over multiple applications to MPD thrusters. Predicted thrust for operating conditions within the range which exhibited no symptoms of the onset phenomenon experimentally also showed agreement between MACH3 and experiment well within the experimental uncertainty. At operating conditions beyond such values , however, there is a discrepancy---up to $\sim20\%$---which implies that certain significant physical processes associated with onset are not currently being modeled. Such processes are also evident in the experimental total voltage data, as is evident by the characteristic ``voltage hash", but not present in predicted plasma voltage. Additionally, analysis of the predicted plasma flow-field shows no breakdown in azimuthal symmetry, which is expected to be associated with onset. This implies that perhaps certain physical processes are modeled by neither MACH2 nor MACH3; the latter indicating that such phenomenon may not be inherently three dimensional and related to the plasma---as suggested by other efforts---but rather a consequence of electrode material processes which have not been incorporated into the current models.
ContributorsParma, Brian (Author) / Mikellides, Pavlos G (Thesis advisor) / Squires, Kyle (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2011
Description
The evolution of single hairpin vortices and multiple interacting hairpin vortices are studied in direct numerical simulations of channel flow at Re-tau=395. The purpose of this study is to observe the effects of increased Reynolds number and varying initial conditions on the growth of hairpins and the conditions under which single hairpins autogenerate hairpin packets. The hairpin vortices are believed to provide a unified picture of wall turbulence and play an important role in the production of Reynolds shear stress which is directly related to turbulent drag. The structures of the initial three-dimensional vortices are extracted from the two-point spatial correlation of the fully turbulent direct numerical simulation of the velocity field by linear stochastic estimation and embedded in a mean flow having the profile of the fully turbulent flow. The Reynolds number of the present simulation is more than twice that of the Re-tau=180 flow from earlier literature and the conditional events used to define the stochastically estimated single vortex initial conditions include a number of new types of events such as quasi-streamwise vorticity and Q4 events. The effects of parameters like strength, asymmetry and position are evaluated and compared with existing results in the literature. This study then attempts to answer questions concerning how vortex mergers produce larger scale structures, a process that may contribute to the growth of length scale with increasing distance from the wall in turbulent wall flows. Multiple vortex interactions are studied in detail.
ContributorsParthasarathy, Praveen Kumar (Author) / Adrian, Ronald (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2011
Description
A new theoretical model was developed utilizing energy conservation methods in order to determine the fully-atomized cross-sectional Sauter mean diameters of pressure-swirl atomizers. A detailed boundary-layer assessment led to the development of a new viscous dissipation model for droplets in the spray. Integral momentum methods were also used to determine the complete velocity history of the droplets and entrained gas in the spray. The model was extensively validated through comparison with experiment and it was found that the model could predict the correct droplet size with high accuracy for a wide range of operating conditions. Based on detailed analysis, it was found that the energy model has a tendency to overestimate the droplet diameters for very low injection velocities, Weber numbers, and cone angles. A full parametric study was also performed in order to unveil some underlying behavior of pressure-swirl atomizers. It was found that at high injection velocities, the kinetic energy in the spray is significantly larger than the surface tension energy, therefore, efforts into improving atomization quality by changing the liquid's surface tension may not be the most productive. From the parametric studies it was also shown how the Sauter mean diameter and entrained velocities vary with increasing ambient gas density. Overall, the present energy model has the potential to provide quick and reasonably accurate solutions for a wide range of operating conditions enabling the user to determine how different injection parameters affect the spray quality.
ContributorsMoradi, Ali (Author) / Lee, Taewoo (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2013
Description
This study performs numerical modeling for the climate of semi-arid regions by running a high-resolution atmospheric model constrained by large-scale climatic boundary conditions, a practice commonly called climate downscaling. These investigations focus especially on precipitation and temperature, quantities that are critical to life in semi-arid regions. Using the Weather Research and Forecast (WRF) model, a non-hydrostatic geophysical fluid dynamical model with a full suite of physical parameterization, a series of numerical sensitivity experiments are conducted to test how the intensity and spatial/temporal distribution of precipitation change with grid resolution, time step size, the resolution of lower boundary topography and surface characteristics. Two regions, Arizona in U.S. and Aral Sea region in Central Asia, are chosen as the test-beds for the numerical experiments: The former for its complex terrain and the latter for the dramatic man-made changes in its lower boundary conditions (the shrinkage of Aral Sea). Sensitivity tests show that the parameterization schemes for rainfall are not resolution-independent, thus a refinement of resolution is no guarantee of a better result. But, simulations (at all resolutions) do capture the inter-annual variability of rainfall over Arizona. Nevertheless, temperature is simulated more accurately with refinement in resolution. Results show that both seasonal mean rainfall and frequency of extreme rainfall events increase with resolution. For Aral Sea, sensitivity tests indicate that while the shrinkage of Aral Sea has a dramatic impact on the precipitation over the confine of (former) Aral Sea itself, its effect on the precipitation over greater Central Asia is not necessarily greater than the inter-annual variability induced by the lateral boundary conditions in the model and large scale warming in the region. The numerical simulations in the study are cross validated with observations to address the realism of the regional climate model. The findings of this sensitivity study are useful for water resource management in semi-arid regions. Such high spatio-temporal resolution gridded-data can be used as an input for hydrological models for regions such as Arizona with complex terrain and sparse observations. Results from simulations of Aral Sea region are expected to contribute to ecosystems management for Central Asia.
ContributorsSharma, Ashish (Author) / Huang, Huei-Ping (Thesis advisor) / Adrian, Ronald (Committee member) / Herrmann, Marcus (Committee member) / Phelan, Patrick E. (Committee member) / Vivoni, Enrique (Committee member) / Arizona State University (Publisher)
Created2012
Description
Passive flow control achieved by surface dimpling can be an effective strategy for reducing drag around bluff bodies - an example of substantial popular interest being the flow around a golf ball. While the general effect of dimples causing a delay of boundary layer separation is well known, the mechanisms contributing to this phenomena are subtle and not thoroughly understood. Numerical models offer a powerful approach for studying drag reduction, however simulation strategies are challenged by complex geometries, and in applications the introduction of ad hoc turbulence models which introduce additional uncertainty. These and other factors provide much of the motivation for the current study, which focused on the numerical simulations of the flow over a simplified configuration consisting of a dimpled flat plate. The principal goals of the work are to understand the performance of the numerical methodology, and gain insight into the underlying physics of the flow. Direct numerical simulation of the incompressible Navier-Stokes equations using a fractional step method was employed, with the dimpled flat plate represented using an immersed boundary method. The dimple geometry utilizes a fixed dimple aspect ratio, with dimples arranged in a single spanwise row. The grid sizes considered ranged from approximately 3 to 99 million grid points. Reynolds numbers of 3000 and 4000 based on the inlet laminar boundary layer thickness were simulated. A turbulent boundary layer was induced downstream of the dimples for Reynolds numbers which did not transition for the flow over an undimpled flat plate. First and second order statistics of the boundary layer that develops agree reasonably well with those for turbulent channel flow and flat plate boundary layers in the sublayer and buffer layers, but differ in the outer layer. Inspection of flow visualizations suggest that early transition is promoted by thinning of the boundary layer, initiation of shear layer instabilities over the dimples, flow separation and reattachment, and tripping of the boundary layer at the trailing edge of the dimples.
ContributorsMode, Jeffrey Michael (Author) / Squires, Kyle (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2010
Description
As the push to develop ever more efficient aircraft increases, the use of lightweight composite materials to meet this push has increased. Traditional aircraft structural component sizing has revolved around the tensile yield strength of materials. Since composite materials excel in tensile strength, these traditional sizing tools provide overly optimistic weight reduction predictions. Furthermore, composite materials, in general, are weak under compression and shear. Thus, proper structural sizing yields heavier-than-expected designs. Nevertheless, a wing using thin, lightweight composites in the primary load-bearing components significantly impacts its static aeroelastic properties. These thin structures have a decreased flexural rigidity, making them more susceptible to bending. The bending of swept wings decreases the design wing twist and dihedral angle, potentially impacting the aerodynamic performance and the lateral stability and control, respectively. This work aims to determine what, if any, are the effects of excessive static aeroelastic properties on the aerodynamic performance of an aircraft. Does the perceived gain in the theoretical reduction in structural weight outweigh the potential reduction in aerodynamic performance?
ContributorsWebb, Benjamin David (Author) / Takahashi, Timothy (Thesis advisor) / Herrmann, Marcus (Committee member) / Perez, Ruben (Committee member) / Arizona State University (Publisher)
Created2022
Description
This study identifies the influence that leading-edge shape has on the aerodynamic characteristics of a wing using surface far-field and near-field analysis. It examines if a wake survey is the appropriate means for measuring profile drag and induced drag. The paper unveils the differences between sharp leading-edge and blunt leading-edge wings with the tools of pressure loop, chordwise pressure distribution, span load plots and with wake integral computations. The analysis was performed using Computational Fluid Dynamics (CFD), vortex lattice potential flow code (VORLAX), and a few wind-tunnels runs to acquire data for comparison. This study found that sharp leading-edge wings have less leading-edge suction and higher drag than blunt leading-edge wings.
The blunt leading-edge wings have less drag because the normal vector of the surface in the front section of the airfoil develops forces at opposed skin friction. The shape of the leading edge, in conjunction with the effect of viscosity, slightly alter the span load; both the magnitude of the lift and the transverse distribution. Another goal in this study is to verify the veracity of wake survey theory; the two different leading-edge shapes reveals the shortcoming of Mclean’s equation which is only applicable to blunt leading-edge wings.
The blunt leading-edge wings have less drag because the normal vector of the surface in the front section of the airfoil develops forces at opposed skin friction. The shape of the leading edge, in conjunction with the effect of viscosity, slightly alter the span load; both the magnitude of the lift and the transverse distribution. Another goal in this study is to verify the veracity of wake survey theory; the two different leading-edge shapes reveals the shortcoming of Mclean’s equation which is only applicable to blunt leading-edge wings.
ContributorsOu, Che Wei (Author) / Takahashi, Timothy (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
Description
The study of deflagration to detonation transition (DDT) in explosives is of prime importance with regards to insensitive munitions (IM). Critical damage owing to thermal or shock stimuli could translate to significant loss of life and material. The present study models detonation and deflagration of a commonly used granular explosive: cyclotetramethylene-tetranitramine, HMX. A robust literature review is followed by computational modeling of gas gun and DDT tube test data using the Sandia National Lab three-dimensional multi-material Eulerian hydrocode CTH. This dissertation proposes new computational practices and models that aid in predicting shock stimulus IM response. CTH was first used to model experimental data sets of DDT tubes from both Naval Surface Weapons Center and Los Alamos National Laboratory which were initiated by pyrogenic material and a piston, respectively. Analytical verification was performed, where possible, for detonation via empirical based equations at the Chapman Jouguet state with errors below 2.1%, and deflagration via pressure dependent burn rate equations. CTH simulations include inert, history variable reactive burn and Arrhenius models. The results are in excellent agreement with published HMX detonation velocities. Novel additions include accurate simulation of the pyrogenic material BKNO3 and the inclusion of porosity in energetic materials. The treatment of compaction is especially important in modeling precursory hotspots, caused by hydrodynamic collapse of void regions or grain interactions, prior to DDT of granular explosives. The CTH compaction model of HMX was verified within 11% error via a five pronged validation approach using gas gun data and employed use of a newly generated set of P-α parameters for granular HMX in a Mie-Gruneisen Equation of State. Next, the additions of compaction were extended to a volumetric surface burning model of HMX and compare well to a set of empirical burn rates. Lastly, the compendium of detonation and deflagration models was applied to the aforementioned DDT tubes and demonstrate working functionalities of all models, albeit at the expense of significant computational resources. A robust hydrocode methodology is proposed to make use of the deflagration, compaction and detonation models as a means to predict IM response to shock stimulus of granular explosive materials.
ContributorsMahon, Kelly Susan (Author) / Lee, Taewoo (Thesis advisor) / Herrmann, Marcus (Committee member) / Chen, Kangping (Committee member) / Jiao, Yang (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2015
Description
A moving overlapping mesh methodology that achieves spectral accuracy in space and up to second-order accuracy in time is developed for solution of unsteady incompressible flow equations in three-dimensional domains. The targeted applications are in aerospace and mechanical engineering domains and involve problems in turbomachinery, rotary aircrafts, wind turbines and others. The methodology is built within the dual-session communication framework initially developed for stationary overlapping meshes. The methodology employs semi-implicit spectral element discretization of equations in each subdomain and explicit treatment of subdomain interfaces with spectrally-accurate spatial interpolation and high-order accurate temporal extrapolation, and requires few, if any, iterations, yet maintains the global accuracy and stability of the underlying flow solver. Mesh movement is enabled through the Arbitrary Lagrangian-Eulerian formulation of the governing equations, which allows for prescription of arbitrary velocity values at discrete mesh points.
The stationary and moving overlapping mesh methodologies are thoroughly validated using two- and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal global convergence, for both methods, is documented and is in agreement with the nominal order of accuracy of the underlying solver.
Stationary overlapping mesh methodology was validated to assess the influence of long integration times and inflow-outflow global boundary conditions on the performance. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics are validated against the available data.
Moving overlapping mesh simulations are validated on the problems of two-dimensional oscillating cylinder and a three-dimensional rotating sphere. The aerodynamic forces acting on these moving rigid bodies are determined, and all results are compared with published data. Scaling tests, with both methodologies, show near linear strong scaling, even for moderately large processor counts.
The moving overlapping mesh methodology is utilized to investigate the effect of an upstream turbulent wake on a three-dimensional oscillating NACA0012 extruded airfoil. A direct numerical simulation (DNS) at Reynolds Number 44,000 is performed for steady inflow incident upon the airfoil oscillating between angle of attack 5.6 and 25 degrees with reduced frequency k=0.16. Results are contrasted with subsequent DNS of the same oscillating airfoil in a turbulent wake generated by a stationary upstream cylinder.
The stationary and moving overlapping mesh methodologies are thoroughly validated using two- and three-dimensional benchmark problems in laminar and turbulent flows. The spatial and temporal global convergence, for both methods, is documented and is in agreement with the nominal order of accuracy of the underlying solver.
Stationary overlapping mesh methodology was validated to assess the influence of long integration times and inflow-outflow global boundary conditions on the performance. In a turbulent benchmark of fully-developed turbulent pipe flow, the turbulent statistics are validated against the available data.
Moving overlapping mesh simulations are validated on the problems of two-dimensional oscillating cylinder and a three-dimensional rotating sphere. The aerodynamic forces acting on these moving rigid bodies are determined, and all results are compared with published data. Scaling tests, with both methodologies, show near linear strong scaling, even for moderately large processor counts.
The moving overlapping mesh methodology is utilized to investigate the effect of an upstream turbulent wake on a three-dimensional oscillating NACA0012 extruded airfoil. A direct numerical simulation (DNS) at Reynolds Number 44,000 is performed for steady inflow incident upon the airfoil oscillating between angle of attack 5.6 and 25 degrees with reduced frequency k=0.16. Results are contrasted with subsequent DNS of the same oscillating airfoil in a turbulent wake generated by a stationary upstream cylinder.
ContributorsMerrill, Brandon Earl (Author) / Peet, Yulia (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Calhoun, Ronald (Committee member) / Arizona State University (Publisher)
Created2016
Description
The effect of reduced frequency on dynamic stall behavior of a pitching NACA0012 airfoil in a turbulent wake using Direct Numerical Simulations is presented in the current study. Upstream turbulence with dynamically oscillating blades and airfoils is associated with ambient flow unsteadiness and is encountered in many operating conditions. Wake turbulence, a more realistic scenario for airfoils in operation, is generated using a small solid cylinder placed upstream, the vortices shed from which interact with the pitching airfoil affecting dynamic stall behavior.
A recently developed moving overlapping grid approach is used using a high-order Spectral Element Method (SEM) for spatial discretization combined with a dynamic time-stepping procedure allowing for up to third order temporal discretization. Two cases of reduced frequency (k = 0:16 and 0:25) for airfoil oscillation are investigated and the change in dynamic stall behavior with change in reduced frequency is studied and documented using flow-fields and aerodynamic coefficients (Drag, Lift and Pitching Moment) with a focus on understanding vortex system dynamics (including formation of secondary vortices) for different reduced frequencies and it’s affect on airfoil aerodynamic characteristics and fatigue life. Transition of the flow over the surface of an airfoil for both undisturbed and disturbed flow cases will also be discussed using Pressure coefficient and Skin Friction coefficient data for a given cycle combined with a wavelet analysis using Morse wavelets in MATLAB.
A recently developed moving overlapping grid approach is used using a high-order Spectral Element Method (SEM) for spatial discretization combined with a dynamic time-stepping procedure allowing for up to third order temporal discretization. Two cases of reduced frequency (k = 0:16 and 0:25) for airfoil oscillation are investigated and the change in dynamic stall behavior with change in reduced frequency is studied and documented using flow-fields and aerodynamic coefficients (Drag, Lift and Pitching Moment) with a focus on understanding vortex system dynamics (including formation of secondary vortices) for different reduced frequencies and it’s affect on airfoil aerodynamic characteristics and fatigue life. Transition of the flow over the surface of an airfoil for both undisturbed and disturbed flow cases will also be discussed using Pressure coefficient and Skin Friction coefficient data for a given cycle combined with a wavelet analysis using Morse wavelets in MATLAB.
ContributorsGandhi, Anurag (Author) / Peet, Yulia (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2017