Search NASA⌕ Search

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.

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

At least 145 records · Page 8

Integrating geophysical monitoring data into multiphase fluid flow reservoir simulation

Simulation of multiphase flow systems are of critical importance in managing hydrological systems. Flow simulations are affected by a number of factors including structure and flow properties including porosity and permeability as well as the anisotropy and heterogeneity of these properties. In many cases traditional hydrological and reservoir data are highly affected by these parameters, but are not directly sensitive to them. As such modellers often adjust these parameters in an ad hoc manner until solutions numerically converge. Simulation models are generally based on structural data from reflection seismics whose physical flow properties are then populated using geostatistical extrapolation techniques utilizing a sparse number of borehole logs and core analysis. In multiphase systems including enhanced oil recovery and carbon capture and sequestration uncertainties regarding phase-dependent physical properties confounds this challenge further. Geophysical methods provide a means by which to gain an improved understanding of phase distributions in the subsurface. In this paper we will look at applications from active carbon capture and sequestration and enhanced oil recovery applications, as well as synthetic examples. Geophysical data including electromagnetic and gravity are inverted using structural constraints from the reservoir model. Inversions are then mapped into flow properties using calibrated relations such as Archie’s Equation. The coupled models can then be used to both verify and improve on the reservoir flow model which improves it’s predictive power and utility as a management tool.

02 PETROLEUM↗

Velocity structure functions in multiphase turbulence: interpreting kinematics of Hα filaments in cool-core clusters

ABSTRACT The central regions of cool-core galaxy clusters harbour multiphase gas, with gas temperatures ranging from $10$ to $10^7\, \mathrm{K}$. Feedback from active galactic nuclei jets prevents the gas from undergoing a catastrophic cooling flow. However, the exact mechanism of this feedback energy input is unknown, mainly due to the lack of velocity measurements of the hot-phase gas. However, recent observations have measured the velocity structure functions (VSFs) of the cooler molecular (${\sim} 10\, \mathrm{K}$) and Hα filaments (${\sim} 10^4\, \mathrm{K}$) and used them to indirectly estimate the motions of the hot phase. In the first part of this study, we conduct high-resolution (3843–15363 resolution elements) simulations of homogeneous isotropic subsonic turbulence, without radiative cooling. We analyse the second-order velocity structure functions (VSF2) in these simulations and study the effects of varying spatial resolution, the introduction of magnetic fields, and the effect of projection along the line of sight (LOS) on it. In the second part of the study, we analyse high-resolution (7683 resolution elements) idealized simulations of multiphase turbulence in the intracluster medium from the companion study Mohapatra et al. We compare the VSF2 for both the hot ($T\sim 10^7\, \mathrm{K}$) and cold ($T\sim 10^4\, \mathrm{K}$) phases and find that their amplitude depends on the density contrast between the phases. They have similar scaling with separation, but introducing magnetic fields steepens the VSF2 of only the cold phase. We also find that projection along the LOS steepens the VSF2 for the hot phase and mostly flattens it for the cold phase.

Mohapatra, Rajsekhar (ORCID:0000000216007552)↗

Modeling Multiphase Effects in CO 2 Compressors at Subcritical Inlet Conditions

An advanced numerical framework to model CO 2 compressors over a wide range of subcritical conditions is presented in this paper. Thermodynamic and transport properties are obtained through a table look-up procedure with specialized features for subcritical conditions. Phase change is triggered by the difference between the local values of pressure and saturation pressure, and both vaporization and condensation can be modeled. Rigorous validation of the framework is presented for condensation in high pressure CO 2 using test data in a De Laval nozzle. The comparisons between computations and test data include: condensation onset locations, Wilson line, and nozzle pressure profiles as a function of inlet pressures. The framework is applied to the Sandia compressor that has been modeled over broad range of conditions spanning the saturation dome including: near critical inlet conditions (305.4 K, and 7.843 MPa), pure liquid inlet conditions (at 295 K), pure vapor inlet conditions (at 302 K), and two-phase inlet conditions (at 290 K). Multiphase effects ranging from cavitation at the liquid line to condensation at the vapor line have been simulated. Here, the role of real fluid effects in enhancing multiphase effects at elevated temperatures closer to the critical point has been identified. The performance of the compressor has been compared with test data; the computational fluid dynamics (CFD) results also show that the head-flow coefficient curve collapses with relatively minor scatter, similar to the test data, when the flow coefficient is defined on the impeller exit meridional velocity.

14 SOLAR ENERGY↗

Effects of Carbonic Acid-Rock Interactions on CO2/Brine Multiphase Flow Properties in the Upper Minnelusa Sandstones

Summary Carbon dioxide (CO2) injection into a deep saline aquifer can dissolve into formation brine and generate carbonic acid. The resulting acid can drive fluid-rock geochemical reactions. The impact of these fluid-rock geochemical reactions on porosity, permeability, and multiphase flow responses is relevant to the determination of CO2 storage capacity of deep saline aquifers. In this research, carbonic acid flooding experiments were performed on core samples consisting of poorly sorted, quartz-rich sand with laminated bedding from a possible CO2 storage target in northwest Wyoming. Complementary pre- and post-injection porosity and permeability, thin-section, Brunauer-Emmett-Teller (BET) surface area, mercury intrusion capillary pressure (MICP), and time-domain nuclear magnetic resonance (TD-NMR) measurements were conducted. Overall, both core porosity and permeability increased after a 7-day carbonic acid injection, from 6.2 to 8.4% and 1.6 to 3.7 md, respectively. We attributed these changes to carbonate mineral dissolution, which was evidenced by the effluent brine geochemistry, pore-throat size distribution (PTSD), and BET surface area. To be more specific, within the more-permeable section of core samples containing larger pore size, the permeability increment is apparent due to dolomite mineral grains and cements dissolution. However, for the lower-permeability section corresponding to the smaller pore size, mineral precipitation possibly lessened dissolution effects, leading to insignificant petrophysical properties changes. Consequently, the observed heterogeneous carbonic acid-rock interactions resulted in alterations of CO2/brine relative permeability (i.e., the initial CO2 saturation decreased and the CO2 flow capacity was enhanced). This research provides a fundamental understanding regarding effects of fluid-rock reactions on changes in static and multiphase flow properties of eolian sandstones, which lays the foundation for more accurate prediction/simulation of CO2 injection into deep saline aquifers.

Engineering↗

Nonequilibrium Physics and Phase-Field Modeling of Multiphase Flow in Porous Media

The overarching goal of this project was to develop a new continuum theory of multiphase flow in porous media. The theory follows a phase-field modeling approach, and therefore has a sound thermodynamical basis. It is a phenomenological theory in the sense that its formulation is driven by macroscopic phenomena, such as viscous instabilities during multifluid displacement. The research agenda was organized around a set of hypothesis on hitherto unexplained behavior of multiphase flow. All these hypothesis are nontrivial, and testable. Indeed, a central aspect of the project was testing each hypothesis by means of carefully-designed laboratory experiments, therefore probing the validity of the proposed theory. The proposed research places an emphasis on the fundamentals of flow physics, but is motivated by important energy-driven applications in earth sciences, as well as microfluidic technology.

42 ENGINEERING↗

Development of Multi-physics (Multiphase CFD + MCNP) simulation for generic solution vessel power calculation

The current study aims to predict the steady state power of a generic solution vessel and to develop a corresponding heat transfer coefficient correlation for a Moly99 production facility by conducting a fully coupled multi-physics simulation. A prediction of steady state power for the current application is inherently interconnected between thermal hydraulic characteristics (i.e. Multiphase computational fluid dynamics solved by ANSYS-Fluent 17.2) and the corresponding neutronic behavior (i.e. particle transport solved by MCNP6.2) in the solution vessel. Thus, the development of a coupling methodology is vital to understand the system behavior at a variety of system design and postulated operating scenarios. In this study, we report on the k-effective (keff) calculation for the baseline solution vessel configuration with a selected solution concentration using MCNP K-code modeling. The associated correlation of thermal properties (e.g. density, viscosity, thermal conductivity, specific heat) at the selected solution concentration are developed based on existing experimental measurements in the open literature. The numerical coupling methodology between multiphase CFD and MCNP is successfully demonstrated, and the detailed coupling procedure is documented. In addition, improved coupling methods capturing realistic physics in the solution vessel thermal-neutronic dynamics are proposed and tested further (i.e. dynamic height adjustment, mull-cell approach). As a key outcome of the current study, a multi-physics coupling methodology between MCFD and MCNP is demonstrated and tested for four different operating conditions. Those different operating conditions are determined based on the neutron source strength at a fixed geometry condition. The steady state powers for the generic solution vessel at various operating conditions are reported, and a generalized correlation of the heat transfer coefficient for the current application is discussed. The assessment of multi-physics methodology and preliminary results from various coupled calculations (power prediction and heat transfer coefficient) can be further utilized for the system code validation and generic solution vessel design improvement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Final Project Report CFA-14-6357: A New Paradigm for Understanding Multiphase Ceramic Waste Form Performance

This project fabricated model multiphase ceramic waste forms with processing-controlled microstructures followed by advanced characterization with synchrotron and electron microscopy-based 3D tomography to provide elemental and chemical state-specific information resulting in compositional phase maps of ceramic composites. Details of 3D microstructural features were incorporated into computer-based simulations using durability data for individual constituent phases as inputs in order to predict the performance of multiphase waste forms with varying microstructure and phase connectivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗

Synthetic Absorption Lines from Simulations of Multiphase Gas in Galactic Winds

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high-resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T$ ~ 10 4 K supernova-driven wind and a cool $T$ ~ 10 7 K cloud of interstellar material, we create mock observables, including the optical depth ($\tau$) and covering fraction ($C$ $f$ ) of six commonly observed ions (Si II , C II , Si IV , C IV , N V , and O VI ) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth and overestimate the "true" value of $\tau$ for ions with low optical depth relative to the simulated data. In conclusion, the empirically estimated covering fractions match our direct calculations very well for the low-ionization ions; for the high-ionization ions, the empirical covering fractions underestimate the directly computed values by up to ~40%.

79 ASTRONOMY AND ASTROPHYSICS↗

Application of Nano-Indentation Test in Estimating Constituent Phase Properties for Microstructure-Based Modeling of Multiphase Steels

For multiphase advanced high strength steels (AHSS), the constituent phase properties play a crucial role in determining the overall mechanical behaviors. Therefore, it is important to accurately measure/estimate the constituent phase properties in the research of AHSS. In this study, a new nanoindentation-based inverse method that we developed was adopted in estimating the phase properties of a low alloy Quenching and Partitioning (Q&P) steel. A microstructure-based Finite Element (FE) model was also generated based on the Electron BackScatter Diffraction (EBSD) and Scanning Electron Microscopy (SEM) images of the Q&P steel. The phase properties estimated from nanoindentation were first compared with those estimated from in-situ High Energy X-Ray Diffraction (HEXRD) test and, then, employed in the generated FE model to examine whether they can be appropriately used as the input properties for the model. The results show that the estimated phase properties from the inverse method are similar to those from HEXRD, and that the Ultimate Tensile Strength (UTS) and Uniform Elongation (UE) predicted from the FE model based on the estimated phase properties are also similar to those of tensile experiment of the Q&P steel. Based on the results in this study, the nanoindentation-based inverse method appears to be a viable way in determining the phase properties of complex multiphase steels with submicrometer grain sizes.

Cheng, Guang↗

Simultaneous real-time measurement of composition, flow, attenuation, density, and pipe-wallthickness in multiphase fluids

Apparatus (10) and methods for making simultaneous measurements of composition (water-cut), fluid flow, and sound attenuation in a multiphase fluid flowing (15) through a pipe (12) in real-time, using the same apparatus (10) are described. Additionally, the apparatus (10) provides real-time pipe wall thickness monitoring for observing pipe corrosion or internal deposition. Knowledge of wall thickness is necessary to correct for water-cut (oil-water composition) automatically by adjusting the liquid path length internal to the pipe (spool). The use of short duration frequency chirp excitation signals (24) enables the apparatus to provide information that can be used to extract multiple levels of information from the same measurement in multiphase fluids including the presence of a significant quantity of gas (˜60% gas volume fraction) in different flow regimes. Besides measuring steady flow, this device is useful for measurements during fast changing flows, such as for a rod-pumped well. Measurements up to about 1000 times a second can reliably be made.

Sinha, Dipen N.↗

Evaluating the Practicability of Multiphase CFD for BWR Fuel Assembly Analysis

Multiphase computational fluid dynamics (M-CFD) models of a boiling water reactor (BWR) fuel assembly have been generated in the commercial CFD code, STAR-CCM+ Ver.11.04.010-R8. The Eulerian multiphase approach was employed, in combination with a baseline ‘Zero Closure model’ for two-phase flow boiling, applicable to high void fraction regimes. The capability and practicability of MCFD for BWR fuel assembly were evaluated against the international OECD/NRC BWR Full-size Finemesh Bundle Test (BFBT) benchmark data. The benchmark selected, BFBT-4101-61, is a steady-state full power case, with high void fraction conditions. Geometrical modeling and mesh construction of BWR assemblies introduce challenges in generating high quality computational meshes, while retaining a low cell count to allow acceptable runtimes. In order to evaluate the scalability of the CFD code, parallel CFD simulations were carried out with the High Performance Computing (HPC) system of the Idaho National Laboratory. M-CFD results display good agreement with the experimental data for both local void fraction distribution and exit quality. The parallel speedup factor is linearly proportional to the number of processors up to 500 cores, and retains over 92% efficiency with 1000 cores. The execution time of the simulation decreased exponentially with increasing number of processors while the cost of the simulation increased linearly. The proposed baseline closure, in combination with the parallel performance of the STAR-CCM+ code demonstrated promising capabilities for analyzing BWR fuel assemblies.

97 MATHEMATICS AND COMPUTING↗

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗

Computational Studies of Thermal Conductivity of Multiphase Ceramics for Inert Matrix Fuel

Nuclear energy is a major contributor to base load electrical power around the world and is predicted to continue to expand in the coming years to meet the increasing global demand for electricity. As new reactors are being constructed to meet these demands, it is imperative to make power plants intrinsically safer and more cost efficient. A major consideration to improving safety and performance is the enhancement of current ceramic nuclear fuels. By using an inert matrix fuel, additional phases can be added to the fuel to increase both the thermal and mechanical properties of the fuel at room and reactor conditions. In this study, multiphase ceramics and the individual single-phase components, will be investigated to understand the effects of multiple phases, various grain sizes, and interfaces on the temperature dependent thermal conductivity. A coupled approach involving both experimental measurements and computational modeling will give insight into the development of future models for understanding the role interfaces play in thermal transport within multiphase systems. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Chapter 10: Experimental Methods for the Determination of Parameters. In the book Multiphase Catalytic Reactors: Theory, Design, Manufacturing, and Applications.

The determination of kinetic parameters is a key element in catalyst research and multiphase catalytic process development. In this chapter, the experimental methods used to determine kinetic parameters, the nature of the parameters, and the areas of application and advantages of the different experimental approaches are discussed. Emphasis is placed on transient techniques and methods that provide intrinsic kinetic parameters, i.e., parameters that can be directly linked to the surface composition/structure of a heterogeneous catalyst.

02 PETROLEUM↗

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗