SEARCH · Engineering Papers

Results for “variational fracture modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Micro-tensile characteristics of As-fabricated and irradiated AGR-2 TRISO fuel particle buffer, IPyC, and buffer-IPyC interlayer regions

A recently developed micro-tensile sample preparation technique was implemented to evaluate the tensile strengths of the buffer, IPyC, and buffer-IPyC interlayer regions of the unirradiated and irradiated AGR-2 TRISO fuel particles. Understanding the mechanical properties of the buffer-IPyC interlayer is essential for developing thermomechanical models of buffer-IPyC separation, yet there is a lack of experimental data on its micro-tensile properties. TEM analysis was conducted on these regions to determine the microstructural changes relevant to the samples' tensile properties. In the unirradiated TRISO particle samples, the buffer layer demonstrated the weakest tensile strength, while the IPyC layer exhibited the highest. Conversely, in the irradiated TRISO particle samples, the buffer-IPyC interlayer region showed the lowest tensile strength, with the IPyC layer being the strongest. Fractures in the samples from the buffer-IPyC region predominantly occurred either in the buffer layer or at the buffer-IPyC interface. However, some buffer-IPyC interlayer samples displayed stress-strain and fracture behaviors more akin to the IPyC layer than the buffer layer. Analysis of diffraction patterns suggests that irradiation may have increased anisotropy in the three regions tested. Despite this suggested increase in anisotropy, there was no evidence that it affected the measured strengths. The irradiated TRISO particles demonstrated a considerable increase in void space and a decrease in ultimate tensile strength within the buffer-IPyC interlayer region due to the densification and contraction of the buffer layer. Minor variations in diffraction ring patterns were also observed. These changes, coupled with a significant reduction in the Weibull modulus/shape parameter, imply that irradiation-induced densification leads to tearing between the buffer and IPyC layers at locations of elevated porosity in the buffer-IPyC interlayer region.

Tristructural isotropic (TRISO)↗

NASA NTRS1994

Patterns of brittle deformation under extension on Venus

The development of fractures at regular length scales is a widespread feature of Venusian tectonics. Models of lithospheric deformation under extension based on non-Newtonian viscous flow and brittle-plastic flow develop localized failure at preferred wavelengths that depend on lithospheric thickness and stratification. The characteristic wavelengths seen in rift zones and tessera can therefore provide constraints on crustal and thermal structure. Analytic solutions were obtained for growth rates in infinitesimal perturbations imposed on a one-dimensional, layered rheology. Brittle layers were approximated by perfectly-plastic, uniform strength, overlying ductile layers exhibiting thermally-activated power-law creep. This study investigates the formation of faults under finite amounts of extension, employing a finite-element approach. Our model incorporates non-linear viscous rheology and a Coulomb failure envelope. An initial perturbation in crustal thickness gives rise to necking instabilities. A small amount of velocity weakening serves to localize deformation into planar regions of high strain rate. Such planes are analogous to normal faults seen in terrestrial rift zones. These 'faults' evolve to low angle under finite extension. Fault spacing, orientation and location, and the depth to the brittle-ductile transition, depend in a complex way on lateral variations in crustal thickness. In general, we find that multiple wavelengths of deformation can arise from the interaction of crustal and mantle lithosphere.

Neumann, G. A.↗

DOE OSTI2022

HFTS-1 Natural Joint Data and Engineering Summary

This report provides an analysis of engineering and geologic data collected from the Hydraulic Fracturing Test Site #1 (HTSF-1) project in the southern Midland Basin, Reagan County, Texas. The site is being studied as part of the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative at the National Energy Technology Laboratory (NETL). The data collected is intended to provide a basis of understanding of the site and to construct simulation models. This report provides a summary analysis of the engineering and geologic data collected from the project to provide a basis of understanding of the site and to construct simulation models. The report examines various possible correlations in engineering properties based on natural fracture data from the program, which consisted of 11 horizontal wells, one vertical well, and one slant well. The horizontal wells are in two horizons: 1) Upper and 2) Middle Wolfcamp formations. Program data include: fracture frequency and fracture orientation data from four core runs in a slant well; fracture spacing and orientation from a vertical pilot well; laboratory triaxial testing and mineralogical determinations; and porosity results from magnetic resonance analyses from various wells, together with observations based on the data collection. In addition, available references were reviewed on the site for additional insights. As the focus of the report is on the natural system, hydraulic fracture data from the site were not examined in detail in this report. Data variability is the chief observation in examination of the database. Fracture frequency in the slant well can range from sections with values as high as five fractures per ft to sections up to 100+ ft in length with no natural observed fractures. Fracture spacing across is typically less than 10 ft, but can range up to hundreds of feet. Apparent fracturing shows the trends in two predominate orientations, E-SW and WNW-ESE, but minor variations exist. The rock units vary across the site from siliceous mudstones to calcareous mudstones, showing a general layering with depth. The laboratory properties such as strength and modulus show no apparent trend with depth, but appear to correlate with rock mineralogy with high strength and modulus values where calcium content is high. In addition, an attempt to examine variability and mineralogy on a larger scale was made using a color-coded system based on gamma ray measurements. A staged colored approach was adopted, presuming that lower gamma ray values indicate higher value of calcium content (blue scale) and that higher gamma ray values indicate higher clay mineral content (orange scale). As provided in report appendices, the system correlated well with visual examination of the slant core and the petrofabric analyses of the vertical pilot well. The results showed large variability in mineral content along the horizontal plane across the site. The change in mineral content was also rapid, on a scale less than that of the average hydraulic fracture stage length of about 180 ft.

42 ENGINEERING↗

DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗

DOE OSTI2020

Molecular scale insight of pore morphology relation with mechanical properties of amorphous silica using ReaxFF

Porous materials are typically heterogeneous and they contain large variations of micro-/nano-pore structures, causing complicated behaviors. In continuum models, most mechanical properties of porous materials are estimated based on porosity, while the variations of micro/nano structures are ignored. That could be problematic as the microscopic heterogeneity may affect the mechanical response of porous materials. Thus, understanding micro/nano heterogeneity impact has been the focus in many scientific and engineering subjects. In the present study, we investigated the effect of nanopore structure (including pore shape and orientation) as well as porosity on mechanical properties of amorphous silica (a-SiO 2 ). The pore sizes in our simulations are comparable to the corresponding ones observed in a-SiO 2 based materials. We found that the existing of nanopores strongly influences Young’s modulus (E) and critical energy release rate (G IC ). These properties decrease with increasing porosity. Importantly, the impact of nanopores was characterized by structural parameters of porous materials. In addition to dependency on porosity, Young’s modulus also was found to vary as a function of potential energy per atom, which highly depends on nanopore shape. Furthermore, critical energy release rate was found to increase with increasing ligament length (also known as pore wall thickness). Finally, the results highlighted the importance of nanopore structures, which must be taken into account when studying fracture mechanisms in porous materials. Based on our findings, it was proposed that mechanical properties of porous materials can be controlled by nano-engineering pore structures.

36 MATERIALS SCIENCE↗

NASA NTRS2025

Erosion Behavior of Ti-hBN Multifunctional Coatings in A Custom-Made Planetary Test Rig at Extreme Lunar Temperatures

Spacecraft landings and takeoffs on the lunar surface, along with extreme temperature variations between day and night (-196 to 150° C), cause high-velocity dust impacts and erosion, resulting in the premature failure of structures. Ti/2 vol% hBN coatings were deposited using atmospheric (APS) and vacuum plasma spray (VPS) using cryo-milled powder feedstock to protect the structural components. The erosion performance of coatings at extreme lunar temperature regimes (-150 to 150° C) was evaluated in a custom-made planetary erosion test rig (PETR) at low (50 mph) and high impact velocities (250 mph). The mass loss of VPS coatings was reduced by 50% compared to the APS coatings and 40% compared to the Ti6Al4V substrate. The cryogenic temperature induces brittleness in the material, rendering it susceptible to extreme conditions of material loss. The particle impact-deformation behavior was captured using a high-speed camera to study the erosion mechanism. This analysis revealed chipping in substrates and brittle APS coatings, while particles rebounding and embedding were observed in VPS coatings. Energy calculations, aided by particle trajectory tracking from the high-speed camera, have conclusively shown that VPS coatings absorb 5–10% more energy than APS coatings during erosion tests. A modified erosion index was developed incorporating the fracture toughness and temperatures. New erosion models for brittle and ductile target materials are proposed for developing erosion-resistant material systems.

Abhijith Kunneparambil Sukumaran↗

DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗

NASA NTRS1991

Lake Hoare, Antarctica: sedimentation through a thick perennial ice cover

Lake Hoare in the Dry Valleys of Antarctica is covered with a perennial ice cover more than 3 m thick, yet there is a complex record of sedimentation and of growth of microbial mats on the lake bottom. Rough topography on the ice covering the lake surface traps sand that is transported by the wind. In late summer, vertical conduits form by melting and fracturing, making the ice permeable to both liquid water and gases. Cross-sections of the ice cover show that sand is able to penetrate into and apparently through it by descending through these conduits. This is the primary sedimentation mechanism in the lake. Sediment traps retrieved from the lake bottom indicate that rates of deposition can vary by large amounts over lateral scales as small as 1 m. This conclusion is supported by cores taken in a 3 x 3 grid with a spacing of 1.5 m. Despite the close spacing of the cores, the poor stratigraphic correlation that is observed indicates substantial lateral variability in sedimentation rate. Apparently, sand descends into the lake from discrete, highly localized sources in the ice that may in some cases deposit a large amount of sand into the lake in a very short time. In some locations on the lake bottom, distinctive sand mounds have been formed by this process. They are primary sedimentary structures and appear unique to the perennially ice-covered lacustrine environment. In some locations they are tens of centimetres high and gently rounded with stable slopes; in others they reach approximately 1 m in height and have a conical shape with slopes at angle of repose. A simple formation model suggests that these differences can be explained by local variations in water depth and sedimentation rate. Rapid colonization of fresh sand surfaces by microbial mats composed of cyanobacteria, eukaryotic algae, and heterotrophic bacteria produces a complex intercalation of organic and sandy layers that are a distinctive form of modern stromatolites.

NASA Discipline Exobiology↗

NASA NTRS2015

Survival Times of Meter-Sized Rock Boulders on the Surface of Airless Bodies

This study considers the survival times of meter-sized rock boulders on the surfaces of several airless bodies. As the starting point, we employ estimates of the survival times of such boulders on the surface of the Moon by[1], then discuss the role of destruction due to day-night temperature cycling, consider the meteorite bombardment environment on the considered bodies in terms of projectile flux and velocities and finally estimate the survival times. Survival times of meter-sized rocks on lunar surface: The survival times of hand specimen-sized rocks exposed to the lunar surface environment were estimated based on experiments modeling the destruction of rocks by meteorite impacts, combined with measurements of the lunar surface meteorite flux, (e.g.,[2]). For estimations of the survival times of meter-sized lunar boulders, [1] suggested a different approach based on analysis of the spatial density of boulders on the rims of small lunar craters of known absolute age. It was found that for a few million years, only a small fraction of the boulders ejected by cratering process are destroyed, for several tens of million years approx.50% are destroyed, and for 200-300 Ma, ~90 to 99% are destroyed. Following [2] and other works, [1] considered that the rocks are mostly destroyed by meteorite impacts. Destruction of rocks by thermal-stress. However, high diurnal temperature variations on the surface of the Moon and other airless bodies imply that thermal stresses may also be a cause of surface rock destruction. Delbo et al. [3] interpreted the observed presence of fine debris on the surface of small asteroids as due to thermal surface cycling. They stated that because of the very low gravity on the surface of these bodies, ejecta from meteorite impacts should leave the body, so formation there of fine debris has to be due to thermal cycling. Based on experiments on heating-cooling of cm-scale pieces of ordinary and carbonaceous chondrites and theoretical modeling of expansion of the cracks formed they concluded that thermal fragmentation breaks up rocks larger than a few centimeters more quickly than do micrometeoroid impacts. According to them at 1 AU distance from the Sun the lifetime of 10 cm rock fragments on asteroids with a period of rotation from 2.2 to 6 hours should be only ~103 to 104 years and the larger the rock the faster it gets destroyed. But although [3] are obviously correct stating that impact ejecta should leave small asteroids, the low-velocity part of escaping ejecta will mostly stay in orbits close this given asteroid and part of them will eventually return to it. Moreover, directly beneath the impact point the target rock should be fractured and crushed but may not leave the body (Figure 1). These two points question the conclusions of [3].

Basilevsky, A. T.↗

NASA NTRS1999

Synergistic Effects of Physical Aging and Damage on Long-Term Behavior of Polymer Matrix Composites

The research consisted of two major parts, first modeling and simulation of the combined effects of aging and damage on polymer composites and secondly an experimental phase examining composite response at elevated temperatures, again activating both aging and damage. For the simulation, a damage model for polymeric composite laminates operating at elevated temperatures was developed. Viscoelastic behavior of the material is accounted for via the correspondence principle and a variational approach is adopted to compute the temporal stresses within the laminate. Also, the effect of physical aging on ply level stress and on overall laminate behavior is included. An important feature of the model is that damage evolution predictions for viscoelastic laminates can be made. This allows us to track the mechanical response of the laminate up to large load levels though within the confines of linear viscoelastic constitutive behavior. An experimental investigation of microcracking and physical aging effects in polymer matrix composites was also pursued. The goal of the study was to assess the impact of aging on damage accumulation, in ten-ns of microcracking, and the impact of damage on aging and viscoelastic behavior. The testing was performed both at room and elevated temperatures on [+/- 45/903](sub s) and [02/903](sub s) laminates, both containing a set of 90 deg plies centrally located to facilitate investigation of microcracking. Edge replication and X-ray-radiography were utilized to quantify damage. Sequenced creep tests were performed to characterize viscoelastic and aging parameters. Results indicate that while the aging times studied have limited ]Influence on damage evolution, elevated temperature and viscoelastic effects have a profound effect on the damage mode seen. Some results are counterintuitive, including the lower strain to failure for elevated temperature tests and the catastrophic failure mode observed for the [+/- 45/9O3](sub s), specimens. The fracture toughness for transverse cracks increases with increasing temperature for both systems: transverse cracking was completely absent prior to failure in [+/- 45/903](sub s), and was suppressed for [02/903](sub s). No significant effect of damage on aging or viscoelastic parameters was observed.

Brinson, L. Cate↗

DOE OSTI2021

Simultaneous x-ray and neutron 4D tomographic study of drying-driven hydro-mechanical behavior of cement-based materials at moderate temperatures

The drying of cement-based materials is intimately related to their durability, which has significant economic, social and environmental repercussions. The evolution of the saturation of the pore network and the associated drying shrinkage are in fact leading causes of cracking and of the ensuing penetration of aggressive chemicals. This process is highly heterogeneous, due to the thermo-hydric spatial gradients developing in the material from the exposed surfaces to its core and because of local effects, driven by the intrinsically heterogeneous micro-structure (e.g., by the distribution of pores and aggregates). It follows that macroscopic, sample-scale measurements cannot fully disclose the complexity of the underlying processes. In the last few decades, significant advances in full-field techniques have allowed an unprecedented insight into these local processes. For cement-based materials, x-ray and neutron tomography lend themselves as ideal, and highly complementary, tools for the study of their thermo-hydro-mechanical behavior. Notably, the high sensitivity to density variations of x-ray imaging gives access to the developments of fractures, in 4D (3D + time). On the other hand, neutron tomography allows the study of the evolution of the moisture field in 4D, thanks to its high hydrogen sensitivity. The combination of these two techniques provides a unique insight in thermo-hydro-mechanical couplings, e.g., the effect of cracks on the water content field. This contribution presents novel 5D datasets (3D tomographies along time, plus truly simultaneous x-ray and Neutron rapid acquisitions) in-operando of a cement paste and of a concrete sample heated at moderate temperatures (up to 140{sup ∘}C). The analysis of this 5D data-set (once aligned in time and across modalities) allows for example a quantification of the 4D moisture profiles which were found to predict an overall water loss at hydric equilibrium coherent with the corresponding analytical analysis. In the cement paste sample, the x-ray dataset captures the evolution of an extensive cracking network, opening and propagation toward the core of the sample. A novel analysis procedure is here proposed which allows the extraction of these fractures and the analysis of their interplay with local drying as captured through neutron imaging. This for example reveals the depth of penetration of drying in the vicinity of the fractures along time, which is essential for the assessment and calibration of hydro-mechanical coupled models.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗

NASA NTRS2003

Structural Reliability Using Probability Density Estimation Methods Within NESSUS

A reliability analysis studies a mathematical model of a physical system taking into account uncertainties of design variables and common results are estimations of a response density, which also implies estimations of its parameters. Some common density parameters include the mean value, the standard deviation, and specific percentile(s) of the response, which are measures of central tendency, variation, and probability regions, respectively. Reliability analyses are important since the results can lead to different designs by calculating the probability of observing safe responses in each of the proposed designs. All of this is done at the expense of added computational time as compared to a single deterministic analysis which will result in one value of the response out of many that make up the density of the response. Sampling methods, such as monte carlo (MC) and latin hypercube sampling (LHS), can be used to perform reliability analyses and can compute nonlinear response density parameters even if the response is dependent on many random variables. Hence, both methods are very robust; however, they are computationally expensive to use in the estimation of the response density parameters. Both methods are 2 of 13 stochastic methods that are contained within the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) program. NESSUS is a probabilistic finite element analysis (FEA) program that was developed through funding from NASA Glenn Research Center (GRC). It has the additional capability of being linked to other analysis programs; therefore, probabilistic fluid dynamics, fracture mechanics, and heat transfer are only a few of what is possible with this software. The LHS method is the newest addition to the stochastic methods within NESSUS. Part of this work was to enhance NESSUS with the LHS method. The new LHS module is complete, has been successfully integrated with NESSUS, and been used to study four different test cases that have been proposed by the Society of Automotive Engineers (SAE). The test cases compare different probabilistic methods within NESSUS because it is important that a user can have confidence that estimates of stochastic parameters of a response will be within an acceptable error limit. For each response, the mean, standard deviation, and 0.99 percentile, are repeatedly estimated which allows confidence statements to be made for each parameter estimated, and for each method. Thus, the ability of several stochastic methods to efficiently and accurately estimate density parameters is compared using four valid test cases. While all of the reliability methods used performed quite well, for the new LHS module within NESSUS it was found that it had a lower estimation error than MC when they were used to estimate the mean, standard deviation, and 0.99 percentile of the four different stochastic responses. Also, LHS required a smaller amount of calculations to obtain low error answers with a high amount of confidence than MC. It can therefore be stated that NESSUS is an important reliability tool that has a variety of sound probabilistic methods a user can employ and the newest LHS module is a valuable new enhancement of the program.

Chamis, Chrisos C.↗

NASA NTRS1994

The Large Impact Process Inferred from the Geology of Lunar Multiring Basins

The study of the geology of multiring impact basins on the Moon over the past ten years has given us a rudimentary understanding of how these large structures have formed and evolved on the Moon and other bodies. Two-ring basins on the Moon begin to form at diameters of about 300 km; the transition diameter at which more than two rings appear is uncertain, but it appears to be between 400 and 500 km in diameter. Inner rings tend to be made up of clusters or aligned segments of massifs and are arranged into a crudely concentric pattern; scarp-like elements may or may not be present. Outer rings are much more scarp-like and massifs are rare to absent. Basins display textured deposits, interpreted as ejecta, extending roughly an apparent basin radius exterior to the main topographic rim. Ejecta may have various morphologies, ranging from wormy and hummocky deposits to knobby surfaces; the causes of these variations are not known, but may be related to the energy regime in which the ejecta are deposited. Outside the limits of the textured ejecta are found both fields of satellitic craters (secondaries) and light plains deposits. Impact melt sheets are observed on the floors of relatively unflooded basins. Samples of impact melts from lunar basins have basaltic major-element chemistry, characterized by K, rare-earth elements (REE), P, and other trace elements of varying concentration (KREEP); ages are between 3.8 and 3.9 Ga. These lithologies cannot be produced through the fusion of known pristine (plutonic) rock types, suggesting the occurrence of unknown lithologies within the Moon. These melts were probably generated at middle to lower crustal levels. Ejecta compositions, preservation of pre-basin topography, and deposit morphologies all indicate that the excavation cavity of multiring basins is between about 0.4 and 0.6 times the diameter of the apparent crater diameter. Basin depths of excavation can be inferred from the composition of basin ejecta. A variety of mechanisms has been proposed to account for the formation of basin rings but none of them are entirely plausible. Mechanisms can be divided into two broad groups: (1) forcible uplift due to fluidization of the target; (2) concentric, brittle, fracturing and failure of the target, on regional (megaterraces) to global scales (lithospheric fracturing). Most basin rings are spaced at a constant factor on all planets. Evidence supports divergent ringforming models, so it may be that the ring-locating mechanism differs from the ring-forming mechanism. Thus, large-scale crustal foundering (megaterracing) could occur along concentric zones of weakness created by some type of resonant wave mechanism (fluidization and uplift); such immediate crustal adjustment could then be followed by long-term adjustment of the fractured lithosphere.

Spudis, Paul D.↗

NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗

NASA NTRSDate not supplied

Organic-rich fluids in Plutoâ€™s interiorâ€”A source of cryovolcanic activity

The highly variegated surface of Pluto revealed by the New Horizons spacecraft is composed of both volatile ices (N2, CH4, CO) and H2O ice, all of which are intrinsically colorless. The colors on Pluto reveal a non-ice component generally acknowledged to be a refractory complex organic material (tholins) produced by photolysis in the planet’s atmosphere [1], and by photolysis and radiolysis of the surface ices [2-4]. These processes have produced a range of colors from pale yellow to red to brown, and result in a variation of a factor of 10in albedo across the surface. Tholins are disordered polymer-like materials made of repeating chains of linked subunits and complex combinations of functional groups containing carbon. Nitrogen, oxygen, and other elements can be substituted in both the aliphatic and aromatic subunits. The detection of the spectral absorption of an ammoniated compound associated with surface exposures of red-colored H2O ice (herein, RAW) in the LEISA spectral mapping data from New Horizons suggests that a third source of a colored material is ejected from Pluto’s interior as a water-rich fluid(cryomagma) that rapidly freezes at the low temperature (~40K) of the surface [5,6]. The ammoniated compound is thought to be an ammonia hydrate or an ammoniated salt; the spectral data cover a limited wavelength range, resulting in ambiguity in the identification of the form of the ammonia. The distribution of RAW in the vicinity of tectonically stressed and deformed geological structures, particularly graben (fossae), is seen as the result of both fluid effusion and fountaining ejection of cryoclastic materials in Pluto’s Virgil Fossae region [5,6]. RAW is seen elsewhere on Pluto, both in association with graben complexes and on broad expanses of old surface, but primarily in zones of tectonic stress. The nature of the putative subsurface fluid is unknown, but it must have existed at relatively shallow depths such that the graben faulting could reach it. Geochemical modeling of fluids in the interiors of small planetary bodies in the outer Solar System[7] shows that ammonia-bearing H2O in contact with rock and metal components produce a rich chemical soup, including gases that can, in principle, help propel the fluid from shallow reservoirs through crustal fractures to reach the surface. An organic component of the fluid cryomagma is predicted to originate from thermal processing of the components [ref] and may be an example of one-pot synthesis of such complexes as amino acid precursors [ref].

Dale P Cruikshank↗

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↗

NASA NTRS2005

Mechanical Properties and Microstructure Investigation of Lead Free Solder

While the electronics industry appears to be focusing on Sn-Ag-Cu as the alloy of choice for lead free electronics assembly, ,the exact composition varies by geographic region, supplier and user. Add to that dissolved copper and silver from the printed circuit board traces and surface finish, and there can be significant variation in the final solder joint composition. A systematic study of the mechanical and microstructural properties of Sn-Ag-Cu alloys with Ag varying from 2wt% to 4wt% and Cu varying from 0.5wt% to lSwt%, was undertaken in this research study. Different sample preparation techniques (water quenched, oil quenched and water quenched followed by reflow) were explored and the resulting microstructure compared to that of a typical reflowed lead free chip scale package (CSP) solder joint. Tensile properties (modulus, 0.2% yield strength and the ultimate tensile strength) and creep behavior of selected alloy compositions (Sn-4Ag-1 X u , Sn-4Ag-OSCu, Sn- 2Ag-1 X u , Sn-2Ag-OSCu, Sn-3.5Ag-O.SCu) were determined for three conditions: as- cast; aged for 100 hours at 125OC; and aged for 250 hours at 125OC. There was no significant difference in Young's Modulus as a function of alloy composition. After an initial decrease in modulus after 100 hours at 125"C, there was an insignificant change with further aging. The distribution of 0.2% strain yield stress and ultimate tensile strength as a function of alloy composition was more significant and decreased with aging time and temperature. The microstructures of these alloys were examined using light and scanning electron microscopy (LM and SEM) respectively and SEM based energy dispersive x-ray spectroscopy (EDS). Fracture surface and cross-section analysis were performed on the specimens after creep testing. The creep testing results and the effect of high temperature aging on mechanical properties is presented for the oil quenched samples. In general the microstructure of oil quenched specimen exhibited a eutectic region of Sn with moderately dispersed Ag3Sn intermetallic, surrounded by a dendritic Sn-rich phase. The SEM images of the fracture surface indicated the presence of a tough shear surface at the initial cavity break area and a break line in the middle of specimen along the failure direction. A hyperbolic-sine creep model was adopted and used to fit the creep experiment data. The effect on the mechanical properties by adding the quaternary element bismuth to the Sn-3.5Ag-0.8Cu alloy was measured and compared with the mechanical properties of the ternary alloys. The results of this research study provide necessary data for the modeling of solder joint reliability for a range of Sn-Ag-Cu compositions and a baseline for evaluating the effects of subsequent quaternary additions.

Wang, Qing↗