Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coarse graining”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microscopic Theory of Fluctuating Hydrodynamics in Nonlinear Lattices

The theory of fluctuating hydrodynamics has been an important tool for analyzing macroscopic behavior in nonlinear lattices. However, despite its practical success, its microscopic derivation is still incomplete. In this work, we provide the microscopic derivation of fluctuating hydrodynamics, using the coarse-graining and projection technique; the equivalence of ensembles turns out to be critical. The Green-Kubo (GK)-like formula for the bare transport coefficients are presented in a numerically computable form. Our numerical simulations show that the bare transport coefficients exist for a sufficiently large but finite coarse-graining length in the infinite lattice within the framework of the GK-like formula. Finally, this demonstrates that the bare transport coefficients uniquely exist for each physical system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Carbonaceous chondrite clasts in the howardites Bholghati and EET87513

Twenty-two carbonaceous chondrite clasts from the two howardites Bholghati and EET87513 were analyzed. Clast N from EET87513 is a fragment classified as CM2 material on the basis of texture, bulk composition, mineralogy, and bulk O isotopic composition. Carbonaceous chondrite clasts from Bholghati, for which less data are available because of their small size, can be divided into two petrologic types: C1 and C2. C1 clasts are composed of opaque matrix with rare coarse-grained silicates as individual mineral fragments; textures resemble CI meteorites and some dark inclusions from CR meteorites. Opaque matrix is predominantly composed of flaky saponite; unlike typical CI and CR meteorites, serpentine is absent in the samples we analyzed. C2 clasts contain chondrules, aggregates, and individual fragments of coarse-grained silicates in an opaque matrix principally composed of saponite and anhydrous ferromagnesian silicates with flaky textures similar to phyllosilicates. These anhydrous ferromagnesian silicates are interpreted as the product of heating of pre-existing serpentine. The carbonaceous chondrite clasts we have studied from these two howardites are, with one notable exception (clast N from EET87513), mineralogically distinct from typical carbonaceous chondrites. However, these clasts have very close affinities to carbonaceous chondrites and have also experienced thermal metamorphism and aqueous alteration, but to different degrees.

Buchanan, P. C.↗
DOE OSTI2024

Beyond nothingness in the formation and functional relevance of voids in polymer films

Abstract Voids—the nothingness—broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

36 MATERIALS SCIENCE↗
NASA NTRS2006

Formation of Metal and Silicate Globules in Gujba: A New Bencubbin-like Meteorite Fall

Gujba is a coarse-grained meteorite fall composed of 41 vol% large kamacite globules, 20 vol% large light-colored silicate globules with cryptocrystalline, barred pyroxene and barred olivine textures, 39 vol% dark-colored, silicate-rich matrix, and rare refractory inclusions. Gujba resembles Bencubbin and Weatherford in texture, oxygen-isotopic composition and in having high bulk delta N-15 values (approximately +685%0). The He-3 cosmic-ray exposure age of Gujba (26 +/- 7 Ma) is essentially identical to that of Bencubbin, suggesting that they were both reduced to meter-size fragments in the same parent-body collision. The Gujba metal globules exhibit metal-troilite quench textures and vary in their abundances of troilite and volatile siderophile elements. We suggest that the metal globules formed as liquid droplets either via condensation in an impact-generated vapor plume or by evaporation of preexisting metal particles in a plume. The lower the abundance of volatile elements in the metal globules, the higher the globule quench temperature. We infer that the large silicate globules also formed from completely molten droplets; their low volatile-element abundances indicate that they also formed at high temperatures, probably by processes analogous to those that formed the metal globules. The coarse-grained Bencubbin-Weatherford-Gujba meteorites may represent a depositional component from the vapor cloud enriched in coarse and dense particles. A second class of Bencubbin-like meteorites (represented by Hammadah a1 Hamra 237 and QUE 94411) may be a finer fraction derived from the same vapor cloud

Rubin, Alan E.↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
DOE OSTI2022

Polymers in Deep Eutectic Solvents

The project investigated the behavior of polymers in ionic liquids such as deep eutectic solvents using an array of techniques. These included the development of atomistic force fields, coarse graining these force fields to the united atom level, large scale molecular dynamics (MD) simulations using new thermostat algorithms, and machine learning (ML) methods for the phase behavior. The project demonstrated the feasibility and accuracy of first principles force fields for ionic liquids, deep-eutectic solvents, and urea-water mixtures. A novel hierarchical coarse graining method was then used to develop accurate and efficient united-atom models, using which microsecond simulations were performed for polymers in ionic liquids. These simulations were in quantitative agreement with experiment, thus resolving previous controversies. Methods were also developed to obtain the potential of mean force between complex ions in solution. Finally, supervised ML methods were developed for the phase behavior of polymers in ionic liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2013

Mineralogical, Chemical, and Isotopic Heterogeneity in Zagami: Evidence for a Complex Petrogenesis

Textural variations in the shergottite Zagami were initially interpreted as evidence that it formed in a heterogeneous lava flow. Variations in initial Sr-87/Sr-86 ratios between a Coarse Grained (CG) and a Fine Grained (FG) lithology and evidence for more extensive fractionation of the Rb/Sr ratio in a Dark Mottled Lithology (DML) are consistent with such an interpretation. More recently, Niihara et al. and Misawa et al. have reported the mineralogy and Sr-isotopic systematics of an Olivine Rich Lithology (ORL) found in association with the coarse-grained DML lithology in the Kanagawa Zagami specimen [6,7]. Here we call this lithology DML(Ka) to maintain a distinction with DML(USNM) as studied. An Ar-Ar study by Park et al. of a late stage K-rich melt enriched in K2O to approx 7% and intruded into ORL yielded an Ar-Ar age of 202+/0 7 Ma. The present work extends the study of Kanagawa Zagami to Nd-isotopes.

Nyquist, L. E.↗
DOE OSTI2023

Classification of four-qubit entangled states via machine learning

We apply the support vector machine (SVM) algorithm to derive a set of entanglement witnesses (EW) to identify entanglement patterns in families of four-qubit states. The effectiveness of SVM for practical EW implementations stems from the coarse-grained description of families of equivalent entangled quantum states. The equivalence criteria in our work is based on the stochastic local operations and classical communication classification and the description of the four-qubit entangled Werner states. We numerically verify that the SVM approach provides an effective tool to address the entanglement witness problem when the coarse-grained description of a given family state is available. Here, we also discuss and demonstrate the efficiency of nonlinear kernel SVM methods as applied to four-qubit entangled state classification.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Trace Elemental Abundances in Calcium-Aluminum-Rich Inclusions in CV Chondrites

Introduction: Calcium-aluminum-rich inclusions (CAIs), are the first formed solids that define the age of the Solar System [1,2]. CAIs are thought to have condensed from nebular gas [3,4] within the first <1 Ma of Solar System formation [5,6]. CAIs have experienced numerous early Solar System processes including condensation, evaporation, melting, recrystallization, and aqueous alteration [e.g., 7]. The chemical, mineralogical, and textural diversity among CAIs results from a range of chemical and physical processes recorded during nebular and parent body epoch. This study aims to explore the mineralogical, textural, and chemical compositions of CAIs including the trace elemental abundances in CAI phases to determine the early Solar System processes recorded in them. Samples and Analytical Methods: We analyzed one CAI each from CV3 chondrites Northwest Africa (NWA) 5508 designated as ‘Saguaro’, and Northwest Africa (NWA) 12772 designated as ‘Hoopoe’. Back-scatter electron (BSE) images were collected using a Phenom XL scanning electron microscope (SEM) at the Lunar and Planetary Institute (LPI) and the JEOL JXA-8530F electron probe microanalyzer (EPMA) at Johnson Space Center (JSC)-NASA. Additionally, energy dispersive X-ray spectrometry (EDS) elemental maps of select areas for these samples were collected using a 15.0kV beam energy and a 40µA emission current. Using the EPMA, wavelength-dispersive X-ray spectroscopy (WDS) quantitative data were collected. In-situ trace element measurements for both CAIs were determined at JSC-NASA using a Photon Machines 193nm laser ablation system and a Thermo-Scientific Element-XR inductively coupled plasma mass spectrometer (ICP-MS). Analyses consisted of 30s ablations at 10Hz, spot sizes of 20-25µm, and a fluence of 6.0 J/cm2 for anorthite and melilite, and a 3.5 J/cm2 fluence for all other phases. NIST612 was used to correct for instrument drift, while BHVO-2g was used as a primary calibration standard. BCR-2g and in-house mineral standards were regularly measured as unknowns to ensure accuracy. Results: Saguaro is a coarse-grained CAI, ~11 x 6 mm in dimensions. Saguaro contains spinel, Al-rich pyroxene, anorthite, Mg-rich melilite, and minor perovskite in its interior and is therefore classified as a Type B CAI. Individual melilite grains shows normal compositional zoning with an Ak content ranging from ~24 to 54 with no apparent trend from the core to the edge of the CAI. The spinel appears euhedral and occurs both as clusters and as spinel palisades [8]. Two rim sequences surround most of the sample: the inner rim being a Wark-Lovering (WL) rim (~10-35 µm) containing pyroxene, spinel, and melilite (or anorthite), and the outer rim is a finer-grained, thicker (~100 µm), accretionary rim (Fig. 1). The mineral phases in Saguaro record an overall flat REE pattern with an average negative Eu anomaly in pyroxene, and an average positive Eu anomaly in anorthite and melilite respectively. Anorthite, melilite, and pyroxene have a minor depletion in Tm (Fig. 2). The Hoopoe CAI is a compact, coarse-grained ~6 × 4 mm in size. The major mineralogy includes hibonite, spinel, melilite, anorthite, and perovskite. Therefore, it is classified as a compact transitional type A and B. (?)zoning was observed in some hibonites. Individual melilite grains show both reverse and normal zoning, where the Ak content ranges from ~6- to 28. Melilite shows two distinct textures. One texture consisted of smooth melilite that appeared homogenous, while the second appeared to consist of many fine fractures. The spinel also often appears clustered. The WL-rim sequence surrounding Hoopoe is ~25 µm thick and composed of spinel, perovskite, hibonite, and melilite/anorthite. It is then partially surrounded by an outer accretionary rim (~75µm). Like before, refractory metal nuggets appeared concentrated near the WL rims. Other metal assemblages rich in Fe and Ni were also observed. All major mineral phases in Hoopoe display relatively flat REE patterns, except for varying Eu and Tm between phases (Fig. 2). There is a prominent negative Eu anomaly in perovskite and an average positive Eu anomaly in spinel, anorthite, and hibonite respectively (Fig. 2, 3). The mixed phases along the rim of the CAI also display a negative Eu anomaly, and all phases the CAI were depleted in Pb. Discussion: The CV3 CAIs analyzed in this study were classified based on their mineralogy and textures into Type A versus Type B CAIs [10]. Hibonite appears to be pseudomorphically replacing the spinel, (i.e., is hibonite in composition, but appears in the shape of spinel). Spinel palisades. The presence of spinel palisades present in Saguaro are consistent with the melting and recrystallization experienced by this CAI. Trace elemental analyses. Saguaro and Hoopoe display similar trace element patterns to each other, with both appearing generally flat, with anomalies in Eu, and Tm. Melilite and anorthite display positive Eu anomalies in both CAIs, in addition to the hibonite in Hoopoe (Fig. 2). The phases that are depleted in Eu are pyroxene and perovskite in both Saguaro and Hoopoe, respectively (Fig. 3). Given that Eu is volatile in reducing environments [11], this could possibly indicate reducing conditions at the time anorthite and melilite crystallized, with the gas they formed from containing Eu. As these CAIs continued to form, this gas as a result would become depleted in Eu. This also could be supported by the propensity of anorthite and melilite to take up Eu from its surroundings and incorporate it into their structure [12]. In addition, analyzing the assemblage of the phases in the Saguaro, melilite and anorthite (Eu enriched) often surround the pyroxene (Eu depleted) as they are crystallized. This intergrowth of phases and the proximity of the phases would support that the Eu is being incorporated into some phases, preventing it from incorporating into other. Trace elemental analyses of the CAI rims will be evaluated in more detail, as they are complicated by the transient signal being composed of a mixture of mineral phases. Broadly, however, the patterns in the rims of both CAIs are comparable to each other, and for Hoopoe, to the mixed phase patterns in the core (Fig. 3). Other studies have found that CAI rims can be depleted in Ce and Yb [13], however we did not observe these anomalies in the two CAIs discussed here. Given their similarity, the trace elemental analyses of the mixed interior (i.e. core) and rim phases could be interpreted as forming from similar, if not the same, reservoirs. The REE abundance between the rim and core of Hoopoe are also similar, indicating they may have formed from a gas of the same or similar composition. Acknowledgments: We thank the ASU Center for Meteorite Studies for loaning the samples used in this work and Tabb Prissel for his assistance with the analysis. Mouti Al-Hashimi thanks Sam Crossley and Cyrena Goodrich for their help with the LPI SEM training. This work was supported by the LPI Summer Intern Program in Planetary Science and the LPI Cooperative Agreement. References: [1] Connelly J.N. (2012) Science, 338, 651-655. [2] MacPherson G. J. (2014) Treatise on Geochem., 2, 139-179. [3] Grossman L. (1972) GCA, 36, 597-619. [4] Ebel, D.S. (2006) Meteorites and the Early Solar System II (D. S. Lauretta & H. Y. McSween, Eds.) 253-277. [5] MacPherson G. J. (2012) Earth Planet. Sci. Lett., 331-332, 43-54. [6] MacPherson G.J. (2017) GCA, 201, 65-82. [7] Krot A.N. (1995) Meteoritics & Planet. Sci., 30, 748-775. [8] Wark and Lovering (1982) GCA, 46, 2595-2607. [9] Palme H. and Jones A. (2003) Treatise on Geochemistry (H. D. Holland and K. K. Turekian Eds.), 1, 41-61. [10] Grossman L. (1980) Ann. Rev. Earth Planet. Sci., 8, 559-608. [11] Floss C. et al. (1996) GCA, 60, 1975-1997. [12] Mason B. and Martin P. M. (1974) Earth Planet. Sci. Lett., 22, 141-144. [13] Wark B. and Boynton W. V. (2001) Meteoritics & Planet. Sci., 36, 1135-1166.

X Mouti↗
DOE OSTI2026

Emergence of vorticity and viscous stress in finite-scale quantum hydrodynamics

The Madelung equations offer a hydrodynamic description of quantum systems, from single particles to quantum fluids. In this formulation, the probability density is mapped onto the fluid density and the phase is treated as a scalar potential generating the velocity field. As examples of potential flows, quantum fluids described in this way are inherently irrotational, but quantum vortices may arise at discrete points where the phase is undefined. In this paper, starting from this irrotational description of a quantum fluid, a coarse-graining procedure is applied to arrive at a macroscopic description of the quantum fluid in terms of a hierarchy of moments in which the role of velocity is played by a Favre average of the microscopic velocity field. This hierarchy is truncated using an explicit closure derived from an expansion in a finite length scale. The resulting coarse-grained fields are shown to allow for finite vorticity at any point in the fluid. Additionally, it is shown that this vorticity obeys a similar equation to the vorticity equation in classical hydrodynamics and includes a vortex-stretching term. The particular closure employed here also gives rise to a novel stress term in the fluid equations, which in the appropriate limit appears analogous to an artificial viscous stress from computational fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Multidisciplinary Optimization Methods for Aircraft Preliminary Design

This paper describes a research program aimed at improved methods for multidisciplinary design and optimization of large-scale aeronautical systems. The research involves new approaches to system decomposition, interdisciplinary communication, and methods of exploiting coarse-grained parallelism for analysis and optimization. A new architecture, that involves a tight coupling between optimization and analysis, is intended to improve efficiency while simplifying the structure of multidisciplinary, computation-intensive design problems involving many analysis disciplines and perhaps hundreds of design variables. Work in two areas is described here: system decomposition using compatibility constraints to simplify the analysis structure and take advantage of coarse-grained parallelism; and collaborative optimization, a decomposition of the optimization process to permit parallel design and to simplify interdisciplinary communication requirements.

Kroo, Ilan↗
NASA NTRS2004

What Multilevel Parallel Programs do when you are not Watching: A Performance Analysis Case Study Comparing MPI/OpenMP, MLP, and Nested OpenMP

With the current trend in parallel computer architectures towards clusters of shared memory symmetric multi-processors, parallel programming techniques have evolved that support parallelism beyond a single level. When comparing the performance of applications based on different programming paradigms, it is important to differentiate between the influence of the programming model itself and other factors, such as implementation specific behavior of the operating system (OS) or architectural issues. Rewriting-a large scientific application in order to employ a new programming paradigms is usually a time consuming and error prone task. Before embarking on such an endeavor it is important to determine that there is really a gain that would not be possible with the current implementation. A detailed performance analysis is crucial to clarify these issues. The multilevel programming paradigms considered in this study are hybrid MPI/OpenMP, MLP, and nested OpenMP. The hybrid MPI/OpenMP approach is based on using MPI [7] for the coarse grained parallelization and OpenMP [9] for fine grained loop level parallelism. The MPI programming paradigm assumes a private address space for each process. Data is transferred by explicitly exchanging messages via calls to the MPI library. This model was originally designed for distributed memory architectures but is also suitable for shared memory systems. The second paradigm under consideration is MLP which was developed by Taft. The approach is similar to MPi/OpenMP, using a mix of coarse grain process level parallelization and loop level OpenMP parallelization. As it is the case with MPI, a private address space is assumed for each process. The MLP approach was developed for ccNUMA architectures and explicitly takes advantage of the availability of shared memory. A shared memory arena which is accessible by all processes is required. Communication is done by reading from and writing to the shared memory.

Jost, Gabriele↗
NASA NTRS1986

Identification of cosmogenic argon components in Allende by laser microprobe

New techniques are presented for using a laser microprobe to determine the spallation argon systematics of calcium-aluminum inclusions. The Ar-38(s) amounts determined for melilite and anorthite in a coarse-grained inclusion from Allende are 2.9 x 10 to the -8th and 1.3 x 10 to the -8th cc/g, respectively. The ratio of the amounts is consistent with the calcium contents of these two minerals. The Ar-38(s) amount determined for a fine-grained inclusion from Allende is 1.1 x 10 to the -8th cc/g. Calcium and potassium amounts were determined from irradiated samples of the same inclusions so that production of Ar-38 from calcium during the cosmic ray exposure of Allende could be determined for these samples. The production observed was 12.4 + or - 2.1 x 10 to the -8th cc STP Ar-38/g Ca for the coarse-grained inclusion and 9.9 + or - 2.4 cc STP Ar-38/g Ca for the fine-grained inclusion. No evidence of unusual exposure was observed in the two inclusions studied.

Kirschbaum, C.↗
DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗
DOE OSTI2022

ExaSGD: 2022 Kernel Thrust Activities

The Kernel Thrust milestone ADSE22-407 covers the development of device-capable optimization algorithms and solvers technologies required by the ExaSGD project’s software stack in order to solve security-constrained alternating current optimal power flow (SC-ACOPF) problems on emerging exascale architectures. To this extent, in FY22 the main objective of the Kernel Thrust was (i) provide sparse optimization solver that runs efficiently on hardware accelerator devices (i.e., NVIDIA and AMD GPUs) to perform intra-node computations, (ii) strengthen the reliability and increase the performance of the mixed-dense sparse (MDS) solver of HiOp for deployment on the FY22 target architectures, Summit and Crusher, and (iii) increase performance by improving the mathematical algorithm and refining the parallel MPI-based implementation of the coarse-grain parallel solver HiOp-PriDec for capabilities deployment on the FY22 target architectures, Summit and Crusher. This document presents the developments and contributions done by the Kernels Thrust Team in FY22 toward completion of the above-mentioned objectives. These contributions progressed along four main development (sub)thrusts: (1) Design and implementation of a sparse optimization solver for use on hardware accelerators; (2) Improvement of the mathematical algorithm and of the parallel implementation of HiOp-PriDec to ensure readiness and efficient coarse-grain parallelism for FY23 target exascale machine; and (3) Support Software and Application Development Thrusts of the exaSGD project in their deployment of the project’s software stack on AMD- and NVIDIA-based architectures. The development of the sparse optimization solver (thrust 1 above) was new in FY22 and resulted in a new sparse solver in HiOp (available as of version 0.6). The second development thrust was a continuation of the efforts from FY21 and improved the mathematical algorithm and the communication strategy of the HiOp-PriDec solver. The last developement thrust is a large collaborative effort. Namely, the project’s teams from multiple labs (LLNL, PNNL, ORNL, and NREL) performed large-scale demonstration of the ExaSGD software stack, namely the optimization solvers of HiOp interfaced with the modeling front-end ExaGO and the stochastic sampler PowerScenarios. These demonstration efforts solved large-scale instances of the SC-ACOPF challenge problem of medium network sizes (10, 000-bus system) and large number of contingencies on Summit (NVIDIA accelerators) and Crusher (AMD accelerators) systems at ORNL.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

Breccia Formation at a Complex Impact Crater: Slate Islands, Lake Superior, Ontario, Canada

The Slate Islands impact structure is the eroded remnant of a approximately 30-32 km-diameter complex impact structure located in northern Lake Superior, Ontario, Canada. Target rocks are Archean supracrustal and igneous rocks and Proterozoic metavolcanics, metasediments, and diabase. A wide variety of breccias occurs on the islands, many of which contain fragments exhibiting shock metamorphic features. Aphanitic, narrow and inclusion-poor pseudotachylite veins, commonly with more or less parallel boundaries and apophyses branching off them, represent the earliest breccias formed during the compression stage of the impact process. Coarse-grained, polymictic elastic matrix breccias form small to very large, inclusion-rich dikes and irregularly shaped bodies that may contain altered glass fragments. These breccias have sharp contacts with their host rocks and include a wide range of fragment types some of which were transported over minimum distances of approximately 2 km away from the center of the structure. They cut across pseudotachylite veins and contain inclusions of them. Field and petrographic evidence indicate that these polymictic breccias formed predominantly during the excavation and central uplift stages of the impact process. Monomictic breccias, characterized by angular fragments and transitional contacts with their host rocks, occur in parautochthonous target rocks, mainly on the outlying islands of the Slate Islands archipelago. A few contain fragmented and disrupted, coarse-grained, polymictic clastic matrix breccia dikes. This is an indication that at least some of these monomictic breccias formed late in the impact process and that they are probably related to a late crater modification stage. A small number of relatively large occurrences of glass-poor, suevitic breccias occur at the flanks of the central uplift and along the inner flank of the outer ring of the Slate Islands complex crater. A coarse, glass-free, allogenic breccia, containing shatter-coned fragments derived from Proterozoic target rocks (upper target strata), observed at two locations may be analogous to the 'Bunt Breccia' of the Ries crater in Germany. At one of these locations this breccia lies close to a crater suevite deposit. At the other, it overlies parautochthonous, monomictic breccia. The State Islands impact breccias are superbly exposed, much better than breccias in most other terrestrial impact structures. Observations, including those indicative of multiple and and sequential processes, provide insight on how impact breccias form and how they relate to the various phases of the impact process. Eventually they will lead to an improved understanding of planetary impact processes.

Dressler, B. O.↗
NASA NTRS1993

Igneous fractionation and subsolidus equilibration of diogenite meteorites

Diogenites are coarse-grained orthopyroxenite breccias of remarkably uniform major element composition. Most diogenites contain homogeneous pyroxene fragments up to 5 cm across of Wo2En74Fs24 composition. Common minor constituents are chromite, olivine, trolite and metal, while silica, plagioclase, merrillite and diopside are trace phases. Diogenites are generally believed to be cumulates from the eucrite parent body, although their relationship with eucrites remains obscure. It has been suggested that some diogenites are residues after partial melting. I have performed EMPA and INAA for major, minor and trace elements on most diogenites, concentrating on coarse-grained mineral and lithic clasts in order to elucidate their igneous formation and subsequent metamorphic history. Major element compositions of diogenites are decoupled from minor and trace element compositions; the latter record an igneous fractionation sequence that is not preserved in the former. Low equilibration temperatures indicate that major element diffusion continued long after crystallization. Diffusion coefficients for trivalent and tetravalent elements in pyroxene are lower than those of divalent elements. Therefore, major element compositions of diogenites may represent means of unknown portions of a cumulate homogenized by diffusion, while minor and trace elements still yield information on their igneous history. The scale of major element equilibration is unknown, but is likely to be on the order of a few cm. Therefore, the diogenite precursors may have consisted largely of cm-sized, igneously zoned orthopyroxene grains, which were subsequently annealed during slow cooling, obliterating major element zoning but preserving minor and trace incompatible element zoning.

Mittlefehldt, David W.↗
NASA NTRS1999

Detection of Crystalline Hematite Mineralization on Mars by the Thermal Emission Spectrometer: Evidence for Near-surface Water

The Thermal Emission Spectrometer (TES) instrument on the Mars Global Surveyor (MGS) mission has discovered a remarkable accumulation of crystalline hematite ((alpha-Fe2O3) that covers an area with very sharp boundaries approximately 350 by 350-750 km in size centered near 2 S latitude between 0 and 5 W longitude (Sinus Meridiani). Crystalline hematite is uniquely identified by the presence of fundamental vibrational absorption features centered near 300, 450, and >525/cm, and by the absence of silicate fundamentals in the 1000/cm region. Spectral features resulting from atmospheric CO2, dust, and water ice were removed using a radiative transfer model. The spectral properties unique to Sinus Meridiani were emphasized by removing the average spectrum of the surrounding region. The depth and shape of the hematite fundamental bands show that the hematite is crystalline and relatively coarse grained (>5-10 micron). Diameters up to and greater than 100s of micrometers are permitted within the instrumental noise and natural variability of hematite spectra. Hematite particles <5-10 micron in diameter (either as an unpacked or hard-packed powders) fail to match the TES spectra. The spectrally-derived areal abundance of hematite varies with particle size from approximately 10% for particles >30 micron in diameter to 40-60% for unpacked 10 micron powders. The hematite in Sinus Meridiani is thus distinct from the fine-grained (diameter <5-10 micron), red, crystalline hematite considered, on the basis of visible, near-IR data, to be a minor spectral component in Martian bright regions like Olympus-Amazonis. Sinus Meridiani hematite is closely associated with a smooth, layered, friable surface that is interpreted to be sedimentary in origin. This material may be the uppermost surface in the region, indicating that it could be a late-stage sedimentary unit, or it could be a layered portion of the heavily cratered plains units. We consider five possible mechanisms for the formation of coarse-grained, crystalline hematite. These processes fall into two classes depending on whether they require a significant amount of near-surface water: (1) chemical precipitation that includes origin by (a) precipitation from oxygenated, Fe-rich water (iron formations), (b) hydrothermal extraction and crystal growth.

Christensen, P. R.↗