Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data similarity”

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

First spectral measurement of deuterium-tritium fusion γ rays in inertial fusion experiments

Here, the R-matrix analysis of A = 5 nuclear systems has been partially validated by applying the technique to the 5 Li system and comparing the predicted γ-ray spectrum with historical data. R-matrix analysis of the similar 5 He system was then used to predict the γ-ray spectral shape for the deuterium-tritium (DT) reaction. The resulting spectra have been used in the analysis of DT implosions on the Omega laser where the γ-ray interaction rate was measured by a gas Cherenkov detector. Comparison of predictions to experiment confirmed the presence of both 16.75 and ≈ 13 MeV γ-ray contributions; analysis, using R-matrix spectra, yielded a ratio of γ-ray emission from a transition to the intermediate excited state to that from a transition to the ground state of (2.1 ± 0.4):1, substantiating the first spectral measurement of the DT fusion γ ray in an inertial fusion environment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Time-Integrated Neutrino Source Searches with 10 Years of IceCube Data

This Letter presents the results from pointlike neutrino source searches using ten years of IceCube data collected between April 6, 2008 and July 10, 2018. We evaluate the significance of an astrophysical signal from a pointlike source looking for an excess of clustered neutrino events with energies typically above ~1 TeV among the background of atmospheric muons and neutrinos. We perform a full-sky scan, a search within a selected source catalog, a catalog population study, and three stacked Galactic catalog searches. The most significant point in the northern hemisphere from scanning the sky is coincident with the Seyfert II galaxy NGC 1068, which was included in the source catalog search. The excess at the coordinates of NGC 1068 is inconsistent with background expectations at the level of 2.9σ after accounting for statistical trials from the entire catalog. The combination of this result along with excesses observed at the coordinates of three other sources, including TXS 0506+056, suggests that, collectively, correlations with sources in the northern catalog are inconsistent with background at 3.3σ significance. The southern catalog is consistent with background. Finally these results, all based on searches for a cumulative neutrino signal integrated over the 10 years of available data, motivate further study of these and similar sources, including time-dependent analyses, multimessenger correlations, and the possibility of stronger evidence with coming upgrades to the detector.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Data-driven approach to parameterize SCAN + U for an accurate description of 3 d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the DFT + U method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. Here we develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-U correction, and (3) DFT/DFT + U compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for 3d transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve SCAN/SCAN + U compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize DFT + U performance for other target physical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electrochemical Formation of Li-M-(M')-Si Phases Using Multivalent Electrolyte Salt Additives

Lithium-rich silicides (Li 15 Si 4 ), formed during the electrochemical lithiation of silicon, show high reactivity with electrolyte components that contribute to capacity decay, formal lithium loss, and low coulombic efficiency. Recently, the reactivity of lithium silicides was found to be suppressed by substituting a multivalent cation (i.e. Mg, Ca) for lithium that results in the room temperature formation of a ternary Li-M-Si phase. In this study, we explored a range of multivalent electrolyte salt additives (M = Ni, Cu, La, Ce, Sr, Ba, and Ca-Mg mixed salt) in a lithium-ion cell configuration and identified a room temperature electrochemical route to the formation of new ternary and quaternary lithium silicides. Using this method, both nickel and copper salts were found to plate onto the silicon electrode surface upon lithiation. Further, based on refined synchrotron XRD data, multivalent cations with an ionic radius similar to Na (~1.03 Å) or smaller can be inserted electrochemically into a formally cation-deficient Li 15 Si 4 host lattice to form new ternary (or quartenary) phases. The electrochemical synthesis of a new quaternary Li-M-M’-Si phase represents a facile route to preparing and scaling materials isostructural to the Heusler phase and electron-precise Li 14 MgSi 4 phase that results in enhanced cycling and calendar life performance.

25 ENERGY STORAGE↗
DOE OSTI2020

Petrographic and spectral study of hydrothermal mineralization in drill core from Hawaii: A potential analog to alteration in the martian subsurface

Abstract Continuous rock core was collected for 1764 m (5786’) on the Pohakuloa Army Training base near the center of the big island of Hawaii. The core traverses basaltic lava flows from the volcano's shield-building phase, and perched aquifers and higher temperature groundwaters were encountered. The collected samples record water-rock interactions of basaltic materials in a setting that may be a model for groundwater interactions on Mars. We collected visible and infrared point spectra of materials in the lowest portion of the core, where alteration was noted to become more prominent. We identified three types of phyllosilicate spectral signatures and three types of zeolites. The phyllosilicates show similarity to those identified on Mars using data from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM). Based on the field survey, 25 depths were selected for sampling and petrographic analysis of thin sections. The spectral data and thin section work have a strong agreement in the types of materials identified by the two different techniques. Both the spectral and petrographic data indicate low to moderate temperature geothermal alteration occurred in the lower half of the core. The field spectra are a useful reconnaissance tool for selecting mineralogic diversity for subsequent higher resolution and more time-consuming laboratory analysis.

Geochemistry & Geophysics↗
DOE OSTI2026

IER-620 CED-3b: Experiment Execution Summary for the Pulsed-Neutron Die-Away Experiments (PNDA) with Propylene Glycol and Mobilmet 423

There is a strong need for new benchmarks to validate neutron thermal scattering laws (TSLs). Lawrence Livermore National Laboratory (LLNL) has designed a Pulsed-Neutron Die Away (PNDA) testbed for this purpose. The experiment has a deuterium-tritium (D-T) neutron generator that impinges a 10 -4 s, mono-energetic pulse of 14.1 MeV neutrons on a target sample. After the pulse, the neutron population moderates and establishes a thermal equilibrium within the sample, with a fundamental spatial mode and characteristic decay-time eigenvalue, ⍺. The ⍺ eigenvalue can be extracted from the experimental measurements of the time-dependent neutron flux coming off the surface of the sample and can then be used as an integral parameter (similar to k eff ) to validate nuclear data involved with neutron migration, thermalization, and absorption. For moderating materials and geometric configurations, the ⍺ eigenvalue is heavily dependent on thermal neutron scattering of the target material. For that reason, a PNDA experiment can have a higher sensitivity to TSLs than is commonly available with the k eff parameter in critical experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multiscale Modeling of Vinyl-Addition Polynorbornenes: The Effect of Stereochemistry

Vinyl-addition polynorbornenes are candidates for designing high-performance polymers due to unique characteristics, which include a high glass transition temperature associated with a rigid backbone. Recent studies have established that the processability and properties of these polymers can be fine-tuned by using targeted substitutions. However, synthesis with different catalysts results in materials with distinct properties, potentially due to the presence of various stereoisomers that are difficult to quantify experimentally. Herein, we develop all-atom models of polynorbornene oligomers based on classical force fields and density functional theory. To establish the relationship between chemical architecture, chain conformations, and melt structure, we perform detailed molecular dynamics simulations with the fine-tuned atomistic force field and propose simpler coarse-grained descriptions to address the high molecular weight limit. All-atom simulations of oligomers suggest high glass transition temperatures in the range of 550–600 K. In the melt state (800 K), meso chains form highly rigid extended coils (C∞≈11) with amorphous structural characteristics similar to the X-ray diffraction data observed in the literature. In contrast, simulations with racemo chains predict highly helical tubular chain conformations that could promote assembly into crystalline structures.

Polymer Science↗
DOE OSTI2020

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

Here we present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature thereby integrating over the thermal seed photon flux from the disk surface. The method either assumes the ion and electron temperatures are equal (1T) or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. This new approach leads to radiative efficiency values sime65% above Novikov–Thorne, with a larger fraction of total cooling in the corona as compared to simulations performed using the original target-temperature cooling function. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach, which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

CLASSY III. The Properties of Starburst-driven Warm Ionized Outflows

We report the results of analyses of galactic outflows in a sample of 45 low-redshift starburst galaxies in the COS Legacy Archive Spectroscopic SurveY (CLASSY), augmented by five additional similar starbursts with Cosmic Origins Spectrograph (COS) data. The outflows are traced by blueshifted absorption lines of metals spanning a wide range of ionization potential. The high quality and broad spectral coverage of CLASSY data enable us to disentangle the absorption due to the static interstellar medium (ISM) from that due to outflows. We further use different line multiplets and doublets to determine the covering fraction, column density, and ionization state as a function of velocity for each outflow. We measure the outflow's mean velocity and velocity width, and find that both correlate in a highly significant way with the star formation rate, galaxy mass, and circular velocity over ranges of four orders of magnitude for the first two properties. We also estimate outflow rates of metals, mass, momentum, and kinetic energy. We find that, at most, only about 20% of silicon created and ejected by supernovae in the starburst is carried out in the warm phase we observe. The outflows' mass-loading factor increases steeply and inversely with both circular and outflow velocity (log–log slope ~–1.6), and reaches ~10 for dwarf galaxies. We find that the outflows typically carry about 10%–100% of the momentum injected by massive stars and about 1%–20% of the kinetic energy. We show that these results place interesting constraints on, and new insights into, models and simulations of galactic winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The S-PLUS Fifth Data Release: Over 4500 Square Degrees of the Southern Sky and a Multicolor View of the Hydra and Antlia Galaxy Clusters

Abstract We present the fifth data release (DR5) of the Southern Photometric Local Universe Survey (S-PLUS), covering 4592 deg 2 across 2491 fields. Observations were conducted with the T80-South, a Brazilian robotic telescope equipped with the Javalambre 12-filter system, containing five broadband and seven narrowband filters. Data products feature FITS images and extensive catalogs containing fluxes, magnitudes, and shape parameters for over 113 million detections. In addition, several value-added catalogs are provided, offering photometric redshifts (photo- zs ), object classifications, masks, and extinction coefficients. For the first time, this release includes coverage of 110 deg 2 along the Galactic disk, facilitating new research into Galactic structure and stellar populations. The release also provides full coverage of the Hydra Supercluster and numerous other nearby clusters with improved data reduction and calibration, enhancing photo- z accuracy, which is vital for large-scale structure studies. A preliminary analysis of the Hydra and Antlia galaxy clusters up to 5 × R 200 yields an updated catalog of 1706 cluster members based on both spectroscopic and high-quality photo- zs . Our photometric data shows that both clusters have a similar proportion of galaxies with an H α excess relative to their clustercentric distance, though Hydra has a higher fraction near its center. Additionally, the spatial distribution of all objects in our sample highlights a bridge connecting both clusters. We verify that S-PLUS DR5 provides a solid foundation for future scientific investigations, ranging from solar system studies to cosmology.

Vinicius Rodrigues de Lima, Erik [Universidade de ↗
DOE OSTI2022

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series

Abstract. Improving the quantification of soil thermal and physical properties is key to achieving a better understanding and prediction of soil hydro-biogeochemical processes and their responses to changes in atmospheric forcing. Obtaining such information at numerous locations and/or over time with conventional soil sampling is challenging. The increasing availability of low-cost, vertically resolved temperature sensor arrays offers promise for improving the estimation of soil thermal properties from temperature time series, and the possible indirect estimation of physical properties. Still, the reliability and limitations of such an approach need to be assessed. In the present study, we develop a parameter estimation approach based on a combination of thermal modeling, sliding time windows, Bayesian inference, and Markov chain Monte Carlo simulation to estimate thermal diffusivity and its uncertainty over time, at numerous locations and at an unprecedented vertical spatial resolution (i.e., down to 5 to 10 cm vertical resolution) from soil temperature time series. We provide the necessary framework to assess under which environmental conditions (soil temperature gradient, fluctuations, and trend), temperature sensor characteristics (bias and level of noise), and deployment geometries (sensor number and position) soil thermal diffusivity can be reliably inferred. We validate the method with synthetic experiments and field studies. The synthetic experiments show that in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d the proposed method provides reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. Finally, the application of the developed approach to field data indicates significant repeatability in results and similarity with independent measurements, as well as promise in using a sliding time window to estimate temporal changes in soil thermal diffusivity, as needed to potentially capture changes in bulk density or water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

On the circulation, water mass distribution, and nutrient concentrations of the western Chukchi Sea

Substantial amounts of nutrients and carbon enter the Arctic Ocean from the Pacific Ocean through the Bering Strait, distributed over three main pathways. Water with low salinities and nutrient concentrations takes an eastern route along the Alaskan coast, as Alaskan Coastal Water. A central pathway exhibits intermediate salinity and nutrient concentrations, while the most nutrient-rich water enters the Bering Strait on its western side. Towards the Arctic Ocean, the flow of these water masses is subject to strong topographic steering within the Chukchi Sea with volume transport modulated by the wind field. In this contribution, we use data from several sections crossing Herald Canyon collected in 2008 and 2014 together with numerical modelling to investigate the circulation and transport in the western part of the Chukchi Sea. We find that a substantial fraction of water from the Chukchi Sea enters the East Siberian Sea south of Wrangel Island and circulates in an anticyclonic direction around the island. This water then contributes to the high-nutrient waters of Herald Canyon. The bottom of the canyon has the highest nutrient concentrations, likely as a result of addition from the degradation of organic matter at the sediment surface in the East Siberian Sea. The flux of nutrients (nitrate, phosphate, and silicate) and dissolved inorganic carbon in Bering Summer Water and Winter Water is computed by combining hydrographic and nutrient observations with geostrophic transport referenced to lowered acoustic Doppler current profiler (LADCP) and surface drift data. Even if there are some general similarities between the years, there are differences in both the temperature–salinity and nutrient characteristics. To assess these differences, and also to get a wider temporal and spatial view, numerical modelling results are applied. According to model results, high-frequency variability dominates the flow in Herald Canyon. This leads us to conclude that this region needs to be monitored over a longer time frame to deduce the temporal variability and potential trends.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
DOE OSTI2021

Data-Driven Learning for the Mori--Zwanzig Formalism: A Generalization of the Koopman Learning Framework

A theoretical framework which unifies the conventional Mori--Zwanzig formalism and the approximate Koopman learning of deterministic dynamical systems from noiseless observation is presented. In this framework, the Mori--Zwanzig formalism, developed in statistical mechanics to tackle the hard problem of construction of reduced-order dynamics for high-dimensional dynamical systems, can be considered as a natural generalization of the Koopman description of the dynamical system. We next show that, similar to the approximate Koopman learning methods, data-driven methods can be developed for the Mori--Zwanzig formalism with Mori's linear projection operator. We have developed two algorithms to extract the key operators, the Markov and the memory kernel, using time series of a reduced set of observables in a dynamical system. We have adopted the Lorenz `96 system as a test problem and solved for the above operators. These operators exhibit complex behaviors, which are unlikely to be captured by traditional modeling approaches in Mori--Zwanzig analysis. The nontrivial generalized fluctuation-dissipation relationship, which relates the memory kernel with the two-time correlation statistics of the orthogonal dynamics, was numerically verified as a validation of the solved operators. Here we present numerical evidence that the generalized Langevin equation, a key construct in the Mori--Zwanzig formalism, is more advantageous in predicting the evolution of the reduced set of observables than the conventional approximate Koopman operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Development of a Representative Molten Chloride Fast Reactor Model to Assess the Impact of Nuclear Data

The SCALE code system was employed to conduct a preliminary investigation of nuclear data impacts for a fast spectrum molten chloride salt reactor. A computationally effective depletion model that is representative of the reactor system was successfully developed and used to conduct fuel depletion simulations. Development of this model draws from the SLICE method that was developed at Oak Ridge National Laboratory to enable generation of fuel compositions for an advanced reactor core at equilibrium operation. Eigenvalue uncertainty calculations using the ENDF/B-VII.1 nuclear data library were performed for the reactor in the fresh fuel state and an irradiated fuel state. It was determined that the primary driver of eigenvalue uncertainty was the uncertainty in the 235U (n, 𝛾) cross section. Uncertainty calculation results from this study were compared to results available for a different fast system, a sodium-cooled fast reactor, to confirm similarities and identify differences with respect to nuclear data impacts between the two fast advanced reactor systems.

Hirji, Rakim [Georgia Institute of Technology]↗
DOE OSTI2020

Locality-Aware Scheduling for Scalable Heterogeneous Environments

Heterogeneous computing promise boost performance of scientific applications by allowing massively parallel execution of computational tasks. However, manually managing extremely heterogeneous, multi-device systems is complicated and may result in sub-optimal performance. Specifically, data management is an extremely challenging problem on multi-device systems. In this work, we introduce two locality-aware schedulers for the Minos Computing Library (MCL), an asynchronous, task-based programming model and runtime for extremely heterogeneous systems. The first scheduler implements a pure locality-aware algorithm to maximize data reuse, though it might incur in ”hot-spots” that limit system utilization. The second scheduler mitigates this drawback by dynamically targeting between locality-awareness and system utilization based on the current workload and available computing devices. Our results show that locality-awareness greatly benefit applications that exhibit data reuse, providing up to 6.9x and 7.9x over the original MCL scheduler and equivalent OpenCL implementations, respectively. Moreover, our schedulers introduce negligible overhead compared with the original MCL scheduler and achieve similar performance for applications that don’t benefit from data locality.

Architecture, co-design, Task-based programming mo↗
DOE OSTI2023

Establishing metrics to quantify spatial similarity in spherical and red blood cell distributions

As computational power increases and systems with millions of red blood cells can be simulated, it is important to note that varying spatial distributions of cells may affect simulation outcomes. Since a single simulation may not represent the ensemble behavior, many different configurations may need to be sampled to adequately assess the entire collection of potential cell arrangements. In order to determine both the number of distributions needed and which ones to run, we must first establish methods to identify well-generated, randomly placed cell distributions and to quantify distinct cell configurations. We utilize metrics to assess (1) the presence of any underlying structure to the initial cell distribution and (2) similarity between cell configurations. We propose the use of the radial distribution function to identify long-range structure in a cell configuration and apply it to a randomly distributed and structured set of red blood cells. To quantify spatial similarity between two configurations, we make use of the Jaccard index, and characterize sets of red blood cell and sphere initializations. As an extension to our work submitted to the International Conference on Computational Science, we significantly increase our data set size from 72 to 1048 cells, include a similar set of studies using spheres, compare the effects of varying sphere size, and utilize the Jaccard index distribution to probe sets of extremely similar configurations. Our results show that the radial distribution function can be used as a metric to determine long-range structure in both distributions of spheres and RBCs. We determine that the ideal case of spheres within a cube versus bi-concave shaped cells within a cylinder affects the shape of the Jaccard index distributions, as well as the range of Jaccard values, showing that both the shape of particle and the domain may play a role. Furthermore, we also find that the distribution is able to capture very similar configurations through Jaccard index values greater than 95% when appending several nearly identical configurations into the data set.

59 BASIC BIOLOGICAL SCIENCES↗