SEARCH · Engineering Papers
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.
Steady state mercury evaporation experiment to measure transport coefficient in liquid-vapor interface boundary condition study
Irreversible thermodynamics of liquid-vapor phase change process
Thermodynamic rate equations for boundary analysis at liquid-vapor interface of mercury
To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized this is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory's Marmot application for phase-field simulation of nuclear materials. Simulations of bubble growth in the high burnup region during steady-state conditions showed that initially overpressurized bubbles decreased in pressure during growth, but still remained above equilibrium pressure. During a LOCA transient, simulations of bubbles in the high burnup region showed that bubble size did not change significantly, and the pressure increase due to the transient was calculated and passed to a phase-field model of fracture.
Boiling is known to be a very efficient mode of heat transfer, and as such, it is employed in component cooling and in various energy-conversion systems. In space, boiling heat transfer may be used in thermal management, fluid handling and control, power systems, and on-orbit storage and supply systems for cryogenic propellants and life-support fluids. Recent interest in the exploration of Mars and other planets and in the concept of in situ resource utilization on the Martian and Lunar surfaces highlights the need to understand how gravity levels varying from the Earth's gravity to microgravity (1g = or > g/g(sub e) = or > 10(exp -6)g) affect boiling heat transfer. Because of the complex nature of the boiling process, no generalized prediction or procedure has been developed to describe the boiling heat transfer coefficient, particularly at reduced gravity levels. Recently, Professor Vijay K. Dhir of the University of California at Los Angeles proposed a novel building-block approach to investigate the boiling phenomena in low-gravity to microgravity environments. This approach experimentally investigates the complete process of bubble inception, growth, and departure for single bubbles formed at a well-defined and controllable nucleation site. Principal investigator Professor Vijay K. Dhir, with support from researchers from the NASA Glenn Research Center at Lewis Field, is performing a series of pool boiling experiments in the low-gravity environments of the KC 135 microgravity aircraft s parabolic flight to investigate the inception, growth, departure, and merger of bubbles from single- and multiple-nucleation sites as a function of the wall superheat and the liquid subcooling. Silicon wafers with single and multiple cavities of known characteristics are being used as test surfaces. Water and PF5060 (an inert liquid) were chosen as test liquids so that the role of surface wettability and the magnitude of the effect of interfacial tension on boiling in reduced gravity can be investigated.
The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.
The previously developed dynamic theory of growth and detachment of vapour bubbles under conditions of nucleate pool boiling is applied to study motion and deformation of a bubble evolving at a single nucleation site. The bubble growth is presumed to be thermally controlled, and two components of heat transfer to the bubble are accounted of: the one from the bulk of surrounding liquid and the one due to heat conduction across a liquid microlayer formed underneath the bubble. Bubble evolution is governed by the buoyancy and an effective surface tension force, both the forces making the bubble centre of mass move away from the wall and, thus, assisting its detachment. Buoyancy-controlled and surface-tension-controlled regimes are considered separately in a meticulous way. The duration of the whole process of bubble evolution till detachment, the rate of growth, and the bubble departure size are found as functions of time and physical and operating parameters. Some repeatedly observed phenomena, such as an influence of gravity on the growth rate, are explained. Inferences of the model agree qualitatively with available experimental evidence, and conclusions pertaining to the dependence on gravity of the bubble radius at detachment and the whole time of the bubble development when being attached to the wall are confirmed quantitatively.
A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.
A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.
The radius-time history of a gas bubble located concentrically within a spherical liquid drop in a space laboratory is analyzed within the framework of the quasi-stationary approximation. Illustrative results are calculated from the theory which demonstrate interesting qualitative features. For instance, when a pure gas bubble dissolves within a liquid drop in an environment containing the same gas and some inert species, the dissolution can be more or less rapid than that in an unbounded liquid depending on the initial relative size of the drop. Further, given a similar growth situation, indefinite growth is not possible, and the bubble will initially grow, but always dissolve in the end.
Our hypothesis is that metabolic gases play a role in the initial explosive growth phase of bubble formation during hypobaric exposures. Models that account for optimal internal tensions of dissolved gases to predict the probability of occurrence of venous gas emboli were statistically fitted to 426 hypobaric exposures from National Aeronautics and Space Administration tests. The presence of venous gas emboli in the pulmonary artery was detected with an ultrasound Doppler detector. The model fit and parameter estimation were done by using the statistical method of maximum likelihood. The analysis results were as follows. 1) For the model without an input of noninert dissolved gas tissue tension, the log likelihood (in absolute value) was 255.01. 2) When an additional parameter was added to the model to account for the dissolved noninert gas tissue tension, the log likelihood was 251.70. The significance of the additional parameter was established based on the likelihood ratio test (P < 0.012). 3) The parameter estimate for the dissolved noninert gas tissue tension participating in bubble formation was 19. 1 kPa (143 mmHg). 4) The additional gas tissue tension, supposedly due to noninert gases, did not show an exponential decay as a function of time during denitrogenation, but it remained constant. 5) The positive sign for this parameter term in the model is characteristic of an outward radial pressure of gases in the bubble. This analysis suggests that dissolved gases other than N2 in tissues may facilitate the initial explosive bubble-growth phase.
Activation growth, and detachment of boiling bubbles in water from artificial nucleation sites of known geometry and size drilled in polished metal surface
In this work a building block type of approach is used so that a basic understanding of the processes that contribute to nucleate boiling heat fluxes under microgravity conditions can be developed. This understanding will lead to development of a mechanistic model for nucleate boiling heat transfer which could eventually be used as a design tool in space applications. Task Description Task 1: Fabrication of the Experimental Setup. Under this task, the test section and liquid holding and viewing chambers will be fabricated. Artificial cylinder cavities will be formed on silicon wafers. A single cavity and two or four cavities, with a prescribed spacing and size, will be formed. The desired nucleation wall superheat will be used to determine the size of the mouth of the cavities. Task 2: Experiments. The basic experiments for flow and temperature field around single and multiple (2 or 4 separated or merged bubbles growing on downward facing or inclined surfaces) will be carried out under normal gravity conditions. The experiments will be conducted at one atmosphere pressure, but liquid subcooling will be varied from 0 to 30C. Water and PF-5050 will be used as test liquids. Task 3: Analytical/Numerical Models. In this task, transient temperature and flow field in vapor and liquid will be determined during growth of a single bubble. Analysis will include the evolution of the vapor-liquid interface and development of microlayer underneath the bubbles. For merged bubbles, detailed calculations of flow and temperature field will be carried out for transient shapes of vapor stems supporting a large bubble and the corresponding evaporation rate. Flow and temperature field for a bubble sliding along a heated wall will also be determined. Microgravity conditions will be simulated and a framework of a numerical tool for prediction of nucleate boiling heat fluxes under microgravity conditions will be developed. Task 4: Experiments in a KC-135. To understand bubble growth and detachment behavior of single or large merged bubbles, boiling experiments will be conducted under low gravity (10-2 g) conditions of the aircraft. In these experiments, 'designed' surfaces will be used. Visual observations and heat transfer data will be taken, but holography will not be used. The apparatus used for laboratory experiments will also be employed for experiments in the aircraft. Task 5: Experiments in the Space Shuttle. Effort will be devoted for defining a boiling experiment to be conducted on a 'designed' surface. The experiment will provide microgravity data on bubble growth and departure. These data are needed for development of a credible model for nucleate boiling heat fluxes under microgravity conditions. The heat transfer data will also be obtained and will be used to validate the models.
Incipience, growth and detachment of boiling bubbles in saturated distilled water from artificial nucleation sites of given geometry and size in heated metal surface
Micro scale boiling behavior in the vicinity of graphite micro-fiber tips on the coppergraphite composite boiling surfaces is investigated. It is discovered that a large number of micro bubbles are formed first at the micro scratches and cavities on the copper matrix in pool boiling. In virtue of the non-wetting property of graphite, once the growing micro bubbles touch the graphite tips, the micro bubbles are sucked by the tips and merged into larger micro bubbles sitting on the tips. The micro bubbles grow rapidly and coalesce to form macro bubbles, each of which sitting on several tips. The growth processes of the micro and macro bubbles are analyzed and formulated followed by an analysis of bubble departure on the composite surfaces. Based on these analyses, the enhancement mechanism of the pool boiling heat transfer on the composite surfaces is clearly revealed. Experimental results of pool boiling heat transfer both for water and Freon-113 on the composite surfaces convincingly demonstrate the enhancement effects of the unique structure of Cu-Gr composite surfaces on boiling heat transfer.
The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.
The formation and growth processes of a bubble in the vicinity of graphite micro-fiber tips on metal-graphite composite boiling surfaces and their effects on boiling behavior are investigated. It is discovered that a large number of micro bubbles are formed first at the micro scratches and cavities on the metal matrix in pool boiling. By virtue of the non-wetting property of graphite, once the growing micro bubbles touch the graphite tips, the micro bubbles are sucked by the tips and merged into larger micro bubbles sitting on the end of the tips. The micro bubbles grow rapidly and coalesce to form macro bubbles, each spanning several tips. The necking process of a detaching macro bubble is analyzed. It is revealed that a liquid jet is produced by sudden break-off of the bubble throat. The composite surfaces not only have higher temperatures in micro- and macrolayers but also make higher frequency of the bubble departure, which increase the average heat fluxes in both the bubble growth stage and in the bubble departure period. Based on these analyses, the enhancement mechanism of pool boiling heat transfer on composite surfaces is clearly revealed.
Alpha-radiation damage is known to degrade materials exposed to harsh radiation environments via swelling and embrittlement by nano-scale helium bubbles. How these bubbles grow under varying temperature is essential to understanding and predicting material failure and optimization. Using in situ X-ray scattering, we probed ≈1012 bubbles before, during and after annealing at different temperatures in He-implanted aluminum to reveal a critical growth activation temperature that depends strongly on the helium implantation conditions. An analysis of the data aided by our molecular dynamics simulations identifies multiple temperature-dependent bubble growth mechanisms above and below the critical growth activation temperature. Moreover, the irreversible bubble growth that occurs after annealing for ≈4 h was evaluated using dislocation dynamics and found to negligibly affect the yield strength but may affect subsequent bubble growth at temperature. These in situ results provide valuable insight into bubble growth never observed before and reveal, for the first time, bubble shrinkage after cooling.