Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “total variation”

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 487 records · Page 27

Measurement of the QE Map

113 MHz SRF gun utilized CsK 2 Sb and Na 2 KSb cathodes, which were irradiated with a pulsed green laser with a wavelength of 0.53 microns. The drive laser generated 78 kHz pulses with an acousto-optic modulator selecting the desired number of pulses reaching the cathode. The operator can select the aperture size, defining the laser spot diameter on the cathode, and the power level, which is controlled by a rotating polarizer. The laser power reaching the cathode can be measured with a meter based on a photodiode. The quantum efficiency (QE) of a photocathode serves as a primary indicator of its state. Small QE renders the cathode unusable when the desired bunch charge cannot be achieved. Strong variations in the quantum efficiency on the cathode surface might lead to the growth of the beam emittance due to the space charge effects. We have created a routine to measure the distribution of the QE on the cathode. The operator selects the aperture size, the central spot location, and the scan ranges for both planes. The scan starts from the central spots and proceeds along a spiral-like trajectory, making N rounds. The number of rounds can be chosen, and the total number of measured points is (2N+1) 2 . The laser spot movement is shown in Fig. 1.

43 PARTICLE ACCELERATORS↗

Genomic and Phenotypic Characterization of Shiga Toxin-Producing Escherichia albertii Strains Isolated from Wild Birds in a Major Agricultural Region in California

Escherichia albertii is an emerging foodborne pathogen. To better understand the pathogenesis and health risk of this pathogen, comparative genomics and phenotypic characterization were applied to assess the pathogenicity potential of E. albertii strains isolated from wild birds in a major agricultural region in California. Shiga toxin genes stx2f were present in all avian strains. Pangenome analyses of 20 complete genomes revealed a total of 11,249 genes, of which nearly 80% were accessory genes. Both core gene-based phylogenetic and accessory gene-based relatedness analyses consistently grouped the three stx2f-positive clinical strains with the five avian strains carrying ST7971. Among the three Stx2f-converting prophage integration sites identified, ssrA was the most common one. Besides the locus of enterocyte effacement and type three secretion system, the high pathogenicity island, OI-122, and type six secretion systems were identified. Substantial strain variation in virulence gene repertoire, Shiga toxin production, and cytotoxicity were revealed. Six avian strains exhibited significantly higher cytotoxicity than that of stx2f-positive E. coli, and three of them exhibited a comparable level of cytotoxicity with that of enterohemorrhagic E. coli outbreak strains, suggesting that wild birds could serve as a reservoir of E. albertii strains with great potential to cause severe diseases in humans.

59 BASIC BIOLOGICAL SCIENCES↗

A Decade of Radial-velocity Monitoring of Vega and New Limits on the Presence of Planets

We present an analysis of 1524 spectra of Vega spanning 10 yr, in which we search for periodic radial-velocity variations. A signal with a periodicity of 0.676 day and a semi-amplitude of ∼10 m s{sup −1} is consistent with the rotation period measured over much shorter time spans by previous spectroscopic and spectropolarimetric studies, confirming the presence of surface features on this A0 star. The activity signal appears to evolve on long timescales, which may indicate the presence of failed fossil magnetic fields on Vega. TESS data reveal Vega’s photometric rotational modulation for the first time, with a total amplitude of only 10 ppm. A comparison of the spectroscopic and photometric amplitudes suggests that the surface features may be dominated by bright plages rather than dark spots. For the shortest orbital periods, transit and radial-velocity injection recovery tests exclude the presence of transiting planets larger than 2 R {sub ⊕} and most non-transiting giant planets. At long periods, we combine our radial velocities with direct imaging from the literature to produce detection limits for Vegan planets and brown dwarfs out to distances of 15 au. Finally, we detect a candidate radial-velocity signal with a period of 2.43 days and a semi-amplitude of 6 m s{sup −1}. If caused by an orbiting companion, its minimum mass would be ∼20 M {sub ⊕}; because of Vega’s pole-on orientation, this would correspond to a Jovian planet if the orbit is aligned with the stellar spin. We discuss the prospects for confirmation of this candidate planet.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hadron-Argon Cross-Section Measurements at Protodune

The Deep Underground Neutrino Experiment (DUNE) is a next-generation experiment designed to measure neutrino oscillations with unprecedented precision, specifically seeking to identify CP violation in the lepton sector. DUNE utilizes the Liquid Argon Time Projection Chamber (LArTPC) as its core detecting technology, offering high-resolution imaging to reconstruct neutrino interactions on argon. A primary challenge in this reconstruction is the modeling of hadronic final-state interactions (FSI), where hadrons produced at the initial interaction vertex undergo further scattering before escaping the nucleus. In this thesis, I utilize the ProtoDUNE-SP detector and its 1~GeV/$c$ hadron beam at the CERN Neutrino Platform to present the first measurements of total inelastic $\pi^+$-argon and proton-argon cross sections in an energy regime critical to DUNE. I detail the development and validation of the ``slicing method'', which utilizes a kiloton-scale LArTPC to extract hadron-argon cross sections. These measurements provide an indispensable benchmark for informing DUNE's FSI modeling, and our results show no significant tension with current model predictions. Additionally, I study the impact of FSI modeling on DUNE’s sensitivity to oscillation parameters. This study demonstrates that variations in FSI models may be degenerate with oscillation features, highlighting the necessity of using direct experimental data to refine interaction models. Finally, I provide an outlook on related studies that will further contribute to this effort, supporting the ambitious physics goals of upcoming neutrino oscillation programs including DUNE.

Yin-Rui, Liu [U. Chicago (main)]↗

Corresponding Standard Reference Material Data used in Partial Least Squares Regression Models for Sugar Composition Estimates in Biomass in: Economic Impact of Yield and Composition Variation in Bioenergy Crops: Populus trichocarpa

Corresponding Standard Reference Material Data used in Partial Least Squares Regression Models for Sugar Composition Estimates in Biomass in: Economic Impact of Yield and Composition Variation in Bioenergy Crops: Populus trichocarpa (for corresponding manuscript: DOI: 10.1002/bbb.2148) PDF Files: Images of 1H NMR spectra for neutralized 2-stage acid hydrolysates of 4 NIST Standard Reference Material biomass samples (Monterey Pine 8493, Sugarcane Bagasse 8491, Wheat Straw 8494, and Eastern Cottonwood/Poplar 8492) and 2 Center for Bioenergy Innovation reference biomass samples (Poplar - Populus trichocarpa and Switchgrass - Panicum Virgatum). Suppression of the water peak was achieved using a NOESY-1D with presaturation, a recycle delay of 5 s, and a total of 64 scans. Spectra were acquired at 298 K and processed with automatic phase correction, baseline correction, and chemical shift referencing to TSP-d4. Images show all 1H data from 10 to 1ppm with inset spectra of region of interest (4.0 to 3.1 ppm). Text Files: Spectra for neutralized 2-stage acid hydrolysates of 4 NIST Standard Reference Material biomass samples (Monterey Pine 8493, Sugarcane Bagasse 8491, Wheat Straw 8494, and Eastern Cottonwood/Poplar 8492) and 2 Center for Bioenergy Innovation reference biomass samples (Poplar - Populus trichocarpa and Switchgrass - Panicum Virgatum) were converted into text files for plotting. Files contain 8192 points of raw spectral data from 12.23 to -2.78 ppm. The text file contains 4 columns of data and includes: Point number, Intensity, Hz, and ppm. Xcel Spreadsheet: HPLC measured monomeric sugar concentrations and bucketed 1H NMR data used to build monomeric sugar composition prediction models. Sugar composition in biomass determined from HPLC analyses are given in mg sugar/mg of biomass. Spectral bucketing was performed using Bruker’s AMIX software. Spectra were divided into 0.005 ppm buckets in the region of 3.10– 4.15 ppm for a total of 210 buckets. Headers for the bucketed data are the chemical shift in ppm of the center of the bucket. Bucketed data was used to build partial least squares models for subsequent predictions in The Unscrambler v. 10.5(CAMO A/S, Trondheim, Norway). The formation of methanol during hydrolysis interferes with the quantitative NMR analysis of sugars, so the methanol peak centered at 3.37 ppm and spanning four buckets (3.2925 – 3.2775 ppm) was set to zero for all spectra.

09 BIOMASS FUELS↗

Enhanced critical axial tensile strain limit of CORC ® wires: FEM and analytical modeling

Abstract Conductor on Round Core (CORC ® ) cables and wires are composed of spiraled high-temperature superconducting (HTS) rare-earth barium copper oxide (REBCO) tapes, wound in multiple layers, and can carry very high currents in background magnetic fields of more than 20 T. They combine isotropic flexibility and high resilience to electromagnetic and thermal loads. The brittle nature of HTS tapes limits the maximum allowable axial tensile strain in superconducting cables. An intrinsic tensile strain above about 0.45% will introduce cracks in the REBCO layer of straight HTS tapes resulting in irreversible damage. The helical fashion at which the REBCO tapes are wound around the central core allows tapes to experience only a fraction of the total axial tensile strain applied to the CORC ® wire. As a result, the critical strain limit of CORC ® wires can be increased by a factor of more than 10 that of REBCO tapes. Finite element (FE) and analytical models are developed to predict the performance of CORC ® wires under axial tensile strain. A parametric analysis is carried out by varying the winding angle, the Poisson’s ratio of the CORC ® wire core, the core diameter, and the tape width. The results show that a small variation in winding angle can have a significant impact on the cable’s axial tensile strain tolerance. While the radial contraction of the helically wound tapes in a CORC ® wire under axial tensile strain depends on its winding angle, it is mostly driven by the Poisson’s ratio of the central core, affecting the tape strain state and thus its performance. Contact pressure from multiple layers within the CORC ® wire also affects the CORC ® wire performance. The FE model can be used to optimize the cable design for specific application conditions, resulting in an irreversible strain limit of CORC ® cables and wires as high as 7%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

ETG turbulent transport in the Mega Ampere Spherical Tokamak (MAST) pedestal

The results of a gyrokinetic analysis of turbulent transport driven by the electron temperature gradient (ETG) in the MAST pedestal are presented. Local nonlinear simulations from the gyrokinetic code, GENE, show that heat flux produced by ETG turbulence is 10%–30% of the total applied heating power in the upper pedestal and pedestal top during both the pre-edge-localized mode (ELM) (80%–99% inter-ELM period) and post-ELM (0%–20%) periods. Increasing strongly with the ETG, the ETG transport appears to be stiff. Considering radiation losses, ion thermal transport, and the strong sensitivity of the transport to the ETG, we propose that ETG transport is a plausible mechanism mediating the inter-ELM temperature profile on MAST. Cognizant of the possibility that sharp variations in the pedestal profiles may violate the local approximation, we conducted global nonlinear simulations; the results are in good agreement with local simulations except near the pedestal top, where extended radial structures and high transport levels (far beyond experimental) develop in the local simulations. We quantify and parameterize the discrepancy between local and global simulations by calculating the ratio of the radial correlation length to a length scale representative of the profile curvature. When this ratio is sufficiently small, local and global simulations agree as expected.

ETG↗

Quantification of the Radiative Effect of Aerosol–Cloud Interactions in Shallow Continental Cumulus Clouds

The indirect radiative effect of aerosol variability on shallow cumulus clouds is realized in nature with considerable concurrent meteorological variability. Large-eddy simulations constrained by observations at a continental site in Oklahoma are performed to represent the variability of different meteorological states on days with different aerosol conditions. The total radiative effect of this natural covariation between aerosol and other meteorological drivers of total cloud amount and albedo is quantified. The changes to these bulk quantities are used to understand the response of the cloud radiative effect to aerosol–cloud interactions (ACI) in the context of concurrent processes, as opposed to attempting to untangle the effect of individual processes on a case-by-case basis. Mutual information (MI) analysis suggests that meteorological variability masks the strength of the relationship between cloud drop number concentration and the cloud radiative effect. This is shown to be mostly due to variation in solar zenith angle and cloud field horizontal heterogeneity masking the relationship between cloud drop number and cloud albedo. By combining MI and more traditional differential analyses, a framework to identify important modes of covariation between aerosol, clouds, and meteorological conditions is developed. This shows that accounting for solar zenith angle variation and implementing an albedo bias correction increases the detectability of the radiative effects of ACI in simulations of shallow cumulus.

54 ENVIRONMENTAL SCIENCES↗

The AGORA High-resolution Galaxy Simulations Comparison Project. VIII. Disk Formation and Evolution of Simulated Milky Way Mass Galaxy Progenitors at 1 < z < 5

We investigate how differences in the stellar feedback produce disks with different morphologies in Milky Way–like progenitors over 1 ≤ z ≤ 5, using eight state-of-the-art cosmological hydrodynamics simulation codes in the AGORA project. In three of the participating codes, a distinct, rotation-dominated inner core emerges with a formation timescale of ≲300 Myr, largely driven by a major merger event, while two other codes exhibit similar signs of wet compaction—gaseous shrinkage into a compact starburst phase—at earlier epochs. The remaining three codes show only weak evidence of wet compaction. Consequently, we divide the simulated galaxies into two groups: those with strong compaction signatures and those with weaker ones. Galaxies in these two groups differ in size, stellar age gradients, and disk-to-total mass ratios. Specifically, codes with strong wet compaction build their outer disks in an inside-out fashion, leading to negative age gradients, whereas codes with weaker compaction feature flat or positive age gradients caused primarily by outward stellar migration. Although the stellar half-mass radii of these two groups diverge at z ∼ 3, the inclusion of dust extinction brings their sizes and shapes in mock observations closer to each other and to observed galaxies. We attribute the observed morphological differences primarily to variations in the stellar feedback implementations—such as delayed cooling timescales, and feedback strengths—that regulate both the onset and duration of compaction. Overall, our results suggest that disk assembly at high redshifts is highly sensitive to the details of the stellar feedback prescriptions in simulations.

Jung, Minyong [Seoul National Univ. (Korea, Republ↗

Processing and properties of hafnium-doped tin-modified lead zirconate titanate 95/5 ceramics

Niobium (Nb)-doped lead-tin-zirconate-titanate (PSZT) ceramics near the lead-zirconate-titanate 95/5 orthorhombic AFE-rhombohedral FE morphotropic phase boundary (PSZT 13.5/81/5.5 -1.6Nb) were prepared with up to 10 mol.% of hafnium (Hf) substituted for zirconium. The ceramics were prepared by a traditional solid-state synthesis route and sintered to near full density at 1150°C for 6 h in sealed alumina crucibles with self-same material as the lead vapor source. All compositions were ~98% dense with no detectable secondary phases by X-ray diffraction. The grain size was ~3 μm for all compositions, consisting of equiaxed grains with intergranular porosity. The compositions exhibited remnant polarization values of ~32 μC/cm 2 . Depolarization by the hydrostatic pressure-induced FE-AFE phase transition occurred at 310 MPa for all compositions, resulting in a total depolarization output of 32.4 μC/cm 2 for the PSZT ceramics. Evaluation of the R3c-R3m and R3m-Pm $\bar{3}$ m phase transition temperatures by impedance spectroscopy showed temperatures on heating ranging from 86 to 92°C and 186 to 182°C, respectively, for increasing nominal Hf content. Thermal hysteresis of the phase transitions was also observed in the ceramics, with the transition temperature on cooling being 1–4°C lower. The study demonstrated that the PSZT ceramics are relatively insensitive to variations in Hf content in the range of 0 to 10 mol.%.

36 MATERIALS SCIENCE↗

Chlorophyll- a dynamics in the lower Amazon River: insights from in situ and hyperspectral remote sensing using OCI-PACE

Chlorophyll-a concentration (Chla) is a key indicator of phytoplankton biomass and aquatic trophic status. However, satellite-derived Chla in sediment-rich waters, such as those found in the Lower Amazon River, remains challenging. The present study characterizes in situ Chla levels and their relationships with geographic, physical, and biogeochemical parameters in the Lower Amazon. Data collected between 2014 and 2017 across four hydrological seasons included measurements of Chla, remote sensing reflectance, and water quality parameters such as total suspended sediment, conductivity, water surface temperature, dissolved oxygen, pH, dissolved organic carbon and coloured dissolved organic matter. An empirical model was developed to estimate Chla using simulated hyperspectral bands from NASA’s PACE mission, achieving high performance (R 2 = 0.76; RMSE = 0.11 μg·L −1 ). Red bands proved particularly effective for Chla retrieval, while the addition of ultraviolet bands further enhanced model accuracy. The application of the developed model to satellite imagery yielded results consistent with in situ observations for the same hydrologic season. Seasonal variation and geographic location were major factors influencing Chla dynamics. This study provides a novel contribution to Chla estimation in optically complex, highly turbid waters and highlights the potential of the PACE mission to enhance global aquatic ecosystem monitoring. In conclusion, by offering freely available hyperspectral data with high radiometric resolution, PACE represents a significant advancement in the realm of remote sensing of aquatic environments.

Amazon River↗

Core-collapse supernova neutrino emission and detection informed by state-of-the-art three-dimensional numerical models

ABSTRACT Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, namely Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral standing accretion shock instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate (${\gtrsim} 50 {{\ \rm per\ cent}}$), but the time-integrated signal has a relatively modest asymmetry (${\lesssim} 20 {{\ \rm per\ cent}}$). Both features are associated with the lepton-number emission self-sustained asymmetry and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavours of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of ∼20 per cent if the distance to the CCSN is ≲6 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗

On Magnitudes and Inferred Yields of the 39 Underground Nuclear Test Explosions at the Novaya Zemlya Test Site

ABSTRACT We review the previously compiled seismic magnitudes for underground nuclear explosions at the Union of the Soviet Socialist Republics Novaya Zemlya Test Site, and analyze these data in the context of published details about individual tests and official Russian reports of annual cumulative explosive energy release at the test site. We derive relationships between teleseismic body wave magnitude mb and explosive yield Y by predicting the yield of each test using mb and minimizing the difference between predicted and reported annual cumulative yield. A linear relationship between mb and logY with a slope of 0.82 provides a good fit to the annual total yield, with a root mean square (rms) misfit of 0.095 logY units. A spectral magnitude relationship based on the explosion source model of Mueller and Murphy (1971) with a single-free-parameter provides a slightly better fit. For the 10 tests for which individual yields can be directly deduced from the official reports, the rms misfit of reported yields is less than 0.1 logY units or about 25%. More complex relationships that differentiate between the northern and southern test areas, or take account of variations in the scaled depth of burial, do not provide substantially better fits to the reported yields.

Geochemistry & Geophysics↗

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗

Local lattice distortions and electronic phases in perovskite manganite Pr 0.5 Sr 0.5 MnO 3

We use variable temperature and magnetic field total x-ray scattering to study the crystal structure of the strongly correlated Pr 0.5 Sr 0.5 MnO 3 perovskite, which is a paramagnetic insulator at room temperature, becomes a ferromagnetic metal at 272 K and, upon further decreasing the temperature, turns into an antiferromagnetic insulator at 105 K. We find that a model featuring a monoclinic symmetry captures the structure and its temperature and field evolution well, eliminating the need to evoke a phase segregation scenario as done in prior studies. It appears that coupled variations in Mn–oxygen bonding distances and angles from their values in an undistorted perovskite lattice, i.e., coupled local lattice distortions, assist the phase transitions in Pr 0.5 Sr 0.5 MnO 3 , contributing to its unique physical properties. Local structural distortions thus emerge as an important degree of freedom in strongly correlated systems, in particular perovskite manganates, and, therefore, they should be fully accounted for when their fascinating physics is considered.

36 MATERIALS SCIENCE↗

Nucleon localization function in rotating nuclei

An electron localization function was originally introduced to visualize in positional space bond structures in molecules. It became a useful tool to describe electron configurations in atoms, molecules, and solids. In nuclear physics, a nucleon localization function (NLF) has been used to characterize cluster structures in light nuclei, formation of fragments in fission, and pasta phases appearing in the inner crust of neutron stars. In the case of a deformed rotating nucleus, several NLFs can be introduced, depending on the definition of the spin-quantization axis, direction of the total angular momentum, and self-consistent symmetries of the system. Contributions to the NLF from the current density, spin-current tensor density, and density gradient terms are negligible in the nuclear interior. The oscillating pattern of the simplified NLF can be explained in terms of a constructive interference between kinetic-energy and particle densities. The characteristic nodal pattern seen in the NLF in the direction of major axis of a rotating nucleus comes from single-particle orbits carrying large aligned angular momentum. The variation of the NLF along the minor axis of the nucleus can be traced back to deformation-aligned orbits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Autoregressive distributed lag-based dynamic uniformity modeling and monitoring approaches for superconductor manufacturing

High-temperature superconductors (HTS), known for their high efficiency and low energy loss, have found profound applications across various fields, driving the demand for long, uniformly performing tapes. However, ensuring uniform performance over extended lengths of HTS tapes, often characterized by the consistency of critical current, remains challenging due to fluctuations in growth conditions during manufacturing. Here, to elucidate the mechanisms underlying variations in tape uniformity and enable real-time monitoring of associated parameters, we propose an Autoregressive Distributed Lag (ADL)-based Dynamic Uniformity Modeling and Monitoring (ADUM2) approach. This method integrates uniformity measurement, the identification of critical process parameters and real-time monitoring within the manufacturing process. The ADUM2 approach is applied to the advanced metal organic chemical vapor deposition (A-MOCVD) process, a pilot-scale method for superconductor manufacturing. Our model demonstrates superior performance compared to benchmark methods, accounting for over 80% of the total variance in the data and identifying 13 key process parameters influencing the uniformity of HTS tapes. This study offers significant insights into the high-temperature superconductor manufacturing process and holds the potential to facilitate the production of cost-effective, uniformly performing long superconducting tapes in the future.

autoregressive distributed lag analysis↗

A MultiGPU Performance-Portable Solution for Array Programming Based on Kokkos

Today, multiGPU nodes are widely used in high-performance computing and data centers. However, current programming models do not provide simple, transparent, and portable support for automatically targeting multiple GPUs within a node on application areas of array programming. In this paper, we describe a new application programming interface based on the Kokkos programming model to enable array computation on multiple GPUs in a transparent and portable way across both NVIDIA and AMD GPUs. We implement different variations of this technique to accommodate the exchange of stencils (array boundaries) among different GPU memory spaces, and we provide autotuning to select the proper number of GPUs, depending on the computational cost of the operations to be computed on arrays, that is completely transparent to the programmer. We evaluate our multiGPU extension on Summit (#5 TOP500), with six NVIDIA V100 Volta GPUs per node, and Crusher that contains identical hardware/software as Frontier (#1 TOP500), with four AMD MI250X GPUs, each with 2 Graphics Compute Dies (GCDs)for a total of 8 GCDs per node. We also compare the performance of this solution against the use of MPI + Kokkos, which is the cur-rent de facto solution for multiple GPUs in Kokkos. Our evaluation shows that the new Kokkos solution provides good scalability for many GPUs and a faster and simpler solution (from a programming productivity perspective) than MPI + Kokkos.

Valero Lara, Pedro↗