Influence of Mechanical and Thermal Boundary Conditions on Stabilizing/Destabilizing Mechanisms in Evaporating Liquid Films

Liquid films, evaporating or non-evaporating, are ubiquitous in nature and technology. The dynamics of evaporating liquid films is a study applicable in several industries such as water recovery, heat exchangers, crystal growth, drug design etc. The theory describing the dynamics of liquid films crosses several fields such as engineering, mathematics, material science, biophysics and volcanology to name a few. Interfacial instabilities typically manifest by the undulation of an interface from a presumed flat state or by the onset of a secondary flow state from a primary quiescent state or both. To study the instabilities affecting liquid films, an evaporating/non-evaporating Newtonian liquid film is subject to a perturbation. Numerical analysis is conducted on configurations of such liquid films being heated on solid surfaces in order to examine the various stabilizing and destabilizing mechanisms that can cause the formation of different convective structures. These convective structures have implications towards heat transfer that occurs via this process. Certain aspects of this research topic have not received attention, as will be obvious from the literature review. Static, horizontal liquid films on solid surfaces are examined for their resistance to long wave type instabilities via linear stability analysis, method of normal modes and finite difference methods. The spatiotemporal evolution equation, available in literature, describing the time evolution of a liquid film heated on a solid surface, is utilized to analyze various stabilizing/destabilizing mechanisms affecting evaporating and non-evaporating liquid films. The impact of these mechanisms on the film stability and structure for both buoyant and non-buoyant films will be examined by the variation of mechanical and thermal boundary conditions. Films evaporating in zero gravity are studied using the evolution equation. It is found that films that are stable to long wave type instabilities in terrestrial gravity are prone to destabilization via long wave instabilities in zero gravity.

USING MOSSES TO ASSESS CONDITIONS IN NORTHERN WHITE CEDAR SWAMPS

Northern white cedar (Thuja occidentalis L.) (NWC) swamps are valuable both commercially and ecologically. Unfortunately, many NWC swamps are degraded and information about them is not abundant. Especially there have been no definitive studies about mosses in northern white cedar swamps and how they react to disturbances. Mosses are sensitive to changes in their environment and thus they could be used to assess ecosystem conditions of NWC swamps. The objective of this study was to determine if mosses could be used to asses conditions in NWC swamps and if there are differences between moss communities in disturbed and undisturbed sites. Seventeen sample plots were taken from 12 disturbed and undisturbed sites around upper Michigan and northern Minnesota in the summer of 2012. All mosses occurring on the plots were identified and several associated environmental parameters were measured. The main environmental conditions affecting moss communities were identified with non-metric multidimensional scaling (NMS). Multiple response permutation procedures (MRPP) were run to ascertain if there were significant differences in community composition between disturbances. Indicator species analysis was then done to identify species that are related to different types of disturbances. A one-way ANOVA was used to check for significant differences between species richness and moss cover of undisturbed and disturbed sites. Over all sixty-two moss species were identified. The results indicate that there was no significant difference in species richness or moss cover between disturbed and undisturbed sites. However, moss community composition was affected by disturbance and strongly divided by a wetness gradient. Dicranum fuscescens was found to indicate undisturbed conditions. Calliergon cordifolium and Climacium dendroides indicated disturbed sites with wet conditions. Brotherella recurvans and Eurhynchium pulchellum indicated swamps with other disturbances.

Thickness measurement of dynamic thin liquid films generated by plug-annular flow in non-wetting microchannels

Surface tension forces are significant at millimeter length-scales, causing profoundly different flow morphologies in microchannels than in macroscale flows. The existence and morphology of thin liquid films is particularly relevant for predicting performance and operational stability of devices containing microscale two phase flows. Analytical, computational, and experimental methods previously employed in the study of thin liquid films are discussed. Thicknesses before and after a novel film morphology, referred to as a `shock,' are measured with a novel film thickness measurement technique that uses confocal microscopy. Film thicknesses predicted by previous work are compared to experimental results. Methods for increasing the accuracy of the confocal film thickness measurement technique are discussed.

DEVELOPMENT OF AN IONIC LIQUID FERROFLUID ELECTROSPRAY SOURCE AND MODE SHAPE STUDIES OF A FERROFLUID IN A NON-UNIFORM MAGNETIC FIELD

An electrospray source has been developed using a novel new fluid that is both magnetic and conductive. Unlike conventional electrospray sources that required microfabricated structures to support the fluid to be electrosprayed, this new electrospray fluid utilizes the Rosensweig instability to create the structures in the magnetic fluid when an external magnetic field was applied. Application of an external electric field caused these magnetic fluid structures to spray. These fluid based structures were found to spray at a lower onset voltage than was predicted for electrospray sources with solid structures of similar geometry. These fluid based structures were also found to be resilient to damage, unlike the solid structures found in traditional electrospray sources. Further, experimental studies of magnetic fluids in non-uniform magnetic fields were conducted. The modes of Rosensweig instabilities have been studied in-depth when created by uniform magnetic fields, but little to no studies have been performed on Rosensweig instabilities formed due to non-uniform magnetic fields. The measured spacing of the cone-like structures of ferrofluid, in a non-uniform magnetic field, were found to agree with a proposed theoretical model.

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

Non-intrusive pressure measurement in microchannels

A non-intrusive interferometric measurement technique has been successfully developed to measure fluid compressibility in both gas and liquid phases via refractive index (RI) changes. The technique, consisting of an unfocused laser beam impinging a glass channel, can be used to separate and quantify cell deflection, fluid flow rates, and pressure variations in microchannels. Currently in fields such as microfluidics, pressure and flow rate measurement devices are orders of magnitude larger than the channel cross-sections making direct pressure and fluid flow rate measurements impossible. Due to the non-intrusive nature of this technique, such measurements are now possible, opening the door for a myriad of new scientific research and experimentation. This technique, adapted from the concept of Micro Interferometric Backscatter Detection (MIBD), boasts the ability to provide comparable sensitivities in a variety of channel types and provides quantification capability not previously demonstrated in backscatter detection techniques. Measurement sensitivity depends heavily on experimental parameters such as beam impingement angle, fluid volume, photodetector sensitivity, and a channel’s dimensional tolerances. The current apparatus readily quantifies fluid RI changes of 10-5 refractive index units (RIU) corresponding to pressures of approximately 14 psi and 1 psi in water and air, respectively. MIBD reports detection capability as low as 10-9 RIU and the newly adapted technique has the potential to meet and exceed this limit providing quantification in the place of detection. Specific device sensitivities are discussed and suggestions are provided on how the technique may be refined to provide optimal quantification capabilities based on experimental conditions.

Water transport in complex, non-wetting porous layers with applications to water management in low temperature fuel cells

An experimental setup was designed to visualize water percolation inside the porous transport layer, PTL, of proton exchange membrane, PEM, fuel cells and identify the relevant characterization parameters. In parallel with the observation of the water movement, the injection pressure (pressure required to transport water through the PTL) was measured. A new scaling for the drainage in porous media has been proposed based on the ratio between the input and the dissipated energies during percolation. A proportional dependency was obtained between the energy ratio and a non-dimensional time and this relationship is not dependent on the flow regime; stable displacement or capillary fingering. Experimental results show that for different PTL samples (from different manufacturers) the proportionality is different. The identification of this proportionality allows a unique characterization of PTLs with respect to water transport. This scaling has relevance in porous media flows ranging far beyond fuel cells. In parallel with the experimental analysis, a two-dimensional numerical model was developed in order to simulate the phenomena observed in the experiments. The stochastic nature of the pore size distribution, the role of the PTL wettability and morphology properties on the water transport were analyzed. The effect of a second porous layer placed between the porous transport layer and the catalyst layer called microporous layer, MPL, was also studied. It was found that the presence of the MPL significantly reduced the water content on the PTL by enhancing fingering formation. Moreover, the presence of small defects (cracks) within the MPL was shown to enhance water management. Finally, a corroboration of the numerical simulation was carried out. A threedimensional version of the network model was developed mimicking the experimental conditions. The morphology and wettability of the PTL are tuned to the experiment data by using the new energy scaling of drainage in porous media. Once the fit between numerical and experimental data is obtained, the computational PTL structure can be used in different types of simulations where the conditions are representative of the fuel cell operating conditions.

Analysis of a Non-Traditional Micro Tube Hydroforming Process

As the demand for miniature products and components continues to increase, the need for manufacturing processes to provide these products and components has also increased. To meet this need, successful macroscale processes are being scaled down and applied at the microscale. Unfortunately, many challenges have been experienced when directly scaling down macro processes. Initially, frictional effects were believed to be the largest challenge encountered. However, in recent studies it has been found that the greatest challenge encountered has been with size effects. Size effect is a broad term that largely refers to the thickness of the material being formed and how this thickness directly affects the product dimensions and manufacturability. At the microscale, the thickness becomes critical due to the reduced number of grains. When surface contact between the forming tools and the material blanks occur at the macroscale, there is enough material (hundreds of layers of material grains) across the blank thickness to compensate for material flow and the effect of grain orientation. At the microscale, there may be under 10 grains across the blank thickness. With a decreased amount of grains across the thickness, the influence of the grain size, shape and orientation is significant. Any material defects (either natural occurring or ones that occur as a result of the material preparation) have a significant role in altering the forming potential. To date, various micro metal forming and micro materials testing equipment setups have been constructed at the Michigan Tech lab. Initially, the research focus was to create a micro deep drawing setup to potentially build micro sensor encapsulation housings. The research focus shifted to micro metal materials testing equipment setups. These include the construction and testing of the following setups: a micro mechanical bulge test, a micro sheet tension test (testing micro tensile bars), a micro strain analysis (with the use of optical lithography and chemical etching) and a micro sheet hydroforming bulge test. Recently, the focus has shifted to study a micro tube hydroforming process. The intent is to target fuel cells, medical, and sensor encapsulation applications. While the tube hydroforming process is widely understood at the macroscale, the microscale process also offers some significant challenges in terms of size effects. Current work is being conducted in applying direct current to enhance micro tube hydroforming formability. Initially, adding direct current to various metal forming operations has shown some phenomenal results. The focus of current research is to determine the validity of this process.

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.