Effect of microgravity on flow boiling heat transfer of liquid hydrogen in transportation pipes

The flow boiling phenomenon of liquid hydrogen (LH2) during transportation in microgravity is very different from that under terrestrial condition. In this study, a saturated flow boiling of LH2 in a horizontal tube has been simulated under microgravity condition using coupled level-set and volume of fluid method. The validation of the developed model shows good agreement with the experimental data from the literature. The changes of heat fluxes and pressure drops under different gravitational accelerations were analyzed. And, the variation of heat fluxes with different wall superheat and contact angle were compared between microgravity (10⁻⁴g) and normal gravity (1g) condition. Also, the influence of surface tension were studied under microgravity. The numerical results indicate that the heat flux decrease with the decrement of gravitational acceleration. And the heat transfer ratio decrease with the increment of wall superheat in the nucleate boiling regime. The heat transfer slightly reduce when considering surface tension. In addition, the changes of contact angle have a more significant impact on heat transfer under microgravity condition.     

Mixed Thermocapillary and Forced Convection Heat Transfer Around a Hemispherical Bubble in a Miniature Channel—A 3D Numerical Study

It has been established that for certain conditions, such as microgravity boiling, thermocapillary Marangoni flow has associated with it a significant enhancement of heat transfer. Typically, this phenomenon was investigated for the idealized case of an isolated and stationary bubble resting atop a heated solid that is immersed in a semi-infinite quiescent fluid or within a two-dimensional cavity. This article presents a three-dimensional numerical study that investigates the influence of thermal Marangoni convection on the fluid dynamics and heat transfer around a bubble during laminar flow of water in a minichannel. This mixed thermocapillary and forced convection problem is investigated for channel liquid inlet velocity of 0.01 m/s to 0.03 m/s and Marangoni numbers in the range of 10 to 300 under microgravity conditions. Three-dimensional effects become particularly important on the side and rear regions of the bubble. The thermocapillary forces accelerate the flow along almost the entire bubble interface. The hot core fluid from the heated bottom wall region is forced inward and propelled upward into the thermocapillary jet above the bubble. It can be quantified that the influence of thermocapillary flow on heat transfer enhancement shows an average increase by 40% at the downstream of the bubble and by 60% at the front and rear regions. This heat transfer enhancement depends mainly on the temperature differential as the driving potential for thermocapillary flow and bulk liquid velocity.     

Thermal Bubble Dynamics Under the Effects of an Acoustic Field

In this research, dynamics of a micro thermal bubble existing in an acoustic field have been studied by a high-speed camera and a micro temperature sensor. The micro thermal bubble was generated by a micro heater, which was fabricated by the MEMS (micro-electro-mechanical system) technique and packed into a transparent mini chamber. The acoustic field inside the chamber was generated by a piezoelectric plate that was attached on the top side of the chamber's wall. Compared with micro thermal bubble dynamics in normal conditions, several different bubble dynamic phenomena in acoustic conditions have been found, such as bubble departure and attraction around the heater, bubble oscillating in the liquid volume, etc. By theoretical analysis, the main mechanism of bubble movements is attributed to the balance between Marangoni force and acoustic force. All these bubble dynamic phenomena improve the liquid convective flow and enhance the heat and mass transfer. Thus, this investigation about acoustic thermal bubble dynamics may find some potential applications in micro fluid devices for different functions, such as heat/mass transfer enhancement, micro electronic cooling, micro heater protection, etc. Temperature measurement in both normal conditions and acoustic conditions confirmed that the heat transfer was enhanced by the acoustic field.     

Interface Shape and Marangoni Effect Around aBubble Within the Thermal Boundary Layer

INTRoDUCTIoNDuetohighheattransferperformancecharacter-izedbysmalltemPeraturedifferencesandhighheatfluxes,transportprocesseswithphasechange,espe-ciallyboilingandcondensationprocessesarewidelyemployedinnumerousenergyconversionandtrans-portsystems,heatingand/orcoolingdevices,andaerospaceaPplications.Priortotheutilizationofboil-ingprocessesinspaceapplications,suchasspacecraftthermalcontrol,additionalunderstandingofboilingheattransferbehaviorisneeded.Becauselargedmer-encesekistinthefiuiddensiti…     

Microscale heterogeneous boiling on smooth surfaces—from bubble nucleation to bubble dynamics

In this investigation, boiling incipience and bubble dynamics on a microheater with a geometry of 100 μm × 100 μm fabricated with MEMS technology are evaluated using a high-speed digital camera. For the purpose of comparison with conventional boiling heat transfer, boiling incipience and bubble dynamics are also studied on a carefully selected microheater with a fabricated defect (i.e., a microcavity on the heater surface). Of industrial interest are the effects of dissolved gases on boiling incipience and bubble dynamics, which are also discussed in detail. The possible nucleation temperature (or incipience temperature) is analyzed and discussed from the perspective of the measured bulk temperature of the microheater and a 3D heat conduction numerical model. The time-resolved bubble dynamics (i.e., the bubble size evolution, interface velocity and interface acceleration) are all presented along with high-speed digital images. Based upon this investigation, it is clear that explosive boiling can take place on a smooth surface no matter how slow the heating rate, and dissolved gases have a significant influence on the incipience temperature and bubble behavior. Furthermore, this study illustrates that the classical kinetics of boiling can explain the explosive boiling occurring on a smooth surface in principle and can provide a useful guide for the design of microscale heat transfer and/or MEMS devices. Although unexpected, due to the gravitational effects, Marangoni flow on the vapor–liquid interface induced by the temperature gradient was also observed.     

Effect of gravity conditions on heat transfer and flow fluctuations of liquid hydrogen in a non-isothermal horizontal circular tube

Liquid hydrogen phase transition is a common phenomenon in space missions for space vehicles using low temperature liquid hydrogen as propellant. In this study, a numerical model with coupled RANS solver and VOF/Level-set method was used to simulate the liquid hydrogen phase transition in a non-isothermal horizontal circular tube under different gravity conditions (1g-10^?4 g). The gas phase hydrogen produced by evaporation of liquid hydrogen was calculated by Lee model. The statistics of the overall volume, heat flux, mass flow rate, mean velocity of gas phase hydrogen was carried out. The data results shown that the flow fluctuations were strongest under the gravity acceleration of 10^?1 g relative to other gravity conditions. The average bubble volume at 10^?1 g was the smallest, which was 11.58% smaller than that at 10^?3 g condition. The intermittent contact with the tube wall, which leaded to intermittent long bubble and flow resistance, was the main reason.     

Prediction of liquid hydrogen flow boiling critical heat flux condition under microgravity based on the wall heat flux partition model

Critical heat flux (CHF) of liquid hydrogen (LH2) flow boiling under microgravity is vital for designing space cryogenic propellant conveying pipe since the excursion of wall temperature may cause system failure. In this study, a two-dimensional axisymmetric model based on the wall heat flux partition (WHFP) model was proposed to predict the CHF condition under microgravity including the wall temperature and the CHF location. The proposed numerical model was validated to demonstrate a good agreement between the simulated and experimentally reported results. Then, the wall temperature distribution and the CHF location under different gravity conditions were compared. In addition, the WHFP and vapor-liquid distribution along the wall under microgravity were predicted and its difference with terrestrial gravity condition was also analysed and reported. Finally, the effects of flow velocity and inlet sub-cooling on the wall temperature distributions were analysed under microgravity and terrestrial gravity conditions, respectively. The results indicate that the CHF location moves upstream about 5.25 m from 1g to 10⁻⁴g since the void fraction near the wall reaches the breakpoint of CHF condition much earlier under the microgravity condition. Furthermore, the increase of the velocity and decrease of the sub-cooling have smaller effects on the CHF location during LH2 flow boiling under microgravity.     

Numerical analysis of the dynamics of moving vapor bubbles

Bubbles have been observed rapidly sweeping along very fine heated wires during subcooled nucleate boiling with jet flows emanating from the tops of the vapor bubbles. This paper analyzes the physical mechanisms driving the bubble and the jet flows from the tops of these moving bubbles. The flows are analyzed by numerically solving the governing equations for the velocity and temperature distributions around the bubble and the heated wire as the bubble moves along the wire. The bubble motion is due to the non-uniform temperature distribution in the liquid and in the wire caused by the bubble as it moves along the wire. The flow is driven by the horizontal Marangoni flow induced by the temperature difference across the bubble which thrusts the bubble forward. Comparisons with experimental observations suggest that the condensation heat transfer at the bubble interface is restricted by non-condensable gases that increases the surface temperature gradient and the resulting Marangoni flow.     

Convection and evaporation of microlayer underneath bubbles under microgravity

The multidimensional heat transfer and fluid flow in the microlayer region below a vapor bubble formed during boiling in microgravity are investigated by numerically solving the Navier–Stokes equations with the energy equation. The flow is driven by Marangoni flow due to the surface tension gradient along the bubble surface that results from the temperature gradient. The model also includes condensation and evaporation at the bubble surface. The flow field and heat transfer are calculated for microlayer thicknesses from 0.01 mm to 10 mm to investigate the effect of microlayer thickness. The results show that the velocities are small and have only a small effect on the temperature distribution as compared to the solution for pure conduction in the liquid. Natural convection is shown to have a negligible effect on heat transfer. For less than ideal evaporative heat transfer at the bubble interface, Marangoni convection caused the heat transfer to increase several percent. The flow in the microlayer is shown to agree with the lubrication analogy only for thin, relatively flat interfaces. © 2000 Scripta Technica, Heat Trans Asian Res, 30(1): 1–10, 2001     

Heat transfer for Marangoni‐driven boundary layer flow

Marangoni convection induced by variation of the surface tension with temperature along a surface influences crystal growth melts and other processes with liquid–vapor interfaces, such as boiling in both microgravity and normal gravity in some cases. This paper presents the Nusselt number for Marangoni flow over a flat surface calculated using a similarity solution for both the momentum equations and the energy equation assuming developing boundary layer flow along a surface. Solutions are presented for the surface velocity, the total flow rate, and the Nusselt number for various temperature profiles, Marangoni numbers, and Prandtl numbers. For large bubbles, the predicted boundary layer thickness would be less than the bubble diameter, so the curvature effects could be neglected and this analysis could be used as a first estimate of the effect of Marangoni flow around a vapor bubble. © 2002 Scripta Technica, Heat Trans Asian Res, 31(2): 105–116, 2002; DOI 10.1002/htj.10019     

Effect of gravity on cryogenic boiling heat transfer during tube quenching

Cryogenic forced convective boiling under terrestrial and microgravity conditions for the development of cryogenic fluid management on orbit is studied. The experiments are conducted in a low mass velocity region (100–300 kg/m² s) that is easily influenced by gravity, and fluid behavior observations and heat transfer measurements are performed simultaneously. These experiments aim at understanding the effect of gravitational acceleration on the relation between the flow behavior and thermal characteristics during the quenching of the tube by a cryogenic fluid. The heat transfer increases under microgravity conditions, and results from an increase in the quench front velocity.     

CFD simulation of heat transfer and phase change characteristics of the cryogenic liquid hydrogen tank under microgravity conditions

The heat transfer and phase change processes of cryogenic liquid hydrogen (LH₂) in the tank have an important influence on the working performance of the liquid hydrogen-liquid oxygen storage and supply system of rockets and spacecrafts. In this study, we use the RANS method coupled with Lee model and VOF (volume of fraction) method to solve Navier-stokes equations. The Lee model is adopted to describe the phase change process of liquid hydrogen, and the VOF method is utilized to calculate free surface by solving the advection equation of volume fraction. The model is used to simulate the heat transfer and phase change processes of the cryogenic liquid hydrogen in the storage tank with the different gravitational accelerations, initial temperature, and liquid fill ratios of liquid hydrogen. Numerical results indicate greater gravitational acceleration enhances buoyancy and convection, enhancing convective heat transfer and evaporation processes in the tank. When the acceleration of gravity increases from 10^?2 g0 to 10^?5 g0, gaseous hydrogen mass increases from 0.0157 kg to 0.0244 kg at 200s. With the increase of initial liquid hydrogen temperature, the heat required to raise the liquid hydrogen to saturation temperature decreases and causes more liquid hydrogen to evaporate and cools the gas hydrogen temperature. More cryogenic liquid hydrogen (i.e., larger the fill ratio) makes the average fluid temperature in the tank lower. A 12.5% reduction in the fill ratio resulted in a decrease in fluid temperature from 20.35 K to 20.15 K (a reduction of about 0.1%, at 200s).     

Studying the effect of the force convection boundary condition and radius magnetic field on fluid flow and heat transfer between coaxial pipes

An investigation of the heat transfer of Newtonian fluid flow through coaxial two pipes with variable radius ratio has been conducted with the boundary conditions of forced convection on the inner pipe walls and a radius magnetic field. This paper presents an exact analytical solution to the momentum equation and a novel semi-analytic collocation method for solving the full-term energy equation that takes Joule heating into account as well as viscous dissipation. Based on the results of the numerical fourth-order Runge–Kutta method, it was found that increasing the magnetic parameter decreased the amount of friction on the surface of the pipe walls and the rate of heat transfer. As the radius ratio of the two pipes increases, so does the skin friction and heat transfer rate on the internal pipe walls. As Eckert (Ec) and Prandtl (Pr) numbers increase, the mean temperature as well as the dimensionless temperature between the two pipes increases. The increase in Biot number (Bi) has the opposite impact on the mean temperature. As Ec, Pr, and Bi increase, so does the rate of heat transfer on the inner wall of the pipe.     

液层厚度比对环形腔内双层液体的热毛细对流稳定性的影响

为了了解微重力条件下、水平温度梯度作用时,上部为固壁的环形腔内双层流体系统中液层厚度比对流动稳定性的影响,采用隐式重启Arnoldi方法(IRAM)对环形池内5cSt硅油/HT-70双层流体的热对流过程进行了线性稳定性分析,获得了不同液层厚度比下的临界Marangoni数、临界波数、临界相速度,并通过计算特征向量,得到了临界Marangoni数附近液-液界面的热流体波形态。     

A Mathematical Model of an Oscillating Heat Pipe

A mathematical model predicting the oscillating motion in an oscillating heat pipe is developed. The model considers the system multidegree oscillation of vapor bubbles and liquid plugs, including the effects of filling ratio, operating temperature, gravitational force, and temperature difference between the evaporator and condenser. The model shows that the average velocity of liquid slugs is determined by the temperature difference between the evaporator and condenser. As the turn number increases, the temperature difference for the system to start the oscillating motion decreases. Increasing the bubble number will make the system more unstable and the system can be easily started up. The existence of gravity at the bottom heating mode will make the system easily produce the oscillating motion and decrease the temperature difference as well. Results presented here will assist in optimizing the heat transfer performance and provide a better understanding of heat transfer mechanisms occurring in the oscillating heat pipe.     

On the dynamics and heat transfer of bubble train in micro-channel flow boiling

The dynamics and heat transfer characteristics of flow boiling bubble train moving in a micro channel is studied numerically. The coupled level set and volume of fluid (CLSVOF) is utilized to track interface and a non-equilibrium phase change model is applied to calculate the interface temperature as well as heat flux jump. The working fluid is R134a and the wall material is aluminum. The fluid enters the channel with a constant mass flux (335 kg/m² 1 s), and the boundary wall is heated with constant heat flux (14 kW/m²). The growth of bubbles and the transition of flow regime are compared to an experimental visualization. Moreover, the bubble evaporation rate and wall heat transfer coefficient have been examined, respectively. Local heat transfer is significantly enhanced by evaporation occurring vicinity of interface of the bubbles. The local wall temperature is found to be dependent on the thickness of the liquid film between the bubble train and the wall.     

Mixed Convective Heat Transfer Due to Forced and Thermocapillary Flow Around Bubbles in a Miniature Channel: A 2D Numerical Study

Marangoni thermocapillary convection and its contribution to heat transfer during boiling has been the subject of some debate in the literature. Currently, for certain conditions, such as microgravity boiling, it has been shown that Marangoni thermocapillary convection has a significant contribution to heat transfer. Typically, this phenomenon is investigated for the idealized case of an isolated and stationary bubble resting on a heated surface, which is immersed in a semi-infinite quiescent fluid or within a two-dimensional cavity. However, little information is available with regard to Marangoni heat transfer in miniature confined channels in the presence of a cross flow. As a result, this article presents a two-dimensional (2D) numerical study that investigates the influence of steady thermal Marangoni convection on the fluid dynamics and heat transfer around a bubble during laminar flow of water in a miniature channel with the view of developing a refined understanding of boiling heat transfer for such a configuration. This mixed convection problem is investigated under microgravity conditions for channel Reynolds numbers in the range of 0 to 500 at liquid inlet velocities between 0.01 m/s and 0.0 5m/s and Marangoni numbers in the range of 0 to 17,114. It is concluded that thermocapillary flow may have a significant impact on heat transfer enhancement. The simulations predict an average increase of 35% in heat flux at the downstream region of the bubble, while an average 60% increase is obtained at the front region of the bubble where mixed convective heat transfer takes place due to forced and thermocapillary flow.     

On the bubble departure diameter and release frequency based on numerical simulation results

Heterogeneous boiling on a horizontal plate in stagnant and slowly flowing fluid is simulated using the lattice Boltzmann approach. The bubble departure diameter and release frequency are determined from the simulation results. It is found that the bubble departure diameter is proportional to g^?1/2 in a stagnant fluid and the release frequency scales with g^3/4, where g is the gravitational acceleration. Simulation results show no dependence between the bubble departure diameter and the static contact angle, but the bubble release frequency increases exponentially with the latter. Considering forced boiling, exponential relation is observed between the bubble departure diameter and the flow driving pressure gradient.     

环形池内具有自由表面的双层流体热毛细对流

采用数值模拟方法研究微重力条件下环形双层液体内存在水平温度梯度时的热毛细对流及其稳定性。流体为5cSt硅油/HT-70,外壁被加热、内壁被冷却,下固壁和上自由表面均绝热。结果表明:当Ma较小时,流动为稳定的轴对称流动;随着Ma和深宽比的增大,流动加强,等温线发生强烈的非线性变形;当Ma超过临界值后,流动转化为非稳定的多胞流动;随着Ma和深宽比的增大,速度振荡增大并向热壁方向运动,多胞流动结构占据区域拓展;流动转变的临界Ma随着深宽比的增大而减小。     

Numerical investigation of the bubble growth in horizontal rectangular microchannels

To explore the mechanism of flow boiling in microchannels, the processes of a single-vapor bubble evaporating and two lateral bubbles merging in a 2D microchannel are investigated. The temperature recovery model based on volume of fluid method is adopted to perform the flow boiling phenomena. The effects of wall superheat, Reynolds number, contact angle, surface tension, and two-bubble merger on heat transfer are discussed. Wall superheat dominates the bubble growth and is roughly proportional to wall heat flux. The update of velocity and temperature fields’ distribution in the channel increases with increasing inflow Reynolds number, which improves the wall heat flux markedly. Besides, the area of thin liquid film between the wall and the bubble is enlarged by reducing the contact angle, thus, expanding the wall heat flux several times compared with the single-phase cross section. However, variation of surface tension (0.0589, 0.1?N/m) is found to be insignificant.