SEARCH · Search NASA

Results for “Two-phase systems”

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.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Testing of a Methane Cryogenic Heat Pipe with a Liquid Trap Turn-Off Feature for use on Space Interferometer Mission (SIM)

Camera cooling for SIM presents three thermal control challenges; stable operation at 163K (110 C), decontamination heating to +20 C, and a long span from the cameras to the radiator. A novel cryogenic cooling system based on a methane heat pipe meets these challenges. The SIM thermal team, with the help of heat pipe vendor ATK, designed and tested a complete, low temperature, cooling system. The system accommodates the two SIM cameras with a double-ended conduction bar, a single methane heat pipe, independent turn-off devices, and a flight-like radiator. The turn ]off devices consist of a liquid trap, for removing the methane from the pipe, and an electrical heater to raise the methane temperature above the critical point thus preventing two-phase operation. This is the first time a cryogenic heat pipe has been tested at JPL and is also the first heat pipe to incorporate the turn-off features. Operation at 163K with a methane heat pipe is an important new thermal control capability for the lab. In addition, the two turn-off technologies enhance the "bag of tricks" available to the JPL thermal community. The successful test program brings this heat pipe to a high level of technology readiness.

cyrogenic heat pipe↗

NASA NTRS2014

Temperature Control with Two Parallel Small Loop Heat Pipes for GLM Program

The concept of temperature control of an electronic component using a single Loop Heat Pipe (LHP) is well established for Aerospace applications. Using two LHPs is often desirable for redundancy/reliability reasons or for increasing the overall heat source-sink thermal conductance. This effort elaborates on temperature controlling operation of a thermal system that includes two small ammonia LHPs thermally coupled together at the evaporator end as well as at the condenser end and operating "in parallel". A transient model of the LHP system was developed on the Thermal Desktop (TradeMark) platform to understand some fundamental details of such parallel operation of the two LHPs. Extensive thermal-vacuum testing was conducted with two thermally coupled LHPs operating simultaneously as well as with only one LHP operating at a time. This paper outlines the temperature control procedures for two LHPs operating simultaneously with widely varying sink temperatures. The test data obtained during the thermal-vacuum testing, with both LHPs running simultaneously in comparison with only one LHP operating at a time, are presented with detailed explanations.

0000↗

NASA NTRS2015

Adaptable Single Active Loop Thermal Control System (TCS) for Future Space Missions

This presentation will examine the development of a thermal control system (TCS) for future space missions utilizing a single active cooling loop. The system architecture enables the TCS to be reconfigured during the various mission phases to respond, not only to varying heat load, but to heat rejection temperature as well. The system will consist of an accumulator, pump, cold plates (evaporators), condenser radiator, and compressor, in addition to control, bypass and throttling valves. For cold environments, the heat will be rejected by radiation, during which the compressor will be bypassed, reducing the system to a simple pumped loop that, depending on heat load, can operate in either a single-phase liquid mode or two-phase mode. For warmer environments, the pump will be bypassed, enabling the TCS to operate as a heat pump. This presentation will focus on recent findings concerning two-phase flow regimes, pressure drop, and heat transfer coefficient trends in the cabin and avionics micro-channel heat exchangers when using the heat pump mode. Also discussed will be practical implications of using micro-channel evaporators for the heat pump.

two phase heat transfer↗

NASA NTRS2015

Capillary Channel Flow (CCF) EU2-02 on the International Space Station (ISS): An Experimental Investigation of Passive Bubble Separations in an Open Capillary Channel

It would be signicantly easier to design fluid systems for spacecraft if the fluid phases behaved similarly to those on earth. In this research an open 15:8 degree wedge-sectioned channel is employed to separate bubbles from a two-phase flow in a microgravity environment. The bubbles appear to rise in the channel and coalesce with the free surface in much the same way as would bubbles in a terrestrial environment, only the combined effects of surface tension, wetting, and conduit geometry replace the role of buoyancy. The host liquid is drawn along the channel by a pump and noncondensible gas bubbles are injected into it near the channel vertex at the channel inlet. Control parameters include bubble volume, bubble frequency, liquid volumetric flow rate, and channel length. The asymmetrically confined bubbles are driven in the cross-flow direction by capillary forces until they at least become inscribed within the section or until they come in contact with the free surface, whereupon they usually coalesce and leave the flow. The merging of bubbles enhances, but does not guarantee, the latter. The experiments are performed aboard the International Space Station as a subset of the Capillary Channel Flow experiments. The flight hardware is commanded remotely and continuously from ground stations during the tests and an extensive array of experiments is conducted identifying numerous bubble flow regimes and regime transitions depending on the ratio and magnitude of the gas and liquid volumetric flow rates. The breadth of the publicly available experiments is conveyed herein primarily by narrative and by regime maps, where transitions are approximated by simple expressions immediately useful for the purposes of design and deeper analysis.

fluid dynamics↗

NASA NTRS2016

Lidar-Based Estimates of Above-Ground Biomass in the Continental US and Mexico Using Ground, Airborne, and Satellite Observations

Existing national forest inventory plots, an airborne lidar scanning (ALS) system, and a space profiling lidar system (ICESat-GLAS) are used to generate circa 2005 estimates of total aboveground dry biomass (AGB) in forest strata, by state, in the continental United States (CONUS) and Mexico. The airborne lidar is used to link ground observations of AGB to space lidar measurements. Two sets of models are generated, the first relating ground estimates of AGB to airborne laser scanning (ALS) measurements and the second set relating ALS estimates of AGB (generated using the first model set) to GLAS measurements. GLAS then, is used as a sampling tool within a hybrid estimation framework to generate stratum-, state-, and national-level AGB estimates. A two-phase variance estimator is employed to quantify GLAS sampling variability and, additively, ALS-GLAS model variability in this current, three-phase (ground-ALS-space lidar) study. The model variance component characterizes the variability of the regression coefficients used to predict ALS-based estimates of biomass as a function of GLAS measurements. Three different types of predictive models are considered in CONUS to determine which produced biomass totals closest to ground-based national forest inventory estimates - (1) linear (LIN), (2) linear-no-intercept (LNI), and (3) log-linear. For CONUS at the national level, the GLAS LNI model estimate (23.95 +/- 0.45 Gt AGB), agreed most closely with the US national forest inventory ground estimate, 24.17 +/- 0.06 Gt, i.e., within 1%. The national biomass total based on linear ground-ALS and ALS-GLAS models (25.87 +/- 0.49 Gt) overestimated the national ground-based estimate by 7.5%. The comparable log-linear model result (63.29 +/-1.36 Gt) overestimated ground results by 261%. All three national biomass GLAS estimates, LIN, LNI, and log-linear, are based on 241,718 pulses collected on 230 orbits. The US national forest inventory (ground) estimates are based on 119,414 ground plots. At the US state level, the average absolute value of the deviation of LNI GLAS estimates from the comparable ground estimate of total biomass was 18.8% (range: Oregon,−40.8% to North Dakota, 128.6%). Log-linear models produced gross overestimates in the continental US, i.e., N2.6x, and the use of this model to predict regional biomass using GLAS data in temperate, western hemisphere forests is not appropriate. The best model form, LNI, is used to produce biomass estimates in Mexico. The average biomass density in Mexican forests is 53.10 +/- 0.88 t/ha, and the total biomass for the country, given a total forest area of 688,096 sq km, is 3.65 +/- 0.06 Gt. In Mexico, our GLAS biomass total underestimated a 2005 FAO estimate (4.152 Gt) by 12% and overestimated a 2007/8 radar study's figure (3.06 Gt) by 19%.

Hybrid 3-Phase sampling↗

NASA NTRS2018

Challenges and Opportunities in Gen3 Embedded Cooling with High-Quality Microgap Flow

Gen3, Embedded Cooling, promises to revolutionize thermal management of advanced microelectronic systems by eliminating the sequential conductive and interfacial thermal resistances which dominate the present 'remote cooling' paradigm. Single-phase interchip microfluidic flow with high thermal conductivity chips and substrates has been used successfully to cool single transistors dissipating more than 40kW/sq cm, but efficient heat removal from transistor arrays, larger chips, and chip stacks operating at these prodigious heat fluxes would require the use of high vapor fraction (quality), two-phase cooling in intra- and inter-chip microgap channels. The motivation, as well as the challenges and opportunities associated with evaporative embedded cooling in realistic form factors, is the focus of this paper. The paper will begin with a brief review of the history of thermal packaging, reflecting the 70-year 'inward migration' of cooling technology from the computer-room, to the rack, and then to the single chip and multichip module with 'remote' or attached air- and liquid-cooled coldplates. Discussion of the limitations of this approach and recent results from single-phase embedded cooling will follow. This will set the stage for discussion of the development challenges associated with application of this Gen3 thermal management paradigm to commercial semiconductor hardware, including dealing with the effects of channel length, orientation, and manifold-driven centrifugal acceleration on the governing behavior.

microchannel↗

NASA NTRS2019

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen↗

NASA NTRSDate not supplied

Modeling Electrolytic O2 Recovery from Metabolic CO2 for Advanced Closed Loop Life Support Systems in Extraterrestrial Human Missions

The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration human missions to the Moon and Mars will necessitate a sustainable and highly efficient metabolic oxygen recovery system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current metabolic O2 recovery efficiency, expand mission sustainability, and reduce complexity of the system. The novel design combines CO2 conversion to O2 along with C2H4 as byproduct and water electrolysis (currently conducted in two separate units) into a single compact unit that runs at standard conditions and is theoretically capable of generating O2 with a theoretical maximum metabolic CO2 conversion of 73% while consuming less than metabolic water. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and electrical current density distributions is coupled in the model with all the other physics phenomena involved in the MFECR’s process, such as two-phase flow, free and porous fluid regimes, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The EC reaction sections of the MFECR consists of two porous gas diffusion electrodes (GDE) and an electrolyte serpentine channel sandwiched in the middle. The CO2 feeds the cathode serpentine channel and part of the O2 product is fed back to the anode serpentine chamber. An alkaline solution feeds the electrolyte serpentine chamber wetting the GDEs of both, the anode and cathode allowing the OH- ionic transport between them. The EC reactions in the cathode’s GDE yield C2H4 from CO2 and H2 from water while the EC reaction in the anode’s GDE yields O2. The authors will present in this paper the validation of the model using experimental data and the utilization of the validated model in building a reliable simulator that will not only assist the authors on the MFECR design but also the optimization of its operation in the ISS and future spatial human missions.

Jesus A Dominguez↗

NASA NTRS2021

Analysis of Arc Welding Process in Space

This work is motivated by NASA plans to conduct welding experiments on the ISS. It is expected that deployment of welding and additive manufacturing technologies in the space environment has the potential to revolutionise how orbiting platforms are designed, manufactured, and assembled. However, the structure, composition and quality of a weld is extremely dependent on the environment and can be difficult to control in space. Shielding gases would also be tough to manage as gases behave differently in zero gravity and airless environments. Additional points of concern are related to the spatter and sparks dynamics in space. Therefore, there is a need for a more basic understanding of welding processes by computational modelling. To provide such an insight we developed state of art models of the ARC torch, droplet detachment and transfer, and the melt pool build up using magnetohydrodynamics approximation and level set method for two-phase liquid metal/-gas flow modelling. Two finite element models were built using 2D axisymmetric geometry in COMSOL Multiphysics®: (i) stationary model of the ARC torch and (ii) dynamical model of droplet detachment and transfer. Both models demonstrate reasonable agreement with earlier experimental observations and high sensitivity to the temperature dependence of the thermophysical parameters on the system materials. The models were used to provide physical insight into ARC welding in various environments and geometries. It was shown, in particular, that there is a significant probability of gas bubbles trapping in the meltpool and a possibility of unbounded wondering of sparks in the welding chamber in zero gravity. In addition, we estimated metal evaporation rate that may be hazardous in the confined environment of the ISS. These issues raise concern of quality of the weld and safety for ISS applications.

ARC welding↗

NASA NTRS2023

Universal Two-Phase Convection Heat Transfer Correlations for Cryogenic Pipe Chilldown

This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.

Universal Correlations↗

NASA NTRSDate not supplied

The ZBOT-NC Experiment - Effects of Non-Condensable Gases on Propellant Tank Pressurization and Pressure Control

Introduction Integral to all phases of NASA’s projected planetary expeditions is affordable and reliable cryogenic fluid storage for use in propellant or life support systems. It is greatly advantageous to develop innovative vent-less pressure control designs based on cooling/mixing of the bulk tank fluid to allow storage of the cryogenic fluid with zero or reduced boil-off. The presence of noncondensable gases, can interfere with the condensation at the interface impacting tank pressure control during subcooled jet mixing especially in microgravity. The Zero-Boil-Off Tank (ZBOT) Experiments are a series of small-scale experiments aboard the International Space Station (ISS) that use a transparent volatile simulant fluid in a transparent sealed tank to delineate various fundamental fluid flow, heat and mass transport, and phase change phenomena associated with storage tank pressurization and pressure control in microgravity. The ZBOT-1 experiment was performed on the ISS 2017-2018 timeframe and collected data to validate a state-of-the-art CFD model for tank pressurization and pressure control for a pure system. The ZBOT-NC Experiment is the second experiment in the series to be performed on the International Space Station (ISS) in 2025. Its goal is to investigate the effects of noncondensable gases on interfacial evaporation and condensation during self -pressurization and jet-mixing pressure control in microgravity for a two-component system. Materials & Methods In this work, we will describe the detailed features of the ZBOT-NC experimental hardware and Diagnostics that includes nonintrusive Quantum Dot Thermometry (QDT) for whole field temperature measurement. All microgravity pressurization and pressure control tests will be performed for both the pure and the two-components systems with Xenon and Neon as the two noncondensable gases spanning small and large molecular weights and sizes. The two-phase CFD model that is developed as part of the project will be also presented and discussed. Results Ground-based pressurization and jet mixing experiments and CFD simulation results are compared to each other to validate both the fidelity of the CFD model predictions and the accuracy of the QDT measurement. Model simulations for noncondensable gas effects will also be compared against large Cryogenic LH2-GHe experiments to indicate the noncondensable gas effects on tank pressure control during jet mixing in 1G. Finally, CFD results will be presented to predict the effects of the noncondensable gas during subcooled jet mixing during the ZBOT-NC in advance of the microgravity experiment in 2015.

Evaporation Condensation↗

NASA NTRSDate not supplied

The Zero-Boil-Off Tank Experiment: the Effects of Non-Condensable Gases on Pressurization and Pressure Control of Propellant Tanks in Microgravity

Evaporation Condensation↗

DOE OSTI2017

SAS4A/SASSYS-1 VERSION 5.2 SAFETY ANALYSIS CODE SYSTEM

SAS4A/SASSYS-1 is a software simulation tool used to perform deterministic safety analysis of anticipated events as well as design basis and beyond design basis accidents for advanced nuclear reactors. This software can be used to assess the safety of a prescribed reactor design, but it cannot be used to configure a design to meet targeted performance objectives. Detailed, mechanistic models of steady-state and transient thermal, hydraulic, kinetic, and mechanical phenomena are employed to describe the response of the reactor core, the reactor primary and secondary coolant loops, the reactor control and protection systems and the balance-of-plant to accidents caused by changes in coolant flow, oss of heat rejection, or reactivity insertion. The consequences of single and double-fault accidents can be modeled, including fuel and coolant heating, fuel and cladding mechanical behavior, core reactivity feedbacks, coolant loop performance including natural circulation, and decay heat removal. Analyses are typically terminated upon demonstration of reactor and plant shutdown to permanently coolable conditions or upon violation of design basis margins. The objective of the analysis is to quantify accident consequences as measured by the transient behavior of system performance parameters, such as fuel and cladding temperatures, reactivity, and cladding strain. Originally developed for analysis of sodium cooled reactors with oxide fuel clad by stainless steel, the models In SAS4A/SASSYS-1 were subsequently extended and specialized to metallic fuel clad with advanced alloys and to several other coolant options including lead and LBE. METHOD OF SOLUTION: Each subassembly in the reactor core is represented with single or multiple-pin models. One channel represents one or more fuel pins and many channels are employed for a whole-core representation. Heat transfer in each pin is modeled with a two-dimensional (r/z) heat conduction equation. Single and two-phase coolant thermal-hydraulics are simulated with a unique, one-dimensional (axial) multiple-bubble liquid metal boiling model. The transient fuel and cladding mechanical behavior models are integrated with fission product production,release, and transport models. UNIQUE FEATURES: The physical models in SAS4AJSASSYS-1 are highly detailed 1umerical representations of reactor accident conditions based on extensive laboratory and test reactor results. The models are specialized to liquid metal (sodium) cooled fast reactors with oxide or metallic fuel clad with stainless steel.«

Fanning, Thomas↗

DOE OSTI2016

Evaluating performance of two-group interfacial area transport equation for vertical small and large diameter pipes

In the two-fluid transport model, the coupling of mass, momentum and energy transfer between phases is highly dependent on interfacial area transfer terms. Several research efforts in the past have been focused on the development of an interfacial area transport equation model (IATE), in order to eliminate the drawbacks of static flow regime maps currently used in best-estimate thermal-hydraulic system codes. The IATE attempts to model the dynamics that are involved in two phase flows by accounting for the different interaction mechanisms affecting bubble transport in the flow. The further development and validation of IATE models has been hindered by the lack of adequate experimental data in regions beyond the bubbly flow regime. At the Helmholtz-Zentrum Dresden-Rossendorf (HZDR) experiments utilizing wire mesh sensors have been performed over all flow regimes, establishing a database of high-resolution (in space and time) data. Two sets of data from small and large diameter pipes fitted with wire mesh sensors is utilized in this work. Analysis of flow conditions in the bubbly, churn turbulent and annular flow regimes is presented. The performance of the Fu-Ishii two-group model is evaluated against small diameter database. Results indicate good performance (< 10% error) for small group 1 bubbles, and poor performance for large group 2 bubbles. The Smith two-group large diameter IATE model is evaluated for the large diameter database. In low void-fraction regimes, the Smith model performs well. In high void-fraction regimes, there is poor group-wise interfacial area prediction - however the total interfacial area is erroneously predicted well as group-wise errors compensate each other. Overall, the study suggests that further efforts and re-evaluation of closure terms are needed in order to extend the range of validity of the IATE models. (authors)

42 ENGINEERING↗

DOE OSTI2016

Improving performance of two-group interfacial area transport equation models using genetic algorithms and high resolution experimental data

In the two-fluid transport model, the coupling of mass, momentum and energy transfer between phases is highly dependent on interfacial area transfer terms. Several research efforts in the past have been focused on the development of an interfacial area transport equation model (IATE), in order to eliminate the drawbacks of static flow regime maps currently used in best-estimate thermal-hydraulic system codes. The IATE attempts to model the dynamics that are involved in two phase flows by accounting for the different interaction mechanisms affecting bubble transport in the flow. In general, current IATE models perform poorly in high void fraction flow regimes. As both models depend on closure relations that require empirically determined coefficients, a genetic algorithm is implemented to study their optimization. The preliminary study indicates that improvements are achievable by application of the genetic algorithm. Optimization of coefficients for an entire set of high-resolution wire mesh sensor (WMS) data provided a single set of optimal coefficients. The application of WMS-optimized coefficients to external tests improved IATE performance for similar hydraulic diameters and provided a neutral change in performance for significantly smaller/larger hydraulic diameters. This is a crucial outcome as an optimized set of coefficients that provide significant improvement for a subset of experimental data - is shown to benefit performance in external data for which the coefficients are not specifically optimized. (authors)

42 ENGINEERING↗

DOE OSTI2016

Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation

Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

Experimental Study for Two Phase Flow Instability in the PCHE Type Steam Generator in SMR

Conventional U-tube type heat exchangers have some limitations when applied to Small Modular Reactor (SMRs). In particular, the size of steam generators is important in the SMR, as it greatly affects the size, and thus the cost, of the reactor vessel.. Therefore, the circuit heat exchanger (PCHE) is suggested as an alternative. PCHE, a compact type heat exchanger manufactured by HEATRIC, has several advantages such as the reduced size, wide operation range until 900? and 900 bar and high power density due to the large heat transfer area over volume ratio. Recently, it is newly suggested as an application of a steam generator. MIT and KAIST have already conceptually suggested PCHE as a steam generator. Under boiling condition, flow instability is one of the critical issues-resulting in unstable mass flow rate, sudden temperature change and system control failure. According to the flow mechanism, there are various types of instabilities in the heat exchanger and each flow tendency is different. Typical examples are Ledinegg instability of the static instability and pressure-drop oscillation and density-wave oscillation of the dynamic instability. The characteristic of Ledinegg instability is that the flow under unstable condition changes into the stable flow. The characteristic of dynamic instability is that the flow fluctuates continuously. Pressure-drop oscillation has the period of tens of seconds and the system pressure also highly oscillates with the period of flow instability. In comparison, the density-wave oscillation has the period of few seconds and the system pressure is relatively stable. Previous studies on flow instability are focused on the single channel flow and helical tube type heat exchanger. There has been very limited experimental research on flow instability using PCHE. In KAIST, the test facility was constructed to study the thermal hydraulics in PCHE. The objectives of our test facility are: (1) to identify the occurrence of instability under two phase and super-heated steam flow in PCHE, and (2) to examine whether it can operate normally using throttling effect. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2016

Pressure Drop-Flow Rate Curves for Single Phase Steam in Steam Generator U-tubes During Severe Accidents

Thermally induced steam generator tube rupture (TI-SGTR) can occur in pressurized water reactors (PWRs) with U-tube steam generators (SGs) during high pressure severe accident sequences resulting in containment bypass. Robust natural circulation patterns including in-vessel natural circulation and countercurrent flow to and from the SGs of superheated steam and other hot gases heats cooler structures in the reactor coolant system (RCS). The buoyancy driven flows return cooler and denser gases back to the core slowing the progression of some core degradation phenomena and failure of the lower head. As much as 50 percent of the decay energy and energy produced from core degradation oxidation reactions can be transported from the core to SGs resulting in the heat up and thermal creep of RCS pressure boundary components including hot leg nozzles, pressurizer surge line, and the SG tubes, many of which contain flaws. The creep failure of an individual tube is dependent on the transitory heating of the tube, a function of the mass flow rate through the tube and local temperature history of steam entering the tube. Buoyancy driven natural circulation flow rates are determined by changes in gravitational head, a combination of changes in fluid density due to heat transfer and change in elevation, and irreversible pressure drops resulting primarily from friction loss in the flow paths and other form losses. The nonuniform flow distribution in U-tubes of varying lengths and heights during single-phase liquid and two-phase natural circulation has been extensively studied in scaled test facilities and verified numerically in previous works. We hypothesize the countercurrent natural circulation of single-phase superheated steam in U-tube SGs is governed by the same characteristic pressure drop-flow rates curves of single-phase liquid and two-phase natural circulations in U-tubes connected to common plena derived in. Secondly, we propose flow distribution is a function of U-tube geometry - height, length, hydraulic diameter, etc., and the larger SGs of Combustion Engineering (CE) style 2 x 4 PWRs, which are the focus of this study, will show greater nonuniformity in tube by tube flow rates and ratio of forward flowing tubes to reverse flowing tubes. This work derives characteristic pressure drop-flow rate curves for single-phase superheated steam in CE-type SGs during a high pressure station blackout (SBO) sequence with failed auxiliary feedwater (AFW) and dry secondary side which is closely related to the TI-SGTR event. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗