Search NASA⌕ Search

SEARCH · Search NASA

Results for “two phase heat transfer”

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 55 records · Page 3

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗

In tube condensation of low global warming potential refrigerants in an axial micro-fin aluminum tube

Environmental regulations have driven the development of refrigerants with low global warming potential (GWP). To design heat exchangers using these new refrigerants, data are needed concerning the heat transfer coefficient and pressure drop in two-phase flow. Another change is the increasing use of aluminum tubes rather than copper tubes to reduce heat exchanger cost. Hence, this study presents an experimental investigation of flow condensation using an expanded axial micro-fin aluminum tube with a fin-tip diameter of 5.96 mm. Here, the experiments included single compounds R-32, R-1234yf, and R-1234ze(E), zeotropic mixtures with low glide (R-454B), and zeotropic mixtures with high-glide (R-454C and R-455A). Experiments were conducted at condensation temperatures ranging from 40 °C to 50 °C, reduced pressures ranging from 0.21 to 0.55, and mass fluxes ranging from 150 to 350 kg/(m 2 s). Data obtained for these refrigerants constitute one of the first reports for high-glide refrigerants using axial micro-fin aluminum tubes. An evaluation of heat transfer degradation of zeotropic mixtures due to mass transfer resistance at the liquid/vapor interface is presented. This information can be used to design heat exchangers for next generation air conditioning and refrigeration systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Experimentation on Finned-Tube Microchannel Heat Exchanger Incorporating Phase Change Material and R-410A

Integrating latent heat thermal energy storage, particularly phase change materials (PCM), in building heating, ventilation, air conditioning, and refrigeration (HVAC&R) systems, offers substantial benefits for peak load shifting and energy efficiency. While numerous studies explore PCM-embedded heat exchangers (PCM-HX) using single-phase heat transfer fluids, there needs to be more investigation regarding PCM-HX utilizing two-phase refrigerant flow. This gap is even more pronounced for microchannel heat exchanger (MCHX) configurations. In this study, an MCHX-based PCM-HX was built and tested to analyze and evaluate the performance of a PCM-HX operating with a two-phase refrigerant as the working fluid. A commercially available MCHX condenser unit was embedded with 2.6 kg of RT35 PCM and tested using R-410A as the working fluid, using an in-house test loop to control the refrigerant inlet conditions precisely. Thin and minimally invasive T-type thermocouples traced the temperature evolution along MCHX channels and tube banks. Additional larger T-type thermocouple probes at the refrigerant inlet and outlet measured the average power of the PCM-HX during condensing tests, ranging from 0.7 kW to 1.6 kW based on mass flow rates varying from 3.5 g/s to 8 g/s. This experimental study on MCHX-based PCM-HX with a two-phase refrigerant has the potential to offer practical design guidelines for the direct integration of PCM-TES in HVAC&R systems without the use of a secondary loop, while caution is advised to avoid a substantial increase in the refrigerant charge amount.

42 ENGINEERING↗

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗

Benchmarking the design of the cryogenics system for the underground argon in DarkSide-20k

DarkSide-20k (DS-20k) is a dark matter detection experimentunder construction at the Laboratori Nazionali del Gran Sasso (LNGS)in Italy. It utilises ∼ 100 t of low radioactivity argon from anunderground source (UAr) in its inner detector, with half serving astarget in a dual-phase time projection chamber (TPC). The UArcryogenics system must maintain stable thermodynamic conditionsthroughout the experiment's lifetime of over10 years. Continuous removal ofimpurities and radon from the UAr is essential for maximising signalyield and mitigating background. We are developing an efficient andpowerful cryogenics system with a gas purification loop with atarget circulation rate of1000 slpm. Central to itsdesign is a condenser operated with liquid nitrogen which is pairedwith a gas heat exchanger cascade, delivering a combined coolingpower of more than 8 kW. Here wepresent the design choices in view of the DS-20k requirements, inparticular the condenser's working principle and the coolingcontrol, and we show test results obtained with a dedicatedbenchmarking platform at CERN and LNGS. We find that the thermalefficiency of the recirculation loop, defined in terms of nitrogenconsumption per argon flow rate, is95 % and the pressure in the testcryostat can be maintained within±(0.1–0.2) mbar. Wefurther detail a 5-day cool-down procedure of the test cryostat,maintaining a cooling rate typically within-2 K/h, as required for theDS-20k inner detector. Additionally, we assess the circuit's flowresistance, and the heat transfer capabilities of two heat exchangergeometries for argon phase change, used to provide gas forrecirculation. We conclude by discussing how our findings influencethe finalisation of the system design, including necessarymodifications to meet requirements and ongoing testing activities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Pulsed laminar falling films in vertical tubes: Maintaining a continuous liquid film with reduced film thickness

Falling-liquid-film heat transfer is important to a wide variety of industrial applications, including two-phase thermosyphons, falling-film evaporators in the food industry, heat exchangers in fuel-ethanol processing, and the chemical service industry. Particularly for laminar falling films, the heat transfer increases as the film thickness decreases, hence it is desirable to keep the films as thin as possible. This article presents a periodically-pulsed flow concept to reduce the average mass flow inside a vertical tube. The technique provides thinner films than can be achieved with continuous flow. A drainage model is adopted to predict the film thickness as a function of flow rate and location along the tube. The concept is confirmed experimentally by measuring the film thickness in real time in a circular tube using a laser-based absorption technique. A new heat-transfer-based film continuity measurement is also integrated into the optical film measurement setup. Water at ambient temperature and pressure is used for this work. Film thickness reductions of up to 60%, and liquid-film volume-flow reductions of up to 80% compared to the minimum continuous-flow film have been obtained. The time-varying film thickness profile can also be controlled to a limited extent, which may be of interest for process industries.

42 ENGINEERING↗

Multi-fidelity electrochemical modeling of thermally activated battery cells

Thermally activated batteries undergo a series of coupled physical changes during activation that influence battery performance. These processes include energetic material burning, heat transfer, electrolyte phase change, capillary-driven two-phase porous flow, ion transport, electrochemical reactions, and electrical transport. Several of these processes are strongly coupled and have a significant effect on battery performance, but others have minimal impact or may be suitably represented by reduced-order models. Additionally, assessing the relative importance of these phenomena must be based on comparisons to a high-fidelity model including all known processes. In this work, we first present and demonstrate a high-fidelity, multi-physics model of electrochemical performance. This novel multi-physics model enables predictions of how competing physical processes affect battery performance and provides unique insights into the difficult-to-measure processes that happen during battery activation. We introduce four categories of model fidelity that include different physical simplifications, assumptions, and reduced-order models to decouple or remove costly elements of the simulation. Using this approach, we show an order-of-magnitude reduction in computational cost while preserving all design-relevant quantities of interest within 5 percent. The validity of this approach and these model reductions is demonstrated by comparison between results from the full fidelity model and the different reduced models.

25 ENERGY STORAGE↗

COMPUTATION FLUID DYNAMICS ANALYSIS FOR GENERIC SMALL MODULAR REACTOR CONTAINMENT SEPARATE EFFECTS TEST

It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

C-FOG: Life of Coastal Fog

C-FOG is a comprehensive bi-national project dealing with the formation, persistence, and dissipation (life cycle) of fog in coastal areas (coastal fog) controlled by land, marine, and atmospheric processes. Given its inherent complexity, coastal-fog literature has mainly focused on case studies, and there is a continuing need for research that integrates across processes (e.g., air–sea–land interactions, environmental flow, aerosol transport, and chemistry), dynamics (two-phase flow and turbulence), microphysics (nucleation, droplet characterization), and thermodynamics (heat transfer and phase changes) through field observations and modeling. Central to C-FOG was a field campaign in eastern Canada from 1 September to 8 October 2018, covering four land sites in Newfoundland and Nova Scotia and an adjacent coastal strip transected by the Research Vessel Hugh R. Sharp . An array of in situ, path-integrating, and remote sensing instruments gathered data across a swath of space–time scales relevant to fog life cycle. Satellite and reanalysis products, routine meteorological observations, numerical weather prediction model (WRF and COAMPS) outputs, large-eddy simulations, and phenomenological modeling underpin the interpretation of field observations in a multiscale and multiplatform framework that helps identify and remedy numerical model deficiencies. An overview of the C-FOG field campaign and some preliminary analysis/findings are presented in this paper.

58 GEOSCIENCES↗

CTF Improved Drag Model and Flow Regime Transition Criteria

The demand for accurate prediction of two-phase flow behavior in a boiling water reactor (BWR) requires a comprehensive understanding of flow regime, void fraction, heat transfer, and pressure drop. The CTF subchannel code, which is used for the Thermal/Hydraulic (T/H) solution in the Consortium for Advanced Simulation of Light Water Reactors (CASL)-developed Virtual Environment for Reactor Application (VERA) core simulator, is being further developed for BWR applications. In support of this goal, the present work highlights some of the two-phase closure model developments towards improving the CTF void fraction prediction, especially for subcooled boiling. The drift-flux approach has been well-developed for upward dispersed two-phase flows and proven to be accurate in predicting void fraction in bubbly and slug flow regimes. In this work, these kinematic constitutive relations for the drift-flux velocity have been implemented into CTF to describe the interfacial drag of bubbly flow as an alternative to the existing model for better void fraction prediction. The success of these constitutive relations also relies on a good flow regime map that accounts for flow conditions and channel geometry. A more reliable flow regime transition criteria that account for the flow condition has also been implemented in this study for modeling the flow regime transition criteria. The newly implemented models are shown to give improved void fraction predictions in comparison to experimental data.

Hizoum, Belgacem↗

Why ionizing radiation enhances surface wettability

Radiation-Induced Surface Activation is an inherent phenomenon where surfaces that are exposed to gamma irradiation are observed to undergo an increase in wettability. This increase in wettability as a result of the ionizing radiation exposure has so far been demonstrated to have a pronounced impact on Leidenfrost temperature and two-phase fluid dynamics. Test results from previous experiments have shown that incorporation of this effect on heat transfer equipment design may increase the thermal-hydraulic margin leading to higher thermal efficiency. However, the mechanism behind the increased wettability is not clearly understood. In the present work, three different materials (Zircaloy-4, 316 stainless steel, and copper) were exposed at two different dose rates with use of two different gamma irradiation facilities. A detailed surface characterization on the post-irradiated samples is carried out to understand the changes in surface chemistry, wettability and surface morphology. It is observed from the experiments that the increase in wettability upon irradiation depended on the total dose and not on the dose rate. Moreover, localized oxidation and porosity induced by radiolysis was seen to be the predominant mechanism behind increased wettability which leads to improved Leidenfrost temperature.

36 MATERIALS SCIENCE↗

Why ionizing radiation enhances surface wettability

Radiation-Induced Surface Activation is an inherent phenomenon where surfaces that are exposed to gamma irradiation are observed to undergo an increase in wettability. This increase in wettability as a result of the ionizing radiation exposure has so far been demonstrated to have a pronounced impact on Leidenfrost temperature and two-phase fluid dynamics. Test results from previous experiments have shown that incorporation of this effect on heat transfer equipment design may increase the thermal-hydraulic margin leading to higher thermal efficiency. However, the mechanism behind the increased wettability is not clearly understood. In the present work, three different materials (Zircaloy-4, 316 stainless steel, and copper) were exposed at two different dose rates with use of two different gamma irradiation facilities. A detailed surface characterization on the post-irradiated samples is carried out to understand the changes in surface chemistry, wettability and surface morphology. It is observed from the experiments that the increase in wettability upon irradiation depended on the total dose and not on the dose rate. Moreover, localized oxidation and porosity induced by radiolysis was seen to be the predominant mechanism behind increased wettability which leads to improved Leidenfrost temperature.

Seshadria, Arunkumar↗

Extension of Clad Damage Propagation Model for Fission Gas Dispersal and Two-Phase Flow Effects in MOOSE SubChannel Module

This report presents an extension of the Clad Damage Propagation (CDAP) model implemented in the MOOSE SubChannel Module (SCM) to capture post-failure fission-gas dispersal and two-phase flow effects in sodium-cooled fast reactor assemblies. The extended model tracks discharged gas axially and radially, computes channel-averaged flow quality and void fraction using a Lockhart–Martinelli framework, evaluates two-phase frictional pressure-drop multipliers, determines inlet mass-flow degradation under fixed core pressures, and applies an intensified-void-based heat-transfer degradation to affected fuel pins. Radial plume expansion is parameterized using mineral-oil jet experiments mapped to sodium conditions via Reynolds–Weber similarity. Implementation details are documented, along with the new methods and user inputs needed to control plume mapping and two-phase behavior. Demonstration simulations for 19- and 37-pin bundles show that breach size and inlet velocity strongly influence propagation potential: small breaches (≤0.5 mm) produce limited degradation while larger breaches (~1 mm) can drive oscillatory temperature spikes and enhanced failure propagation, especially at higher velocities. These results demonstrate that the extended CDAP model provides a more complete framework for quantifying cladding damage propagation and evaluating propagation potential in transient scenarios. The approach remains computationally efficient, consistent with subchannel-level analysis, yet incorporates sufficient physics to bridge localized post-failure effects with bundle- and assembly-scale degradation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗

Chapter 14: Heat and Mass Transfer

Fluidized beds are used in a variety of processes due to their excellent heat and mass transfer properties. Heat transfer in fluidized beds can refer to interphase heat transfer, or more commonly to bed-to-surface heat transfer. Heat transfer between bubble phase and dense phase can be characterized by an interphase heat transfer coefficient. This chapter elaborates the derivation of the interphase heat transfer coefficient by an analogy to mass transfer. Mass transfer in fluidized beds usually refers to gas-to-particle mass transfer, or interphase mass transfer, where a component migrates from particles to fluids, or vice versa. The chapter describes the two major components of the interphase mass transfer: particle and fluid mass transfer in the dense phase; and bubble- to dense-phase mass transfer.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗