Microwave Doppler Sensing of Particulate Flow in a Chemical Looping Reactor
Presented by NETL at 2019 NETL Multiphase Flow Science Workshop, 08/19
Chorpening, Benjamin↗
SEARCH · Engineering Papers
Results for “Multiphase flows”
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.
Burial-driven methane recycling in marine gas hydrate systems
Natural gas hydrate may be buried with sediments until it is no longer stable at a given pressure and temperature, resulting in conversion of hydrate into free gas. This gas may migrate upward and recycle back into the hydrate stability zone to form hydrate. As of yet, however, no quantitative description of the methane recycling process has been developed using multiphase flow simulations to model burial-driven gas hydrate recycling. In this study, we present a series of 1D multiphase transport simulations to investigate the methane recycling process in detail. By invoking the effects of capillary phenomena on hydrate and gas formation in pores of varying size, we find that a free gas phase can migrate a significant distance above the bulk base of hydrate stability. Since the top of the free gas occurrence is often identified as the base of the hydrate stability zone from seismic data, our results demonstrate that not only could this assumption mischaracterize a hydrate system, but that under recycling conditions the highest hydrate saturations can occur beneath the top of the free gas occurrence. We show that the presence of pore size distributions requires a replacement zone through which hydrate saturations progressively decrease with depth and are replaced with free gas. This replacement zone works to buffer against significant gas buildup that could lead to fracturing of overlying sediments. In conclusion, this work provides a framework for simulating flow and transport of methane within the 3-phase stability zone from a mass conservation perspective.
58 GEOSCIENCES↗
Development of a Spouted Bed Reactor for Chemical Looping Combustion
We have investigated a novel gas/solid contacting configuration for chemical looping combustion (CLC) with potential operating benefits. CLC configurations are typically able to achieve high fuel conversion efficiencies at the expense of high operating costs and low system reliability. The spouted fluid bed (SB) was identified as an improved reactor configuration for CLC, since it typically exhibits high heat transfer rates and offers the ability to use lower gas flows for material movement compared to bubbling beds (BB). Multiphase Flow with Interphase eXchanges (MFIX) software was used to establish a spouted fluid bed reactor design. An experimental setup was built to supplement the model. The experimental setup was also modified for testing under high temperature, reacting conditions (1073–1273 K). The setup was operated in either a spouted fluid bed or a bubbling bed regime to compare the performance attributes of each. Results for the reactor configurations are presented for CLC using a mixture of carbon monoxide and hydrogen as fuel. Compared to the bubbling bed, the spouted fluid bed reactor achieved an equivalent or better fuel conversion at a lower pressure drop over the material bed. The spouted fluid bed design represents a viable configuration to improve gas/solid contacting for efficient fuel conversion, lower energy requirements for material movement and increase operational robustness for CLC. The research laid the groundwork for future research into a multi-phase reacting flow CLC system. The system will be developed from computational fluid dynamic modeling and pilot-scale testing to expedite the development of CLC technologies.
Energy & Fuels↗
Microgravity Combustion Science and Fluid Physics Experiments and Facilities for the ISS
At the NASA Glenn Research Center, the Microgravity Science Program supports both ground-based and flight experiment research in the disciplines of Combustion Science and Fluid Physics. Combustion Science research includes the areas of gas jet diffusion flames, laminar flames, burning of droplets and misting fuels, solids and materials flammability, fire and fire suppressants, turbulent combustion, reaction kinetics, materials synthesis, and other combustion systems. The Fluid Physics discipline includes the areas of complex fluids (colloids, gels, foams, magneto-rheological fluids, non-Newtonian fluids, suspensions, granular materials), dynamics and instabilities (bubble and drop dynamics, magneto/electrohydrodynamics, electrochemical transport, geophysical flows), interfacial phenomena (wetting, capillarity, contact line hydrodynamics), and multiphase flows and phase changes (boiling and condensation, heat transfer, flow instabilities). A specialized International Space Station (ISS) facility that provides sophisticated research capabilities for these disciplines is the Fluids and Combustion Facility (FCF). The FCF consists of the Combustion Integrated Rack (CIR), the Fluids Integrated Rack (FIR) and the Shared Accommodations Rack and is designed to accomplish a large number of science investigations over the life of the ISS. The modular, multiuser facility is designed to optimize the science return within the available resources of on-orbit power, uplink/downlink capacity, crew time, upmass/downmass, volume, etc. A suite of diagnostics capabilities, with emphasis on optical techniques, will be provided to complement the capabilities of the subsystem multiuser or principal investigator-specific experiment modules. The paper will discuss the systems concept, technical capabilities, functionality, and the initial science investigations in each discipline.
Lauver, Richard W.↗
Development and Evaluation of a General Drag Model for Gas-Solid Flows via Deep Learning
This project presents the development and evaluation of a general drag model for gas–solid multiphase flows using deep learning techniques. A comprehensive database of more than 4,000 experimental and numerical data points for spherical and non spherical particles was compiled, incorporating geometric features such as sphericity, aspect ratio, and orientation. Several predictive approaches—including traditional em pirical correlations, machine learning, and deep neural networks—were benchmarked, with the proposed Drag Coefficient Correlation-aided Deep Neural Network (DCC DNN) demonstrating superior accuracy. To account for particle–particle interactions, additional drag data were generated using CFD-based simulations of packed and flu idized beds, leading to the development of a retrained model capable of incorporat ing volume fraction effects. Integration of the trained model with the MFiX CFD solver was achieved using FTorch, enabling drag predictions during discrete element method (DEM) simulations. Validation against experimental data for single particles and fluidized beds confirmed the model’s improved predictive ability, particularly for non-spherical geometries. While the model performed strongly under fluidized con ditions, limitations remained in unfluidized regimes, suggesting a need for expanded datasets. Overall, this study demonstrates the feasibility of combining deep learning with physics-informed CFD to improve drag modeling for gas–solid flows, with promis ing implications for scaling multiphase simulations in industrial applications.
42 ENGINEERING↗
Time-resolved 3D imaging of two-phase fluid flow inside a steel fuel injector using synchrotron X-ray tomography
The multiphase flow inside a diesel injection nozzle is imaged using synchrotron X-rays from the Advanced Photon Source at Argonne National Laboratory. Through acquisitions performed at several viewing angles and subsequent tomographic reconstruction, in-situ 3D visualization is achieved for the first time inside a steel injector at engine-like operating conditions. The morphology of the internal flow reveals strong flow separation and vapor-filled cavities (cavitation), the degree of which correlates with the nozzle's asymmetric inlet corner profile. Micron-scale surface features, which are artifacts of manufacturing, are shown to influence the morphology of the resulting liquid-gas interface. The data obtained at 0.1 ms time resolution exposes transient flow features and the flow development timescales are shown to be correlated with in-situ imaging of the fuel injector's hydraulically-actuated valve (needle). As more than 98.5% of the X-ray photon flux is attenuated within the steel injector body itself, we are posed with a unique challenge for imaging the flow within. Time-resolved imaging under these low-light conditions is achieved by exploiting both the refractive and absorptive properties of X-ray photons. The data-processing strategy converted these images with a signal-to-noise ratio of similar to 10 into a meaningful dataset for understanding internal flow and cavitation in a nozzle of diameter 200 mu m enclosed within 1-2 millimeters of steel.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A coupled, multiphase heat flux boundary condition for the discrete element method
Similar to single-phase flows, multiphase systems with a constant heat flux at the wall are common in practice. Unlike single-phase flows, however, the numerical implementation of a constant boundary flux is non-trivial due to the coupling between phases – i.e., the partition of the total flux between the phases can vary in space and time. A numerical technique for modeling such a boundary condition is proposed here and verified via simulations of gas-solid flows. While discrete-particle simulations of monodisperse particles are considered in this report, the technique can be extended to include radiation, polydisperse systems, and/or a continuum representation of the solids phase.
42 ENGINEERING↗
The Effect of Wettability Heterogeneity on Relative Permeability of Two-Phase Flow in Porous Media: A Lattice Boltzmann Study
Relative permeability is a critical parameter characterizing multiphase flow in porous media and it is strongly dependent on the wettability. In many situations, the porous media are nonuniformly wet. In this study, to investigate the effect of wettability heterogeneity on relative permeability of two-phase flow in porous media, a multi-relaxation-time color-gradient lattice Boltzmann model is adopted to simulate oil/water two-phase flow in porous media with different oil-wet solid fractions. For the water phase, when the water saturation is high, the relative permeability of water increases with the increase of oil-wet solid fraction under a constant water saturation. However, as the water saturation decreases to an intermediate value (about 0.4–0.7), the relative permeability of water in fractionally wet porous media could be lower than that in purely water-wet porous media, meaning additional flow resistance exists in the fractionally wet porous media. For the oil phase, similar phenomenon is observed. This phenomenon is mainly caused by the wettability-related microscale fluid distribution. According to both our simulation results and theoretical analysis, it is found that the relative permeability of two-phase flow in porous media is strongly related to three parameters: the fluid saturation, the specific interfacial length of fluid, and the fluid tortuosity in the flow direction. Lastly, the relationship between the relative permeability and these parameters under different capillary numbers is explored in this paper.
58 GEOSCIENCES↗
Synchrotron X-Ray Interrogation of Turbulent Gas–Liquid Mixing in Cryogenic Rocket Sprays
The atomization and vaporization of liquid jets within turbulent gaseous flows are characterized by the mixing phenomena occurring over a wide range of spatiotemporal scales. This creates a complex, turbid medium that is not easily interrogated using conventional optical-measurement techniques. In the current study, the optically dense, multiphase flow created by a cryogenic liquid jet injected into a turbulent gaseous coflow is probed using high-speed (MHz) X-ray radiography from a focused, narrowband synchrotron source to resolve the internal cascade of scales and the evolution to isotropic, homogeneous turbulence. Changes in the spectral characteristics for different flow conditions are furthermore correlated with changes in the spatial distributions of the liquid and gas phases within the spray using simultaneous X-ray radiography and tracer (krypton) fluorescence. It is found that an increase in entrainment and mixing infers an evolution in spectral characteristics toward the well-known -5/3 law of energy dissipation in the context of the classical Kolmogorov theory. Finally, these data demonstrate the utility of the synchrotron-based X-ray radiography and fluorescence for uncovering the internal, turbulent mixing processes in multiphase and optically dense flows.
42 ENGINEERING↗
An integrated approach to derive relative permeability from capillary pressure
Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.
58 GEOSCIENCES↗
Role of corner interfacial area in uniqueness of capillary pressure-satuation - interfacial area relation under transient conditions
Capillary pressure ( P c ) and phase saturation ( S w ) in two-phase flow are well known to be hysteretically related. Thermodynamically-derived multiphase flow theories conjecture that this hysteresis will be lifted if specific interfacial area ( a nw ) is included as a new state variable to create a unique P c -S-a nw surface. Specific interfacial area is defined as the total interfacial area per volume of a porous medium. Several studies have confirmed the existence of a unique P c -S w -a nw surface under equilibrium conditions for a given porous medium. However, there is only one experimental work in the literature, where the unique- ness of this surface under transient conditions was questioned. However, in the data analysis only the terminal menisci were considered to calculate the specific interfacial area. In this paper, we investigate the uniqueness of P c -S-a nw surfaces with and without the inclusion of corner fluid-fluid interfacial area, under different dynamic conditions, in two different micro-models bearing two different (granular vs. triangulated) pore morphologies. We establish a systematic metric to analyze hysteresis under different hydrodynamic conditions.
Godinez-Brizuela, Omar E.↗
Final Technical Report: Using Solid Particles as Heat Transfer Fluid for use in Concentrating Solar Power (CSP) Plants
In today’s industrial economy, energy consumption has never been higher. Over the last 15 years the US alone has consumed an average of nearly 100 quadrillion BTUs per year [21]. A need for clean and renewable energy sources has become quite apparent. The SunShot Initiative is an ambitious effort taken on by the United States Department of Energy that targets the development of solar energy that is cost-competitive with other methods for generating electricity. Specifically, this work is concerned with the development of concentrating solar power plants (CSPs) with granular media as the heat transfer fluid (HTF) from the solar receiver. Unfortunately, the prediction of heat transfer in multiphase flows is not well understood. For this reason, our aim is to fundamentally advance the understanding of multiphase heat transfer, particularly in gas-solid flows, while providing quantitative input for the design of a near black body receiver (NBB) that uses solid grains (like sand) as the HTF. Over the course of this three-year project, a wide variety of contributions have been made to advance the state-of-the art description for non-radiative heat transfer in dense, gas-solid systems. Comparisons between a state-of-the-art continuum heat transfer model and discrete element method (DEM) simulations have been drawn. The results of these comparisons brought to light the limitations of the continuum model due to inherent assumptions in its derivation. A new continuum model was then developed for heat transfer at a solid boundary by rigorously accounting for the most dominant non-radiative heat transfer mechanism (particle-fluid-wall conduction). The new model is shown to be in excellent agreement with DEM data and captures the dependence of heat transfer on particle size, a dependency that previous continuum models were not capable of. DEM and the new continuum model were then employed to model heat transfer in a variety of receiver geometries. The results provided crucial feedback on the efficiency and feasibility of various designs. Namely, a prototype design consisting of an array of heated hexagonal tubes was later supplanted by a vertical conduit with internal baffles. Due to low solids heat transfer on the bottom faces of the hexagonal tubes in the prototype, the predicted wall temperature gradients exceeded the design limitations. By contrast, the vertical conduit can be constructed to continually force particle-wall contacts, and thus, result in more desirable solids heat transfer and wall temperature gradients. Finally, a new heat flux boundary condition was developed for DEM simulations to assess the aforementioned wall temperature gradients. The new boundary condition advances current state-of-the-art techniques by allowing the heat fluxes to each phase to vary with space and time while the total flux remains constant. Simulations with the new boundary condition show that the total boundary heat flux is in good agreement with the imposed total boundary heat flux. While the methods we have utilized here are primarily numerical and fundamental by nature, they offer some key advantages of: (i) being robust and valid over a large range of conditions, (ii) able to quickly explore large parameter spaces, and (iii) aid in the construction of experiments. We have ultimately leveraged our computational capabilities to provide feedback on the design of a CSP which possesses great potential to become a cost effective source of clean and renewable electricity. Overall, ensuring that future energy demands are met in a responsible and efficient manner has far reaching impacts that span both ecologic and economic concerns. Regarding logistics, the project was successfully re-negotiated after the go/no-decisions of Years 1 and 2. All milestones were successfully completed.
14 SOLAR ENERGY↗
Determination of Flow Patterns by Image Analysis of a Rectangular Spouted Bed
Presented by NETL at 2017 NETL Multiphase Flow Workshop, 08/17
Yang, Jingsi↗
Modeling of dense reactive granular flows
Because of the lack of fundamental knowledge of both the fluid mechanical and thermal behavior of multiphase flow systems, the design of these processes has been necessarily based on empirical correlations and experiments performed in laboratory or pilot scale units. Detailed models and simulations of these systems would allow for the optimization of the processes and their scaling, without the need for expensive and time consuming testing. In this paper, we present the proposed hydrodynamic model in detail.
granular↗
Measurements of Shear Lift Force on a Bubble in Channel Flow in Microgravity
Under microgravity conditions, the shear lift force acting on bubbles, droplets or solid particles in multiphase flows becomes important because under normal gravity, this hydrodynamic force is masked by buoyancy. This force plays an important role in furnishing the detachment process of bubbles in a setting where a bubble suspension is needed in microgravity. In this work, measurements of the shear lift force acting on a bubble in channel flow are performed. The shear lift force is deduced from the bubble kinematics using scaling and then compared with predictions from models in literature that address different asymptotic and numerical solutions. Basic trajectory calculations are then performed and the results are compared with experimental data of position of the bubble in the channel. A direct comparison of the lateral velocity of the bubbles is also made with the lateral velocity prediction from investigators, whose work addressed the shear lift on a sphere in different two-dimensional shear flows including Poiseuille flow.
Nahra, Henry K.↗
Assessment of CO 2 Storage Potential in Naturally Fractured Reservoirs With Dual–Porosity Models
Naturally Fractured Reservoirs (NFR’s) have received little attention as potential CO 2 storage sites. Two main facts deter from storage projects in fractured reservoirs: (1) CO 2 tends to be nonwetting in target formations and capillary forces will keep CO 2 in the fractures, which typically have low pore volume; and (2) the high conductivity of the fractures may lead to increased spatial spreading of the CO 2 plume. Numerical simulations are a powerful tool to understand the physics behind brine-CO 2 flow in NFR’s. Dual-porosity models are typically used to simulate multiphase flow in fractured formations. However, existing dual-porosity models are based on crude approximations of the matrix-fracture fluid transfer processes and often fail to capture the dynamics of fluid exchange accurately. Therefore, more accurate transfer functions are needed in order to evaluate the CO 2 transfer to the matrix. This work presents an assessment of CO 2 storage potential in NFR’s using dual-porosity models. We investigate the impact of a system of fractures on storage in a saline aquifer, by analyzing the time scales of brine drainage by CO 2 in the matrix blocks and the maximum CO 2 that can be stored in the rock matrix. A new model to estimate drainage time scales is developed and used in a transfer function for dual-porosity simulations. Furthermore, we then analyze how injection rates should be limited in order to avoid early spill of CO 2 (lost control of the plume) on a conceptual anticline model. Numerical simulations on the anticline show that naturally fractured reservoirs may be used to store CO 2 .
54 ENVIRONMENTAL SCIENCES↗
Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation
We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.
54 ENVIRONMENTAL SCIENCES↗
MFIX DEM Enhancement for Industry-Relevant Flows (Final Report)
The overall goal of this two-phase project is to implement performance improvements of the Multiphase Flow with Interphase Exchanges (MFIX) Discrete Element Model (DEM) code that enable a transformative shift for industrial use. Prior to this effort, the largest simulations performed using MFIX are O(10 7 ) particles. This falls short of the O(10 9 ) particle simulations that must be completed on a timescale of days or weeks (vs. months or years) to enable simulations with physically-relevant domain sizes to be incorporated into industrial design cycles within five years. This was accomplished by tailoring best-in-class practices to bear on the unique challenges posed by the MFIX-DEM algorithm and code base. Scientific simulations (e.g., in cosmology, turbulent combustion) routinely use massively parallel computing to update far more particles in short wall clock times. Results from Phase 1 (1.5 years in duration) indicated significant gains in speed were possible for a wide range of benchmark cases. Moreover, a survey sent to >35 companies indicates that the timing is ideal for such an enhanced tool, with >80% of the respondents indicating that DEM is already value-added or will be within the next 5 years, and >70% of the respondents indicating that improved speed is the top computational priority. In Phase 2 (3.5 years in duration), the two major barriers that hinder industry from effectively using multiphase Computational Fluid Dynamics (CFD) to cut costs and improve performance, namely computational overhead and confidence in predictions, continued to be addressed. Regarding the former, the results from Phase 1 to guide the effort, with enhancements focused on an improved time-stepping algorithm and particle sorting. Four target problems of 1 billion particles each and increasing complexity were identified: homogeneous cooling, tumbler with continuous particle size distribution, discharge from a rectangular hopper and a cylindrical riser. Each of these were successfully simulated for relevant time scales (on order of seconds) using less than 24 hours of wall clock time. These represent the first 1-billion particle DEM simulations performed with MFIX, namely using the MFIX-Exa code. This code is currently under development at NETL in collaboration with Lawrence Berkeley National Laboratory. Regarding the second barrier on predictive uncertainty, experiments from Phase 1 (interacting nozzles - hydrodynamics only) and Phase 2 (very small-scale segregation experiments) were used to demonstrate the ability of two simplified approaches to uncertainty quantification (UQ). By limiting the number of particles, UQ based on the simplified treatment was compared to standard UQ, which was shown to have much higher computational demands. Experiments were also performed on a pilot-scale stripper unit to provide validation data for future CFD-DEM simulations and UQ.
20 FOSSIL-FUELED POWER PLANTS↗