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

Clear all filters

150092-Thumbnail Image.png
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

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
189204-Thumbnail Image.png
Description
Computing the fluid phase interfaces in multiphase flow is a challenging area of research in fluids. The Volume of Fluid andLevel Set methods are a few algorithms that have been developed for reconstructing the multiphase fluid flow interfaces. The thesis work focuses on exploring the ability of neural networks to reconstruct

Computing the fluid phase interfaces in multiphase flow is a challenging area of research in fluids. The Volume of Fluid andLevel Set methods are a few algorithms that have been developed for reconstructing the multiphase fluid flow interfaces. The thesis work focuses on exploring the ability of neural networks to reconstruct the multiphase fluid flow interfaces using a data-driven approach. The neural network model has liquid volume fraction stencils as an input, and it predicts the radius of the circle as an output of the network which represents a phase interface separating two immiscible fluids inside a fluid domain. The liquid volume fraction stencils are generated for randomly varying circle radii within a 1x1 domain using an open-source VOFI library. These datasets are used to train the neural network. Once the model is trained, the predicted circular phase interface from the neural network output is used to generate back the predicted liquid volume fraction stencils. Error norms values are calculated to assess the error in the neural network model’s predicted liquid volume fraction stencils with the actual liquid volume fraction stencils from the VOFI library. The neural network parameters are optimized by testing them for different hyper-parameters to reduce the error norms. So as to minimize the difference between the predicted and the actual liquid volume fraction stencils and errors in reconstructing the fluid phase interface geometry.
ContributorsPawar, Pranav Rajesh (Author) / Herrmann, Marcus (Thesis advisor) / Zhuang, Houlong (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2023
171359-Thumbnail Image.png
Description
Theoretical analyses of liquid atomization (bulk to droplet conversion) and turbulence have potential to advance the computability of these flows. Instead of relying on full computations or models, fundamental conservation equations can be manipulated to generate partial or full solutions. For example, integral form of the mass and energy for

Theoretical analyses of liquid atomization (bulk to droplet conversion) and turbulence have potential to advance the computability of these flows. Instead of relying on full computations or models, fundamental conservation equations can be manipulated to generate partial or full solutions. For example, integral form of the mass and energy for spray flows leads to an explicit relationship between the drop size and liquid velocities. This is an ideal form to integrate with existing computational fluid dynamic (CFD), which is well developed to solve for the liquid velocities, i.e., the momentum equation(s). Theoretical adaption to CFD has been performed for various injection geometries, with results that compare quite well with experimental data. Since the drop size is provided analytically, computational time/cost for simulating spray flows with liquid atomization is no more than single-phase flows. Some advances have also been made on turbulent flows, by using a new set of perspectives on transport, scaling and energy distributions. Conservation equations for turbulence momentum and kinetic energy have been derived in a coordinate frame moving with the local mean velocities, which produce the Reynolds stress components, without modeling. Scaling of the Reynolds stress is also found at the first- and second-gradient levels. Finally, maximum-entropy principle has been used to derive the energy spectra in turbulent flows.
ContributorsPark, Jung Eun (Author) / Lee, Taewoo (Thesis advisor) / Gardner, Carl (Committee member) / Huang, Huei-Ping (Committee member) / Kim, Jeonglae (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2022
171812-Thumbnail Image.png
Description
The propulsion matrix provides a compact description of the locomotion of a single flagella molecular motor in a low Reynolds number environment. The locomotion properties of individual flagellar motors are central to bacterial behavior, including chemotaxis, pathogenesis, and biofilm formation. However, because conventional hydrodynamic measurement approaches require applied forces, torques,

The propulsion matrix provides a compact description of the locomotion of a single flagella molecular motor in a low Reynolds number environment. The locomotion properties of individual flagellar motors are central to bacterial behavior, including chemotaxis, pathogenesis, and biofilm formation. However, because conventional hydrodynamic measurement approaches require applied forces, torques, or fluid flows, it is not possible to directly measure the propulsion matrix for an individual microscale helical filament. Here, the limitations inherent to conventional measurement approaches are overcome using a combination of theoretical, experimental, and computational advancements. First, the relationship between the elements of the propulsion matrix with translational and rotational Brownian motion is derived using the fluctuation-dissipation theorem. Next, a volumetric fluorescent imaging using high resolution oblique plane microscopy with sufficient spatio-temporal resolution is conducted to resolve both translation and rotation of individual helical filaments isolated from E.coli's flagellar motor. Finally, a computational framework is developed to track individual helical filaments across six degrees of freedom, extract diffusion coefficients, and quantify the temporal correlation between translation and rotation. This study computed the maximum propulsion efficiency to be around 1.7%. Direct measurement of propulsion efficiency generally agrees with the ensemble and large-scale measurements previously performed using conventional hydrodynamic measurements. The findings suggest that the approach described here can be extended to more complex in-vitro experiments that evaluate microscale molecular motors. For example, evaluating sperm motility without inducing chemotaxis or utilizing a microfluidic setup.
ContributorsDjutanta, Franky (Author) / Hariadi, Rizal (Thesis advisor) / Wang, Robert (Thesis advisor) / Yurke, Bernard (Committee member) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2022
168292-Thumbnail Image.png
Description
In this dissertation, two types of passive air freshener products from Henkel, the wick-based air freshener and gel-based air freshener, are studied for their wicking mechanisms and evaporation performances.The fibrous pad of the wick-based air freshener is a porous medium that absorbs fragrance by capillary force and releases the fragrance

In this dissertation, two types of passive air freshener products from Henkel, the wick-based air freshener and gel-based air freshener, are studied for their wicking mechanisms and evaporation performances.The fibrous pad of the wick-based air freshener is a porous medium that absorbs fragrance by capillary force and releases the fragrance into the ambient air. To investigate the wicking process, a two-dimensional multiphase flow numerical model using COMSOL Multiphysics is built. Saturation and liquid pressure inside the pad are solved. Comparison between the simulation results and experiments shows that evaporation occurs simultaneously with the wicking process. The evaporation performance on the surface of the wicking pad is analyzed based on the kinetic theory, from which the mass flow rate of molecules passing the interface of each pore of the porous medium is obtained. A 3D model coupling the evaporation model and dynamic wicking on the evaporation pad is built to simulate the entire performance of the air freshener to the environment for a long period of time. Diffusion and natural convection effects are included in the simulation. The simulation results match well with the experiments for both the air fresheners placed in a chamber and in the absent of a chamber, the latter of which is subject to indoor airflow. The gel-based air freshener can be constructed as a porous medium in which the solid network of particles spans the volume of the fragrance liquid. To predict the evaporation performance of the gel, two approaches are tested for gel samples in hemispheric shape. The first approach is the sessile drop model commonly used for the drying process of a pure liquid droplet. It can be used to estimate the weight loss rate and time duration of the evaporation. Another approach is to simulate the concentration profile outside the gel and estimate the evaporation rate from the surface of the gel using the kinetic theory. The evaporation area is updated based on the change of pore size. A 3D simulation using the same analysis is further applied to the cylindrical gel sample. The simulation results match the experimental data well.
ContributorsYuan, Jing (Author) / Chen, Kangping (Thesis advisor) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Wang, Liping (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2021
187669-Thumbnail Image.png
Description
Advancements to a dual scale Large Eddy Simulation (LES) modeling approach for immiscible turbulent phase interfaces are presented. In the dual scale LES approach, a high resolution auxiliary grid, used to capture a fully resolved interface geometry realization, is linked to an LES grid that solves the filtered Navier-Stokes equations.

Advancements to a dual scale Large Eddy Simulation (LES) modeling approach for immiscible turbulent phase interfaces are presented. In the dual scale LES approach, a high resolution auxiliary grid, used to capture a fully resolved interface geometry realization, is linked to an LES grid that solves the filtered Navier-Stokes equations. Exact closure of the sub-filter interface terms is provided by explicitly filtering the fully resolved quantities from the auxiliary grid. Reconstructing a fully resolved velocity field to advance the phase interface requires modeling several sub-filter effects, including shear and accelerational instabilities and phase change. Two sub-filter models were developed to generate these sub-filter hydrodynamic instabilities: an Orr-Sommerfeld model and a Volume-of-Fluid (VoF) vortex sheet method. The Orr-Sommerfeld sub-filter model was found to be incompatible with the dual scale approach, since it is unable to generate interface rollup and a process to separate filtered and sub-filter scales could not be established. A novel VoF vortex sheet method was therefore proposed, since prior vortex methods have demonstrated interface rollup and following the LES methodology, the vortex sheet strength could be decomposed into its filtered and sub-filter components. In the development of the VoF vortex sheet method, it was tested with a variety of classical hydrodynamic instability problems, compared against prior work and linear theory, and verified using Direct Numerical Simulations (DNS). An LES consistent approach to coupling the VoF vortex sheet with the LES filtered equations is presented and compared against DNS. Finally, a sub-filter phase change model is proposed and assessed in the dual scale LES framework with an evaporating interface subjected to decaying homogeneous isotropic turbulence. Results are compared against DNS and the interplay between surface tension forces and evaporation are discussed.
ContributorsGoodrich, Austin Chase (Author) / Herrmann, Marcus (Thesis advisor) / Dahm, Werner (Committee member) / Kim, Jeonglae (Committee member) / Huang, Huei-Ping (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2023
193358-Thumbnail Image.png
Description
This dissertation investigates the complex dynamics of semi-dilute inertial particles suspended in vortices using the Eulerian-Lagrangian method. The study explores the modulation of flow induced by inertial particles, focusing on the characteristics of a single vortex, instability analysis within particle-laden flows, and the merging process of co-rotating vortices. Simulations reveal

This dissertation investigates the complex dynamics of semi-dilute inertial particles suspended in vortices using the Eulerian-Lagrangian method. The study explores the modulation of flow induced by inertial particles, focusing on the characteristics of a single vortex, instability analysis within particle-laden flows, and the merging process of co-rotating vortices. Simulations reveal a preferential concentration mechanism, where inertial particles cluster around a void fraction bubble, accelerating the decay of the vortex tube. Small-scale clusters, arising from particle-trajectory crossings, induce significant gradients in the fluid vorticity field, contributing to rapid vortex breakdown. Within a specific Stokes number range, increased particle inertia results in faster vortex decay and stronger inhomogeneity in the particle phase. The instability mechanism in particle-laden flows is explored using a Rankine vortex model. Two-way coupling triggers azimuthal perturbations, leading to the breakdown of the vortex structure. Linear Stability Analysis and Two-Fluid modeling demonstrate that the dusty vortex flow exhibits unstable modes, with growth rates increasing with wavenumber. Eulerian-Lagrangian simulations validate these results, showing excellent agreement between computed and predicted growth rates. The dissertation also delves into the co-rotating vortex merger in a semi-dilute dusty flow. For weak inertial effects, merger experiences a delay compared to particle-free vortices. Under moderate inertial conditions, the merger process exhibits repulsion, increased separation, and eventual convective merger stages. Highly inertial particles stretch the vortex core, initiating a merger with an outcome of a particle-free vortex core surrounded by a halo of concentrated particles. In conclusion, the feedback force from the dispersed phase induces instability and significantly influences the dynamics of vortices in particle-laden flows. The findings contribute to a deeper understanding of the intricate interactions between inertial particles and vortical structures.
ContributorsShuai, Shuai (Author) / Kasbaoui, Mohamed Houssem (Thesis advisor) / Herrmann, Marcus (Committee member) / Peet, Yulia (Committee member) / Huang, Huei-Ping (Committee member) / Wang, Zhihua (Committee member) / Arizona State University (Publisher)
Created2024
Description
The goal of this paper was to do an analysis of two-dimensional unsplit mass and momentum conserving Finite Volume Methods for Advection for Volume of Fluid Fields with interfaces and validating their rates of convergence. Specifically three unsplit transport methods and one split transport method were amalgamated individually with four

The goal of this paper was to do an analysis of two-dimensional unsplit mass and momentum conserving Finite Volume Methods for Advection for Volume of Fluid Fields with interfaces and validating their rates of convergence. Specifically three unsplit transport methods and one split transport method were amalgamated individually with four Piece-wise Linear Reconstruction Schemes (PLIC) i.e. Unsplit Eulerian Advection (UEA) by Owkes and Desjardins (2014), Unsplit Lagrangian Advection (ULA) by Yang et al. (2010), Split Lagrangian Advection (SLA) by Scardovelli and Zaleski (2003) and Unsplit Averaged Eulerian-Lagrangian Advection (UAELA) with two Finite Difference Methods by Parker and Youngs (1992) and two Error Minimization Methods by Pilliod Jr and Puckett (2004). The observed order of accuracy was first order in all cases except when unsplit methods and error minimization methods were used consecutively in each iteration, which resulted in second-order accuracy on the shape error convergence. The Averaged Unsplit Eulerian-Lagrangian Advection (AUELA) did produce first-order accuracy but that was due to a temporal error in the numerical setup. The main unsplit methods, Unsplit Eulerian Advection (UEA) and Unsplit Lagrangian Advection (ULA), preserve mass and momentum and require geometric clipping to solve two-phase fluid flows. The Unsplit Lagrangian Advection (ULA) can allow for small divergence in the velocity field perhaps saving time on the iterative solver of the variable coefficient Poisson System.
ContributorsAnsari, Adil (M.S.) (Author) / Herrmann, Marcus (Thesis advisor) / Peet, Yulia (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
Description
Rapid expansion of dense beds of fine, spherical particles subjected to rapid depressurization is studied in a vertical shock tube. As the particle bed is unloaded, a high-speed video camera captures the dramatic evolution of the particle bed structure. Pressure transducers are used to measure the dynamic pressure changes during

Rapid expansion of dense beds of fine, spherical particles subjected to rapid depressurization is studied in a vertical shock tube. As the particle bed is unloaded, a high-speed video camera captures the dramatic evolution of the particle bed structure. Pressure transducers are used to measure the dynamic pressure changes during the particle bed expansion process. Image processing, signal processing, and Particle Image Velocimetry techniques, are used to examine the relationships between particle size, initial bed height, bed expansion rate, and gas velocities.

The gas-particle interface and the particle bed as a whole expand and evolve in stages. First, the bed swells nearly homogeneously for a very brief period of time (< 2ms). Shortly afterward, the interface begins to develop instabilities as it continues to rise, with particles nearest the wall rising more quickly. Meanwhile, the bed fractures into layers and then breaks down further into cellular-like structures. The rate at which the structural evolution occurs is shown to be dependent on particle size. Additionally, the rate of the overall bed expansion is shown to be dependent on particle size and initial bed height.

Taller particle beds and beds composed of smaller-diameter particles are found to be associated with faster bed-expansion rates, as measured by the velocity of the gas-particle interface. However, the expansion wave travels more slowly through these same beds. It was also found that higher gas velocities above the the gas-particle interface measured \textit{via} Particle Image Velocimetry or PIV, were associated with particle beds composed of larger-diameter particles. The gas dilation between the shocktube diaphragm and the particle bed interface is more dramatic when the distance between the gas-particle interface and the diaphragm is decreased-as is the case for taller beds.

To further elucidate the complexities of this multiphase compressible flow, simple OpenFOAM (Weller, 1998) simulations of the shocktube experiment were performed and compared to bed expansion rates, pressure fluctuations, and gas velocities. In all cases, the trends and relationships between bed height, particle diameter, with expansion rates, pressure fluctuations and gas velocities matched well between experiments and simulations. In most cases, the experimentally-measured bed rise rates and the simulated bed rise rates matched reasonably well in early times. The trends and overall values of the pressure fluctuations and gas velocities matched well between the experiments and simulations; shedding light on the effects each parameter has on the overall flow.
ContributorsZunino, Heather (Author) / Adrian, Ronald J (Thesis advisor) / Clarke, Amanda (Committee member) / Chen, Kangping (Committee member) / Herrmann, Marcus (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2019
153520-Thumbnail Image.png
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

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