Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase fluid”

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 253 records · Page 14

Numerical Simulation of an Industrial Fluid Catalytic Cracking Regenerator

Fluid catalytic cracking (FCC) is one of the most important conversion processes in petroleum refineries, and the FCC regenerator is a key part of an FCC unit utilized in the recovery of solid catalyst reactivity by burning off the deposited coke on the catalyst surface. Furthermore, a three-dimensional multiphase, multispecies reacting flow computational fluid dynamics (CFD) model was established to simulate the flow and reactions inside an FCC regenerator. The Euler–Euler approach, where the two phases (gas and solid) are considered to be continuous and fully interpenetrating, is employed. The model includes gas–solid momentum exchange, gas–solid heat exchange, gas–solid mass exchange, and chemical reactions. Chemical reactions incorporated into the model simulate the combustion of coke which is present on the catalyst surface. We validated the simulation results using plant data.

42 ENGINEERING↗

EFRC-MUSE: Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials

Phase interactions and fluid properties in geological and other environments are critical in applications ranging from hydrogen production and geologic storage and recovery, carbon dioxide storage and sequestration, and the sustainable use of water resources. The four goals of EFRC-MUSE: Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials were based on the priority directions articulated in the Basic Research Needs documents, and the scientific needs in nanoscience. 1. Develop a fundamental understanding of confinement and surface interactions in mesoscale media with nanometer-sized pores on the phase behavior, thermodynamic and multiphase flow properties of multicomponent fluid mixtures. 2. Examine the impact of mineralogy and material heterogeneity on mechanical properties to better understand chemo-mechanical interactions in material failure. 3. Determine in-operando cross-scale structural and nanostructural material properties with fluids in confinement and flow under realistic condi

58 GEOSCIENCES↗

A multiphase electrokinetic flow model for electrolytes with liquid/liquid interfaces

A numerical model for electrokinetic flow of multiphase systems with deformable interfaces is presented, based on a combined level set-volume of fluid technique. A new feature is a multiphase formulation of the Nernst–Planck transport equation for advection, diffusion and conduction of individual charge carrier species that ensures their conservation in each fluid phase. The numerical model is validated against the analytical results of Zholkovskij et al. (2002) [1], and results for the problem of two drops coalescing in the presence of mobile charge carriers are presented. The time taken for two drops containing ions to coalesce decreases with increasing ion concentration.

42 ENGINEERING↗

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the Volume of Fluid (VOF) and Sharp Interface multiphase 2D axisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. Turbulence modeling; turbulence damping at the interface; and using constant vs. temperature dependent fluid properties were shown to have the most profound influence on predicted tank pressures and temperatures. Current study indicates that including phase change at the interface into the computational model is crucial for accurate prediction of the tank self-pressurization process. The effect of accommodation coefficient was also studied. Tank pressure values predicted by the VOF, and Sharp Interface models are within 4% of the experimental ones.

Computational Fluid Dynamics↗

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗

Simulation of Boiling Two-Phase Flow in a Helical Coil Steam Generator Using the Spectral Element Code Nek-2P

To enable the design of a light water small modular reactor, the boiling flow inside a helical coil steam generator has been simulated with the two-fluid model in Nek-2P. Nek-2P is the multiphase branch of the spectral element code Nek5000. Details of the implementation of the two-fluid model and the included closure models are discussed. The presented closure models include interactions for momentum, heat, and mass transfer between phases. Models for the drag, lift, and turbulent dispersion forces are included. The complete model is fully consistent in the limits of both phasic volume fractions approaching zero and is able to simulate flows of dispersed vapor, continuous liquid, dispersed liquid, continuous vapor, or any combination thereof. The closure models and their implementation in Nek-2P have been validated by comparing to experimental data for a boiling flow, demonstrating excellent agreement. Results from the simulation of the helical coil indicate strong phasic separation driven by the effects of buoyancy and inertia. Finally, significant differences were observed in the results compared to simulations performed using Star-CCM+, although these differences were somewhat expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupled thermal-hydrological-mechanical behavior of rock mass surrounding a high-temperature thermal energy storage cavern at shallow depth

The present study is aimed at numerically examining the thermal-hydrological-mechanical (THM) processes within the rock mass surrounding a cavern used for thermal energy storage (TES). We considered a cylindrical rock cavern with a height of 50 m and a radius of 10 m storing thermal energy of 350ºC as a conceptual TES model and simulated its operation for 30 years using THM coupled numerical modeling. At first, the insulator performance was not considered for the purpose of investigating the possible coupled THM behavior of the surrounding rock mass; then, the effects of an insulator were examined for different insulator thicknesses. The key concerns were focused on the hydro-thermal multiphase flow and heat transport in the rock mass around the thermal storage cavern, the effect of evaporation of rock mass, thermal impact on near the ground surface and the mechanical behavior of the surrounding rock mass. It is shown that the rock temperature around the cavern rapidly increased in the early stage and, consequently, evaporation of groundwater occurred, raising the fluid pressure. However, evaporation and multiphase flow did not have a significant effect on the heat transfer and mechanical behavior in spite of the high-temperature (350ºC) heat source. The simulations showed that large-scale heat flow around a cavern was expected to be conductiondominated for a reasonable value of rock mass permeability. Thermal expansion as a result of the heating of the rock mass from the storage cavern led to a ground surface uplift on the order of a few centimeters and to the development of tensile stress above the storage cavern, increasing the potentials for shear and tensile failures after a few years of the operation. Finally, the analysis showed that high tangential stress in proximity of the storage cavern can some shear failure and local damage, although large rock wall failure could likely be controlled with appropriate insulators and reinforcement.

58 GEOSCIENCES↗

Dense Multiphase Flow Simulation: Continuum Model for Poly-Dispersed Systems Using Kinetic Theory

The overall objective of the project was to verify the applicability of the FCMOM approach to the kinetic equations describing the particle flow dynamics. For monodispersed systems the fundamental equation governing the particle flow dynamics is the Boltzmann equation. During the project, the FCMOM was successfully applied to several homogeneous and in-homogeneous problems in different flow regimes, demonstrating that the FCMOM has the potential to be used to solve efficiently the Boltzmann equation. However, some relevant issues still need to be resolved, i.e. the homogeneous cooling problem (inelastic particles cases) and the transition between different regimes. In this report, the results obtained in homogeneous conditions are discussed first. Then a discussion of the validation results for in-homogeneous conditions is provided. And finally, a discussion will be provided about the transition between different regimes. Alongside the work on development of FCMOM approach studies were undertaken in order to provide insights into anisotropy or particles kinetics in riser hydrodynamics. This report includes results of studies of multiphase flow with unequal granular temperatures and analysis of momentum re-distribution in risers due to particle-particle and fluid-particle interactions. The study of multiphase flow with unequal granular temperatures entailed both simulation and experimental studies of two particles sizes in a riser and, a brief discussion of what was accomplished will be provided. And finally, a discussion of the analysis done on momentum re-distribution of gas-particles flow in risers will be provided. In particular a discussion of the remaining work needed in order to improve accuracy and predictability of riser hydrodynamics based on two-fluid models and how they can be used to model segregation in risers.

42 ENGINEERING↗

Towards a Modern Theory of Multiphase Filtration Flow

An alternative theoretical model of joint filtration flow of immiscible incompressible fluids is presented. The model takes into account relaxation processes due to the interchange of the fluids between pores of difference sizes which is driven by capillary forces. The fluids occupy connected regions in a four-dimensional space formed by three coordinates and the pore length scale. When the fluid exchange between pores of given sizes is effected by way of successive flow through pores of all the intermediate sizes, the pressure within each region is governed by a hyperbolic equation, the role of time being played by the pore linear scale. Pressure jumps across hypersurfaces separating the regions equal corresponding values of the capillary pressure. A supplementary condition at any such hypersurface requires the speed of its displacement in the four-dimensional space to coincide with the normal velocity components of both the adjoining fluids. As a result, a principally new statement of multiphase filtration flow problems is gained with allowance for capillary relaxation in the porous space.

Buyevich, Yu A.↗

CFD-DEM Onset of Motion Analysis for Application to Bed Scour Risk Assessment

This CFD study with DEM was done as a part of the Federal Highway Administration’s (FHWA’s) effort to improve scour design procedures. The Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) model, available in CD-Adapco’s StarCCM+ software, was used to simulate multiphase systems, mainly those which combine fluids and solids. In this method the motion of discrete solids is accounted for by DEM, which applies Newton's laws of motion to every particle. The flow of the fluid is determined by the local averaged Navier–Stokes equations that can be solved using the traditional CFD approach. The interactions between the fluid phase and solids phase are modeled by use of Newton's third law. The inter-particle contact forces are included in the equations of motion. Soft-particle formulation is used, which allows particles to overlap. In this study DEM was used to model separate sediment grains and spherical particles laying on the bed with the aim to analyze their movement due to flow conditions. Critical shear stress causing the incipient movement of the sediment was established and compared to the available experimental data. An example of scour around a cylindrical pier is considered. Various depths of the scoured bed and flow conditions were taken into account to gain a better understanding of the erosion forces existing around bridge foundations. The decay of these forces with increasing scour depth was quantified with a ‘decay function’, which shows that particles become increasingly less likely to be set in motion by flow forces as a scour hole increases in depth. Computational and experimental examples of the scoured bed around a cylindrical pier are presented.

97 MATHEMATICS AND COMPUTING↗

Experimental measurements of gas-particle flows in large-scale strippers

Measurements of gas–solid particle flow in a large-scale stripper unit are reported. The experiments target industrial (pilot) scale measurements of gas-particle flow that can be used to validate numerical methods for modeling multiphase flows, such as computational fluid dynamics coupled to the discrete element method (CFD-DEM). Specifically, experiments were performed with Geldart Group B (533 µm) glass beads in a stripper unit. In general, previous CFD-DEM validation datasets were performed in benchtop experimental systems and limited to O(10 5 ) particles. Here, in this work, the number of particles in the isolated section of the stripper is O(10 7 ) for all tests. Measurements are reported for a ~2 m tall section of the ~1 m diameter stripper. Radial mass flux profiles measured at five axial positions and three axial pressure drop measurements are reported for six operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Numerical Study of Two-Phase Flow and Interfacial Mass Transfer in a Wetted Wall Column for Counter-Current Gas Absorption

We present a numerical study to investigate the effects of surface waves of falling films on the mass transfer from flowing gas phase into the liquid film. The model employs the volume of fluid method to explicitly track the gas-liquid phase interface and a one-fluid formulation to model the multiphase flow coupled with mass transfer between the two phases. The surface waves are generated by the time-varying injection rate of the falling film. The effects of surface wave frequencies and amplitudes on the overall mass transfer also are investigated for a counter-current gas-liquid flow of oxygen and water. Because of the solubility difference in the two phases, the model can capture the concentration discontinuity across the liquid-gas interface. Good agreement is obtained between the numerical results and experimental data reported in the literature. By systematically varying the frequencies and amplitudes of flow rate for the falling film, the surface waves generated on the film surface enhance mass transfer in general. Here, the model prediction shows a proportional increase in mass transfer for small frequencies of surface waves, but saturation at higher frequencies.

42 ENGINEERING↗

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↗

A numerical comparison of precipitating turbulent flows between large‐eddy simulation and one‐dimensional turbulence

This study presents the results of computational fluid dynamics simulations of a multiphase, reacting, turbulent mixing layer in an idealized geometry. The purpose is to compare large‐eddy simulation (LES) to one‐dimensional turbulence (ODT) and examine the trends of the flow under differing mixing conditions. Aqueous streams are mixed together to precipitate polymorphs of calcium carbonate. The polymorphs of calcium carbonate are tracked numerically using population balance equations (PBE). Each PBE contains all of the relevant physical models to track the particle evolution including nucleation, growth, and aggregation. A simple subgrid mixing model that is convenient for use with PBEs was implemented in the LES code. The higher spatial resolution achievable with ODT allowed an investigation on the effect of resolution on the mixing‐model error. The Reynolds number of the flow is varied and is shown to cause a decrease in average particle sizes with higher mixing rates. © 2015 American Institute of Chemical Engineers AIChE J , 61: 3185–3197, 2015

Abboud, Alex W.↗

Unveiling the Essential Parameters Driving Mineral Reactions during CO2 Storage in Carbonate Aquifers through Proxy Models

Numerical simulation is a commonly employed technique for studying carbon dioxide (CO2) storage processes in porous media, particularly saline aquifers. It enables the representation of diverse trapping mechanisms and the assessment of CO2 retention capacity within the subsurface. The intricate physicochemical phenomena involved necessitate the incorporation of multiphase flow, accurate depiction of fluid and rock properties, and their interactions. Among these factors, geochemical reaction rates and mechanisms are pivotal for successful CO2 trapping in carbonate reactive rocks. However, research on kinetic parameters and the influence of lithology on CO2 storage remains limited. This limitation is partly due to the challenges faced in laboratory experiments, where the time scale of the reactions and the lack of in situ conditions hinder accurate measurement of mineral reaction rates. This study employs proxy models constructed using response surfaces calibrated with simulation results to address uncertainties associated with geochemical reactions. Monte Carlo simulation is utilized to explore a broader range of parameters and identify influential factors affecting CO2 mineralization. The findings indicate that an open database containing kinetic parameters can support uncertainty assessment. Additionally, the proxy models effectively represent objective functions related to CO2 injectivity and mineralization, with calcite dissolution playing a predominant role. pH, calcite concentration, and CO2 injection rate significantly impact dolomite precipitation, while quartz content remains unaffected.

58 GEOSCIENCES↗

Simulation of gravity driven free-surface flow in fractured geological media

Abstract—Simulation of unsaturated free-surface flow in fractured geological media represents a challenge due to the highly heterogeneous flow field induced by extensive faults, joints and fissures. Free-surface flow in unsaturated media leads to highly intermittent flow regimes and flow velocities well above those assumed for the bulk volume. However, common modeling approaches relying on volume-averaged effective equations fail to capture this flow feature. In this work we present micro-scale flow simulations using a three-dimensional multiphase SPH code. Pairwise fluid-fluid and solid-fluid interaction forces are used to simulate a wide range of wetting conditions encountered on rock surfaces. It is shown that static contact angles for sessiles droplets are independent of the model discretization, i.e. the total amount of particles. Thus, computation times can be reduced without sacrificing qualitative or quantitative information. Furthermore we show that our model is in accordance with general scaling laws for droplet flow.

Kordilla, Jannes↗

Density and Cavitating Flow Results from a Full-Scale Optical Multiphase Cryogenic Flowmeter

Liquid propulsion systems are hampered by poor flow measurements. The measurement of flow directly impacts safe motor operations, performance parameters as well as providing feedback from ground testing and developmental work. NASA Marshall Space Flight Center, in an effort to improve propulsion sensor technology, has developed an all optical flow meter that directly measures the density of the fluid. The full-scale sensor was tested in a transient, multiphase liquid nitrogen fluid environment. Comparison with traditional density models shows excellent agreement with fluid density with an error of approximately 0.8%. Further evaluation shows the sensor is able to detect cavitation or bubbles in the flow stream and separate out their resulting effects in fluid density.

Korman, Valentin↗

Device-scale CFD study for mass transfer coefficient and effective mass transfer area in packed column

Carbon dioxide (CO 2 ) is a significant contributor of greenhouse gases emissions responsible for global warming. The fossil-fueled power plant is one of the major sources of CO 2 emission worldwide and needs to be mitigated. Post-combustion carbon capture by chemical absorption can be a promising technology to reduce CO 2 emission from power plants. The absorption is carried out via countercurrent flow in the packed columns. Although monoethanolamine (MEA) is a popular solvent, it has also limitations such as high volatility of solvent, corrosivity that leads to product degradation and damage to apparatuses, low CO 2 loading capacity, etc. Therefore, alternate solvents are being explored as prospective materials for carbon capture. A water-lean solvent, CO 2 binding organic liquids (CO2BOL), is currently being developed as a prospective material in carbon capture at the Pacific Northwest National Laboratory. Since the development of the material is at the primitive stage, it needs extensive investigation before industrial deployments. Hydrodynamics and mass transfer phenomena in the packed column using CO 2 BOL is unknown. Computational fluid dynamics (CFD) can be a useful tool to understand these phenomena in the packed column. Before employing the extensive CFD investigations, one needs to carefully validate CFD predictions for accurate flow predictions. This report demonstrates the validation of multiphase flow simulations against the experimental results for both structured and random packed columns. The reactive multiphase flow simulations using volume of fluid method were conducted at a wide range of liquid loads and contacts angles. The effective area is one of the critical factors dictating interphase mass transfer due to chemical reactions. Accordingly, CFD-predicted effective areas were compared with the experimental results of Tsai et al. [1] for Mellapak 250.Y structured packing. Both results matched well at different liquid loads in the presence of fast chemical kinetics using caustic solvents. Further, reactive multiphase flow simulations were also conducted in pall ring random packing to compute effective areas in the presence of fast chemical reactions. The CFD-predicted effective area compared well with correlation developed by Song et al. [2]. With successful validations, the current simulation framework can be further used for prediction of hydrodynamics and mass transfer phenomena using CO 2 BOL solvents in the packed column. Subsequently, it would be helpful in design optimization of a packed column. In addition to validation of the CFD model, this report also explains the mechanism for enhanced liquid side mass transfer for a structured packed column as compared to the wetted wall column.

20 FOSSIL-FUELED POWER PLANTS↗