Search NASA⌕ Search

SEARCH · Search NASA

Results for “reach code”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Clocking the formation of today’s largest galaxies: wide field integral spectroscopy of brightest cluster galaxies and their surroundings

ABSTRACT The formation and evolution of local brightest cluster galaxies (BCGs) is investigated by determining the stellar populations and dynamics from the galaxy core, through the outskirts and into the intracluster light (ICL). Integral spectroscopy of 23 BCGs observed out to $4\, r_{e}$ is collected and high signal-to-noise regions are identified. Stellar population synthesis codes are used to determine the age, metallicity, velocity, and velocity dispersion of stars within each region. The ICL spectra are best modelled with populations that are younger and less metal-rich than those of the BCG cores. The average BCG core age of the sample is $\rm 13.3\pm 2.8\, Gyr$ and the average metallicity is $\rm [Fe/H] = 0.30\pm 0.09$, whereas for the ICL the average age is $\rm 9.2\pm 3.5\, Gyr$ and the average metallicity is $\rm [Fe/H] = 0.18\pm 0.16$. The velocity dispersion profile is seen to be rising or flat in most of the sample (17/23), and those with rising values reach the value of the host cluster’s velocity dispersion in several cases. The most extended BCGs are closest to the peak of the cluster’s X-ray luminosity. The results are consistent with the idea that the BCG cores and inner regions formed quickly and long ago, with the outer regions and ICL forming more recently, and continuing to assemble through minor merging. Any recent star formation in the BCGs is a minor component, and is associated with the cluster cool core status.

Edwards, Louise O. V.↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
NASA NTRS1995

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗
DOE OSTI2024

Elucidating the corrosion mechanism of commercial Ni-based Superalloys in UCl3 containing-chloride Molten Salt Systems

Molten salt reactors (MSRs) have gained renewed interest, providing several advantages over their predecessors, including the capability to consume spent fuels, enhancing the environmental sustainability of the uranium fuel cycle. For example, molten chloride fast reactors (MCFRs) can reach criticality with molten chloride spent fuel containing high concentrations of impurities, such as actinide products like uranium chloride (UCl3). However, the redox potential of chloride molten salt fuels may change in the presence of these impurities, dictating their corrosivity and in turn the corrosion performance of structural components, such as those constructed from nickel (Ni)-based alloys. The purpose of this investigation is to assess the extent of corrosion of Ni-based alloy, Inconel 617, when exposed to UCl3-LiCl-KCl eutectic salt. Inconel 617 one of only six structural materials that are fully qualified by the American Society for Mechanical Engineers (ASME) Boiler and Pressure Vessel Code for high-temperature nuclear reactor components, making it a technologically mature material to consider for constructing MCFRs. Inconel 617 specimens were submerged in a static LiCl-KCl-UCl3 eutectic salt mixture heated at 700 C for 1000 h under an inert atmosphere. Upon completion, the extent of corrosion was analyzed through a multi-modal characterization approach spanning the engineering to nanoscale, employing computed tomography, focused-ion beam, and transmission electron microscopy techniques. Results from this investigation will enhance our understanding of property-to-performance relationships of candidate structural materials for MSRs with respect to corrosion resistance and interactions between the salt and alloy interface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microscale THMC Modeling of Pressure Solution in Salt Rock: Impacts of Geometry and Temperature

Abstract Pressure solution, a mechanism that involves tight coupling between the geometry and thermal-hydro-mechanical-chemical (THMC) processes, plays an important role in diagenesis. In this study, we make the first attempt to conduct microscale THMC modeling to understand and quantify the impacts of geometry and temperature on pressure solution, taking natural salt rock as an example. This modeling capability is achieved by expanding a novel MC code that we developed previously (Hu et al. J Geophys Res: Solid Earth 126:e2021JB023112, 2021) to include temperature effects. We first conduct a simulation of an example that involves a single brine inclusion within a single halite grain and find that the temperature impact is limited for that case. We then extract geometry from an image of a natural salt rock and conduct simulations with different cases: (A) only temperature and no stress, (B) only stress and no temperature, and (C) with both stress and temperature. These different cases result in quite different phenomena. In case A, dissolution and precipitation occur across the entire system due to isolated pore space reaching a localized mass balance between dissolution, precipitation, and diffusion. In case B, intense geometric features (e.g., major asperities, inclusions) in one area undergo stress concentration, thus dominating pressure solution in that area. In case C, pressure solution is spread out at contacting highly stressed geometric features close to the hotter side. We conclude that geometric features dominate stress distribution, thus dominating pressure solution in a natural salt rock that may be affected by the temperature if a sufficient temperature gradient is applied.

58 GEOSCIENCES↗
DOE OSTI2017

Triggering Collapse of the Presolar Dense Cloud Core and Injecting Short-lived Radioisotopes with a Shock Wave. V. Nonisothermal Collapse Regime

Recent meteoritical analyses support an initial abundance of the short-lived radioisotope (SLRI) {sup 60}Fe that may be high enough to require nucleosynthesis in a core-collapse supernova, followed by rapid incorporation into primitive meteoritical components, rather than a scenario where such isotopes were inherited from a well-mixed region of a giant molecular cloud polluted by a variety of supernovae remnants and massive star winds. This paper continues to explore the former scenario, by calculating three-dimensional, adaptive mesh refinement, hydrodynamical code (FLASH 2.5) models of the self-gravitational, dynamical collapse of a molecular cloud core that has been struck by a thin shock front with a speed of 40 km s{sup −1}, leading to the injection of shock front matter into the collapsing cloud through the formation of Rayleigh–Taylor fingers at the shock–cloud intersection. These models extend the previous work into the nonisothermal collapse regime using a polytropic approximation to represent compressional heating in the optically thick protostar. The models show that the injection efficiencies of shock front materials are enhanced compared to previous models, which were not carried into the nonisothermal regime, and so did not reach such high densities. The new models, combined with the recent estimates of initial {sup 60}Fe abundances, imply that the supernova triggering and injection scenario remains a plausible explanation for the origin of the SLRIs involved in the formation of our solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The HUNTER Dynamic Human Reliability Analysis Tool: Coupling an External Plant Code

The Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) is a framework to support dynamic human reliability analysis (HRA) with the aim to develop standalone software to perform the dynamic HRA calculations. Within the HRA, human actions in nuclear power plants (NPPs) are predicated by plant states, and human actions influence the plant. In other words, plant operations are necessarily recursive, and it becomes challenging to model complex human-plant interactions. Consequently, we have linked two software simulations that complement those shortcomings. RELAP5-3D—the Reactor Excursion and Leak Analysis Program (RELAP; Aumiller, Tomlinson, and Bauer 2001) is the foundational thermal-hydraulic software used to model nuclear systems. Using RELAP5-3D, we have simulated the plant operations proceeding according to procedures developed to address emergent situations in NPPs. Plant operations include various actions such as the operator checking plant parameters, as well as actions that are continuously performed over time until a specific parameter reaches certain criteria. That means that HUNTER and RELAP5-3D exchange information with each other and should be carried out simultaneously over time. To simulate plant operations, which represent the actual operator checks of plant parameters and corresponding manual control actions, changes in plant status are identified through simulation and performed according to the criteria and order of the procedure. Thus, the goal of coupling HUNTER with RELAP5-3D is to facilitate synchronous coupling, where human and plant models provide iterative feedback loops that drive the course of actions. The advantage of coupling with RELAP5-3D to serve as the external environment module in HUNTER is the ability to customize the plant model and streamline for particular model applications. In this paper, we will address the key features of the coupling and the coupling structure built to perform the feedback loops.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

A WISER way for simulating the performance of gratings

Soft X-ray monochromators for synchrotron radiation sources have been continuously developed over the years, improving energy resolution and stability. Great effort has been made in improving the surface quality of the optics involved, reaching values permitting diffraction-limited images. Still, one problem has not been solved, nor fully understood, yet: groove placing errors. Nowadays, these are one of the major factors responsible for resolving the power reduction of diffraction-gratings-based X-ray monochromators. Despite decades of use of gratings, there is not yet a universally established method for predicting and simulating the effect of groove placing errors on monochromator performance. This is especially important in the new scenario of high-coherent X-ray sources, i.e. diffraction-limited storage rings and free-electron lasers. To address this problem, in this article an approach based on WISER (Wavefront propagatIon Simulation codE libRary) is presented. WISER is a physical optics simulation package, also available in the user-friendly Orange Synchrotron Radiation Suite – OASYS. Even though it was originally conceived to assess the focusing performance of X-ray mirrors in the presence of height defects, it perfectly simulates the performance of a periodic (or quasi-periodic) structure like a diffraction grating. In this article, the way to use WISER and its application to a specific case, e.g. the design of a monochromator for the upgrade of a beamline at the Advanced Light Source, are shown. A simple rule for estimating how well the grooves are placed on a grating, based on calculation of the Strehl ratio, is also presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1982

A numerical simulation of the magnetospheric gate model for the X-ray bursters

A Lagrangian, stellar interior computer code which assumes radiation transport by diffusion is used to numerically simulate an X-ray burst occurring on a rapidly rotating, 1 solar mass neutron star accreting material from a companion in noncontinuous, episodic fashion. The blob of gas from the companion extends from the surface of the neutron star to a radius of 185 km, is optically thick, has mass of 10 to the -11th solar masses, and is initially at rest with respect to the surface. The gas initially has a temperature of about 10 million K and a surface luminosity of 0.1 that of the sun. The infall results in a burst which lasts about 0.1 second and reaches a peak luminosity and effective temperature of 240,000 suns and nine million K. The burst is followed by a phase of oscillations with a period of 0.2 second. It is found that the burst is too cool and rapid to resemble normal bursts, although the time scale is in agreement with some observations.

Starrfield, S.↗
DOE OSTI2020

Collective mode Brownian dynamics: A method for fast relaxation of statistical ensembles

Sampling equilibrium configurations of correlated systems of particles with long relaxation times (e.g., polymeric solutions) using conventional molecular dynamics and Monte Carlo methods can be challenging. This is especially true for systems with complicated, extended bond network topologies and other interactions that make the use and design of specialized relaxation protocols infeasible. We introduce a method based on Brownian dynamics simulations that can reduce the computational time it takes to reach equilibrium and draw decorrelated samples. Importantly, the method is completely agnostic to the particle configuration and the specifics of interparticle forces. In particular, we develop a mobility matrix that excites non-local, collective motion of N particles and can be computed efficiently in $\mathcal{O}(N)$ time. Particle motion in this scheme is computed by integrating the overdamped Langevin equation with an Euler–Maruyama scheme, in which Brownian displacements are drawn efficiently using a low-rank representation of the mobility matrix in position and wave space. We demonstrate the efficacy of the method with various examples from the realm of soft condensed matter and release a massively parallel implementation of the code as a plugin for the open-source package HOOMD-blue [J. A. Anderson et al., J. Comput. Phys. 227, 5342 (2008) and J. Glaser et al., Comput. Phys. Commun. 192, 97 (2015)] which runs on graphics processing units

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Dynamical Aspects of the Quark Gluon Plasma

Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far-reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons is exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges is being developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients are calculated across beam energies. New techniques are being developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project provides essential guidance to the Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Parthenon—a performance portable block-structured adaptive mesh refinement framework

On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multidimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of 1.7 × 10 13 zone-cycles/s on 9216 nodes (73,728 logical GPUs) at [Formula: see text] weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Transfer Request from Rutgers University to University of Illinois Urbana Champaign

Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons will be exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges must be developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients will be calculated across beam energies. New techniques will be developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project will provide essential guidance to the upcoming Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Report on Assessment of Advanced Reactor Fuel Performance

Three activities have been performed in FY 2019 to support the Advanced Fuel Cycle Campaign, which include assessment of annular fuel performance in an ultra-high burnup sodium-cooled reactor (SFR), modeling of annular fuel behavior during a transient scenario, and feasibility of a sustainable closed fuel cycle utilizing melt-refining process. Annular fuel performance parameters, including temperature profile, swelling strain, and cumulative damage fraction (CDF), in a 3000 MWt Sustainable Sodium-cooled Fast Reactor (SSFR) were assessed using the BISON fuel performance code. The SSFR is once-through ultra-long-life core based on a 34-batch fuel management scheme with 1.5-year cycle length. The fuel residence time in the SSFR is 51 years and the average discharged burnup is about 30%. Annular fuel has potential to irradiate for ultra-high burnup without major fuel-cladding mechanical interaction. This feature is achieved by using a low smeared density (~55%). In the SSFR design, however, 75% smeared density and a thick fuel pin were preferred to maintain criticality for 51 years. Due to the high smeared density and thick fuel pin configuration, the annular fuel was predicted to suffer premature cladding failure before reaching the discharge burnup of the SSFR. The current fuel creep model predicts a prominent creep strain due to gravity over the half-century irradiation, which may significantly affect the evolution of CDF. Thus, in order to use the annular fuel in an ultra-high burnup core such as the SSFR, a low smeared density and thin annular fuel is strongly recommended.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Thermal energy mitigation and toroidal peaking effects in JET disruptions

Previous investigations on JET suggest half or less of plasma stored thermal energy W th is radiated (f rad,th ≲ 0.5) using either massive gas injection (MGI) or shattered pellet injection (SPI) disruption mitigation. We investigate whether the apparent incomplete f rad,th is explained by radiation peaking near the injection plume. High toroidal peaking throughout the pre-thermal quench is found in argon–deuterium MGI on JET, with typically >3× higher radiation near the injector than toroidally distant. Previously unexplained toroidal bolometry measurements in neon–deuterium SPI are reproduced with similar peaking using the Emis3D radiation analysis code. These observations align with results from Alcator C-Mod and KSTAR. This peaking is not captured by previous JET studies that found poor thermal mitigation. Two sets of neon–deuterium SPI and two sets of argon–deuterium MGI are analyzed using Emis3D. In SPI, f rad,th rises from no-plume estimates of 0.31 and 0.66 to lower bounds of 0.84 and 0.92, respectively, and f rad,th ∼ 1 is possible. In MGI, the toroidal spread of the peaking feature is poorly constrained. f rad,th up to 0.85 and 0.65 are possible using the largest possible spread, increasing from 0.42 and 0.28, although f rad,th ∼ 1 does not appear to be reached. Revised mitigation estimates on JET suggest a lower melt risk to the divertor in mitigated disruptions on ITER and SPARC than previously thought. However, peaking near injectors could increase flash melting risk on nearby plasma facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2026

Benchmarking Bayesian Optimization Frameworks and Acquisition Strategies for Materials Discovery and Autonomous Laboratories

Bayesian optimization (BO) can accelerate materials discovery by guiding expensive experiments toward the most promising processing conditions. We systematically compare five BO surrogate and framework combinations (Gaussian processes in Ax, Gaussian processes and Monte-Carlo neural networks in BayBE, random forests in Lolopy, and tree-structured Parzen (TPE) estimators in Hyperopt) on three benchmarks that mimic common materials design tasks (a discrete solid-electrolyte composition space, a hybrid discrete/continuous laminate-composite design problem solved with micromechanics modeling, and the continuous Ishigami analytic function which is a standard optimization benchmark). Each BO surrogate is paired with posterior mean, probability of improvement, and expected improvement acquisition functions and run for 100 trials from randomized initial samples with uniform random search providing a control. Across five random seeds per setting, BayBE’s Gaussian-process surrogate with expected improvement consistently reached ≥95 % of the known optimum in the fewest evaluations, while Lolopy’s random forest matched or exceeded GP performance on purely categorical or mixed spaces at a higher computational cost. Posterior mean alone often stagnated at local optima, underscoring the need for exploration, whereas probability and expected improvement balanced exploration and exploitation leading to better optimization in fewer trials. Execution times ranged from milliseconds for TPE to minutes for neural-network and random-forest surrogates. These results establish baseline expectations for BO in automated materials laboratories and highlight expected improvement with Gaussian processes as a reliable first choice, with random forests offering a strong alternative when categorical variables dominate. The benchmark suite and code are released to facilitate future surrogate, acquisition, and constraint-handling research in data-driven materials optimization.

Bayesian optimization↗
DOE OSTI2022

Portability: A Necessary Approach for Future Scientific Software

Today's world of scientific software for High Energy Physics (HEP) is powered by x86 code, while the future will be much more reliant on accelerators like GPUs and FPGAs. The portable parallelization strategies (PPS) project of the High Energy Physics Center for Computational Excellence (HEP/CCE) is investigating solutions for portability techniques that will allow the coding of an algorithm once, and the ability to execute it on a variety of hardware products from many vendors, especially including accelerators. We think without these solutions, the scientific success of our experiments and endeavors is in danger, as software development could be expert driven and costly to be able to run on available hardware infrastructure. We think the best solution for the community would be an extension to the C++ standard with a very low entry bar for users, supporting all hardware forms and vendors. We are very far from that ideal though. We argue that in the future, as a community, we need to request and work on portability solutions and strive to reach this ideal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

A cost comparison of various hourly-reliable and net-zero hydrogen production pathways in the United States

Hydrogen (H 2 ) as an energy carrier may play a role in various hard-to-abate subsectors, but to maximize emission reductions, supplied hydrogen must be reliable, low-emission, and low-cost. Here, we build a model that enables direct comparison of the cost of producing net-zero, hourly-reliable hydrogen from various pathways. To reach net-zero targets, we assume upstream and residual facility emissions are mitigated using negative emission technologies. For the United States (California, Texas, and New York), model results indicate next-decade hybrid electricity-based solutions are lower cost ($2.02-$2.88/kg) than fossil-based pathways with natural gas leakage greater than 4% ($2.73-$5.94/kg). These results also apply to regions outside of the U.S. with a similar climate and electric grid. However, when omitting the net-zero emission constraint and considering the U.S. regulatory environment, electricity-based production only achieves cost-competitiveness with fossil-based pathways if embodied emissions of electricity inputs are not counted under U.S. Tax Code Section 45V guidance.

08 HYDROGEN↗