Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reproducibility”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Propagation of a Stress Pulse in a Heterogeneous Elastic Bar

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material in the limit of long and moderately long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analytic fits to atom-in-jellium shear modulus predictions

Fits to AJ shear modulus calculations from ρ0 to the limit of the corresponding AJ EOS table (usually 1000 ρ0), adjusted to reproduce STP shear modulus where possible. Fitting minimized the fractional difference, so low-pressure points would be reproduced accurately despite the wide-ranging tabulation. AJ calculations were at T=0. Isochoric variation G(T) has been small for all cases spot-tested so far, and may have either sign. The fitted equation usually matches the AJ data to within a few percent. Between numerical noise in the AJ calculation and probably-physical structure not captured by the equation, the deviation could be up to 20% in some places in most models, and 30% in a few. AJ does not capture structural phase transitions or phases stabilized by directional bonds, where the shear modulus may vary by a greater amount. Where possible, STP ρ,G were used as parameters, and low pressure AJ points were de-weighted or removed if necessary. Otherwise, such as where AJ fails to capture solid phases with a significantly different shear modulus, the AJ data were fitted as far down in pressure as possible, G0 was also fitted if necessary, and ρ0 was also adjusted if needed to keep G0>0. The resulting models are not intended for use at low pressure, though some are probably adequate for practical purposes

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-driven learning of nonlocal models: from high-fidelity simulations to constitutive laws

We show that machine learning can improve the accuracy of simulations of stress waves in one-dimensional composite materials. We propose a data-driven technique to learn nonlocal constitutive laws for stress wave propagation models. The method is an optimization-based technique in which the nonlocal kernel function is approximated via Bernstein polynomials. The kernel, including both its functional form and parameters, is derived so that when used in a nonlocal solver, it generates solutions that closely match high-fidelity data. The optimal kernel therefore acts as a homogenized nonlocal continuum model that accurately reproduces wave motion in a smaller-scale, more detailed model that can include multiple materials. We apply this technique to wave propagation within a heterogeneous bar with a periodic microstructure. Several one-dimensional numerical tests illustrate the accuracy of our algorithm. The optimal kernel is demonstrated to reproduce high-fidelity data for a composite material in applications that are substantially different from the problems used as training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

National Ignition Campaign Data Analysis by an ICF Outsider

A simple DT burn model, with no alpha-particle warming feedback, using a combination of ICF inputs and the measured ion temperature, T, volume compression, C, ablator velocity, v, and burn width time, Δt, can reproduce the national-ignition campaign (NIC) fusion yields with a rms percent difference of ~40%. Analysis suggestion the fusion-yield reproducibility of NIC events with nominally identically measured T, C, v, and Δt is <30%. The highest fusion yields, at a fixed laser energy are obtained with depleted-uranium (DU) hohlraums and high-density carbon (HDC) ablators. Small diameter fit tubes (≤5 μm) enhance fusion yields by an average of 26±15% relative to 10-μm diameter fit tubes. A dependence of fusion yield on the diameter of the capsule relative to the diameter of the hohlraum is inferred for fixed T, C, v, and Δt. No evidence of conditions close to ignition, or alpha-particle warming is found in the existing NIC data with LEH< 2MJ. To obtain a predictive model for extrapolation proposes we map the T, C, and v into power laws as a function of LEH. The above mentioned extrapolatable model has a rms percent difference from the yields of the second half of the NIC of (56 ± 4)%. Higher (and significantly lower) DT temperatures are obtained with HDC ablators (and low-foot drives). The highest compressions are obtained via low-foot drives, but this does not compensate for the corresponding lower temperatures. An extrapolation of the presented model suggests significant alpha-particle heating at an LEH of ~5 MJ. It is hoped that others more familiar it the conditions needed for ignition will use to presented model to estimate the LEH required for ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dip leg plug formation evaluation and abatement testing in support of the H-Canyon accelerated basin de-inventory

As requested in the Technical Task Request (TTR), the Savannah River National Laboratory (SRNL) performed testing on the current typical bubbler dip leg design used in the H-Area Canyon (HCA). The TTR has a functional classification of General Services. Testing was initiated to ascertain if the plugging described in the TTR could be reproduced and, if so, to obtain general metrics on the resulting plugging dynamics, and evaluate possible mitigating methods, including design changes and cleaning/purging protocols. The bubblers were to be suspended in a simulant typical of H-Canyon solutions containing aluminum nitrate and silica at a temperature of 25°C using a nominal bubbler air flowrate of 7.5 SCFH. Emphasis was placed on minimizing equipment changes necessary to accomplish the mitigations, as equipment changes in the existing facility would be costly. Testing evolved into four (4) phases, concluding that plugging could indeed be reproduced (believed to be attributed to localized evaporation and bridging of crystals in the dip legs). During testing, plug formation information was obtained, including time, pressure, and temperature traces leading to plugging were defined; solution compositions leading to plugging were determined, and various physical characteristics (visual and mechanical) of the crystallized plug were obtained.

42 ENGINEERING↗
DOE OSTI2021

Toward Well-Provenanced Computer System Benchmarking: An Update

This report provides an update on the current design, implementation, and usage of bueno, a Python framework enabled by container technology that supports gradations of reproducibility for wellprovenanced benchmarking of sequential and parallel programs. The ultimate goal of the bueno project is to provide convenient access to mechanisms that aid in the automated generation, collection, and dissemination of data relevant for experimental reproducibility in computer system benchmarking.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Integration of Omics into a New Comprehensive Rate Law for Competitive Terminal Electron-Accepting Processes in Reactive Transport Models: Application to N, Fe, S, and Contaminant Transformations in Stream and Wetland Sediments

Surface waters represent important sources of alternate energy and drinking water in the United States, and characterizing the biogeochemical processes that affect surface water quality is relevant to the DOE-BER mission. Sediment biogeochemical processes regulate the release of carbon (C), nutrients, and contaminants to surface waters and thus influence water quality. Sediment biogeochemical processes are dynamic and affected by the deposition and remobilization of solid material and changes in environmental conditions driven by water discharge variations. Wetlands are important natural filters of surface waters which may either trap, metabolize, or mobilize nutrients and contaminants. Despite their importance, biogeochemical processes regulating nutrient and contaminant release and C transformation in stream and wetland sediments cannot be predicted accurately by current mathematical models. These reactive transport models largely rely on detectable changes in geochemical conditions to activate metabolic processes, do not accurately account for the competition between microbial processes, and poorly constrain effects of hydrological perturbations on biogeochemical processes. In this BER-SBR exploratory project, metagenomic and geochemical signatures were combined to identify microbially-mediated redox processes in anaerobic stream and wetland sediments from the Savannah River Site (SRS, ANL SFA) and East Fork Poplar Creek (EFPC, ORNL SFA) that play important roles in C, uranium (U), and mercury (Hg) transformations. In addition, sediment incubations were conducted to examine the competition between anaerobic respiration processes Finally, new rate laws were developed for reactive transport models that rely on complementary metagenomic and geochemical signatures to identify the underlying anaerobic microbial processes in stream and wetland sediments, describe the competition between the dominant metabolic processes involved in nutrient release and U and Hg mobilization, and more accurately quantify carbon transformation and the response of microbial processes to changes in redox conditions associated with hydrological forcing. These rate laws were optimized in batch reactors with SRS wetland sediments, where iron and sulfate reduction dominate. Anaerobic carbon remineralization processes followed the expected thermodynamic sequence of microbial respiration with depth in both sediments, except that geochemical signals indicated that sulfate reduction was inactive in EFPC sediments and moderate in SRS wetland sediments. Estimates indicated that microbial iron reduction contributed to at least half of the production of reduced iron in these sediments. Incubations demonstrated that nitrate reduction, denitrification, and dissimilatory nitrate reduction to ammonium were active in the natural EFPC sediment and activated upon nitrate amendment in these nitrate-rich sediments. In turn, these processes were outcompeted by the addition of either iron oxides or sulfate as alternative terminal electron acceptors. Although geochemical products of sulfate reduction were not detected in the incubations, the abundance of sulfate reduction genes increased with depth in the sediment and was equally more pronounced in treatments amended with either iron oxides or sulfate. Simultaneously, anaerobic sulfide oxidizing bacteria coupling sulfide oxidation to DNRA (and not conventional denitrification) were apparently enriched over time, regardless of the treatments. These findings indicate that sulfate reduction is important in freshwater stream sediments and probably catalyzed by a cryptic sulfur cycle involving nitrogen species, in which the sulfur products from sulfate reduction are immediately removed by side reactions and not detectable by geochemical measurements alone. Similar experiments in SRS sediments, however, demonstrated little interaction between nitrogen and sulfur cycling microorganisms. Sulfate reduction was impacted by the addition of more thermodynamically favorable electron acceptors, suggesting either that iron-reducing microorganisms outcompeted sulfate-reducing microorganisms for organic substrate, depleting the stock of electron donor available for sulfate reduction, or that the cryptic sulfur cycle was shunted by the precipitation of FeS generated as a result of the abiotic reduction of iron oxides by dissolved sulfide. A diagnostic modeling exercise was conducted to further investigate the competition between terminal electron accepting processes. As conventional kinetic models typically do not account for cryptic cycles and used inaccurate formulations to describe competition between microbial communities, new metabolic rate laws were developed that explicitly express the electron acceptor-specific enzyme of each energetically favorable metabolic process based on gene abundance detected in the incubations. The model was tested with the sediment slurry incubation data to determine whether substrate competition could explain the decrease in sulfate reduction observed in the presence of iron oxide competitor. The model was able to reproduce geochemical concentrations really well in each treatment once the model was calibrated with the unamended control, suggesting that microbial competition was indeed driven by thermodynamic considerations. Overall, carbon remineralization processes and rates will be reproduced much more realistically with the new metabolic rate laws.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The tracer-BHR4 model for subgrid mixing layers

A model for the passive evolution of turbulent subgrid-scale mixing layers is developed based on the model of tracer-BHR model of Braun & Gore. The new model, tracer-BHR4, mainly differs from the previous iteration in that it aims to replicate the behavior of BHR4, instead of BHR3.1. The coefficients from BHR4 are used and a new model for mixing layer growth is developed to allow the tracer model to reproduce the de-mixing behavior seen in BHR4. Additionally, an improved timestep size controller is introduced. The new model is tested on problems containing Rayleigh-Taylor or Richtmyer-Meshkov mixing layers, and is able to reproduce the behavior of well resolved BHR4 simulations with reasonable accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Non-Destructive In-Process Assessment of Thermal Spray Repairs

This exploratory project, entitled “Non-Destructive In-process Assessment of Thermal Spray Repairs,” evaluated the capabilities and limitations of four non-destructive testing (NDE) technologies for detection of defects in as-sprayed thermal spray coatings: acousto-ultrasonics, vibro-thermography, thermal wave infrared (IR) imaging, and flash IR thermography. The project goal was to identify a single technology suitable for use as an in-process inspection for two common classes of defects in thermal spray coatings, porosity and disbonds. Two rounds of test coupons were fabricated with intentionally seeded defects in AISI 420 stainless steel coatings deposited by twin wire arc (TWA) spray onto gray cast iron substrates. These standardized test coupons allowed the sensitivity of the four NDE methods to be directly compared. The results indicated that the acousto-ultrasonic method will not meet requirements for accuracy, inspection time, and reproducibility, while the vibrothermography method will not meet requirements for integration within remanufacturing cells and process flows. The IR methods, thermal wave imaging and flash thermography, both exhibited acceptable accuracy and reproducibility; however, the former method does not meet the requirement for inspection time. Flash IR thermography met all requirements in the initial round of testing, and that result was confirmed in a second round of testing on samples with a larger, more complex geometry. A final recommendation for further development of flash IR thermography was therefore made based upon i) detection performance, and ii) feasibility of deployment into the target application of resurfacing CAT engine blocks and headers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The importance of phonon spectrum measurements for accurate prediction of thermal scattering laws [Slides]

Although approximate phonon spectra (e.g. Debye model) are sufficient to reproduce integral quantities (e.g. total INS cross section) in reasonable agreement with the measured data, more accurate phonon spectra are needed to achieve agreement with the differential data. SANS component is essential for reproducing the total scattering cross section; this is the case for nanodiamonds and other materials where SANS components are present (such as porous materials like nuclear graphite). A format proposal for SANS data in ENDF or GNDS is needed. Therefore, differential (e.g. phonon spectrum) and integral (e.g. total cross section) measurements are necessary in the process of creating an ENDF evaluation. We have successfully developed a method by which ENDF evaluations can be directly compared to the relevant experimentally measured quantities. INS experiments at SNS, ILL, ISIS, and, in the future, ESS, as well as total cross section measurements at the VESUVIO instrument at ISIS, are accessible through user programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Energy-dependent optimization of the prompt fission neutron spectrum with CGMF

Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Bridging the length scales on mechanical property evaluation (Final Report)

The development of small-scale mechanical testing in combination with microstructural investigation is of great interest to the nuclear materials community for both materials development and monitoring applications. Dramatically reducing the sample sizes to reduce radioactivity and obtaining mechanical properties of irradiated samples is truly intriguing. Moreover, such studies promise a range of benefits including cost reduction, fundamental insight in structure-property relationships, increased statistics on less sample material, and reinvestigation of prior irradiated and tested reactor samples while simultaneously enabling the generation of mechanical test data on ion beam irradiated materials with limited penetration depths. Small scale materials testing on sub-sized samples has been studied for several decades, though it has only been after the development of micro-testing based on Focused Ion Beam (FIB) sample manufacturing in that orders of magnitude smaller samples could really be investigated in a quantitative manner. In recent years, small scale mechanical testing techniques at a number of length scales has been developed for both unirradiated and irradiated (ion and neutron) materials. Technological advances made in this field have enabled ex-situ and in-situ transmission electron microscopy (TEM) and scanning electron microscopy (SEM) examination, thus leading to more accurate measurements as well as additional mechanistic information. A recent review of the benefits of these techniques show that these techniques are at a stage to tackle multi-scale ranges of materials investigations and can be utilized to obtain fundamental science-based understanding of nuclear materials. Considering the tremendous advances made, one can see how small-scale mechanical testing techniques combined with modeling can enable true small scale to bulk scale mechanical property correlations. However, for the engineering community to adapt this approach fully, one needs to demonstrate that a) that these techniques can produce results with high fidelity and reproducibility, b) generate engineering stress-strain data that one can utilize to understand bulk behavior, and c) generate new insight into relevant phenomena fostering the true understanding of radiation damage and microstructure in materials for nuclear applications. It is the objective of this proposal to bridge the length scale between macro- and micro- scale mechanical testing of unirradiated and irradiated materials. This involves the development and demonstration of procedures for multi-scale mechanical testing that enable high fidelity reproducibility of data and the generation bulk property data from small-scale mechanical tests. Through this, the proposal aims to enhance the confidence in the obtained data at the smaller length scales and enhance the insight provided from these techniques for bulk scale applications on both unirradiated and irradiated materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Consistent Evaluation of the Prompt-fission Neutron Spectrum and Multiplicity for n+ 235,238 U and n+ 239 Pu

This report was written to satisfy a FY20 NCSP milestone on 235,238 U and 239 Pu. The milestone requires to “finalize a report assessing our methodology to evaluate prompt-fission neutron spectrum (PFNS) and multiplicity consistently”. More specifically, we study whether the code CGMF can reproduce ENDF/B-VIII.0 evaluated PFNS and average prompt-fission neutron multiplicities, $\bar{v}$, for 235,238 U and 239 Pu using one joint parameter set per isotope. If CGMF is shown to be able to reasonably reproduce ENDF/B-VIII.0 within its model-parameter space, this code could be used for future consistent evaluations of PFNS and $\bar{v}$. To answer this question, we explore here the parameter space of CGMF and its impact on calculated values and whether they are close to evaluated and experimental data. We also list experimental data that would enter a future evaluation and statistics method that could be used to obtain evaluated data and covariances. We conclude that values of $\bar{v}$ calculated by CGMF are reasonably close to ENDF/B-VIII.0 data, while more work on modeling the PFNS is needed (parameter optimization and model improvements) to reliably use it for evaluations.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Reassessing the MCNP Random Number Generator

Random number generators are integral components to Monte Carlo codes. They provide the pseudorandom number sequence used to actually sample the distributions of interest. As a result, they are one of the most important components to the software. The current recommended MCNP random number generator is a 63-bit linear congruential generator (LCG). This generator is quite fast, but it has some drawbacks. First, it only has a period of 2 63 . Due to the necessarily non-optimal usage of random numbers to ensure parallel reproducibility, this amount is too few to guarantee random number sequences are not reused in all configurations the code runs under. As simulation size increases, users will need to be aware of the limitations of the generator and tune configuration variables to best suit their simulations, or they will need to assume that reuse is not negatively affecting their answers. Neither of these are optimal. Second, small LCGs are fairly weak in bit generation quality, and this can have an unknown impact on the quality of the simulation. This paper is an investigation into whether or not more modern random number generators can supersede the current ones. The goal is to find a generator that is similar or superior in speed to the LCGs, has a state space large enough to make strong guarantees about random number reuse, and passes all modern random number test suites. If such a generator is found, it would eliminate the need for the user to even be aware of the limitations of the random number generator and would simplify the use of the code. This paper will be broken into several parts. Sec. 2 will discuss the evolution of the random number generator within the MCNP code. Sec. 3 will go over what a Monte Carlo code needs from a generator to be reproducible and portable and how the current generator behaves in that light. Sec. 4 goes through how each generator was tested. Finally, Sec. 5 will discuss improvements that could be made to the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Addendum to ICSBEP Logbooks 12r, 13r, 14r and 15r for Polyethylene and Graphite Reflected HEU Metal Critical Experiments

This report documents the experimenters’ data sheets that were recently (2023) discovered for critical experiments with highly enriched uranium metal with polyethylene and graphite reflectors. These data sheets were produced at the time of the measurements, when the dimensional inspection reports for the polyethylene and graphite were available. However, at this writing, those inspection reports are not available but may be in unmarked storage for the Y-12 National Security Complex. However, the data sheets presented herein contain sketches of the experimental configurations and other additional information that is not in the logbooks, such as dimensions and masses of the graphite measurements. This report reproduces the experimental data sheets that complement the logbooks. Some of the data on these sheets and the logbooks can be used to infer the dimensions and masses that are not documented. Some obvious mistakes in the logbook and data sheets have been corrected. In the reproduced data sheets, some of the information is not clearly visible: some information, such as average values, can be obtained from other data in the sheets or logbook. This report should be sent to the International Criticality Safety Benchmark Program at Idaho National Laboratory to complement the existing Oak Ridge Critical Facility logbooks there.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Integrating HPC, AI, and Workflows for Scientific Data Analysis: Report from Dagstuhl Seminar 23352

The Dagstuhl Seminar 23352, titled “Integrating HPC, AI, and Workflows for Scientific Data Analysis,” held from August 27 to September 1, 2023, was a significant event focusing on the synergy between High-Performance Computing (HPC), Artificial Intelligence (AI), and scientific workflow technologies. The seminar recognized that modern Big Data analysis in science rests on three pillars: workflow technologies for reproducibility and steering, AI and Machine Learning (ML) for versatile analysis, and HPC for handling large data sets. These elements, while crucial, have traditionally been researched separately, leading to gaps in their integration. The seminar aimed to bridge these gaps, acknowledging the challenges and opportunities at the intersection of these technologies. The event highlighted the complex interplay between HPC, workflows, and ML, noting how ML has increasingly been integrated into scientific workflows, thereby enhancing resource demands and bringing new requirements to HPC architectures, like support for GPUs and iterative computations. The seminar also addressed the challenges in adapting HPC for large-scale ML tasks, including in areas like deep learning, and the need for workflow systems to evolve to leverage ML in data analysis fully. Moreover, the seminar explored how ML could optimize scientific workflow systems and HPC operations, such as through improved scheduling and fault tolerance. A key focus was on identifying prestigious use cases of ML in HPC and understanding their unique, unmet requirements. The stochastic nature of ML and its impact on the reproducibility of data analysis on HPC systems was also a topic of discussion.

97 MATHEMATICS AND COMPUTING↗