Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “I&C Modernization”

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 613 records · Page 34

Modern Aerocapture Guidance to Enable Reduced-Lift Vehicles at Neptune

Aerocapture is covered extensively in the literature as means of achieving orbital insertion with dramatic mass-saving results compared to fully-propulsive systems. One of the primary obstacles facing aerocapture is the inherent uncertainty associated with passing through a planet’s upper atmosphere. In-flight dispersions due to delivery errors, environment variables, and aerodynamic performance impose a large flight envelope. System studies for aerocapture often select high lift-to-drag ratios to compensate for these uncertainties. However, modern predictor-corrector guidance strategies have shown promise in recent years to provide robust control schemes in-situ. These algorithms do not rely on a pre-calculated reference trajectory and instead employ a numerical optimizer to continuously solve nonlinear equations of motion each guidance cycle. Numerical predictor-corrector strategies may provide considerable accuracy over heritage guidance schemes. The goal of this study is reproduce a landmark study of Neptune aerocapture and apply modern guidance to illustrate relative performance improvements and cost-saving potential. Capture constraints based on the theoretical corridor width are considered. Results indicate that heritage vehicles with moderate lift-to-drag ratios, lower than previous studies have indicated, may prove viable for aerocapture at Neptune.

Heidrich, C. R.↗

Isopropanol dehydration reaction rate kinetics measurement using H 2 O time histories

H 2 O formation during thermal decomposition of isopropanol was measured behind reflected shock waves at temperatures ranging from 1127 to 1621 K at an average pressure of 1.42 atm using a laser absorption technique. Of the five modern chemical kinetics models used to compare the H 2 O time histories, the model from Li et al (Combust. Flame 2019;207:171-185) showed the best overall agreement. Sensitivity and rate of production analyses using the Li et al model (as well as those from AramcoMech 3.0, CRECK, and Togbe et al (Energy Fuel. 2011;25:676-683)) showed unimolecular dehydration of isopropanol, iC 3 H 7 OH ⇌ H 2 O + C 3 H 6 (R1), is nearly the sole reaction controlling H2O production at early times, allowing for an a priori measurement of the forward rate constant k 1 . The Arrhenius expression k1 (s –1 ) = 2.60 × 10 13 exp(-31 120 K/T) was determined to represent best the data from this study. According to the models and previous experimental investigations, the pressure investigated is well within the high-pressure limit (HPL) for this reaction. Additional, higher-pressure experiments also confirmed the HPL assumption. An uncertainty analysis was performed by varying secondary reactions within their uncertainties and examining their effect on the overall prediction, establishing an uncertainty within ±20% for all but the highest temperature cases, which have a maximum uncertainty of ±40%. Experiments conducted with a radical trapper, toluene, showed little influence from radical chemistry, suggesting this estimated uncertainty is fairly conservative. Experimental data from Heyne et al (Z Phys Chem. 2015;229:881-907) were found to be in good agreement with the rate measurements from this study and, therefore, a second Arrhenius expression, k 1 (s –1 ) = 2.11 × 10 13 exp(-30 820 K/T), was found to represent both datasets well. This second expression has a larger temperature range of 976-1621 K. Here the present study provides the first high-temperature data collected for this reaction, adding to the limited data available for isopropanol in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An experimental and computational investigation of the structure and spectroscopic signatures of α -UO 3

α-UO 3 is a common intermediate compound found in the nuclear fuel cycle, yet the exact crystal structure of this material has long been debated. Inconsistent computational and experimental data in previous works has led to varying conclusions between authors. Furthermore, to ensure the validity of our results in this work, the structural and spectroscopic signatures of pure phase α-UO 3 are investigated using powder X-ray diffraction and optical vibrational spectroscopy (infrared and Raman). Rietveld refinement of powder X-ray diffraction data on pure phase α-UO 3 collected in this work allows us to propose an alteration to the currently accepted C2mm structure (a = 3.9705 Å, b = 6.8553 Å, c = 4.15955 Å, α = β = γ = 90°) for α-UO 3 with no uranyl [UO 2 2+ ] bonds. Raman spectra collected using two excitation wavelengths (two instruments using 532 nm and one 785 nm) are presented, and differences with recently published results are discussed. Infrared spectra from two instruments used here agree well with recently published results, but the spectral range encompassed in our data extends past what has been reported with modern techniques. Additionally, we provide tentative vibrational mode assignments based on density functional perturbation theory calculations and resulting phonon eigenvector visualizations. Unexpected features in the optical vibrational spectra of α-UO 3 are explained by unique features in the structure we present.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Hydrodynamic and Radiographic Toolbox (HART)

With a multi-lab and university team, we propose to develop new methods for a Hydrodynamic and Radiographic Toolbox (HART) that will enable a fuller and more extensive use of experimental radiographic data towards better characterizing and reducing uncertainties in predictive modeling of weapons performance. We will achieve this by leveraging recent developments in areas of computational imaging, statistical and machine learning, and reduced order modeling of hydrodynamics. In terms of software practices, by partnering with XCP (RISTRA Project) we will conform to recent XCP standards that are inline with modern software practices and standards and ensure compatibility and ease of inter-operability with existing codes. Three new activities under this proposal include (a) the development and use of deep learning-based surrogates to accelerate reconstruction and variational inference of density fields from radiographs of hydrotests, (b) a model-data fusion strategy that couples deep learning-based density reconstructions with fast hydrodynamics simulators to better constrain the reconstruction, and (c) a method for treating asymmetries using techniques adopted from limited view tomography. All three new activities will be based on improved treatment of scatter, noise, beam spot movement, detector blur, and flat fielding in the forward model, and a use of sophisticated priors to aid in the re construction. The improvements to the forward model and improved algorithmic design of the reconstruction when complete will be contained in the iterative reconstruction code SHIVA—a code project that we have recently initiated. The many ways in which machine learning can be used in the reconstruction work will be contained in a code HERMES that has been initiated with DTRA support. For example, the significant levels of acceleration that will likely be achieved by the use of machine learning techniques will permit us (and are required) to quantify uncertainties in density retrievals. Next, the two-way coupling between density reconstruction and model-based simulation of the hydrodynamics will be contained in code EREBUS, and will permit a fuller realization of the potential of the data to constrain the hydrodynamic model and better address issues related to asymmetries in the problem. Finally, we anticipate that the better consistency with physics achieved in our reconstructions will allow them to be used by X-Division more so than today.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Revisiting Selenium Interactions with Pyrite: From Adsorption to Coprecipitation

Interactions of selenium (Se), a trace element bio-essential at low concentrations but highly toxic at high concentrations, with the most abundant sulfide mineral in the Earth's crust, namely pyrite, was investigated over a wide range of time scales. At the nanosecond scale, selenate Se(VI)O 4 2– adsorption onto the net pyrite surface is shown by ab-initio computations to proceed via the formation of a chemical bond between an oxyanion oxygen atom and a surface Fe atom, weakening the other Se-O bonds and reducing Se atom oxidation state. At the hour-to-day scale, adsorption and coprecipitation of selenate Se(VI)O 4 2– and selenite, Se(IV)O 3 2– , were investigated through wet chemical batch experiments at various pH values at different sulfide concentrations. Selenium removal from solution is slower and weaker for selenate than for selenite. After 24h, only 10% of selenate, against 60% of selenite (up to 100% in the presence of sulfide), is removed by the pyrite surface. Independently of its original oxidation state, adsorbed Se is completely reduced to elemental trigonal selenium via adsorption, precipitation or coprecipitation, as shown by XANES spectroscopy. Here, our EXAFS results, compared to published data on Se-rich pyrite, show a Se to S substitution within the pyrite structure. The reductive coprecipitation mechanism of selenium with pyrite represents valuable new insights for improving our understanding of modern and ancient biogeochemical cycles involving Se. In addition, several industries can benefit from direct applications of our findings, such as water treatment, green technologies and sustainable mining.

Adsorption↗

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m 3 . This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion, Large drops, Wind tunnel↗

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m3. This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion↗

Permafrost Formation in a Meandering River Floodplain

Abstract Permafrost influences 25% of land in the Northern Hemisphere, where it stabilizes the ground beneath communities and infrastructure and sequesters carbon. However, the coevolution of permafrost, river dynamics, and vegetation in Arctic environments remains poorly understood. As rivers meander, they erode the floodplain at cutbanks and build new land through bar deposition, creating sequences of landforms with distinct formation ages. Here we mapped these sequences along the Koyukuk River floodplain, Alaska, analyzing permafrost occurrence, and landform and vegetation types. We used radiocarbon and optically stimulated luminescence (OSL) dating to develop a floodplain age map. Deposit ages ranged from modern to 10 ka, with more younger deposits near the modern channel. Permafrost rapidly reached 50% areal extent in all deposits older than 200 years then gradually increased up to ∼85% extent for deposits greater than 4 Kyr old. Permafrost extent correlated with increases in black spruce and wetland abundance, as well as increases in permafrost extent within wetland, and shrub and scrub vegetation classes. We developed an inverse model to constrain permafrost formation rate as a function of air temperature. Permafrost extent initially increased by ∼25% per century, in pace with vegetation succession, before decelerating to <10% per millennia as insulating overbank mud and moss slowly accumulated. Modern permafrost extent on the Koyukuk floodplain therefore reflects a dynamic balance between widespread, time‐varying permafrost formation and rapid, localized degradation due to cutbank erosion that might trigger a rapid loss of permafrost with climatic warming.

54 ENVIRONMENTAL SCIENCES↗

Future aerospace ground test facility requirements for the Arnold Engineering Development Center

Arnold Engineering Development Center (AEDC) was conceived at the close of World War II, when major new developments in flight technology were presaged by new aerodynamic and propulsion concepts. During the past 40 years, AEDC has played a significant part in the development of many aerospace systems. The original plans were extended through the years by some additional facilities, particularly in the area of propulsion testing. AEDC now has undertaken development of a master plan in an attempt to project requirements and to plan for ground test and computational facilities over the coming 20 to 30 years. This report was prepared in response to an AEDC request that the National Research Council (NRC) assemble a committee to prepare guidance for planning and modernizing AEDC facilities for the development and testing of future classes of aerospace systems as envisaged by the U.S. Air Force.

Kirchner, Mark E.↗

Investigating Hydrometeorological Conditions Associated with Increasing Dust Events in Southwestern United States

A recent study (Tong et al. 2017) shows a rapid intensification of dust storm activity over the southwestern United States in the past decades. For example, the frequency of windblown dust storms has increased 240 percent from 1990s to 2000s. Increasing dust storms can worsen air quality in the region. The study also finds that the intensification of dust events has a close connection with a fast-rising infectious disease (valley fever) caused by inhaling soil-dwelling fungi (Coccidioides immitis and C. posadasii) in the southwestern United States. Hydrometeorological conditions play an important role in dust storm activity, including winds, precipitation, soil moisture, atmospheric boundary stability, land surface types, etc. In this presentation, we describe our preliminary results of linking the dust trends to variations of hydrometeorological conditions at regional and global scales, using a number of Earth observations, including products from NLDAS (North American Land Data Assimilation System), TRMM (Tropical Rainfall Measuring Mission), and MERRA (Modern Era Retrospective-analysis for Research and Applications) datasets, from NASA Goddard Earth Sciences (GES) Data and Information Services Center (DISC). The GES DISC, one of the 12 NASA data centers, is home to multidisciplinary data archives such as precipitation, hydrology, atmospheric chemistry, atmospheric dynamics, etc. We demonstrate that multidisciplinary datasets and services at GES DISC can be used for interdisciplinary investigation to understand the interactions of Earth system components from an observational perspective.

satellite↗

A Generalized, Compactly-Supported Correlation Function for Data Assimilation Applications

Correlation functions play an essential role in modern data assimilation, where they are used to model covariances given a set of tunable parameters or applied as tapering functions to localize covariances in ensemble-based schemes. One of the most widely-used correlation functions in data assimilation is the Gaspari and Cohn (1999) piecewise-rational, compactly-supported parametric correlation function (hereafter referred to as GC99). The GC99 correlation function is useful due to its tunable cut-off parameter c and Gaussian-like shape achieved when the parameter a is set to one-half. These properties are attractive for tapering functions in data assimilation applications. However, the GC99 correlation function is homogeneous over Euclidean 3-space and isotropic when restricted to the sphere, properties that may be less than ideal for some geophysical applications. GC99 is also compactly-supported on a sphere of fixed radius, which requires tuning of the cut-off parameter c that can depend on the specific application. This work presents a generalization of the GC99 correlation function that allows the cut-off parameter c and shape parameter a to vary over space to gain more flexibility in shape while maintaining its compact support property. The function, which we call the Generalized Gaspari Cohn (GenGC) correlation function, introduces inhomogeneity in Euclidean 3-space and anisotropy when restricted to the sphere by allowing both parameters c and a to vary, as functions, over the spatial domain. The GC99 correlation function is a special case of GenGC where the functions c and a are held constant, as fixed parameters rather than functions. The GenGC correlation function also generalizes the follow-on to the work of Gaspari and Cohn (1999) presented in Gaspari et al. (2006), which allowed a to vary while keeping c fixed. We illustrate through simple one- and two-dimensional examples the variety of inhomogeneous and anisotropic correlation functions GenGC can produce by varying c and a over space, and suggest applications where they may be useful in data assimilation, such as covariance modeling or localization. In particular, we describe how the GenGC correlation function can be used to construct covariances using correlation length and variance fields derived from dynamics. For example, the correlation length field for advective dynamics is governed by a partial differential equation (PDE) in N spatial dimensions, where N is the number of space dimensions of the state. Correlation length fields can be determined from this PDE and used with GenGC to construct the corresponding correlations. We can then approximate the full covariance by rescaling by the variance, which also satisfies a PDE in N spatial dimensions for advective dynamics. Thus we can approximate the full covariance without solving the covariance PDE, which is in 2N spatial dimensions, by solving just two PDEs each in only N spatial dimensions. This approach to evolving the correlation length and variance fields, then reconstructing the correlations using GenGC, is suggested as an alternative to current methods of covariance modeling in data assimilation algorithms.

GC99↗

A Tracer Study on sCO2 Corrosion with Multiple Oxygen-Bearing Impurities

Abstract Several modern power production systems utilize supercritical CO 2 (sCO 2 ), which can contain O 2 and H 2 O as impurities. These impurities may degrade the compatibility of structural alloys through accelerated oxidation. However, it remains unclear which of these impurities plays a bigger role in high-temperature reactions taking place in sCO 2 . In this study, various model and commercial Fe‐ and Ni‐based alloys were exposed in 300 bar sCO 2 at 750 °C to low levels (50 ppm) of O 2 and H 2 O for 1,000 h. 18 O-enriched water was used to enable the identification of the oxygen source in the post-exposure characterization of the samples. However, oxygen from the water did not accumulate in the scale, which consisted of Cr 2 O 3 in the cases where a protective oxide formed. A 2wt.% Ti addition to a Ni-22%Cr model alloy resulted in the formation of thicker oxides in sCO 2 , while a 1wt.% Al addition reduced the scale thickness. A synergistic effect of both Al and Ti additions resulted in an even thicker oxide than what was formed solely by Ti, similar to observations for Ni-based alloy 282.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Treasuring trash: Pt/SrTiO 3 catalysts process plastic waste into high-value materials

Single-use plastics are ubiquitous throughout modern society because they have properties that make them desirable in a wide variety of applications, including low cost of production, high thermal and chemical stability, and tunable mechanical properties. Several of the recycling methods currently used to process single-use plastic waste, such as mechanical recycling and pyrolysis, do not selectively process plastic materials into uniform products, as would be advantageous for repurposing these materials in pursuit of a circular economy. Catalytic hydrogenolysis, whereby C–C bonds in polymers are broken over a supported catalyst in the presence of H 2 , is a promising approach to converting plastics into value-added products that can be sold commercially without further purification. Here, in this account, we present some recent advances in polyolefin hydrogenolysis, with a focus on Pt nanoparticles on SrTiO 3 nanocuboid supports (Pt/STO). Current developments in this field, general challenges, and future directions are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unusual Base Pair between Two 2-Thiouridines and Its Implication for Nonenzymatic RNA Copying

2-Thiouridine (s 2 U) is a nucleobase modification that confers enhanced efficiency and fidelity both on modern tRNA codon translation and on nonenzymatic and ribozyme-catalyzed RNA copying. We have discovered an unusual base pair between two 2-thiouridines that stabilizes an RNA duplex to a degree that is comparable to that of a native A:U base pair. High-resolution crystal structures indicate similar base-pairing geometry and stacking interactions in duplexes containing s 2 U:s 2 U compared to those with U:U pairs. Notably, the C=O···H–N hydrogen bond in the U:U pair is replaced with a C=S···H–N hydrogen bond in the s 2 U:s 2 U base pair. The thermodynamic stability of the s 2 U:s 2 U base pair suggested that this self-pairing might lead to an increased error frequency during nonenzymatic RNA copying. However, competition experiments show that s 2 U:s 2 U base-pairing induces only a low level of misincorporation during nonenzymatic RNA template copying because the correct A:s 2 U base pair outcompetes the slightly weaker s 2 U:s 2 U base pair. In addition, even if an s 2 U is incorrectly incorporated, the addition of the next base is greatly hindered. This strong stalling effect would further increase the effective fidelity of nonenzymatic RNA copying with s 2 U. Our findings suggest that s 2 U may enhance the rate and extent of nonenzymatic copying with only a minimal cost in fidelity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synergistically integrated phosphonated poly(pentafluorostyrene) for fuel cells

Modern electrochemical energy conversion devices require more advanced proton conductors for their broad applications. Phosphonated polymers have been proposed as anhydrous proton conductors for fuel cells. However, the anhydride formation of phosphonic acid functional groups lowers proton conductivity and this prevents the use of phosphonated polymers in fuel cell applications. Here, we report a poly(2,3,5,6-tetrafluorostyrene-4-phosphonic acid) that does not undergo anhydride formation and thus maintains protonic conductivity above 200 °C. We use the phosphonated polymer in fuel cell electrodes with an ion-pair coordinated membrane in a membrane electrode assembly. We find that this synergistically integrated fuel cell reached peak power densities of 1,130 mW cm -2 at 160 °C and 1,740 mW cm -2 at 240 °C under H2/O2 conditions, substantially outperforming polybenzimidazole- and metal phosphate-based fuel cells. Our result indicates a pathway towards using phosphonated polymers in high-performance fuel cells under hot and dry operating conditions.

25 ENERGY STORAGE↗

GISS‐E2.1: Configurations and Climatology

This paper describes the GISS‐E2.1 contribution to the Coupled Model Intercomparison Project, Phase 6 (CMIP6). This model version differs from the predecessor model (GISS‐E2) chiefly due to parameterization improvements to the atmospheric and ocean model components, while keeping atmospheric resolution the same. Model skill when compared to modern era climatologies is significantly higher than in previous versions. Additionally, updates in forcings have a material impact on the results. In particular, there have been specific improvements in representations of modes of variability (such as the Madden‐Julian Oscillation and other modes in the Pacific) and significant improvements in the simulation of the climate of the Southern Oceans, including sea ice. The effective climate sensitivity to 2xCO2 is slightly higher than previously at 2.7‐‐3.1°C (depending on version), and is a result of lower CO2 radiative forcing and stronger positive feedbacks.

General Circulation Model↗

From PeleC to PeleACC, to PeleC++

PeleC is an Exascale Computing Project application for simulating compressible combustion in complex geometries. It has been built on top of the popular AMReX library. In the beginning of the Exascale Computing Project, PeleC was focused on KNL. It uses a mixture of C++, C, and kernels written in Fortran to obtain performance by focusing on vectorization. Recently we have taken two approaches in deciding PeleC's future for obtaining performance on exascale GPU machines. In the first programming model, we decorated the Fortran kernels with OpenACC directives. This expedited our ability to run at large scales on Summit's GPUs, where we achieved a significant speedup over the CPUs on Summit. The second programming model involved rewriting the Fortran kernels in C++ and using AMReX's Kokkos-like lambda abstractions for running on the GPU. This resulted in similar speedups on Summit's GPUs over merely utilizing the CPUs. Both approaches involved AMReX's management of memory transfers between the device and host. In this work, we compare and contrast the benefits and pitfalls to both programming approaches regarding performance, performance portability, and productivity. We also discuss advantages we have found in taking the time to modernize our code and why have chosen a specific pathway to prepare our code for the future DOE exascale machines.

exascale computing↗