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

Effects of lithium insertion on thermal conductivity of silicon nanowires

Recently, silicon nanowires (SiNWs) have been applied as high-performance Li battery anodes, since they can overcome the pulverization and mechanical fracture during lithiation. Although thermal stability is one of the most important parameters that determine safety of Li batteries, thermal conductivity of SiNWs with Li insertion remains unclear. In this letter, using molecular dynamics simulations, we study room temperature thermal conductivity of SiNWs with Li insertion. It is found that compared with the pristine SiNW, there is as much as 60% reduction in thermal conductivity with 10% concentration of inserted Li atoms, while under the same impurity concentration the reduction in thermal conductivity of the mass-disordered SiNW is only 30%. With lattice dynamics calculations and normal mode decomposition, it is revealed that the phonon lifetimes in SiNWs decrease greatly due to strong scattering of phonons by vibrational modes of Li atoms, especially for those high frequency phonons. The observed strong phonon scattering phenomenon in Li-inserted SiNWs is similar to the phonon rattling effect. Our study serves as an exploration of thermal properties of SiNWs as Li battery anodes or weakly coupled with impurity atoms.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Evaluation and Optimization of Well Completion Options for the Utah FORGE Site

Orientation and completion for well pairs that have been subjected to multi-zonal stimulation play a critical role in the long-term performance of an Enhanced Geothermal Reservoir. Enhanced geothermal systems often rely on preferential flow along fractures between well injection and production locations. Modeling this preferential flow using discrete fracture networks (DNF) relies on stochastic realizations of the DFN based on geological sampling. Here we present the development of a stochastic optimization methodology to determine well completion options in a discrete fracture network based on using parallel subset simulation. Stochastic optimization will provide insight into regions where placements of the injection and production wells are optimal. An example optimization of well-pair location optimization based on a deterministic-stochastic DFN model representing FORGE follows a discussion of the theory.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Machine learning for graph-based representations of three-dimensional discrete fracture networks

Structural and topological information play a key role in modeling flow and transport through fractured rock in the subsurface. Discrete fracture network (DFN) computational suites such as dfnWorks (Hyman et al. Comput. Geosci. 84, 10–19 2015) are designed to simulate flow and transport in such porous media. Flow and transport calculations reveal that a small backbone of fractures exists, where most flow and transport occurs. Restricting the flowing fracture network to this backbone provides a significant reduction in the network’s effective size. However, the particle-tracking simulations needed to determine this reduction are computationally intensive. Such methods may be impractical for large systems or for robust uncertainty quantification of fracture networks, where thousands of forward simulations are needed to bound system behavior. In this paper, we develop an alternative network reduction approach to characterizing transport in DFNs, by combining graph theoretical and machine learning methods. We consider a graph representation where nodes signify fractures and edges denote their intersections. Using random forest and support vector machines, we rapidly identify a subnetwork that captures the flow patterns of the full DFN, based primarily on node centrality features in the graph. Our supervised learning techniques train on particle-tracking backbone paths found by dfnWorks, but run in negligible time compared to those simulations. In conclusion, we find that our predictions can reduce the network to approximately 20% of its original size, while still generating breakthrough curves consistent with those of the original network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Mesoscale simulations of High Burnup Fuel Fragmentation: Applying phase-field modeling to understand high burnup fuel fragmentation during loss of coolant accident conditions

During a loss-of-coolant accident (LOCA) in a nuclear power plant, a rapid increase in temperature is experienced in the fuel. In commercial light water reactors (LWRs), the high burnup structure (HBS) forms in regions of the UO 2 fuel where local burnup is high, characterized by smaller grain sizes and large, overpressurized bubbles. During a LOCA transient, fuel in the HBS region is susceptible to fine fragmentation, where the fuel breaks up into micron-size fragments (sometimes referred to as pulverization). To better understand the mechanisms behind this phenomenon, mesoscale computational modeling using the Idaho National Laboratory code Marmot has been employed to simulate the process of HBS formation and response to a LOCA transient. The process of complete and partial HBS formation was demonstrated in prototypical LWR conditions. The simulated HBS microstructures from these simulations were passed to a phase-field model of fracture. To provide a more direct comparison with experiments, a set of experiments conducted by Studsvik on fine fragmentation/pulverization during LOCA conditions was simulated using the BISON fuel performance code. The temperature profiles from the BISON simulations were passed to a phase-field model to determine the gas bubble pressure as a function of time. These pressure histories were used in a phase-field model of fracture to determine under what conditions and in which types of bubbles fine fragmentation was likely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2018

Simulating the Effects of Carbon Nanotube Continuity and Interfacial Bonding on Composite Strength and Stiffness

Molecular dynamics simulations of carbon nanotube (CNT) composites, in which the CNTs are continuous across the periodic boundary, overestimate the experimentally measured mechanical properties of CNT composites along the fiber direction. Since the CNTs in these composites are much shorter than the composite dimensions, load must be transferred either directly between CNTs or through the matrix, a mechanism that is absent in simulations of effectively continuous CNTs. In this study, the elastic and fracture properties of high volume fraction discontinuous carbon nanotube/amorphous carbon composite systems were compared to those of otherwise equivalent continuous CNT composites using ReaxFF reactive molecular dynamics simulations. These simulations were used to show how the number of nanotube-matrix interfacial covalent bonds affect composite mechanical properties. Furthermore, the mechanical impact of interfacial bonding was decomposed to reveal its effect on the properties of the CNTs, the interfacial layer of matrix, and the bulk matrix. For the composites with continuous reinforcement, it was found that any degree of interfacial bonding has a negative impact on axial tensile strength and stiffness. This is due to disruption of the structure of the CNTs and interfacial matrix layer by the interfacial bonds. For the discontinuous composites, the modulus was maximized between 4%-7% interfacial bonding and the strength continues to increase up to the highest levels of interfacial bonding studied. Areas of low stress and voids were observed in the simulated discontinuous composites at the ends of the tubes, from which fracture was observed to initiate. Experimental carbon nanotube yarn composites were fabricated and tested. The results were used to illustrate knockdown factors relative to the mechanical performance of the tubes themselves.

Benjamin D Jensen↗
DOE OSTI2026

Emergence of Diverse Failure Patterns in Weathering‐Induced Landslides: Insights From Particle Finite Element Simulations

Weathering is a fundamental driver of landslide evolution over geological timescales. Despite its ubiquity and importance, quantifying how weathering drives the progressive destabilization of rock slopes remains challenging. In this work, we develop a unified computational framework based on the particle finite element method to investigate the evolution of weathering‐induced landslides, from long‐term weathering to short‐term slope failure and runout dynamics. The framework integrates key processes, including weathering front propagation, time‐dependent strength degradation, rupture surface development, and post‐failure runout dynamics. Through numerical simulation experiments, we elucidate how interactions among weathering characteristics (type, intensity, and rate law), bedrock strength, fracture distribution, and slope geometry govern the failure modes and kinematics of weathering‐induced landslides. Simulations show that matrix‐dominated weathering leads to shallow translational failures, whereas fracture‐dominated weathering produces deep‐seated rotational and compound landslides. Pre‐existing fractures and slope morphology also strongly influence the movement of destabilized landmasses, affecting the failure pattern (e.g., kinematic mode and rupture surface geometry) and post‐failure behavior (e.g., runout velocity). We further demonstrate that the failure time and volume of weathered slopes are governed by the competition between gravitational driving forces and cohesive resisting forces during progressive destabilization. These findings provide new insights into the fundamental mechanisms that drive the emergence of diverse failure patterns of weathering‐induced landslides with important implications for landslide hazard assessment.

Wang, Liang [Eidgenoessische Technische Hochschule↗
DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗
DOE OSTI2023

Natural Fracture Pattern Development Using 3D Coupled Numerical Model: Mechanics and Diagenesis

In the earth's crust, while tectonic processes can mechanically create natural fractures, cementation can simultaneously grow on the surfaces of the fractures. This crystallization can affect the growth of the fractures. To more accurately address the modeling of natural fracture development, both mechanical deformation and chemical diagenesis should be carried out in an integrated manner. Here, we achieve this with a three-dimensional (3D) simulator that couples mechanical and chemical processes for natural fracture pattern development in layer confined situations. Tectonic loading over geologic time is mimicked by displacement driven boundary conditions in a 3D displacement discontinuity method (DDM) scheme. Quartz diagenesis competes with mechanical fracture opening, resulting in three different cementation conditions - open, bridging or fully infilled – which define fracture element stiffness that is used as feedback within the 3D-DDM code during propagation. An important milestone of this work is that this coupled structural diagenetic model can reproduce power-law aperture population statistics which cannot be generated using mechanics alone. More accurately estimated mechanical fracture apertures and their cementation condition is a critical step toward better permeability estimation and fluid path analysis for subsurface reservoirs.

58 GEOSCIENCES↗
DOE OSTI2016

Impact behaviour of Napier/polyester composites under different energy levels

The effects of different energy levels on the impact behaviour of Napier fibre/polyester reinforced composites were investigated. Napier fibre was extracted using traditional water retting process to be utilized as reinforcing materials in polyester composite laminates. 25% fibre loading composite laminates were prepared and impacted at three different energy levels; 2.5,5 and 7.5 J using an instrumented drop weight impact testing machine (IMATEK IM10). The outcomes show that peak force and contact time increase with increased impact load. The energy absorption was then calculated from the force displacement curve. The results indicated that the energy absorption decreases with increasing energy levels of the impact. Impacted specimens were observed visually for fragmentation fracture using an optical camera to identify the failure mechanisms. Fracture fragmentation pattern from permanent dent to perforation with radial and circumferential was observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

CT scanning and flow measurements of shale fractures after multiple shearing events

A shearing apparatus was used in conjunction with a Hassler-style core holder to incrementally shear fractured shale cores while maintaining various confining pressures. Computed tomography scans were performed after each shearing event, and were used to obtain information on evolving fracture geometry. Fracture transmissivity was measured after each shearing event to understand the hydrodynamic response to the evolving fracture structure. The digital fracture volumes were used to perform laminar single phase flow simulations (local cubic law with a tapered plate correction model) to qualitatively examine small scale flow path variations within the altered fractures. Fractures were found to generally increase in aperture after several shear slip events, with corresponding transmissivity increases. Lower confining pressure resulted in a fracture more prone to episodic mechanical failure and sudden changes in transmissivity. Conversely, higher confining pressures resulted in a system where, after an initial setting of the fracture surfaces, changes to the fracture geometry and transmissivity occurred gradually. Flow paths within the fractures are largely controlled by the location and evolution of zero aperture locations. Lastly, a reduction in the number of primary flow pathways through the fracture, and an increase in their width, was observed during all shearing tests.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2019

Flow behavior of hydraulic fractured tight formations considering Pre-Darcy flow using EDFM

For tight reservoirs/formations, flow in fracture and matrix obey different flow patterns. Pre-Darcy flow always exists in low permeability subsurface porous media while Darcy flow can be applied in fracture. It is an outstanding challenge to understand the flow behavior coupling different flow mechanisms. As numerical simulation can predict field-scale performance so the development of efficient simulation method becomes a key issue to evaluate the flow behavior. In this paper, we develop an improved discrete fracture model based on the widely used embedded discrete facture model (EDFM) to simulate the transport both in homogenous and heterogeneous porous media in the presence of Pre-Darcy flow. The conservation equations are derived for fracture and matrix firstly and then discretized by fully implicit method. The transmissibility of different connection types is defined and the detailed calculation workflow is presented as well. The proposed model is verified through the comparison with the model using local grid refined (LGR) method. It is discovered that our simulation method exhibits good accuracy when analyzing the effect of Pre-Darcy flow. Some examples are shown including 1D scalar transport, 2D heterogeneous, 3D flow, complex fracture system and multiple well system problems with fully-coupled flow in fracture and matrix. The pressure distribution and well flux are demonstrated to analyze the Pre-Darcy flow characteristics. Results indicate that the new model works robustly for different flow problems which provide further insight into the Pre-Darcy flow in low permeability media.

02 PETROLEUM↗
DOE OSTI2015

Influences of temperature on permeability changes of flow-paths altered by highly alkaline Ca-rich groundwater - 15192

Fractures (or cracks) in host rock around the geological repository are considered as main flow-paths for the migration of radionuclides. Such fractures would be altered by the alkaline components (Ca, Na, K) leaching from cementitious materials for the construction of the repository. While silicate minerals in fractures are dissolved by highly alkaline groundwater, secondary minerals would deposit on the surface of fractures because of the reaction between the alkaline components and silicate species. The authors' previous study showed the generation of CSH-gel and the clogging of micro flow-path under the condition of Ca-rich and high pH by flow experiments. The present study examined the permeability changes with temperature in micro flow-path flowing highly alkaline Ca-rich solution, considering temperature at the depth of geological disposal system. For the simulation of micro flow-path such as fractures, 'micro flow-cell', which nipped a thin Teflon sheet with a slit between granite chip and Teflon basement, were used for the experiments. Highly alkaline Ca-rich solution (pH>12 and 8.5 mM Ca(OH){sub 2}) were injected into micro flow-cell with 2.0 ml/h of flow rate under the given temperature conditions (278 K, 296 K or 313 K). As a result, the permeability became lower with the increase in temperature. In addition, the deposition rates of CSH-gel became larger with the increase in temperature, and the activation energy was estimated to be 32.0 kJ/mol. These results suggest that flow-paths around the geological repository would undergo the clogging with CSH-gel in a shorter time-period under a condition of higher temperature. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

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↗
DOE OSTI2024

Evaluation of sintered bauxite proppant for binary enhanced geothermal systems

Abstract Solid granular proppant particles are widely used in oil and gas development to sustain permeability through fractures after hydraulic stimulation. Similar proppants are of interest for geothermal applications where the goal of sustaining permeability is the same, but the harsh geothermal environment risks rapid proppant degradation that will reduce fracture permeability. Here, we present proppant conductivity experiments using saw cut granite, tensile fractured granite, and aluminum control specimens packed with sintered bauxite ceramic proppants at concentrations of 0.0, 0.1, and 1.0 kg/m 2 . Simulated geothermal conditions included temperatures up to 130 °C and normal closure stresses up to 60 MPa. Compared to unpropped specimens, peak fracture conductivity enhancement was up to 6 orders of magnitude. At simulated geothermal conditions, proppants were able to sustain fracture conductivity over 60 h, but chemical dissolution and decreasing permeability over time were evident. Irreversible conductivity reductions with crushing and embedment of proppants during loading stages were also observed. Overall, sintered bauxite proppant remains a promising option for low-temperature binary-cycle enhanced geothermal systems.

58 GEOSCIENCES↗
DOE OSTI2015

Wide-range simulation of elastoplastic wave fronts and failure of solids under high-speed loading

The aim of this paper is numerical study of deformation processes and failure of vanadium under shock-wave loading. According developed statistical theory of solid with mesoscopic defects the constitutive equations were proposed in terms of two structural variables characterizing behavior of defects ensembles: defect density tensor and structural scaling parameter. On the basis of wide-range constitutive equations the mathematical model of deformation behavior and failure of vanadium was developed taking into account the bond relaxation mechanisms, multistage of fracture and nonlinearity kinetic of defects. Results of numerical simulation allow the description of the major effects of shock wave propagation (elastic precursor decay, grow of spall strength under grow strain rate)

42 ENGINEERING↗
DOE OSTI2022

Wellbore cement fracture permeability as a function of confining stress and pore pressure

Cemented annulus fractures are a major leakage path in a wellbore system, and their permeability plays an important role in the behavior of fluid flow through a leaky wellbore. The permeability of these fractures is affected by changing conditions including the external stresses acting on the fracture and the fluid pressure within the fracture. Laboratory gas flow experiments were conducted in a triaxial cell to evaluate the permeability of a wellbore cement fracture under a wide range of confining stress and pore pressure conditions. For the first time, an effective stress law that considers the simultaneous effect of confining stress and pore pressure was defined for the wellbore cement fracture permeability. Here the results showed that the effective stress coefficient (λ) for permeability increased linearly with the Terzaghi effective stress ( -p) with an average of λ = 1 in the range of applied pressures. The relationship between the effective stress and fracture permeability was examined using two physical-based models widely used for rock fractures. The results from the experimental work were incorporated into numerical simulations to estimate the impact of effective stress on the interpreted hydraulic aperture and leakage behavior through a fractured annular cement. Accounting for effective stress-dependent permeability through the wellbore length significantly increased the leakage rate at the wellhead compared with the assumption of a constant cemented annulus permeability.

58 GEOSCIENCES↗
DOE OSTI2014

Investigation of Possible Wellbore Cement Failures During Hydraulic Fracturing Operations

We model and assess the possibility of shear failure, using the Mohr-Coulomb model ? along the vertical well by employing a rigorous coupled flow-geomechanic analysis. To this end, we vary the values of cohesion between the well casing and the surrounding cement to representing different quality levels of the cementing operation (low cohesion corresponds to low-quality cement and/or incomplete cementing). The simulation results show that there is very little fracturing when the cement is of high quality.. Conversely, incomplete cementing and/or weak cement can causes significant shear failure and the evolution of long fractures/cracks along the vertical well. Specifically, low cohesion between the well and cemented areas can cause significant shear failure along the well, but the same cohesion as the cemented zone does not cause shear failure. When the hydraulic fracturing pressure is high, low cohesion of the cement can causes fast propagation of shear failure and of the resulting fracture/crack, but a high-quality cement with no weak zones exhibits limited shear failure that is concentrated near the bottom of the vertical part of the well. Thus, high-quality cement and complete cementing along the vertical well appears to be the strongest protection against shear failure of the wellbore cement and, consequently, against contamination hazards to drinking water aquifers during hydraulic fracturing operations.

02 PETROLEUM↗