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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Incompressible multiphase flow and encapsulation simulations using the moment-of-fluid method

A moment-of-fluid method is presented for computing solutions to incompressible multiphase flows in which the number of materials can be greater than two. In this work, the multimaterial moment-of-fluid interface representation technique is applied to simulating surface tension effects at points where three materials meet. The advection terms are solved using a directionally split cell integrated semi-Lagrangian algorithm, and the projection method is used to evaluate the pressure gradient force term. The underlying computational grid is a dynamic block-structured adaptive grid. The new method is applied to multiphase problems illustrating contact-line dynamics, triple junctions, and encapsulation in order to demonstrate its capabilities. Examples are given in two-dimensional, three-dimensional axisymmetric (R–Z), and three-dimensional (X–Y–Z) coordinate systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Gradient flow based phase-field modeling using separable neural networks

Allen–Cahn equation is a reaction–diffusion equation and is widely used for modeling phase separation. Machine learning methods for solving the Allen–Cahn equation in its strong form suffer from inaccuracies in collocation techniques, errors in computing higher-order spatial derivatives, and the large system size required by the space–time approach. To overcome these challenges, we propose solving the gradient flow of the Ginzburg–Landau free energy functional, which is equivalent to the Allen–Cahn equation, thereby avoiding the second-order spatial derivatives associated with the Allen–Cahn equation. A minimizing movement scheme is employed to solve the gradient flow problem, eliminating the complexities of a space–time approach. We utilize a separable neural network that efficiently represents the phase field through low-rank tensor decomposition. As we use the minimizing movement scheme to numerically solve the gradient flow problem, we thus, refer to the proposed method as the Separable Deep Minimizing Movement (SDMM) method. The evaluation of the functional in the minimizing movement scheme using the Gauss quadrature technique bypasses the inaccuracies associated with collocation techniques traditionally used to solve partial differential equations. A hyperbolic tangent transformation is introduced on the phase field prior to the evaluation of the functional to ensure that it remains strictly bounded within the values of the two phases. For this transformation, theoretical guarantee for energy stability of the minimizing movement scheme is established. Our results suggest that this transformation helps to improve the accuracy and efficiency significantly. The proposed method resolves the challenges faced by state-of-the-art machine learning techniques, outperforming them in both accuracy and efficiency. It is also the first machine learning method to achieve an order of magnitude speed improvement over the finite element method. In addition to its formulation and computational implementation, several case studies illustrate the applicability of the proposed method.

42 ENGINEERING↗
DOE OSTI2018

Development of tearing instability in a current sheet forming by sheared incompressible flow

Sweet–Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal magnetohydrodynamics, a coordinate transformation reveals that the flow compresses and stretches perturbations. Here, a multiple scale formulation describes how linear tearing mode theory can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Phases of 2D massless QCD with qubit regularization

We investigate the possibility of reproducing the continuum physics of 2D S U ( N ) gauge theory coupled to a single flavor of massless Dirac fermion using qubit regularization. The continuum theory is described by N free fermions in the ultraviolet (UV) and a coset Wess-Zumino-Witten (WZW) model in the infrared (IR). In this work, we first explore how well these features can be reproduced using the Kogut-Susskind (KS) Hamiltonian with a finite-dimensional link Hilbert space and a generalized Hubbard coupling. We do this by analyzing the renormalization group (RG) flow diagram of the continuum theory and identifying important phases of the theory. Using strong coupling expansions, we show that our lattice model exhibits a gapped dimer phase and a spin-chain phase. Furthermore, for N = 2 , using tensor network methods, we show that there is a second-order phase transition between these two phases, which we identify as the critical surface of the continuum theory that connects the IR and UV fixed points. In the IR, we identify the critical theory at the transition as the expected S U ( 2 ) 1 WZW model. Lastly, we argue that modifications of our model may allow the study of the UV physics of free fermions. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Integrated Smart Distribution RD&D Project Final Technical Report

Title XIII of the Energy Independence and Security Act of 2007 (EISA), tasked OE with accelerating the commercialization, deployment and integration of advanced digital systems that are needed to modernize the nation’s electricity delivery network for enhanced interoperability and cyber security. To support this mandated task, a goal of the SGR&D Program is to develop technology with the ability to dynamically optimize grid operations and resources and incorporate demand response and consumer participation. The majority of historical research and development in the domain of distribution management systems (DMS) had focused on classical, radial and passive networks with one-way power flow from transmission to distribution systems and thus to the end customer. With the focus on customers, a bottom-up approach required a completely new thinking in real-time operations, and automation, of smart distribution systems. Thus, the GE team embarked on a project with a plan of two distinct phases. Work during the first phase concentrated on research and development (R&D) for the first three years. The second phase, lasting approximately two years, accomplished the demonstration of the use cases developed during the first phase, with further sharing and dissemination of the research at industry events. The consolidated technical scope of the project can be summarized into the following six areas: Advanced distribution modeling capability to accurately simulate/model smart grid operations. Accurate representation of the distribution system in real- or near real-time (capture real-time topology). Interoperability with, and seamless communication between, other management systems and databases used by the utility Simulation of distribution systems based on real-time operational planning to analyze the benefits of smart grid assets. Integration of network, market, and renewable resource models for next generation DMS. Management and forecasting of DER (Distributed Generation, Storage, Demand Response).

14 SOLAR ENERGY↗
DOE OSTI2019

Development of Tearing Instability in a Current Sheet Forming by Sheared Incompressible Flow

Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvenic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; Rutherford 1973) is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal MHD, a coordinate transformation reveals that the flow compresses and stretches perturbations. A multiple scale formulation describes how linear tearing mode theory (Furth et al. 1963; Coppi et al. 1976) can be applied to an equilibrium changing under flow, showing that tearing eigenfunctions grow slowly due to the flow and undergo exponential growth at a rate given by standard scalings with time dependence added. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Investigation on the two-stage active magnetic regenerative refrigerator for liquefaction of hydrogen

An active magnetic regenerative refrigerator (AMRR) is expected to be useful for hydrogen liquefaction due to its inherent high thermodynamic efficiency. Because the temperature of the cold end of the refrigerator has to be approximately liquid temperature, a large temperature span of the active magnetic regenerator (AMR) is indispensable when the heat sink temperature is liquid nitrogen temperature or higher. Since magnetic refrigerants are only effective in the vicinity of their own transition temperatures, which limit the temperature span of the AMR, an innovative structure is needed to increase the temperature span. The AMR must be a layered structure and the thermophysical matching of magnetic field and flow convection effects is very important. In order to design an AMR for liquefaction of hydrogen, the implementation of multi-layered AMR with different magnetic refrigerants is explored with multi-staging. In this paper, the performance of the multi-layered AMR using four rare-earth compounds (GdNi{sub 2}, Gd{sub 0.1}Dy{sub 0.9}Ni{sub 2}, Dy{sub 0.85}Er{sub 0.15}Al{sub 2}, Dy{sub 0.5}Er{sub 0.5}Al{sub 2}) is investigated. The experimental apparatus includes two-stage active magnetic regenerator containing two different magnetic refrigerants each. A liquid nitrogen reservoir connected to the warm end of the AMR maintains the temperature of the warm end around 77 K. High-pressure helium gas is employed as a heat transfer fluid in the AMR and the maximum magnetic field of 4 T is supplied by the low temperature superconducting (LTS) magnet. The temperature span with the variation of parameters such as phase difference between magnetic field and mass flow rate of magnetic refrigerants in AMR is investigated. The maximum temperature span in the experiment is recorded as 50 K and several performance issues have been discussed in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Flow and thermal modelling of the argon volume in the DarkSide-20k TPC

The DarkSide-20k dark matter experiment, currently under construction at LNGS, features a dual-phase time projection chamber (TPC) with a ∼ 50 t argon target from an underground well. At this scale, it is crucial to optimise the argon flow pattern for efficient target purification and for fast distribution of internal gaseous calibration sources with lifetimes of the order of hours. To this end, we have performed computational fluid dynamics simulations and heat transfer calculations. The residence time distribution shows that the detector is well-mixed on time-scales of the turnover time (∼ 40 d). Notably, simulations show that despite a two-order-of-magnitude difference between the turnover time and the half-life of 83m Kr of 1.83 h, source atoms have the highest probability to reach the centre of the TPC 13 min after their injection, allowing for a homogeneous distribution before undergoing radioactive decay. We further analyse the thermal aspects of dual-phase operation and define the requirements for the formation of a stable gas pocket on top of the liquid. We find a best-estimate value for the heat transfer rate at the liquid-gas interface of 62 W with an upper limit of 144 W and a minimum gas pocket inlet temperature of 89 K to avoid condensation on the acrylic anode. This study also informs the placement of liquid inlets and outlets in the TPC. The presented techniques are widely applicable to other large-scale, noble-liquid detectors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. 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 relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
DOE OSTI2023

Sequence-dependent self-assembly of supramolecular nanofibers in periodic dynamic block copolymers

Block copolymers (BCPs) can spontaneously self-assemble into various nanostructured morphologies which enable their diverse applications in selective membranes, polymer electrolytes, and optoelectronics. To expand the range of nanostructures accessible to block copolymers, we designed dynamic block copolymers (DBCPs) that combine the phase separation of traditional block copolymers with the supramolecular self-assembly of periodic dynamic polymers. We demonstrate that DBCPs synthesized with a periodic block sequence self-assemble into high aspect ratio supramolecular nanofibers with well-ordered PEG and PDMS domains, in contrast to those synthesized with a random sequence which do not form nanofibers but have disordered morphologies. The periodicity of the block sequence ensures regular placement of dynamic bonds along the polymer chain, which enables stacking of the hydrogen bonding units into ordered 1D supramolecular nanofibers and delays the onset of terminal flow by up to 60 °C compared to the random block sequence. The hierarchically assembled supramolecular nanofibers display complex mechanical and thermal phase behavior arising from the interplay between phase separation of the dissimilar polymer backbones and supramolecular interactions between dynamic bonds. Despite identical bulk composition, the periodic DBCPs demonstrate ionic conductivity values over two orders of magnitude higher than their random counterparts due to the formation of well-ordered, interconnected, high aspect ratio ion-transporting PEG domains. In conclusion, these results highlight the potential for DBCPs as an emerging material platform to achieve advanced, self-assembled nanostructured morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Hydrodynamic dispersion in thin channels with micro-structured porous walls

Flowand transport within porous- and microtextured-walled channels is relevant to a number of natural and industrial processes. Designing and optimizing the topology of the pores and/or microstructure to achieve target performance at the system scale (or macroscale) is still an open question. In this work, we study whether hydrodynamic dispersion in microfluidic channels with walls structured by obstacles can be modeled by dispersion in channels with porous walls described as continuous porous media of zero or finite permeability. We perform single phase microfluidic non-reactive flow experiments in channels embedded in micropatterns with different topologies. Specifically, we focus on transverse riblets and arrays of pillars as examples of impermeable and permeable obstructions, respectively. We compare the experimental results with three models: 3D pore-scale simulations which resolve the micropattern geometry explicitly and two upscaled models which treat the micropattern as a continuum of zero or finite permeability. This study demonstrates that polydimethylsiloxane micromodels with appropriately patterned surfaces can be successfully employed to validate various continuum-scale modeling approximations in different physical regimes, identified by the order of magnitude of the Peclet number and the obstruction permeability.

Ling, Bowen↗
DOE OSTI2019

Diurnal Forcing and Phase Locking of Gravity Waves in the Maritime Continent

Abstract An unfiltered zonal Hovmöller depiction of rainfall in the Maritime Continent (MC) reveals remarkable spatiotemporal continuity of zonally propagating disturbances with a diurnal period, which endure over multiple days and propagate faster than the individual convective storms they coupled with. This phenomenon and its sensitivity to the Madden–Julian oscillation (MJO) during the 2011/12 Dynamics of the MJO (DYNAMO) field campaign is examined here through a well-validated, convection-permitting model simulation conducted on a large domain. We find that these disturbances are zonally propagating diurnal gravity waves excited by vigorous nocturnal mesoscale convective systems over Sumatra and Borneo. These gravity waves are diurnally phase locked: their wavelength very closely matches the distance between these two islands (~1500 km), while their particular zonal phase speed (~±17 m s−1) allows them to propagate this distance in one diurnal cycle. We therefore hypothesize that these waves are amplified by resonant interaction due to diurnal phase locking. While these zonal gravity waves decouple from convection once beyond the MC, their divergent flow signature endures well across the Indian Ocean, provoking the notion that they may influence rainfall at far remote locations. The exact controls over this zonal phase speed remain uncertain; we note, however, that it is roughly consistent with diurnal offshore-propagating modes documented previously. Further study is required to tie this down, and more generally, to understand the sensitivity of these modes to background flow strength and the geography of the MC.

Ruppert, Jr., James H.↗
DOE OSTI2016

ACCRETION FLOW DYNAMICS OF MAXI J1836-194 DURING ITS 2011 OUTBURST FROM TCAF SOLUTION

The Galactic transient X-ray binary MAXI J1836-194 was discovered on 2011 August 29. Here we make a detailed study of the spectral and timing properties of its 2011 outburst using archival data from the RXTE Proportional Counter Array instrument. The evolution of accretion flow dynamics of the source during the outburst through spectral analysis with Chakrabarti–Titarchuk’s two-component advective flow (TCAF) solution as a local table model in XSPEC. We also fitted spectra with combined disk blackbody and power-law models and compared it with the TCAF model fitted results. The source is found to be in hard and hard-intermediate spectral states only during the entire phase of this outburst. No soft or soft-intermediate spectral states are observed. This could be due to the fact that this object belongs to a special class of sources (e.g., MAXI J1659-152, Swift J1753.5-0127, etc.) that have very short orbital periods and that the companion is profusely mass-losing or the disk is immersed inside an excretion disk. In these cases, flows in the accretion disk are primarily dominated by low viscous sub-Keplerian flow and the Keplerian rate is not high enough to initiate softer states. Low-frequency quasi-periodic oscillations (QPOs) are observed sporadically although as in normal outbursts of transient black holes, monotonic evolutions of QPO frequency during both rising and declining phases are observed. From the TCAF fits, we find the mass of the black hole in the range of 7.5–11 M{sub ⊙}, and from time differences between peaks of the Keplerian and sub-Keplerian accretion rates we obtain a viscous timescale for this particular outburst, ∼10 days.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Assessing the combined influence of fluid-rock interactions on reservoir properties and injectivity during CO 2 storage in saline aquifers

Complex fluid-rock interactions can occur during the injection of carbon dioxide (CO 2 ) into saline aquifers for sequestration, which may affect CO 2 injectivity and storage capacity. A comprehensive reactive transport model is established to analyze salt precipitation, CO 2 -water-rock geochemical reactions, and their effects on reservoir physical properties and injectivity. In addition, sensitivity analyses are conducted to investigate the main factors that affect fluid-rock interactions and injectivity with relevance for site selection for CO 2 storage. Results show that the back flow of formation water not only affects the salt precipitation but also affects the CO 2 -water-rock geochemical reactions, resulting in salt and calcite precipitations mainly occurring in the dry-out zone. However, most of the mineral dissolution/precipitation caused by CO 2 -water-rock reaction occurs in the two-phase and aqueous-phase zones, and their effect on reservoir porosity and permeability are small. A considerable amount of sodium chloride precipitates in the dry-out zone as brine is drawn by capillary action into this zone, with significant consequences for porosity, permeability and injectivity. The injection rate, salinity, capillary pressure–saturation relationships, and reservoir permeability strongly affect the distribution of salt precipitation. Finally, moderate injection rates, salinities, capillary pressures, and permeabilities all lead to favorable CO 2 injectivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Densification of non-radioactive porous siliceous particles loaded with cesium potassium copper hexacyanoferrate by spark plasma sintering

In the development of a conditioning process for cesium-bearing sorbents from the decontamination of aqueous waste streams, the treatment of a granular material composed of (Cs,K) 2 CuFe(CN) 6 loaded onto mesoporous silica particles has been investigated in an additive-free two-step conditioning process involving (i) calcination at 1000 °C to remove volatile species and (ii) densification by spark plasma sintering. A 92.0 % dense pellet was obtained after 10 min sintering at 850 °C under 100 MPa. Density and microstructural measurements indicate that densification is first driven by the viscous flow of a vitreous phase containing Cs, formed during calcination. Additionally, the density of the sample then increases further during sintering by grain boundary diffusion. However, a limited but non negligible part of cesium is lost during densification.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solution-phase sample-averaged single-particle spectroscopy of quantum emitters with femtosecond resolution

Here, the development of many optical quantum technologies depends on the availability of solid-state single quantum emitters with near-perfect optical coherence. However, a standing issue that limits systematic improvement is the significant sample heterogeneity and lack of mechanistic understanding of microscopic energy flow at the single emitter level and ultrafast timescales. Here we develop solution-phase single-particle pump-probe spectroscopy with photon correlation detection that captures sample-averaged dynamics in single molecules and/or defect states with unprecedented clarity at femtosecond resolution. We apply this technique to single quantum emitters in two-dimensional hexagonal boron nitride, which suffers from significant heterogeneity and low quantum efficiency. From millisecond to nanosecond timescales, the translation diffusion, metastable-state-related bunching shoulders, rotational dynamics, and antibunching features are disentangled by their distinct photon-correlation timescales, which collectively quantify the normalized two-photon emission quantum yield. Leveraging its femtosecond resolution, spectral selectivity and ultralow noise (two orders of magnitude improvement over solid-state methods), we visualize electron-phonon coupling in the time domain at the single defect level, and discover the acceleration of polaronic formation driven by multi-electron excitation. Corroborated with results from a theoretical polaron model, we show how this translates to sample-averaged photon fidelity characterization of cascaded emission efficiency and optical decoherence time. Our work provides a framework for ultrafast spectroscopy in single emitters, molecules, or defects prone to photoluminescence intermittency and heterogeneity, opening new avenues of extreme-scale characterization and synthetic improvements for quantum information applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The kinetic analog of the pressure–strain interaction

Energy transport in weakly collisional plasma systems is often studied with fluid models and diagnostics. However, the applicability of fluid models is limited when collisions are weak or absent, and using a fluid approach can obscure kinetic processes that provide key insights into the physics of energy transport. Kinetic diagnostics retain all of the information in 3D-3V phase space and thereby reach beyond the insights of fluid models to elucidate the mechanisms responsible for collisionless energy transport. In this work, we derive the Kinetic Pressure–Strain (KPS): a kinetic analog of the pressure–strain interaction, which is the channel between flow energy density and internal energy density in fluid models. Through two case studies of electron Landau damping, we demonstrate that the KPS diagnostic can elucidate kinetic mechanisms that are responsible for energy transport in this channel, just as the related field–particle correlation is known to identify kinetic mechanisms of transport between electromagnetic field energy density and kinetic energy density in particle flows. In addition, we show that resonant electrons play a major role in transferring energy between fluid flows and internal energy during the process of Landau damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Three-phase numerical model for subsurface hydrology in permafrost-affected regions (PFLOTRAN-ICE v1.0)

Degradation of near-surface permafrost due to changes in the climate is expected to impact the hydrological, ecological and biogeochemical responses of the Arctic tundra. From a hydrological perspective, it is important to understand the movement of the various phases of water (gas, liquid and ice) during the freezing and thawing of near-surface soils. We present a new non-isothermal, single-component (water), three-phase formulation that treats air as an inactive component. This single component model works well and produces similar results to a more complete and computationally demanding two-component (air, water) formulation, and is able to reproduce results of previously published laboratory experiments. A proof-of-concept implementation in the massively parallel subsurface flow and reactive transport code PFLOTRAN is summarized, and parallel performance of that implementation is demonstrated. When water vapor diffusion is considered, a large effect on soil moisture dynamics is seen, which is due to dependence of thermal conductivity on ice content. A large three-dimensional simulation (with around 6 million degrees of freedom) of seasonal freezing and thawing is also presented.

54 ENVIRONMENTAL SCIENCES↗