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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
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 OSTI2024

Computational Investigation of Fuel Dispersal Phenomena during Large-Break Loss of Coolant Accident in Light-Water Reactors

In the event of cladding rupture, which could occur in light water reactor fuel assemblies during a loss-of-coolant accident (LOCA), fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod. This expulsion of fragmented fuel particles, referred as fuel dispersal, is the subject of investigation to evaluate the safety implications of increasing fuel burnup in light water reactors, with a specific focus on fuel fragmentation, relocation, and dispersal. Particle trajectories and the resulting mass distribution of the settled particles within the reactor pressure vessel can pose a long-term cooling challenge for the reactor core. The fuel dispersal phenomenon is significantly influenced by the ejection characteristics of the fuel fragments, as well as the size and shape of the cladding rupture and fuel rod depressurization history during LOCA transients. In this study, the transport of fuel particles within a scaled 5 × 5 lattice of a pressurized water reactor rod bundle geometry is modeled through a two-fluid Eulerian framework that treats the gas and solid phases as interpenetrating continua. The required boundary conditions are evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considers solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum and energy exchange between the gas and solid phases. The simulation results provide the volume fraction of the solids settled on the bottom surface of the fuel bundle, quantifying the deposition within the bundle geometry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗
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

Velocity-space Origins of the Pressure–Strain Interaction in Multipopulation Distributions and Its Application to Magnetic Reconnection

A forefront research question is how energy evolves in weakly collisional plasmas for which departures from local thermodynamic equilibrium (LTE) are significant. The standard approach is studying the terms in the non-LTE energy evolution equation derived by taking the second moment of the Boltzmann equation, but the resultant fluid metrics do not retain information about which particles at which velocities drive energy evolution. A widely studied channel for internal energy density evolution is the pressure–strain interaction. Here, we employ the kinetic pressure–strain, a phase-space diagnostic whose velocity-space integral recovers the pressure–strain interaction to disambiguate the contributions to the pressure–strain interaction from disparate particle populations in composite phase-space densities. We develop phase-space analogs of the pressure–strain interaction decompositions to provide the phase-space origins of normal versus sheared flow. We introduce the “kinetic strain-rate” tensor, the phase-space analog of the strain-rate tensor, which we argue is needed to interpret the phase-space origins of the pressure–strain interaction. To demonstrate the utility of these quantities, we investigate them for composite electron distributions near the electron diffusion region in two-dimensional particle-in-cell simulations of antiparallel symmetric magnetic reconnection. We find that the phase-space-based diagnostics isolate the roles of distinct populations. These results contribute to a growing body of work providing new methods for quantifying phase-space energy evolution for a broad array of processes, from magnetic reconnection to collisionless shocks and turbulence, opening new pathways for answering longstanding problems of particle energization in weakly collisional plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2024

Self-Focused Pulse Propagation Is Mediated by Spatiotemporal Optical Vortices

We show that the dynamics of high-intensity laser pulses undergoing self-focused propagation in a nonlinear medium can be understood in terms of the topological constraints imposed by the formation and evolution of spatiotemporal optical vortices (STOVs). STOVs are born from point phase defects on the sides of the pulse nucleated by spatiotemporal phase shear. These defects grow into closed loops of spatiotemporal vorticity that initially exclude the pulse propagation axis, but then reconnect to form a pair of toroidal vortex rings that wrap around it. STOVs constrain the intrapulse flow of electromagnetic energy, controlling the focusing-defocusing cycles and pulse splitting inherent to nonlinear pulse propagation. Here, we illustrate this in two widely studied but very different regimes, relativistic self-focusing in plasma and nonrelativistic self-focusing in gas, demonstrating that STOVs mediate nonlinear propagation irrespective of the detailed physics.

42 ENGINEERING↗
DOE OSTI2024

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗
DOE OSTI2024

The Impact of Wettability on the Co-moving Velocity of Two-Fluid Flow in Porous Media

Abstract The impact of wettability on the co-moving velocity of two-fluid flow in porous media is analyzed herein. The co-moving velocity, developed by Roy et al. (Front Phys 8:4, 2022), is a novel representation of the flow behavior of two fluids through porous media. Our study aims to better understand the behavior of the co-moving velocity by analyzing simulation data under various wetting conditions. We analyzed 46 relative permeability curves based on the Lattice–Boltzmann color fluid model and two experimentally determined relative permeability curves. The analysis of the relative permeability data followed the methodology proposed by Roy et al. (Front Phys 8:4, 2022) to reconstruct a constitutive equation for the co-moving velocity. Surprisingly, the coefficients of the constitutive equation were found to be nearly the same for all wetting conditions. On the basis of these results, a simple approach was proposed to reconstruct the relative permeability of the oil phase using only the co-moving velocity relationship and the relative permeability of the water phase. This proposed method provides new information on the interdependence of the relative permeability curves, which has implications for the history matching of production data and the solution of the associated inverse problem. The research findings contribute to a better understanding of the impact of wettability on fluid flow in porous media and provide a practical approach for estimating relative permeability based on the co-moving velocity relationship, which has never been shown before.

Engineering↗
DOE OSTI2024

Self-assembled reconfigurable pump architectures via magnetic colloidal swarms

Self-assembled swarms of interactive active units, which are adaptive and dynamically reconfigurable to accommodate different functionalities, represent a promising platform for the development of next-generation robotics. Here, we utilize the emergent collective behavior of active magnetic colloids confined in quasi-two-dimensional arrays of overlapping wells to demonstrate the self-organization of a colloidal swarm into a dynamic pump architecture capable of controlled transport of passive cargo particles. This dynamic architecture provides a global unidirectional looping flow pattern along the entire length of the system. We show that the flow direction of the dynamic swarm-based pump can be externally controlled by a phase shift of a driving magnetic field energizing the swarm. The experimental observations are supported by computational modeling based on phenomenological coarse-grained particle dynamics coupled to shallow-water Navier-Stokes hydrodynamics. In conclusion, our findings demonstrate how the emergent collective behavior of a swarm can be orchestrated into a desired functionality by exploiting the interplay between activity and confinement potentials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Rheology improvement for silicon nitride and resin slurries for vat photopolymerization printing and sintering

This work presents the formulation, rheological characterization, and sintering of silicon nitride (Si₃N₄) slurries for vat photopolymerization (VPP) using digital light processing (DLP). Bimodal Si₃N₄ powder was dispersed into commercial photopolymer resin using two different dispersants with opposing effects on surface charge, resulting in varied slurry stability and flow behavior. Slurries were engineered to exhibit shear-thinning behavior suitable for VPP, and their flow properties were quantified using a power-law fluid model. The formulations achieved cure depths of approximately 40 µm and enabled printing of green bodies. Post-processing included thermal debinding and liquid-phase sintering, yielding primarily β-Si₃N₄ with a minor Y-Si-Al-O-N glass phase. The sintered parts reached ~85% of theoretical density and demonstrated a flexural strength of ~330 MPa. Microstructural analysis revealed closed porosity along with some defects related to powder agglomeration and interlayer adhesion. These findings provide insights into slurry formulation strategies for additive manufacturing of high-performance non-oxide ceramics.

Aerospace engineering↗
DOE OSTI2024

LCLS-II Helium Refrigeration System Commissioning Results

SLAC National Accelerator Laboratory has upgraded to LCLS-II, featuring a 4 GeV superconducting linear accelerator composed of 37 cryomodules and two large helium refrigeration systems with a cooling capacity of 4 kW at 2.0 K. The LCLS-II Helium Refrigeration System (HRS) consists of two compressor stations, each with a power of approximately 4.5 MW, two 4.5 K cold boxes, each with a power of 18 kW equivalent at 4.5 K, and two sets of cold compressors that can each produce a flow of 230 g/s at 31 mbar (2.0K). Performance tests of the HRS were meticulously planned and successfully carried out, with results demonstrating that it exceeded the process requirements for LCLS-II operations. This paper provides a detailed presentation of the LCLS-II HRS performance and the challenges encountered during the commissioning phase.

42 ENGINEERING↗
DOE OSTI2024

MENT-Flow: maximum-entropy phase space tomography using normalizing flows

Generative models can be trained to reproduce low-dimensional projections of high-dimensional phase space distributions. Normalizing flows are generative models that parameterize invertible transformations, allowing exact probability density evaluation and sampling. Consequently, flows are unbiased entropy estimators and could be used to solve the high-dimensional maximum-entropy tomography (MENT) problem. In this work, we evaluate a flow-based MENT solver (MENT-Flow) against exact maximum-entropy solutions and Minerbo's iterative MENT algorithm in two dimensions.

Hoover, Austin↗
DOE OSTI2026

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗
DOE OSTI2025

Temperature Measurements in Hypersonic Wind Tunnels via Femtosecond Coherent Anti-Stokes Raman Scattering

A femtosecond coherent anti-Stokes Raman scattering (fs CARS) instrument is developed to perform gas-phase thermometry in cold-flow hypersonic wind tunnels. Measurements are reported for Mach 8 and 14 pure-nitrogen flows. The fs CARS instrument includes a 100 fs pump/Stokes pulse and a spectrally narrow probe pulse from a second harmonic bandwidth compressor. Important experimental considerations such as limits on the pump/Stokes pulse energy are discussed. The fs CARS focusing and collimating optics are mounted on a two-axis translation stage system to scan the measurement location during a 30 second wind tunnel run. Single-laser-shot rotational CARS spectra are recorded at the laser repetition rate of 1 kHz in the wind tunnel freestream and near simple cone models. Spectral fitting is used to determine quantitative gas temperatures. Freestream temperatures at Mach 8 and 14 spanned ranges of 40–75 and 35–50 K, respectively, depending on tunnel operating conditions. Temperature variations across the central 100 mm span of the wind tunnel were quantified. Measured temperature jumps across conical bow shocks from various models varied by less than 1% from predicted values. Hypersonic boundary layer measurements were demonstrated. In conclusion, these measurements illustrate the utility and robustness of this instrument for the study of complex fluid flow phenomena in challenging ground test facilities.

Aerodynamics↗