Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗
DOE OSTI2019

A general peridynamics model for multiphase transport of non-Newtonian compressible fluids in porous media

Herein a general state-based peridynamics model is developed to simulate transport of fluids in an arbitrary heterogeneous porous medium. The generality encompasses modeling of multiphase, multi-component flow of non-Newtonian and compressible fluids, which is often encountered in but not limited to subsurface reservoirs. Peridynamic model is especially useful for solving non-local problems, such as crack propagation, since it does not assume spatial continuity of field variables. Thus, the formulation presented here, combined with peridynamics-based damage model, can be used to simulate hydraulic fracturing with complex fluids. To demonstrate its capability to simulate multi-phase flow in porous media, the derived model is verified against the analytical Buckley-Leverett solution for immiscible Newtonian two-phase flow. Further, the non-Newtonian two-phase fluid flow in porous media is verified by simulating the polymer flood process involving immiscible displacement of a Newtonian fluid by a non-Newtonian fluid against a generalized solution obtained by Wu et al. The non-local solutions are shown to be consistent with the corresponding local solutions in limiting cases. Moreover, mass conservation of all the phases is satisfied, irrespective of discretization and extent of non-locality.

42 ENGINEERING↗
DOE OSTI2017

Advanced Analysis of Hydraulic Fracture Propagation and Efficiency Using Natural Gas as an Alternative to Water

Over the past years, new technologies have been presented as an alternative to minimize or eliminate the use of water in hydraulic fracturing for unconventional reservoirs, improve resource recovery, cost management and reduce environmental impact. As part of these technologies, energized fluids have been developed to mitigate the environmental impact caused by the application of large volumes of water during hydraulic fracturing. When a fracture is energized, the amount of water pumped during the fracturing process is reduced or eliminated by adding a gas component to the fracturing fluid. The less water that is pumped, the less damage it will create. A new technology called “Energized Natural Gas (ENG)” has recently been proposed for effective hydraulic fracturing in unconventional reservoirs. The ENG process shows that it may be a viable solution for the oil and gas industry. Due to the complexity of multiphase fluid behavior and fracture mechanics, further research and documentation is required to effectively apply gas and gas-water mixtures as a treatment fluid. The objective of this research is to further develop and demonstrate the use of natural gas as an optimum fracturing fluid for tight shale gas formations. GeoMechanics Technologies proposes to complete advanced analyses of hydraulic fracture propagation efficiency using natural gas as an alternative to water. A comprehensive review on the transport properties of gas and water mixtures was performed to investigate the rheology differences of conventionally-used fresh water and the natural gas mixture. This research helped to determine how rheology affects hydraulic fracturing. Then we developed a coupled 3D geomechanical model and fluid flow model to simulate hydraulic fracturing using various mixtures of natural gas and water to determine an effective combination. Gas production was also simulated to further test the effectiveness of the induced fracture. Finally, the technical and economic feasibility of this novel technology was considered and documented.

02 PETROLEUM↗
DOE OSTI2018

The Development of Wormholes in Laboratory-Scale Fractures: Perspectives From Three-Dimensional Simulations

Here, we investigate the development of wormholes in laboratory–scale fractures using three–dimensional numerical simulations. Well–controlled initial conditions, involving a small perturbation near the inlet of an otherwise flat fracture aperture field, were used to make a systematic study of the effects of flow rate and reaction rate on the aperture evolution. We find at least two characteristic wormhole shapes, which can be grouped within a phase diagram in the space of Péclet and Damköhler numbers. We investigate how this phase diagram depends on fracture geometry, specifically the length ( L ) and width ( W ) in comparison to the initial aperture ( h 0 ). This information is used to determine an effective Damköhler number that leads to a Péclet and length–independent phase boundary between wormholes and uniform dissolution. We relate these observations to experimental studies of wormhole formation in dissolving fractures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Monte Carlo Simulation of Alloy Design Techniques: Fracture and Welding Studied Using the BFS Method for Alloys

Large-scale simulations of dynamic processes at the atomic level have developed into one of the main areas of work in computational materials science. Until recently, severe computational restrictions, as well as the lack of accurate methods for calculating the energetics, resulted in slower growth in the area than that required by current alloy design programs. The Computational Materials Group at the NASA Lewis Research Center is devoted to the development of powerful, accurate, economical tools to aid in alloy design. These include the BFS (Bozzolo, Ferrante, and Smith) method for alloys (ref. 1) and the development of dedicated software for large-scale simulations based on Monte Carlo- Metropolis numerical techniques, as well as state-of-the-art visualization methods. Our previous effort linking theoretical and computational modeling resulted in the successful prediction of the microstructure of a five-element intermetallic alloy, in excellent agreement with experimental results (refs. 2 and 3). This effort also produced a complete description of the role of alloying additions in intermetallic binary, ternary, and higher order alloys (ref. 4).

Bozzolo, Guillermo H.↗
DOE OSTI2020

The Influence of an Interlayer on Dual Hydraulic Fractures Propagation

Multi-cluster hydraulic fracturing of long-range horizontal wells is an approach for enhancing the productivity of low-permeability shale reservoirs. In this study, RFPA-Petrol (rock failure process analysis on petroleum problems) is applied for modeling hydraulic fracture propagation in multilayered formations. RFPA-Petrol based on coupled hydraulic-mechanical-damage (HMD) modeling was first tested by modeling a laboratory scale experiment on a physical (cement) model with a single completion. The modeling demonstrated the capability of RFPA-Petrol for simulating hydraulic fracture propagation. Then, we used RFPA-Petrol to investigate how the difference in material properties between oil-bearing layers and interlayers and the fracturing fluid properties influence the propagation of dual fractures in multilayered laboratory-scale models. In this case, the models with geological discontinuities in the vertical direction are strongly heterogeneous and RFPA-Petrol simulations successfully modeled the fracture configurations.

58 GEOSCIENCES↗
DOE OSTI2023

Utah FORGE Project 3-2417: Simulations for Distributed Acoustic Sensing Strain Signatures as an Indicator of Fracture Connectivity

This dataset encompasses simulations of strain signatures from both hydraulically connected and "near-miss" fractures in enhanced geothermal systems (EGS). The files and results are presented from the perspective of digital acoustic sensing's (DAS) potential to differentiate the two fracture types. This dataset was acquired by the FOGMORE R&D project (Fiber Optic Geophysical MOnitoring of Reservoir Evolution), under Utah FORGE R&D Project 3-2417. Included are simulation and results via MatLab and COMSOL files, as well as a thesis and paper summarizing the results. Some stimulated fractures may be incomplete, approaching but not intersecting the production well. These "near-miss" fractures can be addressed in future stimulation stages or re-stimulated to complete the connection. We propose the use of fiber optic distributed acoustic sensing (DAS) as a method by which near-miss stimulated fractures may be identified and distinguished from hydraulically connected fractures. The low-frequency sub-nanostrain signatures of both complete and near-miss fractures in DAS data are simulated in this study using a hydrogeomechanical discrete fracture network model. The spatial distribution of strain was found to be an accurate indicator. However, this indicator must be evaluated in the context of DAS gauge length and spatial sampling. These simulations are a precursor to tests conducted at FORGE in 2023.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Predicting ductility of Mg/SiCp nanocomposite under multiaxial loading conditions based on unit cell modeling

This article presents an investigation on the ductility of Mg/SiCp under nine different multiaxial loading conditions based on a unit cell (or called representative volumetric element, RVE) finite element model. The constitutive relationship of magnesium matrix is characterized by a model combined Lode angle dependent plasticity model with the CPB06 plasticity model, which is able to describe the tension-compression asymmetry in magnesium matrix. The debonding of the matrix-particle interface and matrix fracture are considered by the modified Mohr-Coulomb (MMC) fracture model. Simulation results show that the interface debonding and subsequent matrix fracture are responsible for the ductility and fracture of nanocomposite. The computational results with parametric studies are compared and well correlated with experimental data, and it finally gives a stress-based full fracture locus of Mg/SiCp nanocomposite.

42 ENGINEERING↗
DOE OSTI2017

Phase field model of fluid-driven fracture in elastic media: Immersed-fracture formulation and validation with analytical solutions

Propagation of fluid-driven fractures plays an important role in natural and engineering processes, including transport of magma in the lithosphere, geologic sequestration of carbon dioxide, and oil and gas recovery from low-permeability formations, among many others. The simulation of fracture propagation poses a computational challenge as a result of the complex physics of fracture and the need to capture disparate length scales. Phase field models represent fractures as a diffuse interface and enjoy the advantage that fracture nucleation, propagation, branching, or twisting can be simulated without ad hoc computational strategies like remeshing or local enrichment of the solution space. Here we propose a new quasi-static phase field formulation for modeling fluid-driven fracturing in elastic media at small strains. The approach fully couples the fluid flow in the fracture (described via the Reynolds lubrication approximation) and the deformation of the surrounding medium. The flow is solved on a lower dimensionality mesh immersed in the elastic medium. This approach leads to accurate coupling of both physics. We assessed the performance of the model extensively by comparing results for the evolution of fracture length, aperture, and fracture fluid pressure against analytical solutions under different fracture propagation regimes. Thus, the excellent performance of the numerical model in all regimes builds confidence in the applicability of phase field approaches to simulate fluid-driven fracture.

58 GEOSCIENCES↗
DOE OSTI2021

"Hybrid fracture/matrix modeling for well completion options evaluation"

Orientation and completion for well pairs that have been subjected to multi-zonal stimulation play a critical role in the long-term performance of an Enhanced Geothermal Reservoir. Here we present the development of a methodology to rapidly and efficiently numerically simulate mixed fracture-matrix flow systems for evaluation of well design and completion options. The methodology is based on a loose coupling framework, allowing the fracture and matrix systems to be meshed separately. The fracture system includes the integration of fracture growth and aperture data from well stimulation simulations of stochastically generated fracture networks. Automatic mesh refinement is used in the matrix simulation to resolve heat transfer near the fracture network. This simulation framework is used to efficiently determine optimal production and injection well placement using adaptive sampling.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Simulation of Microcapsule Transport in Fractured Media Using Coupled CFD‐DEM

Geothermal energy is sustainable and gaining momentum as a solution to energy crises and environmental issues. However, challenges like production temperature and thermal breakthrough can impact geothermal project efficiency. One innovative solution to alleviate the thermal breakthrough is to inject polymer-based materials that are encapsulated in microcapsules into fractures to modify fracture permeability and prevent preferential flow. In our study, we utilized a coupled computational fluid dynamics and discrete element method to simulate the transport of microcapsules under various scenarios controlled by microcapsule size, microcapsule concentration, and fracture roughness. For a smooth fracture, the results indicate that small microcapsules can travel through a smooth fracture regardless of their concentrations. Large microcapsules can transport through a smooth fracture when present in lower concentrations. However, medium and mixed-size microcapsules tend to cause the sealing of a smooth fracture, irrespective of their concentrations. For a rough fracture, the transport of microcapsules is complicated by their interactions with the rough fracture walls. The presence of two sealing positions in a rough fracture adds further complexity to this transport phenomenon. The size and concentration of microcapsules control one sealing location, while the rough fracture walls determine the other sealing location. The rough walls substantially affect microcapsule transport, rendering the role of microcapsule size and concentration less significant. The simulation results suggest that complex fracture surfaces significantly elevate the occurrence of sealing behavior. To mitigate sealing behavior within more complex fractures, it would be beneficial to use smaller and lower concentrations of microcapsules.

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

A Unified Cohesive Zone Approach to Model Ductile Brittle Transition in Reactor Pressure Vessel Steels

In this study, a unified cohesive zone model has been proposed to predict, Ductile to Brittle Transition, DBT, in Reactor Pressure Vessel, RPV, steels. A general procedure is described to obtain the Cohesive Zone Model, CZM, parameters for the different temperatures and fracture probabilities. In order to establish the full master-curve, the procedure requires three calibration points with one at the upper-shelf for ductile fracture and two for the fracture probabilities, Pf, of 5% and 95% at the lower-shelf. In the current study, these calibrations were carried out by utilizing the experimental fracture toughness values and flow curves. After the calibration procedure, the simulations of fracture behavior (ranging from completely unstable to stable crack extension behavior) in one inch thick compact tension specimens at different temperatures yielded values that were comparable to the experimental fracture toughness values, indicating the viability of such unified modeling approach.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Test-Free Fracture Toughness

Computational simulation results can give the prediction of damage growth and progression and fracture toughness of composite structures. The experimental data from literature provide environmental effects on the fracture behavior of metallic or fiber composite structures. However, the traditional experimental methods to analyze the influence of the imposed conditions are expensive and time consuming. This research used the CODSTRAN code to model the temperature effects, scaling effects and the loading effects of fiberbraided composite specimens with and without fiber-optic sensors on the damage initiation and energy release rates. The load-displacement relationship and fracture toughness assessment approach is compared with the test results from literature and it is verified that the computational simulation, with the use of established material modeling and finite element modules, adequately tracks the changes of fracture toughness and subsequent fracture propagation for any fiberbraided composite structure due to the change of fiber orientations, presence of large diameter optical fibers, and any loading conditions.

Minnetyan, Levon↗
NASA NTRS2003

Test-Free Fracture Toughness

Computational simulation results can give the prediction of damage growth and progression and fracture toughness of composite structures. The experimental data from literature provide environmental effects on the fracture behavior of metallic or fiber composite structures. However, the traditional experimental methods to analyze the influence of the imposed conditions are expensive and time consuming. This research used the CODSTRAN code to model the temperature effects, scaling effects and the loading effects of fiber/braided composite specimens with and without fiber-optic sensors on the damage initiation and energy release rates. The load-displacement relationship and fracture toughness assessment approach is compared with the test results from literature and it is verified that the computational simulation, with the use of established material modeling and finite element modules, adequately tracks the changes of fracture toughness and subsequent fracture propagation for any fiber/braided composite structure due to the change of fiber orientations, presence of large diameter optical fibers, and any loading conditions.

Minnetyan, Levon↗
DOE OSTI2021

Nanoscale Interfacial Smoothing and Dissolution during Unconventional Reservoir Stimulation Implications for Hydrocarbon Mobilization and Transport

Hydraulic fracturing of low-permeability rocks significantly enhances hydrocarbon production from unconventional reservoirs. However, fluid transport through low-permeability rocks, and the influence of geochemical transformations on pore networks, are poorly constrained. Mineral reactivity during interactions with injected water may alter the physical nature of the rock, which may affect hydrocarbon mobility. To assess alteration to the rock, we previously conducted a hydrothermal experiment that reacted cubed rock samples (1 cm3) with synthetic hydraulic fracturing fluid to simulate physicochemical reactivity during hydraulic fracturing. We analyzed unreacted and reacted rocks by small-angle neutron scattering and high-pressure mercury intrusion to determine how nanoscale (~1-1400 nm) pore networks of unconventional reservoir rocks are influenced by reaction with hydraulic fracturing injectates. We show that fluid-rock interactions have a two-fold influence on hydrocarbon recovery, promoting both hydrocarbon mobilization and transport. Pore-matrix interfaces smooth via removal of clay mineral surface asperities, reducing available surface area for hydrocarbon adsorption by ca. 12-75%. Additionally, hydraulic fracturing fluid-induced dissolution creates new pores with diameters ranging in size from 800-1400 nm, increasing the permeability of the rocks by a factor of 5-10. These two consequences of mineral dissolution act in concert to release hydrocarbons from the host rock and facilitate transport through the rock during unconventional reservoir production.

Herz-Thyhsen, Ryan↗
DOE OSTI2022

Efficient co-solution of time step size and independent state in simulations of fluid-driven fracture propagation with embedded meshes

Here we present an efficient time-continuation scheme for fluid-driven fracture propagation problems in the extended finite element method framework. The approach applies a monolithic solution strategy to a fully coupled and implicit approximation of hydro-mechanical systems in conjunction with simultaneous linear elastic propagation of multiple fractures. At the end of each time step, the process ensures that the weakest fracture tip is in an equilibrium propagation regime. Furthermore, the solution process provides an initialization procedure for the newly created fracture spaces and an a priori estimate of the stress intensity factor growth rate, improving simulation robustness, and efficiency. The solution process is validated using the Kristianovich-Geertsma-de Klerk analytical solution under the toughness- and viscosity-dominated regimes. It is also extended to and demonstrated on problems with multiple fractures undergoing simultaneous propagation with stress shadow interactions. Numerical examples demonstrate that the solution process can reduce the required computational cost by one order of magnitude compared to other existing methods.

42 ENGINEERING↗
NASA NTRS2006

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized From Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, intersticies, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Glaessgen, Edward H.↗
DOE OSTI2014

Hydroshear Simulation Lab Test 2, Additional Data Sheets

This file contains additional data for this test (hydroshear test 3.2). In this test a sample of granite with a pre-cut (man made fracture) is confined, heated and differential stress is applied. The max temperature in this this system development test is 95C. More test details can be found on the spreadsheets.

15 GEOTHERMAL ENERGY↗