Matching Items (39)
Filtering by
- All Subjects: Mechanical Engineering
- Creators: Herrmann, Marcus
- Member of: Theses and Dissertations
- Resource Type: Text
Description
Conventional fluid dynamics models such as the Navier-Stokes equations are derived for prediction of fluid motion at or near equilibrium, classic examples being the motion of fluids for which inter-molecular collisions are dominant. Flows at equilibrium permit simplifications such as the introduction of viscosity and also lead to solutions that are single-valued. However, many other regimes of interest include "fluids"' far from equilibrium; for example, rarefied gases or particle-laden flows in which the dispersed phase can be comprised of granular solids, droplets, or bubbles. Particle motion in these flows is not typically dominated by collisions and may exhibit significant memory effects; therefore, is often poorly described using continuum, field-based (Eulerian) approaches. Non-equilibrium flows generally lack a straightforward counterpart to viscosity and their multi-valued solutions cannot be represented by most Eulerian methods. This strongly motivates different strategies to address current shortcomings and the novel approach adopted in this work is based on the Conditional Quadrature Method of Moments (CQMOM). In CQMOM, moment equations are derived from the Boltzmann equation using a quadrature approximation of the velocity probability density function (PDF). CQMOM circumvents the drawbacks of current methods and leads to multivariate and multidimensional solutions in an Eulerian frame of reference. In the present work, the discretized PDF is resolved using an adaptive two-point quadrature in three-dimensional velocity space. The method is applied to computation of a series of non-equilibrium flows, ranging from simple two-dimensional test cases to fully-turbulent three-dimensional wall-bounded particle-laden flows. The primary contribution of the present effort is on development, application, and assessment of CQMOM for predicting the key features of dilute particle-laden flows. Statistical descriptors such as mean concentration and mean velocity are in good agreement with previous results, for both collision-less and collisional flows at varying particle Stokes numbers. Turbulent statistics and measures of local accumulation agree less favorably with prior results and identify areas for improvement in the modeling strategy.
ContributorsDunn, Dennis Martin (Author) / Squires, Kyle D. (Thesis advisor) / Calhoun, Ronald J. (Committee member) / Chen, Kangping (Committee member) / Dai, Lenore L. (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2015
Description
The formation and stability of a slowly evolving zonal jet in 2-D flow with beta effect is analyzed using the model developed by Manfroi and Young in which the final governing equation was derived by means of a perturbation analysis of a barotropic vorticity equation with sinusoidal meridional mean flow. However in the original study the term β0, that represents the effect of large-scale Rossby waves, was dropped and was proceeded on a path of finding solutions for a simplified 1-D flow. The idea of this study is to understand the effects of the dropped term on the overall dynamics of the zonal jet evolution. For this purpose the system that is entirely deterministic with no additional forcing is solved by means of a standard finite difference scheme. The Numerical solutions are found for varying β0 and μ values where μ represents the bottom drag. In addition to this the criteria for the formation of zonal jets developed originally for the 1-D system is verified for the 2-D system as well. The study reveals the similarity in some of the results of the 1-D and the 2-D system like the merging of jets in the absence of bottom drag, formation of steady jets in presence of a non-zero bottom drag and the adherence to the boundary criteria for the formation of zonal jets. But when it comes to the formation of steady jets, a finite β0 value is required above which the solution is similar to the 1-D system. Also the jets formed under the presence of non-zero bottom drag seem wavy in nature which is different from the steady horizontal jets produced in the 1-D system.
ContributorsRaghunathan, Girish Nigamanth (Author) / Huang, Huei-Ping (Thesis advisor) / Herrmann, Marcus (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2017
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
Description
Durable, cost-effective, and environmentally friendly anti-icing methods are desired to reduce the icing hazard in many different industrial areas including transportation systems, power plants, power transmission, as well as offshore oil and gas production. In contrast to traditional passive anti-icing surfaces, this thesis work introduces an anti-icing coating that responds to different icing conditions by releasing an antifreeze liquid. It consists of an outer porous superhydrophobic epidermis and a wick-like underlying dermis that is infused with the antifreeze liquid. This bi-layer coating prevents accumulation of frost, freezing fog, and freezing rain, while conventional anti-icing surfaces typically work only in one of these conditions. The bi-layer coating also delays condensation on the exterior surface at least ten times longer than identical system without antifreeze.
It is demonstrated that the significant delay in condensation onset is due to the integral humidity sink effect posed by the hygroscopic antifreeze liquid infused in the porous structure. This effect significantly alters the water vapor concentration field at the coating surface, which delays nucleation of drops and ice. It was demonstrated that with a proper design of the environmental chamber the size of the region of inhibited condensation and condensation frosting around an isolated pore, as well as periodically spaced pores, filled by propylene glycol can be quantitatively predicted from quasi-steady state water vapor concentration field. Theoretical analysis and experiments revealed that the inhibition of nucleation is governed by only two non-dimensional geometrical parameters: the pore size relative to the unit cell size and the ratio of the unit cell size to the thickness of the boundary layer. It is demonstrated that by switching the size of the pores from millimeters to nanometers, a dramatic depression of the nucleation onset temperature, as well as significantly greater delay in nucleation onset can be achieved.
It is demonstrated that the significant delay in condensation onset is due to the integral humidity sink effect posed by the hygroscopic antifreeze liquid infused in the porous structure. This effect significantly alters the water vapor concentration field at the coating surface, which delays nucleation of drops and ice. It was demonstrated that with a proper design of the environmental chamber the size of the region of inhibited condensation and condensation frosting around an isolated pore, as well as periodically spaced pores, filled by propylene glycol can be quantitatively predicted from quasi-steady state water vapor concentration field. Theoretical analysis and experiments revealed that the inhibition of nucleation is governed by only two non-dimensional geometrical parameters: the pore size relative to the unit cell size and the ratio of the unit cell size to the thickness of the boundary layer. It is demonstrated that by switching the size of the pores from millimeters to nanometers, a dramatic depression of the nucleation onset temperature, as well as significantly greater delay in nucleation onset can be achieved.
ContributorsSun, Xiaoda (Author) / Rykaczewski, Konrad (Thesis advisor) / Lin, Jerry (Committee member) / Phelan, Patrick (Committee member) / Wang, Robert (Committee member) / Herrmann, Marcus (Committee member) / Wang, Liping (Committee member) / Arizona State University (Publisher)
Created2017
Description
The central purpose of this work is to investigate the large-scale, coherent structures that exist in turbulent Rayleigh-Bénard convection (RBC) when the domain is large enough for the classical ”wind of turbulence” to break down. The study exclusively focuses on the structures that from when the RBC geometry is a cylinder. A series of visualization studies, Fourier analysis and proper orthogonal decomposition are employed to qualitatively and quantitatively inspect the large-scale structures’ length and time scales, spatial organization, and dynamic properties. The data in this study is generated by direct numerical simulation to resolve all the scales of turbulence in a 6.3 aspect-ratio cylinder at a Rayleigh number of 9.6 × 107 and Prandtl number of 6.7. Single and double point statistics are compared against experiments and several resolution criteria are examined to verify that the simulation has enough spatial and temporal resolution to adequately represent the physical system.
Large-scale structures are found to organize as roll-cells aligned along the cell’s side walls, with rays of vorticity pointing toward the core of the cell. Two different large- scale organizations are observed and these patterns are well described spatially and energetically by azimuthal Fourier modes with frequencies of 2 and 3. These Fourier modes are shown to be dominant throughout the entire domain, and are found to be the primary source for radial inhomogeneity by inspection of the energy spectra. The precision with which the azimuthal Fourier modes describe these large-scale structures shows that these structures influence a large range of length scales. Conversely, the smaller scale structures are found to be more sensitive to radial position within the Fourier modes showing a strong dependence on physical length scales.
Dynamics in the large-scale structures are observed including a transition in the global pattern followed by a net rotation about the central axis. The transition takes place over 10 eddy-turnover times and the subsequent rotation occurs at a rate of approximately 1.1 degrees per eddy-turnover. These time-scales are of the same order of magnitude as those seen in lower aspect-ratio RBC for similar events and suggests a similarity in dynamic events across different aspect-ratios.
Large-scale structures are found to organize as roll-cells aligned along the cell’s side walls, with rays of vorticity pointing toward the core of the cell. Two different large- scale organizations are observed and these patterns are well described spatially and energetically by azimuthal Fourier modes with frequencies of 2 and 3. These Fourier modes are shown to be dominant throughout the entire domain, and are found to be the primary source for radial inhomogeneity by inspection of the energy spectra. The precision with which the azimuthal Fourier modes describe these large-scale structures shows that these structures influence a large range of length scales. Conversely, the smaller scale structures are found to be more sensitive to radial position within the Fourier modes showing a strong dependence on physical length scales.
Dynamics in the large-scale structures are observed including a transition in the global pattern followed by a net rotation about the central axis. The transition takes place over 10 eddy-turnover times and the subsequent rotation occurs at a rate of approximately 1.1 degrees per eddy-turnover. These time-scales are of the same order of magnitude as those seen in lower aspect-ratio RBC for similar events and suggests a similarity in dynamic events across different aspect-ratios.
ContributorsSakievich, Philip Sakievich (Author) / Peet, Yulia (Thesis advisor) / Adrian, Ronald (Committee member) / Squires, Kyle (Committee member) / Herrmann, Marcus (Committee member) / Kostelich, Eric (Committee member) / Arizona State University (Publisher)
Created2017
Description
Many defense, healthcare, and energy applications can benefit from the development of surfaces that easily shed droplets of liquids of interest. Desired wetting properties are typically achieved via altering the surface chemistry or topography or both through surface engineering. Despite many recent advancements, materials modified only on their exterior are still prone to physical degradation and lack durability. In contrast to surface engineering, this thesis focuses on altering the bulk composition and the interior of a material to tune how an exterior surface would interact with liquids. Fundamental and applied aspects of engineering of two material systems with low contact angle hysteresis (i.e. ability to easily shed droplets) are explained. First, water-shedding metal matrix hydrophobic nanoparticle composites with high thermal conductivity for steam condensation rate enhancement are discussed. Despite having static contact angle <90° (not hydrophobic), sustained dropwise steam condensation can be achieved at the exterior surface of the composite due to low contact angle hysteresis (CAH). In order to explain this observation, the effect of varying the length scale of surface wetting heterogeneity over three orders of magnitude on the value of CAH was experimentally investigated. This study revealed that the CAH value is primarily governed by the pinning length which in turn depends on the length scale of wetting heterogeneity. Modifying the heterogeneity size ultimately leads to near isotropic wettability for surfaces with highly anisotropic nanoscale chemical heterogeneities. Next, development of lubricant-swollen polymeric omniphobic protective gear for defense and healthcare applications is described. Specifically, it is shown that the robust and durable protective gear can be made from polymeric material fully saturated with lubricant that can shed all liquids irrespective of their surface tensions even after multiple contact incidences with the foreign objects. Further, a couple of schemes are proposed to improve the rate of lubrication and replenishment of lubricant as well as reduce the total amount of lubricant required in making the polymeric protective gear omniphobic. Overall, this research aims to understand the underlying physics of dynamic surface-liquid interaction and provides simple scalable route to fabricate better materials for condensers and omniphobic protective gear.
ContributorsDamle, Viraj (Author) / Rykaczewski, Konrad (Thesis advisor) / Phelan, Patrick (Committee member) / Lin, Jerry (Committee member) / Herrmann, Marcus (Committee member) / Wang, Robert (Committee member) / Wang, Liping (Committee member) / Arizona State University (Publisher)
Created2017
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
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
Description
Owing to the surge in development of endovascular devices such as coils and flow diverter stents, doctors are inclined to approach surgical cases non-invasively more often than before. Treating brain aneurysms as a bulging of a weakened area of a blood vessel is no exception. Therefore, promoting techniques that can help surgeons have a better idea of treatment outcomes are of invaluable importance.
In order to investigate the effects of these devices on intra-aneurysmal hemodynamics, the conventional computational fluid dynamics (CFD) approach uses the explicit geometry of the device within an aneurysm and discretizes the fluid domain to solve the Navier-Stokes equations. However, since the devices are made of small struts, the number of mesh elements in the boundary layer region would be considerable. This cumbersome task led to the implementation of the porous medium assumption. In this approach, the explicit geometry of the device is eliminated, and relevant porous medium assumptions are applied. Unfortunately, as it will be shown in this research, some of the porous medium approaches used in the literature are over-simplified. For example, considering the porous domain to be homogeneous is one major drawback which leads to significant errors in capturing the intra-aneurysmal flow features. Specifically, since the devices must comply with the complex geometry of an aneurysm, the homogeneity assumption is not valid.
In this research, a novel heterogeneous porous medium approach is introduced. This results in a substantial reduction in the total number of mesh elements required to discretize the flow domain while not sacrificing the accuracy of the method by over-simplifying the utilized assumptions.
In order to investigate the effects of these devices on intra-aneurysmal hemodynamics, the conventional computational fluid dynamics (CFD) approach uses the explicit geometry of the device within an aneurysm and discretizes the fluid domain to solve the Navier-Stokes equations. However, since the devices are made of small struts, the number of mesh elements in the boundary layer region would be considerable. This cumbersome task led to the implementation of the porous medium assumption. In this approach, the explicit geometry of the device is eliminated, and relevant porous medium assumptions are applied. Unfortunately, as it will be shown in this research, some of the porous medium approaches used in the literature are over-simplified. For example, considering the porous domain to be homogeneous is one major drawback which leads to significant errors in capturing the intra-aneurysmal flow features. Specifically, since the devices must comply with the complex geometry of an aneurysm, the homogeneity assumption is not valid.
In this research, a novel heterogeneous porous medium approach is introduced. This results in a substantial reduction in the total number of mesh elements required to discretize the flow domain while not sacrificing the accuracy of the method by over-simplifying the utilized assumptions.
ContributorsYadollahi Farsani, Hooman (Author) / Herrmann, Marcus (Thesis advisor) / Frakes, David (Thesis advisor) / Chong, Brian (Committee member) / Peet, Yulia (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2018
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 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