ASU Electronic Theses and Dissertations
This collection includes most of the ASU Theses and Dissertations from 2011 to present. ASU Theses and Dissertations are available in downloadable PDF format; however, a small percentage of items are under embargo. Information about the dissertations/theses includes degree information, committee members, an abstract, supporting data or media.
In addition to the electronic theses found in the ASU Digital Repository, ASU Theses and Dissertations can be found in the ASU Library Catalog.
Dissertations and Theses granted by Arizona State University are archived and made available through a joint effort of the ASU Graduate College and the ASU Libraries. For more information or questions about this collection contact or visit the Digital Repository ETD Library Guide or contact the ASU Graduate College at gradformat@asu.edu.
Displaying 1 - 9 of 9
Filtering by
- Creators: Herrmann, Marcus
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
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
Description
Structural features of canonical wall-bounded turbulent flows are described using several techniques, including proper orthogonal decomposition (POD). The canonical wall-bounded turbulent flows of channels, pipes, and flat-plate boundary layers include physics important to a wide variety of practical fluid flows with a minimum of geometric complications. Yet, significant questions remain for their turbulent motions' form, organization to compose very long motions, and relationship to vortical structures. POD extracts highly energetic structures from flow fields and is one tool to further understand the turbulence physics. A variety of direct numerical simulations provide velocity fields suitable for detailed analysis. Since POD modes require significant interpretation, this study begins with wall-normal, one-dimensional POD for a set of turbulent channel flows. Important features of the modes and their scaling are interpreted in light of flow physics, also leading to a method of synthesizing one-dimensional POD modes. Properties of a pipe flow simulation are then studied via several methods. The presence of very long streamwise motions is assessed using a number of statistical quantities, including energy spectra, which are compared to experiments. Further properties of energy spectra, including their relation to fictitious forces associated with mean Reynolds stress, are considered in depth. After reviewing salient features of turbulent structures previously observed in relevant experiments, structures in the pipe flow are examined in greater detail. A variety of methods reveal organization patterns of structures in instantaneous fields and their associated vortical structures. Properties of POD modes for a boundary layer flow are considered. Finally, very wide modes that occur when computing POD modes in all three canonical flows are compared. The results demonstrate that POD extracts structures relevant to characterizing wall-bounded turbulent flows. However, significant care is necessary in interpreting POD results, for which modes can be categorized according to their self-similarity. Additional analysis techniques reveal the organization of smaller motions in characteristic patterns to compose very long motions in pipe flows. The very large scale motions are observed to contribute large fractions of turbulent kinetic energy and Reynolds stress. The associated vortical structures possess characteristics of hairpins, but are commonly distorted from pristine hairpin geometries.
ContributorsBaltzer, Jon Ronald (Author) / Adrian, Ronald J (Thesis advisor) / Calhoun, Ronald (Committee member) / Gelb, Anne (Committee member) / Herrmann, Marcus (Committee member) / Squires, Kyle D (Committee member) / Arizona State University (Publisher)
Created2012
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
This study uses the Weather Research and Forecasting (WRF) model to simulate and predict the changes in local climate attributed to the urbanization for five desert cities. The simulations are performed in the fashion of climate downscaling, constrained by the surface boundary conditions generated from high resolution land-use maps. For each city, the land-use maps of 1985 and 2010 from Landsat satellite observation, and a projected land-use map for 2030, are used to represent the past, present, and future. An additional set of simulations for Las Vegas, the largest of the five cities, uses the NLCD 1992 and 2006 land-use maps and an idealized historical land-use map with no urban coverage for 1900.
The study finds that urbanization in Las Vegas produces a classic urban heat island (UHI) at night but a minor cooling during the day. A further analysis of the surface energy balance shows that the decrease in surface Albedo and increase effective emissivity play an important role in shaping the local climate change over urban areas. The emerging urban structures slow down the diurnal wind circulation over the city due to an increased effective surface roughness. This leads to a secondary modification of temperature due to the interaction between the mechanical and thermodynamic effects of urbanization.
The simulations for the five desert cities for 1985 and 2010 further confirm a common pattern of the climatic effect of urbanization with significant nighttime warming and moderate daytime cooling. This effect is confined to the urban area and is not sensitive to the size of the city or the detail of land cover in the surrounding areas. The pattern of nighttime warming and daytime cooling remains robust in the simulations for the future climate of the five cities using the projected 2030 land-use maps. Inter-city differences among the five urban areas are discussed.
The study finds that urbanization in Las Vegas produces a classic urban heat island (UHI) at night but a minor cooling during the day. A further analysis of the surface energy balance shows that the decrease in surface Albedo and increase effective emissivity play an important role in shaping the local climate change over urban areas. The emerging urban structures slow down the diurnal wind circulation over the city due to an increased effective surface roughness. This leads to a secondary modification of temperature due to the interaction between the mechanical and thermodynamic effects of urbanization.
The simulations for the five desert cities for 1985 and 2010 further confirm a common pattern of the climatic effect of urbanization with significant nighttime warming and moderate daytime cooling. This effect is confined to the urban area and is not sensitive to the size of the city or the detail of land cover in the surrounding areas. The pattern of nighttime warming and daytime cooling remains robust in the simulations for the future climate of the five cities using the projected 2030 land-use maps. Inter-city differences among the five urban areas are discussed.
ContributorsKamal, Samy (Author) / Huang, Huei-Ping (Thesis advisor) / Anderson, James (Thesis advisor) / Herrmann, Marcus (Committee member) / Calhoun, Ronald (Committee member) / Myint, Soe (Committee member) / Arizona State University (Publisher)
Created2015
Description
This dissertation studies two outstanding microscale fluid mechanics problems: 1) mechanisms of gas production from the nanopores of shale; 2) enhanced mass flow rate in steady compressible gas flow through a micro-conduit.
The dissertation starts with a study of a volumetric expansion driven drainage flow of a viscous compressible fluid from a small capillary and channel in the low Mach number limit. An analysis based on the linearized compressible Navier-Stokes equations with no-slip condition shows that fluid drainage is controlled by the slow decay of the acoustic wave inside the capillary and the no-slip flow exhibits a slip-like mass flow rate. Numerical simulations are also carried out for drainage from a small capillary to a reservoir or a contraction of finite size. By allowing the density wave to escape the capillary, two wave leakage mechanisms are identified, which are dependent on the capillary length to radius ratio, reservoir size and acoustic Reynolds number. Empirical functions are generated for an effective diffusive coefficient which allows simple calculations of the drainage rate using a diffusion model without the presence of the reservoir or contraction.
In the second part of the dissertation, steady viscous compressible flow through a micro-conduit is studied using compressible Navier-Stokes equations with no-slip condition. The mathematical theory of Klainerman and Majda for low Mach number flow is employed to derive asymptotic equations in the limit of small Mach number. The overall flow, a combination of the Hagen-Poiseuille flow and a diffusive velocity shows a slip-like mass flow rate even through the overall velocity satisfies the no-slip condition. The result indicates that the classical formulation includes self-diffusion effect and it embeds the Extended Navier-Stokes equation theory (ENSE) without the need of introducing additional constitutive hypothesis or assuming slip on the boundary. Contrary to most ENSE publications, the predicted mass flow rate is still significantly below the measured data based on an extensive comparison with thirty-five experiments.
The dissertation starts with a study of a volumetric expansion driven drainage flow of a viscous compressible fluid from a small capillary and channel in the low Mach number limit. An analysis based on the linearized compressible Navier-Stokes equations with no-slip condition shows that fluid drainage is controlled by the slow decay of the acoustic wave inside the capillary and the no-slip flow exhibits a slip-like mass flow rate. Numerical simulations are also carried out for drainage from a small capillary to a reservoir or a contraction of finite size. By allowing the density wave to escape the capillary, two wave leakage mechanisms are identified, which are dependent on the capillary length to radius ratio, reservoir size and acoustic Reynolds number. Empirical functions are generated for an effective diffusive coefficient which allows simple calculations of the drainage rate using a diffusion model without the presence of the reservoir or contraction.
In the second part of the dissertation, steady viscous compressible flow through a micro-conduit is studied using compressible Navier-Stokes equations with no-slip condition. The mathematical theory of Klainerman and Majda for low Mach number flow is employed to derive asymptotic equations in the limit of small Mach number. The overall flow, a combination of the Hagen-Poiseuille flow and a diffusive velocity shows a slip-like mass flow rate even through the overall velocity satisfies the no-slip condition. The result indicates that the classical formulation includes self-diffusion effect and it embeds the Extended Navier-Stokes equation theory (ENSE) without the need of introducing additional constitutive hypothesis or assuming slip on the boundary. Contrary to most ENSE publications, the predicted mass flow rate is still significantly below the measured data based on an extensive comparison with thirty-five experiments.
ContributorsShen, Di (Author) / Chen, Kangping (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Calhoun, Ronald (Committee member) / Lopez, Juan (Committee member) / Arizona State University (Publisher)
Created2019
Description
Autonomic closure is a recently-proposed subgrid closure methodology for large eddy simulation (LES) that replaces the prescribed subgrid models used in traditional LES closure with highly generalized representations of subgrid terms and solution of a local system identification problem that allows the simulation itself to determine the local relation between each subgrid term and the resolved variables at every point and time. The present study demonstrates, for the first time, practical LES based on fully dynamic implementation of autonomic closure for the subgrid stress and the subgrid scalar flux. It leverages the inherent computational efficiency of tensorally-correct generalized representations in terms of parametric quantities, and uses the fundamental representation theory of Smith (1971) to develop complete and minimal tensorally-correct representations for the subgrid stress and scalar flux. It then assesses the accuracy of these representations via a priori tests, and compares with the corresponding accuracy from nonparametric representations and from traditional prescribed subgrid models. It then assesses the computational stability of autonomic closure with these tensorally-correct parametric representations, via forward simulations with a high-order pseudo-spectral code, including the extent to which any added stabilization is needed to ensure computational stability, and compares with the added stabilization needed in traditional closure with prescribed subgrid models. Further, it conducts a posteriori tests based on forward simulations of turbulent conserved scalar mixing with the same pseudo-spectral code, in which velocity and scalar statistics from autonomic closure with these representations are compared with corresponding statistics from traditional closure using prescribed models, and with corresponding statistics of filtered fields from direct numerical simulation (DNS). These comparisons show substantially greater accuracy from autonomic closure than from traditional closure. This study demonstrates that fully dynamic autonomic closure is a practical approach for LES that requires accuracy even at the smallest resolved scales.
ContributorsStallcup, Eric Warren (Author) / Dahm, Werner J.A. (Thesis advisor) / Herrmann, Marcus (Committee member) / Calhoun, Ronald (Committee member) / Kim, Jeonglae (Committee member) / Kostelich, Eric J. (Committee member) / Arizona State University (Publisher)
Created2020
Description
Drainage flow of a viscous compressible gas from a semi-sealed narrow conduit is a pore-scale model for studying the fundamental flow physics of fluid recovery from a porous reservoir without using fluid injection. Thermal effect has been routinely neglected for these flows in the traditional petroleum engineering literature. Since the motion is entirely driven by volumetric expansion, temperature change always accompanies the density change. This thesis examines such thermal effects on the drainage flow.
Thermal drainage flow is first studied by simultaneously solving the linearized continuity, momentum and energy equations for adiabatic walls. It is shown that even in the absence of an imposed temperature drop, gas expansion induces a transient temperature decrease inside the channel, which slows down the drainage process compared to the isothermal model and Lighthill’s model. For a given density drop, gas drains out faster as the initial-to-final temperature ratio increases; and the transient density can undershoot the final equilibrium value. A parametric study is then carried out to explore the influence of various thermal boundary conditions on drainage flow. It is found that as the wall transitions from adiabatic to isothermal condition, the excess density changes from a plane wave solution to a non-plane wave solution and the drainage rate increases. It is shown that when the exit is also cooled and the wall is non-adiabatic, the total recovered fluid mass exceeds the amount based on the isothermal theory which is determined by the initial and final density difference alone. Finally, a full numerical simulation is conducted to mimic the channel-reservoir system using the finite volume method. The Ghost-Cell Navier-Stokes Characteristic Boundary Condition technique is applied at the far end of the truncated reservoir, which is an open boundary. The results confirm the conclusions of the linear theory.
Thermal drainage flow is first studied by simultaneously solving the linearized continuity, momentum and energy equations for adiabatic walls. It is shown that even in the absence of an imposed temperature drop, gas expansion induces a transient temperature decrease inside the channel, which slows down the drainage process compared to the isothermal model and Lighthill’s model. For a given density drop, gas drains out faster as the initial-to-final temperature ratio increases; and the transient density can undershoot the final equilibrium value. A parametric study is then carried out to explore the influence of various thermal boundary conditions on drainage flow. It is found that as the wall transitions from adiabatic to isothermal condition, the excess density changes from a plane wave solution to a non-plane wave solution and the drainage rate increases. It is shown that when the exit is also cooled and the wall is non-adiabatic, the total recovered fluid mass exceeds the amount based on the isothermal theory which is determined by the initial and final density difference alone. Finally, a full numerical simulation is conducted to mimic the channel-reservoir system using the finite volume method. The Ghost-Cell Navier-Stokes Characteristic Boundary Condition technique is applied at the far end of the truncated reservoir, which is an open boundary. The results confirm the conclusions of the linear theory.
ContributorsHuang, Wei (Author) / Chen, Kangping (Thesis advisor) / Huang, Huei-Ping (Committee member) / Herrmann, Marcus (Committee member) / Calhoun, Ronald (Committee member) / Baer, Steven (Committee member) / Arizona State University (Publisher)
Created2020
Description
Large amplitude westward propagating long waves in midlatitudes of Northern Hemisphere occasionally sustain coherent phase propagation over multiple weeks. Owing to the large amplitude and the life cycle of these waves previous studies have speculated their influence on extended-range weather forecasts but have not quantified them. The primary aim of this study is to establish an updated long-term catalog of Retrograde events which can then be used to investigate the statistics and structure of these waves. Guided by the newly created catalog the dynamics of these waves are further explored. A preliminary look into the dynamics of these waves reveal a sequence of poleward extrusion, westward migration and vortex shedding occurring frequently during certain strong Retrograde wave events. A strong connection between the westward moving low PV structures and the East Asian cold air outbreak is uncovered. Also, the initiation of the sequence of low PV extrusion and vortex shedding is found to be linked with the phase of propagating Wave-1 zonal component. Enhanced predictability of global midlatitude Geopotential Height at 500mb is noted during active period of strong Retrograde wave activity in comparison to inactive period. Skilled forecasts were produced almost (on an average) 12 days in advance during the active period of one of the winters (1995/96) as compared to 9 days during the inactive period of the season.
ContributorsRaghunathan, Girish Nigamanth (Author) / Huang, Huei-Ping (Thesis advisor) / Chen, Kangping (Committee member) / Calhoun, Ronald (Committee member) / Herrmann, Marcus (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2021