Matching Items (44)
Filtering by
- All Subjects: Mechanical Engineering
- Creators: Herrmann, Marcus
- Member of: Theses and Dissertations
Numerical Modeling of Cake Formation and Permeate Flux Decline in Membrane Filtration Using OpenFOAM
Description
Membrane filtration is an important technology in industry. In the past few decades, equations have been developed from experimental results to predict cake formation and permeate flux decline in the membrane filtration process. In the current work, the block of particles on membrane surface is achieved by setting surface flux on membrane surface zero. This approach is implemented for both microfiltration and nanofiltration using OpenFOAM. Moreover, a new method to deal with cake resistance for nanofiltration is introduced. Cake resistance is applied to both cake and membrane. To validate the new techniques, results of crossflow microfiltration are compared to theoretical results and results of two crossflow nanofiltration cases are compared to experimental data. In addition, the new techniques are applied to dead end filtration to observe the different structure of the cake and explore the effect of resistance on velocity profile.
ContributorsHu, Jueming (Author) / Herrmann, Marcus (Thesis advisor) / Huang, Huei-Ping (Committee member) / Lee, Taewoo (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
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 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.
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
Description
The traditional understanding of robotics includes mechanisms of rigid structures, which can manipulate surrounding objects, taking advantage of mechanical actuators such as motors and servomechanisms. Although these methods provide the underlying fundamental concepts behind much of modern technological infrastructure, in fields such as manufacturing, automation, and biomedical application, the robotic structures formed by rigid axels on mechanical actuators lack the delicate differential sensors and actuators associated with known biological systems. The rigid structures of traditional robotics also inhibit the use of simple mechanisms in congested and/or fragile environments. By observing a variety of biological systems, it is shown that nature models its structures over millions of years of evolution into a combination of soft structures and rigid skeletal interior supports. Through technological bio-inspired designs, researchers hope to mimic some of the complex behaviors of biological mechanisms using pneumatic actuators coupled with highly compliant materials that exhibit relatively large reversible elastic strain. This paper begins the brief history of soft robotics, the various classifications of pneumatic fluid systems, the associated difficulties that arise with the unpredictable nature of fluid reactions, the methods of pneumatic actuators in use today, the current industrial applications of soft robotics, and focuses in large on the construction of a universally adaptable soft robotic gripper and material application tool. The central objective of this experiment is to compatibly pair traditional rigid robotics with the emerging technologies of sort robotic actuators. This will be done by combining a traditional rigid robotic arm with a soft robotic manipulator bladder for the purposes of object manipulation and excavation of extreme environments.
ContributorsShuster, Eden S. (Author) / Thanga, Jekan (Thesis director) / Asphaug, Erik (Committee member) / Mechanical and Aerospace Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Numerical simulation of environmental flow over urban landscape for applications to renewable energy
Description
Development of renewable energy solutions has become a major interest among environmental organizations and governments around the world due to an increase in energy consumption and global warming. One fast growing renewable energy solution is the application of wind energy in cities. To qualitative and quantitative predict wind turbine performance in urban areas, CFD simulation is performed on real-life urban geometry and wind velocity profiles are evaluated. Two geometries in Arizona is selected in this thesis to demonstrate the influence of building heights; one of the simulation models, ASU campus, is relatively low rise and without significant tall buildings; the other model, the downtown phoenix model, are high-rise and with greater building height difference. The content of this thesis focuses on using RANS computational fluid dynamics approach to simulate wind acceleration phenomenon in two complex geometries, ASU campus and Phoenix downtown model. Additionally, acceleration ratio and locations are predicted, the results are then used to calculate the best location for small wind turbine installments.
ContributorsYing, Xiaoyan (Author) / Huang, Huei-Ping (Thesis advisor) / Peet, Yulia (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2015
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