Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pattern clustering”

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.

734 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Methane emission hotspots in a boreal forest-fen mosaic potentially linked to deep taliks

Permafrost thaw is transforming boreal forests into mosaics of wetlands and drier uplands. Topographic controls on hydrological and ecological conditions impact methane (CH 4 ) fluxes, contributing to uncertainty in local and regional CH 4 budgets and underlying drivers. The objective of this study was to explore CH 4 fluxes and their drivers in a transitioning boreal forest-fen ecosystem (Goldstream Valley, Alaska, USA). This landscape is characterized by thawing discontinuous permafrost and heterogeneous mosaics of fens, collapse-scar channels, and small mounds of permafrost soils. From a survey in July 2021, observed chamber CH4 fluxes included fen areas with intermediate to very high emissions (29.8–635.3 mg CH 4 m −2 d −1 ), clustered locations with CH 4 uptake (−2.11 to −0.7 mg CH 4 m −2 d −1 ), and three anomalous emission hotspots (342.4–772.4 mg CH 4 m −2 d −1 ) that were located near samples with lower emissions. Some surface and near-surface variables partially explained the spatial variation in CH 4 flux. Log-transformed CH 4 flux had a positive linear relationship with soil moisture at 20 cm depth ( R 2 = 0.31, p -value < 1e-5) and negative linear relationships with microtopography ( R 2 = 0.13, p -value < 0.006) and slope ( R 2 = 0.28, p -value < 2e-5). Methane emissions generally occurred in flat, wet, graminoid-dominated fens, whereas CH 4 uptake occurred on permafrost mounds dominated by feather mosses and woody vegetation. However, the CH 4 hotspots occurred on drier, slightly sloped locations with low or undetectable near-surface methanogen abundance, suggesting that CH 4 was produced in deeper soils. When the hotspot samples were omitted, log-transformed CH 4 flux had a positive linear relationship with near-surface methanogen abundance ( R 2 = 0.29, p -value = 0.0023), and stronger linear relationships with soil moisture, slope, and soil macronutrient concentrations. Our findings suggest that some CH 4 emission hotspots could arise from CH 4 in deep taliks. The inference that methanogenesis occurs in deep taliks was strengthened by the identification of intrapermafrost taliks across the study area using low-frequency geophysical induction. This study assesses surface spatial heterogeneity in the context of subsurface permafrost conditions and highlights the complexity of CH 4 flux patterns in transitioning forest-wetland ecosystems. To better inform regional CH 4 budgets, further research is needed to understand the spatial distribution of terrestrial CH 4 hotspots and to resolve their surface, near-surface, and subsurface drivers.

boreal↗
NASA NTRS2016

Scattering and Diffraction of Electromagnetic Radiation: An Effective Probe to Material Structure

Scattered electromagnetic waves from material bodies of different forms contain, in an intricate way, precise information on the intrinsic, geometrical and physical properties of the objects. Scattering theories, ever deepening, aim to provide dependable interpretation and prediction to the complicated interaction of electromagnetic radiation with matter. There are well-established multiple-scattering formulations based on classical electromagnetic theories. An example is the Generalized Multi-particle Mie-solution (GMM), which has recently been extended to a special version ̶ the GMM-PA approach, applicable to finite periodic arrays consisting of a huge number (e.g., >>106) of identical scattering centers [1]. The framework of the GMM-PA is nearly complete. When the size of the constituent unit scatterers becomes considerably small in comparison with incident wavelength, an appropriate array of such small element volumes may well be a satisfactory representation of a material entity having an arbitrary structure. X-ray diffraction is a powerful characterization tool used in a variety of scientific and technical fields, including material science. A diffraction pattern is nothing more than the spatial distribution of scattered intensity, determined by the distribution of scattering matter by way of its Fourier transform [1]. Since all linear dimensions entered into Maxwell's equations are normalized by wavelength, an analogy exists between optical and X-ray diffraction patterns. A large set of optical diffraction patterns experimentally obtained can be found in the literature [e.g., 2,3]. Theoretical results from the GMM-PA have been scrutinized using a large collection of publically accessible, experimentally obtained Fraunhofer diffraction patterns. As far as characteristic structures of the patterns are concerned, theoretical and experimental results are in uniform agreement; no exception has been found so far. Closely connected with the spatial distribution of scattered intensities are cross sections, such as for extinction, scattering, absorption, and radiation pressure, as a critical type of key quantity addressed in most theoretical and experimental studies of radiative scattering. Cross sections predicted from different scattering theories are supposed to be in general agreement. For objects of irregular shape, the GMM-PA solutions can be compared with the highly flexible Discrete Dipole Approximation (DDA) [4,5] when dividing a target to no more than ~106 unit cells. Also, there are different ways to calculate the cross sections in the GMM-PA, providing an additional means to examine the accuracy of the numerical solutions and to unveil potential issues concerning the theoretical formulations and numerical aspects. To solve multiple scattering by an assembly of material volumes through classical theories such as the GMM-PA, the radiative properties of the component scatterers, the complex refractive index in particular, must be provided as input parameters. When using a PA to characterize a material body, this involves the use of an adequate theoretical tool, an effective medium theory, to connect Maxwell's phenomenogical theory with the atomistic theory of matter. In the atomic theory, one regards matter as composed of interacting particles (atoms and molecules) embedded in the vacuum [6]. However, the radiative properties of atomic-scaled particles are known to be substantially different from bulk materials. Intensive research efforts in the fields of cluster science and nanoscience attempt to bridge the gap between bulk and atom and to understand the transition from classical to quantum physics. The GMM-PA calculations, which place virtually no restriction on the component-particle size, might help to gain certain insight into the transition.

Xu, Yu-Lin↗
NASA NTRS2017

Washboard Terrain on Pluto

Washboard texture or patterning consists of fields of parallel to sub-parallel ridges typically spaced ~1-2 km crest to crest and a few 100 m in amplitude (Fig. 4a in Moore et al., 2016, Science, 351, 1284-1293). For the most part, underlying topography can be easily discerned. We will refer to discrete, well-bounded patches of these landforms as Washboard Terrain (WT). WT is observed to occur along the rim, and just beyond the rim, of Sputnik basin from the West to NNW. Where it is seen in high-resolution data, it has clearly defined limits, beyond which it would be able to be seen if it were there. WT doesn't occur at very low latitudes or very high latitudes (ranging from 22degN to 62degN). WT seems to occur most conspicuously on relatively level, gently sloping terrain. It is restricted to elevations between approximately 2 km to less than +1.5 km (i.e. not at high elevations). The most noticeable regional aspect of the area in which WT occurs is the sinuous valley network, which is suspected to have been formed, or at least substantially modified, by glaciation. WT also appears to occur mainly on an intermediate-albedo reddish material, where seen in enhanced color data. Where it occurs in level terrain, WT tends to trend ENE - there doesn't seem to be a strong local control of its orientation in response to valley drainage directions. WT can display a greater range of orientations where it occurs in higher-relief (not higher elevation) settings such as spurs. WT appears superposed on very ancient landscapes, but is itself cratered locally by clusters of small (approximately 1-3 km) craters, which may be secondaries. This implies that WT may be intermediate in age. Of several working hypotheses, we currently provisionally favor that WT may be akin to terrestrial recessional moraines (or de Geer moraines) associated with the retreat of a higher stand of N2 glaciation that once overfilled Sputnik basin. These putative moraine features may owe their spacing to superseasonal retreat on Milankovitch timescales of approximately 1 Ma. If this hypothesis has validity, then perhaps the intermediate-albedo reddish material may be akin to ground moraine deposits.

Moore, Jeffrey M.↗
NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗
NASA NTRS1993

Geochemical evidences for two chondritic-like cometary or asteroidal impacts before and at the K/T boundary

A number of geological and palaeontological evidences support multiple impacts of cometary showers within a short time (approximately 1-3 Ma) and their connection with mass extinctions. Observations include clustered crater ages, stratigraphic horizons of impact ejecta closely spaced in time, and evidence for stepwise mass extinctions spanning intervals of 1-3 Ma. For the K/T boundary, three candidates, Popigai, Manson, and Yucatan, have been proposed as impact craters. Two distinct strata at the K/T boundary in western North America have been interpreted as evidence for two sequential impacts. If multiple impacts occurred within a time span of about 1 Ma then multiple Ir enrichments should be observed. DSDP Hole 577B on the Shatsky Plateau in the northern Pacific at K/T time is the first site. Samples contain approximately greater than 97 percent CaCO3, which exhibit clear chemical signals associated with asteroidal/cometary impact. Ir, Fe, and Cr data are presented. From the Th-normalized data, two satellite peaks below the major peak at 78 cm and 81 cm of 577B-1-4 are clearly shown. The major Ir peak (K/T boundary) is at 72 cm. Fe and Cr, from C1-like impactor ejecta fallout, also show two peaks at the same positions. For hole 738C on the southern Kerguelen Plateau, Ir values reach a peak concentration of 18 ppb in the clay layer at 96.0-96.2 cm in section 20R-5, and gradually tail off. In the sample 115 cm above the boundary, Ir concentrations have still not reached background levels. From the Ir peak downward to the lowermost sample analyzed at 102 cm, the Ir concentration is still as high as 1.7 ppb. From the Th-normalized data, we observe a small Ir/Th peak at 100-101 cm. Though this peak is within the error margin, the trend is clear. Fe and Cr exhibit the same pattern. The third case is Hole 690C on the Queen Maud Ridge. Again, the Ir/Th plot indicates the strong possibility of satellite peaks at approximately 52 cm. The main peak is at 39-40 cm. For the Stevns Klint K/T boundary layers, the stratification of trace elements appears threefold with peak concentrations in sublayers A1, A3, and B2 for different element groups, including Ir. C1 ratios for many siderophile elements found in combined layers III and IV, corresponding to layers A, B, C, and D, strongly support the impact hypothesis. Also, multiple Ir anomalies in the K/T section at Lattengebirge, Bavarian Alps are reported. Recent works on Ni-rich spinels and Ir at the K/T boundaries clearly establish cometary/asteroidal impacts at the K/T boundary. Lastly, cometary showers can explain the enhanced Ir contents over approximately a 1 Ma interval in Gubbio shales.

Liu, Y.-G.↗
NASA NTRS2000

Shape-dependent control of cell growth, differentiation, and apoptosis: switching between attractors in cell regulatory networks

Development of characteristic tissue patterns requires that individual cells be switched locally between different phenotypes or "fates;" while one cell may proliferate, its neighbors may differentiate or die. Recent studies have revealed that local switching between these different gene programs is controlled through interplay between soluble growth factors, insoluble extracellular matrix molecules, and mechanical forces which produce cell shape distortion. Although the precise molecular basis remains unknown, shape-dependent control of cell growth and function appears to be mediated by tension-dependent changes in the actin cytoskeleton. However, the question remains: how can a generalized physical stimulus, such as cell distortion, activate the same set of genes and signaling proteins that are triggered by molecules which bind to specific cell surface receptors. In this article, we use computer simulations based on dynamic Boolean networks to show that the different cell fates that a particular cell can exhibit may represent a preprogrammed set of common end programs or "attractors" which self-organize within the cell's regulatory networks. In this type of dynamic network model of information processing, generalized stimuli (e.g., mechanical forces) and specific molecular cues elicit signals which follow different trajectories, but eventually converge onto one of a small set of common end programs (growth, quiescence, differentiation, apoptosis, etc.). In other words, if cells use this type of information processing system, then control of cell function would involve selection of preexisting (latent) behavioral modes of the cell, rather than instruction by specific binding molecules. Importantly, the results of the computer simulation closely mimic experimental data obtained with living endothelial cells. The major implication of this finding is that current methods used for analysis of cell function that rely on characterization of linear signaling pathways or clusters of genes with common activity profiles may overlook the most critical features of cellular information processing which normally determine how signal specificity is established and maintained in living cells. Copyright 2000 Academic Press.

Review↗
DOE OSTI2022

Building a Framework to Understand Transition Metals' Behavior in Euxinic Conditions (Final Technical Report)

This project focuses first and foremost on metal sulfide geochemistry and mineralogy as controlled by a complex matrix of environmental factors. The principal investigator’s group aim to illuminate the metal-sulfide reaction mechanisms, rates, and pathways through systematic experimentation and data collection and analyzing the relationships between the characteristics of the produced metal sulfide solid-phase/aqueous complexes and the environmental factors. This understanding is essential for obtaining a full picture of the complex cycling patterns of single or multi metal species in sulfidic environments ranging from deep-see basins, hydrothermal vents, inland seas, terrestrial water bodies, to engineered remediation systems. The specific goal of this past project was to illuminate the reaction mechanisms and kinetics of metal anions and sulfide in mixed metal cation-metal anion-sulfide systems under various aqueous conditions (which resembled a range of naturally occurring euxinic settings). For the period of this contract, we investigated the molybdenum-iron-sulfide system, with an emphasis on the conditions that caused solid phase formation. We focused on quantifying the mobility/sequestration of molybdenum under each experimental condition and identified the changes of valence states for each involved element (i.e., Mo, Fe, and S) in the precipitate. We also proposed pathways for the electron transfer that occurred in aqueous chemistry. The major analytical tools used for this study include UV-visible light spectroscopy, transmission electron microscopy, X-ray photoelectron spectroscopy, and synchrotron-based X-ray absorption spectroscopy (access gained through facility proposals to the Canadian Light Source). The biggest finding of this project was that besides pH, the iron-sulfur chemistry has a dominant control of the thiolation kinetics and subsequent reduction of Mo(VI), which are likely prerequisites of molybdenum sequestration in anoxic conditions. The results have been written up as manuscript by the end of this project (see Phillips et al.). The experimental results of this project may be critical for advancing our understanding of (1) basic chemistry involving transition metals and reduced sulfur species, (2) the validity of certain geochemical proxies, and (3) the stability and evolution of euxinic geochemical environments. It is noted that the basic results obtained through this project also have implications for Mo-S cluster-based catalyst development in inorganic chemistry and materials sciences.

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

Luciola Hypertelescope Space Observatory

Luciola is a large (one kilometer) "multi-aperture densified-pupil imaging interferometer", or "hypertelescope" employing many small apertures, rather than a few large ones, for obtaining direct snapshot images with a high information content. A diluted collector mirror, deployed in space as a flotilla of small mirrors, focuses a sky image which is exploited by several beam-combiner spaceships. Each contains a pupil densifier micro-lens array to avoid the diffractive spread and image attenuation caused by the small sub-apertures. The elucidation of hypertelescope imaging properties during the last decade has shown that many small apertures tend to be far more efficient, regarding the science yield, than a few large ones providing a comparable collecting area. For similar underlying physical reasons, radio-astronomy has also evolved in the direction of many-antenna systems such as the proposed Low Frequency Array having hundreds of thousands of individual receivers . With its high limiting magnitude, reaching the mv=30 limit of HST when 100 collectors of 25cm will match its collecting area, high-resolution direct imaging in multiple channels, broad spectral coverage from the 1200 Angstrom ultra-violet to the 20 micron infra-red, apodization, coronagraphic and spectroscopic capabilities, the proposed hypertelescope observatory addresses very broad and innovative science covering different areas of ESA s Cosmic Vision program. In the initial phase, a focal spacecraft covering the UV to near IR spectral range of EMCCD photon-counting cameras ( currently 200 to 1000nm), will image details on the surface of many stars, as well as their environment, including multiple stars and clusters. Spectra will be obtained for each resel. It will also image neutron star, black-hole and micro-quasar candidates, as well as active galactic nuclei, quasars, gravitational lenses, and other Cosmic Vision targets observable with the initial modest crowding limit. With subsequent upgrade missions, the spectral coverage can be extended from 120nm to 20 microns, using four detectors carried by two to four focal spacecraft. The number of collector mirrors in the flotilla can also be increased from 12 to 100 and possibly 1,000. The imaging and spectroscopy of habitable exoplanets in the mid infra-red then becomes feasible once the collecting area reaches 6m2 , using a specialized mid infra-red focal spacecraft. Calculations ( Boccaletti et al., 2000) have shown that hypertelescope coronagraphy has unequalled sensitivity for detecting, at mid infra-red wavelengths, faint exoplanets within the exo-zodiacal glare. Later upgrades will enable the more difficult imaging and spectroscopy of these faint objects at visible wavelengths, using refined techniques of adaptive coronagraphy (Labeyrie. & Le Coroller, 2004). Together, the infra-red and visible spectral data carry rich information on the possible presence of life. The close environment of the central black-hole in the Milky Way will be imageable with unprecedented detail in the near infra-red . Cosmological imaging of remote galaxies at the limit of the known universe is also expected, from the ultra-violet to the near infra-red, following the first upgrade, and with greatly increasing sensitivity through successive upgrades. These areas will indeed greatly benefit from the upgrades, in terms of dynamic range, limiting complexity of the objects to be imaged, size of the elementary Direct Imaging Field , and limiting magnitude, approaching that of an 8-meter space telescope when 1000 apertures of 25cm are installed. Similar gains will occur for addressing fundamental problems in physics and cosmology, particularly when observing neutron stars and black holes, single or binary, including the giant black holes, with accretion disks and jets, in active galactic nuclei beyond the Milky Way. Gravitational lensing and micro-lensing patterns, including time-variable patterns and perhaps millisecond lensing flasheshich may be beamed by diffraction from sub-stellar masses at sub-parsec distances (Labeyrie, 1994) , will also be observable initially in the favourable cases, and upgrades will greatly improve the number of observable objects. The observability of gravitational waves emitted by binary lensing masses, in the form of modulated lensing patterns, is a debated issue ( Ragazzoni et al., 2003) but will also become addressable observationally. The technology readiness of Luciola approaches levels where low-orbit testing and stepwise implementation will become feasible in the 2015-2025 time frame. For the following decades beyond 2020, once accurate formation flying techniques will be mastered, much larger hypertelescopes such as the proposed 100km Exo-Earth Imager and the 100,000 km Neutron Star Imager should also become feasible. Luciola is therefore also seen as a precursor toward such very powerful instruments.

Labeyrie, Antoine↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
NASA NTRS2003

Patterns of global gene expression in rat skeletal muscle during unloading and low-intensity ambulatory activity

Physical inactivity and unloading lead to diverse skeletal muscle alterations. Our goal was to identify the genes in skeletal muscle whose expression is most sensitive to periods of unloading/reduced physical activity and that may be involved in triggering initial responses before phenotypic changes are evident. The ability of short periods of physical activity/loading as an effective countermeasure against changes in gene expression mediated by inactivity was also tested. Affymetrix microarrays were used to compare mRNA levels in the soleus muscle under three experimental treatments (n = 20-29 rats each): 12-h hindlimb unloading (HU), 12-h HU followed by 4 h of intermittent low-intensity ambulatory and postural activity (4-h reloading), and control (with ambulatory and postural activity). Using a combination of criteria, we identified a small set of genes (approximately 1% of 8,738 genes on the array or 4% of significant expressed genes) with the most reproducible and largest responses to altered activity. Analysis revealed a coordinated regulation of transcription for a large number of key signaling proteins and transcription factors involved in protein synthesis/degradation and energy metabolism. Most (21 of 25) of the gene expression changes that were downregulated during HU returned at least to control levels during the reloading. In surprising contrast, 27 of 38 of the genes upregulated during HU remained significantly above control, but most showed trends toward reversal. This introduces a new concept that, in general, genes that are upregulated during unloading/inactivity will be more resistant to periodic reloading than those genes that are downregulated. This study reveals genes that are the most sensitive to loading/activity in rat skeletal muscle and indicates new targets that may initiate muscle alterations during inactivity.

NASA Discipline Musculoskeletal↗
DOE OSTI2024

Acceleration of Graph Neural Network-Based Prediction Models in Chemistry via Co-Design Optimization on Intelligence Processing Units

Atomic structure prediction and associated property calculations are the bedrock of chemical physics. Since high-fidelity ab initio modeling techniques for computing the structure and properties can be prohibitively expensive, this motivates the development of machine-learning (ML) models that make these predictions more efficiently. Training graph neural networks over large atomistic databases introduces unique computational challenges such as the need to process millions of small graphs with variable size and support communication patterns that are distinct from learning over large graphs such as social networks. We demonstrate a novel hardware-software co-design approach to scale up the training of atomistic graph neural networks (GNN) for structure and property prediction. First, to eliminate redundant computation and memory associated with alternative padding techniques and to improve throughput via minimizing communication, we formulate the effective coalescing of the batches of variable-size atomistic graphs as the bin packing problem and introduce a hardware-agnostic algorithm to pack these batches. In addition, we propose hardware-specific optimizations including a planner and vectorization for the gather-scatter operations targeted for Graphcore’s Intelligence Processing Unit (IPU), as well as model-specific optimizations such as merged communication collectives and optimized softplus. Putting these all together, we demonstrate the effectiveness of the proposed co-design approach by providing an implementation of a well-established atomistic GNN on the Graphcore IPUs. We evaluate the training performance on multiple atomistic graph databases with varying degrees of graph counts, sizes and sparsity. Here, we demonstrate that such a co-design approach can reduce the training time of atomistic GNNs and can improve the performance by up to 1.5× compared to the baseline implementation of the model on the IPUs. Additionally, we compare our IPU implementation with a Nvidia GPU-based implementation and show that our atomistic GNN implementation on the IPUs can run 1.8× faster on average compared to the execution time on the GPUs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

SMC 2021 : Analyzing Resource Utilization and User Behavior on Titan Supercomputer

Resource utilization statistics of submitted jobs on a supercomputer can help us understand how users from various scientific domains use HPC platforms and better design a job scheduler. We explore to generate insight regarding workload distribution and usage pattern domains from job scheduler trace, GPU failure information, and project-specific information collected from Titan supercomputer. Furthermore, we want to know how the scheduler performance varies over time and how the users' scheduling behavior changes following a system failure. These observations have the potential to provide valuable insight, which is helpful to prepare for system failures. These practices will help us develop and apply novel machine learning algorithms in understanding system behavior, requirement, and better scheduling of HPC systems. There are two datasets, RUR and GPU. RUR: This dataset is the job scheduler traces collected from the Titan supercomputerfrom 01/01/2015 to 07/31/2019 (2015.csv - 2019.csv). These were collected usingResource Utilization Report (RUR), a Cray-developed resource-usage data collectionand reporting system. It contains the usage information of its critical resources (CPU,Memory, GPU, and I/O) of each running job on Titan during that period [2]. ProjectAreas: Every job is associated with a project ID. TheProjectAreas.csvdatasetprovides a mapping of the project ID to its domain science. GPU: There have been some hardware-related issues in the GPUs in Titan that caused some GPUs to fail, sometimes irrecoverably during some job runs. This dataset provides information regarding these failures during the execution of the submitted jobs. GPUs on Titan are uniquely identified by a serial number (SN), and they are installed in a location. A GPU can be installed in a location, then removed from that location following a failure, and then re-installed in a different location after fixing the problem. If the failure can't be recovered, the GPU might be removed entirely from Titan. There are two prominent types of failures that resulted in the removal of GPUs from Titan: Double Bit Error (DBE) and Out of the Bus (OTB). The dataset (gc_full.csv) has the following fields: 1. SN : Serial number of a GPU 2. location : The location where it is installed 3. insert : The time when it was inserted into that location 4. remove : The time when it was removed from that location 5. duration : Amount of time the GPU spent in this location 6. out : If the device was taken out entirely w/o a re-installment into a new location. 7. event : If the GPU was taken out entirely, the reason for its removal. To learn more about this dataset, please refer to the git repositoryhttps://github.com/olcf/TitanGPULifeand the related publication [1]. References [1] George Ostrouchov, Don Maxwell, Rizwan A Ashraf, Christian Engelmann, MallikarjunShankar, and James H Rogers. Gpu lifetimes on titan supercomputer: Survival analysisand reliability. InSC20: International Conference for High Performance Computing,Networking, Storage and Analysis, pages 1-14. IEEE, 2020. [2] Feiyi Wang, Sarp Oral, Satyabrata Sen, and Neena Imam. Learning from five-yearresource-utilization data of titan system. In2019 IEEE International Conference onCluster Computing (CLUSTER), pages 1-6. IEEE, 2019.

42 ENGINEERING↗
DOE OSTI2020

Evaluating Energy Efficiency Opportunities from Connected and Automated Vehicle Deployments Coupled with Shared Mobility in California

Connected and Automated Vehicles (CAVs) can be considered to be a disruptive transportation technology, with the potential to significantly improve overall transportation system efficiency; however, CAVs may increase induce vehicle miles traveled (VMT) and bring on greater energy consumption. Further, shared mobility is another disruptive transportation event that is reshaping our travel patterns. The primary goal of this project was to extensively collect data from vehicles and associated infrastructure equipped with CAV technologies from both real-world experiments and simulation studies mainly deployed in California, and develop a comprehensive framework for evaluating energy efficiency opportunities from large-scale (e.g., statewide) introduction of CAVs and a wide deployment of shared mobility systems in a variety of scenarios. To quantify the combined impact of CAV and shared mobility on travel behavior, traffic performance, and energy efficiency, a unique mesoscopic simulation-based model was developed for mobility and energy efficiency evaluation considering these disruptive transportation technologies. As a complement to existing studies on nationwide evaluation of CAVs’ energy impacts, this project was focused on data collection efforts and CAV applications under congested traffic environments that are frequently experienced on a massive scale across the major metropolitan areas in California. Extensive real-world data collection supplemented with simulation studies were conducted to cover a variety of CAV and shared mobility scenarios, particularly on scenarios less-explored in the existing research. Another key component of this project was to consider the interaction between different CAV technologies and shared mobility models, and the compound effect on energy efficiency. A comprehensive modeling suite was developed to quantify the impact of new mobility technologies on travel behavior and traffic performance. The developed modeling framework includes an energy intensity module, mode choice module and activity generation module that are integrated into an agent-based BEAM simulation platform to perform impact analysis based on a variety of scenarios. In addition, the RouteE model has been upgraded to incorporate the impact of CAVs on traffic flow, VMT and energy intensity, using micro-simulation data collected from both freeways and urban arterials. A novel fundamental influencing factor (FIF) mode choice model was developed to link CAV and shared mobility components with travel behaviors, and adapted into the BEAM-centered model framework. A statewide energy inventory was constructed under various CAV technology deployment scenarios by incorporating datasets and models for predicting vehicle market share and vehicle usage, which are tightly associated with the penetration of shared mobility systems. Based applying this modeling suite to a calibrated network in Riverside California, it was found that cooperative automated driving in general will improve mobility, but automated vehicles, even when deployed in a shared autonomous fleet, will likely bring an increase of VMT (up to 36%) due to mode shifts and deadheading. Ride-hailing vehicles typically have better energy efficiency and a higher share of electric vehicles, which helps offset the negative impact from VMT increases when estimating the system-level energy consumption. In general, simulation results show a 6% increase in energy consumption for the scenarios with an increasing shift to ride-hailing modes. The statewide analysis based on the National Household Travel Survey (NHTS) sample data is consistent with the findings from the Riverside network and validate the developed clustering-prediction modeling methodology. The outcomes from this project will help close the knowledge gap on recognizing the potential performance and energy impacts of a broad deployment of CAV and shared mobility technologies across a wide range of roadway infrastructure with varying levels of congestion. Results from this project: 1) will support policymakers in steering CAV development and deployment towards an energy favorable direction; 2) reduce uncertainties in estimating energy saving opportunities from new mobility technologies and services; 3) increase the confidence of CAV technology investors both on the infrastructure side (i.e., transportation agencies) and on the vehicle side (i.e., OEMs); and 4) expedite the deployment of energy-efficient CAV and shared mobility applications.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Evidence for autotrophic growth of purple sulfur bacteria using pyrite as electron and sulfur source

ABSTRACT Purple sulfur bacteria (PSB) are capable of anoxygenic photosynthesis via oxidizing reduced sulfur compounds and are considered key drivers of the sulfur cycle in a range of anoxic environments. In this study, we show that Allochromatium vinosum (a PSB species) is capable of autotrophic growth using pyrite as the electron and sulfur source. Comparative growth profile, substrate characterization, and transcriptomic sequencing data provided valuable insight into the molecular mechanisms underlying the bacterial utilization of pyrite and autotrophic growth. Specifically, the pyrite-supported cell cultures (“py”’) demonstrated robust but much slower growth rates and distinct patterns from their sodium sulfide-amended positive controls. Up to ~200-fold upregulation of genes encoding various c - and b -type cytochromes was observed in “py,” pointing to the high relevance of these molecules in scavenging and relaying electrons from pyrite to cytoplasmic metabolisms. Conversely, extensive downregulation of genes related to LH and RC complex components indicates that the electron source may have direct control over the bacterial cells’ photosynthetic activity. In terms of sulfur metabolism, genes encoding periplasmic or membrane-bound proteins (e.g., FccAB and SoxYZ) were largely upregulated, whereas those encoding cytoplasmic proteins (e.g., Dsr and Apr groups) are extensively suppressed. Other notable differentially expressed genes are related to flagella/fimbriae/pilin(+), metal efflux(+), ferrienterochelin(−), and [NiFe] hydrogenases(+). Characterization of the biologically reacted pyrite indicates the presence of polymeric sulfur. These results have, for the first time, put the interplay of PSB and transition metal sulfide chemistry under the spotlight, with the potential to advance multiple fields, including metal and sulfur biogeochemistry, bacterial extracellular electron transfer, and artificial photosynthesis. IMPORTANCE Microbial utilization of solid-phase substrates constitutes a critical area of focus in environmental microbiology, offering valuable insights into microbial metabolic processes and adaptability. Recent advancements in this field have profoundly deepened our knowledge of microbial physiology pertinent to these scenarios and spurred innovations in biosynthesis and energy production. Furthermore, research into interactions between microbes and solid-phase substrates has directly linked microbial activities to the surrounding mineralogical environments, thereby enhancing our understanding of the relevant biogeochemical cycles. Our study represents a significant step forward in this field by demonstrating, for the first time, the autotrophic growth of purple sulfur bacteria using insoluble pyrite (FeS2) as both the electron and sulfur source. The presented comparative growth profiles, substrate characterizations, and transcriptomic sequencing data shed light on the relationships between electron donor types, photosynthetic reaction center activities, and potential extracellular electron transfer in these organisms capable of anoxygenic photosynthesis. Furthermore, the findings of our study may provide new insights into early-Earth biogeochemical evolutions, offering valuable constraints for understanding the environmental conditions and microbial processes that shaped our planet’s history.

59 BASIC BIOLOGICAL SCIENCES↗