Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1986

Infrared heating of Comet Halleys atmosphere

The warm circumnuclear dust in the inner cometary coma reradiates in the IR in the wavelength range of the ground state rotational band of the dominant atmospheric molecule H2O. However, the interaction of this radiation with H2O has hitherto not been taken into account in cometary atmospheric models. Here, an earlier two-phase, multifluid model of the dusty atmosphere is extended by including this effect. Although this IR radiation initially pumps the rotational levels of H2O, frequent intermolecular collisions in the inner coma transfer this energy from rotational modes to translational modes. As a result the temperature in the innermost coma no longer decreases to about 10 K, as predicted by the earlier models, but reaches a minimum of only about 120 K.

Marconi, M. L.↗
NASA NTRS1985

Ignition and Flame Spread Above Liquid Fuel Pools

Phenomena of ignition and flame spread above liquid fuel pools are studied and factors that can improve fire safety are determined. Dominant mechanisms for convective heat transfer and impact on ignition and flame spread above liquid fuel pools will be determined from these efforts. The approach is: (1) analytical and computational studies to evaluate scaling factors; (2) experimental design, development, and operation in laboratory at Earth's gravity; (3) experimental design for drop towers at NASA Lewis Research Center; (4) development and operation of drop tower experiments together with NASA Lewis Research personnel and; (5) develop recommendations for Lear Jet and/or Space Shuttle experiments. Through this research, a mathematical model of two-phase systems is formulated and is being coded for the computer.

Sirignano, W. A.↗
NASA NTRS1978

Heat Transfer and Fluid Mechanics Institute, Meeting, 26th, Washington State University, Pullman, Wash., June 26-28, 1978, Proceedings

Papers are presented on the heat transfer and fluid mechanics in energy systems including diesel combustion analysis using rapid sampling and the numerical simulation of thermal response of architectural structures. Experiments associated with hypersonic boundary layers, the prediction of nozzle-wall pressure distribution, and the effects of nozzle geometry on an axisymmetric thrust augmentor are outlined. Predictions are made of film thickness on a roll-metering device and the thermal characteristics of multiple-surface submerged condenser modules are described. Numerical studies are made of the interaction between natural and forced convection in a complex geometry and of pressure drop and holdup in two-phase stratified flow. A numerical solution is given for viscous flows at high Reynolds numbers and for compressible boundary layer flows near separation or reattachment. The structure and transport properties in turbulent shear flows are reviewed along with the spectra of Reynolds stresses in a developing pipe flow.

Crowe, C. T.↗
NASA NTRS1981

A pumped, two-phase flow heat transport system for orbiting instrument payloads

A pumped two-phase (heat absorption/heat rejection) thermal transport system for orbiting instrument payloads is investigated. The thermofluid characteristics necessary for the system design are discussed. A preliminary design with a series arrangement of four instrument heat stations and six radiators in a single loop is described in detail, and the total mass is estimated to be 134 kg, with the radiators, instrument heat stations, and fluid reservoir accounting for approximately 86, 24, and 12 kg, respectively. The evaluation of preliminary test results shows that the system has potential advantages; however, further research is necessary in the areas of one-g and zero-g heat transfer coefficients/fluid regimes, fluid by-pass temperature control, and reliability of small pumps.

Fowle, A. A.↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
NASA NTRS1986

Evaluation of active thermal control options for Space Station

An analysis of various Space Station (MSS) active thermal control systems options under consideration because of their potential low weight, efficiency and reliability is reported. The study addressed ordinary and diode-action body mounted radiators, thermal storage, the area and pumping power requirements for single-phase cooling of cold plates, and single-phase and two-phase active cooling loops. The base configuration considered was a core MSS formed by four habitable modules on which are mounted heat pipe radiators articulated to be always edge-on to the sun. A simulation was performed which accounted for the available heat sinks, several thermal loads and the heat rejection capability. No benefits were found with diode-action radiators if the solar absorption is 0.1 or less, although diode-action heat pipes will maintain a higher level of performance in the presence of coating degradation. Thermal storage becomes important only with radiator coating degradation. Water can be up to three times as efficient as Freon as a heat transfer medium. Finally, single-phase cooling offers a lower system mass than two-phase cooling if varying temperature heat loads can be accommodated.

Schuster, J. R.↗
NASA NTRS1985

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗
DOE OSTI2026

Holistic energy analysis method for thermal management architectures of data centers

Modern high-performance computing (HPC) data centers (DCs), particularly those supporting energy-intensive artificial intelligence (AI) workloads, face escalating thermal management challenges that degrade performance through thermal throttling and drive up cooling power consumption and operational costs. To address this challenge, many have developed a wide variety of thermal management solutions (single-phase, two-phase, direct, indirect, hybrid, and more) which attempt to cool HPC DCs effectively while attempting to minimize overall system power consumption. However, the analysis of these solutions and methods to effectively compare one with another is lacking. Overall power usage effectiveness (PUE) and total-power usage effectiveness (TUE) provide a metric to quantify power consumption but fail to identify components in the system which require further optimization. To address this, we propose a holistic analytical framework – the waterfall diagram (WFD) – which leverages a waterfall chart methodology, offering a comprehensive visualization of both the thermal management system loop and heat flow pathways from individual server components to the outdoor ambient. Use of the WFD enables graphical estimations of power efficiency and cooling performance across each component of a DC cooling system and complements Sankey-style energy flow visualizations by additionally resolving stage-wise temperature changes and incremental TUE contributions. The framework is used in conjunction with simulation-based approaches, to conduct a detailed pressure drop and flow distribution analysis aimed at identifying the optimal coolant distribution architecture for a single-phase direct-to-chip water-cooled DC, which serves as the baseline for subsequent WFD analysis. Among the evaluated architectures, the 3 U modular coolant distribution architecture is found to demonstrate the best performance, considering minimal pressure drop and uniform flow distribution. In addition, TUE is calculated for each cooling loop component based on its associated pressure drop and corresponding pumping power, which are integrated into the WFD. This correlation between TUE and local temperature offers immediate insight into the power efficiency and thermal performance contributions of individual components, facilitating further development and optimization. Examples of WFD applications are presented under varying thermal loads and ambient conditions, demonstrating reasonable cooling strategies. Notably, the 3 U modular architecture maintains a consistent chip case temperature of 85°C, achieving a TUE of 1.016 at ambient temperature of 47°C, and a TUE of 1.026 at ambient temperature of 52°C. The WFD methodology provides an efficient, holistic, and streamlined framework for DC thermal management architecture assessment and enables design optimization which is important for addressing the thermal-fluidic energy challenges of current and next-generation DCs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for heat pipes

Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in heat pipes, including capillary, phase change, turbulence, and compressibility effects, there are high uncertainties in the predictability of their operational regimes and performance. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. Additional analyses and discussion are provided for phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics. The discussions included the recognizing challenges and proposing future research directions for both modeling and simulation and experimental efforts. Additionally, the report addresses phenomena with medium importance and low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geyser boiling, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. In conclusion, this comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1977

Assessment of analytical techniques for predicting solid propellant exhaust plumes

The calculation of solid propellant exhaust plume flow fields is addressed. Two major areas covered are: (1) the applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size and particle size distributions, and (2) thermochemical modeling of the gaseous phase of the flow field. Comparisons of experimentally measured and analytically predicted data are made. The experimental data were obtained for subscale solid propellant motors with aluminum loadings of 2, 10 and 15%. Analytical predictions were made using a fully coupled two-phase numerical solution. Data comparisons will be presented for radial distributions at plume axial stations of 5, 12, 16 and 20 diameters.

Tevepaugh, J. A.↗
NASA NTRS1977

Conceptual design for spacelab two-phase flow experiments

KC-135 aircraft tests confirmed the gravity sensitivity of two phase flow correlations. The prime component of the apparatus is a 1.5 cm dia by 90 cm fused quartz tube test section selected for visual observation. The water-cabin air system with water recycle was a clear choice for a flow regime-pressure drop test since it was used satisfactorily on KC-135 tests. Freon-11 with either overboard dump or with liquid-recycle will be used for the heat transfer test. The two experiments use common hardware. The experimental plan covers 120 data points in six hours with mass velocities from 10 to 640 kg/sec-sq m and qualities 0.01 to 0.64. The apparatus with pump, separator, storage tank and controls is mounted in a double spacelab rack. Supporting hardware, procedures, measured variables and program costs are defined.

Bradshaw, R. D.↗
NASA NTRS1980

Two-phase working fluids for the temperature range of 50 to 350 deg, phase 2

Several two phase heat transfer fluids were tested in aluminum and carbon steel reflux capsules for over 25,000 hours at temperatures up to 300 C. Several fluids showed very good stability and would be useful for long duration heat transfer applications over the range 100 to 350 C. Instrumentation for the measurement of surface tension and viscosity were constructed for use with heat transfer fluids over the temperature range 0 to 300 C and with pressures from 0 to 10 atmospheres. The surface tension measuring device constructed requires less than a 1.0 cc sample and displays an accuracy of about 5 percent in preliminary tests, while the viscometer constructed for this program requires a 0.05 cc sample and shows an accuracy of about 5 percent in initial tests.

Saaski, E. W.↗
NASA NTRS1985

A swirl flow evaporative cold plate

A forced flow evaporative cold plate is under development for future application to the thermal bus concept being pursued by NASA for Space Station Thermal Control. The vaporizer is a swirl-flow device employing a spiral tube coil geometry sandwiched between conductive metal plates upon which electric components could be mounted. This concept is based on the inherent phase separation that occurs in a two phase stream in curvilinear flow. This is a zero 'g' design with one 'g' all-attitude capability and is capable of high heat transfer coefficients, good isothermality, and the ability to function at heat fluxes approaching 5w/sq cm on the cold plates (10w/sq cm on the tube wall) with Freon 114. The advantages of this design over other two phase evaporator approaches are high heat flux capability, simplified control requirements, insensitivity to micro-gravity oscillations, and inexpensive manufacturability. The program included design, fabrication, and test of such a cold plate utilizing an existing test stand developed for two-phase thermal management system (TPTMS) testing. Test results analysis and conclusions are included.

Niggemann, R. E.↗
NASA NTRS1986

Heat-Pipe Array for Large-Area Cooling

High rates of heat transfer anticipated. Prototype evaporative cold plate gathers waste heat from equipment mounted on it. Plate made by welding together flanges of several sections of heat pipe. Since plate separates liquid and vapor phases at inlet and outlet ports, eliminates complexities and uncertainties of two-phase flow in zero gravity. On earth, inlet valve enables plate to operate at relatively-large height differences with other plates in same system.

Edelstein, F.↗
NASA NTRS1985

Two-phase flows within systems with ambient pressure

In systems where the design inlet and outlet pressures are maintained above the thermodynamic critical pressure, it is often assumed that two phase flows within the system cannot occur. Designers rely on this simple rule of thumb to circumvent problems associated with a highly compressible two phase flow occurring within the supercritical pressure system along with the uncertainties in rotordynamics, load capacity, heat transfer, fluid mechanics, and thermophysical property variations. The simple rule of thumb is adequate in many low power designs but is inadequate for high performance turbomachines and linear systems, where two phase regions can exist even though outlet pressure is greater than critical pressure. Rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two phase zone can differ significantly from those for a single-phase zone. Using the Reynolds equation the angular velocity, eccentricity, geometry, and ambient conditions are varied to determine the point of two phase flow incipience.

Hendricks, R. C.↗