Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

An overview of TOUGH-based geomechanics models

After the initial development of the first TOUGH-based geomechanics model 15 years ago based on linking TOUGH2 multiphase flow simulator to the FLAC3D geomechanics simulator, at least 15 additional TOUGH-based geomechanics models have appeared in the literature. This development has been fueled by a growing demand and interest for modeling coupled multiphase flow and geomechanical processes related to a number of geoengineering applications, such as in geologic CO 2 sequestration, enhanced geothermal systems, unconventional hydrocarbon production, and most recently, related to reservoir stimulation and injection-induced seismicity. This paper provides a short overview of these TOUGH-based geomechanics models, focusing on some of the most frequently applied to a diverse set of problems associated with geomechanics and its couplings to hydraulic, thermal and chemical processes.

58 GEOSCIENCES↗
DOE OSTI2018

Beyond the standard two-film theory: Computational fluid dynamics simulations for carbon dioxide capture in a wetted wall column

The standard two-film theory (STFT) is a diffusion-based mechanism that can be used to describe gas mass transfer across liquid film. Fundamental assumptions of the STFT impose serious limitations on its ability to predict mass transfer coefficients. To better understand gas absorption across liquid film in practical situations, a multiphase computational fluid dynamics (CFD) model fully equipped with mass transport and chemistry capabilities has been developed for solvent-based carbon dioxide (CO2) capture to predict the CO2 mass transfer coefficient in a wetted wall column. The hydrodynamics is modeled using a volume of fluid method, and the diffusive and reactive mass transfer between the two phases is modeled by adopting a one-fluid formulation. We demonstrate that the proposed CFD model can naturally account for the influence of many important factors on the overall mass transfer that cannot be quantitatively explained by the STFT, such as the local variation in fluid velocities and properties, flow instabilities, and complex geometries. The CFD model also can predict the local mass transfer coefficient variation along the column height, which the STFT typically does not consider.

Wang, Chao↗
DOE OSTI2021

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗
DOE OSTI2019

Spatial characterization of heterogeneous nanopore surfaces from XCT scans of Niobrara shale

The multiphase fluid flow through natural and geo-architected nanopore structures is fundamentally controlled by the geometry of such nanopore structures, pore wettability and interfacial tension. Furthermore, multiphase fluid flow through the nanopore structures is compromised when the pore surfaces are heterogeneous, having surface regions of different composition, polarity, and wetting characteristics. In this regard, compositional analysis of the pore channel surface has been a big challenge in the past, since 3D analysis is required. Currently, due to advances in nano X-ray CT, with a voxel resolution of 20 nm, it is possible to identify nanopore network structures and quantify their surface composition in 3D. In this paper, the development of procedures for detailed analysis to quantify the composition of the pore channel surface regions, including the corresponding area, is described using nanoCT images. Image analysis procedures, including preprocessing of the 3D tomographic images, correction for the partial volume effect (PVE), and the use of a new algorithm for determination of pore channel surface composition and area are reported on. Videos of the nanopore structure have been provided as supplementary information to further describe the heterogeneity of the nanopore surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Impact of solutes on the lattice parameters and elastic stiffness coefficients of hcp Fe from first-principles calculations

The hexagonal close-packed (hcp) -martensite phase in steels nucleates from the -austenite parent phase and can undergo further transformation to the -martensite phase or exist as a metastable phase depending on temperature, mechanical loading, and alloy chemistry. The solute-dependent lattice parameters and elastic stiffness coefficients of hcp Fe influence the mechanical properties of steels containing the -martensite phase, as well as the martensitic transformations between the phases. We use density functional theory to calculate the lattice parameters and of single-crystal hcp Fe as functions of solute concentration in the dilute limit for the substitutional solutes Al, B, Cu, Mn, and Si, and the octahedral interstitial solutes C and N. Our computationally efficient methodology separates the solute dependence of the into lattice strain and chemical bonding contributions. The computed data can be used to estimate the effect of solutes on polycrystalline elastic moduli and the strain energy associated with martensitic transformations. The data can also serve as inputs to microstructure-based models of multiphase steels containing the -martensite phase.

Fellinger, Michael R.↗
DOE OSTI2022

Modeling gas breakthrough and flow phenomena through engineered barrier systems using a discrete fracture approach

Compacted clays are being considered to build engineered barrier systems (EBS) intended for the safe isolation of high-level nuclear waste (HLW) and spent nuclear fuel (SNF). The corrosion of the metallic canister containing the HLW/SNF will lead to the generation and buildup of the gas pressure in the more internal part of the clay buffer. This phenomenon would eventually trigger the formation and propagation of fractures in the clay barrier, jeopardizing its safety functions. Here, in this work, we propose to use the fragmentation technique (MFT) to model evolving fractures in clays triggered by gas pressurization. The MFT has been successfully used to model the formation of fractures in concrete, drying cracks in soil, hydraulic and thermo-fractures in rocks. In this work, we extend the MFT to deal with multiphase fluid flow in deformable porous media, and we upgraded a fully coupled computer finite element code using the extended technique. The proposed approach is first verified against analytical solutions and is then applied to model gas breakthrough experiments in clays. A very satisfactory performance of the method is observed in all the analyses, showing the potential of the MFT to tackle multiphase flow problems in deformable porous media with evolving discontinuities.

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

High-speed characterization of two-phase flow and bubble dynamics in titanium felt porous media for hydrogen production

Bubble dynamics and two-phase flow phenomena are closely related to the performance of proton exchange membrane electrolyzer cells (PEMECs). This paper reports an in-situ study of the oxygen bubble behavior and associated multiphase evolutions in the anode side of PEMECs with titanium (Ti) felt liquid gas diffusion layers (LGDLs) by a high-speed visualization system. The micro oxygen bubble dynamics was captured and analyzed at different locations and virous operating conditions. The results show that the bubble detachment frequency and detachment diameter greatly increase with the operating current density. Additionally, they are significantly impacted by the local pore structure and morphology of Ti felt LGDLs. In the flow channels, there exist only several discrete micro bubbles at a low current density (0.04 A/cm 2 ) and a large flow velocity (133 mm/s). At a current density (0.2 A/cm 2 ) and a flow velocity (67 mm/s), a number of gas slugs are formed in the follow channels, in addition to discrete micro bubbles. At a high current density (1 A/cm 2 ) and a flow velocity (67 mm/s), more bubbles appear in the channel, and the flow field is dominated by slug or annular flows. These investigations can help to better understand the two-phase flow and bubble detachment mechanism, and provide a foundation for electrochemical reaction, multiphase flow studies and optimize the design of gas diffusion layers and flow fields for PEMECs in the future.

08 HYDROGEN↗
DOE OSTI2019

Development of High-Pressure Multigrain X-Ray Diffraction for Exploring the Earth’s Interior

The lower mantle comprises >55% by volume of our planet and enigmatic seismic features in the deep lower mantle call for interperations from high pressure-temperature (P-T) petrological experiments. Studying phase changes in a realistic bulk composition or with addition of volatiles requires in situ phase identification and structure analysis in a multiphase assemblage under corresponding high P-T conditions. In a powder X-ray diffraction (XRD) pattern, however, peaks of new structures and minor phases are often overshadowed by the diffraction of coexisting major phases. In the past 5 years, we have overcome the difficulties and successfully studied new structures and minor phases in great detail in situ under deep lower manlte conditions applying the multigrain XRD technique that treats each individual 100-500 nm-size grains as a single-crystal in a sample consisting of multiple phases and multiple crystals. This paper summarized the advantages of the multigrain XRD technique and provided an example of identification of hydrous phases with average ~100 nm grain size in a multiphase assemblage through combination of in situ phase identification and ex situ microscopic analysis.

58 GEOSCIENCES↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗
DOE OSTI2018

Pressure transient analysis during CO 2 push-pull tests into faults for EGS characterization

Here, with the goal of detecting and characterizing faults and fractures in enhanced geothermal systems (EGS), a new technology involving CO 2 push-pull testing, active-source geophysical imaging, and well logging has recently been proposed. This technique takes advantage of (1) the contrasting properties of supercritical CO 2 and water which cause CO 2 to appear distinct from surrounding brine in seismic and other geophysical logging approaches, (2) the non-wetting nature of CO 2 which keeps it localized to the faults and fractures to create contrast potentially sufficient for active seismic and well-logging approaches to image faults and fracture zones at EGS sites. In this study, we evaluate the feasibility of using pressure transient monitoring during CO 2 push-pull tests to complement active seismic and wireline well logging for EGS characterization. For this purpose, we developed a 2D model of a prototypical geothermal site (Desert Peak, NV) that includes a single fault. The fault zone consists of a slip plane, fault gouge, and damage zone, and is bounded by the surrounding matrix of the country rock. Through numerical simulation using iTOUGH2, we found that the pressure transient at the monitoring wells in the fault gouge shows unique traits due to the multiphase flow conditions developed by CO 2 injection, and varies sensitively on the arrival of the CO 2 plume and the degree of CO 2 saturation. A sensitivity analysis shows the pressure transient is most sensitive to the fault gouge permeability, but also depends on multiphase flow parameters and the boundary conditions of the fault. Lastly, an inversion study reveals that the fault gouge permeability can be best estimated with the pressure transient data, whereas additional CO 2 saturation data do not improve the accuracy of the inversion significantly.

58 GEOSCIENCES↗
DOE OSTI2016

Modeling the Impact of Carbon Dioxide Leakage into an Unconfined, Oxidizing Carbonate Aquifer

Multiphase, reactive transport modeling was used to identify the mechanisms controlling trace metal release under elevated CO2 conditions from a well-characterized carbonate aquifer. Modeling was conducted for two experimental scenarios: batch experiments to simulate sudden, fast, and short-lived release of CO2 as would occur in the case of well failure during injection, and column experiments to simulate more gradual leaks such as those occurring along undetected faults, fractures, or well linings. Observed and predicted trace metal concentrations are compared to groundwater concentrations from this aquifer to determine the potential for leaking CO2 to adversely impact drinking water quality. Finally, a three-dimensional multiphase flow and reactive-transport simulation of CO2 leakage from an abandoned wellbore into a generalized model of the shallow, unconfined portion of the aquifer is used to determine potential impacts on groundwater quality. As a measure of adverse impacts on groundwater quality, both the EPA’s MCL limits and the maximum trace metal concentration observed in the aquifer were used as threshold values.

Aquifer↗
DOE OSTI2020

Time-Lapse Gravity Monitoring of CO2 Migration Based on Numerical Modeling of a Faulted Storage Complex

In this study, the performance of both surface and borehole time-lapse gravity monitoring to detect CO2 leakage from a carbon storage site is evaluated. Several hypothetical scenarios of CO2 migration in a leaky fault, and thief zones at different depths at the Kimberlina site (California, USA) constitute the basis of the approach. The CO2 displacement is simulated using the TOUGH2 simulator applied to a detailed geological model of the site. The gravity responses to these CO2 plumes are simulated using forward modeling with sensors at ground surface and in vertical boreholes. Results of inversion on one scenario are also presented. The surface-based gravity responses obtained for the different leakage scenarios demonstrate that leakage can be detected at the surface in all the scenarios but the time to detection is highly variable (10 to 40 years) and dependent on the detection threshold considered. Borehole measurements of the vertical component of gravity provide excellent constrains in depth when they are located in proximity of the density anomaly associated with the presence of CO2, thus discriminating multiple leaks in different thief zones. Joint inversion of surface and borehole data can bring valuable information of the occurrence of leakages and their importance by providing a reasonable estimate of mass of displaced fluids. This study demonstrates the importance of combining multiphase flow simulations with gravity modeling in order to define if and when gravity monitoring would be applicable at a given storage site.

Time-lapse gravity monitoring, leak detection, CCS↗
DOE OSTI2018

Hard X-ray fluorescence spectroscopy of high pressure cavitating fluids in aluminum nozzles

X-rays are frequently used to study the internal geometry of dense objects, and to measure the density of multiphase flows. However, quantitatively measuring fluid density inside a metallic object such as a high pressure spray nozzle is difficult. X-rays of sufficiently high energy to penetrate a metal object are not appreciably absorbed by the fluid inside. This requires the use of plastic or beryllium test sections, which are not suited to high pressure conditions. We present a high-energy X-ray fluorescence technique which can overcome this problem. The experiments were conducted at the 7-BM beamline of the Advanced Photon Source at Argonne National Laboratory. A hydrocarbon fluid was seeded with cerium nanoparticles. The fluid was pumped at high pressure through aluminum nozzles with inner diameters of 0.36-0.90 mm and wall thicknesses of 2-3 mm. A collimated, monochromatic 42.5 keV X-ray beam excited K-edge fluorescence from the cerium. The K-alpha emission lines at 34-35 keV were recorded by a cryogenic germanium detector. Changes in fluid density due to cavitation of the liquid inside the nozzle were measured by raster scanning the nozzle across the beam. A spatial resolution of 20 x 20 mu m(2) was achieved with a slitted beam, which was improved to 5 x 10 mu m(2) with X-ray focusing mirrors. The uncertainty in the path-integrated vapor fraction was 30-40 mu m at 95% confidence. A limitation of this approach is that for low vapor pressure fluids, the nanoparticles increase the vapor pressure of the fluid and act as additional nucleation sites. These experiments demonstrate a path forward for measurements of multiphase flows inside metal components under conditions that are not feasible in optically accessible materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

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

Changes in Liver Volume Observed Following Sorafenib and Liver Radiation Therapy

Purpose: The purpose of this study was to quantify unexpected liver volume reductions in patients treated with sorafenib prior to and during liver radiation therapy (RT). Methods and Materials: Fifteen patients were treated in a phase 1 study of sorafenib for 1 week, followed by concurrent sorafenib-RT (in 6 fractions). Patients had either focal cancer (treated with stereotactic body RT [SBRT]) or diffuse disease (treated with whole-liver RT). Liver volumes were contoured and recorded at planning (day 0) from the exhale CT. After 1 week of sorafenib (day 8), RT image guidance at each fraction was performed using cone beam CT (CBCT). Planning liver contours were propagated and modified on the reconstructed exhale CBCT. This was repeated in 12 patients treated with SBRT alone without sorafenib. Three subsequent patients (2 sorafenib-RT and 1 non-sorafenib) were also assessed with multiphasic helical breath-hold CTs. Results: Liver volume reductions on CBCT were observed in the 15 sorafenib-RT patients (median decrease of 68 cc, P=.02) between day 0 and 8; greater in the focal (P=.025) versus diffuse (P=.52) cancer stratum. Seven patients (47%) had reductions larger than the 95% intraobserver contouring error. Liver reductions were also observed from multiphasic CTs in the 2 additional sorafenib-RT patients between days 0 and 8 (decreases of 232.5 cc and 331.7 cc, respectively) and not in the non-sorafenib patient (increase of 92 cc). There were no significant changes in liver volume between planning and first RT in 12 patients with focal cancer treated with SBRT alone (median increase, 4.8 cc, P=.86). Conclusions: Liver volume reductions were observed after 7 days of sorafenib, prior to RT, most marked in patients with focal liver tumors, suggesting an effect of sorafenib on normal liver. Careful assessment of potential liver volume changes immediately prior to SBRT may be necessary in patients in sorafenib or other targeted therapies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Coupling DAEM and CFD for simulating biomass fast pyrolysis in fluidized beds

We report results from computational simulations of an experimental, lab-scale bubbling bed biomass pyrolysis reactor that include a distributed activation energy model (DAEM) for the kinetics. In this study, we utilized multiphase computational fluid dynamics (CFD) to account for the turbulent hydrodynamics, and this was combined with the DAEM kinetics in a multi-component, multi-step reaction network. Our results indicate that it is possible to numerically integrate the coupled CFD–DAEM system without significantly increasing computational overhead. It is also clear, however, that reactor operating conditions, reaction kinetics, and multiphase flow dynamics all have major impacts on the pyrolysis products exiting the reactor. We find that, with the same pre-exponential factors and mean activation energies, inclusion of distributed activation energies in the kinetics can shift the predicted average value of the exit vapor-phase tar flux and its statistical distribution, compared to single-valued activation-energy kinetics. Perhaps the most interesting observed trend is that increasing the diversity of the DAEM activation energies appears to increase the mean tar yield, all else being equal. As a result, these findings imply that accurate resolution of the reaction activation energy distributions will be important for optimizing biomass pyrolysis processes.

09 BIOMASS FUELS↗