Fracture statistics and integration into computational models.
Abstract not provided.
Portone, Teresa↗
SEARCH · Engineering Papers
Results for “fracture modelling”
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.
Three-Dimensional Evaluation of Sand Particle Fracture Using Discrete-Element Method and Synchrotron Microcomputed Tomography Images
Recent research showed that fracture of sand particles plays a significant role in determining the plastic bulk volumetric changes of granular materials under different loading conditions. One of the major tools used to better understand the influence of particle fracture on the behavior of granular materials is discrete-element modeling (DEM). This paper employed the bonded block model (BBM) to simulate the fracture behavior of sand. Each sand particle is modeled as an agglomerate of rigid blocks bonded at their contacts using the linear-parallel contact model, which can transmit both moment and force. DEM simulated particles closely matched the actual three-dimensional (3D) shape of sand particles acquired using high-resolution 3D synchrotron microcomputed tomography (SMT). Results from unconfined one-dimensional (1D) compression of a single synthetic silica cube were used to calibrate the model parameters. Particle fracture was investigated for specimens composed of three sand particles that were loaded under confined 1D compression. Breakage energy measured from DEM models matched well with that measured experimentally. The paper studied the effects of contact loading condition and particle interaction on the fracture mode of particles using BBM that can closely capture the 3D shape of real sand particles.
58 GEOSCIENCES↗
Simulation of flame acceleration and deflagration-to-detonation transition with heat transfer in HE product and fractured HE
A simple heat transfer model has been developed for high explosives. For the gaseous high explosive (HE) product, the Smargorinsky eddy viscosity model has been adopted for estimating the turbulent thermal conductivity. Convective heat transfer in the fractured solid HE is modeled as the effective thermal conduction in the porous media, which is equivalent to the enthalpy flux carried by the velocity determined by Darcy’s law with Forchheimer term.Other features of the model include a multi-reaction chain, equation-of-state for each species involved, and a two-phase hotspot model. Test calculations are performed with various initial disturbances. Onset and progress of flame acceleration (FA) and chemical reactions are strongly affected by the heat transfer and initial disturbance. Staged development of FA also affirms the necessity of modeling based on a multi-reaction chain. When conditions warrant,e.g., stronger initial disturbance, FA evolves to deflagration-to-detonation transition (DDT),demonstrating automatic capture of the DDT.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
The ductile fracture of anisotropic materials - A theoretical and experimental evaluation of the Dugdale model applied to anisotropic materials
Anisotropic materials ductile fracture involving crack initation and propagation, using Dugdale mathematical model
Brinson, H. F.↗
Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing
Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.
58 GEOSCIENCES↗
Coupled Thermo–Hydro–Mechanical Modeling of Hydro-Shearing Stimulation in an Enhanced Geothermal System in the Raft River Geothermal Field, USA
The Raft River Enhanced Geothermal System (EGS) demonstration project aims to improve the geothermal production by enhancing the reservoir permeability via shear stimulation in well RRG-9. In this paper, we performed a series of 3D thermo– hydro–mechanical (THM) simulations to investigate the influence of hydraulic and thermal effects on the development of this EGS. The model includes synthetic fracture populations based on borehole televiewer images and in-situ stress measurements from well RRG-9. Fracture permeability evolution is determined using an empirical permeability law developed from laboratory experiments. The model was calibrated by comparing the hydraulic response to field observations, including wellhead pressure, injection rate, and well injectivity. Particularly, we analyzed the enhancement of reservoir permeability and the spatial extent of the stimulation zone for the given injection schedule. Our results indicate that the permeability enhancement of fractured geothermal reservoir is caused by the combined effects of injection-induced cooling and fluid pressure increase. The decrease of temperature plays a dominant role in reactivation of natural fractures under the hydro-shearing mechanism, while the higher injection pressure promotes shear failure and enlarges the stimulation zone. Finally, for the specified extensional stress state, the model favors greatest permeability enhancement along the maximum horizontal principal stress, a moderate vertical enhancement, and a smaller gain along the minimum horizontal principal stress.
58 GEOSCIENCES↗
Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction
Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.
Chen, Bin↗
The role of stress and fluid saturation on the acoustic response of fractured rock
Standard rock physics models are formulated to describe the behavior of porous sedimentary reservoirs, with clean sandstones being the archetypal system; however, many situations demand geophysical monitoring of rocks with significantly different structures, such as low porosity, fractured reservoirs. Conventional models also suggest that these “stiff” reservoirs can be challenging to monitor seismically due to small fluid substitution effects, but the presence of fractures leads to stress dependence which may be leveraged for remote monitoring purposes. Using samples from the Duperow Formation (dolostone) obtained from the Danielson test well in Kevin Dome, MT, we conducted ultrasonic and multi-scale structural (profilometry, synchrotron micro-tomography, pressure sensitive film) measurements on naturally fractured core in order to characterize the effects of fluid substitution and effective stress on the acoustic response of fractured reservoir rock with a focus in particular on the textural and seismic characteristics of natural fractures. We find that changes in effective stress can yield changes in velocity of up to 20% and changes in attenuation up to 200%. Measured fluid substitution effects are resolvable, but stress effects dominate. These measurements provide insight into the physical processes controlling acoustic response of fractured rocks in general and can also be used to inform monitoring efforts in fractured reservoirs.
58 GEOSCIENCES↗
A deep learning upscaling framework: Reactive transport and mineral precipitation in fracture-matrix systems
Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes for which effective medium approximations may become inaccurate, e.g., mineral precipitation and clogging during reactive transport. We have developed a deep learning upscaling framework, in which pore-scale modeling is directly employed in macroscale systems, without relying on effective medium approximations. The upscaling framework is first developed for general multiscale systems and then applied to modeling reactive transport with mineral precipitation in the altered layer in fracture-matrix structures. Solute transport from the fractures to the matrix is modeled as a wall boundary condition for the fractures, which, in turn, is predicted by recurrent neural networks using the concentration histories at the fracture-matrix boundary. Specifically, we consider a meter-scale fracture network embedded in sandstones, where the smallest feature is at the micron scale. Here the proposed framework allows us to span five orders of magnitude in length scales by capturing mineral precipitation in the altered layer of the rock matrix at the pore scale across the entire meter-scale fracture network.
42 ENGINEERING↗
Interlaminar fracture toughness: Three-dimensional finite element modeling for end-notch and mixed-mode flexure
A computational procedure is described for evaluating End-Notch-Flexure (ENF) and Mixed-Mode-Flexure (MMF) interlaminar fracture toughness in unidirectional fiber composites. The procedure consists of a three-dimensional finite element analysis in conjunction with the strain energy release rate concept and with composite micromechanics. The procedure is used to analyze select cases of ENF and MMF. The strain energy release rate predicted by this procedure is in good agreement with limited experimental data. The procedure is used to identify significant parameters associated with interlaminar fracture toughness. It is also used to determine the critical strain energy release rate and its attendant crack length in ENF and/or MMF. This computational procedure has considerable versatility/generality and provides extensive information about interlaminar fracture toughness in fiber composites.
Murthy, P. L. N.↗
Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants
Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.
42 ENGINEERING↗
An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem
The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.
36 MATERIALS SCIENCE↗
CabanaPD
CabanaPD is a meshfree peridynamics application built with Cabana and Kokkos. Kokkos enables performance portability across hardware architectures and Cabana provides particle capabilities including multi-node MPI support. The main components of CabanaPD are particle initialization, neighbor list generation, force and energy computation, time integration, and multi-node particle communication. In addition, options for creating pre-cracked regions and particle boundary conditions are available. CabanaPD currently enables two common force models: prototype microelastic brittle (PMB) and linear peridynamic solid (LPS). For both of these models, versions with and without fracture as well as linearized model options are available. CabanaPD is designed to be extensible for addition of more complex force models, boundary conditions, etc.
Reeve, Sam↗
Fracture analysis of stiffened panels under biaxial loading with widespread cracking
An elastic-plastic finite-element analysis with a critical crack-tip opening angle (CTOA) fracture criterion was used to model stable crack growth and fracture of 2024-T3 aluminum alloy (bare and clad) panels for several thicknesses. The panels had either single or multiple-site damage (MSD) cracks subjected to uniaxial or biaxial loading. Analyses were also conducted on cracked stiffened panels with single or MSD cracks. The critical CTOA value for each thickness was determined by matching the failure load on a middle-crack tension specimen. Comparisons were made between the critical angles determined from the finite-element analyses and those measured with photographic methods. Predicted load-against-crack extension and failure loads for panels under biaxial loading, panels with MSD cracks, and panels with various numbers of stiffeners were compared with test data whenever possible. The predicted results agreed well with the test data even for large-scale plastic deformations. The analyses were also able to predict stable tearing behavior of a large lead crack in the presence of MSD cracks. The analyses were then used to study the influence of stiffeners on residual strength in the presence of widespread fatigue cracking. Small MSD cracks were found to greatly reduce the residual strength for large lead cracks even for stiffened panels.
Newman, J. C., Jr.↗
Fracture analysis of stiffened panels under biaxial loading with widespread cracking
An elastic-plastic finite-element analysis with a critical crack-tip-opening angle (CTOA) fracture criterion was used to model stable crack growth and fracture of 2024-T3 aluminum alloy (bare and clad) panels for several thicknesses. The panels had either single or multiple-site damage (MSD) cracks subjected to uniaxial or biaxial loading. Analyses were also conducted on cracked stiffened panels with single or MSD cracks. The critical CTOA value for each thickness was determined by matching the failure load on a middle-crack tension specimen. Comparisons were made between the critical angles determined from the finite-element analyses and those measured with photographic methods. Predicted load-against-crack extension and failure loads for panels under biaxial loading, panels with MSD cracks, and panels with various number of stiffeners were compared with test data, whenever possible. The predicted results agreed well with the test data even for large-scale plastic deformations. The analyses were also able to predict stable tearing behavior of a large lead crack in the presence of MSD cracks. The analyses were then used to study the influence of stiffeners on residual strength in the presence of widespread fatigue cracking. Small MSD cracks were found to greatly reduce the residual strength for large lead cracks even for stiffened panels.
Newman, J. C., Jr.↗
Mechanical characterization and modeling of non-linear deformation and fracture of a fiber reinforced metal matrix composite
The nonlinear anisotropic mechanical behavior of an aluminum alloy metal matrix composite reinforced with continuous alumina fibers was determined experimentally. The mechanical behavior of the composite were modeled by assuming that the composite has a periodical microstructure. The resulting unit cell problem was solved with the finite element method. Excellent agreement was found between theoretically predicted and measured stress-strain responses for various tensile and shear loadings. The stress-strain responses for transverse and inplane shear were found to be identical and this will provide a simplification of the constitutive equations for the composite. The composite has a very low ductility in transverse tension and a limited ductility in transverse shear that was correlated to high hydrostatic stresses that develop in the matrix. The shape of the initial yield surface was calculated and good agreement was found between the calculated shape and the experimentally determined shape.
Jansson, S.↗
Progressive fracture of fiber composites
Refined models and procedures are described for determining progressive composite fracture in graphite/epoxy angleplied laminates. Lewis Research Center capabilities are utilized including the Real Time Ultrasonic C Scan (RUSCAN) experimental facility and the Composite Durability Structural Analysis (CODSTRAN) computer code. The CODSTRAN computer code is used to predict the fracture progression based on composite mechanics, finite element stress analysis, and fracture criteria modules. The RUSCAN facility, CODSTRAN computer code, and scanning electron microscope are used to determine durability and identify failure mechanisms in graphite/epoxy composites.
Irvin, T. B.↗
Progressive fracture of fiber composites
Refined models and procedures are described for determining progressive composite fracture in graphite/epoxy angleplied laminates. Unique Lewis Research Center capabilities are utilized including the Real-Time Ultrasonic C-San (RUSCAN) experimental facility and the Composite Durability Structural Analysis (CODSTRAN) computer code. CODSTRAN is used to predict the fracture progression based on composite mechanics, finite element stress analysis, and fracture criteria modules. The RUSCAN facility, CODSTRAN computer code, and scanning electron microscope are used to determine durability and identify failure mechanisms in graphite/epoxy coomposites. Results indicate that RUSCAN/CODSTRAN is an effective method of Studying progressive fracture of composites.
Irvine, T. B.↗