Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two Phase Flow”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Operando X-ray tomography and sub-second radiography for characterizing transport in polymer electrolyte membrane electrolyzer

Utilizing hydrogen gas as an energy carrier and enabling a hydrogen economy is a promising step towards decarbonization. Water electrolysis is a key process in hydrogen generation, but electrolyzers face a couple technological hurdles before widespread integration could occur. Understanding morphology evolution and transport processes in an operating polymer electrolyte membrane (PEM) electrolyzer is key to reducing cost and increasing efficiency. In this study, combined operando X-ray computed tomography and radiography are used to study transport and degradation in PEM electrolyzers subject to applied current densities. Tomography enables three-dimensional steady-state imaging but does not capture transport subtleties. Radiography is limited to two-dimensions but does capture transient phenomena. The tomography results depict degradation of the catalyst layer on the anode side of the electrolyzer. The rate of degradation was observed to increase as the applied current density increased. The catalyst particle detachment is due to autonomous propulsion during the oxygen evolution reaction. Particles that mechanically rip off the membrane were observed to be redeposited into the porous transport layer. During radiography, sub-second exposure time is required to capture the transient behavior of oxygen evolution, even at the lower current densities (50-200 mAcm -2 ). Various flow regimes were observed in the channel ranging from bubbly to slug flow depending on current density. In conclusion, experimental observations agree with a model that an increase in current density increases oxygen bubble diameter and decreases bubble residence time within the electrolyzer channel.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Experimental study of three-dimensional CO 2 -water drainage and fracture-matrix interactions in fractured porous media

Estimation of CO 2 storage capacity in fractured porous reservoirs requires a better understanding of the CO 2 -water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO 2 -water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. Here, the core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO 2 -water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO 2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO 2 saturations in the matrix increase with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO 2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO 2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO 2 storage efficiency in fractured porous media.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Interfacial gauge methods for incompressible fluid dynamics

Designing numerical methods for incompressible fluid flow involving moving interfaces, for example, in the computational modeling of bubble dynamics, swimming organisms, or surface waves, presents challenges due to the coupling of interfacial forces with incompressibility constraints. A class of methods, denoted interfacial gauge methods, is introduced for computing solutions to the corresponding incompressible Navier-Stokes equations. These methods use a type of "gauge freedom" to reduce the numerical coupling between fluid velocity, pressure, and interface position, allowing high-order accurate numerical methods to be developed more easily. Making use of an implicit mesh discontinuous Galerkin framework, developed in tandem with this work, high-order results are demonstrated, including surface tension dynamics in which fluid velocity, pressure, and interface geometry are computed with fourth-order spatial accuracy in the maximum norm. Applications are demonstrated with two-phase fluid flow displaying fine-scaled capillary wave dynamics, rigid body fluid-structure interaction, and a fluid-jet free surface flow problem exhibiting vortex shedding induced by a type of Plateau-Rayleigh instability. The developed methods can be generalized to other types of interfacial flow and facilitate precise computation of complex fluid interface phenomena.

42 ENGINEERING↗
DOE OSTI2019

Empirical derivation of a Nusselt number in a thermally stratified enclosure

A hydrophilic gelatin suspension possesses thermal properties that may affect the ratio of convective to conductive heat transfer across a boundary. This ratio, also known as Nusselt's number, can represent laminar or turbulent flow regimes. In this study, thermal transitioning empirical values were determined using a gelatin as a polymer gel surrogate. The crux of this study was to (a) quantify heat transfer properties of the gelatin as a suspended semi-solid and (b) understand thermal profiles with respect to particle size and loading. We can adduce the existence of a thermally induced: lower rise velocity, stratified Nusselt number, and a boundary layer.

42 ENGINEERING↗
DOE OSTI2016

The roll-up and merging of coherent structures in shallow mixing layers

The current study seeks a fundamental explanation to the development of two-dimensional coherent structures (2DCSs) in shallow mixing layers. A nonlinear numerical model based on the depth-averaged shallow water equations is used to investigate the temporal evolution of shallow mixing layers, where the mapping from temporal to spatial results is made using the velocity at the center of the mixing layers. The flow is periodic in the streamwise direction. Transmissive boundary conditions are used in the cross-stream boundaries to prevent reflections. Numerical results are compared to linear stability analysis, mean-field theory, and secondary stability analysis. Results suggest that the onset and development of 2DCS in shallow mixing layers are the result of a sequence of instabilities governed by linear theory, mean-field theory, and secondary stability theory. The linear instability of the shearing velocity gradient gives the onset of 2DCS. When the perturbations reach a certain amplitude, the flow field of the perturbations changes from a wavy shape to a vortical (2DCS) structure because of nonlinearity. The development of the vertical 2DCS does not appear to follow weakly nonlinear theory; instead, it follows mean-field theory. After the formation of 2DCS, separate 2DCSs merge to form larger 2DCS. In this way, 2DCSs grow and shallow mixing layers develop and grow in scale. The merging of 2DCS in shallow mixing layers is shown to be caused by the secondary instability of the 2DCS. Eventually 2DCSs are dissipated by bed friction. The sequence of instabilities can cause the upscaling of the turbulent kinetic energy in shallow mixing layers.

42 ENGINEERING↗
DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗
DOE OSTI2016

Roles of bulk viscosity on Rayleigh-Taylor instability: Non-equilibrium thermodynamics due to spatio-temporal pressure fronts

Direct numerical simulations of Rayleigh-Taylor instability (RTI) between two air masses with a temperature difference of 70 K is presented using compressible Navier-Stokes formulation in a non-equilibrium thermodynamic framework. The two-dimensional flow is studied in an isolated box with non-periodic walls in both vertical and horizontal directions. The non-conducting interface separating the two air masses is impulsively removed at t = 0 (depicting a heaviside function). No external perturbation has been used at the interface to instigate the instability at the onset. Computations have been carried out for rectangular and square cross sections. The formulation is free of Boussinesq approximation commonly used in many Navier-Stokes formulations for RTI. Effect of Stokes’ hypothesis is quantified, by using models from acoustic attenuation measurement for the second coefficient of viscosity from two experiments. Effects of Stokes’ hypothesis on growth of mixing layer and evolution of total entropy for the Rayleigh-Taylor system are reported. The initial rate of growth is observed to be independent of Stokes’ hypothesis and the geometry of the box. Following this stage, growth rate is dependent on the geometry of the box and is sensitive to the model used. As a consequence of compressible formulation, we capture pressure wave-packets with associated reflection and rarefaction from the non-periodic walls. The pattern and frequency of reflections of pressure waves noted specifically at the initial stages are reflected in entropy variation of the system.

42 ENGINEERING↗
DOE OSTI2025

First‐Order Empirical Interpolation Method for Real‐Time Solution of Parametric Time‐Dependent Nonlinear PDEs

ABSTRACT We present a model reduction approach for the real‐time solution of time‐dependent nonlinear partial differential equations (PDEs) with parametric dependencies. A major challenge in constructing efficient and accurate reduced‐order models for nonlinear PDEs is the efficient treatment of nonlinear terms. We address this by unifying the implementation of hyperreduction methods to deal with nonlinear terms. Furthermore, we introduce a first‐order empirical interpolation method (EIM) to provide an efficient approximation of the nonlinear terms in time‐dependent PDEs. We demonstrate the effectiveness of our approach on the Allen–Cahn equation, which models phase separation, and the Buckley–Leverett equation, which describes two‐phase fluid flow in porous media. Numerical results highlight the accuracy, efficiency, and stability of the proposed method compared with both the Galerkin–Newton approach and hyper‐reduced models using the standard EIM.

Nguyen, Ngoc Cuong [Center for Computational Engin↗
DOE OSTI2018

An efficient Bayesian data-worth analysis using a multilevel Monte Carlo method

Improving the understanding of subsurface systems and thus reducing prediction uncertainty requires collection of data. As the collection of subsurface data is costly, it is important that the data collection scheme is cost-effective. Design of a cost-effective data collection scheme, i.e., data-worth analysis, requires quantifying model parameter, prediction, and both current and potential data uncertainties. Assessment of these uncertainties in large-scale stochastic subsurface hydrological model simulations using standard Monte Carlo (MC) sampling or surrogate modeling is extremely computationally intensive, sometimes even infeasible. In this paper, we propose an efficient Bayesian data-worth analysis using a multilevel Monte Carlo (MLMC) method. Compared to the standard MC that requires a significantly large number of high-fidelity model executions to achieve a prescribed accuracy in estimating expectations, the MLMC can substantially reduce computational costs using multifidelity approximations. Since the Bayesian data-worth analysis involves a great deal of expectation estimation, the cost saving of the MLMC in the assessment can be outstanding. While the proposed MLMC-based data-worth analysis is broadly applicable, we use it for a highly heterogeneous two-phase subsurface flow simulation to select an optimal candidate data set that gives the largest uncertainty reduction in predicting mass flow rates at four production wells. Finally, the choices made by the MLMC estimation are validated by the actual measurements of the potential data, and consistent with the standard MC estimation. But compared to the standard MC, the MLMC greatly reduces the computational costs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Upgrade of the mechanically pumped CO 2 two-phase cooling system for the alpha magnetic spectrometer on the international space station

This study describes the successful upgrade of a mechanically pumped CO 2 two-phase cooling system in space by designing a new pump module for the Alpha Magnetic Spectrometer-02 on the International Space Station. Key factors for mission success are emphasized, including achieving high coolant filling accuracy within 10% of the target and maintaining system stability within ±1 °C. The impact of adding radiators to improve cooling efficiency is examined, and it is found that operating multi-radiators out-of-phase does not significantly affect system reliability. The centrifugal pump design is shown to allow for better lubricant circulation, while the in-house designed controller incorporates protective measures to prevent cavitation, overheating, and over-current. Further, this research advances the understanding of circulation loop systems and their upgrades in space and demonstrates the potential for extending the lifetime of space-borne mechanically pumped two-phase cooling systems.

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

Numerical investigation and experimental validation of water condensation in the gas diffusion layer with different properties

Liquid water in porous media hinders the transport of reactant to catalyst layer, where the electrochemical reactions occur, which in turn, reduces performance and lifetime of a proton exchange membrane fuel cell (PEMFC). In this study, two types of gas diffusion layer (GDL) materials are used to study the two-phase saturation distribution in a PEMFC. The correlations between the effective and anisotropic transport properties and porosity of the two GDLs are solved by the pore scale model. In addition, the anisotropic liquid water permeability and the relationship between saturation pressure and capillary pressure is determined by using the Lattice-Boltzmann method. Fuel cell performance as well as liquid water distribution in the GDL using neutron radiography are obtained for validation. Under wet conditions, the cell with Freudenberg GDL performs better than that using Toray GDL, especially at high current densities. The results of the two-phase model simulation show that the peak water saturation for Toray GDL occurs in the catalyst layer and can reach 40%, while the peak water saturation for Freudenberg GDL occurs inside the GDL directly under the ribs and only reaches 25% saturation. The combined results provide key insights to enable high power density operation of a PEMFC through GDL material optimization.

42 ENGINEERING↗
DOE OSTI2021

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↗
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 OSTI2018

On the Size Effect for Neutron Diffraction by Two-Phase Structures and Applicability of Composite Standards in Analysis of Residual Austenite in Steels

The possibility of using neutron diffraction in quantitative analysis of residual austenite in steels is investigated. According to the obtained results, some properties of neutron diffraction by two-phase structures remain size-independent, which makes it possible to use large composite standards. It is shown that this approach ensures admissible accuracy without using a priori coefficients.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

RELAP-7 Theory Manual

This document summarizes the physical models and mathematical formulations used in the RELAP-7 code. In summary, the MOOSE based RELAP-7 code development is an ongoing effort. The MOOSE framework enables rapid development of the RELAP-7 code. The developmental efforts and results demonstrate that the RELAP-7 project is on a path to success. This theory manual documents the main features implemented into the RELAP-7 code. Because the code is an ongoing development effort, this RELAP-7 Theory Manual will evolve with periodic updates to keep it current with the state of the development, implementation, and model additions/revisions.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2016

CTF Theory Manual

Coolant-Boiling in Rod Arrays|Two Fluids (COBRA-TF) is a thermal/ hydraulic (T/H) simulation code designed for light water reactor (LWR) vessel analysis. It uses a two-fluid, three-field (i.e. fluid film, fluid drops, and vapor) modeling approach. Both sub-channel and 3D Cartesian forms of 9 conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and had been used and modified by several institutions over the last few decades. COBRA-TF also found use at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Management Group (RDFMG) and has been improved, updated, and subsequently re-branded as CTF. As part of the improvement process, it was necessary to generate sufficient documentation for the open-source code which had lacked such material upon being adopted by RDFMG. This document serves mainly as a theory manual for CTF, detailing the many two-phase heat transfer, drag, and important accident scenario models contained in the code as well as the numerical solution process utilized. Coding of the models is also discussed, all with consideration for updates that have been made when transitioning from COBRA-TF to CTF. Further documentation outside of this manual is also available at RDFMG which focus on code input deck generation and source code global variable and module listings.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

High Fidelity BWR Fuel Simulations

This report describes the Consortium for Advanced Simulation of Light Water Reactors (CASL) work conducted for completion of the Thermal Hydraulics Methods (THM) Level 3 milestone THM.CFD.P13.03: High Fidelity BWR Fuel Simulation. High fidelity computational fluid dynamics (CFD) simulation for Boiling Water Reactor (BWR) was conducted to investigate the applicability and robustness performance of BWR closures. As a preliminary study, a CFD model with simplified Ferrule spacer grid geometry of NUPEC BWR Full-size Fine-mesh Bundle Test (BFBT) benchmark has been implemented. Performance of multiphase segregated solver with baseline boiling closures has been evaluated. Although the mean values of void fraction and exit quality of CFD result for BFBT case 4101-61 agreed with experimental data, the local void distribution was not predicted accurately. The mesh quality was one of the critical factors to obtain converged result. The stability and robustness of the simulation was mainly affected by the mesh quality, combination of BWR closure models. In addition, the CFD modeling of fully-detailed spacer grid geometry with mixing vane is necessary for improving the accuracy of CFD simulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗