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.

Development of one-dimensional and two-dimensional computational tools that simulate steady internal condensing flows in terrestrial and zero-gravity environments

This dissertation presents an effective quasi one-dimensional (1-D) computational simulation tool and a full two-dimensional (2-D) computational simulation methodology for steady annular/stratified internal condensing flows of pure vapor. These simulation tools are used to investigate internal condensing flows in both gravity as well as shear driven environments. Through accurate numerical simulations of the full two dimensional governing equations, results for laminar/laminar condensing flows inside mm-scale ducts are presented. The methodology has been developed using MATLAB/COMSOL platform and is currently capable of simulating film-wise condensation for steady (and unsteady flows). Moreover, a novel 1-D solution technique, capable of simulating condensing flows inside rectangular and circular ducts with different thermal boundary conditions is also presented. The results obtained from the 2-D scientific tool and 1-D engineering tool, are validated and synthesized with experimental results for gravity dominated flows inside vertical tube and inclined channel; and, also, for shear/pressure driven flows inside horizontal channels. Furthermore, these simulation tools are employed to demonstrate key differences of physics between gravity dominated and shear/pressure driven flows. A transition map that distinguishes shear driven, gravity driven, and “mixed” driven flow zones within the non-dimensional parameter space that govern these duct flows is presented along with the film thickness and heat transfer correlations that are valid in these zones. It has also been shown that internal condensing flows in a micro-meter scale duct experiences shear driven flow, even in different gravitational environments. The full 2-D steady computational tool has been employed to investigate the length of annularity. The result for a shear driven flow in a horizontal channel shows that in absence of any noise or pressure fluctuation at the inlet, the onset of non-annularity is partly due to insufficient shear at the liquid-vapor interface. This result is being further corroborated/investigated by R. R. Naik with the help of the unsteady simulation tool. The condensing flow results and flow physics understanding developed through these simulation tools will be instrumental in reliable design of modern micro-scale and spacebased thermal systems.

AN ASSESSMENT OF FLOW REGIME MAPS AND A NUMERICAL HEAT TRANSFER CORRELATION FOR THE STRATIFIED/ANNULAR REGIME OF SHEAR-DRIVEN INTERNAL CONDENSING FLOWS

Several modern-day cooling applications require the incorporation of mini/micro-channel shear-driven flow condensers. There are several design challenges that need to be overcome in order to meet those requirements. The difficulty in developing effective design tools for shear-driven flow condensers is exacerbated due to the lack of a bridge between the physics-based modelling of condensing flows and the current, popular approach based on semi-empirical heat transfer correlations. One of the primary contributors of this disconnect is a lack of understanding caused by the fact that typical heat transfer correlations eliminate the dependence of the heat transfer coefficient on the method of cooling employed on the condenser surface when it may very well not be the case. This is in direct contrast to direct physics-based modeling approaches where the thermal boundary conditions have a direct and huge impact on the heat transfer coefficient values. Typical heat transfer correlations instead introduce vapor quality as one of the variables on which the value of the heat transfer coefficient depends. This study shows how, under certain conditions, a heat transfer correlation from direct physics-based modeling can be equivalent to typical engineering heat transfer correlations without making the same apriori assumptions. Another huge factor that raises doubts on the validity of the heat-transfer correlations is the opacity associated with the application of flow regime maps for internal condensing flows. It is well known that flow regimes influence heat transfer rates strongly. However, several heat transfer correlations ignore flow regimes entirely and present a single heat transfer correlation for all flow regimes. This is believed to be inaccurate since one would expect significant differences in the heat transfer correlations for different flow regimes. Several other studies present a heat transfer correlation for a particular flow regime - however, they ignore the method by which extents of the flow regime is established. This thesis provides a definitive answer (in the context of stratified/annular flows) to: (i) whether a heat transfer correlation can always be independent of the thermal boundary condition and represented as a function of vapor quality, and (ii) whether a heat transfer correlation can be independently obtained for a flow regime without knowing the flow regime boundary (even if the flow regime boundary is represented through a separate and independent correlation). To obtain the results required to arrive at an answer to these questions, this study uses two numerical simulation tools - the approximate but highly efficient Quasi-1D simulation tool and the exact but more expensive 2D Steady Simulation tool. Using these tools and the approximate values of flow regime transitions, a deeper understanding of the current state of knowledge in flow regime maps and heat transfer correlations in shear-driven internal condensing flows is obtained. The ideas presented here can be extended for other flow regimes of shear-driven flows as well. Analogous correlations can also be obtained for internal condensers in the gravity-driven and mixed-driven configuration.

SOLUTE-DERIVED THERMAL STABILITY OF NANOCRYSTALLINE ALUMINUM AND PROCESSING FACTOR INFLUENCE ON THE FORMATION OF AL6MN QUASICRYSTALS IN MELT-SPINNING

Thermal stability of nanograined metals can be difficult to attain due to the large driving force for grain growth that arises from the significant boundary area constituted by the nanostructure. Kinetic approaches for stabilization of the nanostructure effective at low homologous temperatures often fail at higher homologous temperatures. Thermodynamic approaches for thermal stabilization may offer higher temperature stability. In this research, modest alloying of aluminum with solute (1 at.% Sc, Yb, or Sr) was examined as a means to thermodynamically stabilize a bulk nanostructure at elevated temperatures. After using melt-spinning and ball-milling to create an extended solid-solution and nanostructure with average grain size on the order of 30-45 nm, 1 h annealing treatments at 673 K (0.72 Tm) , 773 K (0.83 Tm) , and 873 K (0.94 Tm) were applied. The alloys remain nanocrystalline (<100 nm) as measured by Warren-Averbach Fourier analysis of x-ray diffraction peaks and direct observation of TEM dark field micrographs, with the efficacy of stabilization: Sr>Yb>Sc. Disappearance of intermetallic phases in the Sr and Yb alloys in the x-ray diffraction spectra are observed to occur coincident with the stabilization after annealing, suggesting that precipitates dissolve and the boundaries are enriched with solute. Melt-spinning has also been shown to be an effective process to produce a class of ordered, but non-periodic crystals called quasicrystals. However, many of the factors related to the creation of the quasicrystals through melt-spinning are not optimized for specific chemistries and alloy systems. In a related but separate aspect of this research, meltspinning was utilized to create metastable quasicrystalline Al6Mn in an α-Al matrix through rapid solidification of Al-8Mn (by mol) and Al-10Mn (by mol) alloys. Wheel speed of the melt-spinning wheel and orifice diameter of the tube reservoir were varied to determine their effect on the resulting volume proportions of the resultant phases using integrated areas of collected x-ray diffraction spectra. The data were then used to extrapolate parameters for the Al-10Mn alloy which consistently produced Al6Mn quasicrystal with almost complete suppression of the equilibrium Al6Mn orthorhombic phase.