Matching Items (17)
Filtering by
- All Subjects: Mechanical Engineering
- Genre: Doctoral Dissertation
- Creators: Herrmann, Marcus
- Resource Type: Text
Description
The atomization of a liquid jet by a high speed cross-flowing gas has many applications such as gas turbines and augmentors. The mechanisms by which the liquid jet initially breaks up, however, are not well understood. Experimental studies suggest the dependence of spray properties on operating conditions and nozzle geom- etry. Detailed numerical simulations can offer better understanding of the underlying physical mechanisms that lead to the breakup of the injected liquid jet. In this work, detailed numerical simulation results of turbulent liquid jets injected into turbulent gaseous cross flows for different density ratios is presented. A finite volume, balanced force fractional step flow solver to solve the Navier-Stokes equations is employed and coupled to a Refined Level Set Grid method to follow the phase interface. To enable the simulation of atomization of high density ratio fluids, we ensure discrete consistency between the solution of the conservative momentum equation and the level set based continuity equation by employing the Consistent Rescaled Momentum Transport (CRMT) method. The impact of different inflow jet boundary conditions on different jet properties including jet penetration is analyzed and results are compared to those obtained experimentally by Brown & McDonell(2006). In addition, instability analysis is performed to find the most dominant insta- bility mechanism that causes the liquid jet to breakup. Linear instability analysis is achieved using linear theories for Rayleigh-Taylor and Kelvin- Helmholtz instabilities and non-linear analysis is performed using our flow solver with different inflow jet boundary conditions.
ContributorsGhods, Sina (Author) / Herrmann, Marcus (Thesis advisor) / Squires, Kyle (Committee member) / Chen, Kangping (Committee member) / Huang, Huei-Ping (Committee member) / Tang, Wenbo (Committee member) / Arizona State University (Publisher)
Created2013
Description
The heat transfer enhancements available from expanding the cross-section of a boiling microchannel are explored analytically and experimentally. Evaluation of the literature on critical heat flux in flow boiling and associated pressure drop behavior is presented with predictive critical heat flux (CHF) and pressure drop correlations. An optimum channel configuration allowing maximum CHF while reducing pressure drop is sought. A perturbation of the channel diameter is employed to examine CHF and pressure drop relationships from the literature with the aim of identifying those adequately general and suitable for use in a scenario with an expanding channel. Several CHF criteria are identified which predict an optimizable channel expansion, though many do not. Pressure drop relationships admit improvement with expansion, and no optimum presents itself. The relevant physical phenomena surrounding flow boiling pressure drop are considered, and a balance of dimensionless numbers is presented that may be of qualitative use. The design, fabrication, inspection, and experimental evaluation of four copper microchannel arrays of different channel expansion rates with R-134a refrigerant is presented. Optimum rates of expansion which maximize the critical heat flux are considered at multiple flow rates, and experimental results are presented demonstrating optima. The effect of expansion on the boiling number is considered, and experiments demonstrate that expansion produces a notable increase in the boiling number in the region explored, though no optima are observed. Significant decrease in the pressure drop across the evaporator is observed with the expanding channels, and no optima appear. Discussion of the significance of this finding is presented, along with possible avenues for future work.
ContributorsMiner, Mark (Author) / Phelan, Patrick E (Thesis advisor) / Baer, Steven (Committee member) / Chamberlin, Ralph (Committee member) / Chen, Kangping (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2013
Description
A new method of adaptive mesh generation for the computation of fluid flows is investigated. The method utilizes gradients of the flow solution to adapt the size and stretching of elements or volumes in the computational mesh as is commonly done in the conventional Hessian approach. However, in the new method, higher-order gradients are used in place of the Hessian. The method is applied to the finite element solution of the incompressible Navier-Stokes equations on model problems. Results indicate that a significant efficiency benefit is realized.
ContributorsShortridge, Randall (Author) / Chen, Kang Ping (Thesis advisor) / Herrmann, Marcus (Thesis advisor) / Wells, Valana (Committee member) / Huang, Huei-Ping (Committee member) / Mittelmann, Hans (Committee member) / Arizona State University (Publisher)
Created2011
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
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
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
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