Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flow simulation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Open source implementation of glued sphere discrete element method and nonspherical biomass fast pyrolysis simulation

Here, in this research, a glued-sphere Discrete Element Method (DEM) method was implemented in the open-source, Computational Fluid Dynamics (CFD) software MFiX. The implementation was verified using a cylinder-wall collision and then validated by simulating the packing and fluidization of non-spherical particles. The validated code was applied to simulate fast pyrolysis of non-spherical biomass particles in a bubbling fluidized bed. The glued sphere occupancy ratio was proposed to quantify the quality of shape resolution using glued sphere. Shape resolution showed significant influence on the packing height in the simulation of particle packing and an occupancy ratio of 80 % was recommended. Its influence is minor in fully fluidized bed but can be 8 times higher in fixed fluidized bed. Three tested heat transfer models predicted similar yields of elongated biomass fast pyrolysis. The solver developed in this research can be used to simulate other multiphase reacting flows involving non-spherical particles.

42 ENGINEERING↗
DOE OSTI2021

A micromechanical hydro-mechanical-damage coupled model for layered rocks considering multi-scale structures

An anisotropic Hydro-Mechanical-Damage coupled model was established in the framework of thermodynamics for saturated layered rocks, which contain arbitrarily-distributed cracks of smaller scales and orientated bedding planes of much larger scales. Here, the proposed model considers the anisotropic damage growth and frictional sliding of cracks and shear sliding and dilatancy of bedding planes, as well as the interaction between cracks and bedding planes. Numerical implementation of the nonlinear coupled model was developed based on TOUGHREACT, which is a well-established simulator for multiphase fluid flow and reactive transport analysis. The proposed model was first validated with laboratory water injection test results on rock blocks, and then used to study the effects of injection rates, in situ stresses and bedding planes on fluid injection processes at the field scale, demonstrating the importance to consider the variation of multi-scale structures for better understanding the coupled Hydro-Mechanical responses in layered rocks.

58 GEOSCIENCES↗
DOE OSTI2021

An Investigation of the Effects of Volume Fraction on Drag Coefficient of Non-Spherical Particles Using PR-DNS

Prediction of the drag coefficient is required in gas-particle multiphase flow modeling and simulation. Experimental data and correlations on the fixed-bed system of spherical particles with high volume fractions for various possible arrangements are available in the literature. However, the effect of volume fraction on the drag coefficient of non-spherical particles is not well studied. In solving the momentum equation, the volume fraction plays a vital role in determining the flow resistances. In this paper, we study the impact of volume fraction in the range of 0.069 to 0.65 on the drag coefficient using the computational fluid dynamics (CFD) simulation of air for Reynold number in the range of 10 to 10000 using particle resolved direct numerical solution (PR-DNS). Regular non-spherical particles such as a cube, tetrahedron, and spheroids are used in this study since their single particle’s drag coefficient data are available in the literature for comparison. For this work, the simulations are carried out in the Ansys Fluent using polyhedral mesh, which consumes significantly less computational time and power. The study showed the sphericity and volume fraction have significant impact on the bed pressure drop and average drag coefficient of the particles in the bed especially in high Reynolds number regime. The bed of the spheroid experiences the lowest drag being the most streamlined particle, and the particles with the edges result in a large drag coefficient due to flow separation at the discontinuity. The vector plots verify this behavior where large wake regions are observed behind the tetrahedron particle.

Mahyawansi, Pratik↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. 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 relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
DOE OSTI2022

A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles

This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.

42 ENGINEERING↗
DOE OSTI2022

Geologic Carbon Storage of Anthropogenic CO2 under the Colorado Plateau in Emery County, Utah

Geologic Carbon Storage (GCS) is a promising technology for storing large volumes of anthropogenic CO2 effectively and permanently. Numerical simulations are an integral part of site selection and characterization for any potential GCS site. As part of the DOE-funded CarbonSAFE Rocky Mountains Phase I project, a regional GCS analysis was undertaken to understand the efficacy of storing CO2 emissions from the power generation and heavy industry in central Utah’s favorable geology. In this study, the injection of CO2 for geologic storage was simulated in the Navajo Sandstone Formation in Emery County, Utah. Carbon dioxide was sourced from regional power generation stations and heavy industries throughout Utah, with an emphasis on emissions reduction at the Hunter Power Plant near Castle Dale, Utah. A simulation grid was extracted from the project’s geological model encompassing an area around Price, Huntington, and Castle Dale in central Utah. The Navajo Sandstone Member of the Glen Canyon Group was the target of CO2 injection with the overlying Carmel formation providing the primary seal. A suite of simulations was performed assessing the viability of this area for permanent CO2 storage. Results indicate that the area can not only store 46 million metric tons of anthropogenic CO2, meeting the project goals, but this area has the capacity to securely store at least 1.3 billion tons of CO2, suggesting the injection site and surrounding geology are suitable locations for commercial-scale GCS.

58 GEOSCIENCES↗
DOE OSTI2021

Polysiloxane coatings on microspheres based on Multiphase Flow with Interface Exchange-Discrete Element Modelling

Polysiloxane coatings on yttria stabilized zirconia (YSZ) microspheres of 500 μm were simulated in a spout fluidized bed coater using Multiphase Flow with Interface Exchange-Discrete Element Modelling (MFiX-DEM). In this work, two different coater configurations were developed to study the influence of gas velocity and its distribution on particle dynamics. The presence of the Wurster tube not only enhances the distribution but also increases the overall residence time of the particles. Investigations were also carried out with different Wurster tube positions (normal, 10% and 20% lowered from its initial position). Among these, 20% lowered Wurster tube position demonstrated the most effective coating process. The effects of gas inlet pressure on the average gas velocity and the distribution of particles were analyzed. More than 97% of the particles can be retained. The derived results, including average gas velocity, particle retention percentage, and distribution of particles with gas velocity, are being used to guide the experimental work in obtaining defect-free coatings for YSZ microspheres.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GEOS: A performance portable multi-physics simulation framework for subsurface applications

GEOS is a simulation framework focused on solving tightly coupled multi-physics problems with an initial emphasis on subsurface reservoir applications. Currently, GEOS supports capabilities for studying carbon sequestration, geothermal energy, hydrogen storage, and related subsurface applications. The unique aspect of GEOS that differentiates it from existing reservoir simulators is the ability to simulate tightly coupled compositional flow, poromechanics, fault slip, fracture propagation, and thermal effects, etc. Extensive documentation is available on the GEOS documentation pages (GEOS Documentation, 2024). Note that GEOS, as presented here, is a complete rewrite of the previous incarnation of the GEOS referred to in (Settgast et al., 2017).

58 GEOSCIENCES↗
DOE OSTI2024

Modeling gas migration through clay-based buffer material using coupled multiphase fluid flow and geomechanics with stress-dependent gas permeability

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2022

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗
DOE OSTI2024

Model Development for Thermal-Hydrology Simulations of a Full-Scale Heater Experiment in Opalinus Clay

Disposal of commercial spent nuclear fuel in a geologic repository is studied. In situ heater experiments in underground research laboratories provide a realistic representation of subsurface behavior under disposal conditions. Here, this study describes process model development and modeling analysis for a full-scale heater experiment in opalinus clay host rock. The results of thermal-hydrology simulation, solving coupled nonisothermal multiphase flow, and comparison with experimental data are presented. The modeling results closely match the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Using AMReX's Embedded Boundaries to Support MFIX-Exa’s Geometry Capabilities

MFIX-Exa is an AMReX-based application used for the simulation of chemically reacting multiphase flows inside reactors. The reactor geometry is mostly a combination of primitive shapes such as cylinders, cubes and spheres of varying dimensions connected together with potentially a lot of dead space between them. Additionally, the reactors may have particle-laden jets introduced at different parts of their geometry to enhance reactant mixing. This talk presents how MFIX-Exa uses OpenSCAD’s csg format to allow users to define the embedded boundary (EB) geometry and how it supports specifying inflow through certain sections of the EB. The talk will describe how MFIX-Exa leverages the EB features of grid-pruning, restart files and redistribution. Finally, the talk goes over some capabilities which are currently not available in MFIX-Exa but planned for the future, such as moving parts of the reactor geometry.

embedded boundaries↗
DOE OSTI2021

Sensitivity Analysis of MFiX-PIC Parameters Using Nodeworks, PSUADE, and DAKOTA

The study presented in this report was aimed to demonstrate UQ analysis performed not only with Nodeworks, but also two other well-established UQ software tools from the U.S. DOE’s National Laboratories (PSUADE from Lawrence Livermore National Laboratory and DAKOTA from Sandia National Laboratory). It is important to emphasize that the motivation of this study was not to determine the best UQ software, but to verify if the global sensitivity analyses from the end-to-end workflow in Nodeworks are consistent with the results of other two UQ software. The components of Nodeworks from Python’s ecosystem have been tested as standalone libraries. However, an assessment study for the complete workflow targeting a specific UQ analysis has not been performed for Nodeworks. Hence, this study is expected to serve as an equivalent of solution verification for Nodeworks using other established UQ tools as reference solution. For this purpose, three distinct flow configurations (i.e., settling bed, bubbling fluidized, and circulating fluidized bed) have been used as representative multiphase flow problems of interest. The results of the systematic simulation campaigns performed in an earlier study using the particle-in-cell (PIC) approach in the Multiphase Flow with Interphase eXchanges (MFIX) suite of solvers (i.e., MFiX-PIC) was utilized. The same set of tabulated results was provided as input to the different UQ software for global sensitivity analysis. Results for the three cases indicate that based on the Sobol’ Sensitivity Indices method the order of importance ranking determined by Nodeworks for the Sobol’ Total Sensitivity Indices is consistent with PSUADE and DAKOTA in each case for the five model parameters considered. The input files for Nodeworks for the three cases are also shared through NETL’s Gitlab repository for the reader interested in reproducibility and further analysis (See Section 1.2).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H 2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories. This is particularly the case for bentonite, an often-preferred barrier material. Gas flow in preferential pathways is characterized by localization and spontaneous behavior, which is challenging to simulate in numerical models due to strong hydro-mechanical coupling. To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project, this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations. The proposed model utilized hydro-mechanical spatial distributions, namely Young's modulus and gas entry pressure, and elastoplasticity combined with a linear swelling model. A conceptual strain-dependent permeability approach simulated dilatancy-controlled gas flow based on hydro-mechanical coupling. To test the effectiveness of the presented approach, a gas injection test in a compacted, saturated bentonite sample was simulated using the open-source code OpenGeoSys 5.8 and compared with experimental observations. The presented methodology is capable of simulating localized gas flow in preferential pathways. The spatial distributions of Young's modulus and gas entry pressure affect the swelling pressure, relative permeability and, in combination with the strain-dependent permeability model, also the intrinsic permeability.

42 ENGINEERING↗
DOE OSTI2024

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↗
DOE OSTI2021

SiOC coatings on yttria stabilized zirconia microspheres using a fluidized bed coating process

In this study, defect-free SiOC coatings were prepared on yttria stabilized zirconia (YSZ) microspheres by a fluidized bed coating process. Effects of the rheological properties of the coating solution on the coating process were elucidated. An impact regime diagram was constructed, which demonstrated that the coating mechanisms were collision/impact. During the fluidized bed coating, longer spouted time resulted in wider dispersion, longer residence time, and more circulatory motion of particles; the fluid distributed more uniformly throughout the column, as demonstrated in our Multiphase Flow with Interface eXchange (MFiX) simulations. Two-step pyrolysis in Ar achieved complete coating layers, which were comprised of SiOC, SiO 2 , SiC, and graphite. The two stage mass loss during the pyrolysis corresponded to simultaneous reactions due to depolymerization and hydrocarbon loss from 400 to 600°C. Carbon cluster size in the pyrolyzed samples was calculated to be 25 ± 2 Å. In conclusion, this work provides a new method for producing SiOC coatings on micron spheres, with nuclear TRISO fuel particles as the application.

42 ENGINEERING↗
DOE OSTI2023

PFLOTRAN 5

PFLOTRAN leverages massively parallel, high performance computing to simulate large-scale non-isothermal multiphase flow, multicomponent reactive transport and electrical resistivity tomography (ERT) problems in the subsurface environment. Researchers have employed PFLOTRAN to simulate these Earth system processes on leadership class supercomputers for over two decades. The code is designed to predict the future estate of environmental systems and better inform stakeholders in the regulatory decision making process (e.g., fate of contaminants, long-term stewardship for nuclear waste, impact of climate change, etc.). A diverse team of scientists oversees PFLOTRAN development and maintenance under an open-source licensing agreement and manages contributions from an international community of researchers.

Hammond, Glenn↗