Onset of flow oscillations in forced-flow subcooled boiling
Onset of flow oscillations of subcooled liquid in forced flow conditions in single tube boiler
SEARCH · Search NASA
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.
Onset of flow oscillations of subcooled liquid in forced flow conditions in single tube boiler
Onset of flow oscillations when subcooled liquid undergoes phase change under forced flow in single tube boiler
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.
The two-phase bubbly flow and boiling heat transfer in microgravity represents a substantial challenge to scientists and engineers and yet there is an urgent need to seek fundamental understanding in this area for future spacecraft design and space missions. At Washington State University, we have successfully designed, built and tested a 2.1 second drop tower with an innovation airbag deceleration system. Microgravity boiling experiments performed in our 0.6 second Drop Tower produced data flow visualizations that agree with published results and also provide some new understanding concerning flow boiling and microgravity bubble behavior. On the analytical and numerical work, the edge effects of finite divergent electrode plates on the forces experienced by bubbles were investigated. Boiling in a concentric cylinder microgravity and an electric field was numerically predicted. We also completed a feasibility study for microgravity boiling in an acoustic field.
Pressure drop oscillations in forced convection flow of Freon-11 near onset of boiling, noting flow rate for instability and limit cycles
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.
Available experimental data dealing with critical heat flux (CHF) of liquid hydrogen (LH2), liquid methane (LCH4), and liquid oxygen (LO2) in pool and flow boiling are compiled. The compiled data are compared with widely used correlations. Experimental pool boiling CHF data for the aforementioned cryogens are scarce. Based on only 25 data points found in five independent sources, the correlation of Sun and Lienhard (1970) is recommended for predicting the pool CHF of LH2. Only two experiments with useful CHF data for the pool boiling of LCH4 could be found. Four different correlations including the correlation of Lurie and Noyes (1964) can predict the pool boiling CHF of LCH4 within a factor of two for more than 70% of the data. Furthermore, based on the 19 data points taken from only two available sources, the correlation of Sun and Lienhard (1970) is recommended for the prediction of pool CHF of LO2. Flow boiling CHF data for LH2 could be found in seven experimental studies, five of them from the same source. Based on the 91 data points, it is suggested that the correlation of Katto and Ohno (1984) be used to predict the flow CHF of LH2. No useful data could be found for flow boiling CHF of LCH4 or LO2. The available databases for flow boiling of LCH4 and LO2 are generally deficient in all boiling regimes. This deficiency is particularly serious with respect to flow boiling.
Available experimental data dealing with critical heat flux (CHF) of liquid hydrogen (LH 2 ), liquid methane (LCH 4 ), and liquid oxygen (LO 2 ) in pool and flow boiling are compiled. The compiled data are compared with widely used correlations. Experimental pool boiling CHF data for the aforementioned cryogens are scarce. Based on only 25 data points found in five independent sources, the correlation of Sun and Lienhard (1970) is recommended for predicting the pool CHF of LH 2 . Only two experiments with useful CHF data for the pool boiling of LCH 4 could be found. Four different correlations including the correlation of Lurie and Noyes (1964) can predict the pool boiling CHF of LCH 4 within a factor of two for more than 70% of the data. Furthermore, based on the 19 data points taken from only two available sources, the correlation of Sun and Lienhard (1970) is recommended for the prediction of pool CHF of LO 2 . Flow boiling CHF data for LH 2 could be found in seven experimental studies, five of them from the same source. Based on the 91 data points, it is suggested that the correlation of Katto and Ohno (1984) be used to predict the flow CHF of LH 2 . No useful data could be found for flow boiling CHF of LCH 4 or LO 2 . The available databases for flow boiling of LCH 4 and LO 2 are generally deficient in all boiling regimes. This deficiency is particularly serious with respect to flow boiling.
The Two-phase Heat Transfer International Topical Team consists of researchers and members from various space agencies including ESA, JAXA, CSA, and RSA. This presentation included descriptions various fluid experiments either being conducted by or planned by NASA for the International Space Station in the areas of two-phase flow, flow boiling, capillary flow, and crygenic fluid storage.
This study defined the preliminary designs of flight software for the Space Shuttle Orbiter for three two-phase fluid research experiments: (1) liquid reorientation - to study the motion of liquid in tanks subjected to small accelerations; (2) pool boiling - to study low-gravity boiling from horizontal cylinders; and (3) flow boiling - to study low-gravity forced flow boiling heat transfer and flow phenomena in a heated horizontal tube. The study consisted of eight major tasks: reassessment of the existing experiment designs, assessment of the Spacelab facility approach, assessment of the individual carry-on approach, selection of the preferred approach, preliminary design of flight hardware, safety analysis, preparation of a development plan, estimates of detailed design, fabrication and ground testing costs. The most cost effective design approach for the experiments is individual carry-ons in the Orbiter middeck. The experiments were designed to fit into one or two middeck lockers. Development schedules for the detailed design, fabrication and ground testing ranged from 15 1/2 to 18 months. Minimum costs (in 1981 dollars) ranged from $463K for the liquid reorientation experiment to $998K for the pool boiling experiment.
The high power density of emerging electronic devices is driving the transition from remote cooling, which relies on conduction and spreading, to embedded cooling, which extracts dissipated heat on-site. Two-phase microgap coolers employ the forced flow of dielectric fluids undergoing phase change in a heated channel within or between devices. Such coolers must work reliably in all orientations for a variety of applications (e.g., vehicle-based equipment), as well as in microgravity and high-g for aerospace applications, but the lack of acceptable models and correlations for orientation- and gravity-independent operation has limited their use. Reliable criteria for achieving orientation- and gravity-independent flow boiling would enable emerging systems to exploit this thermal management technique and streamline the technology development process. As a first step toward understanding the effect of gravity in two-phase microgap flow and transport, in an earlier effort, the authors studied the effects of evaporator orientation, mass flux, and heat flux on flow boiling of HFE7100 in a 1.01 mm tall by 13.0 mm wide by 12.7 mm long microgap channel. Orientation-independence, defined as achieving similar critical heat fluxes, heat transfer coefficients, and flow regimes across orientations, was achieved for mass fluxes of 400 kg/sq.m-s and greater (corresponding to a Froude number of about 0.8). In the present effort, the authors have studied the effects of gravity, mass flux, and subcooling on flow boiling of HFE7100 in a 0.17 mm tall by 13.0 mm wide by 12.7 mm long microgap channel. The Flow Boiling in Microgap Coolers payload experienced about three minutes of weightlessness and shorter periods of high-g (up to about 5 g) during two recent flights aboard the Blue Origin New Shepard reusable launch vehicle. The results from the flight experiments will be presented and compared with published criteria for achieving gravity-independence.
Experiments were performed to examine the effects of body force on flow boiling critical heat flux (CHF). FC-72 was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface just prior to CHF. High-speed video imaging techniques were used to identify dominant CHF mechanisms corresponding to different flow orientations and liquid velocities. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed significant sensitivity to orientation for flow velocities below 0.2 m/s, where extremely low CHF values where measured, especially with downward-facing heated wall and downflow orientations. High flow velocities dampened the effects of orientation considerably. The CHF data were used to assess the suitability of previous CHF models and correlations. It is shown the Interfacial Lift-off Model is very effective at predicting CHF for high velocities at all orientations. The flooding limit, on the other hand, is useful at estimating CHF at low velocities and for downflow orientations. A new method consisting of three dimensionless criteria is developed for determining the minimum flow velocity required to overcome body force effects on near-saturated flow boiling CHF. Vertical upflow boiling experiments were performed in pursuit of identifying the trigger mechanism for subcooled flow boiling CHF. While virtually all prior studies on flow boiling CHF concern the prediction or measurement of conditions that lead to CHF, this study was focused on events that take place during the CHF transient. High-speed video imaging and photomicrographic techniques were used to record the transient behavior of interfacial features from the last steady-state power level before CHF until the moment of power cut-off following CHF. The video records show the development of a wavy vapor layer which propagates along the heated wall, permitting cooling prior to CHF only in wetting fronts corresponding to the wave troughs. Image analysis software was developed to estimate void fraction from the individual video images. The void fraction records for subcooled flow boiling show the CHF transient is accompanied by gradual lift-off of wetting fronts culminating in some maximum vapor layer mean thickness, following which the vapor layer begins to thin down as the transition to film boiling ensues. This study proves the Interfacial Lift-off Model, which has been validated for near-saturated flow boiling CHF, is equally valid for subcooled conditions.
Boiling flow instabilities - inlet pressure drop versus overall density ratio relationships
An experimental proceedure for examining flow boiling in low gravity environment is presented. The proceedure involves both ground based and KC-135 flight experiments. Two experimental apparati were employed, one for studying subcooled liquid boiling and another for examining saturated liquid boiling. For the saturated flow experiments, liquid nitrogen was used while freon 113 was used for the subcooled flow experiments. The boiling phenomenon was investigated in both cases using flow visualization techniques as well as tube wall temperature measurements. The flow field in both cases was established by injecting cold liquid in a heated tube whose temperature was set above the saturation values. The tubes were both vertically and horizontally supported with the liquid injected from the lower end of the tube. The results indicate substantial differences in the flow patterns established during boiling between the ground based, (1-g), experiments and the flight experiments, (low-g). These differences in the flow patterns will be discussed and some explanations will be offered.
Future manned missions to Mars are expected to greatly increase the space vehicle's size, weight, and heat dissipation requirements. An effective means to reducing both size and weight is to replace single-phase thermal management systems with two-phase counterparts that capitalize upon both latent and sensible heat of the coolant rather than sensible heat alone. This shift is expected to yield orders of magnitude enhancements in flow boiling and condensation heat transfer coefficients. A major challenge to this shift is a lack of reliable tools for accurate prediction of two-phase pressure drop and heat transfer coefficient in reduced gravity. Developing such tools will require a sophisticated experimental facility to enable investigators to perform both flow boiling and condensation experiments in microgravity in pursuit of reliable databases. This study will discuss the development of the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS), which was initiated in 2012 in collaboration between Purdue University and NASA Glenn Research Center. This facility was recently tested in parabolic flight to acquire condensation data for FC-72 in microgravity, aided by high-speed video analysis of interfacial structure of the condensation film. The condensation is achieved by rejecting heat to a counter flow of water, and experiments were performed at different mass velocities of FC-72 and water and different FC-72 inlet qualities. It is shown that the film flow varies from smooth-laminar to wavy-laminar and ultimately turbulent with increasing FC-72 mass velocity. The heat transfer coefficient is highest near the inlet of the condensation tube, where the film is thinnest, and decreases monotonically along the tube, except for high FC-72 mass velocities, where the heat transfer coefficient is enhanced downstream. This enhancement is attributed to both turbulence and increased interfacial waviness. One-ge correlations are shown to predict the average condensation heat transfer coefficient with varying degrees of success, and a recent correlation is identified for its superior predictive capability, evidenced by a mean absolute error of 21.7%.
This paper presents a case study of using FLOCAD to model the heat transfer and boiling flow behavior of a Liquid Nitrogen (LN2) ground station equipment heat exchanger. The unique aspect of this work is that is one of very few studies documenting the use of FLOCAD in the spacecraft thermal control community. The results of the simulations described herein were used to guide and plan the execution of thermal vacuum testing spacecraft hardware using a ground test heat exchanger such as the one modeled herein. The paper will review the theory of two-phase boiling flow and correlations used by SINDA/FLUINT as well as the terminology and nomenclature used within Thermal Desktop / FLOCAD regarding two-phase boiling flow heat transfer modeling and simulation. Results for using the ground station heat exchanger as a heat sink in a thermal vacuum test set up are shown in order to demonstrate the set-up and execution of a typical FLOCAD model. Results indicate that properly modeling of the two-phase behavior is critical in order to ascertain the time constant associated with the heat exchanger, as well as understanding the temperature distribution across the test equipment and prediction of required LN2 flowrates to be used during testing. The paper is meant to augment the FLOCAD user’s manual and serve as a tutorial for thermal engineers and analysts wishing to learn and apply FLOCAD to industrial problems.
One of the primary motivations behind this work is the attempt to understand the physics of a two-phase jet pump which constitutes part of a flow boiling test facility at NASA-Marshall. The flow boiling apparatus is intended to provide data necessary to design highly efficient two-phase thermal control systems for aerospace applications. The facility will also be capable of testing alternative refrigerants and evaluate their performance using various heat exchangers with enhanced surfaces. The test facility is also intended for use in evaluating single-phase performance of systems currently using CFC refrigerants. Literature dealing with jet pumps is abundant and covers a very wide array of application areas. Example application areas include vacuum pumps which are used in the food industry, power station work, and the chemical industry; ejector systems which have applications in the aircraft industry as cabin ventilators and for purposes of jet thrust augmentation; jet pumps which are used in the oil industry for oil well pumping; and steam-jet ejector refrigeration, to just name a few. Examples of work relevant to this investigation includes those of Fairuzov and Bredikhin (1995). While past researchers have been able to model the two-phase flow jet pump using the one-dimensional assumption with no shock waves and no phase change, there is no research known to the author apart from that of Anand (1992) who was able to account for condensation shocks. Thus, one of the objectives of this work is to model the dynamics of fluid interaction between a two-phase primary fluid and a subcooled liquid secondary fluid which is being injected employing atomizing spray injectors. The model developed accounts for phase transformations due to expansion, compression, and mixing. It also accounts for shock waves developing in the different parts of the jet pump as well as temperature and pressure dependencies of the fluid properties for both the primary two-phase mixture and the secondary subcooled liquid. The research effort on which this document partly reports described a relatively simple model capable of describing the performance of a two-phase flow jet pump. The model is based on the isentropic homogeneous expansion/compression hypothesis and is capable of fully incorporating the effects of shocks in both the mixing chamber and the throat/diffuser parts of the pump. The physical system chosen is identical to that experimentally tested by Fairuzov and Bredikhin (1995) and should therefore be relatively easy to validate.
The high power density of emerging electronic devices is driving the transition from remote cooling, which relies on con-duction and spreading, to embedded cooling, which extracts dissipated heat on-site. Two-phase microgap coolers employ the forced flow of dielectric fluids undergoing phase change in a heated channel within or between devices. Such coolers must work reliably in all orientations for a variety of applications (e.g., vehicle-based equipment), as well as in microgravity and high-g for other applications (e.g., spacecraft and aircraft). The lack of acceptable models and correlations for orientation- and gravity-independent operation has limited the use of two-phase coolers in such applications. Previous research has revealed that gravita-tional acceleration plays a diminishing role in establishing flow regimes and transport rates as the channel size shrinks, but there is considerable variation among the proposed microscale criteria and limited research on two-phase flows in low aspect ratio mi-crogap channels. Reliable criteria for achieving orientation- and gravity-independent flow boiling would enable emerging sys-tems to exploit this thermal management technique and streamline the technology development process. As a first step toward understanding the effect of gravity on two-phase microgap flow and transport, in the present effort the authors have studied the effect of evaporator orientation and mass flux on near-saturated flow boiling of HFE7100 in a 1.01 mm tall by 13.0 mm wide by 12.7 mm long microgap channel. Orientation-independence, defined as achieving similar critical heat fluxes, heat transfer coefficients, and flow regimes across evaporator orientations, was achieved for mass fluxes of 400 kg/m2-s and greater. The present results are compared to pub-lished criteria for achieving gravity-independence.