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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

CTF Void Drift Validation Study

This milestone report is a summary of work performed in support of expansion of the validation and verification (V&V) matrix for the thermal-hydraulic subchannel code, CTF. The focus of this study is on validating the void drift modeling capabilities of CTF and verifying the supporting models that impact the void drift phenomenon. CTF uses a simple turbulent-diffusion approximation to model lateral cross-flow due to turbulent mixing and void drift. The void drift component of the model is based on the Lahey and Moody model. The models are a function of two-phase mass, momentum, and energy distribution in the system; therefore, it is necessary to correctly model the ow distribution in rod bundle geometry as a first step to correctly calculating the void distribution due to void drift.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Pulsar Wind Bubble Blowout from a Supernova

For pulsars born in supernovae, the expansion of the shocked pulsar wind nebula is initially in the freely expanding ejecta of the supernova. While the nebula is in the inner flat part of the ejecta density profile, the swept-up, accelerating shell is subject to the Rayleigh–Taylor instability. We carried out two- and three-dimensional simulations showing that the instability gives rise to filamentary structure during this initial phase but does not greatly change the dynamics of the expanding shell. The flow is effectively self-similar. If the shell is powered into the outer steep part of the density profile, the shell is subject to a robust Rayleigh–Taylor instability in which the shell is fragmented and the shocked pulsar wind breaks out through the shell. The flow is not self-similar in this phase. For a wind nebula to reach this phase requires that the deposited pulsar energy be greater than the supernova energy, or that the initial pulsar period be in the ms range for a typical 10{sup 51} erg supernova. These conditions are satisfied by some magnetar models for Type I superluminous supernovae. We also consider the Crab Nebula, which may be associated with a low energy supernova for which this scenario applies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Experiments and CFD simulation of mass transfer and hydrodynamics in a cylindrical bubble column

Flow visualization, Particle Image Velocimetry (PIV) and dissolved oxygen charging experiments are carried out in a small-scale cylindrical air-water bubble column at different superficial gas velocities. Void fraction profile, total gas holdup, liquid phase velocity field and volumetric mass transfer coefficient are obtained. Details of each experiment are given, with a new procedure for obtaining a 2D void fraction profile from flow visualization by image processing and time averaging which could provide useful information for CFD model validation. The experimental results are analyzed and compared with simulations where good agreement is found for both local and global flow patterns. The averaged governing equations and constitutive relations used in the simulation are discussed briefly, with emphasis on the proper comparison between simulations and experiments based on the derivation of the two-fluid model. Finally, by conducting three sets of experiments related to phase distribution, velocity field and mass transfer in the same apparatus, a validated CFD model for multiphase mass transfer is established and can be used in the future design of a multi-phase reactor.

42 ENGINEERING↗
DOE OSTI2020

Conservation Voltage Reduction (CVR) via Two-Timescale Control in Unbalanced Power Distribution Systems

Voltage control devices are employed in power distribution systems to reduce the power consumption by operating the system closer to the lower acceptable voltage limits; this technique is called conservation voltage reduction (CVR). The different modes of operation for system’s legacy devices (with binary control) and new devices (e.g. smart inverters with continuous control) coupled with variable photovoltaic (PV) generation results in voltage fluctuations which makes it challenging to achieve CVR objective. This paper presents a two-timescale control of feeder’s voltage control devices to achieve CVR that includes (1)a centralized controller operating in a slower time-scale to coordinate voltage control devices across the feeder and (2) local controllers operating in a faster timescale to mitigate voltage fluctuations due to PV variability. The centralized controller utilizes a three-phase optimal power flow model to obtain the decision variables for both legacy devices and smart inverters. The local controllers operate smart inverters to minimize voltage fluctuations and restore nodal voltages to their reference values by adjusting the reactive power support. The proposed approach is validated using the IEEE-123 bus (medium-size) and R3-12.47-2 (large-size) feeders. It is demonstrated that the proposed approach is effective in achieving the CVR objective for unbalanced distribution systems.

optimal control, distributed generation, fluctuati↗
DOE OSTI2022

Vapor–liquid equilibrium estimation of n-alkane/nitrogen mixtures using neural networks

Understanding fluid phase behavior, like VLE, in high P&T conditions is crucial for developing high-fidelity simulations of chemically reacting flows in liquid-fueled combustion systems and also forms an integral part of the design-modeling of the control processes in chemical industries. Two data-driven models have been proposed here in this study, each of which was competent in estimating VLE for the Type III binary systems of C 10 /N 2 and C 12 /N 2 , at pressures ranging up to 50–60 MPa. Both models showed better performance in predicting equilibrium pressure as compared to VLE modeled using PR-EOS. A modified model has also been proposed, capable of estimating the full phase envelope for the binary systems of C 10 /N 2 and C 12 /N 2 across a wide range of temperatures, and thus exhibit the mixture critical pressure at the concerned temperature. The diverse applicability of the proposed network architecture was further exhibited while estimating the VLE of a ternary system of C 1 /C 10 /N 2 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Thermal Gradient Effects on Redox Evolution and Volatility-Driven Fractionation in Ternary U/Ce/Cs Condensates

Understanding how thermal history influences redox evolution and chemical fractionation is essential for characterizing high-temperature condensation in complex materials, including nuclear debris. Here, we tested the hypothesis that distinct thermal regimes in a plasma flow reactor influence redox pathways and elemental partitioning in ternary U/Ce/Cs systems. A configurable plasma flow reactor was modified with an external tube furnace to impose two distinct thermal gradients: continuous ambient cooling and a furnace-assisted thermal hold-up near 1400 K followed by rapid cooling. Transmission electron microscopy characterized phase identity, morphology, and nanoscale element distributions, while inductively coupled plasma-mass spectrometry quantified bulk elemental ratios. Across both thermal regimes, uranium and cerium condensed as UO 2 and CeO 2 as dominant refractory oxide products. Uranium partially oxidized to α-UO 3 during extended ambient cooling, while furnace-assisted hold-up preserved UO 2 and produced partial reduction of cerium to Ce 2 O 3 . Cesium remained volatile upstream and condensed later in the reactor, forming Cs 2 O and Cs-uranate phases with the highest incorporation after thermal hold-up. Bulk ICP-MS measurements supported these observations. U/Ce ratios remained comparatively stable and Cs displayed delayed and apparent transient enrichment that matched the nanoscale measurements. This integrated approach provides a quantitative method for linking thermal gradients to redox evolution and volatility-driven fractionation. These results show how the plasma flow reactor can identify where equilibrium descriptions remain adequate and where kinetic effects from residence time and temperature history must be considered when interpreting condensation behavior in multicomponent systems.

and nuclear chemistry↗
DOE OSTI2020

Mechanistic modeling of CO 2 well leakage in a generic abandoned well through a bridge plug cement-casing gap

Both known and unmapped plugged and abandoned wells are potential leakage pathways for CO 2 from geologic carbon sequestration (GCS) sites. Although many abandoned wells have cement bridge plugs installed to prevent leakage, the seal between the cement and the inner casing wall is subject to failure. In this study, we carried out detailed T2Well simulations of cases of sudden non-Darcy flow of CO 2 and brine leakage up the gap between a cement plug and the inner steel casing wall that becomes a fully connected flow path during the post-injection period. The goal of our study was two-fold: (1) to understand the dynamics, rates, and the characteristic temporal signals associated with the onset of leakage through various gap-aperture sizes, and (2) to suggest potential monitoring strategies based on the findings. Simulation results show that the leakage of CO 2 and brine upward is transient with interesting phase interference behavior. Time-dependent oscillatory flows with varying pressure, temperature, and flow rates of CO 2 and brine show strong dependence on gap aperture. Phase-change and decompression lead to very low temperatures at the top of the well for gap apertures larger than 4 mm suggesting that remote thermal monitoring at the ground surface may be an effective way of monitoring even if well locations are not known a priori. Pressure in the well is also indicative of CO 2 leakage. The temporal patterns of changing temperature and pressure may be useful diagnostic signals for leakage detection. Finally, these transient leakage signals may provide information on the cause of leakage and/or characteristics of the flow path that could inform effective remediation design and execution approaches.

58 GEOSCIENCES↗
DOE OSTI2024

The role of RF-induced E×B flows in the mitigation of scrape-off-layer convective transport during ion cyclotron resonance heating

Abstract While multiple experiments have reported a decrease in intermittent fluctuations in the far Scrape-Off-Layer (SOL) during ion cyclotron resonance heating (Antar et al 2010 Phys. Rev. Lett. 105 165001, Li et al 2022 Nucl. Eng. Technol. 54 207–19, Antar et al 2012 Nucl. Fusion 52 103005), the physical mechanism behind this observation has not been fully established yet. In this work, we demonstrate, for the first time, a direct correlation between the amplitude of RF-induced E×B flows and turbulence suppression in the far SOL. Using the Gas Puff Imaging (GPI) diagnostic on Alcator C-Mod, we show again that Ion Cyclotron Range of Frequencies can significantly alter the flow in the SOL and introduce a shear layer in regions magnetically connected close to the antenna (Cziegler et al 2012 Plasma Phys. Control. Fusion 54 105019). With the 4-strap field-aligned antenna operated in dipole phasing, the ratio of the power coupled by the central two straps to the power coupled by the outer two straps was varied. The resulting RF-induced radial electric field magnitude thus varied, and we show that the impact on the far SOL turbulence correlates with the modified E×B velocity. We then apply a newly-developed blob tracking algorithm (Han et al 2022 Sci. Rep. 12 18142) to higher-resolution GPI videos in order to directly observe the process of blob shearing by RF-induced E×B flows. We show that the radially sheared poloidal flows act as a transport barrier by stretching, stopping, and destroying filaments, which is consistent with the observed difference in turbulence statistics.

Diab, R. (ORCID:0000000342065868)↗
DOE OSTI2015

A CIRCULAR-RIBBON SOLAR FLARE FOLLOWING AN ASYMMETRIC FILAMENT ERUPTION

The dynamic properties of flare ribbons and the often associated filament eruptions can provide crucial information on the flaring coronal magnetic field. This Letter analyzes the GOES-class X1.0 flare on 2014 March 29 (SOL2014-03-29T17:48), in which we found an asymmetric eruption of a sigmoidal filament and an ensuing circular flare ribbon. Initially both EUV images and a preflare nonlinear force-free field model show that the filament is embedded in magnetic fields with a fan-spine-like structure. In the first phase, which is defined by a weak but still increasing X-ray emission, the western portion of the sigmoidal filament arches upward and then remains quasi-static for about five minutes. The western fan-like and the outer spine-like fields display an ascending motion, and several associated ribbons begin to brighten. Also found is a bright EUV flow that streams down along the eastern fan-like field. In the second phase that includes the main peak of hard X-ray (HXR) emission, the filament erupts, leaving behind two major HXR sources formed around its central dip portion and a circular ribbon brightened sequentially. The expanding western fan-like field interacts intensively with the outer spine-like field, as clearly seen in running difference EUV images. We discuss these observations in favor of a scenario where the asymmetric eruption of the sigmoidal filament is initiated due to an MHD instability and further facilitated by reconnection at a quasi-null in corona; the latter is in turn enhanced by the filament eruption and subsequently produces the circular flare ribbon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

In-silico quantification of intersubject variability on aerosol deposition in the oral airway

The extrathoracic oral airway I s not only a major mechanical barrier for pharmaceutical aerosols to reach the lung but also a major source of variability in lung deposition. Using computational fluid dynamics, deposition of 1-30 µm particles was predicted in 11 CT-based models of the oral airways of adults (7 healthy subjects and 4 patients with chronic obstructive pulmonary disease). Simulations were performed for mouth breathing during both inspiration and expiration at two steady-state flow rates representative of resting/light activity (18 L/min) and of that typically used during aerosol inhalation from a dry powder inhaler (45 L/min). For each subject, oral deposition increased with increasing particle size and flow rate. When averaged over all subjects, there was no significant difference in oral deposition efficiency between inspiration and expiration, unlike subregional deposition that show significantly different patterns between the two breathing phases. Consistent with previous in-vitro inspiratory studies, our data showed a large intersubject variability in oral deposition. These results highlight the need for incorporating morphological variation of the upper airway in predictive models of aerosol deposition for accurate predictions of particle dosimetry in the intrathoracic region of the lung.

Computational fluid dynamics (CFD), Inertial impac↗
DOE OSTI2022

Large Magnetic Gap in a Designer Ferromagnet–Topological Insulator–Ferromagnet Heterostructure

Abstract Combining magnetism and nontrivial band topology gives rise to quantum anomalous Hall (QAH) insulators and exotic quantum phases such as the QAH effect where current flows without dissipation along quantized edge states. Inducing magnetic order in topological insulators via proximity to a magnetic material offers a promising pathway toward achieving the QAH effect at a high temperature for lossless transport applications. One promising architecture involves a sandwich structure comprising two single‐septuple layers (1SL) of MnBi 2 Te 4 (a 2D ferromagnetic insulator) with ultrathin few quintuple layer (QL) Bi 2 Te 3 in the middle, and it is predicted to yield a robust QAH insulator phase with a large bandgap greater than 50 meV. Here, the growth of a 1SL MnBi 2 Te 4 /4QL Bi 2 Te 3 /1SL MnBi 2 Te 4 heterostructure via molecular beam epitaxy is demonstrated and the electronic structure probed using angle‐resolved photoelectron spectroscopy. Strong hexagonally warped massive Dirac fermions and a bandgap of 75 ± 15 meV are observed. The magnetic origin of the gap is confirmed by the observation of the exchange‐Rashba effect, as well as the vanishing bandgap above the Curie temperature, in agreement with density functional theory calculations. These findings provide insights into magnetic proximity effects in topological insulators and reveal a promising platform for realizing the QAH effect at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2016

NUMERICAL SIMULATIONS OF KELVIN–HELMHOLTZ INSTABILITY: A TWO-DIMENSIONAL PARAMETRIC STUDY

Using two-dimensional simulations, we numerically explore the dependences of Kelvin–Helmholtz (KH) instability upon various physical parameters, including viscosity, the width of the sheared layer, flow speed, and magnetic field strength. In most cases, a multi-vortex phase exists between the initial growth phase and the final single-vortex phase. The parametric study shows that the evolutionary properties, such as phase duration and vortex dynamics, are generally sensitive to these parameters, except in certain regimes. An interesting result is that for supersonic flows, the phase durations and saturation of velocity growth approach constant values asymptotically as the sonic Mach number increases. We confirm that the linear coupling between magnetic field and KH modes is negligible if the magnetic field is weak enough. The morphological behavior suggests that the multi-vortex coalescence might be driven by the underlying wave–wave interaction. Based on these results, we present a preliminary discussion of several events observed in the solar corona. The numerical models need to be further improved to perform a practical diagnostic of the coronal plasma properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗
DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Comparison of URANS and LES predictions for the open phase of the OECD NEA CSNI fluid structure interaction CFD benchmark

The OECD NEA CSNI WGAMA CFD Task Group ran a benchmark in 2020 and 2021 to assess the predictive capabilities of coupled fluid structure interaction (FSI) CFD analysis methods. This paper presents the predictions made for the open phase of the benchmark using URANS and LES turbulence modelling approaches, and a comparison of the results to the experimental data. The benchmark comprised a channel containing two inline cylinders in cross-flow. The cylinders were fixed at one end, free at the other, and had measured resonant frequencies and damping properties. The URANS modelling used ANSYS Fluent 2-way coupled to ANSYS Mechanical. The LES modelling used Nek5000, 1-way coupled to Diablo. Comparisons with cross-channel velocity profiles are presented, both for the mean flow and its RMS. Comparisons are also made to the frequency spectra for point measurements of fluid velocity and pressure, and for the accelerations of the free end of each cylinder. URANS predicts the average velocity profiles relatively well, and is able to predict the velocity and acceleration spectra at the shedding frequency. However, the frequency content at the 4th harmonic of the shedding frequency is low in the URANS flow fields, and so does not excite accelerations at the resonant frequency of the cylinders. LES makes better predictions of the average profiles, and the velocity spectra agree well at both the shedding frequency and at higher frequencies. In conclusion, the 1-way coupled LES results show good agreement for acceleration spectra.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Combining Distributed Electrochemical Measurements with a Semi-Empirical Model to Identify Local Variations in Overpotentials for PEM Water Electrolyzers

Abstract A segmented cell and PCB approach was used to demonstrate an advanced technique for measuring and understanding the electrochemical behavior of PEM water electrolyzers. By employing a semi-empirical model, discrete modes of overpotential and their spatial distributions were identified in an operating cell. Distributed area-specific resistance (DASR measurements indicate that local dehydration of the polymer electrolyte membrane is the primary driver of increased polarization, behaving similarly to mass transport-limited currents at low flow rates due to ineffective gas-phase bubble removal. To improve quantification, the overall ohmic overpotential is divided into background (η Ohm,b ) and mass transport-coupled (η Ohm,mt ) components. Two types of anode diffusion media, a baseline PTL and a thin foil liquid/gas diffusion layer (LGDL), were examined. The results show significant variations in η Ohm,b and highlight that ineffective bubble removal leads to pronounced local dehydration. The effect of cell temperature on various overpotentials was also examined. Negligible differences in diffusion overpotential were observed when comparing PEWE performance using the PTL and LGDL, despite major differences in two-phase fluid transport within each diffusion media. This underscores the critical role of the catalyst layer's micro-porous structure in diffusion processes. The study highlights the complex nature and intricacies of mass transport limitations within PEWEs.

Roy, Anirban (ORCID:0000000163803386)↗
DOE OSTI2019

Small-scale flow topologies in decaying isotropic turbulence laden with finite-size droplets

The topology of the fine-scale motions in decaying isotropic turbulence laden with droplets of super-Kolmogorov size is investigated using results from direct numerical simulations. The invariants of the velocity-gradient, rate-of-strain, and rate-of-rotation tensors are computed in the carrier phase. The joint probability density functions of the invariants are calculated and conditioned on different distances from the droplet surface. The results show that outside the viscous region near the interface, the flow topologies favor stable focus/stretching and unstable node/saddle/saddle structures, which is in agreement with those found in canonical homogeneous isotropic turbulence. Inside the viscous layer at the droplet surface, the flow topologies shift from a preference for high-enstrophy/low-dissipation motions to favoring low-enstrophy/high-dissipation. At the droplet surface, there is a strong tendency for boundary-layer-like and vortex-sheet flow topologies in which the strain and rotation rates are positively correlated. Furthermore, an interesting observation is that the shapes of the invariant distributions at the droplet surface are remarkably similar to those reported in the viscous sublayer of turbulent wall-bounded flows. Also, the results show that the smallest hydrodynamic length scale of the carrier fluid turbulence is located at the droplet interface and that this length scale is two to three times smaller than that of the surrounding bulk flow. From a computational viewpoint, this suggests a more stringent spatial resolution requirement for the direct numerical simulation of finite-size droplets in isotropic turbulence than its single-phase counterpart.

42 ENGINEERING↗
DOE OSTI2021

Characterizing the Effect of Capillary Heterogeneity on Multiphase Flow Pulsation in an Intermediate–Scale Beadpack Experiment using Time Series Clustering and Frequency Analysis

An intermediate-scale beadpack drainage experiment was conducted to investigate how simple layered lamination heterogeneity affects CO 2 flow. Two simple layers of capillary barriers are manually packed in the tank and slow drainage was carried out using analog fluids to mimic the capillary- and gravity-dominated CO 2 upward migration process in deep saline aquifers. Nonwetting phase saturation time series clustering analysis and frequency analysis have been conducted on the experimental data. Additionally, modified invasion percolation numerical simulations were done on a digital model of the beadpack to compare to experimental results. Results show that capillary barriers can lead to strong pulsation behavior, which in turn can cause unexpected early breaching through other barriers. The inlet pressure is found to be able to respond to saturation changes in far regions of the domain, indicating that the wetting phase can transmit pressure changes from the other phase. Although static simulations were not able to capture all the dynamic behavior observed in the experiment, Monte Carlo composite simulation results combining many different realizations can better illustrate how the nonwetting phase will behave in the heterogeneous domain. Furthermore, our results suggest the need for CO 2 storage site selection with preference given to aquifers with more capillary barriers with finer grain sizes to avoid flow pulsation and to retard plume upward migration.

58 GEOSCIENCES↗