TU Darmstadt / ULB / TUprints

Bubble Nucleation in Pulsating Heat Pipes

Frey, Yannick Fabian (2024)
Bubble Nucleation in Pulsating Heat Pipes.
Technische Universität Darmstadt
doi: 10.26083/tuprints-00028599
Ph.D. Thesis, Primary publication, Publisher's Version

[img] Text
2024_10_27_Dissertation_Frey.pdf
Copyright Information: CC BY 4.0 International - Creative Commons, Attribution.

Download (15MB)
Item Type: Ph.D. Thesis
Type of entry: Primary publication
Title: Bubble Nucleation in Pulsating Heat Pipes
Language: English
Referees: Stephan, Prof. Dr. Peter ; Marengo, Prof. Dr. Marco
Date: 10 December 2024
Place of Publication: Darmstadt
Collation: xiv, 124 Seiten
Date of oral examination: 10 July 2024
DOI: 10.26083/tuprints-00028599
Abstract:

In this work, experimental, and numerical investigations are presented that aim to allow further insights into the operational mechanisms of pulsating heat pipes. The focus is on the link between necessary wall superheat for bubble nucleation and thermal resistance. Therefore, three transparent experimental setups were manufactured for testing that were geometrically identical but had different evaporator roughness manufactured by milling, sand blasting, and glass bead blasting techniques. The overall thermal resistance between evaporator and condenser was measured with R1233zd(E) as a working fluid and used for the evaluation of thermal performance. The results showed that decreasing filling ratio, increasing heat flow, operating temperature, and evaporator roughness influenced the thermal resistance beneficially. The influence of surface roughness was most pronounced for low heat flow. From confocal measurements of the evaporator surface structure and a newly introduced evaluation algorithm, bubble radius (cavity size) distribution profiles could be estimated for each surface structure. This allowed to apply the pool boiling nucleation equation to calculate theoretical required wall superheats for nucleation for each data point of the measurements. Based on the results, a relative nucleation threshold (RNT) was introduced that represents the quotient of theoretical required wall superheat for nucleation to the temperature difference from evaporator to condenser available due to thermal conduction in an inactive PHP. It was found that the thermal resistance of the PHP correlates directly with RNT, and that the correlation is valid for multiple heat flow, surface roughness values, and operating temperatures. The coefficients of the correlation depended on the filling ratio. Using a high-speed camera, the flow inside the channels of the evaporator section was recorded. In combination with an optical bubble detection and motion tracking algorithm, average bubble nucleation rates (number of bubble generations per second) in the evaporator area for the analyzed time intervals were derived. The data at 70% filling ratio showed a connection between RNT and the bubble nucleation rate. The maximum bubble nucleation rate was located at an RNT value of 0.3. At RNT close to zero and high values of RNT, the bubble nucleation rate was approaching zero. For RNT>0.3, the behavior could be explained as higher RNT values lead to a lower probability of bubble nucleation and, thus, a reduced bubble nucleation rate. The decrease in nucleation rate from RNT 0.3 to RNT zero could not be quantitatively explained by the experiments. Flow observations suggested that low values of RNT cause a decrease of liquid mass in the evaporator. As bubbles can only be nucleated if liquid is present, the bubble nucleation rate is decreased here, even though RNT is also reduced. Subsequently, 1D-simulations of the PHP were conducted to gain further insights. To obtain correct physical understanding of PHP function from the simulations, they had to first be validated in the intended parameter range. Consequently, the experimental setups were reconstructed in the simulations. It could be shown that the simulation yielded accurate thermal resistances at 70% filling ratio and various condenser temperatures, heat flow, and surface roughness values. However, at lower filling ratios the simulation did not produce accurate results. Therefore, only simulations at 70% filling ratio were further analyzed. It could be shown that heat transfer coefficients between the fluid and the channel walls, the flow oscillation frequencies, and the liquid mass fractions inside the evaporator area correlated directly with RNT, strengthening the plausibility of the correlation found in the experiments. By presuming that bubbles can only be nucleated, if liquid is present within the evaporator, the liquid mass fractions found in the simulations allowed to quantitatively explain the bubble nucleation rates in the experiments. Compendiously, it could be shown that the required wall superheat for nucleation has a major influence on the PHP performance. To allow for the use of the correlation with different channel geometries, the correlation was empirically extended by further simulation data to provide a feasible PHP design tool. For this, simulations with various numbers of channels and evaporator lengths were conducted and a modified version of the correlation was introduced, accounting for the total channel length within the evaporator area.

Alternative Abstract:
Alternative AbstractLanguage

In dieser Arbeit werden experimentelle und numerische Untersuchungen präsentiert, die weiterführende Einblicke in die Funktionsweise eines pulsierenden Wärmerohrs gewähren. Der Schwerpunkt liegt dabei auf dem Zusammenhang zwischen der für Nukleation notwendigen Wandüberhitzung und dem thermischen Widerstand. Daher wurden drei transparente, geometrisch identische, Versuchsaufbauten angefertigt, die sich lediglich durch deren Oberflächenrauigkeit im Verdampferbereich unterscheiden. Diese wurden durch Fräsen sowie Sand- und Glasperlenstrahlen hergestellt. Um die thermische Performance bewerten zu können, wurde der thermische Widerstand zwischen Verdampfer und Kondensator des pulsierenden Wärmerohrs mit R1233zd(E) als Kältemittel gemessen. Geringer Flüssigfüllstand, hohe Wärmeströme, Fluidtemperatur und Oberflächenrauigkeit verringerten dabei den thermischen Widerstand. Der Einfluss der Oberflächenrauigkeit war besonders signifikant bei kleinen Wärmeströmen. Die Oberflächenstruktur im Verdampferbereich wurde durch Aufnahmen mit einem Konfokal-Mikroskop sowie eines neuartigen Auswertealgorithmus charakterisiert, dadurch konnten Blasenradienverteilungen (Kavitätsgrößenverteilungen) aus den Aufnahmen abgeleitet werden. Auf dieser Basis konnten mittlere notwendige Wandüberhitzungen für das Blasensieden für jeden aufgezeichneten Messpunkt anhand der Siedeverzugsgleichung bestimmt werden. Eine neu eingeführte Größe, die relative Nukleationshürde (RNT), konnte somit bestimmt werden. Sie stellt den Quotienten der theoretisch notwendigen Wandüberhitzung zur Temperaturdifferenz dar, die aufgrund von Wärmeleitung in einem inaktiven pulsierenden Wärmerohr verfügbar ist. Es konnte gezeigt werden, dass der thermische Widerstand des pulsierenden Wärmerohrs direkt mit RNT korreliert. Die Korrelation ist für unterschiedliche Wärmeströme, Fluidtemperaturen und Verdampferrauigkeiten gültig. Die Koeffizienten der Korrelation hängen von dem Flüssigfüllstand des pulsierenden Wärmerohrs ab. Mithilfe einer Hochgeschwindigkeitskamera konnten Aufnahmen der Fluidströmung in den Kanälen des Verdampferbereichs aufgezeichnet werden. In Kombination mit einem neuartigen Blasendetektions- und Tracking-Algorithmus konnten somit mittlere Blasenentstehungsraten (Anzahl neu entstandener Blasen pro Sekunde) im Verdampferbereich in den untersuchten Zeitintervallen bestimmt werden. Bei einem Flüssigkeitsfüllstand von 70% wurde ein Zusammenhang zwischen der Blasenentstehungsrate und RNT entdeckt. Die maximale Blasenentstehungsrate tritt bei RNT=0,3 auf. Für große RNT-Werte oder gegen RNT=0 nähert sich die Blasenentstehungsrate null an. Für RNT>0,3 lässt sich dieser Trend dadurch erklären, dass ein höherer RNT-Wert zu einer geringeren Wahrscheinlichkeit für Blasenbildung und damit zu einer niedrigeren Blasenentstehungsrate führt. In den experimentellen Untersuchungen konnte jedoch nicht quantitativ bestimmt werden, wie der Abfall der Blasenentstehungsrate gegen RNT=0 zu erklären ist. Beobachtungen des Strömungsverhaltens legen jedoch nahe, dass sich für niedrige RNT-Werte weniger Flüssigphase im Verdampfer befindet. Da Blasen nur entstehen können, wenn Flüssigphase im Verdampfer vorhanden ist, entstehen hier weniger Blasen, obwohl RNT reduziert wird. Um weitere Einblicke in das Verhalten eines pulsierenden Wärmerohrs zu bekommen, wurden daraufhin 1D-Simulationen durchgeführt. Um ein korrektes physikalisches Verständnis anhand der Simulationen ableiten zu können, mussten diese zunächst mittels Versuchsdaten im entsprechenden Parameterbereich validiert werden. Daher wurden die Randbedingungen aus den durchgeführten Versuchen in den Simulationen rekonstruiert. Die Simulation zeigte dabei für 70% Flüssigkeitsfüllstand ausreichend genaue Ergebnisse für den thermischen Widerstand bei Variationen des Wärmestroms, der Fluidtemperatur und der Oberflächenrauigkeit. Andere Füllstände lieferten jedoch keine ausreichende Genauigkeit in der Simulation, daher wurden für die weiterführende Analyse der Simulationsergebnisse nur Berechnungen bei 70% Füllstand verwendet. Anhand der Simulationsdaten konnte aufgezeigt werden, dass die Wärmeübergangskoeffizienten zwischen der Kanalwand und dem Fluid, die Oszillationsfrequenz der Strömung, sowie der Anteil an Flüssigphase im Verdampferbereich von RNT abhängen. Dies stärkt die Plausibilität der experimentell ermittelten Korrelation. Unter der Annahme, dass Blasen nur entstehen können, wenn Flüssigphase im Verdampfer vorhanden ist, konnten die simulativ ermittelten Flüssigphasenanteile im Verdampfer genutzt werden, um die Blasenentstehungsraten in den Versuchen quantitativ zu erklären. Zusammenfassend konnte gezeigt werden, dass die für das Blasensieden notwendige Wandüberhitzung einen großen Einfluss auf die Leistung eines pulsierenden Wärmerohrs hat. Um die Korrelation auch für verschiedene Kanalgeometrien nutzen zu können, wurde sie anhand von Simulationsdaten empirisch erweitert, sodass der Auslegungsprozess für ein pulsierendes Wärmerohr vereinfacht werden kann. Dazu wurden Simulationen bei unterschiedlichen Verdampferlängen und unterschiedlicher Anzahl an Kanälen durchgeführt. Anhand der Gesamtkanallänge im Verdampferbereich wurde eine modifizierte Korrelation eingeführt.

German
Status: Publisher's Version
URN: urn:nbn:de:tuda-tuprints-285990
Classification DDC: 500 Science and mathematics > 500 Science
600 Technology, medicine, applied sciences > 600 Technology
Divisions: 16 Department of Mechanical Engineering > Institute for Technical Thermodynamics (TTD)
Date Deposited: 10 Dec 2024 13:51
Last Modified: 11 Dec 2024 10:57
URI: https://tuprints.ulb.tu-darmstadt.de/id/eprint/28599
PPN: 524505829
Export:
Actions (login required)
View Item View Item