SEARCH · Search NASA

Results for “Multiphase”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Direct, efficient and selective capture of low concentration of CO 2 from natural gas flue gas using a high temperature tubular carbon capture membrane

Natural gas (NG) fired power plants emit low concentration (4–5%) of CO 2 , which presents additional technical and economic challenges to the current benchmark amine absorption technology. The newly emerged high-temperature multiphase membranes operated on molten carbonate (MC) chemistry for CO 2 capture/separation/conversion have been demonstrated with great potential to meet this challenge. In this study, we report on the CO 2 capture performance of such a membrane in tubular geometry from a mockup NG flue gas. The membrane is comprised of a mixture of Gd 0.20 Ce 0.80 O 1.95 (GDC) and MC, in which GDC forms a porous skeleton to contain MC. Here, we show that the membrane with a dimension of 6.1 mm in outer diameter, 5.1 mm in inner diameter and 5 cm in effective length (resulting in 4cm 2 effective surface area) can achieve a CO 2 flux density of 0.46–0.55 mL/min·cm 2 at 650°C, capturing 97% pure CO 2 at a rate of 37–42% from 5%CO 2 –N 2 using moistened Ar as the sweep gas. The level of performance demonstrated by this study suites the membrane well for stationary CO 2 capture from NG power plants.

03 NATURAL GAS↗

DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

A novel microstructural approach for TRIP/TWIP β-Ti alloys using intermetallic precipitation

Interest in metastable β-Ti has grown, due to the exceptional work hardening caused by TRansformation and TWinning Induced Plasticity (TRIP/TWIP). A current challenge is increasing the yield strength, while maintaining the benefits of TRIP/TWIP. This study presents the first report of TRIP/TWIP in a metastable β-Ti alloy containing intermetallic compounds, specifically Ti silicide precipitates. Silicon (Si) additions enable the formation of silicide particles, which can be precipitated, dissolved, and controlled through tailored heat treatments. In-situ synchrotron x-ray diffraction during tensile testing reveals primary {332}<113> twinning activates first, followed by TRIP, leading to extensive work hardening and uniform elongation in the solutionized state. In conclusion, these findings open a novel microstructural approach to the design of multiphase TRIP/TWIP titanium alloys, relying on intermetallic precipitation rather than allotropes such as α or ω phase, and lay the groundwork for future studies on intermetallic precipitation strengthening.

Materials science↗

DOE OSTI2017

Development of an integrated four-channel fast avalanche-photodiode detector system with nanosecond time resolution

In this paper, a four-channel nanosecond time-resolved avalanche-photodiode (APD) detector system is developed at Beijing Synchrotron Radiation. It uses a single module for signal processing and readout. This integrated system provides better reliability and flexibility for custom improvement. The detector system consists of three parts: (i) four APD sensors, (ii) four fast preamplifiers and (iii) a time-digital-converter (TDC) readout electronics. The C30703FH silicon APD chips fabricated by Excelitas are used as the sensors of the detectors. It has an effective light-sensitive area of 10 × 10 mm 2 and an absorption layer thickness of 110 μ m. A fast preamplifier with a gain of 59 dB and bandwidth of 2 GHz is designed to readout of the weak signal from the C30703FH APD. The TDC is realized by a Spartan-6 field-programmable-gate-array (FPGA) with multiphase method in a resolution of 1ns. The arrival time of all scattering events between two start triggers can be recorded by the TDC. The detector has been used for nuclear resonant scattering study at both Advanced Photon Source and also at Beijing Synchrotron Radiation Facility. Finally, for the X-ray energy of 14.4 keV, the time resolution, the full width of half maximum (FWHM) of the detector (APD sensor + fast amplifier) is 0.86 ns, and the whole detector system (APD sensors + fast amplifiers + TDC readout electronics) achieves a time resolution of 1.4 ns.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2017

Interface tracking simulations of bubbly flows in PWR relevant geometries

[en] Highlights: • Simulations were performed for turbulent bubbly flows in PWR subchannel geometry. • Liquid turbulence is fully resolved by direct numerical simulation approach. • Bubble behavior is captured using level-set interface tracking method. • Time-averaged single- and two-phase turbulent flow statistical quantities are obtained. - Abstract: The advances in high performance computing (HPC) have allowed direct numerical simulation (DNS) approach coupled with interface tracking methods (ITM) to perform high fidelity simulations of turbulent bubbly flows in various complex geometries. In this work, we have chosen the geometry of the pressurized water reactor (PWR) core subchannel to perform a set of interface tracking simulations (ITS) with fully resolved liquid turbulence. The presented research utilizes a massively parallel finite-element based code, PHASTA, for the subchannel geometry simulations of bubbly flow turbulence. The main objective for this research is to demonstrate the ITS capabilities in gaining new insight into bubble/turbulence interactions and assisting the development of improved closure laws for multiphase computational fluid dynamics (M-CFD). Both single- and two-phase turbulent flows were studied within a single PWR subchannel. The analysis of numerical results includes the mean gas and liquid velocity profiles, void fraction distribution and turbulent kinetic energy profiles. Two sets of flow rates and bubble sizes were used in the simulations. The chosen flow rates corresponded to the Reynolds numbers of 29,079 and 80,775 based on channel hydraulic diameter (D_h) and mean velocity. The finite element unstructured grids utilized for these simulations include 53.8 million and 1.11 billion elements, respectively. This has allowed to fully resolve all the turbulence scales and the deformable interfaces of individual bubbles. For the two-phase flow simulations, a 1% bubble volume fraction was used which resulted in 17 bubbles in the smaller case and 262 bubbles in the larger case. In the larger simulation case the size of the resolved bubbles is 0.65 mm in diameter, and the bulk mesh cell size is about 30 microns. Those large-scale simulations provide new level of details previously unavailable and were enabled by the excellent scaling performance of our two-phase flow solver and access to the state-of-the-art supercomputing resources. The presented simulations used up to 256 thousand processing threads on the IBM BG/Q supercomputer “Mira” (Argonne National Laboratory).

Nuclear Science & Technology↗

DOE OSTI2017

Wall-resolved spectral cascade-transport turbulence model

A spectral cascade-transport model has been developed and applied to turbulent channel flows (Reτ= 550, 950, and 2000 based on friction velocity, uτ ; or ReδΜ= 8,500; 14,800 and 31,000, based on the mean velocity and channel half-width). This model is an extension of a spectral model previously developed for homogeneous single and two-phase decay of isotropic turbulence and uniform shear flows; and a spectral turbulence model for wall-bounded flows without resolving the boundary layer. Data from direct numerical simulation (DNS) of turbulent channel flow was used to help develop this model and to assess its performance in the 1D direction across the channel width. The resultant spectral model is capable of predicting the mean velocity, turbulent kinetic energy and energy spectrum distributions for single-phase wall-bounded flows all the way to the wall, where the model source terms have been developed to account for the wall influence. We implemented the model into the 3D multiphase CFD code NPHASE-CMFD and the latest results are within reasonable error of the 1D predictions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Direct numerical simulation of reactor two-phase flows enabled by high-performance computing

Nuclear reactor two-phase flows remain a great engineering challenge, where the high-resolution two-phase flow database which can inform practical model development is still sparse due to the extreme reactor operation conditions and measurement difficulties. Owing to the rapid growth of computing power, the direct numerical simulation (DNS) is enjoying a renewed interest in investigating the related flow problems. A combination between DNS and an interface tracking method can provide a unique opportunity to study two-phase flows based on first principles calculations. More importantly, state-of-the-art high-performance computing (HPC) facilities are helping unlock this great potential. This paper reviews the recent research progress of two-phase flow DNS related to reactor applications. The progress in large-scale bubbly flow DNS has been focused not only on the sheer size of those simulations in terms of resolved Reynolds number, but also on the associated advanced modeling and analysis techniques. Specifically, the current areas of active research include modeling of sub-cooled boiling, bubble coalescence, as well as the advanced post-processing toolkit for bubbly flow simulations in reactor geometries. A novel bubble tracking method has been developed to track the evolution of bubbles in two-phase bubbly flow. Also, spectral analysis of DNS database in different geometries has been performed to investigate the modulation of the energy spectrum slope due to bubble-induced turbulence. In addition, the single-and two-phase analysis results are presented for turbulent flows within the pressurized water reactor (PWR) core geometries. The related simulations are possible to carry out only with the world leading HPC platforms. Furthermore, these simulations are allowing more complex turbulence model development and validation for use in 3D multiphase computational fluid dynamics (M-CFD) codes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

High-fidelity pool boiling simulations on multiple nucleation sites using interface capturing method

Boiling has proved to be one of the most efficient means for heat transfer and is a very important phenomenon during severe accident scenarios in light water reactors. High-fidelity pool boiling simulations can provide a numerical database for improving mechanistic boiling models by allowing for specific evaluation of interactions among bubbles. Previously published pool boiling simulations investigated two nucleation sites in which bubble growth at one site suppressed nucleation at the other site. Based on previous study results, more complicated interface-capturing simulations on pool boiling were conducted using PHASTA code with locally refined unstructured mesh. First, different boundary conditions (BCs) were assessed to support robustness and reproducibility of the boiling model. Then, a scale study was conducted at a larger domain with nine nucleation sites where either nine or four nucleation sites are activated. Involving more nucleation sites increased the complexity of bubble interactions from surrounding sites. Finally, bubble departure behavior influenced by wall heat flux was investigated. In this work, when heat flux was increased, the order of bubble departure changed, but diagonal bubbles always departed after one another. The departure time interval between the first and second bubble reduced as heat flux increased. The corresponding frequency was almost linearly proportional to the heat flux. In addition, bubble departure behavior was found to be greatly influenced by the nucleation site pattern. Multiple nucleation sites resulted in superimposed inhibitive effects from surrounding sites to each bubble, which extensively delayed the departure. This new observation was not discussed in previously published works. The work presented here provides new insight on the fundamental understanding of boiling phenomena, contributes to the development of a 3D multiphase computational fluid dynamics (M-CFD) model, and provides a more comprehensive database for data-driven pool boiling studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗

DOE OSTI2018

Integrated simulation of vertical fracture propagation induced by water injection and its borehole electromagnetic responses in shale gas systems

We analyze fracture propagation induced by hydraulic fracturing with water injection and examine their detectability with crosswell electromagnetic (EM) geophysical methods. For rigorous 3D coupled flow-geomechanical modeling, we employ a numerical method that can model failure by tensile and shear stresses, dynamic nonlinear permeability, dual continuum approach, and thermo-poro-mechanical effects. From numerical simulation, we find that the fracture propagation is not the same as propagation of the water front, because fracturing is governed by geomechanics whereas water saturation is determined by multiphase flow. At early times, the water front is almost identical to the fracture tip, suggesting that the fracture is mostly filled with the injected water. However, at late times, movement of the water front is retarded compared to fracture propagation, yielding a significant gap between the water front and the fracture top, filled with reservoir gas. During fracture propagation, the coupled flow-geomechanical models are transformed via a rock-physics model into electrical conductivity models. We employ a full 3D finite-element EM geophysical simulator to evaluate the sensitivity of the crosswell EM method to fracture propagation. It is shown that anomalous distribution of electrical conductivity is closely related to the injected water saturation, but not closely related to newly created unsaturated fractures. Our numerical modeling experiments demonstrate that the crosswell EM method can be highly sensitive to electrical conductivity changes that directly indicate the migration pathways of the injected fluid. Accordingly, the EM method can serve as an effective monitoring tool for monitoring the injected fluids during hydraulic fracturing operations.

58 GEOSCIENCES↗

DOE OSTI2018

Above zone pressure interpretation for leaky well characterization and its identification from leaky caprock/fault

Pressure transient analysis has long been used for reservoir characterization. Above-zone (AZ) pressure has been recently investigated for inferring leakage pathway characteristics in leakage events from subsurface injection operations. The recorded pressure in the AZ should be purely related to leakage and therefore it can be safely inverted to deduce leakage characteristics. It is crucial to evaluate fluid leakage through abandoned wells to plan for further measurements of leakage prevention. However, the recorded AZ leakage signal may not be related to leaky well(s). Therefore, identification and spatial investigation of well leakage is required for leakage evaluation. Here in this paper, we propose a pressure interpretation method for early detection of leaky pathways, applying two observation points in the AZ. We distinguish leaky well, fault and caprock based on their corresponding flow regime identification. We show that the pressure difference of the two observation wells can be applied as a proxy for unknown leakage rate, which is crucial for leakage identification as well as characterization. Results show that the estimated location of the leaky well, leakage coefficient, and the leakage rate are in good agreement with the actual values. The estimated leakage coefficient of the leaky well can be used to evaluate well leakage in multiphase systems such as carbon dioxide leakage in deep saline aquifers.

02 PETROLEUM↗

DOE OSTI2017

Chemical freeze-out in relativistic heavy-ion collisions

One surprising result in relativistic heavy-ion collisions is that the abundance of various particles measured in experiments is consistent with the picture that they reach chemical equilibrium at a temperature much higher than the temperature they freeze out kinetically. Using a multiphase transport model to study particle production in these collisions, we find, as an example, that the effective pion to nucleon ratio, which includes those from resonance decays, indeed changes very little during the evolution of the hadronic matter from the chemical to the kinetic freeze-out, and it is also accompanied by an almost constant specific entropy. Finally, we further use a hadron resonance gas model to illustrate the results from the transport model study.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Multiplicity scaling of light nuclei production in relativistic heavy-ion collisions

Using the nucleon coalescence model based on kinetic freeze-out nucleons from the 3D MUSIC+UrQMD and the 2D VISHNU hybrid model with a crossover equation of state, we study the multiplicity dependence of deuteron ( d ) and triton ( t ) production from central to peripheral Au+Au collisions at s NN = 7.7, 14.5, 19.6, 27, 39, 62.4 and 200 GeV and Pb+Pb at s NN = 2.76 TeV, respectively. It is found that the ratio N t N p / N d 2 of the proton yield N p , deuteron yield N d and triton yield N t exhibits a scaling behavior in its multiplicity dependence, i.e., decreasing monotonically with increasing charged-particle multiplicity. A similar multiplicity scaling of this ratio is also found in the nucleon coalescence calculation based on kinetic freeze-out nucleons from a multiphase transport (AMPT) model. The scaling behavior of N t N p / N d 2 can be naturally explained by the interplay between the sizes of light nuclei and the nucleon emission source. We further argue that the multiplicity scaling of N t N p / N d 2 can be used to validate the production mechanism of light nuclei, and the collision energy dependence of this yield ratio can further serve as a baseline in the search for the QCD critical point in relativistic heavy-ion collisions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Impact of globally spin-aligned vector mesons on the search for the chiral magnetic effect in heavy-ion collisions

In high-energy heavy-ion collisions, the chiral magnetic effect (CME) is predicted to arise from the interplay between the chirality imbalance of quarks in the nuclear medium and the intense magnetic field, and will cause a charge separation along the magnetic field direction. While the search for the CME is still ongoing in experiments at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), the CME-sensitive observables need to be scrutinized to exclude the non-CME contributions. In this work, we examine the influence of globally spin-aligned ρ mesons on the γ 112 correlator, the R Ψ 2 (ΔS) correlator, and the signed balance functions, via a toy model and a multiphase transport model (AMPT). The global spin alignment of vector mesons could originate from non-CME mechanisms in heavy-ion collisions, and is characterized by the 00-component of the spin density matrix, ρ 00 . Here, we find that the CME observables show similar dependence on ρ 00 , and could receive a positive (negative) contribution from ρ-decay pions, if the ρ 00 of ρ mesons is larger (smaller) than 1/3. Since pions are the most abundant particles in such collisions, the ρ 00 measurements for ρ mesons become crucial to the interpretation of the CME data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Responsive lignin for shape memory applications

Here, we analyze thermally unstable lignin linkages that generate reactive sites and promote crosslinking reactions during thermal annealing at 180 °C. This results in significantly enhanced glass transition temperature (Tg). We use these characteristics of lignin to enhance interfacial reactivity of a lignin-based multiphase polymer synthesized by a dynamic shear process. An approximately 18 °C increase in Tg and more than 230% increase in storage modulus (E') are achieved by thermally annealing the 50 wt% hardwood lignin-nitrile rubber composite for 5 h at 180 °C. Tunable chemical and physical crosslinks within lignin and rubber resulted in shape programmability demonstrating excellent strain recovery of the synthesized renewable materials. The substantial improvement of the elastic work density of the lignin-based elastomers allows a broader applicable stress window for stress/strain sensing. The renewable materials also exhibit a better strain recovery after thermal annealing.

36 MATERIALS SCIENCE↗

DOE OSTI2017

CFD-DEM modeling the effect of column size and bed height on minimum fluidization velocity in micro fluidized beds with Geldart B particles

The fluidization behavior of Geldart B particles in micro fluidized beds is investigated numerically using Computational Fluid Dynamics coupled with Discrete Element Method (CFD-DEM) available in the open-source Multiphase Flow with Interphase eXchanges (MFIX) code. The effects of different bed inner diameters (D) of 8 mm, 12 mm, 16 mm and various initial static bed heights (H) were examined. It is found that both decreasing the column diameter and increasing the bed height in a micro fluidized bed increases the minimum fluidization velocity (Umf). The observed overshoot in pressure drop that occurs before the onset of fluidization decreases in magnitude with increasing column diameter, however there is less sensitivity to bed height. Overall, the numerical results agree qualitatively with existing theoretical correlations and experimental studies. The simulations show that both column diameter and particle-wall friction contribute to the variation in minimum fluidization velocity. Finally, these two factors are coupled and hard to separate. The detailed influences of wall friction on minimum fluidization velocity are then investigated for a prescribed column diameter of 8 mm by varying the wall friction from 0 to 0.4.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Particles climbing along a vertically vibrating tube: numerical simulation using the Discrete Element Method (DEM)

It has been reported experimentally that granular particles can climb along a vertically vibrating tube partially inserted inside a granular silo. Here, we use the Discrete Element Method (DEM) available in the Multiphase Flow with Interphase eXchanges (MFIX) code to investigate this phenomenon. By tracking the movement of individual particles, the climbing mechanism was illustrated and analyzed. The numerical results show that a sufficiently high vibration strength is needed to form a low solids volume fraction region inside the lower end of the vibrating tube, a dense region in the middle of the tube, and to bring the particles outside from the top layers down to fill in the void. The results also show that particle compaction in the middle section of the tube is the main cause of the climbing. Consequently, varying parameters which influence the compacted region, such as the restitution coefficient, change the climbing height.

42 ENGINEERING↗

DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗