Search NASA⌕ Search

SEARCH · Search NASA

Results for “Data Abstractions and Types”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

High‐Pressure and High‐Temperature Single‐Crystal Elasticity of Cr‐Pyrope: Implications for the Density and Seismic Velocity of Subcontinental Lithospheric Mantle

Abstract Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp 71.0 Alm 12.6 Sps 0.7 Grs 3.5 Uvr 12.2 ) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus ( K S 0 ), shear modulus ( G 0 ), their pressure ((∂ K S /∂ P ) T and (∂ G /∂ P ) T ) and temperature (( ∂K S /∂T ) P and ( ∂G/∂T ) P ) derivatives, K S 0 = 167.7(8) GPa, G 0 = 91.5(5) GPa, (∂ K S /∂ P ) T = 4.3(1), (∂ G /∂ P ) T = 1.4(1), ( ∂K S /∂T ) P = − 0.0175(1) GPa/K and ( ∂G/∂T ) P = − 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity, ρ / V P and ρ / V S are much more sensitive to the compositional change of the SCLM.

Xu, Jingui↗
DOE OSTI2022

Measuring Cosmological Parameters with Type Ia Supernovae in redMaGiC Galaxies

Abstract Current and future cosmological analyses with Type Ia supernovae (SNe Ia) face three critical challenges: (i) measuring the redshifts from the SNe or their host galaxies; (ii) classifying the SNe without spectra; and (iii) accounting for correlations between the properties of SNe Ia and their host galaxies. We present here a novel approach that addresses each of these challenges. In the context of the Dark Energy Survey (DES), we analyze an SN Ia sample with host galaxies in the redMaGiC galaxy catalog, a selection of luminous red galaxies. redMaGiC photo- z estimates are expected to be accurate to σ Δ z /(1+ z ) ∼ 0.02. The DES-5YR photometrically classified SN Ia sample contains approximately 1600 SNe, and 125 of these SNe are in redMaGiC galaxies. We demonstrate that redMaGiC galaxies almost exclusively host SNe Ia, reducing concerns relating to classification uncertainties. With this subsample, we find similar Hubble scatter (to within ∼0.01 mag) using photometric redshifts in place of spectroscopic redshifts. With detailed simulations, we show that the bias due to using redMaGiC photo- z s on the measurement of the dark energy equation of state w is up to Δ w ∼ 0.01–0.02. With real data, we measure a difference in w when using the redMaGiC photo- z s versus the spec- z s of Δ w = 0.005. Finally, we discuss how SNe in redMaGiC galaxies appear to comprise a more standardizable population, due to a weaker relation between color and luminosity ( β ) compared to the DES-3YR population by ∼5 σ . These results establish the feasibility of performing redMaGiC SN cosmology with photometric survey data in the absence of spectroscopic data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Developing a Consistent Travel-Time Framework for Comparing Three-Dimensional Velocity Models for Seismic Location Accuracy

Abstract Location algorithms have historically relied on simple, one-dimensional (1D) velocity models for fast seismic event locations. 1D models are generally used as travel-time lookup tables, one for each seismic phase, with travel-times pre-calculated for event distance and depth. These travel-time lookup tables are extremely fast to use and this fast computational speed makes them the preferred type of velocity model for operational needs. Higher-dimensional (i.e., three-dimensional—3D) seismic velocity models are becoming readily available and provide more accurate event locations over 1D models. The computational requirements of these 3D models tend to make their operational use prohibitive. Additionally, comparing location accuracy for 3D seismic velocity models tends to be problematic, as each model is determined using different ray-tracing algorithms. Attempting to use a different algorithm than the one used to develop a model usually results in poor travel-time prediction. We demonstrate and test a framework to create first-P and first-S 3D travel-time correction surfaces using an open-source framework ( PCalc + GeoTess , https://www.sandia.gov/salsa3d/software/geotess ) that easily stores 3D travel-time and uncertainty data. This framework produces fast travel-time and uncertainty predictions and overcomes the ray-tracing algorithm hurdle because the lookup tables can be generated using the exact ray-tracing algorithm that is preferred for a model.

3D velocity models↗
DOE OSTI2023

Accurate model of the projected velocity distribution of galaxies in dark matter haloes

ABSTRACT We present a per cent-level accurate model of the line-of-sight velocity distribution of galaxies around dark matter haloes as a function of projected radius and halo mass. The model is developed and tested using synthetic galaxy catalogues generated with the UniverseMachine run on the Multi-Dark Planck 2 N-body simulations. The model decomposes the galaxies around a cluster into three kinematically distinct classes: orbiting, infalling, and interloping galaxies. We demonstrate that: (1) we can statistically distinguish between these three types of galaxies using only projected line-of-sight velocity information; (2) the halo edge radius inferred from the line-of-sight velocity dispersion is an excellent proxy for the three-dimensional halo edge radius; and (3) we can accurately recover the full velocity dispersion profile for each of the three populations of galaxies. Importantly, the velocity dispersion profiles of the orbiting and infalling galaxies contain five independent parameters – three distinct radial scales and two velocity dispersion amplitudes – each of which is correlated with mass. Thus, the velocity dispersion profile of galaxy clusters has inherent redundancies that allow us to perform non-trivial systematics checks from a single data set. We discuss several potential applications of our new model for detecting the edge radius and constraining cosmology and astrophysics using upcoming spectroscopic surveys.

Astronomy & Astrophysics↗
NASA NTRSDate not supplied

The Nucleosynthetic Heritage of Bennu

The OSIRIS-REx mission sampled the B-type asteroid Bennu and delivered 121.6 g of material to Earth in September 2023. Nucleosynthetic isotope variations are powerful tracers that can be used to establish relationships between different meteorite groups and the returned samples. They can also constrain the samples’ nucleosynthetic heritage, as well as mixing and reservoir formation in the protoplanetary disk. Nucleosynthetic variations reflect the heterogeneous distribution of isotopically distinct dust in the disk, and as such, they provide evidence for the formation of meteorite parent bodies from isotopically distinct dust reservoirs. They are thus also a means to identify exogenous clasts in Bennu materials. Here, we report Ti and Fe isotope data obtained from Bennu samples and carbonaceous chondrites. We also present Ti isotope data for hand-picked matrix from CV and CR chondrites. An isotopic dataset was obtained on aggregate sample OREX-803015-100 consisting of particles typically <500 μm; therefore, its composition represents an average for Bennu materials. A 20.66 mg aliquot was dissolved at Washington University, St. Louis. From this, 5 mg was analysed for Ti and Fe isotopes on a NeptunePlus MC-ICPMS at ETH Zurich; data from a separate aliquot, analysed at Lawrence Livermore National Laboratory, are reported in a companion abstract. Our results demonstrate a close relationship of Bennu to CI and Ryugu samples and therefore a formation outside the snowline, as hypothesized by the mission. Furthermore, the 50 Ti data from the CV and CR matrix fall consistently below the CI and Bennu values, indicating the presence of a uniform matrix mixed with ⁓0.1–0.7% of refractory inclusions in CIs and Bennu.

M. Schönbächler↗
DOE OSTI2020

The completed SDSS-IV extended baryon oscillation spectroscopic survey: pairwise-inverse probability and angular correction for fibre collisions in clustering measurements

ABSTRACT The completed extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalogues contain redshifts of 344 080 quasars at 0.8 < z < 2.2, 174 816 luminous red galaxies between 0.6 < z < 1.0, and 173 736 emission-line galaxies over 0.6 < z < 1.1 in order to constrain the expansion history of the Universe and the growth rate of structure through clustering measurements. Mechanical limitations of the fibre-fed spectrograph on the Sloan telescope prevent two fibres being placed closer than 62 arcsec in a single pass of the instrument. These ‘fibre collisions’ strongly correlate with the intrinsic clustering of targets and can bias measurements of the two-point correlation function resulting in a systematic error on the inferred values of the cosmological parameters. We combine the new techniques of pairwise-inverse probability and the angular upweighting (PIP+ANG) to correct the clustering measurements for the effect of fibre collisions. Using mock catalogues, we show that our corrections provide unbiased measurements, within data precision, of both the projected $\rm {\mathit{ w}_p}\left(\mathit{ r}_p\right)$ and the redshift-space multipole ξ(ℓ = 0, 2, 4)(s) correlation functions down to $0.1\, h^{-1}{\rm Mpc}$, regardless of the tracer type. We apply the corrections to the eBOSS DR16 catalogues. We find that, on scales $s\gtrsim 20\, h^{-1}{\rm Mpc}$ for ξℓ, as used to make baryon acoustic oscillation and large-scale redshift-space distortion measurements, approximate methods such as nearest-neighbour upweighting are sufficiently accurate given the statistical errors of the data. Using the PIP method, for the first time for a spectroscopic program of the Sloan Digital Sky Survey, we are able to successfully access the one-halo term in the clustering measurements down to $\sim 0.1\, h^{-1}{\rm Mpc}$ scales. Our results will therefore allow studies that use the small-scale clustering to strengthen the constraints on both cosmological parameters and the halo occupation distribution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Image feature extraction and galaxy classification: a novel and efficient approach with automated machine learning

ABSTRACT In this work, we explore the possibility of applying machine learning methods designed for 1D problems to the task of galaxy image classification. The algorithms used for image classification typically rely on multiple costly steps, such as the point spread function deconvolution and the training and application of complex Convolutional Neural Networks of thousands or even millions of parameters. In our approach, we extract features from the galaxy images by analysing the elliptical isophotes in their light distribution and collect the information in a sequence. The sequences obtained with this method present definite features allowing a direct distinction between galaxy types. Then, we train and classify the sequences with machine learning algorithms, designed through the platform Modulos AutoML. As a demonstration of this method, we use the second public release of the Dark Energy Survey (DES DR2). We show that we are able to successfully distinguish between early-type and late-type galaxies, for images with signal-to-noise ratio greater than 300. This yields an accuracy of $86{{\ \rm per\ cent}}$ for the early-type galaxies and $93{{\ \rm per\ cent}}$ for the late-type galaxies, which is on par with most contemporary automated image classification approaches. The data dimensionality reduction of our novel method implies a significant lowering in computational cost of classification. In the perspective of future data sets obtained with e.g. Euclid and the Vera Rubin Observatory, this work represents a path towards using a well-tested and widely used platform from industry in efficiently tackling galaxy classification problems at the peta-byte scale.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Iodine-Xenon Dating: Sensitive Chronometer for Reprocessing in the Primitive Solar System

The I-Xe chronometer is based upon decay of I-129 to Xe-129 in the early Solar System. Recent comparison of I-Xe system in individual mineral separates from twelve different meteorites with independent Pb-Pb data has demonstrated that I-Xe clock is a reliable sensitive chronometer when applied to a single mineral system. Since most iodine hosts are secondary minerals, the I-Xe clock generally records post-formational processing, providing the information on early meteorite evolution. Absolute I-Xe ages can be found by normalization using the measured I-Xe and Pb-Pb ages of Acapulco phosphate (4.557 plus or minus 0.002 Ga). Absolute ages for the I-Xe internal standards Shallow water and Bjurbole, 4.566 plus or minus 0.002 Ga and 4.565 plus or minus 0.003 Ga, respectively, provide absolute I-Xe ages for all other samples. The I-Xe age of bulk meteorite is meaningful and interpretable only when the carrier of primordial iodine is a major mineral phase (e. g., enstatite chondrites). Using the "monomineral" approach, separated phases from the Richardton H5 chondrite provide a case history of post-formational alteration in this object. This work applies the I-Xe chronometer to determine the times of reprocessing of selected minerals in single meteorite types. A preliminary account of this work was recently reported. Additional information is contained in the original extended abstract.

Pravdivtseva, O. V.↗
DOE OSTI2021

Risk assessment of engineering diseases of embankment–bridge transition section for railway in permafrost regions

Abstract The embankment–bridge transition section (EBTS) is one of the zones where railway diseases occur frequently in permafrost regions. Disease risk assessment of EBTSs can provide guidance for maintenance. In this study, considering the engineering geological conditions, climate characteristics, and embankment structure types along the Qinghai–Tibet Railway (QTR) as well as based on the disease inventory of the QTR from 2010 to 2019, the logistic regression (LR), support vector machine (SVM), and combination‐weight‐based gay relation analysis (GRA) were used for disease risk assessment of the EBTSs along the QTR in permafrost regions. The results indicate that the LR and SVM models have a better capability for EBTS disease prediction than the GRA model, and the SVM model can select more disease samples in relatively larger regions than the LR model. Based on the SVM and LR models, the risk level of EBTSs is divided into four classes: low‐ (29.9%), moderate‐ (39.6%), high‐ (22.1%), and very high (8.4%) risk. Finally, we selected 272 EBTSs in high‐ and very‐high‐risk classes for key observation during the maintenance of the QTR in permafrost regions. This study provides a reference for the risk assessment of railways built in permafrost regions using data‐driven methods.

Zhang, Saize↗
DOE OSTI2018

Nature of Decahydro-closo-decaborate Anion Reorientations in an Ordered Alkali-Metal Salt: Rb 2 B 10 H 10

The ordered monoclinic phase of the alkali-metal decahydro-closo-decaborate salt Rb 2 B 10 H 10 was found to be stable from about 250 K all the way up to an order-disorder phase transition temperature of ~762 K. The broad temperature range for this phase allowed for a detailed quasielastic neutron scattering (QENS) and nuclear magnetic resonance (NMR) study of the protypical B 10 H 10 2- anion reorientational dynamics. The QENS and NMR combined results are consistent with an anion reorientational mechanism comprised of two types of rotational jumps expected from the anion geometry and lattice structure, namely, more rapid 90 degree jumps around the anion C 4 symmetry axis (e.g., with correlation frequencies of ~2.6 x 10 10 s -1 at 530 K) combined with order of magnitude slower orthogonal 180 degree reorientational flips (e.g., ~3.1 x 10 9 s -1 at 530 K) resulting in an exchange of the apical H (and apical B) positions. Each latter flip requires a concomitant 45 degree twist around the C 4 symmetry axis to preserve the ordered Rb 2 B 10 H 10 monoclinic structural symmetry. This result is consistent with previous NMR data for ordered monoclinic Na 2 B 10 H 10 , which also pointed to two types of anion reorientational motions. The QENS-derived reorientational activation energies are 197(2) and 288(3) meV for the C 4 fourfold jumps and apical exchanges, respectively, between 400 and 680 K. Below this temperature range, NMR (and QENS) both indicate a shift to significantly larger reorientational barriers, for example, 485(8) meV for the apical exchanges. Finally, subambient diffraction measurements identify a subtle change in the Rb 2 B 10 H 10 structure from monoclinic to triclinic symmetry as the temperature is decreased from around 250 to 210 K.

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

Two years of optical and NIR observations of the superluminous supernova UID 30901 discovered by the UltraVISTA SN survey

ABSTRACT We present deep optical and near-infrared photometry of UID 30901, a superluminous supernova (SLSN) discovered during the UltraVISTA survey. The observations were obtained with VIRCAM (YJHKs) mounted on the VISTA telescope, DECam (griz) on the Blanco telescope, and SUBARU Hyper Suprime-Cam (HSC; grizy). These multiband observations comprise +700 d making UID 30901 one of the best photometrically followed SLSNe to date. The host galaxy of UID 30901 is detected in a deep HST F814W image with an AB magnitude of 27.3 ± 0.2. While no spectra exist for the SN or its host galaxy, we perform our analysis assuming z = 0.37, based on the photometric redshift of a possible host galaxy found at a projected distance of 7 kpc. Fitting a blackbody to the observations, the radius, temperature, and bolometric light curve are computed. We find a maximum bolometric luminosity of 5.4 ± 0.34 × 1043 erg s−1. A flattening in the light curve beyond 600 d is observed and several possible causes are discussed. We find the observations to clearly favour an SLSN type I, and plausible power sources such as the radioactive decay of 56Ni and the spin-down of a magnetar are compared to the data. We find that the magnetar model yields a good fit to the observations with the following parameters: a magnetic field B = 1.4 ± 0.3 × 1014 G, spin period of P = 6.0 ± 0.1 ms, and ejecta mass $M_{\mathrm{ ej}} = 11.9^{+4.8}_{-6.4} \,\mathrm{ M}_{\odot }$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

ESS-DIVE Reporting Format for Dataset Package Metadata

ESS-DIVE’s (Environmental Systems Science Data Infrastructure for a Virtual Ecosystem) dataset metadata reporting format is intended to compile information about a dataset (e.g., title, description, funding sources) that can enable reuse of data submitted to the ESS-DIVE data repository. The files contained in this dataset include instructions (dataset_metadata_guide.md and README.md) that can be used to understand the types of metadata ESS-DIVE collects. The data dictionary (dd.csv) follows ESS-DIVE’s file-level metadata reporting format and includes brief descriptions about each element of the dataset metadata reporting format. This dataset also includes a terminology crosswalk (dataset_metadata_crosswalk.csv) that shows how ESS-DIVE’s metadata reporting format maps onto other existing metadata standards and reporting formats.Data contributors to ESS-DIVE can provide this metadata by manual entry using a web form or programmatically via ESS-DIVE’s API (Application Programming Interface). A metadata template (dataset_metadata_template.docx or dataset_metadata_template.pdf) can be used to collaboratively compile metadata before providing it to ESS-DIVE.Since being incorporated into ESS-DIVE’s data submission user interface, ESS-DIVE’s dataset metadata reporting format, has enabled features like automated metadata quality checks, and dissemination of ESS-DIVE datasets onto other data platforms including Google Dataset Search and DataCite.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Using long‐term data from a whole ecosystem warming experiment to identify best spring and autumn phenology models

Abstract Predicting vegetation phenology in response to changing environmental factors is key in understanding feedbacks between the biosphere and the climate system. Experimental approaches extending the temperature range beyond historic climate variability provide a unique opportunity to identify model structures that are best suited to predicting phenological changes under future climate scenarios. Here, we model spring and autumn phenological transition dates obtained from digital repeat photography in a boreal Picea ‐ Sphagnum bog in response to a gradient of whole ecosystem warming manipulations of up to +9°C, using five years of observational data. In spring, seven equally best‐performing models for Larix utilized the accumulation of growing degree days as a common driver for temperature forcing. For Picea , the best two models were sequential models requiring winter chilling before spring forcing temperature is accumulated. In shrub, parallel models with chilling and forcing requirements occurring simultaneously were identified as the best models. Autumn models were substantially improved when a CO 2 parameter was included. Overall, the combination of experimental manipulations and multiple years of observations combined with variation in weather provided the framework to rule out a large number of candidate models and to identify best spring and autumn models for each plant functional type.

Schädel, Christina↗
DOE OSTI2018

Towards in-situ high precision local material velocity measurements in lattice materials under dynamic compression

The study of dynamic compression of materials requires high precision spatial and temporal resolution to infer the material state at any given time. In ultra-fast imaging experiments such as those performed at the Dynamic Compression Sector at the Advanced Photon Source, the current CCD frame rate is less than the X-ray pulse rate necessitating that 4 CCDs be used to image a single scintillator screen to obtain the same projection of transmitted or diffracted beam at time intervals synchronized with the X-ray beam. The use of multiple-camera setups is often employed for a variety of reasons including enhancing the temporal resolution through interleaving frames or increasing the number of viewpoints, for example. Physical corrections are made to have each CCD approximately viewing the same field of view, however we present corrections that can be made after the data collection to further refine and be able to compare movements and changes in the sample between independent CCDs. These software corrections enable higher precision and higher fidelity comparison of position changes due to sample deformation. Furthermore, we present a method to correct images (or frames) as a demonstration toward high fidelity displacement measurements in these types of ultra-fast multi-camera setup experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Deploy threading in Nalu solver stack

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary code in the ExaWind project is Nalu, which is an unstructured-grid solver for the acousticallyincompressible Navier-Stokes equations, and mass continuity is maintained through pressure projection. The model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as re-initialization of matrices and re-computation of preconditioners is required at every time step In this Milestone, we examine the effect of threading on the solver stack performance against flat-MPI results obtained from previous milestones using Haswell performance data full-turbine simulations. Whereas the momentum equations are solved only with the Trilinos solvers, we investigate two algebraic-multigrid preconditioners for the continuity equations: Trilinos/Muelu and HYPRE/BoomerAMG. These two packages embody smoothed-aggregation and classical Ruge-Stiiben AMG methods, respectively. In our FY18 Q2 report, we described our efforts to improve setup and solve of the continuity equations under flat-MPI parallelism. While significant improvement was demonstrated in the solve phase, setup times remained larger than expected. Starting with the optimized settings described in the Q2 report, we explore here simulation performance where OpenMP threading is employed in the solver stack. For Trilinos, threading is acheived through the Kokkos abstraction where, whereas HYPRE/BoomerAMG employs straight OpenMP. We examined results for our mid-resolution baseline turbine simulation configuration (229M DOF). Simulations on 2048 Haswell cores explored the effect of decreasing the number of MPI ranks while increasing the number of threads. Both HYPRE and Trilinos exhibited similar overal solution times, and both showed dramatic increases in simulation time in the shift from MPI ranks to OpenMP threads. This increase is attributed to the large amount of work per MPI rank starting at the single-thread configuration. Decreasing MPI ranks, while increasing threads, may be increasing simulation time due to thread synchronization and start-up overhead contributing to the latency and serial time in the model. These result showed that an MPI+OpenMP parallel decomposition will be more effective as the amount per MPI rank computation per MPI rank decreases and the communication latency increases. This idea was demonstrated in a strong scaling study of our low-resolution baseline model (29M DOF) with the Trilinos-HYPRE configuration. While MPI-only results showed scaling improvement out to about 1536 cores, engaging threading carried scaling improvements out to 4128 cores — roughly 7000 DOF per core. This is an important result as improved strong scaling is needed for simulations to be executed over sufficiently long simulated durations (i.e., for many timesteps). In addition to threading work described above, the team examined solver-performance improvements by exploring communication-overhead in the HYPRE-GMRES implementation through a communicationoptimal- GMRE algorithm (CO-GMRES), and offloading compute-intensive solver actions to GPUs. To those ends, a HYPRE mini-app was allow us to easily test different solver approaches and HYPRE parameter settings without running the entire Nalu code. With GPU acceleration on the Summitdev supercomputer, a 20x speedup was achieved for the overall preconditioner and solver execution time for the mini-app. A study on Haswell processors showed that CO-GMRES provides benefits as one increases MPI ranks.

17 WIND ENERGY↗