Two-Phase Flow In Microchannels: Morphology And Interface Phenomena

The existence and morphology, as well as the dynamics of micro-scale gas-liquid interfaces is investigated numerically and experimentally. These studies can be used to assess liquid management issues in microsystems such as PEMFC gas flow channels, and are meant to open new research perspectives in two-phase flow, particularly in film deposition on non-wetting surfaces. For example the critical plug volume data can be used to deliver desired length plugs, or to determine the plug formation frequency. The dynamics of gas-liquid interfaces, of interest for applications involving small passages (e.g. heat exchangers, phase separators and filtration systems), was investigated using high-speed microscopy - a method that also proved useful for the study of film deposition processes. The existence limit for a liquid plug forming in a mixed wetting channel is determined by numerical simulations using Surface Evolver. The plug model simulate actual conditions in the gas flow channels of PEM fuel cells, the wetting of the gas diffusion layer (GDL) side of the channel being different from the wetting of the bipolar plate walls. The minimum plug volume, denoted as critical volume is computed for a series of GDL and bipolar plate wetting properties. Critical volume data is meant to assist in the water management of PEMFC, when corroborated with experimental data. The effect of cross section geometry is assessed by computing the critical volume in square and trapezoidal channels. Droplet simulations show that water can be passively removed from the GDL surface towards the bipolar plate if we take advantage on differing wetting properties between the two surfaces, to possibly avoid the gas transport blockage through the GDL. High speed microscopy was employed in two-phase and film deposition experiments with water in round and square capillary tubes. Periodic interface destabilization was observed and the existence of compression waves in the gas phase is discussed by taking into consideration a naturally occurring convergent-divergent nozzle formed by the flowing liquid phase. The effect of channel geometry and wetting properties was investigated through two-phase water-air flow in square and round microchannels, having three static contact angles of 20, 80 and 105 degrees. Four different flow regimes are observed for a fixed flow rate, this being thought to be caused by the wetting behavior of liquid flowing in the corners as well as the liquid film stability. Film deposition experiments in wetting and non-wetting round microchannels show that a thicker film is deposited for wetting conditions departing from the ideal 0 degrees contact angle. A film thickness dependence with the contact angle theta as well as the Capillary number, in the form h_R ~ Ca^(2/3)/ cos(theta) is inferred from scaling arguments, for contact angles smaller than 36 degrees. Non-wetting film deposition experiments reveal that a film significantly thicker than the wetting Bretherton film is deposited. A hydraulic jump occurs if critical conditions are met, as given by a proposed nondimensional parameter similar to the Froude number. Film thickness correlations are also found by matching the measured and the proposed velocity derived in the shock theory. The surface wetting as well as the presence of the shock cause morphological changes in the Taylor bubble flow.

Diesel spray mixing limited vaporization with non-ideal and multi-component fuel thermophysical property effects

Internal combustion engines are, and will continue to be, a primary mode of power generation for ground transportation. Challenges exist in meeting fuel consumption regulations and emission standards while upholding performance, as fuel prices rise, and resource depletion and environmental impacts are of increasing concern. Diesel engines are advantageous due to their inherent efficiency advantage over spark ignition engines; however, their NOx and soot emissions can be difficult to control and reduce due to an inherent tradeoff. Diesel combustion is spray and mixing controlled providing an intrinsic link between spray and emissions, motivating detailed, fundamental studies on spray, vaporization, mixing, and combustion characteristics under engine relevant conditions. An optical combustion vessel facility has been developed at Michigan Technological University for these studies, with detailed tests and analysis being conducted. In this combustion vessel facility a preburn procedure for thermodynamic state generation is used, and validated using chemical kinetics modeling both for the MTU vessel, and institutions comprising the Engine Combustion Network international collaborative research initiative. It is shown that minor species produced are representative of modern diesel engines running exhaust gas recirculation and do not impact the autoignition of n-heptane. Diesel spray testing of a high-pressure (2000 bar) multi-hole injector is undertaken including non-vaporizing, vaporizing, and combusting tests, with sprays characterized using Mie back scatter imaging diagnostics. Liquid phase spray parameter trends agree with literature. Fluctuations in liquid length about a quasi-steady value are quantified, along with plume to plume variations. Hypotheses are developed for their causes including fuel pressure fluctuations, nozzle cavitation, internal injector flow and geometry, chamber temperature gradients, and turbulence. These are explored using a mixing limited vaporization model with an equation of state approach for thermopyhysical properties. This model is also applied to single and multi-component surrogates. Results include the development of the combustion research facility and validated thermodynamic state generation procedure. The developed equation of state approach provides application for improving surrogate fuels, both single and multi-component, in terms of diesel spray liquid length, with knowledge of only critical fuel properties. Experimental studies are coupled with modeling incorporating improved thermodynamic non-ideal gas and fuel

COMPUTATIONAL STUDY OF MICROSTRUCTURE-PROPERTYMECHANISM RELATIONS IN FERROIC COMPOSITES

Ferroic materials, as notable members of smart materials, have been widely used in applications that perform sensing, actuation and control. The macroscopic property change of ferroic materials may become remarkably large during ferroic phase transition, leading to the fact that the macroscopic properties can be tuned by carefully applying a suitable external field (electric, magnetic, stress). To obtain an enhancement in physical and/or mechanical properties, different kinds of ferroic composites have been fabricated. The properties of a ferroic composite are determined not only by the properties and relative amounts of the constituent phases, but also by the microstructure of individual phase such as the phase connectivity, phase size, shape and spatial arrangement. This dissertation mainly focuses on the computational study of microstructure – property – mechanism relations in two representative ferroic composites, i.e., two-phase particulate magnetoelectric (ME) composite and polymer matrix ferroelectric composite. The former is a great example of ferroic composite exhibiting a new property and functionality that neither of the constituent phases possesses individually. The latter well represents the kind of ferroic composites having property combinations that are better than the existing materials. Phase field modeling was employed as the computing tool, and the required models for ferroic composites were developed based on existing models for monolithic materials. Extensive computational simulations were performed to investigate the microstructure-property relations and the underlying mechanism in ferroic composites. In particulate, it is found that for ME composite 0-3 connectivity (isolated magnetostrictive phase) is necessary to exhibit ME effect, and small but finite electrical conductivity of isolated magnetic phase can beneficially enhance ME effect. It is revealed that longitudinal and transverse ME coefficients of isotropic 0-3 particulate composites can be effectively tailored by controlling magnetic domain structures without resort to anisotropic two-phase microstructures. Simulations also show that the macroscopic properties of the ferroelectricpolymer composites critically depend on the ferroelectric phase connectivity while are not sensitive to the sizes and internal grain structures of the ceramic particles. Texturing is found critical to exploit the paraelectric«ferroelectric phase transition and nonlinear polarization behavior in paraelectric polycrystal and its polymer matrix composite. Additionally, a Diffuse Interface Field model was developed to simulate packing and motion in liquid phase which is promising for studying the fabrication of particulatepolymer composites.

MAGNETO-PHOTONIC CRYSTALS FOR OPTICAL SENSING APPLICATIONS

Among the optical structures investigated for optical sensing purpose, a significant amount of research has been conducted on photonic crystal based sensors. A particular advantage of photonic crystal based sensors is that they show superior sensitivity for ultra-small volume sensing. In this study we investigate polarization changes in response to the changes in the cover index of magneto-optic active photonic band gap structures. One-dimensional photonic-band gap structures fabricated on iron garnet materials yield large polarization rotations at the band gap edges. The enhanced polarization effects serve as an excellent tool for chemical sensing showing high degree of sensitivity for photonic crystal cover refractive index changes. The one dimensional waveguide photonic crystals are fabricated on single-layer bismuth-substituted rare earth iron garnet films ((Bi, Y, Lu)3(Fe, Ga)5O12 ) grown by liquid phase epitaxy on gadolinium gallium garnet substrates. Band gaps have been observed where Bragg scattering conditions links forward-going fundamental waveguide modes to backscattered high-order waveguide modes. Large near-band-edge polarization rotations which increase progressively with backscattered-mode order have been experimentally demonstrated for multiple samples with different composition, film thickness and fabrication parameters. Experimental findings are supported by theoretical analysis of Bloch modes polarization states showing that large near stop-band edge rotations are induced by the magneto-photonic crystal. Theoretical and experimental analysis conducted on polarization rotation sensitivity to waveguide photonic crystal cover refractive index changes shows a monotonic enhancement of the rotation with cover index. The sensor is further developed for selective chemical sensing by employing Polypyrrole as the photonic crystal cover layer. Polypyrrole is one of the extensively studied conducting polymers for selective analyte detection. Successful detection of aqueous ammonia and methanol has been achieved with Polypyrrole deposited magneto-photonic crystals.

Dysprosium transport in Nd-Fe-B pellets

The addition of heavy rare earth (RE) elements to Nd2Fe14B based magnets to form (Nd,Dy)2Fe14B is known to increase the coercivity and high temperature performance required for hybrid vehicle electric motors and other extreme temperature applications. Attempts to conserve heavy rare earth elements for high temperature (RE)2Fe14B based magnets have led to the development of a grain boundary diffusion process for bulk magnets. This process relies on transport of a heavy rare earth, such as Dy, into a bulk Nd2Fe14B magnet along pores, a low volume fraction of eutectic liquid along grain boundary grain triple junctions and grain boundaries. This enriches the grain surfaces in Dy through the thickness of the bulk magnet, leading to larger increases coercivity with a smaller Dy concentration than can be achieved with homogeneous alloys. Attempts to carry out the same process during sintering require significant control of Dy transport efficiency. The macroscopic transport of Dy in Nd2.7Fe14B1.4 based powder packs is studied using a 'layered' pellet, where Nd2.7Fe14B1.4powder is an interlayer and Dy source as a center layer. The sintering of this layered pellet provided evidence for very large effective diffusion lengths aided by Dy rich liquid flow through connected porosity. Approaches to controlling Dy transportation include decreasing the liquid phase transport capability of the powder pack by increasing the melting point of the Dy source and the decreasing amount of RE rich liquid in the powder packs. The solid-liquid reaction is studied in which melt spun Nd2.7Fe14B1.4 ribbons are PVD coated with Dy-Fe eutectic composition and then thermally treated. The resulting microstructure from the reaction between Dy-Fe eutectic coating and Nd2.7Fe14B1.4 ribbon is interpreted as support for a proposed dissolution/reprecipitation process between solid and liquid phases. The estimate the diffusion coefficient and the effective diffusion length of Dy sources in Nd2.7Fe14B1.4 layered pellets and melt spun ribbons were obtained from the calculation of Fick's second law combined with EDS results from the experiment. The results indicate that the effective diffusion coefficient of Dy in the layered pellets is higher than the diffusion in ribbons due to its higher porosity than ribbons.

DIFFUSE-INTERFACE FIELD APPROACH TO MODELING SELF-ASSEMBLY OF HETEROGENEOUS COLLOIDAL SYSTEMS AND RELATED DIPOLE-DIPOLE INTERACTION PHENOMENA

Colloid self-assembly under external control is a new route to fabrication of advanced materials with novel microstructures and appealing functionalities. The kinetic processes of colloidal self-assembly have attracted great interests also because they are similar to many atomic level kinetic processes of materials. In the past decades, rapid technological progresses have been achieved on producing shape-anisotropic, patchy, core-shell structured particles and particles with electric/magnetic charges/dipoles, which greatly enriched the self-assembled structures. Multi-phase carrier liquids offer new route to controlling colloidal self-assembly. Therefore, heterogeneity is the essential characteristics of colloid system, while so far there still lacks a model that is able to efficiently incorporate these possible heterogeneities. This thesis is mainly devoted to development of a model and computational study on the complex colloid system through a diffuse-interface field approach (DIFA), recently developed by Wang et al. This meso-scale model is able to describe arbitrary particle shape and arbitrary charge/dipole distribution on the surface or body of particles. Within the framework of DIFA, a Gibbs-Duhem-type formula is introduced to treat Laplace pressure in multi-liquid-phase colloidal system and it obeys Young-Laplace equation. The model is thus capable to quantitatively study important capillarity related phenomena. Extensive computer simulations are performed to study the fundamental behavior of heterogeneous colloidal system. The role of Laplace pressure is revealed in determining the mechanical equilibrium of shape-anisotropic particles at fluid interfaces. In particular, it is found that the Laplace pressure plays a critical role in maintaining the stability of capillary bridges between close particles, which sheds light on a novel route to in situ firming compact but fragile colloidal microstructures via capillary bridges. Simulation results also show that competition between like-charge repulsion, dipole-dipole interaction and Brownian motion dictates the degree of aggregation of heterogeneously charged particles. Assembly and alignment of particles with magnetic dipoles under external field is studied. Finally, extended studies on the role of dipole-dipole interaction are performed for ferromagnetic and ferroelectric domain phenomena. The results reveal that the internal field generated by dipoles competes with external field to determine the dipole-domain evolution in ferroic materials.

DEVELOPMENT OF CONFOCAL IMAGING TECHNIQUES FOR PROBING INTERFACIAL DYNAMICS IN MICROSCALE, GAS-LIQUID, TWO-PHASE FLOW

Micro-scale, two-phase flow is found in a variety of devices such as Lab-on-a-chip, bio-chips, micro-heat exchangers, and fuel cells. Knowledge of the fluid behavior near the dynamic gas-liquid interface is required for developing accurate predictive models. Light is distorted near a curved gas-liquid interface preventing accurate measurement of interfacial shape and internal liquid velocities. This research focused on the development of experimental methods designed to isolate and probe dynamic liquid films and measure velocity fields near a moving gas-liquid interface. A high-speed, reflectance, swept-field confocal (RSFC) imaging system was developed for imaging near curved surfaces. Experimental studies of dynamic gas-liquid interface of micro-scale, two-phase flow were conducted in three phases. Dynamic liquid film thicknesses of segmented, two-phase flow were measured using the RSFC and compared to a classic film thickness deposition model. Flow fields near a steadily moving meniscus were measured using RSFC and particle tracking velocimetry. The RSFC provided high speed imaging near the menisci without distortion caused the gas-liquid interface. Finally, interfacial morphology for internal two-phase flow and droplet evaporation were measured using interferograms produced by the RSFC imaging technique. Each technique can be used independently or simultaneously when.