Search NASA⌕ Search

SEARCH · Search NASA

Results for “2d materials”

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.

208 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Wave Propagation in 2-D Granular Matrix and Dust Mitigation of Fabrics for Space Exploration Mission

Wave Propagation study is essential to exploring the soil on Mars or Moon and Dust Mitigation is a necessity in terms of crew's health in exploration missions. The study of Dust Mitigation has a significant impact on the crew s health when astronauts track dust back into their living space after exploration trips. We are trying to use piezoelectric fiber to create waves and vibrations at certain critical frequencies and amplitudes so that we can shake the particles off from the astronaut s fabrics. By shaking off the dust and removing it, the astronauts no longer have to worry about breathing in small and possibly hazardous materials, when they are back in their living quarters. The Wave Propagation in 2-D Granular Matrix studies how the individual particles interact with each other when a pressure wave travels through the matrix. This experiment allows us to understand how wave propagates through soils and other materials. By knowing the details about the interactions of particles when they act as a medium for waves, we can better understand how wave propagates through soils and other materials. With this experiment, we can study how less gravity effects the wave propagation and hence device a way to study soils in space and on Moon or Mars. Some scientists treat the medium that waves travel through as a "black box", they did not pay much attention to how individual particles act as wave travels through them. With this data, I believe that we can use it to model ways to measure the properties of different materials such as density and composition. In order to study how the particles interact with each other, I have continued Juan Agui's experiment of the effects of impacts on a 2-D matrix. By controlling the inputs and measuring the outputs of the system, I will be able to study now the particles in that system interact with each other. I will also try to model this with the software called PFC2D in order to obtain theoretical data to compare with the experiment. PFC2D is a program that allows the user to control the number of particle's characteristic, and the environment of the particle. With this I can run simulations that mimic the impulse test. This software uses a language called FISH, probably created by the creator of the software. This means that in order to model anything, one must use the command terminal instead of GUI's. I will also use this program to simulate the Moon/Mars simulate adhering to the fabric for the Dust Mitigation project. My goals for this summer are just to complete preliminary studies of the feasibility of the Shaking Fabric, learn the PFC-2D program, and to complete building and testing the wave propagation experiment.

Thanh, Phi Hung X.↗
NASA NTRS2002

Dynamics of Sheared Granular Materials

This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Kondic, Lou↗
NASA NTRS2017

Fusion-Enabled Pluto Orbiter and Lander

The Pluto orbiter mission proposed here is credible and exciting. The benefits to this and all outer-planet and interstellar-probe missions are difficult to overstate. The enabling technology, Direct Fusion Drive, is a unique fusion engine concept based on the Princeton Field-Reversed Configuration (PFRC) fusion reactor under development at the Princeton Plasma Physics Laboratory. The truly game-changing levels of thrust and power in a modestly sized package could integrate with our current launch infrastructure while radically expanding the science capability of these missions. During this Phase I effort, we made great strides in modeling the engine efficiency, thrust, and specific impulse and analyzing feasible trajectories. Based on 2D fluid modeling of the fusion reactors outer stratum, its scrape-off-layer (SOL), we estimate achieving 2.5 to 5 N of thrust for each megawatt of fusion power, reaching a specific impulse, Isp, of about 10,000 s. Supporting this model are particle-in-cell calculations of energy transfer from the fusion products to the SOL electrons. Subsequently, this energy is transferred to the ions as they expand through the magnetic nozzle and beyond. Our point solution for the Pluto mission now delivers 1000 kg of payload to Pluto orbit in 3.75 years using 7.5 N constant thrust. This could potentially be achieved with a single 1 MW engine. The departure spiral from Earth orbit and insertion spiral to Pluto orbit require only a small portion of the total delta-V. Departing from low Earth orbit reduces mission cost while increasing available mission mass. The payload includes a lander, which utilizes a standard green propellant engine for the landing sequence. The lander has about 4 square meters of solar panels mounted on a gimbal that allows it to track the orbiter, which beams 30 to 50 kW of power using a 1080 nm laser. Optical communication provides dramatically high data rates back to Earth. Our mass modeling investigations revealed that if current high-temperature superconductors are utilized at liquid nitrogen temperatures, they drive the mass of the engine, partly because of the shielding required to maintain their critical temperature. Second generation materials are thinner but the superconductor is a very thin layer deposited on a substrate with additional layers of metallic classing. Tremendous research is being performed on a variety of these superconducting materials, and new irradiation data is now available. This raises the possibility of operating nearfuture high-temperature superconductors at a moderately low temperature to dramatically reduce the amount of shielding required. At the same time, a first generation space engine may require low-temperature superconductors, which are higher TRL and have been designed for space coils before (AMS-02 experiment for the ISS). We performed detailed analysis of the startup system and thermal conversion system components. The ideal working fluid was determined to be a blend of Helium and Xenon. No significant problems were identified with these subsystems. For the RF system, we conceived of a new, more efficient design using state-of-the-art switch amplifiers, which have the potential for 100% efficiency. This report presents details of our engine and trajectory analyses, mass modeling efforts, and updated vehicle designs.

temperature↗
NASA NTRS2022

Global geologic map of asteroid (101955) Bennu indicates heterogeneous resurfacing in the past 500,000 years

Global geologic maps are useful tools for efficient interpretation of a planetary body, and they provide global context for the diversity and evolution of the surface. We used data acquired by the OSIRIS-REx spacecraft to create the first global geologic map of the near-Earth asteroid (101955) Bennu. As this is the first geologic map of a small, non-spherical, rubble-pile asteroid, we discuss the distinctive mapping challenges and best practices that may be useful for future exploration of similar asteroids, such as those to be visited with the Hera and Janus missions. By mapping on two centimeter-scale global image mosaics (2D projected space) and a centimeter-scale global shape model (3D space), we generated three input maps respectively describing Bennu’s shape features, geologic features, and surface texture. Based on these input maps, we defined two geologic units: the Smooth Unit and the Rugged Unit. The units are differentiated primarily on the basis of surface texture, concentrations of boulders, and the distributions of lineaments, mass movement features, and craters. They are bounded by several scarps. The Rugged Unit contains abundant boulders and signs of recent mass movement. It also has fewer small(<20 m), putatively fresh craters than the Smooth Unit, suggesting that such craters have been erased in the former. Based on these geologic indicators, we infer that the Rugged Unit has the younger surface of the two. Differential crater size-frequency distributions and the distribution of the freshest craters suggest that both unit surfaces formed ~10–65 million years ago, when Bennu was located in the Main Asteroid Belt, and the Smooth Unit has not been significantly resurfaced in the past 2 million years. Meanwhile, the Rugged Unit has experienced resurfacing within the past ~500,000years during Bennu’s lifetime as a near-Earth asteroid. The geologic units are consistent with global diversity in slope, surface roughness, normal albedo, and thermal emission spectral characteristics. The site on Bennu where the OSIRIS-REx mission collected a regolith sample is located in the Smooth Unit, in a small crater nested within a larger one. So although the Smooth Unit is an older surface than the Rugged Unit, the impact-crater settingindicates that the material sampled was recently exposed. Several similarities are apparent between Bennu and asteroid (162173) Ryugu from a global geologic perspective, including two geologic units distinguishable by variations in the number density of boulders, as well as in other datasets such as brightness.

E R Jawin↗
NASA NTRSDate not supplied

Pink Spinel in Apollo Impact Melt Rock 68815: Implications for Mg-Suite Magmatism

Introduction: Magnesian rocks from the lunar highlands are collectively termed the Mg-suite. Characterized by high (>60) Mg# (molar 100×Mg/[Mg+Fe]) in mafic phases and calcic plagioclase, these rocks are plutonic to hypabyssal in origin, and include a range of bulk mineralogies such as troctolites, dunites, norites, gabbronorites, and spinel troctolites [1]. These Mg-suite lithologies have distinct trace element concentrations and ratios that differentiate them from other lunar rock types. These rocks are ancient, generally dated to between 4.5 and 4.1 Ga, although it is unknown if this represents the full range of Mg-suite ages [1–2]. Among the Mg-suite lithologies, the spinel troctolites are relatively rare, to date only found in polymict breccias [3]. Spinel troctolites, as their name suggests, consist of calcic plagioclase and forsteritic olivine, with minor amounts of spinel (MgAl2O4), ± pyroxene and cordierite [1,4]. This form of spinel is often called ‘pink’ spinel because of its appearance in thin section under plane polarized light (PPL; Fig. 1), due to minor amounts of Cr. Spinel troctolites are generally plutonic or hypabyssal in origin (subsequently exhumed and incorporated into polymict breccias), or formed through impact processes (e.g., crystalline impact melt) [5]. A spinel-rich lithology has also been found in the Moscoviense region of the Moon via the Moon Mineralogy Mapper (M3) and lacks other mafic phases [6]. Finally, while the Mg-suite sampled thus far consists of plutonic (or hypabyssal) rocks, the question remains if such magmas could have erupted on the surface of the Moon [7]. These magmas have much lower density than mare basalts, but little sample or remote sensing evidence has been found to support the idea that extrusive Mg-suite volcanism occurred [7]. Here, we present a coordinated microanalytical study of spinel-bearing lithic and mineral clasts found in Apollo sample 68815. These data will be used to understand their petrogenesis (magmatic or impact) and modification histories, and to shed light on the existence of volcanic Mg-suite rocks. Sample Description: Apollo sample 68815 is a polymict impact melt breccia containing a variety of lithic and mineral fragments embedded in devitrified impact melt. This sample was chipped off the top of a boulder at Station 8 during the Apollo 16 mission and had an original weight of nearly 1.8 kg. In this study, we investigated two polished thin sections of 68815: 68815,17 and 68815,148, both containing spinel. Methods: The thin sections of 68815 were studied using optical light microscopy (PPL, cross-polarized light, and reflected light) with a Keyence VHX-7100 Digital Microscope. Each section was then X-ray mapped for 13–14 elements using a Cameca SX100 electron probe microanalyzer (EPMA) located in the Kuiper Materials Imaging and Characterization Facility (KMICF) at the University of Arizona. We have obtained geochemical information about the phases (olivine, plagioclase, spinel, pyroxene) in the thin sections also using the EPMA. In addition, we have used ThermoScientific Helios NanoLab 660 Focused-Ion-Beam Scanning-Electron Microscope (FIB-SEM) and a Hitachi S-4800 SEM (both in KMICF) to obtain backscattered electron (BSE) images and energy dispersive Xray spectrometry (EDS) maps of areas of interest. Using a JEOL 7900F SEM at the Astromaterials Research & Exploration Science (ARES) at NASA Johnson Space Center (JSC), we have obtained electron backscatter diffraction (EBSD) maps of the spinel-bearing portions of the thin sections. The EBSD data were collected under beam conditions of 20 kV, and ~90 μA, with step sizes varying from 0.05 to 2 μm. Following EBSD data collection, we processed the data using AZtecCrystal and MTEX, a free MATLAB toolbox. Results: We have found clasts with subophitic textures, that consist of primarily olivine and plagioclase, with minor amounts of pink Mg-Al spinel and pyroxene (Fig. 1). These clasts are up to ~1 mm in length and contain spinels up to 50 μm across. We have additionally identified pink Mg-Al spinels within the impact melt (i.e., not contained in lithic clasts) in both thin sections. In one instance, a single spinel grain is approximately 300 μm across (Fig. 1b, 2). The spinel fragments embedded in impact melt have varying compositions, typically distinct from the compositions of spinels in the lithic clasts. Spinel-Bearing Clasts: Ten lithic clasts with similar textures and mineral compositions were identified between 68815,17 (two clasts) and ,148 (eight clasts). These clasts fall into two groups. The first has skeletal olivine with intergranular plagioclase, with minor amounts of pyroxene and spinel (Fig. 1a, 2c, 2d). The spinel in these clasts are found amid the plagioclase. The second group have an intergranular texture of olivine and plagioclase, again with minor spinel and pyroxene. The second group may contain spinels surrounded by plagioclase, and spinels enclosed in olivine. Spinels located within both clast types range from no apparent Cr zoning, to reverse zoning (Cr-enrichment inward; Fig. 2d), to normal zoning (Cr-enrichment outward). In the clasts thus far investigated with EPMA, plagioclase compositions range from An# (molar 100×Ca/[Ca+Na+K]) 92–96. Olivine Mg# ranged from 77 to 94, while pyroxene had Mg# from 54–84. Spinel in the clasts have Cr# (molar 100×Cr/[Cr+Al]) 2–4 and Mg# 88–91, which is within the range of pristine and plutonic spinel troctolites [8]. Isolated Spinels: These crystals are generally euhedral to subhedral, and can exhibit reverse Cr zoning (Cr enrichment inward) or no apparent Cr zoning. The spinels thus far investigated via EPMA have Cr# 9–14 and Mg# 65–82. The Cr# for these spinels is within the range reported by [8], but have lower Mg#. Future Work: We will continue to process the EBSD data for these lithic and mineral clasts. We will also continue to characterize these clasts using EPMA and SEM. By thoroughly characterizing the various spinels and spinel-bearing clasts, we aim to constrain the petrogenesis of these minerals and rock fragments. Acknowledgments: We thank NASA for the loan of these thin sections. We thank Ken Domanik and Jerry Chang for their support with data collection. Work was supported by a University of Arizona RII Core Facilities Pilot Program grant and start-up funds to JJB. We acknowledge support from NASA’s Planetary Science Research program for analysis performed at JSC. References: [1] Shearer C. K. et al. (2015) Am. Min. 100, 294–325. [2] Borg L. E. et al. (2020) GCA 290, 312–332. [3] Warren P. H. (1993) Am. Min. 78, 360–376. [4] Dymek R. F. et al. (1976) LPS VII, 2335–2378. [5] Treiman et al. (2019) Am. Min. 104, 370–384. [6] Pieters et al. (2011) JGR: Plan. 116:E00G08. [7] Prissel et al. (2016) Icarus 277, 319–329. [8] Prissel et al. (2016) Am. Min. 101, 1624–1635.

spinel↗
NASA NTRS2014

AE (Acoustic Emission) for Flip-Chip CGA/FCBGA Defect Detection

C-mode scanning acoustic microscopy (C-SAM) is a nondestructive inspection technique that uses ultrasound to show the internal feature of a specimen. A very high or ultra-high-frequency ultrasound passes through a specimen to produce a visible acoustic microimage (AMI) of its inner features. As ultrasound travels into a specimen, the wave is absorbed, scattered or reflected. The response is highly sensitive to the elastic properties of the materials and is especially sensitive to air gaps. This specific characteristic makes AMI the preferred method for finding "air gaps" such as delamination, cracks, voids, and porosity. C-SAM analysis, which is a type of AMI, was widely used in the past for evaluation of plastic microelectronic circuits, especially for detecting delamination of direct die bonding. With the introduction of the flip-chip die attachment in a package; its use has been expanded to nondestructive characterization of the flip-chip solder bumps and underfill. Figure 1.1 compares visual and C-SAM inspection approaches for defect detection, especially for solder joint interconnections and hidden defects. C-SAM is specifically useful for package features like internal cracks and delamination. C-SAM not only allows for the visualization of the interior features, it has the ability to produce images on layer-by-layer basis. Visual inspection; however, is only superior to C-SAM for the exposed features including solder dewetting, microcracks, and contamination. Ideally, a combination of various inspection techniques - visual, optical and SEM microscopy, C-SAM, and X-ray - need to be performed in order to assure quality at part, package, and system levels. This reports presents evaluations performed on various advanced packages/assemblies, especially the flip-chip die version of ball grid array/column grid array (BGA/CGA) using C-SAM equipment. Both external and internal equipment was used for evaluation. The outside facility provided images of the key features that could be detected using the most advanced C-SAM equipment with a skilled operator. Investigation continued using in-house equipment with its limitations. For comparison, representative X-rays of the assemblies were also gathered to show key defect detection features of these non-destructive techniques. Key images gathered and compared are: Compared the images of 2D X-ray and C-SAM for a plastic LGA assembly showing features that could be detected by either NDE technique. For this specific case, X-ray was a clear winner. Evaluated flip-chip CGA and FCBGA assemblies with and without heat sink by C-SAM. Only the FCCGA package that had no heat sink could be fully analyzed for underfill and bump quality. Cross-sectional microscopy did not revealed peripheral delamination features detected by C-SAM. Analyzed a number of fine pitch PBGA assemblies by C-SAM. Even though the internal features of the package assemblies could be detected, C-SAM was unable to detect solder joint failure at either the package or board level. Twenty times touch ups by solder iron with 700degF tip temperature, each with about 5 second duration, did not induce defects to be detected by C-SAM images. Other techniques need to be considered to induce known defects for characterization. Given NASA's emphasis on the use of microelectronic packages and assemblies and quality assurance on workmanship defect detection, understanding key features of various inspection systems that detect defects in the early stages of package and assembly is critical to developing approaches that will minimize future failures. Additional specific, tailored non-destructive inspection approaches could enable low-risk insertion of these advanced electronic packages having hidden and fine features.

Ghaffarian, Reza↗
NASA NTRS1994

Fabric geometry distortion during composites processing

Waviness and tow misalignment are often cited as possible causes of data scatter and lower compression stiffness and strength in textile composites. Strength differences of as much as 40 percent have been seen in composites that appear to have the same basic material and structural properties -- i.e., yarn orientation, yarn size, interlacing geometry. Fabric geometry distortion has been suggested as a possible reason for this discrepancy, but little quantitative data or substantial evidence exists. The focus of this research is to contribute to the present understanding of the causes and effects of geometric distortion in textile composites. The initial part of the study was an attempt to gather qualitative information on a variety of textile structures. Existing and new samples confirmed that structures with a significant direction presence would be more susceptible to distortion due to the compaction process. Thus, uniweaves (fiber vol frac: 54-72 percent) biaxial braids (vf: 34-58 percent) demonstrated very little fabric geometry distortion. In stitched panels, only slight buckling of z-direction stitches was observed, primarily near the surface. In contrast, for structures with high compaction ratios -- e.g., large cylindrical yarns (2.5:1) orpowder towpreg (4:1) -- there were visible distortions where previously smooth and periodic undulations were transformed to abrupt changes in direction. A controlled study of the effect of forming pressure on distortion was conducted on type 162 glass plain weave fabrics. Panels (6 x 6 in) were produced via a resin infusion type setup, but with an EPON 815 epoxy resin. Pressures ranging from hand layup to 200 psi were used (vf: 34-54 percent). Photomicrographs indicated that at pressures up to 50 psi, large changes in thickness were due primarily to resin squeeze out. At higher pressures, when intimate contact was made between the layers, there was some tow flattening and in-plane shifting to optimize nesting. However, even at 200 psi the period and amplitude of the tow undulation remained constant, suggesting that for this relatively fine fabric, distortions from compaction were not a problem. Because of the interest in using larger tows (to reduce cost) and more complex structures, tests were also run on 2D triaxial glass braid (113 yd/lb at 0, 225 yd/lb at +/- 45). Forming pressures of 20, 50, 200, and 500 psi were used, and short block compression tests were run. The 500 psi specimen had a 10 percent decrease in modulus and an almost 50 percent decrease in strength (vs. 20 psi). Because the total fiber wgt/panel was kept constant, the thickness varied from 0.32 to 0.22 in (49-70 percent vf). Yet, the strength value is clearly below what would be expected, even with the decrease in thickness. Photomicrographs of these samples will be taken to determine if more fabric distortion exists in the 500 psi specimens. Finally, because the ultimate goal is to be able to predict and control distortion in a variety of textile structures, a model compaction test was developed to directly measure the deformation of the tows during compaction. Layers of dry glass fabric were placed in a mold with a clear plexiglass window. The yarn amplitude and period was then calculated using image analysis of the videotaped deformation. Preliminary tests demonstrated the feasibility of this technique for simple fabrics with large tows.

Chen, Julie↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗
NASA NTRS2021

Petrology, Geochemistry, and Pairing of Lunar Meteorites from the Dominion Range

Introduction: During the 2018-2019 Antarctic Search for Meteorites (ANSMET) field season in the Dominion Range (DOM), 7 lunar meteorite stones were collected: DOM 18242 (15.1 g), DOM 18244 (25.1 g), DOM 18262 (6.8 g), DOM 18509 (16.5 g), DOM 18543 (13.6 g), DOM 18666 (45.9 g), and DOM 18678 (11.6 g). Here we present the initial results of electron microprobe and X-ray computed tomography (XCT) studies of these stones and look at the details of their petrography and mineral chemistry, as well as investigate possible pairing relationships, both with each other and with previously described lunar meteorites. Most of the work presented here is on the DOM 18509, 18543, and 18678 stones; subsamples of the other stones are in hand and similar measurements will be made on them by the time of the meeting. Methods: Textures in DOM 18509, DOM 18543, and DOM 18678 were characterized in 2D by optical microscopy, backscattered electron (BSE) imaging and elemental X-ray images on thin sections, as well as in 3D by X-ray computed tomography (XCT) on sample chips. Mineral compositions were assessed through a combination of wavelength dispersive spectroscopy EPMA (electron probe microanalysis) and x-ray mapping on the JEOL 8530 at NASA JSC. The bulk composition of all three meteorites was determined based on analyses of the fusion crust glass. XCT analyses were done on the Nikon XTH 320 at NASA JSC. ICP-MS data on bulk rock subsamples for each meteorite will be carried out in the near future. Results: The stones are all similar in macroscopic appearance with a dark aphanitic matrix hosting a variety of small- to medium-sized angular mineral and lithic fragments (often light colored in nature) [1,2]. Based on EPMA and XCT results, the three meteorites are polymict regolith breccias comprised of mineral, glass, and lithic clasts ranging up to several mm in length. Melt veins run through all three meteorite samples. Mineral clasts in all 3 stones are dominated by pyroxene and plagioclase (An82-96), with minor amounts of SiO2, olivine (Fo1-52), and FeTi-oxides. Pyroxene grains are mostly Fe-rich pigeonite and augite, and larger clasts are normally zoned and have fine exsolution lamellae. The lithic clasts in all stones consists of: (1) basalt clasts that contain zoned pyroxene, plagioclase laths, and ilmenite, with minor silica and Fe-rich olivine; (2) granulitic clasts; (3) anorthosite clasts; (4) Si-rich clasts that also contain ilmenite, troilite, high-Ca pyroxene, fayalite, and K-feldspar likely mesostasis from late stage basalts). All three meteorites contain glassy fusion crust that is highly vesicular, high in FeO and Al2O3 (15-17 wt% each), ferroan (Mg# of 23-24), and moderately rich in TiO2 (1.4-1.8 wt%). The composition of the fusion crust can serve as a proxy for the bulk meteorite composition and is identical within error for all three meteorites. Implications: The lithic and mineral clasts in all three stones are similar in clast population and assemblages as well as mineral chemistry. In addition, the fusion crust composition, a proxy for bulk composition, is within error of each other for all three stones. Thus DOM 18509, DOM 18543, and DOM 18678 are almost certainly paired. Based on similarities in macroscopic description as well as preliminary classification data [1,2], all 7 lunar stones from DOM are likely paired, though additional quantitative analyses are needed to confirm this. The presence of spherules and vesicular fusion crusts indicates that the DOM pairing group is a regolith breccia. The presence of basaltic and gabbroic clasts as well as more feldspathic materials, suggest that the regolith from which these meteorites formed contained a mixture of feldspathic highland material and mare material, suggesting a possible provenance near a marehighlands boundary. No evidence of KREEPy lithologies have been observed so far in these meteorites, however, future ICP-MS data on bulk rock chips for the stones will reveal any KREEP component if present. The DOM pairing group has many similarities to previously described lunar breccia meteorite MET 01210, however more detailed compositional data will be needed to make a definitive comparison.

R.A. Zeigler↗
NASA NTRSDate not supplied

Characterizing Detailed Grain Shape and Size Distribution Properties of Lunar Regolith

Introduction: As the nation prepares to return to the Moon, there is an increasing need for testing tools, instruments, and equipment in simulated environments on Earth to ensure successful operations during lunar missions. Regolith will affect all aspects of future lunar missions, from plume interactions during landing to space suit and tool design [1]. Because of this, it is important to understand the grain shape and size properties of lunar regolith and how those influence regolith behavior in order to prepare for these missions. This knowledge is also vital to create more accurate lunar regolith simulants for testing equipment in a lunar environment. While particle size analyses have been performed on most Apollo soils using simple sieving, shape has only been crudely addressed [2]. This work analyzes 4 lunar regolith samples to provide a better understanding of these size and shape parameters and will provide a more accurate baseline of data to create high-fidelity lunar regolith simulants. Methods: New technologies exist today that are capable of measuring size and shape simultaneously for hundreds of thousands of particles in a single measurement. We conducted a rigorous analysis of the particle size distribution (PSD) as well as the size-dependent 2D and 3D shape parameters of lunar regolith samples of different compositions and maturity levels. This analysis was done using a Microtrac SYNC which provides a unique combination of tri-laser diffraction and Dynamic Image Analysis (DIA). Sample Selection. Four regolith samples were selected for analysis based on maturity level and lunar terrain type: • 10084 – Mature high-Ti mare regolith • 15601 – Immature low-Ti mare regolith • 64501 – Mature highland regolith • 67461 – Immature highland regolith Sample Analysis. After receiving the samples, each sample was imaged with an optical microscope (Figure 1). To obtain 2D and 3D particle size and shape, 0.1 g of each sample was then added to the SYNC for analysis, which outputs more than 30 size and shape parameters for each individual grain as well as the complete size distribution from 0.01-2000 μm by blending laser diffraction and DIA together. For the 0.1 g sample masses requested, ~100,000 grains per sample were captured by DIA. Results: From the PSD analysis, it can be seen that the average particle size of 10084 is ~24.5 µm which is much smaller than the average particle size of the Apollo sample collection (~72 µm) [3], while the average particle size of samples 15601, 64501, and 67461 are larger than the Apollo sample average at ~106 µm, ~103 µm, and ~118.5 µm respectively. The size and shape measurements for the samples output ~30 parameters, four of which were focused on for this study: sphericity, aspect ratio, roundness, and concavity (Figure 3). However, after analyzing the plots it was found that only sphericity and aspect ratio showed differences between the samples. These two parameters are measured on a scale from 0 to 1, with 1 being a perfect sphere with equal dimensions. The results show that the sphericity values are slightly high-er for the mature samples 10084 and 64501 (~0.96) than they are for samples 15601 and 67461 (~0.95) (Figure 4). The aspect ratio values are slightly lower for samples 64501 (~0.7) and 67461 (~0.74) than they are for samples 10084 (~0.76) and 15601 (~0.8) (Figure 5). Discussion: Mature regoliths are those that have been exposed to micrometeorite and solar wind bombardment for long periods of time, breaking up particles and causing them to become more rounded [4]. Therefore, the smaller particle sizes and higher sphericity values for samples 10084 and 64501 are expected due to their higher maturity level compared to 15601 and 67461. However, the aspect ratio values are not dependent on maturity and are instead dependent on terrain type. The lower aspect ratio values for the high-land samples is potentially due to the higher plagioclase content, which occurs in elongated particles and does not break down as easily as the pyroxenes or olivines that are present in the mare. Conclusion and Future Work: The results of this work provide a baseline of high-quality data that will contribute to the creation of high-fidelity lunar simulants and will greatly benefit NASA’s efforts of establishing a human presence on the Moon. Future work includes performing unsupervised image classification on the ~105 particle images per sample in order to identify different classes of grains. These grain classes can then be linked to detailed shape properties, and the relative abundance of each class in the samples can be compared. Acknowledgments: We would like to thank the Extraterrestrial Materials Analysis Group (ExMAG) and the Astromaterials Allocation Review Board (AARB) for allocating lunar regolith samples 10084,9010, 15601,365, 64501,249, and 67461,171 for this work. References: [1] Taylor, L. A. et al. (2005) AIAA #2510. [2] Katagiri et al. (2015) ASCE. [3] Carrier III, W. D. (2005) Lunar Geotech Institute Tech Report. [4] McKay, D. S. et al. (1991) The Lunar Sourcebook, Chapter 7.

S R Deitrick↗