Search NASA⌕ Search

SEARCH · Search NASA

Results for “scale facto”

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 19 records

Dither Gyro Scale Factor Calibration: GOES-16 Flight Experience

This poster is a sequel to a paper presented at the 34th Annual AAS Guidance and Control Conference in 2011, which first introduced dither-based calibration of gyro scale factors. The dither approach uses very small excitations, avoiding the need to take instruments offline during gyro scale factor calibration. In 2017, the dither calibration technique was successfully used to estimate gyro scale factors on the GOES-16 satellite. On-orbit dither calibration results were compared to more traditional methods using large angle spacecraft slews about each gyro axis, requiring interruption of science. The results demonstrate that the dither technique can estimate gyro scale factors to better than 2000 ppm during normal science observations.

calibration↗

Phase-Shifting Liquid Crystal Point-Diffraction Interferometry

Microgravity fluid physics experiments frequently measure concentration and temperature. Interferometers such as the Twyman Green illustrated have performed full-field measurement of these quantities. As with most such devices, this interferometer uses a reference path that is not common with the path through the test section. Recombination of the test and reference wavefronts produces interference fringes. Unfortunately, in order to obtain stable fringes, the alignment of both the test and reference paths must be maintained to within a fraction of the wavelength of the light being used for the measurement. Otherwise, the fringes will shift and may disappear. Because these interferometers are extremely sensitive to bumping, jarring and transmitted vibration, they are typically mounted on optical isolation tables. Schlieren deflectometers or the more recent Shack-Hartmann wavefront sensors also measure concentration and temperature in laboratory fluid flows. Ray optics describe the operation of both devices. In a schlieren system, an expanded, collimated beam passes through a test section where refractive index gradients deflect rays. A lens focuses the beam to a filter placed in the rear focal plane of the decollimating lens. In a quantitative color schlieren system, gradients in the index of refraction appear as colors in the field of view due to the action of the color filter. Since sensitivity is a function of the focal length of the decollimating lens, these systems are rather long and filter fabrication and calibration is rather difficult. A Shack-Hartmann wavefront sensor is an array of small lenslets. Typical diameters are on the order of a few hundred microns. Since these lenslets divide the test section into resolution elements, the spatial resolution can be no smaller than an individual lenslet. Such a device was recently used to perform high-speed tomography of heated air exiting a 1.27 cm diameter nozzle. While these wavefront sensors are very compact, the limited spatial resolution and the methods required for data reduction suggest that a more useful instrument needs to be developed. The category of interferometers known as common path interferometers can eliminate much of the vibration sensitivity associated with traditional interferometry as described above. In these devices, division of the amplitude of the wavefront following the test section produces the reference beam. Examples of these instruments include shearing and point diffraction interferometers. In the latter case, shown schematically, a lens focuses light passing through the test section onto a small diffracting object. Such objects are typically either a circle of material on a high quality glass plate or a small sphere in a glass cell. The size of the focused spot is several times larger than the object so that the light not intercepted by the diffracting object forms the test beam while the diffracted light generates a spherical reference beam. While this configuration is mechanically stable, phase shifting one beam with respect to the other is difficult due to the common path. Phase shifting enables extremely accurate measurements of the phase of the interferogram using only gray scale intensity measurements and is the de facto standard of industry. Mercer and Creath 2 demonstrated phase shifting in a point diffraction interferometer using a spherical spacer in a liquid crystal cell as the diffracting object. By changing the voltage across the cell, they were able to shift the phase of the undiffracted beam relative to the reference beam generated by diffraction from the sphere. While they applied this technology to fluid measurements, the device shifted phase so slowly that it was not useful for studying transient phenomena. We have identified several technical problems that precluded operation of the device at video frame rates and intend to solve them to produce a phase-shifting liquid crystal point-diffraction interferometer operating at video frame rates. The first task is to produce high contrast fringes. Since the diffracted beam is much weaker than the transmitted beam, interferograms have poor contrast unless a dye is added to the liquid crystal to reduce the intensity of the undiffracted light. Dyes previously used were not rigorously characterized and suffered from hysteresis in both the initial alignment state of the device and the electro-optic switching characteristics. Hence, our initial effort will identify and characterize dyes that do not suffer from these difficulties and are readily soluble in the liquid crystal host. Since the ultimate goal of this research is to produce interferometers capable of phase shifting at video frame rates, we will quantify the difference in switching times between ferroelectric and nematic liquid crystals. While we have more experience with nematic crystals, they typically switch more slowly than ferroelectric cells. As part of that effort, we will investigate the difference in the modulation of the interferograms as a function of the type of liquid crystal in the cell. Because the temporal switching response of a liquid crystal cell is directly related its thickness, we intend to explore techniques required to produce cells that are as thin as possible. However, the cells must still produce a total phase shift of two pi radians.

Griffin, DeVon W.↗

A new interatomic potential for mixed Mg-Al-Ga-In spinels

While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. Here, in this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO 4 structure while correctly predicting that neither MgAlInO 4 nor MgGaInO 4 are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.

36 MATERIALS SCIENCE↗

The PetscSF Scalable Communication Layer

PetscSF, the communication component of the Portable, Extensible Toolkit for Scientific Computation (PETSc), is designed to provide PETSc's communication infrastructure suitable for exascale computers that utilize GPUs and other accelerators. PetscSF provides a simple application programming interface (API) for managing common communication patterns in scientific computations by using a star-forest graph representation. PetscSF supports several implementations based on MPI and NVSHMEM, whose selection is based on the characteristics of the application or the target architecture. An efficient and portable model for network and intra-node communication is essential for implementing large-scale applications. The Message Passing Interface, which has been the de facto standard for distributed memory systems, has developed into a large complex API that does not yet provide high performance on the emerging heterogeneous CPU-GPU-based exascale systems. Here, we discuss the design of PetscSF, how it can overcome some difficulties of working directly with MPI on GPUs, and we demonstrate its performance, scalability, and novel features.

97 MATHEMATICS AND COMPUTING↗

IDAES-PSE Software Tools for Optimizing Energy Systems and Market Interactions

Modern power grids coordinate electricity production and consumption via multi-scale wholesale energy markets. Historically, levelized cost metrics were the de facto standard for techno-eco-nomic analyses of energy systems and comparison of technology options. However, these metrics neglect the complexity of energy infrastructure including the time-varying value of electricity. An emerging alternative is multi-period optimization, which considers the locational marginal price of electricity as input data (parameters). In this work, we present a general interface for multi-period optimization with time-varying energy prices to facilitate rapid analysis and comparison of potential energy systems models. The PriceTakerModel class is written in the IDAES-PSE platform and allows users to generate a multi-period, price-taker model instance, as well as automatically generate common operational constraints for their model, such as start-up and shutdown. We show this interface successfully generates multi-period price-taker models, facilitates model discrimination, and aids in analyzing various technologies for deployment in unique energy markets.

Laky, Daniel↗

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Offload NIC for NASA, NLR, and Grid Computing

This work addresses distributed data management and access dynamically configurable high-speed access to data distributed and shared over wide-area high-speed network environments. An offload engine NIC (network interface card) is proposed that scales at nX10-Gbps increments through 100-Gbps full duplex. The Globus de facto standard was used in projects requiring secure, robust, high-speed bulk data transport. Novel extension mechanisms were derived that will combine these technologies for use by GridFTP, bandwidth management resources, and host CPU (central processing unit) acceleration. The result will be wire-rate encrypted Globus grid data transactions through offload for splintering, encryption, and compression. As the need for greater network bandwidth increases, there is an inherent need for faster CPUs. The best way to accelerate CPUs is through a network acceleration engine. Grid computing data transfers for the Globus tool set did not have wire-rate encryption or compression. Existing technology cannot keep pace with the greater bandwidths of backplane and network connections. Present offload engines with ports to Ethernet are 32 to 40 Gbps f-d at best. The best of ultra-high-speed offload engines use expensive ASICs (application specific integrated circuits) or NPUs (network processing units). The present state of the art also includes bonding and the use of multiple NICs that are also in the planning stages for future portability to ASICs and software to accommodate data rates at 100 Gbps. The remaining industry solutions are for carrier-grade equipment manufacturers, with costly line cards having multiples of 10-Gbps ports, or 100-Gbps ports such as CFP modules that interface to costly ASICs and related circuitry. All of the existing solutions vary in configuration based on requirements of the host, motherboard, or carriergrade equipment. The purpose of the innovation is to eliminate data bottlenecks within cluster, grid, and cloud computing systems, and to add several more capabilities while reducing space consumption and cost. Provisions were designed for interoperability with systems used in the NASA HEC (High-End Computing) program. The new acceleration engine consists of state-ofthe- art FPGA (field-programmable gate array) core IP, C, and Verilog code; novel communication protocol; and extensions to the Globus structure. The engine provides the functions of network acceleration, encryption, compression, packet-ordering, and security added to Globus grid or for cloud data transfer. This system is scalable in nX10-Gbps increments through 100-Gbps f-d. It can be interfaced to industry-standard system-side or network-side devices or core IP in increments of 10 GigE, scaling to provide IEEE 40/100 GigE compliance.

Awrach, James↗

Parametric study of statistical bias in laser Doppler velocimetry

Analytical studies have often assumed that LDV velocity bias depends on turbulence intensity in conjunction with one or more characteristic time scales, such as the time between validated signals, the time between data samples, and the integral turbulence time-scale. These parameters are presently varied independently, in an effort to quantify the biasing effect. Neither of the post facto correction methods employed is entirely accurate. The mean velocity bias error is found to be nearly independent of data validation rate.

Gould, Richard D.↗

Testing and Expertise for Marine Energy (TEAMER) Program Support - Numerical Modeling Assistance for iProTech's "PIP" WEC Device: Cooperative Research and Development Final Report, CRADA Number CRD-20-17303

Develop a time-domain, 3 degree of freedom numerical model of IProTech's PIP device, including power take-off (PTO), that can be used to develop control algorithms and establish a baseline technology performance level (TPL). This is important to IProTech because having a good numerical model allows IProTech to optimize design and control algorithms on an ongoing basis. A baseline TPL enables the PIP technology to be benchmarked against other wave energy converter (WEC) concepts, providing (i) justification (or not) for investment in wave tank model testing, (ii) identifying optimum design parameters and control algorithms for the test model, (iii) specifying test conditions and (iv) interpreting wave tank test results. An accredited WEC-Sim model of the PIP device is an essential tool for further development of the technology. 1) The development of numerical models of the PIP concept will enable iProTech to evaluate the potential performance of this concept, and to investigate how sensitive the device is to different variables within the system (e.g. hydraulic components, geometry dimensions, etc.) 2) The numerical models were developed with WEC-Sim and PTO-Sim -- the de-facto industry standard WEC numerical modeling tools. A detailed model of the system's hydraulic circuit was developed with PTO-Sim. However, WEC-Sim is based on linear hydrodynamics and has several limitations; validation against high fidelity models/physical scale models would help to build confidence in the results.

16 TIDAL AND WAVE POWER↗

Formation of 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 during Thermochemical Reduction of 12R-Ba 4 CeMn 3 O 12 : Identification of a Polytype in the Ba(Ce,Mn)O 3 Family

The resurgence of interest in a hydrogen economy and the development of hydrogen-related technologies has initiated numerous research and development efforts aimed at making the generation, storage, and transportation of hydrogen more efficient and affordable. Solar thermochemical hydrogen production (STCH) is a process that potentially exhibits numerous benefits such as high reaction efficiencies, tunable thermodynamics, and continued performance over extended cycling. Although CeO 2 has been the de facto standard STCH material for many years, more recently 12R-Ba 4 CeMn 3 O 12 (BCM) has demonstrated enhanced hydrogen production at intermediate H 2 /H 2 O conditions compared to CeO 2 , making it a contender for large-scale hydrogen production. However, the thermo-reduction stability of 12R-BCM dictates the oxygen partial pressure (pO 2 ) and temperature conditions optimal for cycling. In this study, we identify the formation of a 6H-BCM polytype at high temperature and reducing conditions, experimentally and computationally, as a mechanism and pathway for 12R-BCM decomposition. 12R-BCM was synthesized with high purity and then controllably reduced using thermogravimetric analysis (TGA). Synchrotron X-ray diffraction (XRD) data is used to identify the formation of a 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 (6H-BCM) polytype that is formed at 1350 degrees C under strongly reducing pO 2 . Density functional theory (DFT) total energy and defect calculations show a window of thermodynamic stability for the 6H-polytype consistent with the XRD results. These data provide the first evidence of the 6H-BCM polytype and could provide a mechanistic explanation for the superior water-splitting behaviors of 12R-BCM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗

Reducing the Cost of CCSD Basis Set Extrapolation in Ab Initio Computational Thermochemistry

Here, a series of approximations to CCSD contributions in computational model chemistries is presented in the context of kcal mol –1 , kJ mol –1 , and 20 cm –1 theoretical predictions of total atomization energies, benchmarked within the HEAT+CH 4 test suite. A specific set of circumstances where MP2, without empirical scaling, may be used as an effective intermediate in the first two of these accuracy ranges was determined. However, SDQ-MP4, a method long used in pursuit of kcal mol –1 accuracy but relatively unstudied in the subchemical accuracy community, offers significant improvement over the quality of MP2 as a basis-set intermediate at significantly reduced cost compared to CCSD. Given this, we argue for SDQ-MP4 as the de facto CCSD basis-set intermediate in sub-chemical accuracy calculations when CCSD in a desired basis set becomes unaffordable. We additionally report on a “CBS-like” scheme, where MP2 and SDQ-MP4 are used in conjunction to create a “cheap” three-part approximation of large CCSD basis set limits. The data for the CCSD approximation schemes are organized in such a way that model chemistry developers can locate an analog of their current approach for the CCSD basis set limit and explore alternative intermediates that either decrease computational cost or increase computational accuracy. We also show, for a handful of molecules, that SDQ-MP4 shows promise as an effective basis-set intermediate for harmonic and fundamental frequency computations, allowing for zero-point corrections of nearly CCSD(T)/ANO1 quality using simple composite methods that only require CCSD(T)/ANO0.

Thorpe, James H. [Argonne National Laboratory (ANL↗

Patch2Self2: Self-supervised Denoising on Coresets via Matrix Sketching

Diffusion MRI (dMRI) non-invasively maps brain white matter yet necessitates denoising due to low signal-to-noise ratios. Patch2Self (P2S) employing self-supervised techniques and regression on a Casorati matrix effectively denoises dMRI images and has become the new de-facto standard in this field. P2S however is resource intensive both in terms of running time and memory usage as it uses all voxels (n) from all-but-one held-in volumes (d-1) to learn a linear mapping Phi : \mathbb R ^ n x(d-1) \mapsto \mathbb R ^ n for denoising the held-out volume. The increasing size and dimensionality of higher resolution dMRI acquisitions can make P2S infeasible for large-scale analyses. This work exploits the redundancy imposed by P2S to alleviate its performance issues and inspect regions that influence the noise disproportionately. Specifically this study makes a three-fold contribution: (1) We present Patch2Self2 (P2S2) a method that uses matrix sketching to perform self-supervised denoising. By solving a sub-problem on a smaller sub-space so called coreset we show how P2S2 can yield a significant speedup in training time while using less memory. (2) We present a theoretical analysis of P2S2 focusing on determining the optimal sketch size through rank estimation a key step in achieving a balance between denoising accuracy and computational efficiency. (3) We show how the so-called statistical leverage scores can be used to interpret the denoising of dMRI data a process that was traditionally treated as a black-box. Experimental results on both simulated and real data affirm that P2S2 maintains denoising quality while significantly enhancing speed and memory efficiency achieved by training on a reduced data subset.

Fadnavis, Shreyas↗

DeepUQ: Assessing the Aleatoric Uncertainties from two Deep Learning Methods

Assessing the quality of aleatoric uncertainty estimates from uncertainty quantification (UQ) deep learning methods is important in scientific contexts, where uncertainty is physically meaningful and important to characterize and interpret exactly. We systematically compare aleatoric uncertainty measured by two UQ techniques, Deep Ensembles (DE) and Deep Evidential Regression (DER). Our method focuses on both zero-dimensional (0D) and two-dimensional (2D) data, to explore how the UQ methods function for different data dimensionalities. We investigate uncertainty injected on the input and output variables and include a method to propagate uncertainty in the case of input uncertainty so that we can compare the predicted aleatoric uncertainty to the known values. We experiment with three levels of noise. The aleatoric uncertainty predicted across all models and experiments scales with the injected noise level. However, the predicted uncertainty is miscalibrated to $\rm{std}(\sigma_{\rm al})$ with the true uncertainty for half of the DE experiments and almost all of the DER experiments. The predicted uncertainty is the least accurate for both UQ methods for the 2D input uncertainty experiment and the high-noise level. While these results do not apply to more complex data, they highlight that further research on post-facto calibration for these methods would be beneficial, particularly for high-noise and high-dimensional settings.

Nevin, Rebecca↗

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗

What is flat ΛCDM, and may we choose it?

The Universe is neither homogeneous nor isotropic, but it is close enough that we can reasonably approximate it as such on suitably large scales. The inflationary-Λ-Cold Dark Matter (ΛCDM) concordance cosmology builds on these assumptions to describe the origin and evolution of fluctuations. With standard assumptions about stress-energy sources, this system is specified by just seven phenomenological parameters, whose precise relations to underlying fundamental theories are complicated and may depend on details of those fields. Nevertheless, it is common practice to set the parameter that characterizes the spatial curvature, Ω K , exactly to zero. This parameter-fixed ΛCDM is awarded distinguished status as separate model, "flat ΛCDM." Ipso facto this places the onus on proponents of "curved ΛCDM" to present sufficient evidence that Ω K ≠ 0, and is needed as a parameter. While certain inflationary model Lagrangians, with certain values of their parameters, and certain initial conditions, will lead to a present-day universe well-described as containing zero curvature, this does not justify distinguishing that subset of Lagrangians, parameters and initial conditions into a separate model. Absent any theoretical arguments, we cannot use observations that suggest small Ω K to enforce Ω K = 0. Our track record in picking inflationary models and their parameters a priori makes such a choice dubious, and concerns about tensions in cosmological parameters and large-angle cosmic-microwave-background anomalies strengthens arguments against this choice. We argue that Ω K must not be set to zero, and that ΛCDM remains a phenomenological model with at least 7 parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization of the Aerosol Source Term in Dry Storage Canisters

As long-term dry storage of used nuclear fuel at independent spent fuel storage installations (ISFSI) trends toward the de facto back end of the US fuel cycle, it becomes appropriate to investigate potential degradation and dispersion scenarios for suitable risk mitigation purposes. Pitting and subsequent stress corrosion cracking of the canister wall is currently viewed as a potential scenario leading to a through-wall pathway for contamination to be transferred from within the storage container to the surrounding environment. While stress corrosion cracking measurements are currently underway to further characterize this scenario, a parallel effort endeavors to perform a consequence analysis of conditions in which through wall cracks are indeed formed. This effort consists of engineering scale modeling using the GOTHIC and MELCOR software packages along with experimental depletion and penetration tests.

spent fuel storage, aerosol deposition, Stress Cor↗

Mixed-precision iterative refinement using tensor cores on GPUs to accelerate solution of linear systems

Double-precision floating-point arithmetic (FP64) has been the de facto standard for engineering and scientific simulations for several decades. Problem complexity and the sheer volume of data coming from various instruments and sensors motivate researchers to mix and match various approaches to optimize compute resources, including different levels of floating-point precision. In recent years, machine learning has motivated hardware support for half-precision floating-point arithmetic. A primary challenge in high-performance computing is to leverage reduced-precision and mixed-precision hardware. We show how the FP16/FP32 Tensor Cores on NVIDIA GPUs can be exploited to accelerate the solution of linear systems of equations Ax = b without sacrificing numerical stability. The techniques we employ include multiprecision LU factorization, the preconditioned generalized minimal residual algorithm (GMRES), and scaling and auto-adaptive rounding to avoid overflow. We also show how to efficiently handle systems with multiple right-hand sides. On the NVIDIA Quadro GV100 (Volta) GPU, we achieve a 4×-5× performance increase and 5× better energy efficiency versus the standard FP64 implementation while maintaining an FP64 level of numerical stability.

GMRES↗