SEARCH · Engineering Papers

Results for “Fracture Models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Immobile Pore-Water Storage Enhancement and Retardation of Gas Transport in Fractured Rock

The effect of immobile pore-water on gas transport in fractured rock has implications for numerical modeling of soil vapor extraction, methane leakage, gaseous CO 2 leakage from sequestration operations, radionuclide transport from underground nuclear explosions, and nuclear waste disposal. While the ability for immobile pore-water storage to effect gas transport has been recognized in the past, the details and specific scenarios leading to enhanced, retarded, or unaffected gas transport have not been explored. In this paper, we performed numerical investigations into the enhancement and retardation of gas transport due to immobile pore-water storage in order to identify implications for gas transport applications. To do this, we developed a numerical approach to model gas transport with a single-phase flow solution coupled to the advection–dispersion equation modified to account for immobile pore-water storage. Other than the immobility of pore water, the formulation contains all other physics included in two-phase formulations (advective and diffusive gas transport in fractures and rock matrix and dissolution in immobile pore water). The assumption of immobile pore water is valid here since for many applications involving transport of soluble gases in fractured rock, the rate of aqueous transport is insignificant compared to gas transport. We verify our modeling approach with analytical solutions of: (1) 1D gas diffusion, (2) 1D gas advection, (3) barometric pumping of a fracture, and (4) gas transport with uniform fracture flow and diffusion into the matrix. We account for pore-water storage in our model by implementing a kinetic formulation of gas dissolution wherein the dissolved (aqueous) phase is considered an immobile constituent. Using this formulation, we model the effect of dissolution rate and saturation on the retardation of gas transport during pure diffusion and pure advection. We also demonstrate that although it is commonly believed that pore-water storage will always enhance gas transport in fractures during oscillatory flow (e.g., during reversing pressure gradients such as barometric pumping cycles), our simulations indicate that this may not always be the case. Our numerical investigations indicate that scenarios with lower effective diffusion coefficients ( ⪅ 10 − 5  m 2 / s ) and lower dissolution coefficients (i.e., the dissolution diffusion coefficient) ( ⪅ 10 − 11  m 2 / s ) will result in enhanced gas transport. Other combinations of gaseous diffusion and dissolution coefficients result in delayed gas transport or insignificant effect on gas transport. Based on these results, and given the fact that the free-air diffusion coefficients for many gases are slightly above 10 − 5  m 2 / s (near the boundary of enhanced/non-enhanced gas transport), tortuosity would have to be significant for gas transport enhancement to occur. Similarly, if we consider the free-water diffusion coefficient of gases (around 10 − 9  m 2 / s for many gases of interest) to be a maximum bound for the effective dissolution coefficients (i.e., the rate of diffusion into the water away from the air/water interface will limit the dissolution rate), the effective dissolution coefficient would have to be at least 2 orders of magnitude less for enhanced gas transport to occur. Otherwise, pore-water storage will delay or have negligible effect on the gas transport. The results also indicate that the assumption of instantaneous equilibrium, often invoked in numerical codes for dissolution processes (i.e., dissolution coefficient is effectively infinite), will fail to accurately capture important details of soluble gas transport. Finally, the implication of our results is that gas breakthrough times in fractured rock during oscillatory flow may be highly dependent on the effective gaseous diffusion and dissolution coefficients.

58 GEOSCIENCES↗

DOE OSTI2020

Wellbore injectivity response to step-rate CO 2 injection: Coupled thermo-poro-elastic analysis in a vertically heterogeneous formation

Safe and permanent carbon geological storages require an elaborate analysis of geomechanical stability. Subsurface injection of CO 2 changes the local temperature and pore pressure and, further, alters the stress due to thermo-poro-elastic responses. A permanent CO 2 storage testing was conducted in the water leg of Cranfield reservoir in Mississippi, USA, where hosted CO 2 enhanced-oil recovery activities. During CO 2 injection, the injection rate was ramped up twice but the bottom-hole pressure did not increase with the imposed injection rates as expected. This unexpected field observation suggests the possibility of an open-mode fracture development at the injector. However, the injector response and potential fracture development have not been rigorously interpreted with detailed near-wellbore temperature and pore pressure change upon the CO 2 injection. In this study, we performed history matching of CO 2 injection using coupled thermo-poro-elastic reservoir simulation to examine the possibility of fracturing. We built a reservoir model including vertical heterogeneity in both petrophysical and geomechanical properties estimated from well-logging analysis and laboratory experiments. The simulation results show that CO 2 injection changes stresses and support the hypothesis of development of an open-mode fracture at the injector during the Cranfield test. The near-injector region exhibits a large temperature reduction up to 55 °C with ensuing effective horizontal stress reduction up to 9.1 MPa. However, a caprock integrity issue is unlikely because of the horizontal stress contrast within layers and low hydraulic communication with the injection zone. This work indicates that the injectant temperature should be considered in the design of high-rate CO 2 injector.

58 GEOSCIENCES↗

DOE OSTI2023

Simulation of Chilled-water Injection at EGS Collab Testbed 2 using the GEOS simulation framework

This work seeks to model the chilled-water injection that occurred at EGS Collab Testbed 2 as part of experiment 3, using the open-source GEOS simulation framework. Modeling of this process, previously conducted using STOMP-GT, confirmed that thermal breakthrough was expected to occur during the chill water injection and may have been masked by strong Joule-Thompsons effects induced by large pressure drops between the fractures and the production wells. Building a GEOS numerical model of this process will allow to include additional physics. In fact, GEOS can model fully coupled thermo-poromechanical processes in fractured porous media and can be employed to model the process of artificially enhancing rock permeability by hydraulically fracturing the rocks. For example, the GEOS model could help understanding the dynamic nature of the EGS collab experiment 3 flow system. Since the limited time and resources do not allow for an extensive study, the focus is on building a baseline model of this process. As such, the main outcomes of this modelling efforts are listed below.

15 GEOTHERMAL ENERGY↗

DOE OSTI2018

An exploration of pressure dynamics using differential equations defined on a fractal geometry

Pressure tests have been used to understand the behavior of naturally fractured reservoirs (NFRs) by establishing parameters such as permeability, connectivity of fractures, system compressibility, fracture network parameters, and others associated with the physical properties of the reservoir. The complex fracture network of a NFR produces anomalous fluid flow, which is shown by the pressure behavior at the well. A plausible model to explain this behavior is to assume that it arises from a fractal fracture system. However, to estimate the parameters that define this fractal fracture system, fluid flow models based on the modification of Darcy’s law are used. On the other hand, in this work, a fractal structure with known fractal dimension: Sierpinski gasket (SG) is proposed. The objective of this work is to solve numerically the diffusion equation on SG, and establish a comparison with flow models used in the oil industry to investigate the fractality assumption. This study is focused on the qualitative behavior of the solution of the diffusion equation defined on a SG supported by the Kigami’s theory of differential equations on fractals.

42 ENGINEERING↗

DOE OSTI2019

Prediction of the Shear Failure of Opened Rock Fractures and Implications for Rock Slope Stability Evaluation

Rock slope commonly fails due to the shear failure of rock fractures. Shear strength of rock fractures are reduced substantially once the fracture surfaces are mismatched or opened. We propose a new criterion to predict the shear strength of rock fractures in different opening states. The degree of interlocking representing the true asperity contact area is incorporated into the modified model of Saeb and Amadei. The effect of fracture opening on asperity dilation and degradation is separately considered. The transitional stress that is a critical parameter involved in the model is analytically determined based on energy consideration. The new model is validated with experimental results from direct shear tests on synthetic fractures with regular‐shaped asperities. Good agreement between the analytical solution and the experimental data confirms the capacity of the proposed model. Therefore, the model has great potential for assessing the stability of rock slopes where fractures are often opened due to stress relief and engineering disturbances.

Li, Yingchun (ORCID:0000000201083704)↗

DOE OSTI2019

Scalable En Echelon Shear‐Fracture Aperture‐Roughness Mechanism: Theory, Validation, and Implications

Abstract Shear fractures can facilitate fluid conductivity through rock. Aperture and roughness are controlling characteristics for a fracture's fluid conductivity. Inspired by en echelon fractures, we develop a shear “fracturelet” model that predicts anisotropic aperture with respect to the direction of shearing, rougher (nonplanar) rather than smoother (planar) fractures, and the bounds of this roughness for a coalesced fracture. This tendency for rougher fracture creation is validated by in situ X‐ray images and fluid conductivity measurements from triaxial direct shear experiments on anhydrite and shale. These experiments were conducted at confining stresses from 4 to 30 MPa and shear displacement magnitudes from 0 to 2 mm on initially intact rock specimens. Hydraulic, dilatational, and local fracture apertures were measured in the experiments. Apertures exhibited strong anisotropy with more conductive flow paths forming perpendicular to the direction of shearing. Local and dilatational aperture were found to be positively correlated with increasing shear displacement but hydraulic aperture was found to vary significantly, always having values smaller than the other aperture measures at factors ranging from 0.6 to 0.0. An implication of these results is that shear fractures have a mechanism for simultaneously exhibiting very low fluid conductivity and high fluid storage volume.

Frash, Luke P.↗

DOE OSTI2023

Earth Sciences Applications for Energy and Global Security: A Few Vignettes [Slides]

Understanding the physics of flow and transport in the subsurface environment is crucial to several National Security applications. In the Energy Security realm, some examples include characterizing radioactive waste disposal, exploring natural gas extraction and sequestering CO2. An example in the Global Security arena is to be able to distinguish between anthropogenic and naturally occurring seismic disturbances, the former being of particular concern as underground nuclear testing activity has seen an increase in recent years. A common theme across these applications is the ability to accurately model fluid flow in porous and/or fractured subsurface media. For decades, researchers have grappled with how to include the effects of uncertainties in these models, including but not limited to heterogeneity, unknown initial and boundary conditions and issues with scalability. More recently advances in computing and machine learning methodologies have allowed a broader exploration of the uncertainty space as we strive for real time decision making. I will present an overview of the field and discuss a few examples of how my own research in these topics has evolved over nearly two decades.

58 GEOSCIENCES↗

DOE OSTI2026

Oil and Gas Recovery Capability

Oil and gas recovery increasingly relies on stimulation and production strategies that involve complex interactions between rock deformation, fractures, and fluid flow underground. Conventional modeling approaches often simplify or separate these effects, which can limit their ability to accurately predict performance in fractured or geologically complex reservoirs.

02 PETROLEUM↗

DOE OSTI2021

A comprehensive review of proppant embedment in shale reservoirs: Experimentation, modeling and future prospects

This paper provides a comprehensive review on the application of proppants to maintain fracture permeability over the lifetime of a well based on published observations from experiments and modeling. The review identifies and describes important processes occurring during proppant embedment, during hydraulic fracturing, laboratory testing of fracture conductivity, proppant embedment and modeling of proppant embedment. Lastly, this paper identifies the challenges and knowledge gaps that also provide future avenues of research and opportunities for collaborative technological development which requires an interdisciplinary approach of science, engineering in academia, government, and private sector.

03 NATURAL GAS↗

DOE OSTI2020

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗

DOE OSTI2024

Optimization of a Lethal, Combat-Relevant Model of Sterile Inflammation in Mice for Drug Candidate Screening

ABSTRACT Introduction Extensive trauma, commonly seen in wounded military Service Members, often leads to a severe sterile inflammation termed systemic inflammatory response syndrome (SIRS), which can progress to multiple organ dysfunction syndrome (MODS) and death. MODS is a serious threat to wounded Service Members, historically causing 10% of all deaths in trauma admissions at a forward deployed combat hospital. The importance of this problem will be exacerbated in large-scale combat operations, in which evacuation will be delayed and care of complex injuries at lower echelons of care may be prolonged. The main goal of this study was to optimize an existing mouse model of lethal SIRS/MODS as a therapeutic screening platform for the evaluation of immunomodulatory drugs. Materials and Methods Male C57BL/6 mice were euthanized, and the bones and muscles were collected and blended into a paste termed tissue–bone matrix (TBX). The TBX at 12.5%–20% relative to body weight of each recipient mouse was implanted into subcutaneous pouches created on the dorsum of anesthetized animals. Mice were observed for clinical scores for up to 48 hours postimplantation and euthanized at the preset point of moribundity. To test effects of anesthetics on TBX-induced mortality, animals received isoflurane or ketamine/xylazine (K/X). In a separate set of studies, mice received TBX followed by intraperitoneal injection with 20 mg/kg or 40 mg/kg Eritoran or a placebo carrier. All Eritoran studies were performed in a blinded fashion. Results We observed that K/X anesthesia significantly increased the lethality of the implanted TBX in comparison to inhaled anesthetics. Although all the mice anesthetized with isoflurane and implanted with 12.5% TBX survived for 24 hours, 60% of mice anesthetized with K/X were moribund by 24 hours postimplantation. To mimic more closely the timing of lethal SIRS/MODS following polytrauma in human patients, we extended observation to 48 hours. We performed TBX dose–response studies and found that as low as 15%, 17.5%, and 20% TBX caused moribundity/mortality in 50%, 80%, and 100% mice, respectively, over a 48-hour time period. With 17.5% TBX, we tested if moribundity/mortality could be rescued by anti-inflammatory drug Eritoran, a toll-like receptor 4 antagonist. Neither 20 mg/kg nor 40 mg/kg doses of Eritoran were found to be effective in this model. Conclusions We optimized a TBX mouse model of SIRS/MODS for the purpose of evaluating novel therapeutic interventions to prevent trauma-related pathophysiologies in wounded Service Members. Negative effects of K/X on lethality of TBX should be further evaluated, particularly in the light of widespread use of ketamine in treatment of pain. By mimicking muscle crush, bone fracture, and necrosis, the TBX model has pleiotropic effects on physiology and immunology that make it uniquely valuable as a screening tool for the evaluation of novel therapeutics against trauma-induced SIRS/MODS.

General & Internal Medicine↗

DOE OSTI2016

Closure of Fracture Due to Cover Stress Re-establishment After Coal Mining

In situ measurements of deformations, stresses, and closure of fractures, affecting water inflow following coal mining, are challenging due to the inaccessibility of fractured rock. Here, the authors studied the closure process of the fractured rock mass with the cover stress re-establishment based on a theoretical analysis and a scale model testing. A quantitative analysis is used to study the fracture distribution in the fractured zone. A function to describe a fracture aperture distribution in the fractured zone is proposed, which takes into account the curvature and thickness of the fractured rock. The theoretical analysis and a scale model testing both indicate that the cover stress re-establishment with mining distance increasing and the relationship between the fracture closure and cover stress re-establishment both satisfy a logarithmic function. The scale model test also shows the following features: (1) the fracture ratio (which is the fracture area divided by the total area of fracture and intact rock with a unit width in the vertical or horizontal direction) in the lower part of the fractured rock mass is greater than that in the upper part; (2) the initially fast decreased of fracture ratios is then followed by a slower decrease during the cover stress re-establishment process; (3) in the upper part of the rock mass, the vertical directional fractures with small apertures are being closed with cover stress re-establishment, which indicates an increase in the water resistance reducing the seepage from these parts of the fractured zone. This study improves the general understanding of the fracture closure process and cover stress re-establishment in the fractured rock mass after coal mining ceased, and provides a theoretical basis for water resource protection in case of underground coal mining.

58 GEOSCIENCES↗

DOE OSTI2019

TOUGH-UDEC: A simulator for coupled multiphase fluid flows, heat transfers and discontinuous deformations in fractured porous media

A numerical simulator entitled TOUGH-UDEC is presented for the assessment of coupled thermal-hydraulic-mechanical processes in fractured porous media. Two existing well-established codes, TOUGH2 and UDEC, are coupled to model multiphase fluid flows, heat transfers, and discontinuous deformations in fractured porous media by means of discrete fracture representation. TOUGH2 is widely used for the modeling of heat transfers and multiphase multicomponent fluid flows, and UDEC is a well-known distinct element code for rock mechanics. The two codes are solved sequentially, with coupling parameters passed to each equation at specific intervals. After solving thermal-hydraulic equations within the TOUGH2 code, pressure and temperature information is imported into the UDEC model. Upon solving the mechanical equation within the UDEC code the calculated fracture apertures are converted to the equivalent permeability and porosity values for a TOUGH2 flow analysis. The solution is calculated by iteratively following an explicit sequence for numerical efficiency. Verifications are presented to demonstrate the capabilities of the coupled TOUGH-UDEC simulator. Three application examples of (1) shear dilation due to increased pore pressure, (2) thermal stress and (3) CO 2 injection, show that the new simulator can be an effective tool for geoengineering applications involving shear activation of fractures and faults.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Gas Flow Tightly Coupled to Elastoplastic Geomechanics for Tight- and Shale-Gas Reservoirs: Material Failure and Enhanced Permeability

We investigate coupled flow and geomechanics in gas production from extremely low permeability reservoirs such as tight and shale gas reservoirs, using dynamic porosity and permeability during numerical simulation. In particular, we take the intrinsic permeability as a step function of the status of material failure, and the permeability is updated every time step. We consider gas reservoirs with the vertical and horizontal primary fractures, employing the single and dynamic double porosity (dual continuum) models. We modify the multiple porosity constitutive relations for modeling the double porous continua for flow and geomechanics. The numerical results indicate that production of gas causes redistribution of the effective stress fields, increasing the effective shear stress and resulting in plasticity. Shear failure occurs not only near the fracture tips but also away from the primary fractures, which indicates generation of secondary fractures. These secondary fractures increase the permeability significantly, and change the flow pattern, which in turn causes a change in distribution of geomechanical variables. From various numerical tests, we find that shear failure is enhanced by a large pressure drop at the production well, high Biot's coefficient, low frictional and dilation angles. Smaller spacing between the horizontal wells also contributes to faster secondary fracturing. When the dynamic double porosity model is used, we observe a faster evolution of the enhanced permeability areas than that obtained from the single porosity model, mainly due to a higher permeability of the fractures in the double porosity model. These complicated physics for stress sensitive reservoirs cannot properly be captured by the uncoupled or flow-only simulation, and thus tightly coupled flow and geomechanical models are highly recommended to accurately describe the reservoir behavior during gas production in tight and shale gas reservoirs and to smartly design production scenarios.

02 PETROLEUM↗

DOE OSTI2016

Recovery Act. Development and Validation of an Advanced Stimulation Prediction Model for Enhanced Geothermal System

The research project aims to develop and validate an advanced computer model that can be used in the planning and design of stimulation techniques to create engineered reservoirs for Enhanced Geothermal Systems. The specific objectives of the proposal are to: 1) Develop a true three-dimensional hydro-thermal fracturing simulator that is particularly suited for EGS reservoir creation. 2) Perform laboratory scale model tests of hydraulic fracturing and proppant flow/transport using a polyaxial loading device, and use the laboratory results to test and validate the 3D simulator. 3) Perform discrete element/particulate modeling of proppant transport in hydraulic fractures, and use the results to improve understand of proppant flow and transport. 4) Test and validate the 3D hydro-thermal fracturing simulator against case histories of EGS energy production. 5) Develop a plan to commercialize the 3D fracturing and proppant flow/transport simulator. The project is expected to yield several specific results and benefits. Major technical products from the proposal include: 1) A true-3D hydro-thermal fracturing computer code that is particularly suited to EGS, 2) Documented results of scale model tests on hydro-thermal fracturing and fracture propping in an analogue crystalline rock, 3) Documented procedures and results of discrete element/particulate modeling of flow and transport of proppants for EGS applications, and 4) Database of monitoring data, with focus of Acoustic Emissions (AE) from lab scale modeling and field case histories of EGS reservoir creation.

15 GEOTHERMAL ENERGY↗

DOE OSTI2016

Implications of Grain-Scale Mineralogical Heterogeneity for Radionuclide Transport in Fractured Media

The geological disposal of nuclear waste is based on the multi-barrier concept, comprising various engineered and natural barriers, to confine the radioactive waste and isolate it from the biosphere. Some of the planned repositories for high-level nuclear waste will be hosted in fractured crystalline rock formations. The potential of these formations to act as natural transport barriers is related to two coupled processes: diffusion into the rock matrix and sorption onto the mineral surfaces available in the rock matrix. Different in situ and laboratory experiments have pointed out the ubiquitous heterogeneous nature of the rock matrix: mineral surfaces and pore space are distributed in complex microstructures and their distribution is far from being homogeneous (as typically assumed by Darcy-scale coarse reactive transport models). We use a synthetically generated fracture–matrix system to assess the implications of grain-scale physical and mineralogical heterogeneity on cesium transport and retention. The resulting grain-scale reactive transport model is solved using high-performance computing technologies, and the results are compared with those derived from two alternative models, denoted as upscaled models, where mineral abundance is averaged over the matrix volume. In the grain-scale model, the penetration of cesium into the matrix is faster and the penetration front is uneven and finger-shaped. The analysis of the cesium breakthrough curves computed at two different points in the fracture shows that the upscaled models provide later first-arrival time estimates compared to the grain-scale model. The breakthrough curves computed with the three models converge at late times. These results suggest that spatially averaged upscaled parameters of sorption site distribution can be used to predict the late-time behavior of breakthrough curves but could be inadequate to simulate the early behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

Matrix Diffusion in Fractured Media: New Insights Into Power Law Scaling of Breakthrough Curves

Here, we develop a theoretical model for power law tailing behavior of transport in fractured rock based on the relative dominance of the decay rate of the advective travel time distribution, modeled using a Pareto distribution (with tail decaying as ~ time -(1+α) ), versus matrix diffusion, modeled using a Lévy distribution. The theory predicts that when the advective travel time distribution decays sufficiently slowly (α<1), the late-time decay rate of the breakthrough curve is -(1+α/2) rather than the classical -3/2. Yet, if α>1, the -3/2 decay rate is recovered. For weak matrix diffusion or short advective first breakthrough times, we identify an early-time regime where the breakthrough curve follows the Pareto distribution, before transitioning to the late-time decay rate. The theoretical predictions are validated against particle tracking simulations in the three-dimensional discrete fracture network simulator dfnWorks, where matrix diffusion is incorporated using a time domain random walk.

58 GEOSCIENCES↗

NASA NTRS1985

Fracture and crack growth in orthotropic laminates

A mathematical model based on the classical shear-lag assumptions is used to study the residual strength and fracture behavior of composite laminates with symmetrically placed buffer strips. The laminate is loaded by a uniform remote longitudinal tensile strain and has initial damage in the form of a transverse crack in the parent laminate between buffer strips. The crack growth behavior as a function of material properties, number of buffer-strip plies, spacing, width of buffer strips, longitudinal matrix splitting, and debonding at the interface is studied. Buffer-strip laminates are shown to arrest fracture and increase the residual strengths significantly over those of one material laminates, with S-glass being a more effective buffer strip material than Kevlar in increasing the damage tolerance of graphite/epoxy panels. For a typical graphite/epoxy laminate with S-glass buffer-strips, the residual strength is about 2.4 times the residual strength of an all graphite/epoxy panel with the same crack length. Approximately 50% of this increase is due to the S-glass/epoxy buffer-strips, 40% due to longitudinal splitting of the buffer strip interface and 10% due to bonding.

Goree, James G.↗