Matching Items (43)
Filtering by
- All Subjects: Mechanical Engineering
- Creators: Herrmann, Marcus
Description
Recent literature indicates potential benefits in microchannel cooling if an inlet orifice is used to suppress pressure oscillations that develop under two-phase conditions. This study investigates the costs and benefits of using an adjustable microchannel inlet orifice. The focus is on orifice effect during steady-state boiling and critical heat flux (CHF) in the channels using R134a in a pumped refrigerant loop (PRL). To change orifice size, a dam controlled with a micrometer was placed in front of 31 parallel microchannels. Each channel had a hydraulic diameter of 0.235 mm and a length of 1.33 cm. For steady state two-phase conditions, mass fluxes of 300 kg m-2 s-1 and 600 kg m-2 s-1were investigated. For orifice sizes with a hydraulic diameter to unrestricted hydraulic diameter (Dh:Dh,ur) ratio less than 35 percent, oscillations were reduced and wall temperatures fell up to 1.5 °C. Critical heat flux data were obtained for 7 orifice sizes with mass fluxes from 186 kg m-2 s-1 to 847 kg m-2 s-1. For all mass fluxes and inlet conditions tested, CHF values for a Dh:Dh,ur ratio of 1.8 percent became increasingly lower (up to 37 W cm-2 less) than those obtained with larger orifices. An optimum orifice size with Dh:Dh,ur of 35 percent emerged, offering up to 5 W cm-2 increase in CHF over unrestricted conditions at the highest mass flux tested, 847 kg m-2 s-1. These improvements in cooling ability with inlet orifices in place under both steady-state and impending CHF conditions are modest, leading to the conclusion that inlet orifices are only mildly effective at improving heat transfer coefficients. Stability of the PRL used for experimentation was also studied and improved. A vapor compression cycle's (VCC) proportional, integral, and derivative controller was found to adversely affect stability within the PRL and cause premature CHF. Replacing the VCC with an ice water heat sink maintained steady pumped loop system pressures and mass flow rates. The ice water heat sink was shown to have energy cost savings over the use of a directly coupled VCC for removing heat from the PRL.
ContributorsOdom, Brent A (Author) / Phelan, Patrick E (Thesis advisor) / Herrmann, Marcus (Committee member) / Trimble, Steve (Committee member) / Tasooji, Amaneh (Committee member) / Holcomb, Don (Committee member) / Arizona State University (Publisher)
Created2012
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
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
Multiphase flows are relevant to various industrial processes and are also a ubiquitous feature of nature. Atomization is a Gas-Liquid class of multiphase flow in which the liquid bulk disintegrates into a spectrum of drops. The final drop size distribution of fragmenting liquids is important and is crucial to quantifying the performance of atomizers. This thesis implements two models of ligament breakup. The first model provides a method to determine the droplet size distribution of fragmenting ligaments. The second model provides a relation between ligament stretching, aspect ratio and dimensionless properties like Ohnesorge and Weber numbers for ligaments being stretched by aerodynamic force. The first model by Villermaux et.al considers a ligament as a linear succession of liquid blobs which undergo continuous interplay during destabilization. The evolution of their size distribution ultimately rules the droplet size distribution which follow a gamma distribution [14]. The results show that the Direct Numerical Simulations (DNS) of ligaments with different perturbations fragmented into very few drops and cannot be used to confirm that they follow the predicted gamma distribution. The second model considers a ligament breakup due to Rayleigh-Plateau Instability and provides an equation for ligament stretching. Through test runs the proportionality constant in the equation is determined by a least square fit. The theoretical number of drops is compared with the number of drops resulting from the Direct Numerical Simulation of ligament with a sinusoidal perturbation. It is found that the wavelength of the initial perturbation does not determine the number of drops obtained by ligament breakup
ContributorsRama Krishna, Prathyush (Author) / Herrmann, Marcus (Thesis advisor) / Takahashi, Timothy (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2021
Description
Spray flows are important in a myriad of practical applications including fuel injection, ink-jet printing, agricultural sprays, and industrial processes. Two-phase sprays find particular use for spot cooling applications with high heat fluxes as in casting processes and power electronics. Computability of sprays in a cost-effective manner provides a path to optimize the design of nozzles to tune the spray characteristics for the needs of a particular application. Significant research has so far been devoted to understand and characterize spray flows better, be it from a theoretical, experimental or computational standpoint. The current thesis discusses a methodology for modeling primary atomization using the Quadratic Formula which is derived from an integral formulation of the governing equations. The framework is then applied to different examples of flat-fan hydraulic sprays. For each case, the spray is first resolved as a continuous fluid using the volume of fluid method. Atomization criterion is then applied to the velocity flow-field to determine the sites for primary atomization. At each site, local diameters for particle injection is determined using the quadratic formula. The trajectory of injected particles are then monitored through a particle tracking algorithm. The results from the numerical analysis are compared with experimental data to validate the computational framework.
ContributorsBhardwaj, Angshuman (Author) / Lee, T.-W. (Thesis advisor) / Herrmann, Marcus (Committee member) / Rykaczewski, Konrad (Committee member) / Arizona State University (Publisher)
Created2022
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
Identifying and tracking the location of the fluid interface is a fundamental aspect of multiphase flows. The Volume of Fluid (VOF) and Level Set methods are widely used to track the interface accurately. Analyzing the liquid structures such as sheets, ligaments, and droplets helps understand the flow physics and fluid breakup mechanism, aids in predicting droplet formation, improves atomization modeling and spray combustion. The thesis focuses on developing a new method to identify these liquid structures and devise a sphere model for droplet size prediction by augmenting concepts of linear algebra, rigid body dynamics, computational fluid mechanics, scientific computing, and visualization. The first part of the thesis presents a new approach to classify the fluid structures based on their length scales along their principal axes. This approach provides a smooth tracking of the structures' generation history instead of relying on high-speed video imaging of the experiment. A droplet is observed to have three equal length scales, while a ligament has one and a sheet has two significantly larger length scales. The subsequent breakup of ligaments and droplets depends on the atomizer geometry, operating conditions, and fluid physical properties. While it's straightforward to apply DNS and estimate this breakup, it is proven to be computationally expensive. The second part of the thesis deals with developing a sphere model that would essentially reduce this computational cost. After identifying a liquid structure, the sphere model utilizes the level set data in the domain to quantify the structure using spheres. By using the evolution information of these spheres as they separate from each other, the subsequent droplet size distribution can be evaluated.
ContributorsKashetty, Sindhuja (Author) / Herrmann, Marcus (Thesis advisor) / Wells, Valana (Committee member) / Kim, Jeonglae (Committee member) / Arizona State University (Publisher)
Created2021
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
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
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