Investigations on the Effect of Heater Surface Characteristics on Bubble Dynamics in Subcooled Nucleate Boiling

Nucleating boiling is a repeating cycle of bubble initiation, growth and departure at many nucleation sites at the heated wall. Thereby, the bubble growth process significantly affects the dynamics of bubble departure. Experiments were performed to study the influence of heater surface characteristics, such as wettability and roughness, on single bubble growth and departure dynamics for natural circulation and upward flow boiling conditions. Self-assembled monolayer (SAM) coating, wet-etching and femtosecond pulsed laser treatment were used to alter the surface wettability and produce nano- and microstructures on stainless steel surfaces with a roughness in the range of micrometers. These surface preparation techniques allowed to separately quantify the effect of surface wettability and roughness on the bubble dynamics. The surface wettability and roughness are represented by the liquid contact angle hysteresis (θhys) and root mean square roughness of the surface (Sq). Boiling experiments were conducted at atmospheric pressure with degassed deionized water at low-subcooling. Stainless steel heater surfaces were vertically oriented during natural circulation boiling. In the experiments, bubbles were generated from an artificial nucleation cavity on the treated stainless steel heater surfaces. High-resolution optical shadowgraphy has been used to record the bubble generation, departure, sliding, detachment and inception of the next bubble. Higher bulk liquid velocity yielded smaller bubble departure diameters and slower bubble growth rates for all heater surface types. The effect of surface wettability on single bubble dynamics was studied for smooth surfaces with different liquid contact angle hysteresis. Low wetting surfaces yielded a greater bubble growth rate and departure diameter. The bubble growth rate and departure diameter were found maximum for an intermediate surface roughness Sq between 0.108 and 0.218 m. The corresponding roughness height is referred to as the ‘optimal roughness height’ in this work. Surface roughness was found very influential to the bubble growth and departure, which can be explained by considering its interaction with the microlayer underneath a bubble. The role of the heater surface parameters for the bubble growth was qualitatively assessed by evaluating the microlayer thickness constant C2. Hence, an improved bubble growth model was derived in this work. The bubble growth model was formulated on the basis of the evaporation of the microlayer beneath a bubble with the dryout area, inertia and heat diffusion controlled bubble growth and condensation at the bubble cap. The model can also predict the superheated liquid layer around a bubble which helps to determine the portion of a bubble that is in contact with the subcooled liquid. As bubble growth Abstract is highly dependent on the effective interactions of heater surface roughness and microlayer, a term Ceff was introduced in the bubble growth model. The effective microlayer thickness constant Ceff incorporates the impact of heater surface characteristics on the bubble growth process until the departure of a bubble. The bubble growth model was utilized in the analysis of high-resolution experimental data of steam bubble growth and the values of Ceff were calculated for different heater surface characteristics. The value of Ceff was found to decrease with the increase of bubble growth rate. A simplified model for the bubble departure criterion was derived from the expressions of forces which act on a nucleating bubble throughout its growth cycle. It was found that 90% of the departing bubbles satisfy the bubble departure criterion model with ±25% deviation. The knowledge gained from this work shall be particularly useful to improve nucleate boiling models for numerical simulations. The findings are also useful for designing heater surfaces in the future.:Abstract v Kurzfassung vii Acknowledgements xiii Abbreviations and Symbols xv Chapter 1: Introduction and Motivation 1 1.1 General overview 1 1.2 Theoretical background 3 1.3 Objectives and outline of the thesis 7 Chapter 2: Fundamentals of Bubble Dynamics in Nucleate Boiling 9 2.1 Bubble growth in nucleate boiling 9 2.2 Bubble growth models 12 2.3 The physical process of bubble departure 16 2.4 Experimental investigations of bubble dynamics 20 2.4.1 Effects of heater surface characteristics 21 2.4.2 Effects of bulk liquid velocity 24 2.5 Chapter conclusion 26 Chapter 3: Heater Surface Preparation and Characterization 27 3.1 Surface properties 27 3.2 Surface preparation 29 3.2.1 Self-assembled monolayer coating 30 3.2.2 High-power pulsed laser irradiation 31 3.2.3 Wet-etching 32 3.3 Surface cleaning 32 3.4 Surface characterization 32 3.4.1 Wettability measurement 32 3.4.2 Roughness measurement 33 3.4.3 Analysis of surface characteristics 34 3.4.4 Uncertainty of surface parameters 38 3.5 Artificial cavity preparation 38 Chapter 4: Experimental Setup and Procedure 41 4.1 Natural circulation boiling (NCB 41 4.1.1 Experimental procedure and measurement techniques 41 4.1.2 Uncertainty analysis 44 4.2 Upward flow boiling (UFB) 45 4.2.1 Experimental procedure and measurement techniques 45 4.2.2 Uncertainty analysis 48 4.3 Image processing 50 Chapter 5: Experimental Results 53 5.1 Introduction to the analysis of the bubble dynamics 53 5.1.1 The bubble life cycle 53 5.1.2 Calculation of the bubble equivalent diameter 55 5.1.3 Bubble dynamics with the increase of heat flux 57 5.1.4 Qualitative assessment of the bubble dynamics for different parameters 60 5.2 Bubble dynamics 61 5.2.1 Effect of heater surface wettability 61 5.2.2 Effect of heater surface roughness 65 5.2.3 Effect of bulk liquid velocity 70 5.3 Bubble departure 76 5.3.1 Effect of heater surface wettablity 76 5.3.2 Effect of heater surface roughness 76 5.3.3 Effect of bulk liquid velocity 78 5.4 Chapter conclusion 79 Chapter 6: Analysis and Model Development 81 6.1 Numerical evaluation of the role of heater surface characteristics 81 6.1.1 Derivation of an improved bubble growth model 86 6.1.2 Calculation of Ceff 82 6.2 Effect of liquid velocity on the bubble growth 93 6.3 Improved modeling of bubble departure 95 6.3.1 Analysis of important parameters 95 6.3.2 Formulation of a bubble departure criterion 100 6.4 Chapter conclusion 102 Chapter 7: Summary and Outlook 105 Bibliography 109 List of Figures 121 List of Tables 127 Appendix: Surface Parameters and Profile 129Der Blasenabriss von einer Keimstellenkavität ist ein komplexer Ablösemechanismus und spielt eine wichtige Rolle beim Wärmetransport. Zur Beschreibung der Blasendynamik sind Kenntnisse über den Blasenwachstumsprozess sowie die Vorhersage eines Kriteriums für die Blasenablösung erforderlich. In den existierenden Blasenwachstums- und Blasenablösungsmodellen wird die Oberflächencharakteristik des Heizers bisher nicht berücksichtigt. Im Rahmen dieser Promotion wurden Experimente durchgeführt, um den Einfluss der Heizeroberfläche und der Hauptströmungsgeschwindigkeit auf diese Parameter für eine vertikale Heizfläche zu untersuchen. Hierbei wurden das Naturkonvektionssieden und das aufwärtsgerichtete Strömungssieden betrachtet. Die Experimente wurden mit vollentsalztem Wasser bei einer Unterkühlung zwischen 1,68 und 4,00 K bei Atmosphärendruck und einem aus Edelstahl gefertigten Heizer durchgeführt, dessen Oberfläche anhand der Parameter Oberflächenrauigkeit und Benetzbarkeit charakterisiert ist. Unterschiedliche Oberflächenbearbeitungstechniken, wie Beschichtung durch Self-Assembled Monolayer (SAM), Nass-Ätzen und Hochleistungspuls-Laserbestrahlung wurden genutzt, um die Oberflächenbenetzung und –rauigkeit zu modifizieren. Der Unterschied zwischen dem gemessenen Fortschritts- (θadv) und Rückzugskontaktwinkel (θrec) der Flüssigkeit wird als Flüssigkeitskontaktwinkelhysterese (θhys) bezeichnet und beschreibt die Oberflächenbenetzbarkeit. Die Oberflächenrauigkeit wurde durch ein Konfokal-Mikroskop bestimmt und durch das gemittelte Quadrat der Rauigkeit (Sq) und den Maximalwert der Rauigkeit (St) definiert. Insgesamt wurden 18 unterschiedliche Heizoberflächen mit einer Größe von 130 x 20 mm² untersucht. Davon kamen jeweils die Hälfte für das Naturkonvektionssieden bzw. aufwärtsgerichtetes Strömungssieden zur Anwendung. Der Einfluss der Oberflächenbenetzbarkeit auf die Blasendynamik wurde für polierte Oberflächen (Sq 0,01 μm) analysiert. Die Wirkung der Oberflächenrauigkeit auf die Blasendynamik wurde für konstante Flüssigkeitskontaktwinkelhysteresen von 40,05°±1,5° und 59,97°±1,5° für Naturzirkulation und Strömungssieden untersucht. Eine künstliche zylindrische Kavität mit einer Fläche von 1963,5 m² und einer Tiefe von 50 m wurde mittels Mikrolaser in die Heizoberflächen eingebracht, um die Blasen in einer spezifischen Position zu erzeugen. Während des Naturkonvektionssiedens betrug die Wärmestromdichte 19,22 bis 30,29 kW/m². Bei den Experimenten mit aufwärtsgerichtetem Strömungssieden wurde die Hauptströmungsgeschwindigkeit im Bereich von 0,052 bis 0,183 m/s variiert und eine Appendix: Surface Parameters and Profile Wärmestromdichte zwischen 39,41 und 45,47 kW/m² aufgeprägt. Daraus resultierten insgesamt 87 Experimentalserien. Um den Blasenlebenszyklus zu erfassen, wurde hochauflösende Bildgebungstechnik verwendet. Mit der Bildverarbeitungssoftware ImageJ wurden die erfassten Videos weiterverarbeitet. Die Temperatur der Hauptströmung wurde mit Typ-K Thermoelementen gemessen. Die zeit- und ortsgemittelten Heizerwandtemperaturen wurden für die Naturzirkulation durch Infrarotthermografie und für das aufwärtsgerichtete Strömungssieden durch Typ-K Thermoelemente erfasst. Die mittlere Flüssigkeitsgeschwindigkeit wurde bei der Naturzirkulation mittels Particle Image Velocimetry (PIV) und beim Strömungssieden mittels Coriolis-Durchflusszähler bestimmt. Eine hochauflösende optische Schattenbildtechnik diente zur Aufzeichnung der Hauptphasen des Blasenlebenszyklus: Blasenerzeugung, Blasenwachstum, Blasenablösung, Blasengleiten und Blasenabriss. In dieser Arbeit wurden die der Blasenablösung vorrausgehenden Phasen untersucht. Blasenhöhe, Blasenbreite, Blasenbasisdurchmesser und Schwerpunkt der Blase wurden mit Hilfe der Bildverarbeitung ermittelt. Der blasenäquivalente Durchmesser wurde mittels des geometrischen Mittelwertes, der Blasenbreite und der Blasenhöhe berechnet. Basierend auf den Messdaten können folgende Erkenntnisse für das Blasenwachstum und den Blasenablösemechanismus postuliert werden: (i) Eine höhere Wärmeströmedichte führen zu größen Blasen und kürzeren Wachstumsperioden. Der Einfluss der Oberflächenbenetzbarkeit und der Oberflächenrauigkeit auf die Blasendynamik zeigt ähnliche Tendenzen für Naturkonvektion und aufwärtsgerichtetes Strömungssieden. (ii) Eine höhere Flüssigkeitskontaktwinkelhysterese führt zu einer schnelleren Expansion der Blasenbasis und zu einem schnellern Blasenwachstum. Für gut benetzbare Oberflächen bewegt sich der Blasenschwerpunkt schneller entlang der Strömungsrichtung. Für Oberflächen mit geringer Benetzbarkeit ist die Blasengröße vor der Blasenablösung größer und die Ablöseperiode länger. Der mittlere Blasenablösedurchmesser für unterschiedliche Hauptströmungsgeschwindigkeiten der Flüssigkeit erhöht sich von 0,75 auf 1,75 mm bei zunehmender Flüssigkeitskontaktwinkelhysterese von 42,32° auf 62,30°. (iii) Eine, bezogen auf die Mikrogrenzschichtdicke, optimale Oberflächenrauigkeit erhöht die Blasenwachstumsrate und die Blasengröße. Dieses Ergebnis ist bisher einzigartig bei der Untersuchung der Einzelblasendynamik beim Blasensieden. Die Expansion der Blasenbasis und der Blasenwachstumsrate erreicht ein Maximum für das gemittelte Quadrat der Rauigkeit (Sq) im Bereich zwischen 0,156 und 0,202 m für Naturzirkulation. Für aufwärtsgerichtetes Strömungssieden war die Expansion der Blasenbasis und die Blasenwachstumsrate für Sq-Werte zwischen 0,108 und 0,218 m maximal. Der Blasenablösedurchmesser wurde für einen großen Bereich der Hauptströmungsgeschwindigkeiten und Wärmestromedichte gemittelt. Das Maximum des mittleren Ablösedurchmessers wurde für die Oberfläche mit einem Wert von Sq = 0,218 m erreicht. Die Oberflächenrauigkeit erweitert die Wärmeübertragungsoberfläche neben der Blasenbasis. Der Einfluss der Oberflächenrauigkeitshöhe auf die Blasen hängt von der Mikrogrenzschichtdicke sowie vom Blasenbasisradius ab. Das Modell der Mikrogrenzschichtdicke von Cooper und Lloyd [1] und die konzeptionelle Idee zur Störung der Mikrogrenzschicht durch die Rautiefe von Sriraman [2] wurden analysiert. Es wurde nachgewiesen, dass die Oberflächenrauigkeit die effektive Mikrogrenzschichtdicke und die dazugehörige Wärmeübertragung beeinflusst. (iv) Es wurden geringere Blasenwachstumsraten für höhere Hauptströmungs-geschwindigkeiten gemessen. Weiterhin reduzieren sich der Blasenablösedurchmesser sowie Ablöseperioden mit zunehmender Hauptströmungsgeschwindigkeit bei unterschiedlichen Wärmeoberflächencharakteristiken. Bei niedrigen Hauptströmungs-geschwindigkeiten im Bereich zwischen ungefähr 0,052 und 0,16 m/s reduziert sich der durchschnittliche Blasenablösedurchmesser deutlich. Die experimentellen Ergebnisse zeigen einen wesentlichen Einfluss der Oberflächenbeschaffenheit auf das Blasenwachstum und den Ablöseprozess beim Blasensieden. Um diesen Einfluss numerisch zu charakterisieren, wurde ein neues Blasenwachstumsmodel entwickelt. Existierende Blasenwachstumsmodelle berücksichtigen den umfangreichen Einfluss der Oberfläche des Heizers bisher nicht. Das vorgeschlagene Model bezieht die plausibelsten Mechanismen des Blasensiedens mit ein. Dazu zählen: Mikrogrenzschichtverdampfung im Bereich der Austrocknung, trägheits- und wärmediffusionskontrolliertes Blasenwachstum und Kondensation an der Blasenoberseite. Das Modell berücksichtigt, dass die überhitzte Flüssigkeitsschicht an der Heizerwand durch die wachsende Blase nach außen verdrängt wird und die so gestreckte Flüssigkeitsschicht einen Teil der Blase einhüllt. Kondensation erfolgt an der Blasengrenze, die in Kontakt mit der unterkühlten Flüssigkeit steht, und demzufolge mit der überhitzen Flüssigkeitsschicht nicht in Kontakt kommt. Das vorgeschlagene Blasenwachstumsmodel arbeitet mit drei Konstanten für die beschriebenen Wärmeübertragungsmechanismen beim Blasenwachstum. Dabei handelt es sich um eine Konstante für die effektive Mikrogrenzschichtdicke (Ceff ), eine weitere Konstante ́ für die Wärmediffusion hin zur Blase und der Trägheit sowie letztendlich einer Konstante S zur Abbildung des Kondensationswärmeübergangs, anhand der Beschreibung des Anteils der Blase, welcher in Kontakt mit der unterkühlten Flüssigkeit steht. Die effektive Mikrogrenzschichtdickenkonstante (Ceff) definiert den Einfluss der Heizoberflächencharakteristik auf die Verdampfung der Mikrogrenzschicht und somit die Blasenwachstumsrate beim Blasensieden. Die numerisch berechnete und experimentell gemessene Blasengröße wurde verglichen, um die Mikrogrenzschichtdickenkonstante Ceff zu definieren. Der Einfluss der Kondensation auf Ceff wurde geprüft.:Abstract v Kurzfassung vii Acknowledgements xiii Abbreviations and Symbols xv Chapter 1: Introduction and Motivation 1 1.1 General overview 1 1.2 Theoretical background 3 1.3 Objectives and outline of the thesis 7 Chapter 2: Fundamentals of Bubble Dynamics in Nucleate Boiling 9 2.1 Bubble growth in nucleate boiling 9 2.2 Bubble growth models 12 2.3 The physical process of bubble departure 16 2.4 Experimental investigations of bubble dynamics 20 2.4.1 Effects of heater surface characteristics 21 2.4.2 Effects of bulk liquid velocity 24 2.5 Chapter conclusion 26 Chapter 3: Heater Surface Preparation and Characterization 27 3.1 Surface properties 27 3.2 Surface preparation 29 3.2.1 Self-assembled monolayer coating 30 3.2.2 High-power pulsed laser irradiation 31 3.2.3 Wet-etching 32 3.3 Surface cleaning 32 3.4 Surface characterization 32 3.4.1 Wettability measurement 32 3.4.2 Roughness measurement 33 3.4.3 Analysis of surface characteristics 34 3.4.4 Uncertainty of surface parameters 38 3.5 Artificial cavity preparation 38 Chapter 4: Experimental Setup and Procedure 41 4.1 Natural circulation boiling (NCB 41 4.1.1 Experimental procedure and measurement techniques 41 4.1.2 Uncertainty analysis 44 4.2 Upward flow boiling (UFB) 45 4.2.1 Experimental procedure and measurement techniques 45 4.2.2 Uncertainty analysis 48 4.3 Image processing 50 Chapter 5: Experimental Results 53 5.1 Introduction to the analysis of the bubble dynamics 53 5.1.1 The bubble life cycle 53 5.1.2 Calculation of the bubble equivalent diameter 55 5.1.3 Bubble dynamics with the increase of heat flux 57 5.1.4 Qualitative assessment of the bubble dynamics for different parameters 60 5.2 Bubble dynamics 61 5.2.1 Effect of heater surface wettability 61 5.2.2 Effect of heater surface roughness 65 5.2.3 Effect of bulk liquid velocity 70 5.3 Bubble departure 76 5.3.1 Effect of heater surface wettablity 76 5.3.2 Effect of heater surface roughness 76 5.3.3 Effect of bulk liquid velocity 78 5.4 Chapter conclusion 79 Chapter 6: Analysis and Model Development 81 6.1 Numerical evaluation of the role of heater surface characteristics 81 6.1.1 Derivation of an improved bubble growth model 86 6.1.2 Calculation of Ceff 82 6.2 Effect of liquid velocity on the bubble growth 93 6.3 Improved modeling of bubble departure 95 6.3.1 Analysis of important parameters 95 6.3.2 Formulation of a bubble departure criterion 100 6.4 Chapter conclusion 102 Chapter 7: Summary and Outlook 105 Bibliography 109 List of Figures 121 List of Tables 127 Appendix: Surface Parameters and Profile 12

Dynamic Scaling, Data-collapse and Self-Similarity in Mediation-Driven Attachment Networks

Recently, we have shown that if the $i$th node of the Barab\'{a}si-Albert (BA) network is characterized by the generalized degree $q_i(t)=k_i(t)t_i^\beta/m$, where $k_i(t)\sim t^\beta$ and $m$ are its degree at current time $t$ and at birth time $t_i$, then the corresponding distribution function $F(q,t)$ exhibits dynamic scaling. Applying the same idea to our recently proposed mediation-driven attachment (MDA) network, we find that it too exhibits dynamic scaling but, unlike the BA model, the exponent $\beta$ of the MDA model assumes a spectrum of value $1/2\leq \beta \leq 1$. Moreover, we find that the scaling curves for small $m$ are significantly different from those of the larger $m$ and the same is true for the BA networks albeit in a lesser extent. We use the idea of the distribution of inverse harmonic mean (IHM) of the neighbours of each node and show that the number of data points that follow the power-law degree distribution increases as the skewness of the IHM distribution decreases. Finally, we show that both MDA and BA models become almost identical for large $m$.Comment: 8 pages, 8 captioned figure

Exploring Work-Life Conflict in Global Software Development (GSD) Contexts: A Survey of IT Professionals based in India

Global Software Development (GSD) is now a mega-trend. While there is a rich literature exploring various facets of the GSD phenomenon, few (if any) studies have focused on the working conditions of IT professionals, specifically their work-life conflict. In this paper, we discuss our research-in-progress on this issue, wherein we empirically examine the effects of four categories of relevant antecedents (individual factors, organizational factors, IS project-based factors, and the nature of personnel distribution in teams) on work-life conflict, and the effect of work-life conflict on outcome variables such as organizational commitment and individual’s performance. Analysis of data collected as part of an on-going study show that the measurement instruments are valid and reliable, and many of the hypothesized relationships hold. The aspiration of this study is to be among the first to empirically examine work-life conflict (WLC) issues in a GSD setting

The Role Of Individual, Family-Related, And Organizational Factors In Shaping WLC In Offshoring Contexts: A Study Of European And Indian IT Professionals

Today, we operate in a networked world, where organizations frequently resort to offshoring, such that work gets accomplished by a globally distributed workforce, whether inside or outside the organizational boundaries. Much of the past research on offshoring has focused on economic rationale and benefits as well as risks associated with offshoring. Offshoring leverages human capital in different parts of the globe, and the issues of WLC (WLC) faced by offshoring workforce issues can have substantial impact on the effectiveness of the offshoring arrangements. In spite of this, WLC in the context of offshoring has not received the attention it deserves. In this paper, based on responses from global software development (GSD) professionals based in Europe and India, we examine the key individual-based, family-based, and organization-based factors that can potentially have an impact on IT workers’ WLC. Further, we investigate the impact on WLC on outcome variables such as the workers’ job satisfaction as well as valence regarding working in offshoring environments in the future. Some interesting differences in factors relevant to workers in Europe and India emerge from our analysis

Legumes in conservation agriculture: A sustainable approach in rice-based ecology of the Eastern Indo-Gangetic Plain of South Asia − an overview

Legumes in Conservation Agriculture (CA) have the potential to increase crop productivity and sustainability of the rice-based system. However, there is limited information available on the importance of legume crops in CA in the rice-based system of the Eastern Indo-Gangetic Plain (EIGP). Rice-based cropping is the mainstay of the people in Bangladesh, on the EIGP. These systems are the major food supplier to the region. In addition, they provide income and employment opportunities to the majority of people in the region. However, the system is facing unprecedented challenges and increased risk due to water, energy, labour and capital scarcity, which are exaggerated due to the effects of climate change. This could be further aggravated by deteriorating soil health, depleting underground water, and reduced land and water productivity which ultimately threaten sustainable food production and food security of the EIGP. Hence, sustainable crop intensification is essential, but increasing cropping intensity has reduced the yield of single crops by degrading soil properties. To address these constraints, Conservation Agriculture (CA), with a minimum tillage system, residue retention and crop diversification with legumes, could be an effective approach for improving crop productivity while sustaining the natural resources in intensive rice-based systems of EIGP. The addition of legumes in crop rotation is a suitable technique for crop diversification due to its multiple benefits related to soil health and natural resources. Diverse legume crops involving rice-based cropping exist at different agro-ecological zones in Bangladesh, and their rotation definitely could act a major role in promoting the CA in rice-based systems. Legume-based rotation offers multiple benefits, such as biological nitrogen fixation, improves soil pores through the deep root system, P-availability, soil fertility and enhanced nutrient cycling, and reducing the use of external input and thereby minimizing greenhouse gas emission and groundwater pollution, improving water productivity, and minimizes diseases and pest incidence. As a result, crop rotation with legumes has a high potential for CA and sustainable rice-based cropping systems in Bangladesh. The gaps between legume and non-legume crops in CA for each parameter suggest a noteworthy possibility for the improvement of rice-based systems in EIGP. This review suggests further sustainability improvements can be achieved through future field research focused on the inclusion of legume crops in the diverse rice-based systems under CA

Investigations on the Effect of Heater Surface Characteristics on Bubble Dynamics in Subcooled Nucleate Boiling

Nucleating boiling is a repeating cycle of bubble initiation, growth and departure at many nucleation sites at the heated wall. Thereby, the bubble growth process significantly affects the dynamics of bubble departure. Experiments were performed to study the influence of heater surface characteristics, such as wettability and roughness, on single bubble growth and departure dynamics for natural circulation and upward flow boiling conditions. Self-assembled monolayer (SAM) coating, wet-etching and femtosecond pulsed laser treatment were used to alter the surface wettability and produce nano- and microstructures on stainless steel surfaces with a roughness in the range of micrometers. These surface preparation techniques allowed to separately quantify the effect of surface wettability and roughness on the bubble dynamics. The surface wettability and roughness are represented by the liquid contact angle hysteresis (θhys) and root mean square roughness of the surface (Sq). Boiling experiments were conducted at atmospheric pressure with degassed deionized water at low-subcooling. Stainless steel heater surfaces were vertically oriented during natural circulation boiling. In the experiments, bubbles were generated from an artificial nucleation cavity on the treated stainless steel heater surfaces. High-resolution optical shadowgraphy has been used to record the bubble generation, departure, sliding, detachment and inception of the next bubble. Higher bulk liquid velocity yielded smaller bubble departure diameters and slower bubble growth rates for all heater surface types. The effect of surface wettability on single bubble dynamics was studied for smooth surfaces with different liquid contact angle hysteresis. Low wetting surfaces yielded a greater bubble growth rate and departure diameter. The bubble growth rate and departure diameter were found maximum for an intermediate surface roughness Sq between 0.108 and 0.218 m. The corresponding roughness height is referred to as the ‘optimal roughness height’ in this work. Surface roughness was found very influential to the bubble growth and departure, which can be explained by considering its interaction with the microlayer underneath a bubble. The role of the heater surface parameters for the bubble growth was qualitatively assessed by evaluating the microlayer thickness constant C2. Hence, an improved bubble growth model was derived in this work. The bubble growth model was formulated on the basis of the evaporation of the microlayer beneath a bubble with the dryout area, inertia and heat diffusion controlled bubble growth and condensation at the bubble cap. The model can also predict the superheated liquid layer around a bubble which helps to determine the portion of a bubble that is in contact with the subcooled liquid. As bubble growth Abstract is highly dependent on the effective interactions of heater surface roughness and microlayer, a term Ceff was introduced in the bubble growth model. The effective microlayer thickness constant Ceff incorporates the impact of heater surface characteristics on the bubble growth process until the departure of a bubble. The bubble growth model was utilized in the analysis of high-resolution experimental data of steam bubble growth and the values of Ceff were calculated for different heater surface characteristics. The value of Ceff was found to decrease with the increase of bubble growth rate. A simplified model for the bubble departure criterion was derived from the expressions of forces which act on a nucleating bubble throughout its growth cycle. It was found that 90% of the departing bubbles satisfy the bubble departure criterion model with ±25% deviation. The knowledge gained from this work shall be particularly useful to improve nucleate boiling models for numerical simulations. The findings are also useful for designing heater surfaces in the future.:Abstract v Kurzfassung vii Acknowledgements xiii Abbreviations and Symbols xv Chapter 1: Introduction and Motivation 1 1.1 General overview 1 1.2 Theoretical background 3 1.3 Objectives and outline of the thesis 7 Chapter 2: Fundamentals of Bubble Dynamics in Nucleate Boiling 9 2.1 Bubble growth in nucleate boiling 9 2.2 Bubble growth models 12 2.3 The physical process of bubble departure 16 2.4 Experimental investigations of bubble dynamics 20 2.4.1 Effects of heater surface characteristics 21 2.4.2 Effects of bulk liquid velocity 24 2.5 Chapter conclusion 26 Chapter 3: Heater Surface Preparation and Characterization 27 3.1 Surface properties 27 3.2 Surface preparation 29 3.2.1 Self-assembled monolayer coating 30 3.2.2 High-power pulsed laser irradiation 31 3.2.3 Wet-etching 32 3.3 Surface cleaning 32 3.4 Surface characterization 32 3.4.1 Wettability measurement 32 3.4.2 Roughness measurement 33 3.4.3 Analysis of surface characteristics 34 3.4.4 Uncertainty of surface parameters 38 3.5 Artificial cavity preparation 38 Chapter 4: Experimental Setup and Procedure 41 4.1 Natural circulation boiling (NCB 41 4.1.1 Experimental procedure and measurement techniques 41 4.1.2 Uncertainty analysis 44 4.2 Upward flow boiling (UFB) 45 4.2.1 Experimental procedure and measurement techniques 45 4.2.2 Uncertainty analysis 48 4.3 Image processing 50 Chapter 5: Experimental Results 53 5.1 Introduction to the analysis of the bubble dynamics 53 5.1.1 The bubble life cycle 53 5.1.2 Calculation of the bubble equivalent diameter 55 5.1.3 Bubble dynamics with the increase of heat flux 57 5.1.4 Qualitative assessment of the bubble dynamics for different parameters 60 5.2 Bubble dynamics 61 5.2.1 Effect of heater surface wettability 61 5.2.2 Effect of heater surface roughness 65 5.2.3 Effect of bulk liquid velocity 70 5.3 Bubble departure 76 5.3.1 Effect of heater surface wettablity 76 5.3.2 Effect of heater surface roughness 76 5.3.3 Effect of bulk liquid velocity 78 5.4 Chapter conclusion 79 Chapter 6: Analysis and Model Development 81 6.1 Numerical evaluation of the role of heater surface characteristics 81 6.1.1 Derivation of an improved bubble growth model 86 6.1.2 Calculation of Ceff 82 6.2 Effect of liquid velocity on the bubble growth 93 6.3 Improved modeling of bubble departure 95 6.3.1 Analysis of important parameters 95 6.3.2 Formulation of a bubble departure criterion 100 6.4 Chapter conclusion 102 Chapter 7: Summary and Outlook 105 Bibliography 109 List of Figures 121 List of Tables 127 Appendix: Surface Parameters and Profile 129Der Blasenabriss von einer Keimstellenkavität ist ein komplexer Ablösemechanismus und spielt eine wichtige Rolle beim Wärmetransport. Zur Beschreibung der Blasendynamik sind Kenntnisse über den Blasenwachstumsprozess sowie die Vorhersage eines Kriteriums für die Blasenablösung erforderlich. In den existierenden Blasenwachstums- und Blasenablösungsmodellen wird die Oberflächencharakteristik des Heizers bisher nicht berücksichtigt. Im Rahmen dieser Promotion wurden Experimente durchgeführt, um den Einfluss der Heizeroberfläche und der Hauptströmungsgeschwindigkeit auf diese Parameter für eine vertikale Heizfläche zu untersuchen. Hierbei wurden das Naturkonvektionssieden und das aufwärtsgerichtete Strömungssieden betrachtet. Die Experimente wurden mit vollentsalztem Wasser bei einer Unterkühlung zwischen 1,68 und 4,00 K bei Atmosphärendruck und einem aus Edelstahl gefertigten Heizer durchgeführt, dessen Oberfläche anhand der Parameter Oberflächenrauigkeit und Benetzbarkeit charakterisiert ist. Unterschiedliche Oberflächenbearbeitungstechniken, wie Beschichtung durch Self-Assembled Monolayer (SAM), Nass-Ätzen und Hochleistungspuls-Laserbestrahlung wurden genutzt, um die Oberflächenbenetzung und –rauigkeit zu modifizieren. Der Unterschied zwischen dem gemessenen Fortschritts- (θadv) und Rückzugskontaktwinkel (θrec) der Flüssigkeit wird als Flüssigkeitskontaktwinkelhysterese (θhys) bezeichnet und beschreibt die Oberflächenbenetzbarkeit. Die Oberflächenrauigkeit wurde durch ein Konfokal-Mikroskop bestimmt und durch das gemittelte Quadrat der Rauigkeit (Sq) und den Maximalwert der Rauigkeit (St) definiert. Insgesamt wurden 18 unterschiedliche Heizoberflächen mit einer Größe von 130 x 20 mm² untersucht. Davon kamen jeweils die Hälfte für das Naturkonvektionssieden bzw. aufwärtsgerichtetes Strömungssieden zur Anwendung. Der Einfluss der Oberflächenbenetzbarkeit auf die Blasendynamik wurde für polierte Oberflächen (Sq 0,01 μm) analysiert. Die Wirkung der Oberflächenrauigkeit auf die Blasendynamik wurde für konstante Flüssigkeitskontaktwinkelhysteresen von 40,05°±1,5° und 59,97°±1,5° für Naturzirkulation und Strömungssieden untersucht. Eine künstliche zylindrische Kavität mit einer Fläche von 1963,5 m² und einer Tiefe von 50 m wurde mittels Mikrolaser in die Heizoberflächen eingebracht, um die Blasen in einer spezifischen Position zu erzeugen. Während des Naturkonvektionssiedens betrug die Wärmestromdichte 19,22 bis 30,29 kW/m². Bei den Experimenten mit aufwärtsgerichtetem Strömungssieden wurde die Hauptströmungsgeschwindigkeit im Bereich von 0,052 bis 0,183 m/s variiert und eine Appendix: Surface Parameters and Profile Wärmestromdichte zwischen 39,41 und 45,47 kW/m² aufgeprägt. Daraus resultierten insgesamt 87 Experimentalserien. Um den Blasenlebenszyklus zu erfassen, wurde hochauflösende Bildgebungstechnik verwendet. Mit der Bildverarbeitungssoftware ImageJ wurden die erfassten Videos weiterverarbeitet. Die Temperatur der Hauptströmung wurde mit Typ-K Thermoelementen gemessen. Die zeit- und ortsgemittelten Heizerwandtemperaturen wurden für die Naturzirkulation durch Infrarotthermografie und für das aufwärtsgerichtete Strömungssieden durch Typ-K Thermoelemente erfasst. Die mittlere Flüssigkeitsgeschwindigkeit wurde bei der Naturzirkulation mittels Particle Image Velocimetry (PIV) und beim Strömungssieden mittels Coriolis-Durchflusszähler bestimmt. Eine hochauflösende optische Schattenbildtechnik diente zur Aufzeichnung der Hauptphasen des Blasenlebenszyklus: Blasenerzeugung, Blasenwachstum, Blasenablösung, Blasengleiten und Blasenabriss. In dieser Arbeit wurden die der Blasenablösung vorrausgehenden Phasen untersucht. Blasenhöhe, Blasenbreite, Blasenbasisdurchmesser und Schwerpunkt der Blase wurden mit Hilfe der Bildverarbeitung ermittelt. Der blasenäquivalente Durchmesser wurde mittels des geometrischen Mittelwertes, der Blasenbreite und der Blasenhöhe berechnet. Basierend auf den Messdaten können folgende Erkenntnisse für das Blasenwachstum und den Blasenablösemechanismus postuliert werden: (i) Eine höhere Wärmeströmedichte führen zu größen Blasen und kürzeren Wachstumsperioden. Der Einfluss der Oberflächenbenetzbarkeit und der Oberflächenrauigkeit auf die Blasendynamik zeigt ähnliche Tendenzen für Naturkonvektion und aufwärtsgerichtetes Strömungssieden. (ii) Eine höhere Flüssigkeitskontaktwinkelhysterese führt zu einer schnelleren Expansion der Blasenbasis und zu einem schnellern Blasenwachstum. Für gut benetzbare Oberflächen bewegt sich der Blasenschwerpunkt schneller entlang der Strömungsrichtung. Für Oberflächen mit geringer Benetzbarkeit ist die Blasengröße vor der Blasenablösung größer und die Ablöseperiode länger. Der mittlere Blasenablösedurchmesser für unterschiedliche Hauptströmungsgeschwindigkeiten der Flüssigkeit erhöht sich von 0,75 auf 1,75 mm bei zunehmender Flüssigkeitskontaktwinkelhysterese von 42,32° auf 62,30°. (iii) Eine, bezogen auf die Mikrogrenzschichtdicke, optimale Oberflächenrauigkeit erhöht die Blasenwachstumsrate und die Blasengröße. Dieses Ergebnis ist bisher einzigartig bei der Untersuchung der Einzelblasendynamik beim Blasensieden. Die Expansion der Blasenbasis und der Blasenwachstumsrate erreicht ein Maximum für das gemittelte Quadrat der Rauigkeit (Sq) im Bereich zwischen 0,156 und 0,202 m für Naturzirkulation. Für aufwärtsgerichtetes Strömungssieden war die Expansion der Blasenbasis und die Blasenwachstumsrate für Sq-Werte zwischen 0,108 und 0,218 m maximal. Der Blasenablösedurchmesser wurde für einen großen Bereich der Hauptströmungsgeschwindigkeiten und Wärmestromedichte gemittelt. Das Maximum des mittleren Ablösedurchmessers wurde für die Oberfläche mit einem Wert von Sq = 0,218 m erreicht. Die Oberflächenrauigkeit erweitert die Wärmeübertragungsoberfläche neben der Blasenbasis. Der Einfluss der Oberflächenrauigkeitshöhe auf die Blasen hängt von der Mikrogrenzschichtdicke sowie vom Blasenbasisradius ab. Das Modell der Mikrogrenzschichtdicke von Cooper und Lloyd [1] und die konzeptionelle Idee zur Störung der Mikrogrenzschicht durch die Rautiefe von Sriraman [2] wurden analysiert. Es wurde nachgewiesen, dass die Oberflächenrauigkeit die effektive Mikrogrenzschichtdicke und die dazugehörige Wärmeübertragung beeinflusst. (iv) Es wurden geringere Blasenwachstumsraten für höhere Hauptströmungs-geschwindigkeiten gemessen. Weiterhin reduzieren sich der Blasenablösedurchmesser sowie Ablöseperioden mit zunehmender Hauptströmungsgeschwindigkeit bei unterschiedlichen Wärmeoberflächencharakteristiken. Bei niedrigen Hauptströmungs-geschwindigkeiten im Bereich zwischen ungefähr 0,052 und 0,16 m/s reduziert sich der durchschnittliche Blasenablösedurchmesser deutlich. Die experimentellen Ergebnisse zeigen einen wesentlichen Einfluss der Oberflächenbeschaffenheit auf das Blasenwachstum und den Ablöseprozess beim Blasensieden. Um diesen Einfluss numerisch zu charakterisieren, wurde ein neues Blasenwachstumsmodel entwickelt. Existierende Blasenwachstumsmodelle berücksichtigen den umfangreichen Einfluss der Oberfläche des Heizers bisher nicht. Das vorgeschlagene Model bezieht die plausibelsten Mechanismen des Blasensiedens mit ein. Dazu zählen: Mikrogrenzschichtverdampfung im Bereich der Austrocknung, trägheits- und wärmediffusionskontrolliertes Blasenwachstum und Kondensation an der Blasenoberseite. Das Modell berücksichtigt, dass die überhitzte Flüssigkeitsschicht an der Heizerwand durch die wachsende Blase nach außen verdrängt wird und die so gestreckte Flüssigkeitsschicht einen Teil der Blase einhüllt. Kondensation erfolgt an der Blasengrenze, die in Kontakt mit der unterkühlten Flüssigkeit steht, und demzufolge mit der überhitzen Flüssigkeitsschicht nicht in Kontakt kommt. Das vorgeschlagene Blasenwachstumsmodel arbeitet mit drei Konstanten für die beschriebenen Wärmeübertragungsmechanismen beim Blasenwachstum. Dabei handelt es sich um eine Konstante für die effektive Mikrogrenzschichtdicke (Ceff ), eine weitere Konstante ́ für die Wärmediffusion hin zur Blase und der Trägheit sowie letztendlich einer Konstante S zur Abbildung des Kondensationswärmeübergangs, anhand der Beschreibung des Anteils der Blase, welcher in Kontakt mit der unterkühlten Flüssigkeit steht. Die effektive Mikrogrenzschichtdickenkonstante (Ceff) definiert den Einfluss der Heizoberflächencharakteristik auf die Verdampfung der Mikrogrenzschicht und somit die Blasenwachstumsrate beim Blasensieden. Die numerisch berechnete und experimentell gemessene Blasengröße wurde verglichen, um die Mikrogrenzschichtdickenkonstante Ceff zu definieren. Der Einfluss der Kondensation auf Ceff wurde geprüft.:Abstract v Kurzfassung vii Acknowledgements xiii Abbreviations and Symbols xv Chapter 1: Introduction and Motivation 1 1.1 General overview 1 1.2 Theoretical background 3 1.3 Objectives and outline of the thesis 7 Chapter 2: Fundamentals of Bubble Dynamics in Nucleate Boiling 9 2.1 Bubble growth in nucleate boiling 9 2.2 Bubble growth models 12 2.3 The physical process of bubble departure 16 2.4 Experimental investigations of bubble dynamics 20 2.4.1 Effects of heater surface characteristics 21 2.4.2 Effects of bulk liquid velocity 24 2.5 Chapter conclusion 26 Chapter 3: Heater Surface Preparation and Characterization 27 3.1 Surface properties 27 3.2 Surface preparation 29 3.2.1 Self-assembled monolayer coating 30 3.2.2 High-power pulsed laser irradiation 31 3.2.3 Wet-etching 32 3.3 Surface cleaning 32 3.4 Surface characterization 32 3.4.1 Wettability measurement 32 3.4.2 Roughness measurement 33 3.4.3 Analysis of surface characteristics 34 3.4.4 Uncertainty of surface parameters 38 3.5 Artificial cavity preparation 38 Chapter 4: Experimental Setup and Procedure 41 4.1 Natural circulation boiling (NCB 41 4.1.1 Experimental procedure and measurement techniques 41 4.1.2 Uncertainty analysis 44 4.2 Upward flow boiling (UFB) 45 4.2.1 Experimental procedure and measurement techniques 45 4.2.2 Uncertainty analysis 48 4.3 Image processing 50 Chapter 5: Experimental Results 53 5.1 Introduction to the analysis of the bubble dynamics 53 5.1.1 The bubble life cycle 53 5.1.2 Calculation of the bubble equivalent diameter 55 5.1.3 Bubble dynamics with the increase of heat flux 57 5.1.4 Qualitative assessment of the bubble dynamics for different parameters 60 5.2 Bubble dynamics 61 5.2.1 Effect of heater surface wettability 61 5.2.2 Effect of heater surface roughness 65 5.2.3 Effect of bulk liquid velocity 70 5.3 Bubble departure 76 5.3.1 Effect of heater surface wettablity 76 5.3.2 Effect of heater surface roughness 76 5.3.3 Effect of bulk liquid velocity 78 5.4 Chapter conclusion 79 Chapter 6: Analysis and Model Development 81 6.1 Numerical evaluation of the role of heater surface characteristics 81 6.1.1 Derivation of an improved bubble growth model 86 6.1.2 Calculation of Ceff 82 6.2 Effect of liquid velocity on the bubble growth 93 6.3 Improved modeling of bubble departure 95 6.3.1 Analysis of important parameters 95 6.3.2 Formulation of a bubble departure criterion 100 6.4 Chapter conclusion 102 Chapter 7: Summary and Outlook 105 Bibliography 109 List of Figures 121 List of Tables 127 Appendix: Surface Parameters and Profile 12

Investigations on the Effect of Heater Surface Characteristics on Bubble Dynamics in Subcooled Nucleate Boiling

Nucleating boiling is a repeating cycle of bubble initiation, growth and departure at many nucleation sites at the heated wall. Thereby, the bubble growth process significantly affects the dynamics of bubble departure. Experiments were performed to study the influence of heater surface characteristics, such as wettability and roughness, on single bubble growth and departure dynamics for natural circulation and upward flow boiling conditions. Self-assembled monolayer (SAM) coating, wet-etching and femtosecond pulsed laser treatment were used to alter the surface wettability and produce nano- and microstructures on stainless steel surfaces with a roughness in the range of micrometers. These surface preparation techniques allowed to separately quantify the effect of surface wettability and roughness on the bubble dynamics. The surface wettability and roughness are represented by the liquid contact angle hysteresis (θhys) and root mean square roughness of the surface (Sq). Boiling experiments were conducted at atmospheric pressure with degassed deionized water at low-subcooling. Stainless steel heater surfaces were vertically oriented during natural circulation boiling. In the experiments, bubbles were generated from an artificial nucleation cavity on the treated stainless steel heater surfaces. High-resolution optical shadowgraphy has been used to record the bubble generation, departure, sliding, detachment and inception of the next bubble. Higher bulk liquid velocity yielded smaller bubble departure diameters and slower bubble growth rates for all heater surface types. The effect of surface wettability on single bubble dynamics was studied for smooth surfaces with different liquid contact angle hysteresis. Low wetting surfaces yielded a greater bubble growth rate and departure diameter. The bubble growth rate and departure diameter were found maximum for an intermediate surface roughness Sq between 0.108 and 0.218 m. The corresponding roughness height is referred to as the ‘optimal roughness height’ in this work. Surface roughness was found very influential to the bubble growth and departure, which can be explained by considering its interaction with the microlayer underneath a bubble. The role of the heater surface parameters for the bubble growth was qualitatively assessed by evaluating the microlayer thickness constant C2. Hence, an improved bubble growth model was derived in this work. The bubble growth model was formulated on the basis of the evaporation of the microlayer beneath a bubble with the dryout area, inertia and heat diffusion controlled bubble growth and condensation at the bubble cap. The model can also predict the superheated liquid layer around a bubble which helps to determine the portion of a bubble that is in contact with the subcooled liquid. As bubble growth Abstract is highly dependent on the effective interactions of heater surface roughness and microlayer, a term Ceff was introduced in the bubble growth model. The effective microlayer thickness constant Ceff incorporates the impact of heater surface characteristics on the bubble growth process until the departure of a bubble. The bubble growth model was utilized in the analysis of high-resolution experimental data of steam bubble growth and the values of Ceff were calculated for different heater surface characteristics. The value of Ceff was found to decrease with the increase of bubble growth rate. A simplified model for the bubble departure criterion was derived from the expressions of forces which act on a nucleating bubble throughout its growth cycle. It was found that 90% of the departing bubbles satisfy the bubble departure criterion model with ±25% deviation. The knowledge gained from this work shall be particularly useful to improve nucleate boiling models for numerical simulations. The findings are also useful for designing heater surfaces in the future.:Abstract v Kurzfassung vii Acknowledgements xiii Abbreviations and Symbols xv Chapter 1: Introduction and Motivation 1 1.1 General overview 1 1.2 Theoretical background 3 1.3 Objectives and outline of the thesis 7 Chapter 2: Fundamentals of Bubble Dynamics in Nucleate Boiling 9 2.1 Bubble growth in nucleate boiling 9 2.2 Bubble growth models 12 2.3 The physical process of bubble departure 16 2.4 Experimental investigations of bubble dynamics 20 2.4.1 Effects of heater surface characteristics 21 2.4.2 Effects of bulk liquid velocity 24 2.5 Chapter conclusion 26 Chapter 3: Heater Surface Preparation and Characterization 27 3.1 Surface properties 27 3.2 Surface preparation 29 3.2.1 Self-assembled monolayer coating 30 3.2.2 High-power pulsed laser irradiation 31 3.2.3 Wet-etching 32 3.3 Surface cleaning 32 3.4 Surface characterization 32 3.4.1 Wettability measurement 32 3.4.2 Roughness measurement 33 3.4.3 Analysis of surface characteristics 34 3.4.4 Uncertainty of surface parameters 38 3.5 Artificial cavity preparation 38 Chapter 4: Experimental Setup and Procedure 41 4.1 Natural circulation boiling (NCB 41 4.1.1 Experimental procedure and measurement techniques 41 4.1.2 Uncertainty analysis 44 4.2 Upward flow boiling (UFB) 45 4.2.1 Experimental procedure and measurement techniques 45 4.2.2 Uncertainty analysis 48 4.3 Image processing 50 Chapter 5: Experimental Results 53 5.1 Introduction to the analysis of the bubble dynamics 53 5.1.1 The bubble life cycle 53 5.1.2 Calculation of the bubble equivalent diameter 55 5.1.3 Bubble dynamics with the increase of heat flux 57 5.1.4 Qualitative assessment of the bubble dynamics for different parameters 60 5.2 Bubble dynamics 61 5.2.1 Effect of heater surface wettability 61 5.2.2 Effect of heater surface roughness 65 5.2.3 Effect of bulk liquid velocity 70 5.3 Bubble departure 76 5.3.1 Effect of heater surface wettablity 76 5.3.2 Effect of heater surface roughness 76 5.3.3 Effect of bulk liquid velocity 78 5.4 Chapter conclusion 79 Chapter 6: Analysis and Model Development 81 6.1 Numerical evaluation of the role of heater surface characteristics 81 6.1.1 Derivation of an improved bubble growth model 86 6.1.2 Calculation of Ceff 82 6.2 Effect of liquid velocity on the bubble growth 93 6.3 Improved modeling of bubble departure 95 6.3.1 Analysis of important parameters 95 6.3.2 Formulation of a bubble departure criterion 100 6.4 Chapter conclusion 102 Chapter 7: Summary and Outlook 105 Bibliography 109 List of Figures 121 List of Tables 127 Appendix: Surface Parameters and Profile 129Der Blasenabriss von einer Keimstellenkavität ist ein komplexer Ablösemechanismus und spielt eine wichtige Rolle beim Wärmetransport. Zur Beschreibung der Blasendynamik sind Kenntnisse über den Blasenwachstumsprozess sowie die Vorhersage eines Kriteriums für die Blasenablösung erforderlich. In den existierenden Blasenwachstums- und Blasenablösungsmodellen wird die Oberflächencharakteristik des Heizers bisher nicht berücksichtigt. Im Rahmen dieser Promotion wurden Experimente durchgeführt, um den Einfluss der Heizeroberfläche und der Hauptströmungsgeschwindigkeit auf diese Parameter für eine vertikale Heizfläche zu untersuchen. Hierbei wurden das Naturkonvektionssieden und das aufwärtsgerichtete Strömungssieden betrachtet. Die Experimente wurden mit vollentsalztem Wasser bei einer Unterkühlung zwischen 1,68 und 4,00 K bei Atmosphärendruck und einem aus Edelstahl gefertigten Heizer durchgeführt, dessen Oberfläche anhand der Parameter Oberflächenrauigkeit und Benetzbarkeit charakterisiert ist. Unterschiedliche Oberflächenbearbeitungstechniken, wie Beschichtung durch Self-Assembled Monolayer (SAM), Nass-Ätzen und Hochleistungspuls-Laserbestrahlung wurden genutzt, um die Oberflächenbenetzung und –rauigkeit zu modifizieren. Der Unterschied zwischen dem gemessenen Fortschritts- (θadv) und Rückzugskontaktwinkel (θrec) der Flüssigkeit wird als Flüssigkeitskontaktwinkelhysterese (θhys) bezeichnet und beschreibt die Oberflächenbenetzbarkeit. Die Oberflächenrauigkeit wurde durch ein Konfokal-Mikroskop bestimmt und durch das gemittelte Quadrat der Rauigkeit (Sq) und den Maximalwert der Rauigkeit (St) definiert. Insgesamt wurden 18 unterschiedliche Heizoberflächen mit einer Größe von 130 x 20 mm² untersucht. Davon kamen jeweils die Hälfte für das Naturkonvektionssieden bzw. aufwärtsgerichtetes Strömungssieden zur Anwendung. Der Einfluss der Oberflächenbenetzbarkeit auf die Blasendynamik wurde für polierte Oberflächen (Sq 0,01 μm) analysiert. Die Wirkung der Oberflächenrauigkeit auf die Blasendynamik wurde für konstante Flüssigkeitskontaktwinkelhysteresen von 40,05°±1,5° und 59,97°±1,5° für Naturzirkulation und Strömungssieden untersucht. Eine künstliche zylindrische Kavität mit einer Fläche von 1963,5 m² und einer Tiefe von 50 m wurde mittels Mikrolaser in die Heizoberflächen eingebracht, um die Blasen in einer spezifischen Position zu erzeugen. Während des Naturkonvektionssiedens betrug die Wärmestromdichte 19,22 bis 30,29 kW/m². Bei den Experimenten mit aufwärtsgerichtetem Strömungssieden wurde die Hauptströmungsgeschwindigkeit im Bereich von 0,052 bis 0,183 m/s variiert und eine Appendix: Surface Parameters and Profile Wärmestromdichte zwischen 39,41 und 45,47 kW/m² aufgeprägt. Daraus resultierten insgesamt 87 Experimentalserien. Um den Blasenlebenszyklus zu erfassen, wurde hochauflösende Bildgebungstechnik verwendet. Mit der Bildverarbeitungssoftware ImageJ wurden die erfassten Videos weiterverarbeitet. Die Temperatur der Hauptströmung wurde mit Typ-K Thermoelementen gemessen. Die zeit- und ortsgemittelten Heizerwandtemperaturen wurden für die Naturzirkulation durch Infrarotthermografie und für das aufwärtsgerichtete Strömungssieden durch Typ-K Thermoelemente erfasst. Die mittlere Flüssigkeitsgeschwindigkeit wurde bei der Naturzirkulation mittels Particle Image Velocimetry (PIV) und beim Strömungssieden mittels Coriolis-Durchflusszähler bestimmt. Eine hochauflösende optische Schattenbildtechnik diente zur Aufzeichnung der Hauptphasen des Blasenlebenszyklus: Blasenerzeugung, Blasenwachstum, Blasenablösung, Blasengleiten und Blasenabriss. In dieser Arbeit wurden die der Blasenablösung vorrausgehenden Phasen untersucht. Blasenhöhe, Blasenbreite, Blasenbasisdurchmesser und Schwerpunkt der Blase wurden mit Hilfe der Bildverarbeitung ermittelt. Der blasenäquivalente Durchmesser wurde mittels des geometrischen Mittelwertes, der Blasenbreite und der Blasenhöhe berechnet. Basierend auf den Messdaten können folgende Erkenntnisse für das Blasenwachstum und den Blasenablösemechanismus postuliert werden: (i) Eine höhere Wärmeströmedichte führen zu größen Blasen und kürzeren Wachstumsperioden. Der Einfluss der Oberflächenbenetzbarkeit und der Oberflächenrauigkeit auf die Blasendynamik zeigt ähnliche Tendenzen für Naturkonvektion und aufwärtsgerichtetes Strömungssieden. (ii) Eine höhere Flüssigkeitskontaktwinkelhysterese führt zu einer schnelleren Expansion der Blasenbasis und zu einem schnellern Blasenwachstum. Für gut benetzbare Oberflächen bewegt sich der Blasenschwerpunkt schneller entlang der Strömungsrichtung. Für Oberflächen mit geringer Benetzbarkeit ist die Blasengröße vor der Blasenablösung größer und die Ablöseperiode länger. Der mittlere Blasenablösedurchmesser für unterschiedliche Hauptströmungsgeschwindigkeiten der Flüssigkeit erhöht sich von 0,75 auf 1,75 mm bei zunehmender Flüssigkeitskontaktwinkelhysterese von 42,32° auf 62,30°. (iii) Eine, bezogen auf die Mikrogrenzschichtdicke, optimale Oberflächenrauigkeit erhöht die Blasenwachstumsrate und die Blasengröße. Dieses Ergebnis ist bisher einzigartig bei der Untersuchung der Einzelblasendynamik beim Blasensieden. Die Expansion der Blasenbasis und der Blasenwachstumsrate erreicht ein Maximum für das gemittelte Quadrat der Rauigkeit (Sq) im Bereich zwischen 0,156 und 0,202 m für Naturzirkulation. Für aufwärtsgerichtetes Strömungssieden war die Expansion der Blasenbasis und die Blasenwachstumsrate für Sq-Werte zwischen 0,108 und 0,218 m maximal. Der Blasenablösedurchmesser wurde für einen großen Bereich der Hauptströmungsgeschwindigkeiten und Wärmestromedichte gemittelt. Das Maximum des mittleren Ablösedurchmessers wurde für die Oberfläche mit einem Wert von Sq = 0,218 m erreicht. Die Oberflächenrauigkeit erweitert die Wärmeübertragungsoberfläche neben der Blasenbasis. Der Einfluss der Oberflächenrauigkeitshöhe auf die Blasen hängt von der Mikrogrenzschichtdicke sowie vom Blasenbasisradius ab. Das Modell der Mikrogrenzschichtdicke von Cooper und Lloyd [1] und die konzeptionelle Idee zur Störung der Mikrogrenzschicht durch die Rautiefe von Sriraman [2] wurden analysiert. Es wurde nachgewiesen, dass die Oberflächenrauigkeit die effektive Mikrogrenzschichtdicke und die dazugehörige Wärmeübertragung beeinflusst. (iv) Es wurden geringere Blasenwachstumsraten für höhere Hauptströmungs-geschwindigkeiten gemessen. Weiterhin reduzieren sich der Blasenablösedurchmesser sowie Ablöseperioden mit zunehmender Hauptströmungsgeschwindigkeit bei unterschiedlichen Wärmeoberflächencharakteristiken. Bei niedrigen Hauptströmungs-geschwindigkeiten im Bereich zwischen ungefähr 0,052 und 0,16 m/s reduziert sich der durchschnittliche Blasenablösedurchmesser deutlich. Die experimentellen Ergebnisse zeigen einen wesentlichen Einfluss der Oberflächenbeschaffenheit auf das Blasenwachstum und den Ablöseprozess beim Blasensieden. Um diesen Einfluss numerisch zu charakterisieren, wurde ein neues Blasenwachstumsmodel entwickelt. Existierende Blasenwachstumsmodelle berücksichtigen den umfangreichen Einfluss der Oberfläche des Heizers bisher nicht. Das vorgeschlagene Model bezieht die plausibelsten Mechanismen des Blasensiedens mit ein. Dazu zählen: Mikrogrenzschichtverdampfung im Bereich der Austrocknung, trägheits- und wärmediffusionskontrolliertes Blasenwachstum und Kondensation an der Blasenoberseite. Das Modell berücksichtigt, dass die überhitzte Flüssigkeitsschicht an der Heizerwand durch die wachsende Blase nach außen verdrängt wird und die so gestreckte Flüssigkeitsschicht einen Teil der Blase einhüllt. Kondensation erfolgt an der Blasengrenze, die in Kontakt mit der unterkühlten Flüssigkeit steht, und demzufolge mit der überhitzen Flüssigkeitsschicht nicht in Kontakt kommt. Das vorgeschlagene Blasenwachstumsmodel arbeitet mit drei Konstanten für die beschriebenen Wärmeübertragungsmechanismen beim Blasenwachstum. Dabei handelt es sich um eine Konstante für die effektive Mikrogrenzschichtdicke (Ceff ), eine weitere Konstante ́ für die Wärmediffusion hin zur Blase und der Trägheit sowie letztendlich einer Konstante S zur Abbildung des Kondensationswärmeübergangs, anhand der Beschreibung des Anteils der Blase, welcher in Kontakt mit der unterkühlten Flüssigkeit steht. Die effektive Mikrogrenzschichtdickenkonstante (Ceff) definiert den Einfluss der Heizoberflächencharakteristik auf die Verdampfung der Mikrogrenzschicht und somit die Blasenwachstumsrate beim Blasensieden. Die numerisch berechnete und experimentell gemessene Blasengröße wurde verglichen, um die Mikrogrenzschichtdickenkonstante Ceff zu definieren. Der Einfluss der Kondensation auf Ceff wurde geprüft.:Abstract v Kurzfassung vii Acknowledgements xiii Abbreviations and Symbols xv Chapter 1: Introduction and Motivation 1 1.1 General overview 1 1.2 Theoretical background 3 1.3 Objectives and outline of the thesis 7 Chapter 2: Fundamentals of Bubble Dynamics in Nucleate Boiling 9 2.1 Bubble growth in nucleate boiling 9 2.2 Bubble growth models 12 2.3 The physical process of bubble departure 16 2.4 Experimental investigations of bubble dynamics 20 2.4.1 Effects of heater surface characteristics 21 2.4.2 Effects of bulk liquid velocity 24 2.5 Chapter conclusion 26 Chapter 3: Heater Surface Preparation and Characterization 27 3.1 Surface properties 27 3.2 Surface preparation 29 3.2.1 Self-assembled monolayer coating 30 3.2.2 High-power pulsed laser irradiation 31 3.2.3 Wet-etching 32 3.3 Surface cleaning 32 3.4 Surface characterization 32 3.4.1 Wettability measurement 32 3.4.2 Roughness measurement 33 3.4.3 Analysis of surface characteristics 34 3.4.4 Uncertainty of surface parameters 38 3.5 Artificial cavity preparation 38 Chapter 4: Experimental Setup and Procedure 41 4.1 Natural circulation boiling (NCB 41 4.1.1 Experimental procedure and measurement techniques 41 4.1.2 Uncertainty analysis 44 4.2 Upward flow boiling (UFB) 45 4.2.1 Experimental procedure and measurement techniques 45 4.2.2 Uncertainty analysis 48 4.3 Image processing 50 Chapter 5: Experimental Results 53 5.1 Introduction to the analysis of the bubble dynamics 53 5.1.1 The bubble life cycle 53 5.1.2 Calculation of the bubble equivalent diameter 55 5.1.3 Bubble dynamics with the increase of heat flux 57 5.1.4 Qualitative assessment of the bubble dynamics for different parameters 60 5.2 Bubble dynamics 61 5.2.1 Effect of heater surface wettability 61 5.2.2 Effect of heater surface roughness 65 5.2.3 Effect of bulk liquid velocity 70 5.3 Bubble departure 76 5.3.1 Effect of heater surface wettablity 76 5.3.2 Effect of heater surface roughness 76 5.3.3 Effect of bulk liquid velocity 78 5.4 Chapter conclusion 79 Chapter 6: Analysis and Model Development 81 6.1 Numerical evaluation of the role of heater surface characteristics 81 6.1.1 Derivation of an improved bubble growth model 86 6.1.2 Calculation of Ceff 82 6.2 Effect of liquid velocity on the bubble growth 93 6.3 Improved modeling of bubble departure 95 6.3.1 Analysis of important parameters 95 6.3.2 Formulation of a bubble departure criterion 100 6.4 Chapter conclusion 102 Chapter 7: Summary and Outlook 105 Bibliography 109 List of Figures 121 List of Tables 127 Appendix: Surface Parameters and Profile 12