Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fracturing”

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 289 records · Page 16

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗

Shifting groundwater fluxes in bedrock fractures: Evidence from stream water radon and water isotopes

Geologic features (e.g., fractures and alluvial fans) can play an important role in the locations and volumes of groundwater discharge and degree of groundwater-surface water (GW-SW) interactions. However, the role of these features in controlling GW-SW dynamics and streamflow generation processes are not well constrained. GW-SW interactions and streamflow generation processes are further complicated by variability in precipitation inputs from summer and fall monsoon rains, as well as declines in snowpack and changing melt dynamics driven by warming temperatures. Using high spatial and temporal resolution radon and water stable isotope sampling and a 1D groundwater flux model, we evaluated how groundwater contributions and GW-SW interactions varied along a stream reach impacted by fractures (fractured-zone) and below the fractured hillslope (non-fractured zone) in Coal Creek, a Colorado River headwater stream affected by summer monsoons. During early summer, groundwater contributions from the fractured zone dominated, but declined throughout the summer. Groundwater contributions from the non-fractured zone were constant throughout the summer and became proportionally more important later in the summer. We hypothesize that groundwater in the non-fractured zone is dominantly sourced from a high-storage alluvial fan at the base of a tributary that is connected to Coal Creek throughout the summer and provides consistent groundwater influx. Water isotope data revealed that Coal Creek responds quickly to incoming precipitation early in the summer, and summer precipitation becomes more important for streamflow generation later in the summer. We quantified the change in catchment dynamic storage and found it negatively related to stream water isotope values, and positively related to modeled groundwater discharge and the ratio of fractured zone to non-fractured zone groundwater. Importantly, we interpret these relationships as declining hydrologic connectivity throughout the summer leading to late summer streamflow supported predominantly by shallow flow paths, with variable response to drying from geologic features based on their storage. As groundwater becomes more important for sustaining summer flows, quantifying local geologic controls on groundwater inputs and their response to variable moisture conditions may become critical for accurate predictions of streamflow.

54 ENVIRONMENTAL SCIENCES↗

Biofilm Feedbacks Alter Hydrological Characteristics of Fractured Rock Impacting Sulfidogenesis and Treatment

Flow-through fractures dominate the movement of fluids in a variety of natural as well as engineered subsurface systems. Microbial activities in fractured rock impact subsurface energy recovery, storage, and waste disposal. It has been recognized that understanding how the contrasting permeability between fracture and matrix interacts with microbial metabolism under thermal and hydrological gradients is key to effective utilization of the subsurface, yet such studies are sparse. Microorganisms mediate the production of hydrogen sulfide (also known as souring) in oil-bearing geological formations. We conducted a comprehensive experimental study of a novel 2D fractured rock system to understand these complex interactions and demonstrated how biofilm development can impact fracture flow, which subsequently feedbacks to moderate sulfidogenesis. Elevated temperature relevant to reservoir conditions interacted with the injection of cold fluid and formed a thermal gradient away from fractures, creating thermal niches for microbial activities in the fractured rock. Finally, results showed that while fracture flows were dominant in the beginning, with time, growth of the biofilm in the fractures reduced permeability, effectively moderating the initial fracture-matrix contrast, and limited microbial accessibility to nutrients and subsequent reactions rates.

59 BASIC BIOLOGICAL SCIENCES↗

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗

Hydro-Mechanical Measurements of Sheared Crystalline Rock Fractures With Applications for EGS Collab Experiments 1 and 2

We present hydro-mechanical measurements that characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens from the enhanced geothermal system (EGS) Collab Project's Experiment 1 and 2 sites (E1 and E2) at the Sanford Underground Research Facility. We employed a triaxial direct shear method augmented with X-ray imaging to perform hydroshearing (injection-induced shearing) and mechanical shearing on naturally fractured specimens at in situ stress conditions. Measurements included fracture permeability, strength, stress-dependent aperture, shear dilation, and frictional strength. Results reveal that in situ natural fractures must be permeable, weak, and shear-oriented to be hydrosheared, and only a subset of the observable in situ fractures were suitable. When sheared, the fracture permeability typically increased by a factor of 10 or more and this increase was retained over time. However, shear slip did not always result in permeability increase. High phyllosilicate content associated with exceptionally weak fractures exhibited poor or even decreased permeability after stimulation. These measurements in combination with site data were used to conduct a slip-tendency analysis for different fracture sets, and we selected the top candidate natural fractures for hydroshearing at the EGS Collab sites. We also found that the lower in situ shear stress and stronger fractures at the E2 site make hydroshearing more challenging than at the E1 site. Overall, the methods and analysis used in our work can be applied to any geothermal project to identify in situ joint sets that are best suited for hydroshearing, which in turn can help to optimize well placement and energy production.

58 GEOSCIENCES↗

Fracture Characterization Via AI‐Assisted Analysis of Temperature Logs

Abstract Fractures control fluid flow, mass transport, and heat transfer in a geothermal reservoir. This makes accurate characterization of fracture networks a prerequisite for optimal design and control of a reservoir's exploitation. We develop a deep‐learning procedure to identify fracture locations via interpretation of temporally and spatially continuous downhole temperature measurements. A long short‐term memory fully convolutional network (LSTM‐FCN) is used both to capture long‐term dependencies in sequential temperature data and to distill local features around fractures. A wellbore and fractured‐reservoir thermal model is established to generate temperature data for network training. The trained LSTM‐FCN exhibits a unique ability to detect multiple fractures intersecting a borehole. We use the LSTM‐FCN algorithm to evaluate the effectiveness of different‐stage wellbore temperature measurements on fracture detection in a complex fractured system. Our experiments reveal that the use of various‐stage temperature information as an input feature set improves the robustness of fracture detection to noise interference. This study indicates the practical feasibility of obtaining accurate fracture‐network reconstructions from temperature signals, at reasonable computational cost.

Yang, Xiaoyu↗

Gravity-driven controls on fluid and carbonate precipitation distributions in fractures

Abstract Many challenges related to carbon-dioxide ( $$\hbox {CO}_2$$ CO 2 ) sequestration in subsurface rock are linked to the injection of fluids through induced or existing fracture networks and how these fluids are altered through geochemical interactions. Here, we demonstrate that fluid mixing and carbonate mineral distributions in fractures are controlled by gravity-driven chemical dynamics. Using optical imaging and numerical simulations, we show that a density contrast between two miscible fluids causes the formation of a low-density fluid runlet that increases in areal extent as the fracture inclination decreases from 90 $$^\circ$$ ∘ (vertical fracture plane) to 30 $$^\circ$$ ∘ . The runlet is sustained over time and the stability of the runlet is controlled by the gravity-driven formation of 3D vortices that arise in a laminar flow regime. When homogeneous precipitation was induced, calcium carbonate covered the entire surface for horizontal fractures (0 $$^\circ$$ ∘ ). However, for fracture inclinations greater than 10 $$^\circ$$ ∘ , the runlet formation limited the areal extent of the precipitation to less than 15% of the fracture surface. These insights suggest that the ability to sequester $$\hbox {CO}_2$$ CO 2 through mineralization along fractures will depend on the fracture orientation relative to gravity, with horizontal fractures more likely to seal uniformly.

54 ENVIRONMENTAL SCIENCES↗

A study of explosive-induced fracture in polymethyl methacrylate (PMMA)

The fracture response of geologic materials is of interest for applications, including geothermal energy harnessing and containment of underground explosions. To better understand the explosively induced fracture response of geomaterials, polymethyl methacrylate (PMMA) was used as a transparent rock surrogate to allow imaging of internal shock propagation and fracture growth processes. Experiments were conducted using high-speed shadowgraphy and photon Doppler velocimetry (PDV), which were compared to numerical simulations. Experiments measured fractures produced in 304.8mm × 304.8mm × 304.8mm PMMA cubes with two simultaneously initiated detonators. The cubes were subjected to varying amounts and directions of externally applied uniaxial stresses, including no stress, 2 MPa stress, and 20 MPa stress. The fracture radius as a function of time was extracted from the high-speed videos. Post-test images of the PMMA cubes aided in the determination of three-dimensional effects not directly imaged by the cameras. The surface velocity history and the shock response captured in PDV and the high-speed videos were compared to the simulated explosive-induced shock response. The simulation results indicate that the shock drives the fracture for the first 20 μs corresponding to a fracture radius of approximately 15 mm in the experiments. The gas-driven fracture extent was estimated analytically using an equilibrium stress distribution calculated after the shock wave propagation through the sample. Reduction in the gas pressure due to the leakage of the explosive products through the crack as a function of time was accounted for. In conclusion, the estimated fracture lengths were in agreement with the experimentally observed fracture lengths.

15 GEOTHERMAL ENERGY↗

Interaction of a normally-incident plane wave with a nonlinear poroelastic fracture

While it has been recognized that a large amplitude incident wave upon a dry fracture can exhibit nonlinear seismic wave scattering due to its stress-dependent mechanical compliance, the impact of pore fluid in the fracture and a fluid-filled poroelastic background medium—features common for fractures in the Earth—are not well understood. As a first step toward an understanding of the nonlinear poroelastic response of elastic waves in fractured media, analytical approximate formulas are used for the amplitude and phase of a normally incident plane wave using a perturbation method, assuming a fluid-filled, highly compliant nonlinear interface embedded in a linear poroelastic solid. The stress-closure behavior of the fracture is modeled by nonlinear, poroelastic displacement-discontinuity boundary conditions (a linear–slip interface). The theory predicts that the static (“Direct current,” or DC) and higher-order-harmonic waves produced by the nonlinear scattering can be greatly reduced by the presence of fluid in the fracture. This, however, depends upon a number of parameters, including fracture compliance, fluid properties (compressibility and viscosity), and the permeability of the background medium, as well as environmental parameters such as the initial fluid pressure and stress acting on the fracture. The static effect produces low-frequency fluid pressure pulses when a finite-duration wave is incident upon the fracture—behavior unique to fluid-filled fractures within a poroelastic medium.

58 GEOSCIENCES↗

Alcohol exposure decreases osteopontin expression during fracture healing and osteopontin-mediated mesenchymal stem cell migration in vitro

Background Alcohol consumption is a risk factor for impaired fracture healing, though the mechanism(s) by which this occurs are not well understood. Our laboratory has previously shown that episodic alcohol exposure of rodents negatively affects fracture callus development, callus biomechanics, and cellular signaling which regulates stem cell differentiation. Here, we examine whether alcohol alters chemokine expression and/or signaling activity in the mouse fracture callus during early fracture healing. Methods A mouse model for alcohol-impaired tibia fracture healing was utilized. Early fracture callus was examined for alcohol-effects on tissue composition, expression of chemokines involved in MSC migration to the fracture site, and biomechanics. The effects of alcohol on MSC migration and cell adhesion receptors were examined in an in vitro system. Results Mice exposed to alcohol showed decreased evidence of external callus formation, decreased callus-related osteopontin (OPN) expression levels, and decreased biomechanical stiffness. Alcohol exposure decreased rOPN-mediated MSC migration and integrin β1 receptor expression in vitro. Conclusions The effects of alcohol exposure demonstrated here on fracture callus-associated OPN expression, rOPN-mediated MSC migration in vitro, and MSC integrin β1 receptor expression in vitro have not been previously reported. Understanding the effects of alcohol exposure on the early stages of fracture repair may allow timely initiation of treatment to mitigate the long-term complications of delayed healing and/or fracture non-union.

Orthopedics↗

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗

Femoral Fracture in Primary Soft-Tissue Sarcoma of the Thigh and Groin Treated with Intensity-Modulated Radiation Therapy: Observed versus Expected Risk

Purpose: This study was designed to compare the observed risk of femoral fracture in primary soft-tissue sarcoma (STS) of the thigh/groin treated with intensity-modulated radiation therapy (IMRT) to expected risk calculated using the Princess Margaret Hospital (PMH) nomogram. Methods: Expected femoral fracture risk was calculated by using the PMH nomogram. Cumulative risk of fracture was estimated by using Kaplan–Meier statistics. Prognostic factors were assessed with univariate and multivariate analysis using Cox’s stepwise regression. Results: Between February 2002 and December 2010, 92 consecutive eligible patients were assessed. Median follow-up was 73 months (106 months in surviving patients). IMRT was delivered preoperatively (50 Gy) in 13 (14%) patients and postoperatively in 79 (86%) patients (median dose, 63 Gy; range, 59.4–66.6 Gy). The observed crude risk of fractures was 6.5% compared with 25.6% expected risk from the nomogram; the cumulative risk of fracture using IMRT at 5 years was 6.7% (95% CI 2.8–16.0%). The median time to fracture was 23 months (range, 6.9–88.6). Significant predictors of fracture on univariate analysis were age ≥ 60 years (p = 0.03), tumor location in the anterior thigh (p = 0.008), and periosteal stripping to > 20 cm (p < 0.0001). On multivariate analysis, age ≥ 60 years and periosteal stripping > 20 cm retained significance (p = 0.04 and p = 0.009, respectively). Conclusions: In this study, the cumulative risk of femur fracture in patients treated with IMRT (6.7%) is less than the expected risk using the PMH nomogram (25.6%). Established predictors of femur fracture, such as gender, tumor size, and dose of RT, seem to have less impact on fracture risk when using IMRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Fracture modes of high modulus graphite/epoxy angleplied laminates subjected to off-axis tensile loads

Angleplied laminates of high modulus graphite fiber/epoxy were examined in several ply configurations at various tensile loading angles to the zero ply direction to determine the effects of ply orientations on tensile properties, fracture modes, and fracture surface characteristics of the various plies. Experimental results consist of stress-strain data, selected plots, fracture stresses and strains, and scanning electron microscope (SEM) photographs of fracture surfaces. It was found that the stress-strain curves were linear to fracture, and that although fracture surface characteristics for a given fracture mode are similar to those for the same fracture mode in uniaxial specimens, no simple load angle range can be associated with a given fracture mode. It was also concluded that SEM results must be supplemented with ply stress calculations in order to identify ranges of fracture modes occurring as a function of ply orientation with respect to the load direction.

Sinclair, J. H.↗

Mechanics of fracture - Fundamentals and some recent developments

An overview is presented of the fundamental aspects of and recent developments in fracture mechanics. Reference is made to linear elastic fracture mechanics including the state of stresses and displacements in the vicinity of cracks, effects of crack geometry and orientation on stress intensity factors, energy balance of Griffith, Irwin's stress intensity concept, and linear elastic fracture mechanics testing for fracture toughness. Other aspects of this paper include the non-linear behavior of materials and their influence on fracture mechanics parameters, consideration of viscoelasticity and plasticity, non-linear fracture toughness parameters as C.O.D., R-curve and J-integral, and a non-linear energy method, proposed by Liebowitz. Finite element methods applied to fracture mechanics problems are indicated. Also, consideration has been given to slow crack growth, dynamic effects on K(IC), Sih's criterion for fracture, Lee and Liebowitz's criterion relating crack growth with plastic energy, and applications of fracture mechanics to aircraft design. Suggestions are offered for future research efforts to be undertaken in fracture mechanics.

Liebowitz, H.↗

Deformation and fracture of thin sheet aluminum-lithium alloys: The effect of cryogenic temperatures

The objective is to characterize the fracture behavior and to define the fracture mechanisms for new Al-Li-Cu alloys, with emphasis on the role of indium additions and cryogenic temperatures. Three alloys were investigated in rolled product form: 2090 baseline and 2090 + indium produced by Reynolds Metals, and commercial AA 2090-T81 produced by Alcoa. The experimental 2090 + In alloy exhibited increases in hardness and ultimate strength, but no change in tensile yield strength, compared to the baseline 2090 composition in the unstretched T6 condition. The reason for this behavior is not understood. Based on hardness and preliminary Kahn Tear fracture experiments, a nominally peak-aged condition was employed for detailed fracture studies. Crack initiation and growth fracture toughness were examined as a function of stress state and microstructure using J(delta a) methods applied to precracked compact tension specimens in the LT orientation. To date, J(delta a) experiments have been limited to 23 C. Alcoa 2090-T81 exhibited the highest toughness regardless of stress state. Fracture was accompanied by extensive delamination associated with high angle grain boundaries normal to the fatigue precrack surface and progressed microscopically by a transgranular shear mechanism. In contrast the two peak-aged Reynolds alloys had lower toughness and fracture was intersubgranular without substantial delamination. The influences of cryogenic temperature, microstructure, boundary precipitate structure, and deformation mode in governing the competing fracture mechanisms will be determined in future experiments. Results contribute to the development of predictive micromechanical models for fracture modes in Al-Li alloys, and to fracture resistant materials.

Wagner, John A.↗

Fracture mechanics; Proceedings of the 22nd National Symposium, Atlanta, GA, June 26-28, 1990. Vols. 1 & 2

Current research on fracture mechanics is reviewed, focusing on ductile fracture; high-temperature and time-dependent fracture; 3D problems; interface fracture; microstructural aspects of fatigue and fracture; and fracture predictions and applications. Particular attention is given to the determination and comparison of crack resistance curves from wide plates and fracture mechanics specimens; a relationship between R-curves in contained and uncontained yield; the creep crack growth behavior of titanium alloy Ti-6242; a crack growth response in three heat resistant materials at elevated temperature; a crack-surface-contact model for determining effective-stress-intensity factors; interfacial dislocations in anisotropic bimaterials; an effect of intergranular crack branching on fracture toughness evaluation; the fracture toughness behavior of exservice chromium-molybdenum steels; the application of fracture mechanics to assess the significance of proof loading; and a load ratio method for estimating crack extension.

Ernst, Hugo A.↗

Tensile Creep Fracture of Polycrystalline Near-Stoichiometric NiAl

Tensile creep fracture behavior of polycrystalline near-stoichiometric NiAl has been studied between 700 and 1200 K under initial applied stresses varying between 10 and 200 MPa. The stress exponent for fracture varied between 5.0 and 10.7 while the activation energy for fracture was 250 +/- 22 kJ/mol. The fracture life was inversely proportional to the secondary creep rate in accordance with the Monkman-Grant relation although there was extensive scatter in the data. This observation suggests that the fracture life for near-stoichiometric NiAl was influenced by creep under these stress and temperature conditions. Several different fracture morphologies were observed. Transgranular ductile cleavage fracture occurs at 700 K and at the higher stresses at 800 K. The fracture mode transitions to transgranular creep fracture at 900 and 1000 K and at lower stresses at 800 K, while plastic rupture and grain boundary cavitation occur at 1100 and 1200 K. An experimental fracture mechanism map is constructed for near-stoichiometric NiAl.

Raj, Sai V.↗