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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Influence of injection mode on transport properties in kilometer-scale three-dimensional discrete fracture networks

We investigate how the choice of injection mode impacts transport properties in kilometer-scale three-dimensional discrete fracture networks (DFN). The choice of injection mode, resident and flux-weighted, is designed to mimic different physical phenomena. It has been hypothesized that solute plumes injected under resident conditions evolve to behave similarly to solutes injected under flux-weighted conditions. Previously, computational limitations have prohibited the large-scale simulations required to investigate this hypothesis. We investigate this hypothesis by using a high-performance DFN suite, dfnWorks, to simulate flow in kilometer-scale three-dimensional DFNs based on fractured granite at the Forsmark site in Sweden, and adopt a Lagrangian approach to simulate transport therein. Results show that after traveling through a pre-equilibrium region, both injection methods exhibit linear scaling of the first moment of travel time and power law scaling of the breakthrough curve with similar exponents, slightly larger than 2. Lastly, the physical mechanisms behind this evolution appear to be the combination of in-network channeling of mass into larger fractures, which offer reduced resistance to flow, and in-fracture channeling, which results from the topology of the DFN.

42 ENGINEERING↗
DOE OSTI2023

Clinical effects of 3D printing-assisted posterolateral incision in the treatment of ankle fractures involving the posterior malleolus

Objective To explore the clinical outcomes of a 3D printing-assisted posterolateral approach for the treatment of ankle fractures involving the posterior malleolus. Methods A total of 51 patients with ankle fractures involving the posterior malleolus admitted to our hospital from January 2018 to December 2019 were selected. The patients were divided into 3D printing group (28 cases) and control group (23 cases). 3D printing was performed for ankle fractures, followed by printing of a solid model and simulation of the operation on the 3D model. The operation was then performed according to the preoperative plan, including open reduction and internal fixation via the posterolateral approach with the patient in the prone position. Routine x-ray and CT examinations of the ankle joint were performed, and ankle function was evaluated using the American Foot and Ankle Surgery Association (AOFAS) ankle-hindfoot score. Results All patients underwent x-ray and CT examinations. All fractures healed clinically, without loss of reduction or failure of internal fixation. Good clinical effects were achieved in both groups of patients. The operation time, intraoperative blood loss and intraoperative fluoroscopy frequency in the 3D printing group were significantly less than those in the control group ( p < 0.05). There was no significant difference between the two groups in the anatomical reduction rate of fractures or the incidence of surgical complications ( p > 0.05). Conclusion The 3D printing-assisted posterolateral approach is effective in the treatment of ankle fractures involving the posterior malleolus. The approach can be well planned before the operation, is simple to perform, yields good fracture reduction and fixation, and has good prospects for clinical application.

Zheng, Hongming↗
DOE OSTI2021

Field Experiments and Model Validation: The EGS Collab Project

Implementing enhanced geothermal systems (EGS) will require improvements in understanding stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The US Department of Energy (DOE) is addressing these and other challenges at multiple scales. The EGS Collab project, addressed here, is performing tests and modeling at the 10 m scale. The FORGE project is performing tests at the full reservoir scale. The EGS Collab team created an underground testbed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently building a second testbed at SURF at a depth of about 1.25 km aimed at investigating shear stimulation (Experiment 2). In Experiment 1, we drilled eight boreholes in a well-characterized phyllite and installed geophysical sensors in six of them to create a well-instrumented testbed to allow careful monitoring of stimulation events and flow tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results both in near-real-time and in detailed process studies. Stimulations performed in this testbed allowed quantification of processes occurring during stimulation and the examination of dynamic flow occurrences. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Our second testbed, targeted at shear stimulation, is currently being built at the SURF Facility at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1.

artificial intelligence↗
DOE OSTI2019

Finite element analyses of single particle crushing tests incorporating computed tomography imaging and damage mechanics

Fracture of individual grains within a granular assembly affect the strength and deformation behavior of granular materials under large stress states. Interparticle forces, leading to contact stresses and ultimately fracture initiation, are influenced by particle morphology. A numerical method that incorporates both particle fracture and morphology can provide a more accurate model of a granular material’s macro-scale response. This research investigates an approach to modeling fracture initiation and propagation within individual grains, utilizing high resolution X-ray computed tomography (CT) imaging to consider grain morphology and an explicit finite element code optimized for high performance computing which incorporates damage mechanics. Using precisely measured force-displacement curves from single particle crushing tests of manufactured quartz spheres, the method is validated, and the effects of fracture properties on grain fragmentation are determined. This approach is then used to simulate single particle crushing of Ottawa sand, with grain morphology incorporated through discretized CT images. Here the effect of grain orientation on the force and applied displacement at catastrophic splitting is also explored. The proposed discrete particle based finite element framework can be applied to modeling granular assemblies, considering the effect of individual particle shape and fracture on the assembly’s deformation response through well calibrated numerical simulations.

42 ENGINEERING↗
DOE OSTI2014

Solute embrittlement of SiC

The energies and stresses associated with the decohesion of β-SiC in the presence of mobile Pd and Ag impurities are studied from first principles. Density functional theory calculations are parameterized with a generalized cohesive zone model and are analyzed within a thermodynamic framework that accounts for realistic boundary conditions in the presence of mobile impurities. We find that Pd impurities will embrittle SiC when Pd is in equilibrium with metallic Pd precipitates. Our thermodynamic analysis predicts that Pd embrittles SiC by substantially reducing the maximum stress of decohesion as a result of a phase transition between decohering planes involving an influx of Pd atoms. The methods presented in this work can be applied to study the thermodynamics of decohesion of SiC in other aggressive environments containing oxygen and water, for example, and yield environment dependent cohesive zone models for use in continuum approaches to study crack propagation and fracture.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

A Stochastic Study of Flow Anisotropy and Channelling in Open Rough Fractures

The quantification of fluid flow in rough fractures is of high interest for reservoir engineering, especially for deep geothermal applications. Herein, rough self-affine fractures are stochastically generated with incremental shear displacement and geometrically described by two aperture definitions, the vertical aperture a vert and the effective aperture a eff . In order to compare their effect on fracture flow, such as anisotropy and channelling, Local Cubic Law (LCL) model-based 2D fluid flow is simulated. The particularity of this approach is the combination of a stochastic generation of self-affine fractures with a statistical analysis (560 individual realizations) of the impact of the LCL’s aperture constraint on fracture flow. The results show that aperture definition affects the quantitative interpretation of flow anisotropy and channeling as well as the aperture distribution of the fractures with shearing. Higher values of mean aperture for a given fracture are found using a vert , whereas the aperture standard deviation is larger with a eff . In addition, flow anisotropy is significantly sensitive to aperture definition for small shear displacements and shows a relative higher dispersion with a eff . Thus, LCL prediction models based on a vert are expected to lead to higher dispersion of anisotropy results with a higher uncertainty (factor ~ 2). Realizations based on avert lead to an enhanced clustering of high flow rates for higher shearing displacements. This channeling development results in higher total flow rates for these simulations. Finally, these findings support the direct calibration of pre-existing LCL anisotropy simulations based on avert towards more representative results using a eff .

58 GEOSCIENCES↗
DOE OSTI2021

Enabling Particulate Materials Processing Science for High-Consequence, Small-Lot Precision Manufacturing

This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Facilitating the exploitation of ERTS imagery using snow enhancement techniques

The author has identified the following significant results. Analysis of all available (Gemini, Apollo, Nimbus, NASA aircraft) small scale snow covered imagery has been conducted to develop and refine snow enhancement techniques. A detailed photographic interpretation of ERTS-simulation imagery covering the Feather River/Lake Tahoe area was completed and the 580-680nm. band was determined to be the optimum band for fracture detection. ERTS-1 MSS bands 5 and 7 are best suited for detailed fracture mapping. The two bands should provide more complete fracture detail when utilized in combination. Analysis of early ERTS-1 data along with U-2 ERTS simulation imagery indicates that snow enhancement is a viable technique for geological fracture mapping. A wealth of fracture detail on snow-free terrain was noted during preliminary analysis of ERTS-1 images 1077-15005-6 and 7, 1077-15011-5 and 7, and 1079-15124-5 and 7. A direct comparison of data yield on snow-free versus snow-covered terrain will be conducted within these areas following receipt of snow-covered ERTS-1 imagery.

Wobber, F. J.↗
DOE OSTI2016

A phase-field approach to model multi-axial and microstructure dependent fracture in nuclear grade graphite

The fracture behavior of nuclear grade graphites is strongly influenced by underlying microstructural features such as the character of filler particles, and the distribution of pores and voids. These microstructural features influence the crack nucleation and propagation behavior, resulting in quasi-brittle fracture with a tortuous crack path and significant scatter in measured bulk strength. This paper uses a phase-field method to model the microstructural and multi-axial fracture in H-451, a historic variant of nuclear graphite that provides the basis for an idealized study on a legacy grade. The representative volume elements are constructed from randomly located pores with random size obtained from experimentally determined log-normal distribution. The representative volume elements are then subjected to simulated multi-axial loading, and a reasonable agreement of the resulting fracture stress with experiments is obtained. Finally, quasi-brittle stress-strain evolution with a tortuous crack path is also observed from the simulations and is consistent with experimental results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Carbide Particle Redistribution in U-10Mo Alloy during Hot and Cold Rolling Processes: Synthetic Microstructure-Based Finite-Element Approach

The influence of initial carbide particle distribution features on particle fracture and the evolution of stringer distribution within an alloy of uranium with 10 weight percent molybdenum during cold rolling was investigated in this study. This work will provide guidance to casting optimization work and identify actionable process changes to eliminate the surface features directly attributed to carbides. In a previous study, various synthetic microstructure-based finite-element (FE) models were generated with consideration of different parameters (i.e., particle shape, particle volume fraction, and particle distribution), and used to examine the effects of hot rolling on the microstructure evolution. In this study, microstructures after ~80% hot rolling were used as the initial microstructures for the subsequent cold rolling simulation. Plane-strain compression loading was employed in the FE models to simulate cold rolling with reductions from ~80% to ~94%. Stringer analyses with the consideration of particle fracture were then performed with the simulated cold rolling results. The results of void evolution analyses show that microstructures with elliptical particles incorporate a much larger (almost double) void volume fraction than those with circular particles, indicating that the elliptical particles are fractured more easily than the circular ones. The results from the fracture mode analyses show that, under the adopted simple compression loading assumption, almost all the particles (i.e., >99%) were fractured by tensile separation. The results of stringer analyses show that, during cold rolling, the stringer volume fractions increase substantially due to the particle fracture, especially for the microstructures with the elliptical particles, and that the stringer volume fractions from stringers comprising many particles tend to increase faster. The maximum stringer lengths for the microstructures with circular particles tend to remain constant or decrease during cold rolling, whereas those for elliptical particles tend to slightly increase after a slight decrease. The number of particles at the surface increases very rapidly due to particle fracture during cold rolling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hydraulic and Thermal Stimulation Program at Raft River Idaho, A DOE EGS

A Department of Energy Enhanced Geothermal System (EGS) stimulation program has injected over 254 million gallons of water into the well RRG-9 ST1 since the summer of 2013. Three major stimulations have been conducted during the program increasing injection flow rates from less than 20 gpm to 550 gpm. Geologic, water chemistry, microseismic activity, and borehole imaging data have been used to develop a conceptual model describing possible flow paths of this injected water. This model contains two major fracture zones one of which intersects the RRG-9 ST1 wellbore. Modified Hall and injectivity index plots constructed using injection flow rates, surface temperatures, and wellhead pressures show steady improvement in the injectivity of the well. Here, the injectivity index has risen from 0.15 gpm/psi to 2.0 gpm/ psi. A pressure falloff test conducted on April 28, 2015 indicates a reservoir permeability of 1,220 md and -5.38 skin factor. The well stimulation program was simulated numerically using an Idaho National Laboratory reservoir simulation code, FALCON. These simulations show a significant increase in the permeability of connecting fracture pathways after each stimulation event.

Enhanced geothermal system↗
DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Accurate Simulation of Acoustic Emission Sources in Composite Plates

Acoustic emission (AE) signals propagate as the extensional and flexural plate modes in thin composite plates and plate-like geometries such as shells, pipes, and tubes. The relative amplitude of the two modes depends on the directionality of the source motion. For source motions with large out-of-plane components such as delaminations or particle impact, the flexural or bending plate mode dominates the AE signal with only a small extensional mode detected. A signal from such a source is well simulated with the standard pencil lead break (Hsu-Neilsen source) on the surface of the plate. For other sources such as matrix cracking or fiber breakage in which the source motion is primarily in-plane, the resulting AE signal has a large extensional mode component with little or no flexural mode observed. Signals from these type sources can also be simulated with pencil lead breaks. However, the lead must be fractured on the edge of the plate to generate an in-plane source motion rather than on the surface of the plate. In many applications such as testing of pressure vessels and piping or aircraft structures, a free edge is either not available or not in a desired location for simulation of in-plane type sources. In this research, a method was developed which allows the simulation of AE signals with a predominant extensional mode component in composite plates requiring access to only the surface of the plate.

Prosser, W. H.↗
NASA NTRS2010

Modeling Progressive Damage Using Local Displacement Discontinuities Within the FEAMAC Multiscale Modeling Framework

A method for performing progressive damage modeling in composite materials and structures based on continuum level interfacial displacement discontinuities is presented. The proposed method enables the exponential evolution of the interfacial compliance, resulting in unloading of the tractions at the interface after delamination or failure occurs. In this paper, the proposed continuum displacement discontinuity model has been used to simulate failure within both isotropic and orthotropic materials efficiently and to explore the possibility of predicting the crack path, therein. Simulation results obtained from Mode-I and Mode-II fracture compare the proposed approach with the cohesive element approach and Virtual Crack Closure Techniques (VCCT) available within the ABAQUS (ABAQUS, Inc.) finite element software. Furthermore, an eccentrically loaded 3-point bend test has been simulated with the displacement discontinuity model, and the resulting crack path prediction has been compared with a prediction based on the extended finite element model (XFEM) approach.

Ranatunga, Vipul↗
DOE OSTI2017

Caprock integrity susceptibility to permeable fracture creation

Caprock leakage is of crucial concern for environmentally and economically sustainable development of carbon dioxide sequestration and utilization operations. One potential leakage pathway is through fractures or faults that penetrate the caprock. In this study, we investigate the permeability induced by fracturing initially intact Marcellus shale outcrop specimens at stressed conditions using a triaxial direct-shear method. Measurements of induced permeability, fracture geometry, displacement, and applied stresses were all obtained at stressed conditions to investigate the coupled processes of fracturing and fluid flow as may occur in the subsurface. Fracture geometry was directly observed at stressed conditions using X-ray radiography video. Numerical simulation was performed to evaluate the stress distribution developed in the experiments. Our experiments show that permeability induced by fracturing is strongly dependent on the stresses at which the fractures are created, the magnitude of shearing displacement, and the duration of flow. The strongest permeability contrast was observed when comparing specimens fractured at low stress to others fractured at higher stress. Measureable fracture permeability decreased by up to 7 orders of magnitude over a corresponding triaxial confining stress range of 3.5 MPa to 30 MPa. These results show that increasing stress, depth, and time are all significant permeability inhibitors that may limit potential leakage through fractured caprock.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

iMatTOUGH: An open-source Matlab-based graphical user interface for pre- and post-processing of TOUGH2 and iTOUGH2 models

TOUGH2 and iTOUGH2 are powerful models that simulate the heat and fluid flows in porous and fracture media, and perform parameter estimation, sensitivity analysis and uncertainty propagation analysis. However, setting up the input files is not only tedious, but error prone, and processing output files is time consuming. Here, we present an open source Matlab-based tool (iMatTOUGH) that supports the generation of all necessary inputs for both TOUGH2 and iTOUGH2 and visualize their outputs. The tool links the inputs of TOUGH2 and iTOUGH2, making sure the two input files are consistent. It supports the generation of rectangular computational mesh, i.e., it automatically generates the elements and connections as well as their properties as required by TOUGH2. The tool also allows the specification of initial and time-dependent boundary conditions for better subsurface heat and water flow simulations. The effectiveness of the tool is illustrated by an example that uses TOUGH2 and iTOUGH2 to estimate soil hydrological and thermal properties from soil temperature data and simulate the heat and water flows at the Rifle site in Colorado.

58 GEOSCIENCES↗
DOE OSTI2019

Effect of Pre-Weld Heat Treatment on Microstructure and Creep Strength of ICHAZ in Grade 91 Steel

In this work, Grade 91 base metal was normalized with two different cooling conditions, water quenching and air cooling. The intercritical welding thermal cycle simulated by using the Gleeble system were applied to the heat-treated base metals to reproduce the intercritical heat-affected zone (ICHAZ). Microstructure, including precipitate and tempered martensite, of the base metal and ICHAZ was carefully characterized with advanced microscopy techniques. Creep strength of the simulated ICHAZs was evaluated at a high creep temperature of 650 °C with a stress of 100 MPa. A correlation between the microstructure and creep behavior of the ICHAZ is built to understand creep rupture mechanisms in the ICHAZ. The results show that more coarse carbides precipitated along the prior austenite grain boundaries after tempering in the water-quenched base metal. These carbides cannot be fully dissolved by the intercritical welding thermal cycle. The simulated ICHAZ generated from the water-quenched base metal underwent a higher grain recovery and growth during the typical post-weld heat treatment at 760 °C. Both simulated ICHAZs from two pre-weld heat treated base metals exhibit an extremely low creep resistance with a typical ductile fracture during creep testing. The faster grain growth and precipitate coarsening in the ICHAZ, simulated from the water-quenched base metal, made it even worse.

Wang, Yiyu↗
NASA NTRS1991

Microfracture in high temperature metal matrix crossply laminates

Microfracture (fiber/matrix fracture, interphase debonding, and inter-ply delamination) in high temperature metal matrix composites (HTMMC), subjected to both mechanical and thermal loading, is computationally simulated. A crossply 0.3 fiber volume ratio SiC/Ti15 composite with 0/90/0 lay-up is evaluated for microfracture using a multicell finite element model. A computational simulation procedure based on strain energy release rates is used to predict the fracture process and establish the hierarchy of fracture modes. Microfracture results for various loading cases are presented and discussed.

Mital, Subodh K.↗