SEARCH · Engineering Papers
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.
Theoretical principles are outlined for the design of complex cutting inserts in conditions of three-dimensional chip formation. Simulation of the formation and fracture of continuous chip provides the basis for developing a system of complex turning inserts with improved performance. Production experience with the new cutting tool is satisfactory.
Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.
A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.
Coupled solid-fluid numerical simulations of proppant and liquid flowing in realistic rock fracture geometries. The plugging of the flow path due to solid particle interactions is shown to be dependent on fracture properties and solid particle sizes.
Interaction between hydraulic fractures (HF) and natural fractures (NF) may lead to complex fracture networks as a consequence of branching and merging of the natural and hydraulically-induced/reopened fractures. Conventional simulation methods, with assumptions of homogeneous and continuous media, are not able to capture realistic fracture patterns in the presence of pre-existing natural fractures (open or healed). These natural discontinuities alter the local stress field and impact the fracturing fluid pressure distribution when hydraulically-driven fractures approach natural discontinuities. Modeling the propagation of one or more propagating hydraulic fractures in naturally fractured reservoirs should carefully consider the complex mechanical and hydraulic interactions between hydraulically-driven and natural fractures. A rigorous hydraulic fracturing model based on a coupled discrete element method (DEM) and conjugate network flow model is described in this work. The coupled DEM-network flow model conveniently accounts for both continuous poroelasticity and discrete natural fractures simultaneously and provides more realistic fracture growth patterns in low-permeability rocks. This dual lattice DEM-network flow model can capture both tension-induced rock failure and shear-induced natural discontinuity sliding and dilation. Simulation results revealed that reactivation of natural fractures becomes more difficult when the following conditions are present: larger magnitude of stress difference, higher intercepting angle between the natural fracture and the propagating hydraulic fracture (higher implies more orthogonal to the natural fracture), larger natural fracture cohesion and angle of internal friction and/or lower natural fracture permeability - leading to continued growth of the hydraulic fracture. The term stress difference was colloquially used in the preceding sentence. For clarity, we will define this as ΔS = S H,max – S h,min ). The quantitative results of numerous simulations reinforced the intuitive understanding of the underlying phenomena. Other simulations reinforced insights from laboratory experimentation and field measurements that lower injection rates and treating fluid viscosities facilitate preferential fluid intake into natural fractures. Finally, this simulation capability for predicting the nature of complex fracture networks potentially enables optimizing treatment design and completion strategy in naturally fractured reservoirs.
Fracture tests were conducted on 76.2mm wide, 2.3mm thick middle crack tension (M(T)) specimens machined from 2024-T3 aluminum sheet. The specimens were tested on the T-L orientation and comparisons were made to similar tests conducted in the L-T orientation. Measurement of critical crack tip opening angle (CTOA), applied stress, and crack front shape were made as a function of crack extension. A two-dimensional, elastic-plastic finite element analysis was used to simulate the fracture behavior for both orientations. The results indicate that the T-L orientation had a 10 percent lower stress at fracture than similar tests conducted in the L-T orientation. Correspondingly, the critical CTOA in the T-L tests reached a constant value of 4.7 degrees after 2-3mm of crack extension and the L-T tests reached a value of 6 degrees. The fracture surfaces of the T-L specimens were observed to remain flat, while those of the L-T specimens transitioned to a 45 degree slant fracture after about 2-3mm of crack extension. The tunneling behavior of the two orientations also differed; the T-L specimens reached a deeply tunneled stabilized crack front shape while, the L-T specimens were observed to have only a small amount of tunneling once the crack began to grow on the 45 degree slant. The two-dimensional, elastic-plastic finite element analysis was able to simulate the fracture behavior for both the T-L and L-T orientations.
In this study, we examine the use of a linear softening cohesive fracture model (LCFM) to predict single-trace fracture growth in short-rod (SR) and notched 3-point-bend (N3PB) test configurations in Indiana Limestone. The broad goal of this work is to (a) understand the underlying assumptions of LCFM and (b) use experimental similarities and deviations from the LCFM to understand the role of loading paths of tensile fracture propagation. Cohesive fracture models are being applied in prediction of structural and subsurface fracture propagation in geomaterials. They lump the inelastic processes occurring during fracture propagation into a thin zone between elastic subdomains. LCFM assumes that the cohesive zone initially deforms elastically to a maximum tensile stress ( σ max ) and then softens linearly from the crack opening width at σ max to zero stress at a critical crack opening width w 1 . Using commercial finite element software, we developed LCFMs for the SR and N3PB configurations. After fixing σ max with results from cylinder splitting tests and finding an initial Young's modulus ( E ) with unconfined compressive strength tests, we manually calibrate E and w 1 in the SR model against an envelope of experimental data. We apply the calibrated LCFM parameters in the N3PB geometry and compare the model against an envelope of N3PB experiments. For accurate simulation of fracture propagation, simulated off-crack stresses are high enough to require inclusion of damage. Different elastic moduli are needed in tension and compression. We hypothesize that the timing and location of shear versus extensional micromechanical failures control the qualitative macroscopic force-versus-displacement response in different tests. For accurate prediction, the LCFM requires a constant style of failure, which the SR configuration maintains until very late in deformation. The N3PB configuration does not maintain this constancy. To be broadly applicable between geometries and failure styles, the LCFM would require additional physics, possibly including elastoplastic damage in the bulk material and more complicated cohesive softening models.
The emergence of hydrocarbons within shale as a major recoverable resource has sparked interest in fluid transport through these tight mudstones. Recent studies suggest the importance to recovery of microfracture networks that connect localized zones with large organic content to the inorganic matrix. This paper presents a joint modeling and experimental study to examine the onset, formation, and evolution of microfracture networks as shale matures. Both the stress field and fractures are simulated and imaged. A novel laboratory-scale, phase-field fracture propagation model was developed to characterize the material failure mechanisms that play a significant role during the shale maturation process. The numerical model developed consists of coupled solid deformation, pore pressure, and fracture propagation mechanisms. Benchmark tests were conducted to validate model accuracy. Laboratory-grade gelatins with varying Young’s modulus were used as scaled-rock analogs in two dimensional Hele-Shaw cell apparatus. Yeast within the gelatin generates gas in a fashion analogous to hydrocarbon formation as shale matures. These setups allow study and visualization of host rock elastic-brittle fracture and fracture network propagation mechanisms. The experimental setup was fitted to utilize photoelasticity principles coupled with birefringence properties of gelatin to explore visually the stress field of the gelatin as the fracture network developed. Stress optics image analysis and Linear Elastic Fracture Mechanics (LEFM) principles for crack propagation were used to monitor fracture growth for each gelatin type. Observed and simulated responses suggest gas diffusion within and deformation of the gelatin matrix as predominant mechanisms for energy dissipation depending on gelatin strength. LEFM, an experimental estimation of principal stress development with fracture growth, at different stages was determined for each gelatin rheology. Here, the interplay of gas diffusion and material deformation determines the resulting frequency and pattern of fractures. Results correlate with Young’s modulus. Experimental and computed stress fields reveal that fractures resulting from internal gas generation are similar to, but not identical to, type 1 opening mode.
It has been shown experimentally that under mixed tensile and compressive stress states, a corresponding mixed-mode fracture will occur. Here, the formation of mixed-mode fractures is investigated using the combined finite–discrete element method. A series of simulations with confining pressures ranging from 7.5 to 150 MPa generate a spectrum of mixed-mode failure conditions. These stress states at failure span the transition from the tensile to compressive failure. The models reproduce previous experiments on Carrara marble using a dog-bone geometry, and a comparison shows good qualitative agreement with the experimental observations. This paper demonstrates that mixed-mode fracture can be captured via numerical simulations and identifies areas where further research is required to better understand mixed-mode fracture processes.
In this study, all grain boundaries are not equal in their predisposition for fracture due to the complex coupling between lattice geometry, interfacial structure, and mechanical properties. The ability to understand these relationships is crucial to engineer materials resilient to grain boundary fracture. Here, a methodology is presented to isolate the role of grain boundary structure on interfacial fracture properties, such as the tensile strength and work of separation, using atomistic simulations. Instead of constructing sets of grain boundary models within the misorientation/structure space by simply varying the misorientation angle around a fixed misorientation axis, the proposed method creates sets of grain boundary models by means of isocurves associated with important fracture-related properties of the adjoining lattices. Such properties may include anisotropic elastic moduli, the Schmid factor for primary slip, and the propensity for simultaneous slip on multiple slip systems. This approach eliminates the effect of lattice properties from the comparative analysis of interfacial fracture properties and thus enables the identification of structure-property relationships for grain boundaries. As an example, this methodology is implemented to study crack propagation along Ni grain boundaries. Segregated H is used as a means to emphasize differences in the selected grain boundary structures while keeping lattice properties fixed.
Not Available
Throughout U.S. Department of Energy (DOE) complexes, safety engineers employ the five-factor formula to calculate the source term (ST) that includes parameters of airborne release fraction (ARF), respirable fraction (RF) and damage ratio (DR). Limited experimental data on fragmentation of solids, such as ceramic pellets (i.e., PuO2), and container breach due to mechanical insults (i.e., drop and forklift impact), can be supplemented by modeling and simulation using high fidelity computational tools to estimate these parameters. This paper presents the use of Sandia National Laboratories' SIERRA solid mechanics (SM) finite element code to investigate the behavior of the widely utilized waste container (such as 7A Drum) subject to a range of free fall impact and puncture scenarios. The resulting behavior of the container is assessed, and the estimates are presented for bounding DRs from calculated breach areas for the various accident conditions considered. This paper also describes a novel multiscale constitutive model recently implemented in SIERRA/SM that simulates the fracture of brittle materials such as PuO 2 and determines ARF during the fracture process. Furthermore, comparisons are made between model predictions and simple bench-top experiments.
Due to the complexity of welding, not all of the important factors are always properly considered and controlled. An automatic system is required. This report outlines a simulation method and all the important considerations to do this. As in many situations where a defect or failure has occurrred, it is freqently necessary to trouble shoot the system and eventually identify those factors that were neglected. This is expensive and time consuming. Very frequently the causes are materials-related that might have been anticipated. Computer simulation can automatically consider all important variables. The major goal of this presentation is to identify the proper relationship to design, processing and materials variables to welding.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
In this project we are developing a multiscale analysis of the evolution of fracture permeability, using numerical simulations and linear stability analysis. Our simulations include fully three-dimensional simulations of the fracture topography, fluid flow, and reactant transport, two-dimensional simulations based on aperture models, and linear stability analysis.