Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bubble growth”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Keyhole fluctuation and pore formation mechanisms during laser powder bed fusion additive manufacturing

Keyhole porosity is a key concern in laser powder-bed fusion (LPBF), potentially impacting component fatigue life. However, some keyhole porosity formation mechanisms, e.g., keyhole fluctuation, collapse and bubble growth and shrinkage, remain unclear. Using synchrotron X-ray imaging we reveal keyhole and bubble behaviour, quantifying their formation dynamics. The findings support the hypotheses that: (i) keyhole porosity can initiate not only in unstable, but also in the transition keyhole regimes created by high laser power-velocity conditions, causing fast radial keyhole fluctuations (2.5-10 kHz); (ii) transition regime collapse tends to occur part way up the rear-wall; and (iii) immediately after keyhole collapse, bubbles undergo rapid growth due to pressure equilibration, then shrink due to metal-vapour condensation. Concurrent with condensation, hydrogen diffusion into the bubble slows the shrinkage and stabilises the bubble size. The keyhole fluctuation and bubble evolution mechanisms revealed here may guide the development of control systems for minimising porosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗

Time and Space Resolved Wall Temperature Measurements during Nucleate Boiling with Constant Heat Flux Boundary Conditions

The lack of temporally and spatially resolved measurements under nucleate bubbles has complicated efforts to fully explain pool-boiling phenomena. The objective of this current work is to acquire time and space resolved temperature distributions under nucleate bubbles on a constant heat flux surface. This was performed using a microheater array with 100 micron resolution that allowed effectively simultaneous measurements of surface temperature while supplying a constant dissipative heat flux. This data is then correlated with high speed (> 1000Hz) visual recordings of the bubble growth and departure from the heater surface acquired from below and from the side of the heater. The data indicate that a significant source of energy during bubble nucleation and initial growth is the superheated layer around the bubble. Bubble coalescence was not observed to decrease surface temperature as significantly as bubble departure from the surface. Since bubble departure is typically followed by a sharp increase in the heater surface temperature, it is surmised that the departing bubble effectively removes the superheated layer, allowing a high local heat transfer rate with the bulk fluid through transient conduction/micro-convection during rewetting.

Myers, Jerry G.↗

Behavior of bubbles in glassmelts. III - Dissolution and growth of a rising bubble containing a single gas

Finite difference solutions of the mass transport equations governing the dissolution (growth) of a rising gas bubble, containing a single gas, in a glassmelt were obtained. These solutions were compared with those obtained from an approximate procedure for a range of the controlling parameters. Applications were made to describe various aspects of O2 and CO2 gas-bubble behavior in a soda-lime-silicate melt.

Onorato, P. I. K.↗

The bubbly-slug transition in a boiling two-phase flow under microgravity

A theory is presented to describe, in reduced gravity flow boiling, the transition from bubbly two-phase flow to slug flow. It is shown that characteristics of the bubbly flow and the transition were controlled by the mechanism of vapor bubble growth dynamics. By considering in nucleate boiling, behavior of vapor bubbles at departure from a heated surface a condition required for transition was determined. Although required, this condition alone could not ensure coalescence of bubbles to cause the transition to slug two-phase flow. The condition leading to coalescence, therefore, was obtained by examining oscillations of vapor bubbles following their departure from the heated surface. The predicted transition conditions were compared with the prediction and test data reported for adiabatic reduced gravity two-phase flow, and good qualitative agreement was found.

Kiper, Ali M.↗

BISON microstructure-based pulverization criterion in high burnup structure

To improve the economics of commercial nuclear power production, utilities are seeking to increase the allowable burnup limit of UO$_2$ fuel. One of the main factors that contributes to the current burnup limit of 62 GWd/MTU in commercial light water reactors (LWRs) is the risk of fine fragmentation or pulverization during a loss of coolant accident (LOCA). Pulverization primarily occurs at high burnups, especially when the high burnup structure (HBS) has formed. To allow the industry to pursue increased burnup and develop mitigation strategies, it is essential to have improved capability to predict the onset of pulverization. However, the mechanism of pulverization is not well understood, and the existing predictive capabilities implemented in the BISON fuel performance code are empirical in nature. In this report, mesoscale simulations are used to improve understanding of the formation mechanism of the HBS and how it responds during a LOCA transient, and inform development of a BISON pulverization criterion. A phase-field model was used to simulate the evolution of bubble pressure as a result of HBS formation. The simulations showed that gas atoms diffuse from grain interiors to the new grain boundaries created during HBS formation, and diffuse rapidly along these grain boundaries to reach existing bubbles. This causes an increase in bubble pressure in existing bubbles, leading to bubble growth during steady-state operation. To simulate the response of HBS bubbles to a LOCA transient, a newly developed phase-field model was used; in agreement with preliminary results from FY20, bubble size did not change significantly during the duration of the transient. A phase-field fracture model was used to study fragmentation patterns in the HBS, including using input from the phase-field model as initial conditions. Phase-field fracture simulations were used to determine a pulverization criterion for BISON. A function for the critical pressure for pulverization to occur was fit to data from the phase-field fracture simulations, and this function was implemented as a material property in BISON. For comparison, an analytical criterion for pulverization was developed and implemented within the same material property.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Synergistic foam stabilization and transport improvement in simulated fractures with polyelectrolyte complex nanoparticles: Microscale observation using laser etched glass micromodels

Inaccessibility to direct pore scale observation in hydrocarbon recovery of tight shale formations poses a great challenge to water-energy nexus initiatives and necessitates the use of high throughput technologies to emulate environmentally friendly processes. Herein, we employ a precise glass micromodel fabrication and visualization method to isolate the supercritical CO 2 bubbles surrounded by CO 2 -water lamella prepared in saline produced water stabilized with molecular complexation of zwitterionic surfactants (ZS) and polyelectrolyte complex nanoparticles (PECNP). The Selective Laser Enhanced Etching (SLE) technique was selected for micromodel simulation of high-pressure flow. Two representative designs, (1) fracture/micro-crack network 28 and (2) fracture/matrix were etched on fused silica glass with a laser printing machine and scCO 2 foam was injected to study the foamability, propagation, stability, and fluid loss properties. The highly monodispersed and uniformly distributed array of scCO 2 bubbles were detected in flow of scCO 2 foam in highly saline brine containing ionic complexes of positively charged PECNPs and ZS, whereas foam flow with the lamella containing ZS in fractures offered a noticeably large and polydisperse array of scCO 2 bubbles. scCO 2 bubble motion and deformation were traced, and local description of foam flow was visually examined. The confined array of scCO 2 bubbles stabilized by ZS in microcracks was affected by bubble growth and coalescence, whereas the super-populated array of monodispersed scC O2 bubbles with lamella containing complexes of PECNP and ZS were able to fill the channels with stable configurations within the timeframe of comparative stability measurements. The ability of complex fluid to prevent the formation damage was evaluated through fluid loss visualization in micromodels. Probing scCO 2 foam transport in homogenous porous media revealed smaller volume leak-off for scCO 2 foam containing PECNP-ZS ionic complexes.

04 OIL SHALES AND TAR SANDS↗

Hydrate Formation on Marine Seep Bubbles and the Implications for Water Column Methane Dissolution

Abstract Methane released from seafloor seeps contributes to a number of benthic, water column, and atmospheric processes. At seafloor seeps within the methane hydrate stability zone, crystalline gas hydrate shells can form on methane bubbles while the bubbles are still in contact with the seafloor or as the bubbles begin ascending through the water column. These shells reduce methane dissolution rates, allowing hydrate‐coated bubbles to deliver methane to shallower depths in the water column than hydrate‐free bubbles. Here, we analyze seafloor videos from six deepwater seep sites associated with a diverse range of bubble‐release processes involving hydrate formation. Bubbles that grow rapidly are often hydrate‐free when released from the seafloor. As bubble growth slows and seafloor residence time increases, a hydrate coating can form on the bubble's gas‐water interface, fully coating most bubbles within ∼10 s of the onset of hydrate formation at the seafloor. This finding agrees with water‐column observations that most bubbles become hydrate‐coated after their initial ∼150 cm of rise, which takes about 10 s. Whether a bubble is coated or not at the seafloor affects how much methane a bubble contains and how quickly that methane dissolves during the bubble's rise through the water column. A simplified model shows that, after rising 150 cm above the seafloor, a bubble that grew a hydrate shell before releasing from the seafloor will have ∼5% more methane than a bubble of initial equal volume that did not grow a hydrate shell after it traveled to the same height.

54 ENVIRONMENTAL SCIENCES↗

Bubble Formation at a Submerged Orifice in Reduced Gravity

The dynamic regime of gas injection through a circular plate orifice into an ideally wetting liquid is considered, when successively detached bubbles may be regarded as separate identities. In normal gravity and at relatively low gas flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole time of its evolution. If the flow rate exceeds a certain threshold value, another stage of the detachment process takes place in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical stem that lengthens as the bubble rises above the plate. The bubble shape resembles then that of a mushroom and the upper envelope continues to grow until the gas supply through the stem is completely cut off. Such a stage is always present under conditions of sufficiently low gravity, irrespective of the flow rate. Two major reasons make for bubble detachment: the buoyancy force and the force due to the momentum inflow into the bubble with the injected gas. The former force dominates the process at normal gravity whereas the second one plays a key role under negligible gravity conditions. It is precisely this fundamental factor that conditions the drastic influence on bubble growth and detachment that changes in gravity are able to cause. The frequency of bubble formation is proportional to and the volume of detached bubbles is independent of the gas flow rate in sufficiently low gravity, while at normal and moderately reduced gravity conditions the first variable slightly decreases and the second one almost linearly increases as the flow rate grows. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Buyevich, Yu A.↗

Effects of Gravity on Bubble Formation at a Plate Orifice

A model of the dynamic regime of gas injection through a submerged plate orifice into an ideally wetting liquid is developed in the circumstance when successively detached bubbles may be regarded as independent objects. Two major factors favor bubble detachment: the buoyancy force and a force due to the momentum inflow into the bubble with injected gas. In normal and moderately reduced gravity, the first factor dominates. At relatively low flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole period of its evolution till detachment. If the flow rate exceeds a critical value, another stage of bubble evolution occurs in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical gaseous stem that lengthens as the bubble rises above the plate. The bubble continues to grow until the gas supply through the stem is completely cut off, after which back flow of gas into the stem from the bubble becomes possible. In microgravity, the second factor prevails, and the latter stage is always present irrespective of the flow rate. However, the gas envelope rises and the stem lengthens very slowly. This difference in the underlying physical mechanisms provides for key properties of bubble growth and detachment being drastically different in appreciable and sufficiently reduced gravity. The frequency of bubble formation slightly decreases with and the detachment volume is almost proportional to the gas flow rate in the first case, in accordance with familiar relations. In the second case, the first variable is proportional to the flow rate whereas the second one is independent of it. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.

Webbon, Bruce W.↗

Potential Fifty Percent Reduction in Saturation Diving Decompression Time Using a Combination of Intermittent Recompression and Exercise

Conventional saturation decompression protocols use linear decompression rates that become progressively slower at shallower depths, consistent with free gas phase control vs. dissolved gas elimination kinetics. If decompression is limited by control of free gas phase, linear decompression is an inefficient strategy. The NASA prebreathe reduction program demonstrated that exercise during O2 prebreathe resulted in a 50% reduction (2 h vs. 4 h) in the saturation decompression time from 14.7 to 4.3 psi and a significant reduction in decompression sickness (DCS: 0 vs. 23.7%). Combining exercise with intermittent recompression, which controls gas phase growth and eliminates supersaturation before exercising, may enable more efficient saturation decompression schedules. A tissue bubble dynamics model (TBDM) was used in conjunction with a NASA exercise prebreathe model (NEPM) that relates tissue inert gas exchange rate constants to exercise (ml O2/kg-min), to develop a schedule for decompression from helium saturation at 400 fsw. The models provide significant prediction (p < 0.001) and goodness of fit with 430 cases of DCS in 6437 laboratory dives for TBDM (p = 0.77) and with 22 cases of DCS in 159 altitude exposures for NEPM (p = 0.70). The models have also been used operationally in over 25,000 dives (TBDM) and 40 spacewalks (NEPM). The standard U.S. Navy (USN) linear saturation decompression schedule from saturation at 400 fsw required 114.5 h with a maximum Bubble Growth Index (BGI(sub max)) of 17.5. Decompression using intermittent recompression combined with two 10 min exercise periods (75% VO2 (sub peak)) per day required 54.25 h (BGI(sub max): 14.7). Combined intermittent recompression and exercise resulted in a theoretical 53% (2.5 day) reduction in decompression time and theoretically lower DCS risk compared to the standard USN decompression schedule. These results warrant future decompression trials to evaluate the efficacy of this approach.

Gernhardt, Michael I.↗

Application of Dynamical System Scaling to Bubble Dynamics

Dynamical System Scaling (DSS) provides a useful method for analyzing, categorizing and scaling time-dependent processes. A key feature of DSS is the temporal displacement rate, D, which relates the natural process time to the reference clock time. Because of its property of being invariant under a two-parameter affine transformation, it provides an underlying basis for process scaling. Application of DSS to bubble dynamics serves to elucidate the physical significance of the temporal displacement rate and its role in scaling a variety of bubble dynamics processes. This paper shows that the temporal displacement rate consists of the sum of dimensionless groups that govern the bubble dynamics. Preserving the temporal displacement rate for a prototype and a scaled model results in similitude of the time-dependent normalized bubble size, growth rates and interfacial accelerations for different fluid conditions. For vapor bubble growth in a superheated liquid, it is shown that inertial controlled bubble growth occurs when D=0 and thermally controlled bubble growth occurs when D=1 .

dynamical system scaling, SMR, bubble dynamics, DS↗

Implications of Microstructure in Helium-Implanted Nanocrystalline Metals

Helium bubbles are known to form in nuclear reactor structural components when displacement damage occurs in conjunction with helium exposure and/or transmutation. If left unchecked, bubble production can cause swelling, blistering, and embrittlement, all of which substantially degrade materials and—moreover—diminish mechanical properties. On the mission to produce more robust materials, nanocrystalline (NC) metals show great potential and are postulated to exhibit superior radiation resistance due to their high defect and particle sink densities; however, much is still unknown about the mechanisms of defect evolution in these systems under extreme conditions. Here, the performances of NC nickel (Ni) and iron (Fe) are investigated under helium bombardment via transmission electron microscopy (TEM). Bubble density statistics are measured as a function of grain size in specimens implanted under similar conditions. While the overall trends revealed an increase in bubble density up to saturation in both samples, bubble density in Fe was over 300% greater than in Ni. To interrogate the kinetics of helium diffusion and trapping, a rate theory model is developed that substantiates that helium is more readily captured within grains in helium-vacancy complexes in NC Fe, whereas helium is more prone to traversing the grain matrices and migrating to GBs in NC Ni. Our results suggest that (1) grain boundaries can affect bubble swelling in grain matrices significantly and can have a dominant effect over crystal structure, and (2) an NC-Ni-based material can yield superior resistance to irradiation-induced bubble growth compared to an NC-Fe-based material and exhibits high potential for use in extreme environments where swelling due to He bubble formation is of significant concern.

36 MATERIALS SCIENCE↗

Diffusion in Melts and Magmas

Diffusion results from random motion of particles and entities. Diffusion in melts and magmas is due to thermally excited random motion of atoms, ions, and clusters, and plays a critical role in magmatic and volcanic processes. In melts and magmas, diffusion is one of the two mechanisms of mass transfer; the other being bulk flow (referred to as convection or advection). When both are present, diffusion refers to the dispersive motion relative to the mean bulk flow in a given reference frame (Richter et al. 1998). Diffusion plays critical roles in controlling magma mixing (Watson 1982; Koyaguchi 1985, 1989; Lesher 1994; Huber et al. 2009; Guo and Zhang 2020), mineral growth and dissolution rates in magmas (e.g., Watson 1982; Harrison and Watson 1983; Zhang et al. 1989; Newcombe et al. 2014; Macris et al. 2018), bubble growth and dissolution rate in magmas (Sparks 1978; Proussevitch and Sahagian 1998; Liu and Zhang 2000; Zhang 2013), and elemental and isotope fractionation during mineral 2 growth and dissolution (Jambon 1980; Richter et al. 1999, 2003; Watson and Muller 2009; Chopra et al. 2012; Watkins et al. 2014, 2017; Holycross et al. 2016, 2018). As a result, diffusion also plays an essential role in explosive volcanic eruptions and magma crystallization. Furthermore, diffusion has important applications in geospeedometry (Lasaga 1983, 1998; Zhang 1994, 2008; Trail et al. 2016; Zhang and Xu 2016). Experimental investigation of diffusion in geologically relevant silicate melts began to flourish in the 1970’s when micro-analytical measurements of diffusion profiles became available. In addition to the vast number of papers published since then, numerous books and reviews are available for diffusion in silicate melts. Hofmann et al. (1974) edited a book titled “Geochemical Transport and Kinetics” published by Carnegie Institution of Washington. This was the first landmark book summarizing the field. Lasaga and Kirkpatrick (1981) edited a book “Kinetics of Geochemical Processes” as volume 8 of the Reviews in Mineralogy (later becoming Reviews in Mineralogy and Geochemistry) series. Zhang and Cherniak (2010) edited “Diffusion in Minerals and Melts” as volume 72 of Reviews of Mineralogy and Geochemistry series, in which one chapter focused on diffusion theory, five chapters on diffusion in silicate melts (Behrens 2010; Lesher 2010; Liang 2010; Zhang and Ni 2010; Zhang et al. 2010), and other chapters were on experimental, analytical, and computational methods, and diffusion in minerals. Several textbooks covered the principles and applications of diffusion theories (Kirkaldy and Young 1987; Shewmon 1989; Cussler 1997; Lasaga 1998; Zhang 2008; Vrentas and Vrentas 2016), and two classic books covered the mathematics of diffusion (Carslaw and Jaeger 1959; Crank 1975). In preparing for this review chapter, we thought carefully about what to cover for this vast field, and decided to briefly go through the fundamentals of diffusion (more complete review can be found in Chakraborty 1995; Zhang 2008, 2010) and solutions to often-encountered 3 diffusion problems, and then focus on post-2010 diffusion studies on silicate melts and magmas. Here, melts refer to (mostly natural) silicate liquid, and magmas refer to crystal-bearing and/or bubble-bearing melts in which the continuous phase is the melt. There is a large body of work on diffusion in glasses, especially in the materials science literature, which is not covered in this review.

Youxue Zhang↗

Evolution of Vapor Bubbles Nucleation Sites in Low Gravity

When liquid is expelled by a vapor bubble growing at a nucleation site on a superheated surface, a thin microlayer underneath the bubble is left behind. It is evaporated from the free microlayer surface that provides for bubble growth. The average thickness of the microlayer determining the evaporation rate increases with time if the latter does not exceed a threshold value associated with the burn-out crisis. The bubble is described as a spherical segment with its flattened part adjoining the microlayer. This introduces two independent variables - the radius of the spherical part of the bubble surface and the polar angle that defines the relative area of the flattened part. They are to be found out from a set of two strongly nonlinear equations resulting from mass and momentum conservation laws. The first one depends on both microlayer thickness and nonmonotonously changing bubble base area. The second involves two major factors favoring bubble detachment - the buoyancy and a force due to the initial momentum of vapor input into the bubble. The former force depends on gravity whereas the latter one does not. It is why the limiting regimes of bubble evolution that correspond to normal or moderately reduced gravity and to microgravity feature drastically different properties. In the first case, the buoyancy dominates and the bubble evolves in such a manner as to become a full sphere at a moment that can be viewed as that of detachment. The detachment volume grows as gravity decreases. In the second case, the buoyancy is negligible and the bubble stays near the surface, while its volume continues to increase for a sufficiently long time. The findings are discussed in connection with experimental data obtained under different gravity conditions, some unpublished experiments being included. They help to understand why the pool boiling heat transfer coefficient frequently increases as gravity falls down and eventually vanishes.

Buyevich, Yu A.↗

Studies on pressure response of gas bubbles contributions of condensed droplets in bubbles generated by a uniform nucleation

The response of a tiny gas bubble under reduced pressure is investigated in its relation to cavitation. Equations of motion are formulated for gas mixtures inside the bubble and numerical calculations performed for several examples. The conclusions are as follows: (1) at the onset of bubble growth, the gas mixture inside it adiabatically expands and the temperature decreases. Condensed droplets appear inside the gas mixture due to a uniform nucleation and the temperature recovers, thus the motion of the bubble is apparently isothermal; (2) the evaporation and condensation coefficient largely affects bubble motions (maximum radius, period and rate of attenuation of the bubble oscillation) including the uniform contraction; (3) the oscillation period of the bubble is longer as the equilibrium bubble radius is larger when the surrounding pressure decreases stepwise. In this circumstance the temperature inside the bubble is kept constant due to condensation evaporation phenomena and is nearly isothermal; and (4) when the surrounding pressure decreases in a stepwise fashion, the critical pressure bubble radius relation becomes closer to that for the isothermal process if the bubble radius is larger than 8 microns.

Matsumoto, Y.↗

Design of an Improved Heater Array to Measure Microscale Wall Heat Transfer

An improved array of microscale heaters is being developed to measure the heat transfer coefficient at many points underneath individual bubbles during boiling as a function of space and time. This heater array enables the local heat transfer from a surface during the bubble growth and departure process to be measured with very high temporal and spatial resolution, and should allow better understanding of the boiling heat transfer mechanisms by pin-pointing when and where in the bubble departure cycle large amounts of wall heat transfer occur. Such information can provide much needed data regarding the important heat transfer mechanisms during the bubble departure cycle, and can serve as benchmarks to validate many of the analytical and numerical models used to simulate boiling. The improvements to the heater array include using a silicon-on-quartz substrate to reduce thermal cross-talk between the heaters, decreased space between the heaters, increased pad sizes on the heaters, and progressive heater sizes. Some results using the present heater array are discussed.

Kim, Jungho↗