CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

C 3-symmetric discotic liquid crystalline materials for molecular electronics: versatile synthesis and self-organization

This thesis presents the versatile synthesis and self-organization of C3-symmetric discotic nanographene molecules as well as their potential applications as materials in molecular electronics. The details can be described as follows: 1) A novel synthetic strategy towards properly designed C3 symmetric 1,3,5-tris-2’arylbenzene precursors has been developed. After the final planarization by treatment with FeCl3 under mild conditions, for the first time, it became possible to access a variety of new C3-symmetric hexa-peri-hexabenzocoronenes (HBCs) and a series of triangle-shaped nanographenes. D3 symmetric HBC with three alkyl substituents and C2 symmetric HBC with two alkyl substituents were synthesized and found to show the surprising decrease of isotropic points., the self-assembly at the liquid-solid interface displayed a unique zigzag and flower patterns. 2) Triangle-shaped discotics revealed a unique self-assembly behavior in solution, solid state as well as at the solution-substrate interface. A mesophase stability over the broad temperature range with helical supramoelcular arrangement were observed in the bulk state. The honeycomb pattern as the result of novel self-assembly was presented. Triangle-shaped discotics with swallow alkyl tails were fabricated into photovoltaic devices, the supramolecular arrangement upon thermal treatment was found to play a key role in the improvement of solar efficiency. 3) A novel class of C3 symmetric HBCs with alternating polar/apolar substituents was synthesized. Their peculiar self-assembly in solution, in the bulk and on the surface were investigated by NMR techniques, X-ray diffraction as well as different electron microscope techniques. 4) A novel concept for manipulating the intracolumnar stacking of discotics and thus for controlling the helical pitch was presented. A unique staggered stacking in the column was achieved for the first time. Theoretical simulations confirmed this self-organization and predicted that this packing should show the highest charge carrier mobility for all discotics.

Molecular organization of n-cyanobiphenyl liquid crystals on a molybdenite surface

The alignement and anchoring of liquid crystals on solid surfaces is a key problem for modern device technology that until now has been treated empirically, but that can now be tackled by atomistic computer simulations. Molecular dynamics (MD) simulations were used in this thesis work to study two films of 7 and 8 n-alkyl-4’cyanobiphenyl (7CB and 8CB) liquid crystals , with a thickness of 15 nm, confined between two (001) surfaces of MoS2 (molybdenite). The isotropic and nematic phases of both liquid crystals were simulated, and the resulting structures characterized structurally. A new force field was designed to model the interactions between the liquid crystal (LC) molecules and the surface of molybdenite, while an accurate force field developed previously was used to model the 7CB and 8CB molecules. The results show that the (001) molybdenite surface induces a planar orientation in both the liquid crystals. For the nematic phase of 8CB, one of the two solid/LC interfaces is composed of a first layer of molecules aligned parallel to the surface, followed by a second layer of molecules aligned perpendicular to the surface (also called, homeotropic). The effect of the surface appears to be local in nature as it is confined to the first 15 Angström of the LC film. Conversely, for the nematic phase of 7CB, a planar ordering is established into the LC film. The LC molecules at the interface with the molybdenite appear to align preferentially their alkyl chains toward the solid substrate. The resulting tilt angle of molecules was found to be in good agreement with experimental measurements available in literature. Despite the fact that the MD simulations spanned a time range of more than 100 ns, the nematic phases of both 7CB and 8CB were found not to be completely formed. In order to confirm the findings presented in this thesis, we propose to extend the current study.

Water-soluble organic compounds in air particulate matter analyzed by HPLC-DAD-MS and HPLC-MS/MS: abundance, sources and transformation of carboxylic acids, nitrophenols and nitrated proteins

Wasserlösliche organische Verbindungen (WSOCs) sind Hauptbestandteile atmosphärischer Aerosole, die bis zu ~ 50% und mehr der organischen Aerosolfraktion ausmachen. Sie können die optischen Eigenschaften sowie die Hygroskopizität von Aerosolpartikeln und damit deren Auswirkungen auf das Klima beeinflussen. Darüber hinaus können sie zur Toxizität und Allergenität atmosphärischer Aerosole beitragen.In dieser Studie wurde Hochleistungsflüssigchromatographie gekoppelt mit optischen Diodenarraydetektion und Massenspektrometrie (HPLC-DAD-MS und HPLC-MS/MS) angewandt, um WSOCs zu analysieren, die für verschiedene Aerosolquellen und -prozesse charakteristisch sind. Niedermolekulare Carbonsäuren und Nitrophenole wurden als Indikatoren für die Verbrennung fossiler Brennstoffe und die Entstehung sowie Alterung sekundärer organischer Aerosole (SOA) aus biogenen Vorläufern untersucht. Protein-Makromoleküle wurden mit Blick auf den Einfluss von Luftverschmutzung und Nitrierungsreaktionen auf die Allergenität primärer biologischer Aerosolpartikel – wie Pollen und Pilzsporen – untersucht.rnFilterproben von Grob- und Feinstaubwurden über ein Jahr hinweg gesammelt und auf folgende WSOCs untersucht: die Pinen-Oxidationsprodukte Pinsäure, Pinonsäure und 3-Methyl-1,2,3-Butantricarbonsäure (3-MBTCA) sowie eine Vielzahl anderer Dicarbonsäuren und Nitrophenole. Saisonale Schwankungen und andere charakteristische Merkmale werden mit Blick auf Aerosolquellen und -senken im Vergleich zu Daten anderen Studien und Regionen diskutiert. Die Verhätlnisse von Adipinsäure und Phthalsäure zu Azelainsäure deuten darauf hin, dass die untersuchten Aerosolproben hauptsächlich durch biogene Quellen beeinflusst werden. Eine ausgeprägte Arrhenius-artige Korrelation wurde zwischen der 3-MBTCA-¬Konzentration und der inversen Temperatur beobachtet (R2 = 0.79, Ea = 126±10 kJ mol-1, Temperaturbereich 275–300 K). Modellrechnungen zeigen, dass die Temperaturabhängigkeit auf eine Steigerung der photochemischen Produktionsraten von 3-MBTCA durch erhöhte OH-Radikal-Konzentrationen bei erhöhten Temperaturen zurückgeführt werden kann. Im Vergleich zur chemischen Reaktionskinetik scheint der Einfluss von Gas-Partikel-Partitionierungseffekten nur eine untergeordnete Rolle zu spielen. Die Ergebnisse zeigen, dass die OH-initiierte Oxidation von Pinosäure der geschwindigkeitsbestimmende Schritt der Bildung von 3-MBTCA ist. 3-MBTCA erscheint somit als Indikator für die chemische Alterung von biogener sekundärer organischer Aerosole (SOA) durch OH-Radikale geeignet. Eine Arrhenius-artige Temperaturabhängigkeit wurde auch für Pinäure beobachtet und kann durch die Temperaturabhängigkeit der biogenen Pinen-Emissionen als geschwindigkeitsbestimmender Schritt der Pinsäure-Bildung erklärt werden (R2 = 0.60, Ea = 84±9 kJ mol-1).rn rnFür die Untersuchung von Proteinnitrierungreaktionen wurde nitrierte Protein¬standards durch Flüssigphasenreaktion von Rinderserumalbumin (BSA) und Ovalbumin (OVA) mit Tetranitromethan (TNM) synthetisiert.Proteinnitrierung erfolgt vorrangig an den Resten der aromatischen Aminosäure Tyrosin auf, und mittels UV-Vis-Photometrie wurde der Proteinnnitrierungsgrad (ND) bestimmt. Dieser ist definiert als Verhältnis der mittleren Anzahl von Nitrotyrosinresten zur Tyrosinrest-Gesamtzahl in den Proteinmolekülen. BSA und OVA zeigten verschiedene Relationen zwischen ND und TNM/Tyrosin-Verhältnis im Reaktionsgemisch, was vermutlich auf Unterschiede in den Löslichkeiten und den molekularen Strukturen der beiden Proteine zurück zu führen ist.rnDie Nitrierung von BSA und OVA durch Exposition mit einem Gasgemisch aus Stickstoffdioxid (NO2) und Ozon (O3) wurde mit einer neu entwickelten HPLC-DAD-¬Analysemethode untersucht. Diese einfache und robuste Methode erlaubt die Bestimmung des ND ohne Hydrolyse oder Verdau der untersuchten Proteine und ernöglicht somit eine effiziente Untersuchung der Kinetik von Protein¬nitrierungs-Reaktionen. Für eine detaillierte Produktstudien wurden die nitrierten Proteine enzymatisch verdaut, und die erhaltenen Oligopeptide wurden mittels HPLC-MS/MS und Datenbankabgleich mit hoher Sequenzübereinstimmung analysiert. Die Nitrierungsgrade individueller Nitrotyrosin-Reste (NDY) korrelierten gut mit dem Gesamt-Proteinnitrierungsgrad (ND), und unterschiedliche Verhältnisse von NDY zu ND geben Aufschluss über die Regioselektivität der Reaktion. Die Nitrierungmuster von BSA und OVA nach Beahndlung mit TNM deuten darauf hin, dass die Nachbarschaft eines negativ geladenen Aminosäurerestes die Tyrosinnitrierung fördert. Die Behandlung von BSA durch NO2 und O3 führte zu anderend Nitrierungemustern als die Behandlung mit TNM, was darauf hindeutet, dass die Regioselektivität der Nitrierung vom Nitrierungsmittel abhängt. Es zeigt sich jedoch, dass Tyrosinreste in Loop-Strukturen bevorzugt und unabhängig vom Reagens nitriert werden.Die Methoden und Ergebnisse dieser Studie bilden eine Grundlage für weitere, detaillierte Untersuchungen der Reaktionskinetik sowie der Produkte und Mechanismen von Proteinnitrierungreaktionen. Sie sollen helfen, die Zusammenhänge zwischen verkehrsbedingten Luftschadstoffen wie Stickoxiden und Ozon und der Allergenität von Luftstaub aufzuklären.rn

Atomistic simulations of liquid crystals in the bulk and at their interfaces

In this thesis, we have dealt with several problems concerning liquid crystals (LC) phases, either in the bulk or at their interfaces, by the use of atomistic molecular dynamics (MD) simulations. We first focused our attention on simulating and characterizing the bulk smectic phase of 4-n-octyl-4'-cyanobiphenyl (8CB), allowing us to investigate the antiparallel molecular arrangement typical of SmAd smectic phases. A second topic of study was the characterization of the 8CB interface with vacuum by simulating freely suspended thin films, which allowed us to determine the influence of the interface on the orientational and positional order. Then we investigated the LC-water and LC-electrolyte water solution interface. This interface has recently found application in the development of sensors for several compounds, including biological molecules, and here we tried to understand the re-orientation mechanism of LC molecules at the interface which is behind the functioning of these sensors. The characterization of this peculiar interface has incidentally led us to develop a polarizable force field for the pentyl-cyanobiphenyl mesogen, whose process of parametrization and validation is reported here in detail. We have shown that this force field is a significant improvement over its previous, static charge non polarizable version in terms of density, orientational order parameter and translational diffusion.

Total OH reactivity measurements at the biosphere-atmosphere interface

The biosphere emits copiously volatile organic compounds (VOCs) into the atmosphere, which are removed again depending on the oxidative capacity of the atmosphere and physical processes such as mixing, transport and deposition. Biogenic VOCs react with the primary oxidant of the atmosphere, the hydroxyl radical (OH), and potentially lead to the formation tropospheric ozone and aerosol, which impact regional climate and air quality. The rate of OH decay in the atmosphere, the total OH reactivity is a function of the atmospheric, reactive compound's concentration and reaction velocity with OH. One way to measure the total OH reactivity, the total OH sink, is with the Comparative Reactivity Method - CRM. Basically, the reaction of OH with a reagent (here pyrrole) in clean air and in the presence of atmospheric, reactive molecules is compared. This thesis presents measurements of the total OH reactivity at the biosphere-atmosphere interface to analyze various influences and driving forces. For measurements in natural environment the instrument was automated and a direct, undisturbed sampling method developed. Additionally, an alternative detection system was tested and compared to the originally used detector (Proton Transfer Reaction-Mass Spectrometer, PTR-MS). The GC-PID (Gas Chromatographic Photo-Ionization Detector) was found as a smaller, less expensive, and robust alternative for total OH reactivity measurements. The HUMPPA-COPEC 2010 measurement campaign in the Finish forest was impacted by normal boreal forest emissions as well as prolonged heat and biomass burning emissions. The measurement of total OH reactivity was compared with a comprehensive set of monitored individual species ambient concentration levels. A significant discrepancy between those individually measured OH sinks and the total OH reactivity was observed, which was characterized in detail by the comparison of within and above the forest canopy detected OH reactivity. Direct impact of biogenic emissions on total OH reactivity was examined on Kleiner Feldberg, Germany, 2011. Trans-seasonal measurements of an enclosed Norway spruce branch were conducted via PTR-MS, for individual compound's emission rates, and CRM, for total OH reactivity emission fluxes. Especially during summertime, the individually monitored OH sink terms could not account for the measured total OH reactivity. A controlled oxidation experiment in a low NOx environment was conducted in the EUPHORE reaction chamber (CHEERS, Spain 2011). The concentration levels of the reactant isoprene and its major products were monitored and compared to total OH reactivity measurements as well as to the results of two models. The individually measured compounds could account for the total OH reactivity during this experiment as well as the traditional model-degradation scheme for isoprene (MCM 3.2). Due to previous observations of high OH levels in the isoprene-rich environment of the tropics, a novel isoprene mechanism was recently suggested. In this mechanism (MIME v4) additional OH is generated during isoprene oxidation, which could not be verified in the conditions of the CHEERS experiment.

Structure and dynamics of ionic liquid,electrode interfaces

In the early 20th century, Gouy, Chapman, and Stern developed a theory to describe the capacitance and the spatial ion distribution of diluted electrolytes near an electrode. After a century of research, considerable progress has been made in the understanding of the electrolyte/electrode interface. However, its molecular-scale structure and its variation with an applied potential is still under debate. In particular for room-temperature ionic liquids, a new class of solventless electrolytes, the classical theories for the electrical double layer are not applicable. Recently, molecular dynamics simulations and phenomenological theories have attempted to explain the capacitance of the ionic liquid/electrode interface with the molecular-scale structure and dynamics of the ionic liquid near the electrode. rnHowever, experimental evidence is very limited. rnrnIn the presented study, the ion distribution of an ionic liquid near an electrode and its response to applied potentials was examined with sub-molecular resolution. For this purpose, a new sample chamber was constructed, allowing in situ high energy X-ray reflectivity experiments under potential control, as well as impedance spectroscopy measurements. The combination of structural information and electrochmical data provided a comprehensive picture of the electric double layer in ionic liquids. Oscillatory charge density profiles were found, consisting of alternating anion- and cation-enriched layers at both, cathodic and anodic, potentials. This structure was shown to arise from the same ion-ion correlations dominating the liquid bulk structure that were observed as a distinct X-ray diffraction peak. Therefore, existing physically motivated models were refined and verified by comparison with independent measurements. rnrnThe relaxation dynamics of the interfacial structure upon potential variation were studied by time resolved X-ray reflectivity experiments with sub-millisecond resolution. The observed relaxation times during charging/discharging are consistent with the impedance spectroscopy data revealing three processes of vastly different characteristic time-scales. Initially, the ion transport normal to the interface happens on a millisecond-scale. Another 100-millisecond-scale process is associated with molecular reorientation of electrode-adsorbed cations. Further, a minute-scale relaxation was observed, which is tentatively assigned to lateral ordering within the first layer.

Spatial and functional relationships between air conduits and blood capillaries in the pulmonary gas exchange tissue of adult and developing chickens

The documented data regarding the three-dimensional structure of the air capillaries (ACs), the ultimate sites of gas exchange in the avian lung is contradictory. Further, the mode of gas exchange, described as cross-current has not been clearly elucidated. We studied the temporal and spatial arrangement of the terminal air conduits of the chicken lung and their relationship with the blood capillaries (BCs) in embryos as well as the definitive architecture in adults. Several visualization techniques that included corrosion casting, light microscopy as well as scanning and transmission electron microscopy were used. Two to six infundibulae extend from each atrium and give rise to numerous ACs that spread centrifugally. Majority of the ACs are tubular structures that give off branches, which anastomose with their neighboring cognates. Some ACs have globular shapes and a few are blind-ending tapering tubes. During inauguration, the luminal aspects of the ACs are characterized by numerous microvillus-like microplicae, which are formed during the complex processes of cell attenuation and canalization of the ACs. The parabronchial exchange BCs, initially inaugurated as disorganized meshworks, are reoriented via pillar formation to lie predominantly orthogonal to the long axes of the ACs. The remodeling of the retiform meshworks by intussusceptive angiogenesis essentially accomplishes a cross-current system at the gas exchange interface in the adults, where BCs form ring-like patterns around the ACs, thus establishing a cross-current system. Our findings clarify the mode of gas exchange in the parabronchial mantle and illuminate the basis for the functional efficiency of the avian lung.

DESIGN OF A LIQUID FUEL INJECTOR FOR ALTERNATIVE FUEL STUDIES IN AN ATMOSPHERIC MODEL GAS TURBINE COMBUSTOR

A new liquid-fuel injector was designed for use in the atmospheric-pressure, model gas turbine combustor in Bucknell University’s Combustion Research Laboratory during alternative fuel testing. The current liquid-fuel injector requires a higher-than-desired pressure drop and volumetric flow rate to provide proper atomization of liquid fuels. An air-blast atomizer type of fuel injector was chosen and an experiment utilizing water as the working fluid was performed on a variable-geometry prototype. Visualization of the spray pattern was achieved through photography and the pressure drop was measured as a function of the required operating parameters. Experimental correlations were used to estimate droplet sizes over flow conditions similar to that which would be experienced in the actual combustor. The results of this experiment were used to select the desired geometric parameters for the proposed final injector design and a CAD model was generated. Eventually, the new injector will be fabricated and tested to provide final validation of the design prior to use in the combustion test apparatus.

Efficient Liquid Liquid Extraction By Emulsion Formation and Separation in Robust Milli-Fluidic Devices

Conventional liquid liquid extraction (LLE) methods require large volumes of fluids to achieve the desired mass transfer of a solute, which is unsuitable for systems dealing with a low volume or high value product. An alternative to these methods is to scale down the process. Millifluidic devices share many of the benefits of microfluidic systems, including low fluid volumes, increased interfacial area-to-volume ratio, and predictability. A robust millifluidic device was created from acrylic, glass, and aluminum. The channel is lined with a hydrogel cured in the bottom half of the device channel. This hydrogel stabilizes co-current laminar flow of immiscible organic and aqueous phases. Mass transfer of the solute occurs across the interface of these contacting phases. Using a y-junction, an aqueous emulsion is created in an organic phase. The emulsion travels through a length of tubing and then enters the co-current laminar flow device, where the emulsion is broken and each phase can be collected separately. The inclusion of this emulsion formation and separation increases the contact area between the organic and aqueous phases, therefore increasing the area over which mass transfer can occur. Using this design, 95% extraction efficiency was obtained, where 100% is represented by equilibrium. By continuing to explore this LLE process, the process can be optimized and with better understanding may be more accurately modeled. This system has the potential to scale up to the industrial level and provide the efficient extraction required with low fluid volumes and a well-behaved system.

Implementation of a phenomenological evaporation model into a porous network simulation for water management in low temperature fuel cells

A phenomenological transition film evaporation model was introduced to a pore network model with the consideration of pore radius, contact angle, non-isothermal interface temperature, microscale fluid flows and heat and mass transfers. This was achieved by modeling the transition film region of the menisci in each pore throughout the porous transport layer of a half-cell polymer electrolyte membrane (PEM) fuel cell. The model presented in this research is compared with the standard diffusive fuel cell modeling approach to evaporation and shown to surpass the conventional modeling approach in terms of predicting the evaporation rates in porous media. The current diffusive evaporation models used in many fuel cell transport models assumes a constant evaporation rate across the entire liquid-air interface. The transition film model was implemented into the pore network model to address this issue and create a pore size dependency on the evaporation rates. This is accomplished by evaluating the transition film evaporation rates determined by the kinetic model for every pore containing liquid water in the porous transport layer (PTL). The comparison of a transition film and diffusive evaporation model shows an increase in predicted evaporation rates for smaller pore sizes with the transition film model. This is an important parameter when considering the micro-scaled pore sizes seen in the PTL and becomes even more substantial when considering transport in fuel cells containing an MPL, or a large variance in pore size. Experimentation was performed to validate the transition film model by monitoring evaporation rates from a non-zero contact angle water droplet on a heated substrate. The substrate was a glass plate with a hydrophobic coating to reduce wettability. The tests were performed at a constant substrate temperature and relative humidity. The transition film model was able to accurately predict the drop volume as time elapsed. By implementing the transition film model to a pore network model the evaporation rates present in the PTL can be more accurately modeled. This improves the ability of a pore network model to predict the distribution of liquid water and ultimately the level of flooding exhibited in a PTL for various operating conditions.

Integrated solar energy and absorption cooling model for HVAC (heating, ventilating, and air conditioning) applications in buildings

The demands in production and associate costs at power generation through non renewable resources are increasing at an alarming rate. Solar energy is one of the renewable resource that has the potential to minimize this increase. Utilization of solar energy have been concentrated mainly on heating application. The use of solar energy in cooling systems in building would benefit greatly achieving the goal of non-renewable energy minimization. The approaches of solar energy heating system research done by initiation such as University of Wisconsin at Madison and building heat flow model research conducted by Oklahoma State University can be used to develop and optimize solar cooling building system. The research uses two approaches to develop a Graphical User Interface (GUI) software for an integrated solar absorption cooling building model, which is capable of simulating and optimizing the absorption cooling system using solar energy as the main energy source to drive the cycle. The software was then put through a number of litmus test to verify its integrity. The litmus test was conducted on various building cooling system data sets of similar applications around the world. The output obtained from the software developed were identical with established experimental results from the data sets used. Software developed by other research are catered for advanced users. The software developed by this research is not only reliable in its code integrity but also through its integrated approach which is catered for new entry users. Hence, this dissertation aims to correctly model a complete building with the absorption cooling system in appropriate climate as a cost effective alternative to conventional vapor compression system.

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.

INVESTIGATING THE ROLE OF THE CONTACT LINE IN HETEROGENEOUS NUCLEATION WITH HIGH SPEED IMAGING

While nucleation of solids in supercooled liquids is ubiquitous [15, 65, 66], surface crystallization, the tendency for freezing to begin preferentially at the liquid-gas interface, has remained puzzling [74, 18, 68, 69, 51, 64, 72, 16]. Here we employ high-speed imaging of supercooled water drops to study the phenomenon of heterogeneous surface crystallization. Our geometry avoids the "point-like contact" of prior experiments by providing a simple, symmetric contact line (triple line defined by the substrate-liquid-air interface) for a drop resting on a homogeneous silicon substrate. We examine three possible mechanisms that might explain these laboratory observations: (i) Line Tension at the triple line, (ii) Thermal Gradients within the droplets and (iii) Surface Texture. In our first study we record nearly perfect spatial uniformity in the immersed (liquid-substrate) region and, thereby, no preference for nucleation at the triple line. In our second study, no influence of thermal gradients on the preference for freezing at the triple line was observed. Motivated by the conjectured importance of line tension (τ) [1, 66] for heterogeneous nucleation, we also searched for evidence of a transition to surface crystallization at length scales on the order of δ ∼ τ/σ, where σ is the surface tension [14]; poorly constrained τ [49] leads to δ ranging from microns to nanometers. We demonstrate that nano-scale texture causes a shift in the nucleation to the three-phase contact line, while micro-scale texture does not. The possibility of a critical length scale has implications for the effectiveness of nucleation catalysts, including formation of ice in atmospheric clouds [7].

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.