Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “different dimensions”

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

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

At least 613 records · Page 34

Pair density wave and s ± id superconductivity in a strongly coupled lightly doped Kondo insulator

Here, we investigate the large Kondo coupling limit of the Kondo-Heisenberg model on one- and two-dimensional lattices. Focusing on the possible superconducting states when slightly doping the Kondo insulator state, we identify different pairing modes to be most stable in different parameter regimes. Possibilities include uniform s-wave, pair-density-wave with momentum π (in both one and two dimensions) and uniform s ± id x 2 -y 2 wave (in two dimensions). We attribute these exotic pairing states to the presence of various pair-hopping terms with a "wrong'' sign in the effective model, a mechanism that is likely universal for inducing pairing states with spatially modulated pair wavefunctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Gauge-invariant renormalization of four-quark operators

We study the renormalization of four-quark operators in one-loop perturbation theory. We employ a coordinate-space gauge-invariant renormalization scheme (GIRS), which can be advantageous compared to other schemes, especially in nonperturbative lattice investigations. From our perturbative calculations, we extract the conversion factors between GIRS and the modified minimal subtraction scheme ( MS ¯ ) at the next-to-leading order. As a by-product, we also obtain the relevant anomalous dimensions in the GIRS scheme. A formidable issue in the study of the four-quark operators is that operators with different Dirac matrices mix among themselves upon renormalization. We focus on both parity-conserving and parity-violating four-quark operators, which change flavor numbers by two units ( Δ F = 2 ). The extraction of the conversion factors entails the calculation of two-point Green’s functions involving products of two four-quark operators, as well as three-point Green’s functions with one four-quark and two bilinear operators. The significance of our results lies in their potential to refine our understanding of QCD phenomena, offering insights into the precision of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and shedding light on the nonperturbative treatment of complex mixing patterns associated with four-quark operators. Published by the American Physical Society 2024

Constantinou, M. (ORCID:0000000269881745)↗

Cell dynamics simulations of coupled charge and magnetic phase transformation in correlated oxides

In this work, we present a comprehensive numerical study on the kinetics of phase transition that is characterized by two nonconserved scalar order parameters coupled by a special linear-quadratic interaction. This particular Ginzburg-Landau theory has been proposed to describe the coupled charge and magnetic transition in nickelates and the collinear stripe phase in cuprates. The inhomogeneous state of such systems at low temperatures consists of magnetic domains separated by quasimetallic domain walls where the charge order is reduced. By performing large-scale cell dynamics simulations, we find a two-stage phase-ordering process in which a short period of independent evolution of the two order parameters is followed by a correlated coarsening process. The long-time growth and coarsening of magnetic domains is shown to follow the Allen-Cahn power law. We further show that the nucleation-and-growth dynamics during phase transformation to the ordered states is well described by the Kolmogorov-Johnson-Mehl-Avrami theory in two dimensions. On the other hand, the presence of quasimetallic magnetic domain walls in the ordered states gives rise to a very different kinetics for transformation to the high-temperature paramagnetic phase. In this scenario, the phase transformation is initiated by the decay of magnetic domain walls into two insulator-metal boundaries, which subsequently move away from each other. Implications of our findings to recent nano-imaging experiments on nickelates are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗

Data augmentation for disruption prediction via robust surrogate models

The goal of this work is to generate large statistically representative datasets to train machine learning models for disruption prediction provided by data from few existing discharges. Such a comprehensive training database is important to achieve satisfying and reliable prediction results in artificial neural network classifiers. Here, we aim for a robust augmentation of the training database for multivariate time series data using Student-t process regression. We apply Student-t process regression in a state space formulation via Bayesian filtering to tackle challenges imposed by outliers and noise in the training data set and to reduce the computational complexity. Thus, the method can also be used if the time resolution is high. We use an uncorrelated model for each dimension and impose correlations afterwards via coloring transformations. We demonstrate the efficacy of our approach on plasma diagnostics data of three different disruption classes from the DIII-D tokamak. To evaluate if the distribution of the generated data is similar to the training data, we additionally perform statistical analyses using methods from time series analysis, descriptive statistics, and classic machine learning clustering algorithms.

97 MATHEMATICS AND COMPUTING↗

An adaptive computation mesh for the solution of singular perturbation problems

In singular perturbation problems, control of zone size variation can affect the effort required to obtain accurate, numerical solutions of finite difference equations. The mesh is generated by the solution of potential equations. Numerical results for a singular perturbation problem in two dimensions are presented. The mesh was used in calculations of resistive magnetohydrodynamic flow in two dimensions.

Brackbill, J. U.↗

Tunnel constraint for a jet in crossflow

A facet of a unified tunnel correction scheme which uses wall pressures to determine tunnel induced blockage and upwash is described. With this method, there is usually no need to use data concerning model forces or power settings to find the interference; it follows directly from the pressures and tunnel dimensions. However, highly inclined jets do not produce good pressure signatures and are highly three dimensional, so they must be treated differently. Flow modeling is also discussed.

Wilsden, D. J.↗

Sensitivity analysis for the dynamic response of thermoviscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric, fully-coupled, thermoviscoplastic response of shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of large-strain, transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be viscoplasticity with strain hardening and thermal hardening, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used for the shell equations with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. The energy-balance equation is solved using a Galerkin procedure, with the temperature as the fundamental unknown. Spatial discretization is performed in one dimension (meridional direction) for the momentum and constitutive equations of the shell, and in two dimensions (meridional and thickness directions) for the energy-balance equation. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) for the momentum equation; a predictor-corrector version of the trapezoidal rule is used for the energy-balance equation; and an explicit scheme consistent with the central difference method is used to integrate the constitutive equations. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading.

Kulkarni, Makarand↗

Finite Difference Time Domain Electromagnetic Scattering from Frequency-Dependent Lossy Materials

During this effort the tasks specified in the Statement of Work have been successfully completed. The extension of Finite Difference Time Domain (FDTD) to more complicated materials has been made. A three-dimensional FDTD code capable of modeling interactions with both dispersive dielectric and magnetic materials has been written, validated, and documented. This code is efficient and is capable of modeling interesting targets using a modest computer work station platform. However, in addition to the tasks in the Statement of Work, a significant number of other FDTD extensions and calculations have been made. RCS results for two different plate geometries have been reported. The FDTD method has been extended to computing far zone time domain results in two dimensions. Finally, the capability to model nonlinear materials has been incorporated into FDTD and validated. The FDTD computer codes developed have been supplied, along with documentation, and preprints describing the other FDTD advances have been included with this report as attachments.

Luebbers, Raymond J.↗

Andreyivanovite: A Second New Phosphide from the Kaidun Meteorite

Andreyivanovite (ideally FeCrP) is another new phosphide species from the Kaidun meteorite, which fell in South Yemen in 1980. Kaidun is a unique breccia containing an unprecedented variety of fragments of different chondritic as well as achondritic lithologies. Andreyivanovite was found as individual grains and linear arrays of grains with a maximum dimension of 8 m within two masses of Fe-rich serpentine. In one sample it is associated with Fe-Ni-Cr sulfides and florenskyite (FeTiP). Andreyivanovite is creamy white in reflected light, and its luster is metallic. The average of nine electron microprobe analyses yielded the formula Fe(Cr0.587 Fe0.150 V0.109 Ti0.081 Ni0.060 Co0.002)P. Examination of single grains of andreyivanovite using Laue patterns collected by in-situ synchrotron X-ray diffraction (XRD), and by electron backscattered diffraction revealed that it is isostructural with florenskyite; we were unable to find single crystals of sufficient quality to perform a complete structure analysis. Andreyivanovite crystallizes in the space group Pnma, and has the anti-PbCl2 structure. Previously-determined cell constants of synthetic material [a = 5.833(1), b = 3.569(1), c = 6.658(1) A] were consistent with our XRD work. We used the XPOW program to calculate a powder XRD pattern; the 5 most intense reflections are d = 2.247 (I = 100), 2.074 (81), 2.258 (46), 1.785 (43), and 1.885 A (34). Andreyivanovite is the second new phosphide to be described from the Kaidun meteorite. Andreyivanovite could have formed as a result of cooling and crystallization of a melted precursor consisting mainly of Fe-Ni metal enriched in P, Ti, and Cr. Serpentine associated with andreyivanovite would then have formed during aqueous alteration on the parent asteroid. It is also possible that the andreyivanovite could have formed during aqueous alteration, however, artificial FeTiP has been synthesized only during melting experiments, at low oxygen fugacity, and there is no evidence that a hydrothermal genesis is reasonable.

Zolensky, Michael↗

Fabrication and Modification of Nanoporous Silicon Particles

Silicon-based nanoporous particles as biodegradable drug carriers are advantageous in permeation, controlled release, and targeting. The use of biodegradable nanoporous silicon and silicon dioxide, with proper surface treatments, allows sustained drug release within the target site over a period of days, or even weeks, due to selective surface coating. A variety of surface treatment protocols are available for silicon-based particles to be stabilized, functionalized, or modified as required. Coated polyethylene glycol (PEG) chains showed the effective depression of both plasma protein adsorption and cell attachment to the modified surfaces, as well as the advantage of long circulating. Porous silicon particles are micromachined by lithography. Compared to the synthesis route of the nanomaterials, the advantages include: (1) the capability to make different shapes, not only spherical particles but also square, rectangular, or ellipse cross sections, etc.; (2) the capability for very precise dimension control; (3) the capacity for porosity and pore profile control; and (4) allowance of complex surface modification. The particle patterns as small as 60 nm can be fabricated using the state-of-the-art photolithography. The pores in silicon can be fabricated by exposing the silicon in an HF/ethanol solution and then subjecting the pores to an electrical current. The size and shape of the pores inside silicon can be adjusted by the doping of the silicon, electrical current application, the composition of the electrolyte solution, and etching time. The surface of the silicon particles can be modified by many means to provide targeted delivery and on-site permanence for extended release. Multiple active agents can be co-loaded into the particles. Because the surface modification of particles can be done on wafers before the mechanical release, asymmetrical surface modification is feasible. Starting from silicon wafers, a treatment, such as KOH dipping or reactive ion etching (RIE), may be applied to make the surface rough. This helps remove the nucleation layer. A protective layer is then deposited on the wafer. The protective layer, such as silicon nitride film or photoresist film, protects the wafer from electrochemical etching in an HF-based solution. A lithography technique is applied to pattern the particles onto the protective film. The undesired area of the protective film is removed, and the protective film on the back side of the wafer is also removed. Then the pattern is exposed to HF/surfactant solution, and a larger DC electrical current is applied to the wafers for a selected time. This step removes the nucleation layer. Then a DC current is applied to generate the nanopores. Next, a large electrical current is applied to generate a release layer. The particles are mechanically suspended in the solvent and collected by filtration or centrifuge.

Ferrari, Mauro↗

Stochastic parametric skeletal dosimetry model for humans: Anatomical-morphological basis and parameter evaluation

Radiation exposure of the hematopoietic system that results in a radiation dose to bone marrow of more than 100 mGy leads to an increase in the risk of leukemia in humans. Excess relative risk of leukemia was observed in cohorts whose members lived in the territories of the Southern Urals that were radioactively contaminated in the 1950s. As part of the dosimetric support of epidemiological studies of these cohorts, an original methodology for stochastic bone dosimetric modeling was developed, termed the Stochastic Parametric Skeletal Dosimetry (SPSD) model. The purpose of this work was to present the anatomical and morphological bases of the SPSD model, which includes an assessment of the parameters of the microstructure of the trabecular bone in the hematopoietic areas of the human skeleton, as well as a description of the macrostructural division (segmentation) of hematopoietic areas into simple bone segments. As a result, an anatomical-morphological basis of the SPSD model was created based on published data. Data collection work included the analysis of original articles, atlases, manuals, monographs and the formation of primary data files. Data on the duration of hematopoiesis in various parts of the skeleton; and data on age-related changes in the microstructure and linear dimensions of human bones and their segments were analyzed. The paper describes the full set of parameters to be used for the dosimetric model for newborns, children aged 1, 5 and 10 years, as well as for adolescents aged 15 years and adults; for the latter, sex differences in bone size were considered. In total, the SPSD model includes 289 unique basic (bone) phantom segments, each of which is described by 7 or more parameters describing the microstructure, thickness of the cortical layer and linear dimensions. Population variability was estimated for each parameter. The approach to SPSD modeling, i.e., the use of simple geometric shapes, was successfully verified using independent datasets on bone masses and volumes.

Science & Technology - Other Topics↗

A regional comparison of sub-daily flow variability in regulated and unregulated rivers in the United States

Regulating rivers for hydropower or other purposes can dramatically alter river flow patterns, including creating substantial changes in flow over short, minutes-to-hours-long timespans known as sub-daily flow variability (SDFV). The impacts of flexible hydropower production on flow and aquatic organisms are increasingly documented in research. However, the degree to which flow alteration relates to different hydropower operational modes in distinct geographical regions and seasons is not well understood. This study offers a methodology for regional- and species-appropriate evaluations of potential impacts of flow on fish based on sub-daily flow characteristics of hydropower operational modes. We analyzed 15-min discharge data between 2018 and 2021 from 69 USGS stream gages to compare SDFV in hydropeaking, run-of-river, and unregulated systems in the US Southeast and Pacific Northwest. Regulated systems exhibited significant SDFV downstream from hydropower facilities relative to unregulated systems, but specific impacts differed between regions. Regulated systems in the Southeast were characterized by high flow coefficients of variation and ratios (hydropeaking only) and extended durations of daily upramping flow phases. Regulated systems in the Pacific Northwest were characterized by many short flow phases per day and large portions of the day spent upramping. Pacific Northwest unregulated systems displayed the strongest seasonal flow patterns while Southeastern hydropeaking systems displayed the greatest SDFV. Given that SDFV impacts multiple dimensions of fish ecology, region-specific sub-daily flow signatures have important implications for understanding and mitigating potential community-, species-, and age-specific effects on fish in different parts of the country.

Fish↗

Nonlinear Wakes Behind a Row of Elongated Roughness Elements

This paper is concerned with the high Reynolds number flow over a spanwise periodic array of roughness elements with inter-element spacing of the order of the local boundary-layer thickness. While earlier work by Goldstein, Sescu, Duck and Choudhari (2010) and Goldstein, Sescu, Duck and Choudhari (2011) was mainly concerned with smaller roughness heights that produced relatively weak distortions of the downstream flow, the focus here is on extending the analysis to larger roughness heights and streamwise elongated planform shapes that together produce a qualitatively different, nonlinear behavior of the downstream wakes. The roughness scale flow now has a novel triple-deck structure that is somewhat different from related studies that have previously appeared in the literature. The resulting flow is formally nonlinear in the intermediate wake region, where the streamwise distance is large compared to the roughness dimensions but small compared to the downstream distance from the leading edge, as well as in the far wake region where the streamwise length scale is of the order of the downstream distance from the leading edge. In contrast, the flow perturbations in both of these wake regions were strictly linear in the earlier work by Goldstein et al (2010, 2011). This is an important difference because the nonlinear wake flow in the present case provides an appropriate basic state for studying the secondary instability and eventual breakdown into turbulence.

Boundary layer control↗

Automatic Segmentation of Building Envelope Point Cloud Data Using Machine Learning

About 50% of buildings in the US were constructed before energy codes were introduced. Modular overclad panel retrofits, in which a new envelope is constructed over the existing building, are a promising solution given that it minimizes occupant disruption and shortens construction time at the jobsite. Current state-of-the-art retrofit panel layout and dimensioning consists of three steps: 1) 3D point cloud data generation of the building envelope using commonly available surveying equipment, 2) manual segmentation of 3D point cloud data by a trained professional to identify and dimension window openings, door openings, and other architectural features, and 3) modular panel layout optimization and dimensioning by an architect or engineer. Among these steps, the second one remains the most difficult and costly because it is very labor-intensive. We propose a methodology to automatically label 3D point cloud data to reduce the time and expense spent in manual segmentation. Machine learning methods were employed to classify the point cloud data into distinct groups, each of which corresponds to different features of the building envelope. After classification, a segmentation algorithm was developed to perform boundary detection and separate the components of the façade. Finally, the algorithm returns the relative positions and dimensions of the features in the building envelope. The measurements obtained with the proposed automated method were compared against the actual dimensions to determine the overall algorithm accuracy. The proposed algorithm can then be used to reduce manual efforts for 3D point cloud labeling before modular panel layout optimization is performed.

Maldonado Puente, Bryan↗

Investigating Mechanisms Driving Differences in the Characteristics of Precipitation in the E3SM Multiscale Modeling Framework With 2D Versus 3D Cloud Resolving Model Configurations

Abstract In this study, we compare the Energy Exascale Earth Systems Model (E3SM) multiscale modeling framework (MMF) with the cloud resolving model (CRM) configured in two (2dMMF) and three (3dMMF) dimensions. We explore how CRM dimensionality impacts the representation of mean and extreme precipitation characteristics. Our results show that tropical mean precipitation patterns are better represented in 3dMMF compared to observations (Integrated Multi‐satellitE Retrivals for GPM and Global Precipitation Climatology Project One Degree Daily products), while 2dMMF better captures extreme precipitation intensity, with systematic land‐ocean differences in precipitation and cloud‐associated variables. These differences are attributed to the co‐occurrence of CRM throttling (i.e., suppressed convection in due to smaller numbers of CRM columns and domain size) and dilution (i.e., 3‐D cloud circulations with increased entrainment and lower precipitation efficiency) effects. Overall, throttling results in more low‐level humidity in 2dMMF and dilution contributes to more high clouds with less precipitation efficiency in 3dMMF. Since throttling occurs more strongly over the ocean than land, the 3dMMF tends to have less cloud liquid and precipitation over the ocean and more cloud ice and precipitation over land. These results may serve as a guide for choosing the CRM structure to reduce precipitation and cloud‐related biases.

54 ENVIRONMENTAL SCIENCES↗

Numerical Simulations of Single and Multiple Scattering by Fractal Ice Clusters

We consider the scattering model in the form of a vertically and horizontally homogeneous particulate slab of an arbitrary optical thickness composed of widely separated fractal aggregates built of small spherical ice monomers. The aggregates are generated by applying three different approaches, including simulated cluster-cluster aggregation (CCA) and diffusion-limited aggregation (DLA) procedures. Having in mind radar remote-sensing applications, we report and analyze the results of computations of the backscattering circular polarization ratio obtained using efficient superposition T-matrix and vector radiative-transfer codes. The computations have been performed at a wavelength of 12.6 cm for fractal aggregates with the following characteristics: monomer refractive index m=1.78+i0.003, monomer radius r=1 cm, monomer packing density p=0.2, overall aggregate radii R in the range 4<=R<=10 cm and fractal dimensions D(sub f) 2.5 and 3. We show that for aggregates generated with simulated CCA and DLA procedures, the respective values of the backscattering circular polarization ratio differ weakly for D(sub f) 2.5, but the differences can increase somewhat for D(sub f)3, especially in case of an optically semi-infinite medium. For aggregates with a spheroidal overall shape, the dependence of the circular polarization ratio on the cluster morphology can be quite significant and increases with increasing the aspect ratio of the circumscribing spheroid.

Dlugach, Janna M.↗

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↗