Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “two’s complement”

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 415 records · Page 23

A Novel Approach to Simulate Realistic Concrete Microstructures under Irradiation

The concrete biological shield of light water reactors is exposed to high neutron and gamma irradiation doses in the long term. Irradiation deteriorates the physical and mechanical properties of concrete. Such effects need to be investigated to predict the concrete’s performance in the event of a lifetime extension of a nuclear power plant. This work combines high-resolution characterization techniques with fast-Fourier transform (FFT)-based 2-D simulations to evaluate the radiation-induced volumetric expansion (RIVE) and damage in concrete microstructures under neutron irradiation. Two concrete microstructures from samples provided by the Japan Concrete Aging Management Program (JCAMP) were characterized using micro x-ray fluorescence (mXRF) to obtain elemental intensity maps, and energy-dispersive x-ray spectroscopy to complement mXRF with local elemental information for Na. Minerals and cement paste are then identified based on the elemental composition to produce high-resolution phase maps, resulting in a more accurate representation of the microstructures compared to previous work. Simulations of radiation-induced volumetric expansion (RIVE), creep, and damage in JCAMP concrete use the fast Fourier transform (FFT)-based code Microstructure Oriented Scientific Analysis of Irradiated Concrete (MOSAIC) combined with the irradiated minerals, aggregates, and concrete (IMAC) database, which contains mineral-specific RIVE models. Overall, the simulation results are in fair agreement with experimental data.

Cheniour, Amani↗

Lunar LIGO: A new concept in gravitational wave astronomy

For three decades, physicists have been in search of an elusive phenomenon predicted by Einstein's general theory of relativity; gravitational radiation. These weak vibrations of spacetime have, thus far, eluded conclusive Earth-based detection due in part to insufficient detector sensitivity and noise isolation. The detection of gravitational waves is crucial for two reasons. It would provide further evidence for the validity of Einstein's theory of relativity, the presently accepted theory of gravitation. Furthermore, the ability to identify the location of a source of a detected gravitational wave event would yield a radical new type of astronomy based on non-electromagnetic emissions. We continue our study of a lunar-based system which can provide an important complement to Earth-based analysis because it is completely independent of the geophysical sources of noise on Earth, while providing an Earth-Moon baseline for pin-pointing burst sources in the Universe. We also propose for the first time that a simplified version of the LIGO beam detector optical system, which we will call LLIGO (Lunar LIGO), could be emplaced on the Moon as part of NASA's robotic lander program now under study (Artemis). The Earth-based investigation has two major programs underway. Both involve large interferometer-type gravitational wave antennas.

Lafave, Norman↗

Pan-tropical daily L-band microwave land surface emissivity retrieval from GNSS-R observations

The uncertainties in microwave land surface emissivity (MLSE) measurements have long limited the use of spaceborne microwave radiometer data. As an emerging observation method, Global Navigation Satellite System Reflectometry (GNSS-R) has demonstrated great potential in several land and ocean applications. In this study, a method for obtaining daily MLSE dataset in the pan-tropical region from Cyclone GNSS (CYGNSS) observations is presented and evaluated. The CYGNSS observations are first aggregated into the Equal-Area-Scalable-Earth (EASE) 2.0 36 km grid by a combined weight function of distance, time, and signal-to-noise ratio variance. Then, the method employs a pixel-by-pixel regression algorithm to conduct the daily MLSE retrieval using reference emissivity derived from the Soil Moisture Active Passive (SMAP) brightness temperature. The CYGNSS MLSE shows good agreement with SMAP MLSE, delivering an overall root-mean-square error (RMSE) of 0.022 and 0.017 for horizontal and vertical polarization, respectively, during the training set spanning the whole year of 2018. Furthermore, on the test set from January 2019 to May 2019, the RMSE values amounted to 0.030 and 0.023 for horizontal and vertical polarization, respectively. Temperature records from the International Soil Moisture Network are employed to calculate the emissivity and for in-situ validation, which yield an RMSE of 0.034 and 0.026 for the two polarizations, respectively. The proposed algorithm provides an encouraging approach to obtain accurate daily MLSE dataset for microwave remote sensing. Compared to the SMAP MLSE, the CYGNSS MLSE has a remarkable improvement of 86% in temporal resolution, greatly complementing the existing microwave emissivity datasets.

54 ENVIRONMENTAL SCIENCES↗

Quantifying uncertainties due to optical potentials in one-neutron knockout reactions

One-neutron knockout reactions have been widely used to extract information about the single-particle structure of nuclei from the valley of stability to the drip lines. The interpretation of knockout data relies on reaction models, where the uncertainties are typically not accounted for. In this work, we quantify uncertainties of optical potentials used in these reaction models and propagate them, for the first time, to knockout observables using a Bayesian analysis. We study two reactions in the present paper, the first of which involves a loosely bound halo projectile, 11 Be, and the second a tightly bound projectile, 12 C. We first quantify the parametric uncertainties associated with phenomenological optical potentials. Complementing this approach, we also quantify the model uncertainties associated with the chiral forces that can be used to construct microscopic optical potentials. For the phenomenological study, we investigate the impact of the imaginary terms of the optical potential on the breakup and stripping components of the knockout cross sections as well as the impact of the angular range. For the 11 Be case, the theoretical uncertainty from the phenomenological method is on the order of the experimental uncertainty of the knockout observables; however, for the 12 C case, the theoretical uncertainty is significantly larger. Additionally, the widths of the uncertainty bands for the knockout observables obtained for the microscopic study and the phenomenological approach are of similar orders of magnitude. Based on this work we conclude that structure information inferred from the ratio of the knockout cross sections will carry a theoretical uncertainty of at least 20% for halo nuclei and at least 40% for tightly bound nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Maximum Entropy Principle in Deep Thermalization and in Hilbert-Space Ergodicity

We report universal statistical properties displayed by ensembles of pure states that naturally emerge in quantum many-body systems. Specifically, two classes of state ensembles are considered: those formed by (i) the temporal trajectory of a quantum state under unitary evolution or (ii) the quantum states of small subsystems obtained by partial, local projective measurements performed on their complements. These cases, respectively, exemplify the phenomena of “Hilbert-space ergodicity” and “deep thermalization.” In both cases, the resultant ensembles are defined by a simple principle: The distributions of pure states have maximum entropy, subject to constraints such as energy conservation, and effective constraints imposed by thermalization. We present and numerically verify quantifiable signatures of this principle by deriving explicit formulas for all statistical moments of the ensembles, proving the necessary and sufficient conditions for such universality under widely accepted assumptions, and describing their measurable consequences in experiments. We further discuss information-theoretic implications of the universality: Our ensembles have maximal information content while being maximally difficult to interrogate, establishing that generic quantum state ensembles that occur in nature hide (scramble) information as strongly as possible. Our results generalize the notions of Hilbert-space ergodicity to time-independent Hamiltonian dynamics and deep thermalization from infinite to finite effective temperature. Our work presents new perspectives to characterize and understand universal behaviors of quantum dynamics using statistical and information-theoretic tools.

Eigenstate thermalization↗

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nondetection of Helium in the Upper Atmospheres of TRAPPIST-1b, e, and f

We obtained high-resolution spectra of the ultracool M-dwarf TRAPPIST-1 during the transit of its planet “b” using two high-dispersion near-infrared spectrographs, the Infrared Doppler (IRD) instrument on the Subaru 8.2m telescope, and the Habitable Zone Planet Finder (HPF) instrument on the 10 m Hobby–Eberly Telescope. These spectroscopic observations are complemented by a photometric transit observation for planet “b” using the APO/ARCTIC, which assisted us in capturing the correct transit times for our transit spectroscopy. Using the data obtained by the new IRD and HPF observations, as well as the prior transit observations of planets “b,” “e” and “f” from IRD, we attempt to constrain the atmospheric escape of the planet using the He i triplet 10830 Å absorption line. We do not detect evidence for any primordial extended H-He atmospheres in all three planets. To limit any planet-related absorption, we place an upper limit on the equivalent widths of <7.754 mÅ for planet “b,” <10.458 mÅ for planet “e,” <4.143 mÅ for planet “f” at 95% confidence from the IRD data, and <3.467 mÅ for planet “b” at 95% confidence from HPF data. Using these limits along with a solar-like composition isothermal Parker wind model, we attempt to constrain the mass-loss rates for the three planets. For TRAPPIST-1b, our models exclude the highest possible energy-limited rate for a wind temperature <5000 K. This nondetection of extended atmospheres with low mean-molecular weights in all three planets aids in further constraining their atmospheric composition by steering the focus toward the search of high-molecular-weight species in their atmospheres.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Geos 3 Orbit determination experiment

The purpose of this experiment was to investigate the value of altimetry data in high-precision satellite orbit determination. To accomplish this, software was developed to process laser, C-band, doppler and altimeter data singly or jointly. Initially, orbit determination studies were undertaken using synthetic data to validate the software. As data became available, preliminary experiments were carried out. When all the data became available, an intensive study was made covering a 4-day span in 1976. The results showed that even with sparse altimeter data it was possible to accurately determine the semimajor axis and eccentricity with altimeter data only. When altimeter data was supplemented with (as few as) two C-band passes, high-precision ephemerides were obtained. Using two laser passes to supplement the altimetry data did not achieve that same high precision. This is probably because the geographic location (mid-Atlantic) of the highly accurate laser data were such that they did not ideally complement the available (south Atlantic and Indian Ocean) altimeter data.

Pisacane, V. L.↗

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

In the last two decades, there have been extensive developments in computational aerodynamics, which constitutes a major part of the general area of computational fluid dynamics. Such developments are essential to advance the understanding of the physics of complex flows, to complement expensive wind-tunnel tests, and to reduce the overall design cost of an aircraft, particularly in the area of aeroelasticity. Aeroelasticity plays an important role in the design and development of aircraft, particularly modern aircraft, which tend to be more flexible. Several phenomena that can be dangerous and limit the performance of an aircraft occur because of the interaction of the flow with flexible components. For example, an aircraft with highly swept wings may experience vortex-induced aeroelastic oscillations. Also, undesirable aeroelastic phenomena due to the presence and movement of shock waves occur in the transonic range. Aeroelastically critical phenomena, such as a low transonic flutter speed, have been known to occur through limited wind-tunnel tests and flight tests. Aeroelastic tests require extensive cost and risk. An aeroelastic wind-tunnel experiment is an order of magnitude more expensive than a parallel experiment involving only aerodynamics. By complementing the wind-tunnel experiments with numerical simulations the overall cost of the development of aircraft can be considerably reduced. In order to accurately compute aeroelastic phenomenon it is necessary to solve the unsteady Euler/Navier-Stokes equations simultaneously with the structural equations of motion. These equations accurately describe the flow phenomena for aeroelastic applications. At Ames a code, ENSAERO, is being developed for computing the unsteady aerodynamics and aeroelasticity of aircraft and it solves the Euler/Navier-Stokes equations. The purpose of this contract is to continue the algorithm enhancements of ENSAERO and to apply the code to complicated geometries. During the last year, the geometric capability of the code was extended to simulate transonic flows, a wing with oscillating control surface. Single-grid and zonal approaches were tested. For the zonal approach, a new interpolation technique was introduced. The key development of the algorithm was an interface treatment between moving zones for a control surface using the virtual-zone concept. The work performed during the period, 1 Apr. 1992 through 31 Mar. 1993 is summarized. Additional details on the various aspects of the study are given in the Appendices.

Obayashi, Shigeru↗

Effects of Above Real Time Training (ARTT) On Individual Skills and Contributions to Crew/Team Performance

Above Real Time Training (ARTT) is the training acquired on a real time simulator when it is modified to present events at a faster pace than normal. The experiments on training of pilots performed by NASA engineers and others have indicated that real time training (RTT) reinforced with ARTT would offer an effective training strategy for such tasks which require significant effort at time and workload management. A study was conducted to find how ARTT and RTT complement each other for training of novice pilot-navigator teams to fly on a required route. In the experiment, each of the participating pilot-navigator teams was required to conduct simulator flights on a prescribed two-legged ground track while maintaining required air speed and altitude. At any instant in a flight, the distance between the actual spatial point location of the airplane and the required spatial point was used as a measure of deviation from the required route. A smaller deviation represented better performance. Over a segment of flight or over complete flight, an average value of the deviation represented consolidated performance. The deviations were computed from the information on latitude, longitude, and altitude. In the combined ARTT and RTT program, ARTT at intermediate training intervals was beneficial in improving the real time performance of the trainees. It was observed that the team interaction between pilot and navigator resulted in maintaining high motivation and active participation throughout the training program.

Ali, Syed Firasat↗

Molecular chaperone RAP interacts with LRP1 in a dynamic bivalent mode and enhances folding of ligand-binding regions of other LDLR family receptors

The low-density lipoprotein receptor (LDLR) family of receptors are cell-surface receptors that internalize numerous ligands and play crucial role in various processes, such as lipoprotein metabolism, hemostasis, fetal development, etc. Previously, receptor-associated protein (RAP) was described as a molecular chaperone for LDLR-related protein 1 (LRP1), a prominent member of the LDLR family. We aimed to verify this role of RAP for LRP1 and two other LDLR family receptors, LDLR and vLDLR, and to investigate the mechanisms of respective interactions using a cell culture model system, purified system, and in silico modelling. Upon coexpression of RAP with clusters of the ligand-binding complement repeats (CRs) of the receptors in secreted form in insect cells culture, the isolated proteins had increased yield, enhanced folding, and improved binding properties compared with proteins expressed without RAP, as determined by circular dichroism and surface plasmon resonance. Within LRP1 CR-clusters II and IV, we identified multiple sites comprised of adjacent CR doublets, which provide alternative bivalent binding combinations with specific pairs of lysines on RAP. Mutational analysis of these lysines within each of isolated RAP D1/D2 and D3 domains having high affinity to LRP1 and of conserved tryptophans on selected CR-doublets of LRP1, as well as in silico docking of a model LRP1 CR-triplet with RAP, indicated a universal role for these residues in interaction of RAP and LRP1. Consequently, we propose a new model of RAP interaction with LDLR family receptors based on switching of the bivalent contacts between molecules over time in a dynamic mode.

59 BASIC BIOLOGICAL SCIENCES↗

Sodium azide mutagenesis induces a unique pattern of mutations

The nature and effect of mutations are of fundamental importance to the evolutionary process. The generation of mutations with mutagens has also played important roles in genetics. Applications of mutagens include dissecting the genetic basis of trait variation, inducing desirable traits in crops, and understanding the nature of genetic load. Previous studies of sodium azide-induced mutations have reported single nucleotide variants (SNVs) found in individual genes. To characterize the nature of mutations induced by sodium azide, we analyze whole-genome sequencing (WGS) of 11 barley lines derived from sodium azide mutagenesis, where all lines were selected for diminution of plant fitness owing to induced mutations. We contrast observed mutagen-induced variants with those found in standing variation in WGS of 13 barley landraces. Here, we report indels that are two orders of magnitude more abundant than expected based on nominal mutation rates. We found induced SNVs are very specific, with C → T changes occurring in a context followed by another C on the same strand (or the reverse complement). The codons most affected by the mutagen include the sodium azide-specific CC motif (or the reverse complement), resulting in a handful of amino acid changes and few stop codons. The specific nature of induced mutations suggests that mutagens could be chosen based on experimental goals. Sodium azide would not be ideal for gene knockouts but will create many missense mutations with more subtle effects on protein function.

Genetics & Heredity↗

FY 2021 Isolated Grids and Grid-Connected Turbine Reference Systems; Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL)

For individuals, businesses, and communities focused on building resilient electrical grid infrastructure, wind energy can provide an affordable, accessible, and compatible distributed energy resource option that also enhances the capabilities of local grid operations. The Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) project is a multi-year distributed wind research effort, driven through a partnership between four Department of Energy National Laboratories and industry to develop and improve the planning, design, and operation of wind-centered microgrids to complement solar, energy storage, and other distributed energy resources for grid-tied and isolated operation (U.S. Department of Energy, 2021). This report documents the application of methods developed through the initial three years of the MIRACL project to two real-world distributed wind reference systems. Specifically, the methods demonstrated in this report include 1) a market valuation framework to comprehensively value the services distributed wind can provide and 2) a resilience framework that enables stakeholders to characterize distribution system resilience and compare grid investment decisions from a resilience perspective. Additional methods mentioned in this report include distributed hybrid system design methods for grid resilience, advanced control co-simulation platforms, and power hardware-in-the-loop (PHIL) models. Preliminary results from these additional methods are presented in this report and will be demonstrated and/or applied to the reference systems in the coming year. The purpose of applying these methods to reference systems is to drive technology transfer of the theories, methodologies, and technologies developed under the MIRACL project and increase the number of referenceable case studies available to stakeholders interested in additional value-added capabilities of wind systems beyond bulk energy supply (i.e. kilowatt-hours).

17 WIND ENERGY↗

Reconnaissance and deep-drill site selection on Taylor Dome, Antarctica

Taylor Dome is a small ice dome near the head of Taylor Valley, Southern Victoria Land. The location of the dome, just west of the Transantarctic Mountains, is expected to make the composition of the accumulating snow sensitive to changes in the extent of the Ross Ice Shelf. Thus, it is linked to the discharge of the West Antarctic Ice Sheet but protected against direct influences of glacial-interglacial sea-level rise. The record of past climatic and environmental changes in the ice provides a valuable complement to the radiocarbon-dated proxy record of climate derived from perched deltas, strandlines, and moraines that have been obtained in the nearby Dry Valleys. We carried out a reconnaissance of the Taylor Dome area over the past two field seasons to determine the most favorable location to obtain a deep core to bedrock. A stake network has been established with an 80-km line roughly along the crest of Taylor Dome, and 40-km lines parallel to it and offset by 10 km. These lines have been surveyed 1990/91, and the positions of 9 grid points have been determined with geoceivers. A higher density stake network was placed and surveyed around the most likely drill area in the second year. Ground-based radar soundings in both years provided details on bedrock topography and internal layering of the ice in the drill area. An airborne radar survey in January 1992, completed the radar coverage of the Taylor Dome field area.

Grootes, Pieter M.↗

Options for flight testing rocket-based combined-cycle (RBCC) engines

While NASA's current next-generation launch vehicle research has largely focused on advanced all-rocket single-stage-to-orbit vehicles (i.e. the X-33 and it's RLV operational follow-on), some attention is being given to advanced propulsion concepts suitable for 'next-generation-and-a-half' vehicles. Rocket-based combined-cycle (RBCC) engines combining rocket and airbreathing elements are one candidate concept. Preliminary RBCC engine development was undertaken by the United States in the 1960's. However, additional ground and flight research is required to bring the engine to technological maturity. This paper presents two options for flight testing early versions of the RBCC ejector scramjet engine. The first option mounts a single RBCC engine module to the X-34 air-launched technology testbed for test flights up to about Mach 6.4. The second option links RBCC engine testing to the simultaneous development of a small-payload (220 lb.) two-stage-to-orbit operational vehicle in the Bantam payload class. This launcher/testbed concept has been dubbed the W vehicle. The W vehicle can also serve as an early ejector ramjet RBCC launcher (albeit at a lower payload). To complement current RBCC ground testing efforts, both flight test engines will use earth-storable propellants for their RBCC rocket primaries and hydrocarbon fuel for their airbreathing modes. Performance and vehicle sizing results are presented for both options.

Olds, John↗

A Study of the Structure of the Source Region of the Solar Wind in Support of a Solar Probe Mission

Despite the richness of the information about the physical properties and the structure of the solar wind provided by the Ulysses and SOHO (Solar and Heliospheric Observatory) observations, fundamental questions regarding the nature of the coronal heating mechanisms, their source, and the manifestations of the fast and slow solar wind, still remain unanswered. The last unexplored frontier to establish the connection between the structure and dynamics of the solar atmosphere, its extension into interplanetary space, and the mechanisms responsible for the evolution of the solar wind, is the corona between 1 and 30 R(sub s). A Solar Probe mission offers an unprecedented opportunity to explore this frontier. Its uniqueness stems from its trajectory in a plane perpendicular to the ecliptic which reaches within 9 R(sub s) of the solar surface over the poles and 3 - 9 R(sub s) at the equator. With a complement of simultaneous in situ and remote sensing observations, this mission is destined to detect remnants and signatures of the processes which heat the corona and accelerate the solar wind. In support of this mission, we fulfilled the following two long-term projects: (1) Study of the evolution of waves and turbulence in the solar wind (2) Exploration of signatures of physical processes and structures in the corona. A summary of the tasks achieved in support of these projects are given below. In addition, funds were provided to support the Solar Wind 9 International Conference which was held in October 1998. A brief report on the conference is also described in what follows.

Habbal, Shadia R.↗

Emissions Relationships in Western Forest Fire Plumes: I. Reducing the Effect of Mixing Errors on Emission Factors

Studies of emission factors from biomass burning using aircraft data complement the results of lab studies and extend them to conditions of immense hot conflagrations. We illustrate and discuss emission relationships for 422 individual samples from many forest-fire plumes in the Western US. The samples are from two NASA investigations: ARCTAS (Arctic Research of the Composition of the Troposphere from Aircraft and Satellites) and SEAC4RS (Studies of Emissions and Atmospheric Composition, Clouds, and Climate Coupling by Regional Surveys). This work provides sample-by-sample enhancement ratios (EnRs) for 23 gases and particulate properties. Many EnRs provide candidates for emission ratios (ERs, corresponding to the EnR at the source) when the origin and degree of transformation is understood and appropriate. From these, emission factors (EFs) can be estimated when the fuel dry mass consumed is known or can be estimated using the carbon mass budget approach. This analysis requires understanding the interplay of mixing of the plume with surrounding air. Some initial examples emphasize that measured C(tot) = CO2 + CO in a fire plume does not necessarily describe the emissions of the total carbon liberated in the flames, C(burn). Rather, it represents C(tot) = C(burn) + C(bkgd), which includes possibly varying background concentrations for entrained air. Consequently, we present a simple theoretical description for plume entrainment for multiple tracers from flame to hundreds of kilometers downwind and illustrate some intrinsic linear behaviors. The analysis suggests a Mixed Effects Regression Emission Technique (MERET), which can eliminate occasional strong biases associated with the commonly used normalized excess mixing ratio (NEMR) method. MERET splits C(tot) to reveal C(burn) by exploiting the fact that C(burn) and all tracers respond linearly to dilution, while each tracer has consistent EnR behavior (slope of tracer concentration with respect to C(burn)). The two effects are separable. Two or three or preferably more emission indicators are required as a minimum; here we used ten. Limited variations in the EnRs for each tracer can be incorporated and the variations and co-variations analyzed. The percentage CO yield (or the modified combustion efficiency) plays some role. Other co-relationships involving nitrogen and organic classes are more prominent; these have strong relationships to the C(burn) to O3 emission relationship. In summary, MERET allows fine spatial resolution (EnRs for individual observations) and comparison of similar plumes distant in time and space. Alkene ratios provide us with an approximate photochemical timescale. This allows discrimination and definition, by fire situation, of ERs, allowing us to estimate emission factors.

Chatfield, Robert B.↗

Towards Enhancing Coding Productivity for GPU Programming Using Static Graphs

The main contribution of this work is to increase the coding productivity of GPU programming by using the concept of Static Graphs. GPU capabilities have been increasing significantly in terms of performance and memory capacity. However, there are still some problems in terms of scalability and limitations to the amount of work that a GPU can perform at a time. To minimize the overhead associated with the launch of GPU kernels, as well as to maximize the use of GPU capacity, we have combined the new CUDA Graph API with the CUDA programming model (including CUDA math libraries) and the OpenACC programming model. We use as test cases two different, well-known and widely used problems in HPC and AI: the Conjugate Gradient method and the Particle Swarm Optimization. In the first test case (Conjugate Gradient) we focus on the integration of Static Graphs with CUDA. In this case, we are able to significantly outperform the NVIDIA reference code, reaching an acceleration of up to 11x thanks to a better implementation, which can benefit from the new CUDA Graph capabilities. In the second test case (Particle Swarm Optimization), we complement the OpenACC functionality with the use of CUDA Graph, achieving again accelerations of up to one order of magnitude, with average speedups ranging from 2x to 4x, and performance very close to a reference and optimized CUDA code. Our main target is to achieve a higher coding productivity model for GPU programming by using Static Graphs, which provides, in a very transparent way, a better exploitation of the GPU capacity. The combination of using Static Graphs with two of the current most important GPU programming models (CUDA and OpenACC) is able to reduce considerably the execution time w.r.t. the use of CUDA and OpenACC only, achieving accelerations of up to more than one order of magnitude. Finally, we propose an interface to incorporate the concept of Static Graphs into the OpenACC Specifications.

58 GEOSCIENCES↗