Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Mathematical Foundations for Nonlocal Interface Problems: Multiscale Simulations of Heterogeneous Materials (Final LDRD Report)

Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗
NASA NTRS1987

Age of fracturing and mesa development in the Elysium Area, northern Martian plains

One of the fundamental questions of Martian crustal history is the origin of the crustal dichotomy between northern plains and southern highlands. Hypotheses for the origin of the dichotomy may be constrained by global scale, geophysical considerations, or by detailed geological studies of the genesis and relative ages of materials and landforms in the northern plains of Mars and along the boundary between the plains and the highlands. This abstract summarizes progress on one aspect of a long-range geological study intended to constrain hypotheses for the dichotomy by tracing the history of the northern plains from the most recent events backward -- essentially the same approach used to understand old events in Earth history. Both the giant impact and the mantle convection models for the crustal dichotomy imply a major coeval fracturing event. As has been known for some time, the present dichotomy boundary lies well south of its original position, at least in some places. Nevertheless, there is severe fracturing of this present dichotomy boundary in many places. The specific objectives of this study are to determine the ages of fracturing along the present dichotomy boundary.

Mcgill, George E.↗
NASA NTRS2010

Experimental Evaluation of Fatigue Crack Initiation from Corroded Hemispherical Notches in Aerospace Structural Materials

A test program was developed and executed to evaluate the influence of corroded hemispherical notches on the fatigue crack initiation and propagation in aluminum 7075-T7351, 4340 steel, and D6AC steel. Surface enhancements such as shot peening and laser shock peening were also incorporated as part of the test effort with the intent of improving fatigue performance. In addition to the testing, fracture mechanics and endurance limit based analysis methods were evaluated to characterize the results with the objective of challenging typical assumptions used in modeling fatigue cracks from corrosion pits. The results specifically demonstrate that the aluminum and steel alloys behave differently with respect to fatigue crack initiation from hemispherical corrosion pits. The aluminum test results were bounded by the fracture mechanics and endurance limit models while exhibiting a general insensitivity to the residual stress field generated by shot peening. The steel specimens were better characterized by the endurance limit fatigue properties and did exhibit sensitivities to residual stresses from the shot peening and laser shock peening

Garcia, Daniel B.↗
DOE OSTI2022

Technical and Economic Evaluation of the First Ever Polymer Flood Field Pilot to Enhance the Recovery of Heavy Oils on Alaska's North Slope via Machine Assisted History Matching

Polymer flooding has become globally established as a potential enhanced oil recovery method for heavy oils. To determine whether this technology may be useful in developing the substantial heavy oil resources on the Alaska North Slope, a polymer flood field pilot commenced at the Milne Point Unit in August 2018. This study seeks to evaluate the results of the field pilot on a technical and economic basis. A reservoir simulation model is constructed and calibrated to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding. Transmissibility contrasts are later reinstated to replicate fracture overextension interpreted in one of the producing wells. The calibrated simulation models produced at each stage of the history matching process are used to forecast oil recovery. These forecasts are used as input for economic analysis, incremental to waterflooding expectations. The simulation forecasts indicate that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone, yielding attractive project economics. However, meaningful variations between simulation scenarios demonstrate that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the technical and economic benefits of an enhanced oil recovery process such as polymer flooding. Subsequently, the simulation model and economic model are used in conjunction to conduct a sensitivity analysis for polymer flood design parameters, from which recommendations are provided for both the continued operation of the current field pilot and future polymer flood designs. The results demonstrate that a higher polymer concentration can be injected due to the development of fractures in the reservoir. The throughput rate should remain high without exceeding operating constraints. A calculated point-forward polymer utilization parameter demonstrates the decreasing efficiency of the polymer flood at later times in the pattern life. Future projects will benefit from starting polymer injection earlier in the pattern life. A pattern with tighter horizontal well spacing will observe a greater incremental benefit from polymer flooding.

Keith, Cody↗
NASA NTRS1993

Strength and deformation properties of basaltic lava flows on planetary surfaces

Basaltic rocks are thought to constitute a volumetrically significant rock type on the Moon, Mercury, Mars, and Venus, in addition to the Earth. Spacecraft images of surfaces with known or suspected basaltic composition on these bodies, particularly on Venus, indicate that these rocks have been deformed in the brittle regime to form faults and perhaps dilatant cracks, in addition to folding and more distributed types of deformation. Predictions of brittle fracture or other types of deformation are made by comparing calculated stresses from a tectonic model to some criterion for rock strength. Common strength criteria used in the planetary science literature for near-surface deformation include a Griffith tensile-strength criterion for intact rock, a Mohr envelope for intact basalt, and a brittle strength envelope based on Byerlee's law of rock frictional resistance. However, planetary terrains of basaltic composition consist of much more than just intact basaltic rock. The aggregate basaltic material, termed the 'rock mass,' consists of both the intact rock and the associated fracture, faults, lithologic contacts, and other discontinuous surfaces. A basaltic rock mass is the relevant material for which strength properties must be defined and calculated model stresses must be compared to in order to more accurately predict brittle deformation. For example, the various strengths of a rock mass are less than that of intact material of the same composition. This means that tectonic models which compare stresses to intact failure strengths overestimate the stresses required for fracture and so underestimate the extent and magnitude of brittle deformation predicted in these models. On the other hand, rock mass shear strength can be greater than that predicted from Byerlee's law. The concept of rock mass strength is central to many engineering design studies in which calculated stresses are used to predict brittle fracture, and this experience indicates that brittle strength envelopes which assume properties for intact rock (Griffith parabolas) or sliding along a single, continuous surface (Byerlee's law) inadequately characterize the tensile, compressive, and shear strengths of rock masses. The criterion adopted here to relate stresses to rock mass fracture is based on a Griffith-type curve for tensile normal stress and a concave downward curve for compressive normal stress. It is the only available criterion that explicitly considers the weakening effects of discontinuities within the rock mass on the stress state required for fracture.

Schultz, Richard A.↗
NASA NTRS1991

Preliminary results on the fracture analysis of multi-site cracking of lap joints in aircraft skins

This paper presents results of a fracture mechanics analysis relevant to fatigue crack growth at rivets in lap joints of aircraft skins are presented. Multi-site damage (MSD) is receiving increased attention within the context of problems of aging aircraft. Fracture analyses previously carried out include small-scale modeling of rivet/skin interactions, larger-scale two-dimensional models of lap joints similar to that developed here, and full scale three-dimensional models of large portions of the aircraft fuselage. Fatigue testing efforts have included flat coupon specimens, two-dimensional lap joint tests, and full scale tests on specimens designed to closely duplicate aircraft sections. Most of this work is documented in the proceedings of previous symposia on the aging aircraft problem. The effect MSD has on the ability of skin stiffeners to arrest the growth of long skin cracks is a particularly important topic that remains to be addressed. One of the most striking features of MSD observed in joints of some test sections and in the joints of some of the older aircraft fuselages is the relative uniformity of the fatigue cracks from rivet to rivet along an extended row of rivets. This regularity suggests that nucleation of the cracks must not be overly difficult. Moreover, it indicates that there is some mechanism which keeps longer cracks from running away from shorter ones, or, equivalently, a mechanism for shorter cracks to catch-up with longer cracks. This basic mechanism has not been identified, and one of the objectives of the work is to see to what extent the mechanism is revealed by a fracture analysis of the MSD cracks. Another related aim is to present accurate stress intensity factor variations with crack length which can be used to estimate fatigue crack growth lifetimes once cracks have been initiated. Results are presented which illustrate the influence of load shedding from rivets with long cracks to neighboring rivets with shorter cracks. Results will also be included for the effect of residual stress due to the riveting process itself. All the results presented in this paper should be regarded as preliminary in the sense that further work by ourselves and others will be required to validate them and to establish their sensitivity to more embellished modelling assumptions.

Aging Aircraft↗
NASA NTRS1984

Fracto-emission - The role of charge separation

Fracto-emission is the emission of particles (e.g., electrons, ions, ground state and excited neutrals, and photons) during and following fracture. It is found that during fracture in vacuum of adhesive bonds and crystalline materials involving large amounts of charge separation on the surface the emission of charged particles, excited neutrals, light, and radio waves occurs with unique and revealing time dependencies. Simultaneous fracto-emission measurements on several systems are reported. The results are interpreted in terms of a conceptual model involving the following steps: (1) charge separation due to fracture, (2) desorption of gases from the material into the crack tip, (3) a gas discharge in the crack, (4) energetic bombardment of the freshly created crack walls, and (5) thermally stimulated electron emission, accompanied by electron stimulated desorption of ions and excited neutrals. In addition to evidence from fracture experiments, results from studies of electron bombardment of a polymer surface are presented.

Dickinson, J. T.↗
NASA NTRS1997

Probabilistic Mesomechanical Fatigue Model

A probabilistic mesomechanical fatigue life model is proposed to link the microstructural material heterogeneities to the statistical scatter in the macrostructural response. The macrostructure is modeled as an ensemble of microelements. Cracks nucleation within the microelements and grow from the microelements to final fracture. Variations of the microelement properties are defined using statistical parameters. A micromechanical slip band decohesion model is used to determine the crack nucleation life and size. A crack tip opening displacement model is used to determine the small crack growth life and size. Paris law is used to determine the long crack growth life. The models are combined in a Monte Carlo simulation to determine the statistical distribution of total fatigue life for the macrostructure. The modeled response is compared to trends in experimental observations from the literature.

Tryon, Robert G.↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2021

Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1973

On the breakup of tectonic plates by polar wandering

The observed boundary system of the major tectonic plates on the surface of the earth lends fresh support to the hypothesis of polar wandering. A dynamic model of the outer shell of the earth under the influence of polar shift is developed. The analysis falls into two parts: (1) deriving equations for stresses caused by polar shifting; and (2) deducing the pattern according to which the fracture of the shell can be expected. For stress analysis, the theory of plates and shells is the dominant feature of this model. In order to determine the fracture pattern, the existence of a mathematical theorem of plasticity is recalled: it says that the plastic flow begins to occur when a function in terms of the differences of the three principal stresses surpasses a certain critical value. By introducing the figures for the geophysical constants, this model generates stresses which could produce an initial break in the lithosphere.

Liu, H. S.↗
NASA NTRS1988

Reliability based analysis of contact problems

The Batdorf model is modified to include the reduction in shear due to the effect of compressive stresses on the crack face. This new formulation was used to obtain the probability of failure of ceramic components under contact stress conditions. The combined effect of the surface and volume flaws are included in the analysis. Due to the nature of the fracture of brittle materials under compressive loading, the component is modeled as a series system in order to establish bounds on the probability of failure.

Powers, Lynn M.↗
NASA NTRSDate not supplied

Experimental Characterization of Non-Associative Plasticity Flow Rule Coefficients and Post-Peak Stress Degradation for the LS-DYNA MAT213 Model

This project is focused on developing an experimental framework for characterizing non-associative plasticity flow rule coefficients through coupon-scale tests for the LS-DYNA MAT213 model. The main objective is to characterize these coefficients based on the multi-scale (i.e., both microscopic and macroscopic) full-field measurement of the evolution of strain and stress fields. This paper focuses on presenting the experimental work on characterizing the full-scale stress-strain curves of T700/LM-PAEK composites under tension, compression, and shear loads. The experimental data set was intended to build a deformation sub-model in the MAT213 model for the material. The strain data were collected using both microscopic and macroscopic digital image correlation techniques. The microscopic technique was particularly useful for fracture cases under small strains. A preliminary simulation result obtained from the MAT213 model is also presented in the paper. The experimental framework herein will be extended to characterize post-peak stress degradation in the composite material and to develop a damage sub-model for the material. This project will contribute to developing a simulation tool based on the MAT213 model for simulating the rate-dependent impact damage in composites under multi-axial loading.

MAT213↗