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.

SourceResearch connectionFrom yearThrough year

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

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Improvements in the reactor core isolation cooling (RCIC) pump model

As part of the efforts to understand the unexpected “self-regulating” mode of the RCIC (Reactor Core Isolation Cooling) systems that was observed during the Fukushima accidents and to extend BWR RCIC and PWR AFW (Auxiliary Feed Water) operational range and flexibility, mechanistic models for the Terry turbine are being developed and tested with the RELAP-7 code to simulate the RCIC system. In the previous work, we developed a set of analytical models for the normal working conditions, as well as the two-phase off-design conditions, of the RCIC system, based on Sandia’s original work. An under-expanded jet model was developed to replace the CFD-based reduced-order model in the Sandia work to obtain the velocity and thermodynamic conditions of the turbine stator inlet. The models include both an adiabatic expansion process inside the nozzle and a free expansion process outside of the nozzle to ambient pressure. For the two-phase off-design conditions, two well-established choking models were used - the Isentropic Homogenous Equilibrium Model (IHEM) and Moody’s model. The new twophase Terry turbine model used the choking models to calculate the mass flow rate, the critical pressure at the nozzle throat, and steam quality. The two-phase nozzle expansion model was also developed by considering the vapor phase with a similar model for the single-phase case with the assumption that the liquid phase would slip along the wall with a much slower speed and will not contribute to the impulse on the rotor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Flow boiling of R134a in a large surface area microchannel array for high-flux laser diode cooling

Packaging high average power laser diode arrays that generate heat at an area average flux in excess of 1 kW·cm -2 is a significant engineering challenge. While liquid microchannel coolers have demonstrated up to 11.9 kW·cm -2 , two-phase microchannel array coolers have not achieved 1 kW cm -2 due to critical heat flux and flow instabilities. In the current study, flow boiling heat transfer was characterized by a 1 × 10 mm heated zone centered over a 5 × 10 mm array of 125 very small channels (45 × 200 μm) with R134a as the phase change fluid. The high aspect ratio channels (4.4:1) were manufactured using MEMS fabrication techniques, which yielded a large heat transfer surface area to volume ratio. A test facility was used to characterize the heat transfer performance of boiling R134a over a range of saturation temperatures (15°C to 25°C), mass fluxes (735-2230 kg·m -2 ·s -1 ), and heat duties (< 110.3 W). During the tests, the calculated outlet vapor quality exceeded 61%, and the base heat flux at the heater reached a maximum of 1.1 kW·cm -2 . The resulting average experimental flow boiling heat transfer coefficients are found to be as large a 13.4 kW·m -2 ·K -1 over the approximately 3 mm two-phase region, with an average uncertainty of ± 2.72%. A substantial amount of heat was spread downstream via the low thermal resistance silicon floor. Specifically, between 29.5% and 55.1% of the heat dissipated in the two-phase region was dissipated over the heater. Finally, the remaining heat dissipated in the two-phase region was dissipated in the 2 mm of channel downstream of the heater. This suggests that heat spreading from the hotspot played a vital role in dissipating the heat load.

42 ENGINEERING↗

DOE OSTI2019

A FFT-based numerical implementation of mesoscale field dislocation mechanics: Application to two-phase laminates

In this study, we present an enhanced crystal plasticity elasto-viscoplastic fast Fourier transform (EVPFFT) formulation coupled with a phenomenological Mesoscale Field Dislocation Mechanics (MFDM) theory here named MFDM-EVPFFT formulation. In contrast with classic CP-EVPFFT, the model is able to tackle plastic flow and hardening due to polar dislocation density distributions or geometrically necessary dislocations (GNDs) in addition to statistically stored dislocations (SSDs). The model also considers GND mobility through a GND density evolution law numerically solved with a recently developed filtered spectral approach, which is here coupled with stress equilibrium. The discrete Fourier transform method combined with finite differences is applied to solve both lattice incompatibility and Lippmann-Schwinger equations in an augmented Lagrangian numerical scheme. Numerical results are presented for two-phase laminate composites with plastic channels and elastic second phase. It is shown that both GND densities and slip constraint at phase boundaries influence the overall and local hardening behavior. In contrast with the CP-EVPFFT formulation, a channel size effect is predicted on the shear flow stress with the present MFDM-EVPFFT formulation. The size effect originates from the progressive formation of continuous screw GND pile-ups from phase boundaries to the channel center. The effect of GND mean free path on local and global responses is also examined for the two-phase composite.

36 MATERIALS SCIENCE↗

DOE OSTI2021

On Two-Phase Pressure and Temperature Equilibration with Mie-Grüneisen Equations of State

In the course of reactive flow modeling, the phases of high explosive (HE) reactants and products introduce additional degrees of freedom to the underlying hydrodynamic equations of compressible reactive fluid mechanics. Hydrodynamic simulations track the bulk (mixture) conservation properties of mass, momentum and energy. When multiple phases are introduced, one must decide on appropriate closure rules to determine the relative phase volumes and energy distribution between the phases. One typical closure is to assume that the phases are in pressure and temperature equilibrium. Here, the requisite equations of p-T equilibration are explored for Mie-Grüneisen Equations of State (EOS). Examples with Davis reactants and Davis products EOSs will be given. The Davis forms are used in both the WSD and AWSD reactive flow models.

42 ENGINEERING↗

DOE OSTI2022

Geologic Heterogeneity Controls on Trapping and Migration of CO 2

When CO 2 is injected into sedimentary basins, the fate of the buoyant CO 2 will be determined by how it interacts with the underlying geologic heterogeneities. Here we geologically print a cross-bedded pattern using two different bead sizes, and then conduct flow experiments to observe buoyant flow pathways, migration speed, and immobilized volumes. The amount of buoyant phase trapped under the heterogeneities is observed to vary by two orders of magnitude, with the controlling parameter being the size contrast between the coarse and fine grains that make up the structure. Furthermore, the results allow an estimate of heterogeneity trapping, and show how CO 2 trapping in aquifers can be experimentally evaluated before injection.

58 GEOSCIENCES↗

DOE OSTI2018

Fluid Flow, Dissolution, and Mixing Phenomena in Argon-Stirred Steel Ladles

In the current study, the multiphase fluid flow in argon-stirred steel ladles is simulated using an Eulerian–Lagrangian two-phase approach. The momentum source and the turbulent kinetic energy source due to the motion of the bubble are considered for the liquid phase. Argon bubbles are treated as discrete phase particles, and the interfacial forces between the bubbles and the liquid phase; the dependence of the gas density and the bubble diameter on the temperature and the static pressure; and the bubble size distribution are considered. When the fluid flow reaches the quasi-steady state, the ferroalloy melting and mixing phenomena is also modeled. The melting time and the trajectory length of each ferroalloy particle are recorded using a user-defined function (UDF). Local mixing time is predicted in the entire computational domain by checking the mixing criteria in every cell. The effects of gas flow rate, porous plug location, and separation angle of two porous plugs on the fluid flow and the mixing phenomena are investigated. The results show that the flow intensity increases, and the mixing time decreases with the increasing gas flow rate. The optimal porous plug’s radial position with one porous plug is 0.50R for its best mixing condition. When two porous plugs are adopted, the separation angle of 90 deg is recommended to improve the flow field and mixing phenomena.

42 ENGINEERING↗

DOE OSTI2018

Renormalization group flow of the Aharonov–Bohm scattering amplitude

The Aharonov–Bohm elastic scattering with incident particles described by plane waves is revisited by using the phase-shifts method. The formal equivalence between the cylindrical Schrödinger equation and the one-dimensional Calogero problem allows us to show that up to two scattering phase-shifts modes in the cylindrical waves expansion must be renormalized. The renormalization procedure introduces new length scales giving rise to spontaneous breaking of the conformal symmetry. The knowledge of the exact beta function permits us to describe the renormalization group flows within the two-parametric family of renormalized Aharonov–Bohm scattering amplitudes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Experimental Investigation of Uranium and Iron Condensation from High-Temperature Plasma Conditions

We used a plasma flow reactor (PFR) to generate synthetic fallout nanoparticles from vapor-phase condensation of two different input concentrations of uranium and iron analytes (U/Fe = 1:1 and 1:2). Synthetic fallout from complex chemical matrices (e.g., mixtures of U and Fe) has not been generated in this setup before and allows for the observation of relative condensation and fractionation of nuclear debris. Experiments were conducted under two different temperature histories with variations in particle flow patterns along the PFR. Transmission electron microscopy (TEM) observation and analysis of the nanoparticles revealed variations in speciation of uranium oxides (UO 2 and α-UO 3 ) depending on the competition between flow mixing and oxygen sequestration by iron. A ternary metal oxide, UFeO 4 , was observed in addition to iron oxides (e.g., FeO and Fe 3 O 4 ), which suggests that fallout models should account for chemical speciation of ternary metal oxides (i.e., UFeO 4 ) and their relative condensation behaviors in addition to those of singular metal oxides (e.g., FeO and UO 2 ). X-ray Energy Dispersive Spectroscopy (EDS) elemental maps showed that some particles had Fe-rich cores surrounded by U-rich regions. This suggests that either condensed U oxides coagulate onto molten Fe oxides or that the increase in iron analyte concentration might drive the system toward a higher degree of supersaturation, leading to earlier formation of iron oxide particles and providing an energetically favored pathway for nucleation of uranium oxides around iron oxide particles.

and nuclear chemistry↗

DOE OSTI2020

effect of ambient pressure on the behavior of single-component fuels in a gasoline multi-hole injector

The injection characteristics of neat ethanol and pure isooctane are studied under different ambient pressure and temperature conditions. Injection underflash-boiling conditions can enhance liquid atomization and evaporation, providing the possibility of improvement in the fuel/air mixing. These superheated conditions often introduce phenomena that are not taken into account in the standard modeling of sprays for engine applications. The present work proposes a numerical investigation of the behavior of Engine Combustion Network's 8-hole spray-G injector, starting at the subcooled nominal condition and reducing the ambient pressure at constant low temperature to reach the flare flash-boiling condition. To initialize the properties of the injected fuel, the flow in the nozzle is simulated with a Eulerian approach, handling the two phases with a mixture model and the phase change, due to cavitation and flash boiling, with the Homogenous Relaxation Model. A map of the mixture's kinematic and thermal behavior is obtained at the interface between the injector and the chamber to initialize the Lagrangian simulations. A literature-based vaporization model is implemented to obtain the proper description of the characteristic features of a multi-hole spray under super-heated conditions, like plume-plume interaction. The numerical representation of the spray is validated in terms of penetration and radial spreading on DBI images, reproducing the light attenuation profiles caused by the presence of the liquid spray. Simulations show that coupled nozzle flow and spray calculations capture the spray morphology and shape better compared to calculations performed without considering the nozzle flow simulation details, especially under flare flash conditions.

Nocivelli, Lorenzo↗

DOE OSTI2020

Quantum interference of currents in an atomtronic SQUID

Quantum interference of currents is the most important and well known quantum phenomenon in a conventional superconducting quantum interference device (SQUID). Here, we report the observation of quantum interference of currents in an atomtronic SQUID. Analogous to a conventional SQUID, currents flowing through two junctions in an atomtronic SQUID interfere due to the phase difference from rotation. This interference results in modulation of critical currents. This modulation was observed for three different radii with clear modulation periods which were measured to be consistent with fundamental rotation rates. This observation shows the possibility of studying various interesting SQUID physics with an atomtronic SQUID and especially, macroscopic quantum phenomena with currents may be realized with an atomtronic SQUID toward the goal of quantum metrology of rotation sensing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Development of an Efficient Meso- scale Multi-phase Flow Solver in Nuclear Applications

The proposed research aims at formulating a predictive high-order Lattice Boltzmann Equation for multi-phase flows relevant to nuclear energy related application - namely, saturated and sub-cooled boiling in reactors, and liquid- liquid mixing and extraction for fuel cycle separation. An efficient flow solver will be developed based on the Finite Element based Lattice Boltzmann Method (FE- LBM), accounting for phase-change heat transfer and capable of treating multiple phases over length scales from the submicron to the meter. A thermal LBM will be developed in order to handle adjustable Prandtl number, arbitrary specific heat ratio, a wide range of temperature variations, better numerical stability during liquid-vapor phase change, and full thermo-hydrodynamic consistency. Two-phase FE-LBM will be extended to liquid–liquid–gas multi-phase flows for application to high-fidelity simulations building up from the meso-scale up to the equipment sub-component scale. While several relevant applications exist, the initial applications for demonstration of the efficient methods to be developed as part of this project include numerical investigations of Critical Heat Flux (CHF) phenomena in nuclear reactor fuel bundles, and liquid-liquid mixing and interfacial area generation for liquid-liquid separations. In addition, targeted experiments will be conducted for validation of this advanced multi-phase model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Formation of Magnetically Truncated Accretion Disks in 3D Radiation-transport Two-temperature GRMHD Simulations

Multiwavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nucleus transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole (BH). However, the formation, structure, and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semiopaque corona reach, respectively, $T_i$ ≳ 10 10 K and $T_e$ ≳ 5 × 10 8 K. The system produces radiation, powerful collimated jets, and broader winds at the total energy efficiency exceeding 90%, the highest ever energy extraction efficiency from a spinning BH by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

A Numerical Study of Factors Affecting Fracture-Fluid Cleanup and Produced Gas/Water in Marcellus Shale: Part II

Horizontal wells combined with successful multi-stage hydraulic fracture treatments are currently the most established method for effectively stimulating and enabling economic development of gas bearing organic-rich shale formations. Fracture cleanup in the Stimulated Reservoir Volume (SRV) is critical to stimulation effectiveness and long-term well performance. However, fluid cleanup is often hampered by formation damage, and post-fracture well performance frequently falls below expectations. A systematic study of the factors that hinder fracture fluid cleanup in shale formations can help optimize fracture treatments and better quantify long term volumes of produced water and gas. Fracture fluid cleanup is a complex process influenced by multi-phase flow through porous media (relative permeability hysteresis, capillary pressure etc.), reservoir rock and fluid properties, fracture fluid properties, proppant placement, fracture treatment parameters, and subsequent flowback and field operations. Changing SRV and fracture conductivity as production progresses further adds to the complexity of this problem. Numerical simulation is the best, and most practical approach to investigate such a complicated blend of mechanisms, parameters, their interactions, and subsequent impact on fracture fluid cleanup and well deliverability. In this paper, a 3-dimensional, 2-phase, dual-porosity model was used to investigate the impact of multiphase flow, proppant crushing, proppant diagenesis, shut-in time, reservoir rock compaction, gas slippage, and gas desorption on fracture fluid cleanup, and well performance in Marcellus shale. The research findings have shed light on the factors that substantially constrains efficient fracture fluid cleanup in gas shales, and provided guidelines for improved fracture treatment designs and water management.

58 GEOSCIENCES↗

DOE OSTI2017

Toward High Fidelity LES Simulations of 5 x 5 PWR Fuel Bundle with Mixing Vane

The present work illustrates one example of a computational fluid dynamics (CFD) application in the nuclear engineering field. Since CFD is one of the most powerful tools available at the actual state of the art that could in principle help and assist the design process, a compromise should be established between the accuracy of the CFD numerical model (RANS - LES - DNS) and the computational cost of the simulations. In this spirit, this work will present the computational results from two different CFD model categories, Large Eddy Simulations (LES) and Reynolds-averaged Navier-Stokes (RANS) models. Both treat the flow turbulence issue in two separate ways, which differ mostly in the resolution of the turbulence statistics and the computational cost. The single-phase hydraulic problem presented involves the steady-state turbulent flow field modelling and resolution around a scaled PWR mixing vane grid, designed by Westinghouse and tested by Texas A and M University. The role of the mixing vane/spacer grid is to enhance the turbulent mixing and hence heat transfer of the core coolant around fuel bundles, thus increasing the departure of nucleate boiling-critical heat flux (DNB-CHF) point, as well as important structural functions for the rods bundles. The authors numerically studied the flow field characteristics at the region downstream of the spacer grids and validated their results to particle imagery velocimetry (PIV) measurements. Dominguez and Hassan performed two-dimensional time resolved PIV (TR-PIV) measurements within a 5 x 5 PWR rod bundle with spacer-grids that utilized a matching index of refraction (MIR) approach. Their averaged data showed a rather uniform direction velocity profile after the spacer grid while the downstream illustrated a no slip wall profile. The available PIV flow field experimental results from Texas A and M show complex structures of the flow field downstream the mixing vanes with the formation of secondary flows that tend to enhance the turbulent heat transfer, due to the mixing vanes geometry itself. Computational results will be validated with the available experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2018

A two-phase mixing layer between parallel gas and liquid streams: multiphase turbulence statistics and influence of interfacial instability

The two-phase mixing layer formed between parallel gas and liquid streams is an important fundamental problem in turbulent multiphase flows. The problem is relevant to many industrial applications and natural phenomena, such as air-blast atomizers in fuel injection systems and breaking waves in the ocean. The velocity difference between the gas and liquid streams triggers an interfacial instability which can be convective or absolute depending on the stream properties and injection parameters. In the present study, a direct numerical simulation of a two-phase gas–liquid mixing layer that lie in the absolute instability regime is conducted. A dominant frequency is observed in the simulation and the numerical result agrees well with the prediction from viscous stability theory. As the interfacial wave plays a critical role in turbulence transition and development, the temporal evolution of turbulent fluctuations (such as the enstrophy) also exhibits a similar frequency. To investigate the statistical response of the multiphase turbulence flow, the simulation has been run for a long physical time so that time-averaging can be performed to yield the statistically converged results for Reynolds stresses and the turbulent kinetic energy (TKE) budget. An extensive mesh refinement study using from 8 million to about 4 billions cells has been performed. The turbulent dissipation is shown to be highly demanding on mesh resolution compared with other terms in TKE budget. The results obtained with the finest mesh are shown to be close to converged results of turbulent dissipation which allow us to obtain estimations of the Kolmogorov and Hinze scales. The estimated Kolmogorov scale is found to be similar to the cell size of the finest mesh used here. The computed Hinze scale is significantly larger than the size of droplets observed and does not seem to be a relevant length scale to describe the smallest size of droplets formed in atomization.

Mechanics↗

DOE OSTI2020

Continuous and discontinuous transitions in the depinning of two-dimensional dusty plasmas on a one-dimensional periodic substrate

We report Langevin dynamical simulations are performed to study the depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate. From the diagnostics of the sixfold coordinated particles P 6 and the collective drift velocity V x , three different states appear, which are the pinning, disordered plastic flow, and moving ordered states. It is found that the depth of the substrate is able to modulate the properties of the depinning phase transition, based on the results of P 6 and V x , as well as the observation of hysteresis of V x while increasing and decreasing the driving force monotonically. When the depth of the substrate is shallow, there are two continuous phase transitions. When the potential well depth slightly increases, the phase transition from the pinned to the disordered plastic flow states is continuous; however, the phase transition from the disordered plastic flow to the moving ordered states is discontinuous. When the substrate is even deeper, the phase transition from the pinned to the disordered plastic flow states changes to discontinuous. When the depth of the substrate further increases, as the driving force increases, the pinned state changes to the moving ordered state directly, so that the disordered plastic flow state disappears completely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Parallel Fokker–Planck-DSMC algorithm for rarefied gas flow simulation in complex domains at all Knudsen numbers

A major challenge for the conventional Direct Simulation Monte Carlo (DSMC) technique lies in the fact that its computational cost becomes prohibitive in the near continuum regime, where the Knudsen number (Kn)—characterizing the degree of rarefaction—becomes small. In contrast, the Fokker–Planck (FP) based particle Monte Carlo scheme allows for computationally efficient simulations of rarefied gas flows in the low and intermediate Kn regime. The Fokker–Planck collision operator—instead of performing binary collisions employed by the DSMC method—integrates continuous stochastic processes for the phase space evolution in time. This allows for time step and grid cell sizes larger than the respective collisional scales required by DSMC. Dynamically switching between the FP and the DSMC collision operators in each computational cell is the basis of the combined FP-DSMC method, which has been proven successful in simulating flows covering the whole Kn range. Until recently, this algorithm had only been applied to two-dimensional test cases. In this contribution, we present the first general purpose implementation of the combined FP-DSMC method. Utilizing both shared- and distributed-memory parallelization, this implementation provides the capability for simulations involving many particles and complex geometries by exploiting state of the art computer cluster technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗