Matching Items (16)
Filtering by
- All Subjects: Mechanical Engineering
- All Subjects: Computational Fluid Dynamics
- Genre: Doctoral Dissertation
- Creators: Chen, Kangping
Description
Microfluidics is the study of fluid flow at very small scales (micro -- one millionth of a meter) and is prevalent in many areas of science and engineering. Typical applications include lab-on-a-chip devices, microfluidic fuel cells, and DNA separation technologies. Many of these microfluidic devices rely on micron-resolution velocimetry measurements to improve microchannel design and characterize existing devices. Methods such as micro particle imaging velocimetry (microPIV) and micro particle tracking velocimetry (microPTV) are mature and established methods for characterization of steady 2D flow fields. Increasingly complex microdevices require techniques that measure unsteady and/or three dimensional velocity fields. This dissertation presents a method for three-dimensional velocimetry of unsteady microflows based on spinning disk confocal microscopy and depth scanning of a microvolume. High-speed 2D unsteady velocity fields are resolved by acquiring images of particle motion using a high-speed CMOS camera and confocal microscope. The confocal microscope spatially filters out of focus light using a rotating disk of pinholes placed in the imaging path, improving the ability of the system to resolve unsteady microPIV measurements by improving the image and correlation signal to noise ratio. For 3D3C measurements, a piezo-actuated objective positioner quickly scans the depth of the microvolume and collects 2D image slices, which are stacked into 3D images. Super resolution microPIV interrogates these 3D images using microPIV as a predictor field for tracking individual particles with microPTV. The 3D3C diagnostic is demonstrated by measuring a pressure driven flow in a three-dimensional expanding microchannel. The experimental velocimetry data acquired at 30 Hz with instantaneous spatial resolution of 4.5 by 4.5 by 4.5 microns agrees well with a computational model of the flow field. The technique allows for isosurface visualization of time resolved 3D3C particle motion and high spatial resolution velocity measurements without requiring a calibration step or reconstruction algorithms. Several applications are investigated, including 3D quantitative fluorescence imaging of isotachophoresis plugs advecting through a microchannel and the dynamics of reaction induced colloidal crystal deposition.
ContributorsKlein, Steven Adam (Author) / Posner, Jonathan D (Thesis advisor) / Adrian, Ronald (Committee member) / Chen, Kangping (Committee member) / Devasenathipathy, Shankar (Committee member) / Frakes, David (Committee member) / Arizona State University (Publisher)
Created2011
Description
Locomotion of microorganisms is commonly observed in nature. Although microorganism locomotion is commonly attributed to mechanical deformation of solid appendages, in 1956 Nobel Laureate Peter Mitchell proposed that an asymmetric ion flux on a bacterium's surface could generate electric fields that drive locomotion via self-electrophoresis. Recent advances in nanofabrication have enabled the engineering of synthetic analogues, bimetallic colloidal particles, that swim due to asymmetric ion flux originally proposed by Mitchell. Bimetallic colloidal particles swim through aqueous solutions by converting chemical fuel to fluid motion through asymmetric electrochemical reactions. This dissertation presents novel bimetallic motor fabrication strategies, motor functionality, and a study of the motor collective behavior in chemical concentration gradients. Brownian dynamics simulations and experiments show that the motors exhibit chemokinesis, a motile response to chemical gradients that results in net migration and concentration of particles. Chemokinesis is typically observed in living organisms and distinct from chemotaxis in that there is no particle directional sensing. The synthetic motor chemokinesis observed in this work is due to variation in the motor's velocity and effective diffusivity as a function of the fuel and salt concentration. Static concentration fields are generated in microfluidic devices fabricated with porous walls. The development of nanoscale particles that swim autonomously and collectively in chemical concentration gradients can be leveraged for a wide range of applications such as directed drug delivery, self-healing materials, and environmental remediation.
ContributorsWheat, Philip Matthew (Author) / Posner, Jonathan D (Thesis advisor) / Phelan, Patrick (Committee member) / Chen, Kangping (Committee member) / Buttry, Daniel (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
Wind measurements are fundamental inputs for the evaluation of potential energy yield and performance of wind farms. Three-dimensional scanning coherent Doppler lidar (CDL) may provide a new basis for wind farm site selection, design, and control. In this research, CDL measurements obtained from multiple wind energy developments are analyzed and a novel wind farm control approach has been modeled. The possibility of using lidar measurements to more fully characterize the wind field is discussed, specifically, terrain effects, spatial variation of winds, power density, and the effect of shear at different layers within the rotor swept area. Various vector retrieval methods have been applied to the lidar data, and results are presented on an elevated terrain-following surface at hub height. The vector retrieval estimates are compared with tower measurements, after interpolation to the appropriate level. CDL data is used to estimate the spatial power density at hub height. Since CDL can measure winds at different vertical levels, an approach for estimating wind power density over the wind turbine rotor-swept area is explored. Sample optimized layouts of wind farm using lidar data and global optimization algorithms, accounting for wake interaction effects, have been explored. An approach to evaluate spatial wind speed and direction estimates from a standard nested Coupled Ocean and Atmosphere Mesoscale Prediction System (COAMPS) model and CDL is presented. The magnitude of spatial difference between observations and simulation for wind energy assessment is researched. Diurnal effects and ramp events as estimated by CDL and COAMPS were inter-compared. Novel wind farm control based on incoming winds and direction input from CDL's is developed. Both yaw and pitch control using scanning CDL for efficient wind farm control is analyzed. The wind farm control optimizes power production and reduces loads on wind turbines for various lidar wind speed and direction inputs, accounting for wind farm wake losses and wind speed evolution. Several wind farm control configurations were developed, for enhanced integrability into the electrical grid. Finally, the value proposition of CDL for a wind farm development, based on uncertainty reduction and return of investment is analyzed.
ContributorsKrishnamurthy, Raghavendra (Author) / Calhoun, Ronald J (Thesis advisor) / Chen, Kangping (Committee member) / Huang, Huei-Ping (Committee member) / Fraser, Matthew (Committee member) / Phelan, Patrick (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
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
Portable devices rely on battery systems that contribute largely to the overall device form factor and delay portability due to recharging. Membraneless microfluidic fuel cells are considered as the next generation of portable power sources for their compatibility with higher energy density reactants. Microfluidic fuel cells are potentially cost effective and robust because they use low Reynolds number flow to maintain fuel and oxidant separation instead of ion exchange membranes. However, membraneless fuel cells suffer from poor efficiency due to poor mass transport and Ohmic losses. Current microfluidic fuel cell designs suffer from reactant cross-diffusion and thick boundary layers at the electrode surfaces, which result in a compromise between the cell's power output and fuel utilization. This dissertation presents novel flow field architectures aimed at alleviating the mass transport limitations. The first architecture provides a reactant interface where the reactant diffusive concentration gradients are aligned with the bulk flow, mitigating reactant mixing through diffusion and thus crossover. This cell also uses porous electro-catalysts to improve electrode mass transport which results in higher extraction of reactant energy. The second architecture uses porous electrodes and an inert conductive electrolyte stream between the reactants to enhance the interfacial electrical conductivity and maintain complete reactant separation. This design is stacked hydrodynamically and electrically, analogous to membrane based systems, providing increased reactant utilization and power. These fuel cell architectures decouple the fuel cell's power output from its fuel utilization. The fuel cells are tested over a wide range of conditions including variation of the loads, reactant concentrations, background electrolytes, flow rates, and fuel cell geometries. These experiments show that increasing the fuel cell power output is accomplished by increasing reactant flow rates, electrolyte conductivity, and ionic exchange areas, and by decreasing the spacing between the electrodes. The experimental and theoretical observations presented in this dissertation will aid in the future design and commercialization of a new portable power source, which has the desired attributes of high power output per weight and volume and instant rechargeability.
ContributorsSalloum, Kamil S (Author) / Posner, Jonathan D (Thesis advisor) / Adrian, Ronald (Committee member) / Christen, Jennifer (Committee member) / Phelan, Patrick (Committee member) / Chen, Kangping (Committee member) / Arizona State University (Publisher)
Created2010
Description
As miniature and high-heat-dissipation equipment became major manufacture and operation trends, heat-rejecting and heat-transport solutions faced increasing challenges. In the 1970s, researchers showed that particle suspensions can enhance the heat transfer efficiency of their base fluids. However, their work was hindered by the sedimentation and erosion issues caused by the relatively large particle sizes in their suspensions. More recently, nanofluids--suspensions of nanoparticles in liquids-were proposed to be applied as heat transfer fluids, because of the enhanced thermal conductivity that has generally been observed. However, in practical applications, a heat conduction mechanism may not be sufficient for cooling high-heat-dissipation devices such as microelectronics or powerful optical equipment. Thus, the thermal performance under convective, i.e., flowing heat transfer conditions becomes of primary interest. In addition, with the presence of nanoparticles, the viscosity of a nanofluid is greater than its base fluid and deviates from Einstein's classical prediction. Through the use of a test rig designed and assembled as part of this dissertation, the viscosity and heat transfer coefficient of nanofluids can be simultaneously determined by pressure drop and temperature difference measurements under laminar flow conditions. An extensive characterization of the nanofluid samples, including pH, electrical conductivity, particle sizing and zeta potential, is also documented. Results indicate that with constant wall heat flux, the relative viscosities of nanofluid decrease with increasing volume flow rate. The results also show, based on Brenner's model, that the nanofluid viscosity can be explained in part by the aspect ratio of the aggregates. The measured heat transfer coefficient values for nanofluids are generally higher than those for base fluids. In the developing region, this can be at least partially explained by Prandtl number effects. The Nusselt number ( Nu ) results for nanofluid show that Nu increases with increasing nanofluid volume fraction and volume flow rate. However, only DI-H2O (deionized water) and 5/95 PG/H2O (PG = propylene glycol) based nanofluids with 1 vol% nanoparticle loading have Nu greater than the theoretical prediction, 4.364. It is suggested that the nanofluid has potential to be applied within the thermally developing region when utilizing the nanofluid as a heat transfer liquid in a circular tube. The suggested Reynold's number is greater than 100.
ContributorsLai, Wei-Yun (Author) / Phelan, Patrick E (Thesis advisor) / Chen, Kangping (Committee member) / Hayes, Mark (Committee member) / Prasher, Ravi S (Committee member) / Sieradzki, Karl (Committee member) / Arizona State University (Publisher)
Created2010
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
This dissertation introduces FARCOM (Fortran Adaptive Refiner for Cartesian Orthogonal Meshes), a new general library for adaptive mesh refinement (AMR) based on an unstructured hexahedral mesh framework. As a result of the underlying unstructured formulation, the refinement and coarsening operators of the library operate on a single-cell basis and perform in-situ replacement of old mesh elements. This approach allows for h-refinement without the memory and computational expense of calculating masked coarse grid cells, as is done in traditional patch-based AMR approaches, and enables unstructured flow solvers to have access to the automated domain generation capabilities usually only found in tree AMR formulations.
The library is written to let the user determine where to refine and coarsen through custom refinement selector functions for static mesh generation and dynamic mesh refinement, and can handle smooth fields (such as level sets) or localized markers (e.g. density gradients). The library was parallelized with the use of the Zoltan graph-partitioning library, which provides interfaces to both a graph partitioner (PT-Scotch) and a partitioner based on Hilbert space-filling curves. The partitioned adjacency graph, mesh data, and solution variable data is then packed and distributed across all MPI ranks in the simulation, which then regenerate the mesh, generate domain decomposition ghost cells, and create communication caches.
Scalability runs were performed using a Leveque wave propagation scheme for solving the Euler equations. The results of simulations on up to 1536 cores indicate that the parallel performance is highly dependent on the graph partitioner being used, and differences between the partitioners were analyzed. FARCOM is found to have better performance if each MPI rank has more than 60,000 cells.
The library is written to let the user determine where to refine and coarsen through custom refinement selector functions for static mesh generation and dynamic mesh refinement, and can handle smooth fields (such as level sets) or localized markers (e.g. density gradients). The library was parallelized with the use of the Zoltan graph-partitioning library, which provides interfaces to both a graph partitioner (PT-Scotch) and a partitioner based on Hilbert space-filling curves. The partitioned adjacency graph, mesh data, and solution variable data is then packed and distributed across all MPI ranks in the simulation, which then regenerate the mesh, generate domain decomposition ghost cells, and create communication caches.
Scalability runs were performed using a Leveque wave propagation scheme for solving the Euler equations. The results of simulations on up to 1536 cores indicate that the parallel performance is highly dependent on the graph partitioner being used, and differences between the partitioners were analyzed. FARCOM is found to have better performance if each MPI rank has more than 60,000 cells.
ContributorsBallesteros, Carlos Alberto (Author) / Herrmann, Marcus (Thesis advisor) / Adrian, Ronald (Committee member) / Chen, Kangping (Committee member) / Huang, Huei-Ping (Committee member) / Lopez, Juan (Committee member) / Arizona State University (Publisher)
Created2019