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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Biocatalytic Yarn for Peroxide Decomposition with Controlled Liquid Transport

In this study, a robust biocatalytic yarn with controllable liquid transport properties is created by coating thin layers of chitosan containing catalase onto a cellulosic yarn. The resulting material integrates enzyme catalytic functionality with protective coating properties of chitosan and structural functionality of the textile. Mild immobilization conditions and good affinity between the two polysaccharides minimize enzyme inactivation during the preparation steps and prevent enzyme from leaching during peroxide decomposition testing and washing, providing a novel and versatile enzyme immobilization strategy. The catalytic efficiency of enzymes in a reaction containing solid, liquid, and gas phases is facilitated when dissolved enzyme substrate is transported by liquid flowing through the coated textile structure. The flow-through configuration decomposes at least two times more peroxide in a twenty-times smaller reaction zone volume compared to a stirred tank configuration. Liquid transport through the yarn and liquid spatial distribution within the yarn are investigated by in situ neutron radiography and neutron computed tomography, revealing a constrained wicking mechanism that benefits biocatalytic yarn performance. This new class of sustainable and flexible biocatalytic textile matrices has beneficial multifunctional properties, not previously described, that are applicable for numerous small- and large-scale applications including controlled flow reactors and reactive filtration.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Online data-driven changepoint detection for high-dimensional dynamical systems

In this study, the detection of anomalies or transitions in complex dynamical systems is of critical importance to various applications. In this study, we propose the use of machine learning to detect changepoints for high-dimensional dynamical systems. Here, changepoints indicate instances in time when the underlying dynamical system has a fundamentally different characteristic—which may be due to a change in the model parameters or due to intermittent phenomena arising from the same model. We propose two complementary approaches to achieve this, with the first devised using arguments from probabilistic unsupervised learning and the latter devised using supervised deep learning. To accelerate the deployment of transition detection algorithms in high-dimensional dynamical systems, we introduce dimensionality reduction techniques. Our experiments demonstrate that transitions can be detected efficiently, in real-time, for the two-dimensional forced Kolmogorov flow and the Rössler dynamical system, which are characterized by anomalous regimes in phase space where dynamics are perturbed off the attractor at potentially uneven intervals. Finally, we also demonstrate how variations in the frequency of detected changepoints may be utilized to detect a significant modification to the underlying model parameters by utilizing the Lorenz-63 dynamical system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Nonequilibrium renormalization group fixed points of the quantum clock chain and the quantum Potts chain

In this work, we derive an exact renormalization group recursion relation for the Loschmidt amplitude of the quantum Q-state clock model and the quantum Q-state Potts model in one dimension. The renormalization group flow is discussed in detail. The fixed points of the renormalization group flow are found to be complex in general. These fixed points control the dynamical phases of the two models, giving rise to nonanalyticities in its Loschmidt rate function, for both the pure and the disordered system. For the quench protocols studied, dynamical quantum phase transitions are found to occur in the clock model for all Qs considered, while in the Potts model, they only occur when Q < 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Helium recovery system at IB3A

The growing demand for sustainable cryogenic operations at Fermilab has underscored the need to improve helium management, particularly at the Industrial Building 3A (IB3A) test facility. IB3A characterizes and tests superconductors, cables, and coils for projects such as the HL-LHC AUP and Mu2e, yet currently relies on 500 L Dewars whose boil-off is vented to atmosphere, wasting a critical, non-renewable resource and increasing the cost of testing. A project is therefore under way to link IB3A to an existing purification and liquefaction station in a neighboring building through a dedicated pipeline. Captured helium will be transferred, purified, reliquefied, and returned for reuse, cutting losses and operating costs. This paper details the first two project phases: "Design and Engineering" and "Procurement and Installation." The design phase finalizes pipeline specifications, establishes flow-control requirements, and resolves integration challenges with existing cryogenic infrastructure. The procurement and installation phase covers material sourcing, pipeline construction, and deployment of control and monitoring systems to assure reliable, efficient operation. Key technical hurdles-route optimization, pressure drop mitigation, and interface compatibility-are discussed alongside implemented solutions. Implementing the pipeline and upgrading IB3A will dramatically reduce helium consumption and therefore testing cost, strengthening Fermilab's capacity to support frontier science far into the future.

Porwisiak, D. [Fermilab; Wroclaw Tech. U.]↗
DOE OSTI2021

Methods and systems for detection of man-in-the-middle attacks for SCADA communication networks and applications of same

A system for detecting MITM for SCADA communication networks includes secure substation-substation communication links for providing secure and reliable paths to exchange OT data between substations for OT data consistency check; a SIB in each substation for sampling CT and PT measurements to calculate voltage magnitude and phase angle thereof; a S&C server in each substation coupled to the SIB for receiving the voltage magnitude and phase angle from the SIB and obtaining a packet carrying active power flow in transmission lines between two substations and a time stamp; an IDS server placed in a SCADA center for collecting the packet of each substation sent by the S&C server; analyzing the received packet from every adjacent substation; inspecting the payload of the received packet; and triggering an intrusion alarm to a SCADA operator when the power flow is not the same as the payload of the packets.

McCann, Roy A.↗
DOE OSTI2019

Multiscale dendritic needle network model of alloy solidification with fluid flow

We present a mathematical formulation of a multiscale model for solidification with convective flow in the liquid phase. The model is an extension of the dendritic needle network approach for crystal growth in a binary alloy. We propose a simple numerical implementation based on finite differences and step-wise approximations of parabolic dendritic branches of arbitrary orientation. Results of the two-dimensional model are verified against reference benchmark solutions for steady, unsteady, and buoyant flow, as well as steady-state dendritic growth in the diffusive regime. Simulations of equiaxed growth under forced flow yield dendrite tip velocities within 10% of quantitative phase-field results from the literature. Finally, we perform illustrative simulations of polycrystalline solidification using physical parameters for an aluminum-10wt% copper alloy. Resulting microstructures show notable differences when taking into account natural buoyancy in comparison to a purely diffusive transport regime. The resulting model opens new avenues for computationally and quantitatively investigating the influence of fluid flow and gravity-induced buoyancy upon the selection of dendritic microstructures. Further ongoing developments include an equivalent formulation for directional solidification conditions and the implementation of the model in three dimensions, which is critical for quantitative comparison to experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TU-H-CAMPUS-IeP2-03: Development of 3D Printed Coronary Phantoms for In-Vitro CT-FFR Validation Using Data from 320- Detector Row Coronary CT Angiography

Purpose: To validate Computed Tomography Fractional Flow Reserve (CT-FFR) measurements with accurate 3D printed coronary phantoms. Methods: DICOM data from four phases in two patients imaged with a standard 320 × 0.5mm coronary CT acquisition (70–80% cardiac cycle) underwent semi-automated segmentation using a research workstation. Both patients had a >50% stenosis from the clinical image interpretation. Each volume was saved as a Stereo Lithographic (STL) file with 250 micron resolution. The 3D geometries were qualitatively assessed; the best of the four phases was 3D printed using a Stratasys Eden260V printer in Tango+, a rubber-like material that roughly emulates mechanical properties of human vasculature. We connected the model to a programmable pump and measured the pressure drop using pressure sensors embedded proximal and distal to the arterial stenosis. Next, the STL files used for the 3D printed models were uploaded in the ANSYS meshing tool (ICEM CFD 16.1). A standard meshing process was applied and the meshed geometry was directly imported in the ANSYS Fluent for Computational Flow Dynamics simulations. The CFD simulations were used to calculate the CT-FFR and compared to the bench top FFR measured in the 3D printed phantoms. Results: FFR-CT measurements and phantoms were completed in within an hour after the segmentation. Patient 1 had a 60% stenosis that resulted in a CT-FFR of 0.68. The second case had a 50% stenosis and a CT-FFR of 0.75. The average bench top FFR measurements were 0.72 and 0.80, respectively. Conclusion: This pilot investigation demonstrated the use of a bench-top coronary model for CT-FFR validation. The measurements and the CFD simulations agreed within 6%. Project supported by Support: Toshiba America Medical Systems Corp.and NIH grant R01-EB002873. Project supported by Toshiba America Medical Systems Corp.and partial support from NIH grant R01-EB002873.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A giant Lyα nebula in the core of an x-ray cluster at z = 1.99: implications for early energy injection

We present the discovery of a giant ≳100 kpc Lyα nebula detected in the core of the X-ray emitting cluster CL J1449+0856 at z = 1.99 through Keck/LRIS narrow-band imaging. This detection extends the known relation between Lyα nebulae and overdense regions of the universe to the dense core of a 5–7 × 10{sup 13} M {sub ⊙} cluster. The most plausible candidates to power the nebula are two Chandra-detected AGN host cluster members, while cooling from the X-ray phase and cosmological cold flows are disfavored primarily because of the high Lyα to X-ray luminosity ratio (L{sub Lyα}/L{sub X}≈0.3, ≳10–1000 times higher than in local cool-core clusters) and by current modeling. Given the physical conditions of the Lyα-emitting gas and the possible interplay with the X-ray phase, we argue that the Lyα nebula would be short-lived (≲10 Myr) if not continuously replenished with cold gas at a rate of ≳1000 M {sub ⊙} yr{sup −1}. We investigate the possibility that cluster galaxies supply the required gas through outflows and we show that their total mass outflow rate matches the replenishment necessary to sustain the nebula. This scenario directly implies the extraction of energy from galaxies and its deposition in the surrounding intracluster medium (ICM), as required to explain the thermodynamic properties of local clusters. We estimate an energy injection of the order of ≈2 keV per particle in the ICM over a 2 Gyr interval. In our baseline calculation, AGNs provide up to 85% of the injected energy and two-thirds of the mass, while the rest is supplied by supernovae-driven winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Blade Designs For Improved Multi-Phase Performance In sCO 2 Compressors; Part II - Optical Diagnostics In sCO 2 And Experimental Evaluation With Particle Image Velocimetry

This paper presents the second part of a study in which the leading-edge and suction surface of a compressor blade was modified to delay onset of phase change for sCO 2 compressors operating near the critical point. Using a first-of-its-kind apparatus for the measurement of sCO 2 flow fields, Particle Image Velocimetry (PIV) is used for local flow field measurements of two compressor blade geometries: the modified “biased-wedge,” and a conventional constant thickness blade. Utilizing the developed hardware, the feasibility of a simple, laser-based diagnostic for qualitatively measuring liquid phase regions, is also presented. The design of the optical diagnostics rig, a discussion of numerous challenges, and necessary considerations involved in performing optical-based measurements like PIV, in sCO 2 , are discussed. Velocity field measurements for the modified compressor profile show a much lower suction peak compared to a conventional blade. Furthermore, these results validate numerical results at the tested conditions, where the suction side profile of the biased wedge works to minimize the local pressure gradient.

14 SOLAR ENERGY↗
DOE OSTI2021

Robust Distributed State Estimator for Interconnected Transmission and Distribution Networks (Final Report RPPR-1)

This project’s objective is to develop a combined transmission and distribution state estimator which accounts for very large system size and model complexity (by way of distributing the computations) and large number of solar PV units connected to the distribution system on multiple feeders. The project not only provides a robust formulation and solution to this problem but also tests the solution by implementing it on a well-established large utility system. It introduces several improvements with respect to the state of the art in existing state estimation software: (a) The developed state estimator (SE) allows robust and accurate monitoring of bidirectional flows in distribution systems which result due to the distributed energy sources which are not observable and thus not incorporated in generation dispatch; (b) Large utility systems with tens of thousands of transmission buses and hundreds of thousands of distribution nodes are difficult to model as a single integrated system. This shortcoming is addressed by developing a “scalable distributed computational framework” which allows splitting the ultra large system models into several small subsystems and coordinating their solution by a robust and practical state estimation formulation; (c) Measurement errors irrespective of their locations are detected and removed by the developed state estimator. Historically, transmission and distribution systems were analyzed and operated as two independent systems. Given the non-transposed short feeder sections, unevenly loaded phases, strictly radial configuration and unidirectional power flows in the absence of remote generation, distribution system analysis was customized to account for these characteristics. However, some of these assumptions are no longer valid (non-radial configuration, bidirectional power flows) and thus distribution system analysis should be revisited. Furthermore, in the past, the interaction between the transmission and distribution systems was quite passive, where distribution substations were modeled as lumped loads in the transmission system model. With substantial generation injected by renewable generation located in the distribution systems, such modeling will no longer be accurate. The developed state estimator facilitates proper monitoring of the interactions between the transmission and distribution systems and enables smart dispatch of these units which are made observable by the state estimator.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Characterization of the external and internal flow structure of an aerated-liquid injector using X-ray radiography and fluorescence

In the present study, the internal flowfield of aerated-liquid fuel injectors is examined through x-ray radiography and x-ray fluorescence. An inside-out injector, consisting of a perforated aerating tube within an annular liquid stream, sprays into a quiescent environment at a fixed mass flow rate of water and nitrogen gas. The liquid is doped with bromine (in the form of NaBr) to create an x-ray fluorescence signal. This allows for reasonable absorption and fluorescence signals, and one or both diagnostics can be used to track the liquid distribution. The injector housing is fabricated from beryllium (Be), which allows the internal flowfield to be examined (as Be has relatively low x-ray attenuation coefficient). Two injector geometries are compared, illustrating the effects of aerating orifice size and location on the flow evolution. Time-averaged equivalent pathlength (EPL) and line-of-sight averaged density ρ(y) reveal the formation of the two-phase mixture, showing that the liquid film thickness along the injector walls is a function of the aerating tube geometry, though only upstream of the nozzle. These differences in gas and liquid distribution (between injectors with different aerating tube designs) are suppressed as the mixture traverses the nozzle contraction. The averaged liquid velocity (computed from the density and liquid mass flow rate) reveal a similar trend. This suggests that at least for the current configurations, the plume width, liquid mass distribution, and averaged liquid velocity for the time-averaged external spray are insensitive to the aerating tube geometry.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows

Here, the present paper addresses important fundamental issues of inter-phase heat transfer and energy coupling in turbulent dispersed multiphase flows through scaling analysis. In typical point-particle or two-fluid approaches, the fluid motion and convective heat transfer at the particle scale are not resolved and the momentum and energy coupling between fluid and particles are provided by proper closure models. By examining the kinetic energy transfer due to the coupling forces from the macroscale to microscale fluid motion, closure models are obtained for the contributions of the coupling forces to the energy coupling. Due to the inviscid origin of the added-mass force, its contribution to the microscale kinetic energy does not contribute to dissipative transfer to fluid internal energy as was done by the quasi-steady force. Time scale analysis shows that when the particle is larger than a critical diameter, the diffusive-unsteady kernel decays at a time scale that is smaller than the Kolmogorov time scale. As a result, the computationally costly Basset-like integral form of diffusive-unsteady heat transfer can be simplified to a non-integral form. Conventionally, the fluid-to-particle volumetric heat capacity ratio is used to evaluate the relative importance of the unsteady heat transfer to the energy balance of the particles. Therefore, for gas-particle flows, where the fluid-to-particle volumetric heat capacity ratio is small, unsteady heat transfer is usually ignored. However, the present scaling analysis shows that for small fluid-to-particle volumetric heat capacity ratio, the importance of the unsteady heat transfer actually depends on the ratio between the particle size and the Kolmogorov scale. Furthermore, the particle mass loading multiplied by the heat capacity ratio is usually used to estimate the importance of the thermal two-way coupling effect. Through scaling argument, improved estimates are established for the energy coupling parameters of each energy exchange mechanism between the phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Three‐Phase CO 2 Flow in a Basalt Fracture Network

Abstract Geologic CO 2 sequestration in basalt reservoirs is predicated on permanent CO 2 isolation via rapid mineralization reactions. This process is supported by a substantial body of evidence, including laboratory experiments documenting rapid mineralization rates, regional storage estimates indicating large, accessible storage reservoirs, and two successful pilot‐scale studies. Nevertheless, there remains significant uncertainty in the behavior of CO 2 flow within basalt fracture networks, particularly in the context estimating physical trapping potential in early time and as CO 2 undergoes phase change. In this study, a Monte Carlo numerical model is designed to simulate a supercritical CO 2 plume infiltrating a low‐permeability flood basalt entablature. The fracture network model is based on outcrop‐scale LiDAR mapping of Columbia River Basalt, and CO 2 flow is simulated within fifty equally probable realizations of the fracture network. The spatial distribution of fracture permeability for each realization is randomly drawn from a basalt aperture distribution, and ensemble results are analyzed with e‐type estimates to compute mean and standard deviation of fluid pressure and CO 2 saturation. Results of this model after 10 years of simulation suggest that (1) CO 2 flow converges on a single dominant flow path, (2) CO 2 accumulates at fracture intersections, and (3) variability in permeability can account for a 1.6 m depth interval within which free CO 2 may change phase from supercritical fluid to subcritical liquid or gas. In the context of CO 2 sequestration in basalt, these results suggest that physical CO 2 trapping may be substantially enhanced as carbonate minerals precipitate within the basalt fracture network.

Gierzynski, Alec O.↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2020

Investigation of countercurrent flow profile and liquid holdup in random packed column with local CFD data

Liquid holdup and mass transfer area are critical parameters used in packed column design and CO 2 capture efficiency prediction. Compared to experimental methods, 3D computational fluid dynamics (CFD) can provide high-resolution data for studying countercurrent flow hydrodynamics in a packed column. In this paper, a framework was established for modeling the liquid–gas countercurrent flow in a random packed column with pall rings. Besides the column-averaged information, the radial pall ring distribution, velocity, and liquid holdup profiles are obtained in the simulations. The radial two-phase velocity and liquid holdup profiles show significant fluctuations caused by the non-uniform pall ring distribution. The entrance effect and wall influence are also observed in the packed column. By using local countercurrent flow information, the packing specific area $a_p$ and liquid velocity $u_L$ range are significantly expanded for liquid holdup hL prediction. This enables us to develop the liquid holdup correlation with one packing size and six CFD runs only. The proposed power law of the liquid holdup correlation $h_L~u_L^{0.44}$ indicates the random packed column falls in a viscous to turbulent transition regime with a Reynolds Number range of [6.7–40.2]. The derived liquid holdup correlation is in good agreement with existing correlations that are developed using the column-averaged experimental data.

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

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model: Improvement and Quality Assurance

The multisegmented wellbore model (MSW) semi-analytically estimates the amount of CO 2 and brine leakage from a leaking legacy well by segmenting it into intervals to simulate site-specific stratigraphic and hydrogeologic properties. The model is a component of the National Risk Assessment Partnership Open-Source Integrated Assessment Model (NRAP-Open-IAM), which was developed to perform risk assessment for geologic CO 2 storage. The new wellbore leakage model, which uses deep learning networks for a caprock segment, was developed to enhance the analytical MSW. The model was trained and validated using a synthetic data set of Subsurface Transport Over Multiple Phases (STOMP) multiphase flow simulations from various geological, well attribute, and operational conditions to ensure its quality. The results demonstrate that the model is more accurate than the existing model in predicting the transport of two-phase fluids (brine and injected CO 2 ) through the well. This report provides a detailed explanation of the model development and quality assurance.

58 GEOSCIENCES↗
DOE OSTI2020

Probe Chiral Magnetic Effect with Signed Balance Function

In this paper a pair of observables are proposed as alternative ways, by examining the fluctuation of net momentum-ordering of charged pairs, to study the charge separation induced by the Chiral Magnetic Effect (CME) in relativistic heavy ion collisions. They are, the out-of-plane to in-plane ratio of fluctuation of the difference between signed balance functions measured in pair’s rest frame, and the ratio of it to similar measurement made in the laboratory frame. In this work, both observables have been studied with simulations including flow-related backgrounds, and for the first time, backgrounds that are related to resonance's global spin alignment. The two observables have similar positive responses to signal, and opposite, limited responses to identifiable backgrounds arising from resonance flow and spin alignment. Both observables have also been tested with two realistic models, namely, a multi-phase transport (AMPT) model and the anomalous-viscous fluid dynamics (AVFD) model. These two observables, when cross examined, will provide useful insights in the study of CME-induced charge separation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Symmetric Mass Generation with Four SU(2) Doublet Fermions

We study a single exactly massless staggered fermion in the fundamental representation of an SU(2) gauge group. We utilize an nHYP-smeared fermion action supplemented with additional heavy Pauli-Villars fields that serve to decrease lattice artifacts. The phase diagram exhibits a clear two-phase structure with a conformal phase at weak coupling and a novel new phase, the symmetric mass generation (SMG) phase, appearing at strong coupling. The SMG phase is confining with all states gapped and chiral symmetry unbroken. Our finite size scaling analysis provides strong evidence that the phase transition between these two phases is continuous, which would allow for the existence of a continuum SMG phase. Furthermore, the renormalization group flows are consistent with a 𝛽 function that vanishes quadratically at the new fixed point, suggesting that the 𝑁 𝑓 = 4 flavor SU(2) gauge theory lies at the opening of the conformal window.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗