Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Results from the DOE Advanced Gas Reactor Fuel Development and Qualification Program

Modular HTGR designs were developed to provide natural safety, which prevents core damage under all design basis accidents and presently envisioned severe accidents. The principle that guides their design concepts is to passively maintain core temperatures below fission product release thresholds under all accident scenarios. This level of fuel performance and fission product retention reduces the radioactive source term by many orders of magnitude and allows potential elimination of the need for evacuation and sheltering beyond a small exclusion area. This level, however, is predicated on exceptionally high fuel fabrication quality and performance under normal operation and accident conditions. Germany produced and demonstrated high quality fuel for their pebble bed HTGRs in the 1980s, but no U.S. manufactured fuel had exhibited equivalent performance prior to the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program. The design goal of the modular HTGRs is to allow elimination of an exclusion zone and an emergency planning zone outside the plant boundary fence, typically interpreted as being about 400 meters from the reactor. To achieve this, the reactor design concepts require a level of fuel integrity that is better than that claimed for all prior US manufactured TRISO fuel, by a few orders of magnitude. The improved performance level is about a factor of three better than qualified for German TRISO fuel in the 1980’s. At the start of the AGR program, without a reactor design concept selected, the AGR fuel program selected to qualify fuel to an operating envelope that would bound both pebble bed and prismatic options. This resulted in needing a fuel form that could survive at peak fuel temperatures of 1250°C on a time-averaged basis and high burnups in the range of 150 to 200 GWd/MTHM (metric tons of heavy metal) or 16.4 to 21.8% fissions per initial metal atom (FIMA). Although Germany has demonstrated excellent performance of TRISO-coated UO2 particle fuel up to about 10% FIMA and 1150°C, UO2 fuel is known to have limitations because of CO formation and kernel migration at the high burnups, power densities, temperatures, and temperature gradients that may be encountered in the prismatic modular HTGRs. With uranium oxycarbide (UCO) fuel, the kernel composition is engineered to prevent CO formation and kernel migration, which are key threats to fuel integrity at higher burnups, temperatures, and temperature gradients. Furthermore, the recent poor fuel performance of UO2 TRISO fuel pebbles measured in Chinese irradiation testing in Russia and in German pebbles irradiated at 1250°C, and historic data on poorer fuel performance in safety testing of German pebbles that experienced burnups in excess of 10% FIMA [1] have each raised concern about the use of UO2 TRISO above 10% FIMA and 1150°C and the degree of margin available in the fuel system. This continues to be an active area of study internationally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Response approach to the squeezed-limit bispectrum: application to the correlation of quasar and Lyman-α forest power spectrum

We present the squeezed-limit bispectrum, which is generated by nonlinear gravitational evolution as well as inflationary physics, measures the correlation of three wavenumbers, in the configuration where one wavenumber is much smaller than the other two. Since the squeezed-limit bispectrum encodes the impact of a large-scale fluctuation on the small-scale power spectrum, it can be understood as how the small-scale power spectrum ``responds'' to the large-scale fluctuation. Viewed in this way, the squeezed-limit bispectrum can be calculated using the response approach even in the cases which do not submit to perturbative treatment. To illustrate this point, we apply this approach to the cross-correlation between the large-scale quasar density field and small-scale Lyman-α forest flux power spectrum. In particular, using separate universe simulations which implement changes in the large-scale density, velocity gradient, and primordial power spectrum amplitude, we measure how the Lyman-α forest flux power spectrum responds to the local, long-wavelength quasar overdensity, and equivalently their squeezed-limit bispectrum. We perform a Fisher forecast for the ability of future experiments to constrain local non-Gaussianity using the bispectrum of quasars and the Lyman-α forest. Combining with quasar and Lyman-α forest power spectra to constrain the biases, we find that for DESI the expected 1-σ constraint is err[f NL ]~60. Ability for DESI to measure f NL through this channel is limited primarily by the aliasing and instrumental noise of the Lyman-α forest flux power spectrum. Lastly, the combination of response approach and separate universe simulations provides a novel technique to explore the constraints from the squeezed-limit bispectrum between different observables.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Coast-Urban-Rural Atmospheric Gradient Experiment

Accurately representing the climate and weather variability within cities and providing well-tested representations of the impacts of urban systems on the atmospheric environment are challenges for current earth system models (ESMs). Field data for testing and developing these models are critical to addressing these limitations. The Atmospheric Radiation Measurement (ARM) user facility will support a field campaign aimed at understanding how different surface-atmosphere interactions around a city are influencing its climate. The Coast-Urban-Rural Atmospheric Gradient Experiment (CoURAGE) is expected to operate from December 2024 to November 2025 in and around Baltimore, Maryland. With its aging infrastructure, growing susceptibility to heat and flooding, and ongoing issues with air and water pollution, Baltimore is characteristic of many large industrial cities in the Eastern United States. ARM, a U.S. Department of Energy (DOE) Office of Science user facility, plans to deploy one of its three mobile observatories at Morgan State University’s Clifton Park site in downtown Baltimore. Ancillary ARM sites are also planned for rural Maryland, northwest of the city; and to the southern end of Kent Island in Chesapeake Bay. Data from CoURAGE will complement measurements from the Baltimore Social-Environmental Collaborative (BSEC), one of four DOE Urban Integrated Field Laboratories that will operate in cities across the United States to study urban climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Distributed hydrological, chemical, and microbiological measurements around Meander A of East River, Colorado.

In July-August 2018, we collected 88 filtered pore water samples at 20 cm depth and analyzed them for anions, cations, dissolved organic carbon, dissolved organic matter (DOM) quality, and standard water quality parameters. Vertical hydraulic head gradients were measured to assess the potential for upward or downward water flow and sediment grab samples were taken for grain size analyses to constrain the relationship between streambed morphology and water flow. Lastly, DNA was extracted from pore water filters for a 16S rRNA gene sequencing analysis. This allowed for a high-resolution characterization of streambed hydrology and biogeochemistry during base flow. Data.csv includes coordinates, time of sampling, grain size analysis D10 values, vertical hydraulic head gradients, microbiological diversity metrics, and multiple chemical analyses. ReadMe.txt contains all information regarding sampling procedure, analysis instrumentation, and reporting limits. Microbiological sequences analyzed here have been deposited in the NCBI Sequence Read Archive and may be accessed using accession number PRJNA527829.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

Experimental and theoretical investigation of three-dimensional turbulent boundary layers and turbulence characteristics inside an axial flow inducer passage

Analytical and experimental investigations of the characteristics of three dimensional turbulent boundary layers in a rotating helical passage of an inducer rotor are reported. Expressions are developed for the velocity profiles in the inner layer, where the viscous effects dominate, in the outer layer, where the viscous effects are small, and in the interference layer, where the end walls influence the flow. The prediction of boundary layer growth is based on the momentum integral technique. The equations derived are general enough to be valid for all turbomachinery rotors with arbitrary pressure gradients. The experimental investigations are carried out in a flat plate inducer 3 feet in diameter. The mean velocity profiles, turbulence intensities and shear stresses, wall shear stress, and limiting streamline angles are measured at various radial and chordwise locations by using rotating probes. The measurements are in general agreement with the predictions. The radial flows are well represented by an expression which includes the effect of stagger angle and radial pressure gradient. The radial flows in the rotor channel are higher than those on a single blade. The collateral region exists only very near the blade surface. The radial component of turbulence intensity is higher than the streamwise component because of the effect of rotation.

Anand, A. K.↗
DOE OSTI2018

Collapse of the electron gas from three to two dimensions in Kohn-Sham density functional theory

Here under pressure, a quasi-2D electron gas can collapse toward the true 2D limit. In this limit, the exact exchange-correlation energy per electron has a known finite limit, but general-purpose semilocal approximate density functionals, like the local density approximation (LDA) and Perdew-Burke-Ernzerhof generalized gradient approximation (PBE GGA), are known to diverge to minus infinity. Here we consider a model density for a non-interacting electron gas confined to a thickness $L$ by infinite-barrier walls, with a fixed 2D density $1/(\pi (r_s^{\text{2D}})^2)$, and $r_s^{\text{2D}}=4$ Bohr. We estimate that LDA, PBE, and the strongly constrained and appropriately normed (SCAN) meta-GGA are accurate for the exchange-correlation energy over a wide quasi-2D range, $1.5< L/r_s^{\text{2D}} < 3.85$, but not for smaller $L$. Of these functionals, only SCAN tends to a finite limit when $L$ tends to 0. Since the non-interacting kinetic energy, treated exactly in Kohn-Sham theory, dominates in this limit within a deformable jellium model, all of the general-purpose functionals can estimate the pressure required to achieve any thickness (with SCAN and LDA better than PBE). This pressure vanishes around $L/r_s^{\text{2D}} = 3.85$, where the 3D electron density is roughly that of the valence electrons in metallic potassium, and reaches about 20 GPa at $L/r_s^{\text{2D}} = 1.5$ and 400 GPa at $L/r_s^{\text{2D}} = 0.6$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1972

Rocket ascent G-limited moment-balanced optimization program (RAGMOP)

This document describes the RAGMOP (Rocket Ascent G-limited Momentbalanced Optimization Program) computer program for parametric ascent trajectory optimization. RAGMOP computes optimum polynomial-form attitude control histories, launch azimuth, engine burn-time, and gross liftoff weight for space shuttle type vehicles using a search-accelerated, gradient projection parameter optimization technique. The trajectory model available in RAGMOP includes a rotating oblate earth model, the option of input wind tables, discrete and/or continuous throttling for the purposes of limiting the thrust acceleration and/or the maximum dynamic pressure, limitation of the structural load indicators (the product of dynamic pressure with angle-of-attack and sideslip angle), and a wide selection of intermediate and terminal equality constraints.

Lyons, J. T.↗
DOE OSTI2018

Genome-centric view of carbon processing in thawing permafrost

As global temperatures rise, large amounts of carbon sequestered in permafrost are becoming available for microbial degradation. However, accurate prediction of carbon gas emissions produced from thawing permafrost is limited by our understanding of the resident microbial communities and their associated carbon metabolism. Here, metagenomic sequencing of 214 samples from intact, thawing and thawed sites collected over three years enabled the recovery of 1,529 metagenome-assembled genomes (MAGs), including many from phyla with poor genomic representation. These genomes were shown to broadly reflect the diversity of this complex ecosystem with genus-level representatives recovered for >60% of the community, constituting a two orders of magnitude increase in the number of genomes available for understanding carbon processing in this environment. Metabolic reconstruction, supported by metatranscriptomic and metaproteomic data, revealed key populations involved in organic matter degradation, including bacteria encoding a pathway for xylose degradation only previously identified in fungi. Combined analysis of the microbial communities and geochemical data highlighted lineages correlated with the production of greenhouse gases and suggest novel syntrophic relationships. Our findings link changing biogeochemistry to specific microbial lineages involved in each stage of carbon processing, providing key information for predicting the impact of climate change on these systems.

Woodcroft, Ben J.↗
DOE OSTI2017

Kink Waves in Non-isothermal Stratified Solar Waveguides: Effect of the External Magnetic Field

We study the effect of an external magnetic field on the properties of kink waves, propagating along a thin non-isothermal stratified and diverging magnetic flux tube. A wave equation, governing the propagation of kink waves under the adopted model is derived. It is shown that the vertical gradient of temperature introduces a spatially local cut-off frequency ω {sub c}. The vertical distribution of the cut-off frequency is calculated for the reference VAL-C model of the solar atmosphere and for different values of a ratio of external to internal magnetic fields. The results show that the cut-off frequency is negative below the temperature minimum due to the negative temperature gradient. In the chromosphere the cut-off frequency at a given height is smaller for a stronger external magnetic field. For the appropriate range of a ratio B{sub e} / B{sub i} ≈ 0–0.8, the cutoff lies in the range ω{sub c} ≈ 0.003–0.010 s{sup −1} (periods 600 < P{sub c} < 2000 s). The estimate of the cut-off frequency in the transition region is provided as well. In the propagating wave regime, the effective wave energy flux in the non-isothermal diverging flux tubes is the same as in the straight and homogeneous cylindrical waveguides. The obtained wave equation in the limit β = 0 is used to study the kink oscillations of non-isothermal coronal loops. It is found that the gradient of temperature along the coronal loops reduces the frequency ratio of the first overtone to the fundamental mode, i.e., ω{sub 2}/ ω{sub 1} < 2. This reduction grows for a larger ratio of temperature at the loop top to the temperature at the footpoints. Moreover, the effect of reduction is most pronounced for the steeper temperature profiles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Coherent coupling between a quantum dot and a donor in silicon

Individual donors in silicon chips are used as quantum bits with extremely low error rates. However, physical realizations have been limited to one donor because their atomic size causes fabrication challenges. Quantum dot qubits, in contrast, are highly adjustable using electrical gate voltages. This adjustability could be leveraged to deterministically couple donors to quantum dots in arrays of qubits. In this work, we demonstrate the coherent interaction of a 31P donor electron with the electron of a metal-oxide-semiconductor quantum dot. We form a logical qubit encoded in the spin singlet and triplet states of the two-electron system. We show that the donor nuclear spin drives coherent rotations between the electronic qubit states through the contact hyperfine interaction. This provides every key element for compact two-electron spin qubits requiring only a single dot and no additional magnetic field gradients, as well as a means to interact with the nuclear spin qubit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Understanding Operando Water Management in Hydroxide‐Exchange‐Membrane Fuel Cells

The water balance in hydroxide-exchange-membrane fuel cells (HEMFCs) is a key challenge for improved performance and durability, intimately linked with the various interfaces and coupled phenomena. For every 4 electrons produced, 4 water molecules are generated in the anode and 2 consumed in the cathode, while electroosmosis transports water across the HEM from the cathode to the anode. Consequently, a concentration gradient drives water back, from anode to cathode. Ineffective water management could lead to cathode dry-out, limiting reaction rate and causing ionomer degradation, or to anode flooding. To address these concerns, it is critical to measure the water transport operando . Herein, a home-built water-flux station is used to measure total water flux during cell operation with different inlet relative humidities and back pressures. Increasing the HEM thickness fourfold decreases the water flux at high current density, and utilizing microporous layers on both the anode and cathode decreases the water flux from the anode to the cathode. However, the most significant variable in changing the water flux was found by increasing the anode back pressure. Furthermore, humidity cycling significantly changed electrochemical performance without affecting the overall water fluxes. These findings can be translated to other devices utilizing an HEM.

AEMFC↗
DOE OSTI2015

Orthogonal Injection Ion Funnel Interface Providing Enhanced Performance for Selected Reaction Monitoring-Triple Quadrupole Mass Spectrometry

The electrodynamic ion funnel facilitates efficient focusing and transfer of charged particles in the higher pressure regions (e.g. ion source interfaces) of mass spectrometers, and thus providing increased sensitivity. An “off-axis” ion funnel design has been developed to reduce the source contamination and interferences from, e.g. ESI droplet residue and other poorly focused neutral or charged particles with very high mass-to charge ratios. In this study a dual ion funnel interface consisting of an orthogonal higher pressure electrodynamic ion funnel (HPIF) and an ion funnel trap combined with a triple quadruple mass spectrometer was developed and characterized. An orthogonal ion injection inlet and a repeller plate electrode was used to direct ions to an ion funnel HPIF at 9-10 Torr pressure. Several critical factors for the HPIF were characterized, including the effects of RF amplitude, DC gradient and operating pressure. Compared to the triple quadrupole standard interface more than 4-fold improvement in the limit of detection for the direct quantitative MS analysis of low abundance peptides was observed. Lastly, the sensitivity enhancement in liquid chromatography selected reaction monitoring (SRM) analyses of low abundance peptides spiked into a highly complex mixture was also compared with that obtained using a both commercial s-lens interface and a in-line dual ion funnel interface.

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

Approximate models for the ion-kinetic regime in inertial-confinement-fusion capsule implosions

“Reduced” (i.e., simplified or approximate) ion-kinetic (RIK) models in radiation-hydrodynamic simulations permit a useful description of inertial-confinement-fusion (ICF) implosions where kinetic deviations from hydrodynamic behavior are important. For implosions in or near the kinetic regime (i.e., when ion mean free paths are comparable to the capsule size), simulations using a RIK model give a detailed picture of the time- and space-dependent structure of imploding capsules, allow an assessment of the relative importance of various kinetic processes during the implosion, enable explanations of past and current observations, and permit predictions of the results of future experiments. The RIK simulation method described here uses moment-based reduced kinetic models for transport of mass, momentum, and energy by long-mean-free-path ions, a model for the decrease of fusion reactivity owing to the associated modification of the ion distribution function, and a model of hydrodynamic turbulent mixing. The transport models are based on local gradient-diffusion approximations for the transport of moments of the ion distribution functions, with coefficients to impose flux limiting or account for transport modification. After calibration against a reference set of ICF implosions spanning the hydrodynamic-to-kinetic transition, the method has useful, quantifiable predictive ability over a broad range of capsule parameter space. Calibrated RIK simulations show that an important contributor to ion species separation in ICF capsule implosions is the preferential flux of longer-mean-free-path species out of the fuel and into the shell, leaving the fuel relatively enriched in species with shorter mean free paths. Also, the transport of ion thermal energy is enhanced in the kinetic regime, causing the fuel region to have a more uniform, lower ion temperature, extending over a larger volume, than implied by clean simulations. We expect that the success of our simple approach will motivate continued theoretical research into the development of first-principles-based, comprehensive, self-consistent, yet useable models of kinetic multispecies ion behavior in ICF plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Perdew-Zunger self-interaction correction: How wrong for uniform densities and large- Z atoms?

Semilocal density functionals for the exchange-correlation energy of a many-electron system cannot be exact for all one-electron densities. In 1981, Perdew and Zunger (PZ) subtracted the fully nonlocal self-interaction error orbital-by-orbital, making the corrected functional exact for all collections of separated one-electron densities and making no correction to the exact functional. Although the PZ self-interaction correction (SIC) eliminates many errors of semilocal functionals, it is often worse for equilibrium properties of sp-bonded molecules and solids. Nonempirical semilocal functionals are usually designed to be exact for electron gases of uniform density and, thus, also make 0% error for neutral atoms in the limit of large atomic number Z, but PZ SIC is not so designed. For localized SIC orbitals, we show analytically that the local spin density approximation (LSDA)-SIC correlation energy per electron of the uniform gas in the high-density limit makes an error of -50% in the spin-unpolarized case and -100% in the fully spin-polarized case. Then we extrapolate from the Ne, Ar, Kr, and Xe atoms to estimate the relative errors of the PZ SIC exchange-correlation energies (with localized SIC orbitals) in the limit of large atomic number: about +5.5% for the LSDA-SIC and about -3.5% for nonempirical generalized gradient [Perdew-Burke-Ernzerhof (PBE)-SIC] and meta-generalized gradient strongly constrained and appropriately normed (SCAN)-SIC approximations. Finally, the SIC errors are considerably larger than those that have been estimated for LSDA-SIC by approximating the localized SIC orbitals for the uniform gas and may explain the errors of PZ SIC for equilibrium properties, opening the door to a generalized SIC that is more widely accurate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

THE BROWN DWARF KINEMATICS PROJECT (BDKP). IV. RADIAL VELOCITIES OF 85 LATE-M AND L DWARFS WITH MagE

Radial velocity measurements are presented for 85 late M- and L-type very low-mass stars and brown dwarfs obtained with the Magellan Echellette spectrograph. Targets primarily have distances within 20 pc of the Sun, with more distant sources selected for their unusual spectral energy distributions. We achieved precisions of 2–3 km s{sup −1}, and combined these with astrometric and spectrophotometric data to calculate UVW velocities. Most are members of the thin disk of the Galaxy, and velocity dispersions indicate a mean age of 5.2 ± 0.2 Gyr for sources within 20 pc. We find signficantly different kinematic ages between late-M dwarfs (4.0 ± 0.2 Gyr) and L dwarfs (6.5 ± 0.4 Gyr) in our sample that are contrary to predictions from prior simulations. This difference appears to be driven by a dispersed population of unusually blue L dwarfs which may be more prevalent in our local volume-limited sample than in deeper magnitude-limited surveys. The L dwarfs exhibit an asymmetric U velocity distribution with a net inward flow, similar to gradients recently detected in local stellar samples. Simulations incorporating brown dwarf evolution and Galactic orbital dynamics are unable to reproduce the velocity asymmetry, suggesting non-axisymmetric perturbations or two distinct L dwarf populations. We also find the L dwarfs to have a kinematic age-activity correlation similar to more massive stars. We identify several sources with low surface gravities, and two new substellar candidate members of nearby young moving groups: the astrometric binary DENIS J08230313–4912012AB, a low-probability member of the β Pictoris Moving Group; and 2MASS J15104786–2818174, a moderate-probability member of the 30–50 Myr Argus Association.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Connection between plasma response and resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression in DIII-D [Connection between plasma response and RMP ELM suppression in DIII-D]

Calculations of the plasma response to applied non-axisymmetric fields in several DIII-D discharges show that predicted displacements depend strongly on the edge current density. This result is found using both a linear two-fluid-MHD model (M3D-C1) and a nonlinear ideal-MHD model (VMEC). Furthermore, it is observed that the probability of a discharge being edge localized mode (ELM)-suppressed is most closely related to the edge current density, as opposed to the pressure gradient. It is found that discharges with a stronger kink response are closer to the peeling–ballooning stability limit in ELITE simulations and eventually cross into the unstable region, causing ELMs to reappear. Thus for effective ELM suppression, the RMP has to prevent the plasma from generating a large kink response, associated with ELM instability. Experimental observations are in agreement with the finding; discharges which have a strong kink response in the MHD simulations show ELMs or ELM mitigation during the RMP phase of the experiment, while discharges with a small kink response in the MHD simulations are fully ELM suppressed in the experiment by the applied resonant magnetic perturbation. The results are cross-checked against modeled 3D ideal MHD equilibria using the VMEC code. The procedure of constructing optimal 3D equilibria for diverted H-mode discharges using VMEC is presented. As a result, kink displacements in VMEC are found to scale with the edge current density, similar to M3D-C1, but the displacements are smaller. A direct correlation in the flux surface displacements to the bootstrap current is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anomalous thermal effect in Zr Te 5 observed via photothermal measurements

In this study, we explore the magneto-thermoelectric power (MTP) of $ZrTe$ 5 , a canonical Dirac semimetal, through a novel photothermal technique. Unlike conventional thermoelectric studies that rely on on-chip heaters and are limited by fabrication processes, especially for stress-sensitive materials, our approach utilizes photothermal effects to induce temperature gradients. Our experiments, applying a magnetic field approximately parallel and transverse to the photocurrent detection direction, reveal that the photothermal method efficiently and reliably extracts both diagonal and off-diagonal components of the thermoelectric coefficient of $ZrTe$ 5 . Here, we observe that the longitudinal MTP reproduces features previously reported in thermal transport studies, while the photoinduced transverse MTP confirms the anomalous Nernst effect. This photothermal measurement technique opens new avenues for investigating transport properties in a wide range of quantum materials, both in 3D and 2D systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

On the performance and energy efficiency of sparse linear algebra on GPUs

In this paper we unveil some performance and energy efficiency frontiers for sparse computations on GPU-based supercomputers. We compare the resource efficiency of different sparse matrix–vector products (SpMV) taken from libraries such as cuSPARSE and MAGMA for GPU and Intel’s MKL for multicore CPUs, and develop a GPU sparse matrix–matrix product (SpMM) implementation that handles the simultaneous multiplication of a sparse matrix with a set of vectors in block-wise fashion. While a typical sparse computation such as the SpMV reaches only a fraction of the peak of current GPUs, we show that the SpMM succeeds in exceeding the memory-bound limitations of the SpMV. We integrate this kernel into a GPU-accelerated Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) eigensolver. LOBPCG is chosen as a benchmark algorithm for this study as it combines an interesting mix of sparse and dense linear algebra operations that is typical for complex simulation applications, and allows for hardware-aware optimizations. In a detailed analysis we compare the performance and energy efficiency against a multi-threaded CPU counterpart. The reported performance and energy efficiency results are indicative of sparse computations on supercomputers.

Anzt, Hartwig↗