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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Anomalous Transport in Three-Dimensional Discrete Fracture Networks: Interplay Between Aperture Heterogeneity and Injection Modes

In this work, we study how the interplay between fracture aperture heterogeneity and tracer injection mode controls fluid flow and tracer transport in three-dimensional (3D) discrete fracture networks (DFNs). The direct 3-D DFN simulations show that tracer injection mode has substantial effects on tracer spreading across all levels of aperture heterogeneity. The key controlling factor for effective transport is the initial Lagrangian velocity distribution, which is determined by the interplay between injection mode and aperture heterogeneity. The fundamental difference between initial Lagrangian velocity distribution and domain-scale Eulerian velocity distribution plays a vital role in determining anomalous transport. We effectively capture the observed anomalous transport using an upscaled transport model that incorporates initial velocity distribution, stationary velocity distribution, velocity correlation length, and average advective tortuosity. With the upscaled transport model, we accurately capture the evolution of Lagrangian velocity distribution and predict longitudinal spreading in 3-D DFN.

58 GEOSCIENCES↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Tip radius preservation for high resolution imaging in amplitude modulation atomic force microscopy

The acquisition of high resolution images in atomic force microscopy (AFM) is correlated to the cantilever's tip shape, size, and imaging conditions. In this work, relative tip wear is quantified based on the evolution of a direct experimental observable in amplitude modulation atomic force microscopy, i.e., the critical amplitude. We further show that the scanning parameters required to guarantee a maximum compressive stress that is lower than the yield/fracture stress of the tip can be estimated via experimental observables. In both counts, the optimized parameters to acquire AFM images while preserving the tip are discussed. The results are validated experimentally by employing IgG antibodies as a model system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

This is a final technical report for the project: Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems (EGS). The goal is to demonstrate the feasibility of foam fracturing in EGS applications. The project, led by Oak Ridge National Laboratory (ORNL), was conducted in collaboration with Temple University. The report describes the research activities with Task 1 at ORNL: foam fracturing testing system development and experimental study on foam fracturing, and Task 2 at Temple University: foam testing and foam characterization. Main findings are: 1. A foam fracturing test system has been developed at ORNL, which can be used to perform foam fracturing under pressure up to 6,000 psi. The system monitors foam density during fracturing online and is capable of testing materials in both monotonic and cyclic (up to 50 Hz) injections. 2. Foam fracturing tests were carried out on Charcoal black granite specimens with a blind borehole to the middle length. Two diameters of blind borehole were tested; G2 series: 9.53 mm and G3 series: 4.76 mm. N2-in-water foam was used with AOS as a surfactant. 3. There was a hole-size effect on fracture initiation pressure. The effect is smaller in the case of foam, which was influenced by the high penetrability of gas in foam. Breakdown pressure showed a behavior just as that of fracture pressure; namely an increased value for small hole samples, while the effect in water fracture was more impressive than in foam fracture. 4. Water mass was reduced in foam fracturing within similar range of breakdown pressures. In G2 series, it was decreased from 10.44 g for water fracturing to 5.17 g, representing more than 50% water reduction. Therefore, there is the potential to reduce water use in EGS stimulation through foam fracturing. 5. Use of cyclic injection has the potential to reduce the breakdown pressure and seismicity in EGS application. Experiments using 4-s cycle period found that specimens can be fractured with a low number of cycles. The fatigue pressure was approximately 64 - 77% of monotonic breakdown pressure for water fracturing and 58 - 94% of the breakdown pressure for foam fracturing. 6. A foam stability testing system has been developed that can test foam at 220 Deg C to 2,000 psi. Tested components of candidate foams included two gases: N2 and CO2; 4 surfactants: AOS, SDS, NP-40 and CTAC; 5 stabilizing agents: guar, bentonite clay, borate salt, silica NPs, and GO. 7. N2 and AOS provided the most stable performance over the tested ranges. Furthermore, the AOS foam with stabilizing agents of guar and borate salt (crosslinker) offered the highest half-life of 20 minutes at 200 Deg C and 1,000 psi. 8. Arrhenius equation and modified power law have been demonstrated to fit well the half-time vs. temperature and pressure data, respectively. These relations can be useful to provide the suggestion for future foam stability study. This submission contains the supporting data developed during the project: 1) A final technical report 2) Granite fracturing data in monotonic and cyclic injections with water and N2 foam Foam performance data in various temperatures and pressures, including half-time, is submitted separately.

15 GEOTHERMAL ENERGY↗
DOE OSTI2017

Modeling radionuclide migration from underground nuclear explosions

The travel time of radionuclide gases to the ground surface in fracture rock depends on many complex factors. Numerical simulators are the most complete repositories of knowledge of the complex processes governing radionuclide gas migration to the ground surface allowing us to verify conceptualizations of physical processes against observations and forecast radionuclide gas travel times to the ground surface and isotopic ratios

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1976

Program definition and preliminary results of a long-term evaluation program of advanced composites for supersonic cruise aircraft applications

Advanced composites will play a key role in the development of the technology required for the design and fabrication of future supersonic vehicles. A program is in progress to determine the time-temperature-stress capabilities of several high-temperature composite materials. Tests included in this study are thermal aging, ambient aging, fatigue, creep, fracture, tensile, and real-time flight-simulation exposure. The program has two parts. The first includes all the material-property determinations and aging and simulation exposures up through 10,000 h. The second continues these tests up to 50,000 cumulative h. The program is currently in the midst of the 10,000-h phase. Some of the results obtained to date and the status of the various tests are discussed.

Kerr, J. R.↗
NASA NTRS1987

Energy dissipation associated with crack extension in an elastic-plastic material

Crack extension in elastic-plastic material involves energy dissipation through the creation of new crack surfaces and additional yielding around the crack front. An analytical procedure, using a two-dimensional elastic-plastic finite element method, was developed to calculate the energy dissipation components during a quasi-static crack extension. The fracture of an isotropic compact specimen was numerically simulated using the critical crack-tip-opening-displacement (CTOD) growth criterion. Two specimen sizes were analyzed for three values of critical CTOD. Results from the analyses showed that the total energy dissipation rate consisted of three components: the crack separation energy rate, the plastic energy dissipation rate, and the residual strain energy rate. All three energy dissipation components and the total energy dissipation rate initially increased with crack extension and finally reached constant values.

Shivakumar, K. N.↗
NASA NTRS1987

Energy dissipation associated with crack extension in an elastic-plastic material

Crack extension in elastic-plastic material involves energy dissipation through the creation of new crack surfaces and additional yielding around the crack front. An analytical procedure, using a two-dimensional elastic-plastic finite element method, was developed to calculate the energy dissipation components during a quasi-static crack extension. The fracture of an isotropic compact specimen was numerically simulated using the critical crack-tip-opening-displacement (CTOD) growth criterion. Two specimen sizes were analyzed for three values of critical CTOD. Results from the analysis showed that the total energy dissipation rate consisted of three components: the crack separation energy rate, the plastic energy dissipation rate, and the residual strain energy rate. All three energy dissipation components and the total energy dissipation rate initially increased with crack extension and finally reached constant values.

Shivakumar, K. N.↗
NASA NTRS1988

An elastic-plastic finite-element analysis of the J-resistance curve using a CTOD criterion

The fracture process in compact and bend specimens was simulated using a two-dimensional finite-element analysis of the J-resistance (J-R) curve and a CTOD criterion. The J-R curves were calculated from the numerical results for each specimen type using several different methods. In general, the J-R curves obtained for the bend specimens were found to be higher than those for the compact specimens, especially beyond maximum load. However, below the maximum load, the modified deformation theory of plasticity and the contour-integral J(r) method resulted in very similar J-R curves for both specimen types.

Newman, J. C., Jr.↗
NASA NTRS1992

High-performance computing for flight vehicles; Proceedings of the Symposium, Washington, Dec. 7-9, 1992

The present conference discusses high-performance computing systems for flight vehicles, large-scale simulations on high-performance flight computers and software, multidisciplinary and design/optimization applications of computers, computational electromagnetics and acoustics, the simulation of aircraft powerplant turbomachinery and reacting flows, and flow calculations on parallel machines. Also discussed are direct flow simulation Monte Carlo methods, structural mechanics sensitivity and fracture calculations on parallel machines, grid-generation and advanced algorithms for CFD, advanced solid-mechanics and structures applications, and advancements in flow visualization technology and neural networks.

Noor, Ahmed K.↗
NASA NTRS1999

Damage Tolerance of Large Shell Structures

Progressive damage and fracture of large shell structures is investigated. A computer model is used for the assessment of structural response, progressive fracture resistance, and defect/damage tolerance characteristics. Critical locations of a stiffened conical shell segment are identified. Defective and defect-free computer models are simulated to evaluate structural damage/defect tolerance. Safe pressurization levels are assessed for the retention of structural integrity at the presence of damage/ defects. Damage initiation, growth, accumulation, and propagation to fracture are included in the simulations. Damage propagation and burst pressures for defective and defect-free shells are compared to evaluate damage tolerance. Design implications with regard to defect and damage tolerance of a large steel pressure vessel are examined.

Minnetyan, L.↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗
NASA NTRS2022

A Cohesive Element Formulation for Simulation of General Mode I, II and III Delamination

Cohesive element formulations proposed thus far have assumed that Mode II and Mode III fracture processes, and hence traction-separation laws, are identical. However, experimental observations suggest that such an assumption may not hold in general. The approach described herein consists of a cohesive element formulation in which the assumption of identical Mode II and Mode III fracture processes is not required. This assumption is alleviated by estimating the normal to the delamination front based on the gradient of the displacement jumps determined within each element. The estimated normal is used to rotate the cohesive element coordinate system. Determining the displacement jumps in the rotated coordinate system enables the distinction between Mode II sliding, and Mode III tearing. This information is subsequently used to extend a recently proposed cohesive formulation, resulting in an approach that can simulate mixed-mode I/II/III fracture via prescribing piecewise-linear traction-separation laws combined with a 3D mixed-mode fracture criteria.

Delamination↗
DOE OSTI2016

Complex fracture nucleation and evolution with nonlocal elastodynamics

In this study, a mechanical model is introduced for predicting the initiation and evolution of complex fracture patterns without the need for damage variable or law. The model, a continuum variant of Newton's second law, uses integral rather than partial differential operators where the region of integration is over finite domain. The force interaction is derived from a novel nonconvex strain energy density function, resulting in a nonmonotomic material model. The resulting equation of motion is proved to be mathematically well-posed. The model has the capacity to simulate nucleation and growth of multiple, mutually interacting dynamic fractures. In the limit of zero region of integration, the model reproduces the classic Griffith model of brittle fracture. The simplicity of the formulation avoids the need for supplemental kinetic relations that dictate crack growth or the need for an explicit damage evolution law.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Toward better hydraulic fracturing fluids and their application in energy production: A review of sustainable technologies and reduction of potential environmental impacts

Recent advances in hydraulic fracturing, in conjunction with horizontal drilling, have enabled large-scale extraction of natural gas and oil from shale formations. Despite its advances and enormous economic benefits, opportunities remain to increase hydraulic fracturing efficiency and minimize potential environmental impacts. This review specifically examines three key themes associated with development and utilization of hydraulic fracturing fluids: 1) characteristics and behavior of fracturing fluids, 2) understanding and predicting migration and fate of fracturing fluids, 3) technologies to reduce environmental impact of fracturing fluids. Furthermore, the paper discusses key and new techniques and findings on rheology of hydrogel-based fluids, high fidelity simulation of propagation transport, potential environmental impacts, geosynthetics in mitigating contamination, and greener fracturing fluids. It is indicated that future development relies on advances in understanding of physical processes, modeling capabilities, and monitoring techniques.

02 PETROLEUM↗
DOE OSTI2018

Assessing the Geomechanical Response of CO 2 Disposal in Flood Basalt Reservoirs

In this work, we complete a multi-scale investigation of the physical processes governing geologic CO 2 sequestration in basalt reservoirs. This work integrates (i) laboratory measurements of relative permeability in a basalt fracture, (ii) theoretical modeling to understand how relative permeability uncertainty affects geomechanical performance attributes of basalt reservoirs, (iii) outcrop-scale simulation to identify the characteristics of vertical CO 2 flow in a basalt fracture network, and (iv) regional-scale assessment of permeability architecture in the Columbia River Basalt Group (CRBG). These individual studies are combined into a site-scale stochastic model of geologic CO 2 sequestration in the Columbia River Plateau, which shows how bulk and relative permeability variability affects geomechanical reservoir integrity and storage capacity of flood basalt reservoirs. This project resulted in a novel laboratory technique that combines X-ray CT and micro-PET scanning to calculate the relative permeability curves for a horizontal fracture in basalt. We also discovered that permeability on the Columbia River Plateau exhibits directionally-dependent spatial correlation, which our analysis suggests is a result of loading-induced subsidence. This analysis further discovered that permeability in the CRBG may increase at depth due to bending moment stresses. In the context of CO 2 sequestration, this result implies that storage capacity and injectivity may actually increase at depths beyond ~1 km. Our outcrop-scale investigation found that buoyant, free-phase CO 2 tends to accumulate at fracture intersections, which inhibits vertical flow over 10+ year timescales and may facilitate both physical and chemical trapping as fluid mobility decreases. At the site-scale, we find that in flood basalt reservoirs the spatial configuration of reservoir permeability exhibits substantial control on CO 2 injectivity. We utilized known borehole geology from the Wallula Basalt Sequestration Pilot Project to develop 50 equally probable permeability configurations and found that a single injection well operating at 95% of the borehole breakout pressure can deliver 0.12 to 2.0 million metric tons of CO 2 per year (MMT CO 2 per year). This maximum injection rate can support a 1,000 MW gas-fired power plant with a single injection well.

58 GEOSCIENCES↗
DOE OSTI2019

Complex Fracture Nucleation and Evolution with Nonlocal Elastodynamics

A mechanical model is introduced for predicting the initiation and evolution of complex fracture patterns without the need for a damage variable or law. The model, a continuum variant of Newton’s second law, uses integral rather than partial differential operators where the region of integration is over finite domain. The force interaction is derived from a novel nonconvex strain energy density function, resulting in a nonmonotonic material model. The resulting equation of motion is proved to be mathematically well-posed. The model has the capacity to simulate nucleation and growth of multiple, mutually interacting dynamic fractures. In the limit of zero region of integration, the model reproduces the classic Griffith model of brittle fracture. As a result, the simplicity of the formulation avoids the need for supplemental kinetic relations that dictate crack growth or the need for an explicit damage evolution law.

36 MATERIALS SCIENCE↗