Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RuN”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Biochemical Conversion of Lignocellulosic Biomass to Hydrocarbon Fuels and Products: 2022 State of Technology and Future Research

The annual State of Technology (SOT) assessment is a key activity for biochemical platform research. It allows the impact of research progress (both directly achieved in-house at the National Renewable Energy Laboratory [NREL] and indirectly extrapolated from available public data) to be quantified in terms of economic improvements in the overall cellulosic biofuel production process for a particular conversion pathway. As such, initial benchmarks can be established for currently demonstrated performance and progress can be tracked toward out-year goals to ultimately demonstrate cost-competitive biofuel technology. Building upon efforts to progress NREL's biochemical platform R&D work toward ultimate 2030 goals to demonstrate less than $2.50/gallon gasoline equivalent (GGE) fuel selling prices, experimental and techno-economic analysis (TEA) activities have primarily focused on "advanced" biochemical processing strategies to fuels and coproducts, guided by TEA modeling to highlight key barriers and priorities toward achieving this goal across a number of potential bioconversion pathways. The purpose of the present effort is to benchmark the latest experimental developments for these pathways as quantified by modeled minimum fuel selling prices (MFSPs), as a measure of current status relative to those final targets. For this SOT, TEA models were run for two separate biological conversion pathways to fuels, based on available data for integrated biomass deconstruction and hydrolysate processing; namely carboxylic acids (primarily butyric acid) and diols (2,3-butanediol [BDO]), reflecting NREL's recently published 2018 biochemical design report focused on those two pathways. The models were run across two scenarios for lignin utilization, namely combustion and conversion to value-added coproducts.

09 BIOMASS FUELS↗
DOE OSTI2023

Dark Matter and Track Triggering with the CMS Experiment (Final Report)

This report summarizes the progress from DOE Early Career Award entitled “Dark Matter and Track Triggering with the CMS Experiment”. The first major goal of the project was to establish a sensitive Dark Matter (DM) search program for the CMS experiment in Run-2 of the LHC. Two related strategies were developed to accomplish this goal. First, we designed and executed searches for DM produced in association with heavy flavor quark pairs, top/anti-top (ttbar) and bottom/anti-bottom (bbbar). In addition, we sought to maximize the power of Run-2 DM searches through the development of a statistical combination of all major search channels within a consistent theoretical framework. The results of these aspects of the project are detailed in Section 1. The second major goal of the project was to develop a real-time “Level-1” tracking trigger system for the high-luminosity LHC (HL-LHC) CMS upgrade. Charged particle tracking in the first stage of the CMS trigger will be crucial for surviving the high-pileup environment expected at the HL-LHC. The L1 tracking trigger must process all front-end hit data sent from the inner Tracker and will have just 4us to output track primitives to the downstream L1 trigger without data loss. Our development of the CMS tracking trigger is described in Section 2. The team supported by the award is given in Section 3. The project resulted in 10 peer reviewed publications, a graduate student thesis, and contributions to the CMS Phase-2 Tracker Technical Design Report, as is detailed in Section 4. Our development work for the CMS tracking trigger has been integrated in the backend system architecture for the Tracker upgrade.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Short-Depth QAOA circuits and Quantum Annealing on Higher-Order Ising Models (Rev.2)

The Quantum Alternating Operator Ansatz (QAOA) and Quantum Annealing (QA) are quantum algorithms that are both based on the adiabatic theorem and both have the goal of sampling the optimal solution(s) of combinatorial optimization problems. Quantum annealing has been physically instantiated on D-Wave devices using superconducting flux qubits, and QAOA can be programmed on digital gate-model quantum computers such as the programmable superconducting transmon qubits devices of the IBMQ series, for instance ibm washington. QAOA and QA address the same types of problems, but it is unclear how they will scale to large problem sizes and to larger and higher-fidelity quantum computers. In this article, we present a direct comparison between QAOA, one and two rounds, run on all 127 qubits of ibm washington and QA run on D-Wave Advantage system4.1 and Advantage system6.1. The problems which allow for this comparison are random Ising model problems whose connectivity matches the heavy hexagonal lattice topology of ibm washington and the Pegasus graph connectivity of the two D-Wave devices. We create two classes of problem instances for this comparison: one with higher order terms (ZZZ variable interactions), linear terms, and quadratic terms, and a separate problem type with only linear and quadratic terms. Our QAOA circuits are novel and extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of ibm washington’s heavy hexagonal lattice and can be applied to future heavy-hex chips. We also test the effectiveness of the error suppression technique digital dynamical decoupling on the QAOA circuits. The QAOA circuits compiled to ibm washington are composed of several thousand circuit instructions, approximately 3, 000 depending on the details of the circuit, making these some the largest quantum circuits ever executed on a digital quantum processor. QAOA and QA are compared against the classical heuristic algorithm of simulated annealing and all problem instances are exactly solved using CPLEX in order to evaluate which samplers, if any, correctly found the ground state solution(s) of the problem instances. We find that (i) QA outperforms QAOA on all problem instances, (ii) QAOA samples the problems better than random sampling, and (iii) QAOA angle computation exhibits clear parameter concentration across the ensemble of Ising models.

127 Qubits↗
DOE OSTI2022

Vitrification of Hanford Tank Wastes for Condensate Recycle and Feed Composition Changeover Testing (Rev. 1)

During the vitrification of Hanford Site nuclear waste at the Waste Treatment and Immobilization Plant (WTP), which is a part of the safe and efficient retrieval, treatment, and disposal mission of the U.S. Department of Energy Office of River Protection, the offgas condensate generated from the waste-to-glass conversion is currently planned to be concentrated by evaporation in the Effluent Management Facility (EMF). This concentrated condensate can then be recycled back to the incoming waste and vitrified. To test the recycle process, an apparatus was designed and built to mimic the EMF evaporator and was then used to concentrate a volume of condensate produced during the vitrification of a sample of Hanford tank 241-AP-107 (referred to herein as AP-107) waste in a continuous laboratory-scale melter (CLSM). The concentrated condensate was added to an additional sample of AP-107 waste, to mimic one round of the recycle process, and the combined solution was vitrified, producing a second round of recycle condensate. In the current study, the EMF test apparatus was used to concentrate the second-round recycle condensate under evaporation conditions (at 45 °C and 1.4 psia) designed to emulate EMF operation. The condensate was successfully concentrated by a factor of ~10 while retaining over 95 % of the technetium-99 (99Tc), Cs, and I inventories in the concentrate. Another portion of AP-107 waste was retrieved by Washington River Protection Solutions, LLC (WRPS) and transferred to Pacific Northwest National Laboratory (PNNL), where it was pretreated and then combined with the second-round recycle AP-107 condensate concentrate and glass-forming chemicals (GFCs) to form the two-time recycle AP-107 melter feed, approximating a second round to the recycling action to be performed at the WTP. A portion of AP-105 waste was also retrieved by WRPS and provided to PNNL for pretreatment and combining with GFCs to form AP-105 melter feed. The two-time recycle AP-107 and AP-105 melter feeds were processed consecutively in the CLSM. The CLSM run proceeded for 13.63 hours, producing 9.70 kg of glass for an average glass production rate of 1464 kg m 2 d -1 during the two-time recycle AP-107 feed charging and 1568 kg m 2 d -1 during the AP-105 feed charging. The rate during AP-107 charging was essentially equivalent to the rate when processing no-recycle AP-107 feed and lower than that achieved when processing one-time recycle AP-107 feed. However, all rates were within the potential range of variability when processing melter feeds with similar composition in the CLSM. Likewise, the rate during AP-105 charging was higher than the previous rate processing AP-105, but within the potential CLSM range. The cold-cap characteristics changed from the typically thin AP-107 cold cap to a foamy-edged cold cap as previously seen with AP-105 shortly after transitioning to the AP-105 melter feed. The glass produced during the CLSM run was within 10 % of its target composition for the primary glass components. The CaO and Li 2 O targets varied by more than 1 wt% between the two-time recycle AP-107 and AP-105 glass targets and it took about 2 turnovers of the CLSM glass inventory to reach a relative chemical steady state in the glass for CaO and Li 2 O after the melter feed inputs were switched.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

β- and γ-Counting for Pre-detonation Nuclear Forensics on Eu-155

Post-detonation nuclear forensics was performed at Los Alamos National Laboratory (LANL) on 155Eu, a fission product on the wing of the fission product production curve whose yield is sensitive to fission fuel and neutron energy. With a half-life of 4.753 years, 155Eu provides a longer-lived option for these measurements than other fission products with similar mass numbers. The Chemistry Division Group - Nuclear and Radiochemistry, at LANL routinely measures a suite of fission products from 235U fissions in thermal neutron flux experiments known as thermal calibration exercises, using a mixture of gas proportional β-decay counting and γ-spectrometry on HPGe detectors. The fission products of interest are reported relative to a high-yield reference fission product from the same sample to create a running average ratio specific to neutron energy and fuel type; Equation 1 below shows the ratio-of-ratios R-value measured fission products are reported in: R i x = [A( i X)/A( 99 Mo)] unknown [A( i X)/A( 99 Mo)] 235 U n,th where A denotes activity, i X is the nuclide of interest, and 235 U n,th denotes irradiations of 235 U with thermal neutrons. Measurements of unknown fission spectra are ratioed to a running-average of thermal calibration results, and the resulting R-value can be referred to as against a library of irradiation conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Optimizing Performance on Trinity Utilizing Machine Learning, Proxy Applications and Scheduling Priorities

The sheer number of nodes continues to increase in today’s supercomputers, the first half of Trinity alone contains more than 9400 compute nodes. Since the speed of today’s clusters are limited by the slowest nodes, it more important than ever to identify slow nodes, improve their performance if it can be done, and assure minimal usage of slower nodes during performance critical runs. This is an ongoing maintenance task that occurs on a regular basis and, therefore, it is important to minimize the impact upon its users by assessing and addressing slow performing nodes and mitigating their consequences while minimizing down time. These issues can be solved, in large part, through a systematic application of fast running hardware assessment tests, the application of Machine Learning, and making use of performance data to increase efficiency of large clusters. Proxy applications utilizing both MPI and OpenMP were developed to produce data as a substitute for long runtime applications to evaluate node performance. Machine learning is applied to identify underperforming nodes, and policies are being discussed to both minimize the impact of underperforming nodes and increase the efficiency of the system. In this paper, I will describe the process used to produce quickly performing proxy tests, consider various methods to isolate the outliers, and produce ordered lists for use in scheduling to accomplish this task.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Nb SRF cavities applied to searches for beyond the standard model physics

The SRF technology developed for accelerators can be successfully applied to new applications, including quantum computing, dark matter searches and beyond the standard model physics. The ultra-high quality factor of SRF cavities can al-low to achieve unprecedented sensitivity in fields outside of the usual accelerator applications, for examples fundamental and beyond the standard model physics. Applications of SRF cavities for gravitational waves searches are also being investigated. The SQMS Physics and Sensing thrust is working on this effort as it strives to combine SRF cavities with QIS, with a focus on BSM and fundamental physics. Using Nb SRF cavities, we run a proof of principle search for the dark photon dark matter candidate at fixed frequency demonstrating the deepest sensitivity to wavelike dark mater [1]. The Dark SRF experiment uses two 1.3GHz SRF cavities to conduct a light shining through wall (LSW) experiment to attempt producing and detecting dark photons in the laboratory. The results of this first superconducting LSW search were recently published [2] and we are now working on the second phase of the experiment, which will be run in a dilution refrigerator using 2.6GHz SRF cavities. In addition to the limits we have achieved, we will also present the experiments we are currently developing, such as the heterodyne axion DM search [3, 4], the SERAPH series of experiments, and gravitational waves searches based on SRF cavities as proposed in [5-8].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Carbon, Nitrogen, and Sulfur Analysis of the Cetama Viognier Standard Reference Material

The Viognier sample was analyzed for C, N, and S concentrations and stable isotope compositions using an Elementar Vario Isotope Cube Elemental Analyzer (EA) that is coupled to an IsoPrime PrecisION IRMS. Powder samples were placed in tin capsules and then loaded onto a rotary autosampler. The autosampler dropped samples into the EA, and samples were combusted at 1175 °C over tungsten oxide in a continuous stream of helium carrier gas. A pulse of oxygen is added to the gas stream resulting in flash combustion of the tin capsule containing the sample, which raises the reaction temperature to approximately 1800 °C for a few seconds. The combustion reaction produces SO 2 , N 2 and CO 2 from any sulfur, nitrogen and carbon present in the sample. The resulting gases were then passed through a reduced copper reactor that was heated to 850 °C, to reduce NO x to N 2 , reduce SO 3 to SO 2 , and trap any volatile halogen compounds on silver wool. Following water removal using an adsorption tube, the N 2 , CO 2 , and SO 2 analyte gases were separated and purified using purge-trap columns. The purified gases were then carried through a thermal conductivity detector. The detector signal was passed to software that calculated elemental abundances based on integrated peak areas. The sample gases were then passed to the IRMS and stable isotope ratios were measured. Raw sample peak areas were corrected by subtracting the average peak area from blanks consisting of empty tin capsules run using the same EA analysis method. Blanks were also run following each replicate to verify that all the material was combusted. The Viognier sample was analyzed in triplicate on two days (target sample masses: 30, 60 and 90 mg). The following standards were analyzed to calibrate EA-IRMS measurements: IAEA-C6, USGS-40, USGS-41, IAEA-S1, and IAEA-S2.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Bench-Scale Testing of a High Efficiency, Ultra-Compact Process for Pre-Combustion CO 2 Capture (Final Report)

The overarching objective of this bench-scale study is to field-validate the technical feasibility of a membrane- and adsorption-enhanced water gas shift (WGS) reaction process that employs a carbon molecular sieve (CMS) membrane reactor (MR) followed by an adsorption reactor (AR) for pre-combustion CO 2 capture while demonstrating progress towards achievement of the overall performance goals of CO 2 capture with 95% CO 2 purity at a cost of electricity (COE) of 30% less than the baseline capture approaches. The main goal of this bench-scale project is to advance the proposed MR-AR transformative CO 2 capture technology to a sufficient maturity level so that our present industrial partners (M&PT and Worley) can potentially adapt the technology for further small pilot-scale testing. The project began at TRL 4, as the system prototype had already been validated in the laboratory on simulated syngas (as part of project DE-FE0026423). The project ended at TRL 5, via scaling-up of the prototype system and testing it on actual syngas at a host site (CAER at the University of Kentucky (UKy)). The project was carried-out in two different phases: In Phase I (during BP1 of the project), the team designed, constructed, and assembled the bench-scale experimental MR-AR system, prepared the membranes, adsorbents, and catalysts, tested the bench-scale unit with simulated syngas to validate functionality, and prepared a preliminary TEA of the technology. In Phase II (during BP2 of the project), the team installed the unit at the test site at the UKy site and completed all utility connections and hookups, field-tested the novel MR-AR process in the bench-scale system using real syngas, collected and analyzed experimental data, and completed a detailed TEA of the technology. A key aim of the project was to identify and address the technical and process risks and to generate information to advance the technology to the next stage of development. All project milestones and success criteria were met. Specifically, the team: (i) designed, constructed and assembled a bench scale system that handled 1 – 5 scfm syngas and tested its functionality with simulated syngas; (ii) prepared high performance CMS membrane tubes (ID: 3.5 mm, OD: 5.7 mm, 30” long) that meet the target H 2 permeance ( > 1 m 3 /(m 2 .hr.bar) or (> 370.3 GPU)) and a target H 2 /CO selectivity >80 at the relevant temperature (up to 300oC) and pressure conditions (up to 25 bar) with a <10% decline in performance over each 250 hr. testing period (iii) procured commercial sour-shift catalyst in sufficient quantity and prepared up to 10 kg of pelletized adsorbent for use at relevant conditions (250oC 2.5wt.% and target sorbent attrition rate < 0.2; (iv) installed the bench-scale unit at the test site (UKy) and tested the MR-AR technology using real syngas for over 250 hr. for both static and flow experiments; (v) updated the TEA analysis based on bench-scale data and met CO 2 capture goals of 95% CO 2 purity at a cost of electricity 30% less than baseline capture approaches. The CMS membranes, HTC adsorbents and catalysts employed all exhibited very robust and stable performance during the long-term run (over a >250 hr live syngas run). Furthermore, the proposed MR-AR IGCC system achieved a LCOE with a N 2 sale price of $\$$ 30/ton, with N 2 compression (no N 2 compression) of 101.2 $\$$/MWh (95.3 $\$$/MWh), which represents a 29.8% (33.9%) LCOE reduction in the baseline IGCC with carbon capture of 144.2 $\$$/MWh. The proposed MR-AR IGCC delivers a CO 2 capture cost of 44.8 $\$$/tonne (39.9 $\$$/tonne) vs 98.1 $\$$/tonne of the baseline capture case and a net power production of 586 MWe (623 MWe) vs 556 MWe of the baseline capture case.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

JIRIAF: JLAB Integrated Research Infrastructure Acros Facilities

The JIRIAF project aims to combine geographically diverse computing facilities into an integrated science infrastructure. This project starts by dynamically evaluating temporarily unallocated or idled compute resources from multiple providers. These resources are integrated to handle additional workloads without affecting local running jobs. This paper describes our approach to launch best-effort batch tasks which exploit these underutilized resources. Our system measures the real-time behavior of jobs running on a machine and learns to distinguish typical performance from outliers. Unsupervised ML techniques are used to analyze hardware-level performance measures, followed by a real-time cross-correlation analysis to determine which applications cause performance degradation. We then ameliorate bad behavior by throttling these processes. We demonstrate that problematic performance interference can be detected and acted on, which makes it possible to continue to share resources between applications and simultaneously maintain high utilization levels in a computing cluster. We relocate the CLAS12 data processing workflow to a remote data center for a case study, preventing file migration and temporal data persistency.

Lawrence, David↗
DOE OSTI2023

Batched Sparse Linear Algebra (Final Report for Subcontract B648960)

This report finalizes design specifications for developing batched kernels for small tensor operations for unassembled matrix-free iterative solvers, batched solvers for partially assembled operators, and batched solvers with support for various sparse formats. The outcome of the project milestones is a set of interfaces to Batched Sparse LA solvers running on hardware accelerators for use in ECP Libraries and Applications. It is part of the development of sparse batched kernels, solvers/preconditioners as well as creating interoperability in xSDK libraries with sparse and dense batched functions to benefit ECP applications. The participants included representatives from ECP libraries (not limited to the xSDK project), applications, and vendors (AMD, Intel, and NVIDIA). Batched sparse linear algebra solvers form the new frontier for algorithmic development and performance engineering. Many applications (ECP and non-ECP alike) require simultaneous solutions of small linear systems of equations that are structurally sparse. To move towards high hardware utilization, it is important to provide these applications with appropriate interfaces to efficient batched sparse solvers running on modern hardware accelerators. We present interface designs in use by HPC software libraries supporting batched sparse linear algebra and the development of sparse batched kernel codes for solvers and preconditioners. We also address the potential interoperability opportunities to keep the software portable between the major hardware accelerators from AMD, Intel, and NVIDIA. The presented interface specifications includes batched band, sparse iterative, and sparse direct solvers. This report summarizes progress in Kokkos Kernels and the xSDK libraries MAGMA, Ginkgo, hypre, SUNDIALS, and SuperLU_dist.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Celeritas R&D Report: Accelerating Geant4

Celeritas is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a Graphics Processing Unit (GPU)-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.4, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling straightforward integration of Celeritas into the high energy physics (HEP) experimental frameworks CMSSW and ATLAS FullSimLight. On the Perlmutter supercomputer, Celeritas performs EM physics between 3× and 18× faster using the machine’s Nvidia GPUs compared to using only CPUs, corresponding to an electrical power efficiency up to a factor of 5. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.7–2.2× on GPU and 1.2× on CPU. In a CMS test application using tt¯ events and the prototype Run 4 configuration, compared to Geant4 CPU, Celeritas with a Nvidia A100 improves overall throughput up to a factor of 2.7× but cannot be efficiently shared with more than 8 cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

MultiPEM Toolbox: User Manual [Rev. 2]

This document explains use of the Multi-Phenomenology Explosion Monitoring (Multi PEM) Toolbox, a collection of R scripts for estimating the unknown device parameters of a new event with uncertainty quantification. The methodology and application used for illustration in this user manual are fully documented in a Los Alamos National Laboratory technical report hereafter designated “WPA” for reference. Additional details on the application are found in a recent journal article. Two assessment types are available: rapid and complete. Rapid assessments are conducted in two stages, as described in Section 2. In the first stage, calibration data are used to estimate forward and error model parameters (WPA, §5.1) and (if relevant) errors-in-variables yield values for calibration sources (WPA, §3, Equation (3)). In the second stage, new event data are used to estimate the unknown new event device parameters (WPA, §5.2) with uncertainty quantification. Two options for treating the inferred first stage parameters in second stage Bayesian analysis are available: fixing them at their maximum likelihood estimate (default), or multiple imputation. Multiple imputation involves utilizing several posterior samples (imputations) of the first stage parameters as fixed values in the second stage posterior sampling of the new event device parameters. Second stage sampling is conducted across imputations in parallel to improve computational efficiency. This method produces improved uncertainty quantification of the new event device parameters compared with the default treatment of the first stage parameters, at the expense of additional computation. Complete assessments are conducted in a single stage, as described in Section 3. Calibration and (if relevant) new event data are used simultaneously to estimate all forward model, error model, and (if relevant) new event device parameters with uncertainty quantification on the latter. As the name suggests, rapid assessments generally run substantially faster than complete assessments (even with multiple imputation), because the results of first stage analysis can be stored and incorporated into estimating a relatively low-dimensional space of new event device parameters whenever relevant new event data becomes available. On the other hand, complete assessments must be run on the full set of model and device parameters with calibration and new event data every time the latter becomes available.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Enhanced Component Performance Study: Emergency Diesel Generators 1998–2022

This report presents an enhanced performance evaluation of the emergency power system (EPS) and high-pressure core spray (HPCS) emergency diesel generators (EDGs) at U.S. commercial nuclear power plants. This report evaluates component performance over time using (1) Institute of Nuclear Power Operations (INPO) Industry Reporting and Information System (IRIS) data from 1998 through 2022 and (2) maintenance unavailability performance data from Mitigating Systems Performance Index (MSPI) Basis Document data from 2002 through 2022. The objective is to show estimates of current failure probabilities and rates related to EDGs, trend these data on an annual basis, determine if the current data are consistent with the probability distributions currently recommended for use in Nuclear Regulatory Commission (NRC) probabilistic risk assessments, show how the reliability data differ for different EDG manufacturers and for EDGs with different ratings; and summarize the subcomponents, causes, detection methods, and recovery associated with each EDG failure mode. The EDG failure modes considered are fail to start (FTS), fail to load and run (FTLR), and fail to run after one hour of operation (FTR>1H). Engineering analyses were performed with respect to time-period and failure mode without regard to the actual number of EDGs at each plant. The factors analyzed include subcomponent, failure cause, detection method, recovery, manufacturer, and EDG rating. The following increasing trends were identified for EDGs for the most recent 10-year period: (1) HPCS EDG FTS failure probability, (2) EPS and HPCS EDG frequency of start demands (demands per reactor year), and (3) EPS and HPCS EDG frequency of FTLR demands. The following decreasing trends were identified for EDGs for the most recent 10-year period: (1) EPS EDG FTS failure probability, (2) EPS EDG FTLR failure probability, (3) EPS EDG unavailability, and (4) EPS and HPCS EDG frequency of FTLR events (failures per reactor year).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Excerpt from Oxygen Monitor Clean Testing Report Performed in 2010

The Rosemount X-STREAM Series X2GP is a paramagnetic oxygen analyzer with a 0-25% oxygen measurement range. It claims a drift less than 2% of span (0.10% volume oxygen) per week, and a response time of less than 5 seconds (for the analyzer to report 90% of a step change). Our test equipment was set-up such that the same gas was running through each oxygen analyzer at the same time. This allowed Calibration Error, Calibration Drift, and Loss of Power Testing to be run on all monitors at the same time, eliminating any variability between the analyzers from tolerances in the gas cylinder, ambient temperature, ambient pressure, test duration, etc. An APTech pressure regulator was placed on the gas cylinder manifold and a 23 psia Paroscientific pressure transducer measured the regulator’s outlet pressure. This did not give the inlet pressure on any one analyzer but outlet pressure from the manifold before the stream was split. The APTech pressure regulator was used to keep operating conditions consistent. A DWYER rotameter was installed on each oxygen monitor’s outlet piping. With the exception of Variation of Flow Rate Testing, the rotameter was left fully open for all testing. The flow rate through each oxygen monitor was controlled by a Swagelok needle valve that was installed on each analyzer’s inlet piping. Additionally, a filter was installed in the each monitor’s inlet piping to keep particles from the measurement cells and a Swagelok ball valve was installed on each monitor’s outlet piping to isolate each monitor between tests. No pumps were installed. Each of the cylinders used was at a pressure greater than ambient room barometric pressure. This pressure differential was used to move gas through the testing equipment. The flow rate through each analyzer was controlled by the Swagelok needle valve installed on each analyzer’s inlet piping. Finally, for all tests, a Paroscientific Digiquartz Barometer was used to measure ambient room barometric pressure and a Fluke Thermo-hygrometer was used to measure ambient room temperature.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Characterization of Non-Science Grade DESI CCDs and Creating Additional ICARUS Monitoring Metrics

DESI, otherwise known as the Dark Energy Spectroscopic Instrument, is a tool to measure optical spectra from the effects of dark energy on the universe and has been running for 3 years. My work delves into the data taking/analysis side of the DESI CCDs which will help in swapping out broken CCDs on the instrument. Additionally, I worked on making new metrics for ICARUS, Imaging Cosmic and Rare Underground Signals for their monitoring of run data. This paper will talk about the types of characterizations done for DESI and some metrics for data quality and tools used for ICARUS monitoring.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗