Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “two-phase flows”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

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↗
DOE OSTI2025

Non-dimensional performance and safety parameters for heat pipes

The use of heat pipes in safety-critical systems such as nuclear microreactors dictates the development of generalized, practical, scalable performance and safety parameters. Traditional dimensional metrics, while informative, lack the universality required for comparative analysis across varying designs and operating regimes. Here, this work introduces a comprehensive set of non-dimensional parameters to characterize heat pipe performance and safety, including capillary performance, effective thermal conductivity, response time, exergetic efficiency, allowable temperature gradients, allowable rate of temperature change, priming coefficients, and factor of safety. A reference heat pipe design representative of microreactor applications was analyzed via the developed parameters using both traditional analytical models and Sockeye simulations under transient and steady-state conditions. Sodium, potassium, and water were evaluated as working fluids to demonstrate the applicability of the framework across a broad temperature range. The proposed non-dimensional parameters effectively captured key thermal-hydraulic behaviors and safety concerns, as was demonstrated via Sockeye simulations. This framework supports the development of design optimization strategies, operational protocols, and safety assurance practices for advanced reactor systems and other high-reliability applications.

42 - ENGINEERING↗
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 OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
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 NTRS2026

Center for High-Efficiency Electrical Technologies for Aircraft: Phase I Final Report

Under the University Leadership Initiative (ULI), the Center for High-Efficiency for Electrical Technologies (CHEETA) was established to develop and mature early-stage technologies pertaining to hydrogen-electric power and energy systems for aircraft. In particular, integration of these technologies on an aircraft system is envisioned to leverage the high specific energy content of liquid hydrogen (LH2) with fuel cell energy conversion and an electrically driven ducted fan system to provide an ultra-efficient propulsion drivetrain. For this concept, the LH2 system is not just used as an energy storage mechanism, but also as a cryogen to enable highly efficient superconducting electric systems. The end result of this concept is an integrated aircraft system with a quiet, efficient propulsion architecture that produces zero CO2, NOx, SOx, and particulate matter emissions at the vehicle level.

Hydrogen↗
DOE OSTI2025

Evaporator Temperature Transient Testing of a High-Performance Sodium Filled Heat Pipe

Heat pipes are two-phase heat transfer devices that enable passive removal of heat from the reactor core to the power conversion system in heat pipe-cooled microreactor designs. Experimental investigations of heat pipe transients are needed for technology demonstration, verification and validation of numerical codes, and the establishment of regulatory requirements. The Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL) serves as a platform for evaluating the dynamic response of high-temperature heat pipes under a variety of operating conditions. The present work details the experimental investigation of a high-performance, defined as over 2 kW sodium heat pipe subjected to rapid input power fluctuations induced by sudden changes in the evaporator temperature setpoint. In addition, the heat pipe was subjected to an asymmetrical heat load where a subset of heaters operated at 30% and 70% below their nominal power. These experimental conditions were chosen to simulate thermal and operational stresses expected to be encountered in microreactors to provide data on heat pipe behavior during such important transient events. Key data and performance metrics, including time series of temperatures and strains, axial temperature profiles, thermal response times, and heat transfer capabilities, the thermal output over thermal input, were reported and discussed. The results highlight the resilience of heat pipes, revealing their potential to maintain thermal stability and efficiency under varying power loads. Lastly, the paper concludes with a discussion on the significance of the results and their implications for future research.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1979

Proceedings of the 25th International Instrumentation Symposium, May 7-10, 1979, Anaheim, California

The papers deal essentially with the instrumentation needs of the transportation industry, machinery diagnostics, and energy research. The sessions covered include: acoustic emission, reentry vehicle ground and flight testing, performance of nose tip materials, instrumentation for solar heating and cooling systems, vibration measurement, data acquisition and analysis, system for aerospace pressure measurements, wind tunnel instrumentation and control, energy source instrumentation, nondestructive testing, pressure and flow measurements, automated test facilities, strain measurements, electrooptical instrumentation, thermal measurements, microcomputer applications in instrumentation and automation, flight testing, blast pressure measurement, two-phase liquid/gas mass flow measurement, and machinery monitoring and instrumentation.

Conferences↗
NASA NTRS2002

Spike-Nosed Bodies and Forward Injected Jets in Supersonic Flow

The paper contains new numerical simulation and experimental test results of blunt body drag reduction using thin spikes mounted in front of a body and one- or two-phase jets injected against a supersonic flow. Numerical simulations utilizing the NASA CFL3D code were conducted at the Hampton University Fluid Mechanics and Acoustics Laboratory (FM&AL) and experimental tests were conducted using the facilities of the IM/MSU Aeromechanics and Gas Dynamics Laboratory. Previous results were presented at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference. Those results were based on some experimental and numerical simulation tests for supersonic flow around spike-nosed or shell-nosed bodies, and numerical simulations were conducted only for a single spike-nosed or shell-nosed body at zero attack angle, alpha=0. In this paper, experimental test results of gas, liquid and solid particle jet injection against a supersonic flow are presented. In addition, numerical simulation results for supersonic flow around a multiple spike-nosed body with non-zero attack angles and with a gas and solid particle forward jet injection are included. Aerodynamic coefficients: drag, C(sub D), lift, C(sub L), and longitudinal momentum, M(sub z), obtained by numerical simulation and experimental tests are compared and show good agreement.

Gilinsky, M.↗
DOE OSTI2025

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2026

Flow Boiling Pressure Drop Characteristics of Next-generation Refrigerants in a Micro‑fin Copper Tube

This paper presents experimental frictional pressure-drop data for flow boiling of R-410A, R-134a, and next-generation alternatives R-454C, R-455A, R-1234yf, and R-1234ze(E) in a horizontal micro-fin copper tube. Tests were conducted over a range of mass fluxes and evaporation temperatures to characterize refrigerant-dependent two-phase pressure-drop behavior. The results show that frictional pressure gradient increased with mass flux and vapor quality and generally increased as evaporation temperature decreased, with liquid viscosity strongly affecting the observed trends. Among the evaluated correlations, the Goto et al. (2001) model gave the best overall agreement with the measurements before optimization. Further optimization of the Kuo and Wang (1996) and Goto et al. (2001) models reduced the overall mean absolute deviation to below 15%, with the optimized Goto (2001) model providing the most consistent predictions across all six refrigerants. The results support improved pressure-drop prediction and evaporator design for next-generation refrigerants in micro-fin tubes.

Hu, Yifeng [ORNL] (ORCID:0000000242875185)↗
DOE OSTI2025

Probing Nanorod Assembly and Dynamics in Polymer Nanocomposites in Equilibrium and Shear

Coarse-grained molecular dynamics simulations are used to examine the structure and dynamics of nanorod assemblies in polymer melts under equilibrium and simple shear. We show that as the concentration of nanorods increases, there is a transition from an isotropic phase to a two-phase region in which the nanorods phase separate into a dilute phase and dense bundles of hexagonally packed nanorods. The onset of the two-phase region is below that predicted by Onsager theory, which we attribute to an effective increase in the diameter of nanorods due to a layer of polymer bound to the rod surfaces. Equilibrium simulations show that increasing polymer chain length N enhances nanorod bundling at fixed nanorod concentration. Simulations of systems undergoing simple shear show that flow enhances nanorod alignment and bundling relative to those of equilibrium systems with similar properties. Finally, simulations reveal that increasing nanorod length enhances nanorod alignment under shear at equivalent shear rates but not at equivalent Péclet numbers. Overall, our simulations highlight that polymer matrix-nanorod attraction (i.e., chemistry), shear rate, and matrix chain length are desirable design variables to control the structure of nanorod-containing soft materials under simple shear.

diffusion↗
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 OSTI2026

Prototype Modeling for a Light-Trapping Planar-Cavity Enclosed Particle Solar Receiver

Concentrating solar thermal (CST) systems present a promising avenue for affordable and reliable energy production. Solar receivers are key components that determine the efficiency and longevity of these systems. Particle-based solar receivers have emerged as a compelling alternative to traditional technologies, offering several advantages that address limitations in current CST systems. This is especially true as next-generation CST technologies target applications including electricity generation, thermochemical processes, and industrial process heat, many of which necessitate higher operating temperatures than current commercial molten salt systems. Molten-salt thermal energy storage (TES) systems, commonly used in CSP, face challenges related to freezing and corrosion. Particle-based TES systems, in contrast, do not experience these issues, as particles are stable at high temperatures, exceeding 1000 degrees Celsius. This capability allows for a wider range of applications, including those requiring higher temperatures for industrial processes and efficient electricity generation. A novel innovation in particle-based solar receiver technology is the light-trapping planar cavity receiver (LTPCR) configuration developed by NREL. The LTPCR design consists of small cavity-like structures using opaque planar surfaces, enabling efficient capture and absorption of solar energy. A high incident flux concentration at the cavity aperture is absorbed on the receiver walls, and subsequently transferred to particles on the inside of cavities. The particles flow through the system, forming a fluidized bed inside of the receiver panels, effectively capturing the absorbed solar heat. Air is used as a fluidizing medium in this process to enhance particle heat transfer and mixing. The effectiveness of this design lies in its ability to manage solar flux conditions and ensure high solar-to-thermal receiver efficiency. A 100-kW prototype is currently being tested at the King Saud University in Saudi Arabia to assess the receiver performance. A range of modeling analyses for the optical, thermal, and mechanical effects were conducted to assess the performance of the receiver under on-sun conditions. The solar flux resulting from the KSU heliostat field was modeled using NREL SolTrace software and produced up to 600 kW/m2 at the receiver aperture. The solar flux absorbed on the receiver walls was then used within a computational fluid dynamics (CFD) model to predict wall temperature distributions along with radiation and convection loss. A two-phase CFD model was developed for the fluidized bed of silica sand inside the receiver panels to predict local wall-to-particle heat transfer coefficients, particle temperature distributions, and outlet temperature of the particles. We have also conducted analyses to understand the thermomechanical behavior of these innovative enclosed light-trapping solar receivers optimized for particle heating. We used finite element analysis (FEA) to predict the receiver's performance using temperature distributions obtained from CFD and based on the resulting stress profiles, evaluated creep-fatigue damage with a goal of achieving a 30-year service life. Analysis showed a significant impact of the particle-to-wall heat transfer coefficients (HTCs) on receiver performance, with higher HTCs resulting in reduced stress and increased lifespan. For instance, when using Inconel 740H, increasing the HTC from 800 W/m2 K to 1400 W/m2 K increased the creep life from 4,000 hours to over 100,000 hours. This highlights the importance of understanding and optimizing heat transfer in the design of high-efficiency receivers.

14 SOLAR ENERGY↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗
DOE OSTI2026

A Kinetic Model of the Long-Term Corrosion of Glass-Ceramic Materials

Multiphase waste forms show promise for increased waste loading and for the ability to dispose of contaminated solid and particulate waste through direct densification. However, achieving predictive capability for long-term durability of multiphase waste forms, and thus assessing their possible deployment, requires expanding the current, limited knowledge base. Here, we describe the development of a corrosion model of a two-phase waste form consisting of crystals of known volume fraction embedded in a glass matrix. This model accounts for the dissolution of both the crystalline and glass phases as well as the hydration of the glass phase through an ion exchange reaction. Because of the large difference in solubility between the two phases, the reactive surface of the crystalline phase is a function of the extent of dissolution of the glass phase in this model. Model parameterization was performed using corrosion data, such as from single-pass flow-through tests, for the individual phases. The parameterized corrosion model was evaluated against static dissolution test data for a glass-ceramic multiphase waste form. This evaluation demonstrated the model’s ability to reproduce the time-dependent release of key tracers of glass and crystalline phase dissolution. Hence, the development of a kinetic model provides a pathway for long-term durability predictions and thus the use of multiphase waste forms in nuclear cleanup missions.

Kerisit, Sebastien N.↗