Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Error Incidents”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Progress Report of Alloy 709 Performance in Liquid Sodium

With improved high-temperature properties, Alloy 709 shows great potential for advanced reactor applications. To understand the effects of sodium on A709, small tensile specimens were machined from two commercial heats and subject to sodium exposure in two forced convection loops operated by Argonne. Thermal aging experiments were also conducted in parallel to isolate any sodium-specific effects. An increasing trend in yield strength and a decrease trend in ultimate tensile strength can be seen with the increasing exposure time, suggesting microstructural changes at these temperatures. Nonetheless, the thermally aged and sodium exposed samples behaved similarly, implying no sodium-specific effect at these temperatures. In contrast, 316H SS exposed to sodium at 650°C showed a strong sodium effect. The different responses of A709 and 316H SS to sodium exposure is of interest and requires additional effort to understand. The cold traps of the sodium loop were replaced early in FY23. During a maintenance operation, an unexpected sodium leak incident occurred at one of the sodium loops, triggering an emergency response. An investigation into the incident revealed that the leak happened at a VCR fitting on the recently replaced cold trap. The metal gasket of the leaked VCR fitting was not compressed sufficiently. A dimensional mismatch between the new component and old piping configuration was responsible for a significant misalignment. A new cold trap with improved design is in production. The new design provides additional flexibility and an access point for post-installation helium leak check and for venting of trapped gas. These design features will make the system more error-tolerant and the installation quality verifiable. A thermodynamic analysis was performed to understand the carburization-decarburization behavior of A709 in sodium environment. The equilibrium carbon activity in the alloy was evaluated with M 23 C 6 , TiC, and NbC. The carburization-decarburization boundary was calculated with different methods developed previously. The results suggests that A709 will be carburized in SMT-1 and SMT-2, and can be decarburized in low-carbon-activity environments such as EBR-II at temperatures above 620-650°C. Nonetheless, A709 should be less likely to be decarburized than 316H SS when exposed to the same low-carbon-activity environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Comparison of active interrogation methods for source location in a scattering and absorbing medium, consisting of PGNAA, and an AmBe quasi-forward biased directional source

Source location of Special Nuclear Material (SNM) encompassing 95% 235 U and 239 Pu is identified by utilizing a directional source from Patent No: US20190013109A1 using Prompt Gamma Neutron Activation Analysis (PGNAA) and neutron spectroscopy simulated with Monte-Carlo N-Particle transport 6.2 (MCNP). BC-408, HPGe, LaBr 3 detector arrays were used to identify the location of the SNM using total counts incident on each detector, and PGNAA photopeaks from HPGe and LaBr 3 detector arrays in a polyethylene shield. The conducted simulations varied the volume and location of the SNM in the MCNP input files to observe how the source location method behaved. PGNAA photopeaks used for source identification include 61 keV from fission, 2.223 MeV prompt gamma from hydrogen, 511 keV annihilation, and a single and double escape peaks from the prompt gamma interaction from hydrogen. The capabilities of each detector systems to acquire well resolved photopeaks with a 1% relative error or less, and total relative error for F4 and F8 tallies were less than 0.015% relative error. Source predictions of the SNM with uneven amounts of polyethylene shielding between the source and detectors was observed to overpredict and give invalid source location predictions. Source locations of the SNM with even amounts of polyethylene material between the source and each detector were found to be valid. With a 1 Ci 241 Am source activity, it was determined that 1630 s were needed to obtain the results for each detector system with the quasi-forward directional AmBe source. Coupling source and material identification together would increase acquisition time but would only require one system to determine.

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

Advanced Photovoltaic Module Characterization: Using Image Transformers for Current–Voltage Curve Prediction From Electroluminescence Images

Individual photovoltaic (PV) module health monitoring can be a daunting task for operation and maintenance of solar farms. Modules can be inspected through luminescence, thermal imaging, and current–voltage (I–V) curve analyzes for identification of damage and power loss. I–V curves provide easily interpretable data to determine module health as they directly provide electrical performance metrics. However, in order to obtain these curves, modules must be disconnected from the array and either removed to a solar simulator or characterized in situ with corrections for module temperature, the incident solar spectrum, and intensity. Luminescence or thermal images of a module are relatively easy to acquire in situ. Electroluminescence (EL) images highlight physical defects in the modules but do not provide easily interpretable features to correlate with electrical performance. This work presents a SWin transformer network to predict I–V curves for PV modules from their corresponding EL images. The predicted I–V curves allow the accurate prediction of the maximum power point (MPP), short-circuit current I sc , and open-circuit voltage V oc with a mean error less of than 1%. Comparing single diode model (SDM) parameters extracted from the predicted curves to those extracted from the true curves, the series resistance R s demonstrates a mean error of 5.19%, and the photocurrent I a mean error of 0.197%. The shunt resistance R sh and dark current Io parameters are predicted with larger errors because of their sensitivity to small changes in the I–V curve.

Byford, Brandon K. [New Mexico State Univ., Las Cr↗
DOE OSTI2022

Pulsed Sphere Time of Flight Neutron Spectrum Measurements: Blank Run and 1.8 Mean-Free-Path thick Polyethylene Moderating Sphere

From the late 1960s to 1985, Lawrence Livermore National Laboratory (LLNL) undertook a series of ‘pulsed sphere’ experiments for 32 materials involving 148 different experiments using 75 different spheres to measure the neutron leakage spectra for different materials from 14-MeV neutrons generated by 3 H(d, n) 4 He reactions induced by an incident D+ beam from the Insulated Core Transformer (ICT) accelerator at LLNL. This deuteron beam was focused to impinge on tritium loaded onto a titanium substrate within a low mass target assembly placed within the center of a spherical shell constructed from the material of interest. The neutron time of flight was measured at a detector outside a collimator. The program was intended to be comprehensive with numerous targets of various thicknesses spanning the periodic table from H to 239 Pu. Five different neutron detectors were utilized to avoid possible systematic errors due to differences in detector responses. Further, these measurements were designed with the aim that the simple geometry of the measurements could be easily simulated to validate Monte Carlo transport codes and nuclear data libraries. As a result, these benchmarks have been extensively used in the validation of ENDL, ENDF and other libraries. Two experiments selected as benchmark cases in this evaluation are: 1) target assembly without material sphere (blank run), and 2) target assembly with a hollow polyethylene sphere, 1.8 mean free path (mfp) or 16.56 cm thick.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A Comparison of Infectious Disease Forecasting Methods across Locations, Diseases, and Time

Accurate infectious disease forecasting can inform efforts to prevent outbreaks and mitigate adverse impacts. This study compares the performance of statistical, machine learning (ML), and deep learning (DL) approaches in forecasting infectious disease incidences across different countries and time intervals. We forecasted three diverse diseases: campylobacteriosis, typhoid, and Q-fever, using a wide variety of features (n = 46) from public datasets, e.g., landscape, climate, and socioeconomic factors. We compared autoregressive statistical models to two tree-based ML models (extreme gradient boosted trees [XGB] and random forest [RF]) and two DL models (multi-layer perceptron and encoder–decoder model). The disease models were trained on data from seven different countries at the region-level between 2009–2017. Forecasting performance of all models was assessed using mean absolute error, root mean square error, and Poisson deviance across Australia, Israel, and the United States for the months of January through August of 2018. The overall model results were compared across diseases as well as various data splits, including country, regions with highest and lowest cases, and the forecasted months out (i.e., nowcasting, short-term, and long-term forecasting). Overall, the XGB models performed the best for all diseases and, in general, tree-based ML models performed the best when looking at data splits. There were a few instances where the statistical or DL models had minutely smaller error metrics for specific subsets of typhoid, which is a disease with very low case counts. Feature importance per disease was measured by using four tree-based ML models (i.e., XGB and RF with and without region name as a feature). The most important feature groups included previous case counts, region name, population counts and density, mortality causes of neonatal to under 5 years of age, sanitation factors, and elevation. This study demonstrates the power of ML approaches to incorporate a wide range of factors to forecast various diseases, regardless of location, more accurately than traditional statistical approaches.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Fluorescence strobo-stereoscopy for specular reflection-suppressed full field of view imaging

This report introduces fluorescence strobo-stereoscopy (FSS) to suppress strong specular reflection and enable the full field of view (FFOV) 3D surface imaging while the part is rotating. Specular reflection off the target surface significantly degrades the image quality and becomes critical for highly reflective surface measurements. In FSS, the fluorescent dye-doped fluid applied on the machined surface is excited upon incident ultra-violet light and becomes a new light source by Stokes' Law. Thus, specular reflection off of smooth surface can be suppressed by separating the fluorescent light from the excitation light. The developed FSS comprises a pair of imaging cameras, spatial filters, and an excitation light source. As a result, FSS effectively rejected the specular reflection and improved the FFOV 3D surface image quality of the machined part by enhancing contrast in the rotating target surface. Such enhancements in 3D imaging allowed to identify manufacturing tolerance of the part and to detect the surface features. The axial and lateral accuracy errors of FSS were 2.3% and 1.4% with the target size of 4.07 mm and 0.215 mm, respectively. A whole view reconstruction of the cylindrical target sample was performed, and the corresponding cylindricity and diameter deviation were assessed. The fluid media effect and the target surface quality effect were discussed.

3D surface imaging↗
DOE OSTI2020

First Trials with EMPHASIS

This report presents a series of three EMPHASIS trial problems we ran in March 2019, to demonstrate (1) that we can set up and run an SGEMP calculation, and (2) that we understand the nature and units of the input specifications and the output results. As these are the first EMPHASIS problems we have set up and run, we present some of the input decks and some of the lessons learned, to better orient a reader who might also be a new user of EMPHASIS. The EMPHASIS manual and the ITS manual are both good, but not complete enough to set up and run these applications. Substantial trial and error was still required. The trial problems are all based on an evacuated aluminum cylinder with 1 m inner radius and 1 m inner height. The cylinder is 100 mil thick. We take the axis of symmetry to be the z axis. The aluminum cylinder will be treated as perfectly conducting. The trial problems are (1) the cavity driven with a short pulse on an antenna at the center of the cavity, (2) the cavity driven by a prescribed short pulse of photocurrent injected normally from the bottom inside of the cylinder, and (3) the cavity driven by a prescribed short pulse of x-rays normally incident on the bottom outside of the cylinder.

42 ENGINEERING↗
DOE OSTI2022

Measurement of Compton Scattering at MAMI for the Extraction of the Electric and Magnetic Polarizabilities of the Proton

A precise measurement of the differential cross sections dσ/dΩ and the linearly polarized photon beam asymmetry Σ 3 for Compton scattering on the proton below pion threshold has been performed with a tagged photon beam and almost 4π detector at the Mainz Microtron. The incident photons were produced by the recently upgraded Glasgow-Mainz photon tagging facility and impinged on a cryogenic liquid hydrogen target, with the scattered photons detected in the Crystal Ball/TAPS setup. Using the highest statistics Compton scattering data ever measured on the proton along with two effective field theories (both covariant baryon and heavy-baryon) and one fixed-t dispersion relation model, constraining the fits with the Baldin sum rule, we have obtained the proton electric and magnetic polarizabilities with unprecedented precision: α E1 = 10.99 ± 0.16 ± 0.47 ± 0.17 ± 0.34, β M1 = 3.14 ± 0.21 ± 0.24 ± 0.20 ± 0.35; in units of 10 -4 fm 3 where the errors are statistical, systematic, spin polarizability dependent, and model dependent.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Accounting for electron-beam-induced warping of molecular nanocrystals in MicroED structure determination

High-energy electrons induce sample damage and motion at the nanoscale to fundamentally limit the determination of molecular structures by electron diffraction. Using a fast event-based electron counting (EBEC) detector, we characterize beam-induced, dynamic, molecular crystal lattice reorientations (BIRs). These changes are sufficiently large to bring reciprocal lattice points entirely in or out of intersection with the sphere of reflection, occur as early events in the decay of diffracted signal due to radiolytic damage, and coincide with beam-induced migrations of crystal bend contours within the same fluence regime and at the same illuminated location on a crystal. These effects are observed in crystals of biotin, a series of amino acid metal chelates, and a six-residue peptide, suggesting that incident electrons inevitably warp molecular lattices. The precise orientation changes experienced by a given microcrystal are unpredictable but are measurable by indexing individual diffraction patterns during beam-induced decay. Reorientations can often tilt a crystal lattice several degrees away from its initial position before irradiation, and for an especially beam-sensitive Zn(II)-methionine chelate, are associated with dramatic crystal quakes prior to 1 e − Å −2 electron beam fluence accumulates. Since BIR coincides with the early stages of beam-induced damage, it echoes the beam-induced motion observed in single-particle cryoEM. As with motion correction for cryoEM imaging experiments, accounting for BIR-induced errors during data processing could improve the accuracy of MicroED data.

Vlahakis, Niko (ORCID:0000000250920265)↗
DOE OSTI2020

Deterministic incident-wave elevation prediction in intermediate water depth

Potential performance gains from optimal (non-causal) impedance-matching control of wave energy devices in irregular ocean waves are dependent on deterministic wave elevation prediction techniques that work well in practical applications. Although a number of devices are designed for operation in intermediate water depths, little work has been reported on deterministic wave prediction in such depths. Here, we investigated a deterministic wave-prediction technique based on an approximate propagation model that leads to an analytical formulation, which may be convenient to implement in practice. To improve accuracy, an approach to combine predictions based on multiple up-wave measurement points is evaluated. The overall method is tested using experimental time-series measurements recorded in the U.S. Navy MASK basin in Carderock, MD, USA. For comparison, an alternative prediction approach based on Fourier coefficients is also tested with the same data. Comparison of prediction approaches with direct measurements suggest room for improvement. Possible sources of error including tank reflections are estimated, and potential mitigation approaches are discussed.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

A new route to obtain fluorescence X-ray absorption spectra of compounds and to remove the self-absorption induced nonlinearity in the spectra

A new route to obtain fluorescence X-ray absorption spectra of compounds and to remove the self-absorption induced nonlinearity in the spectra is described. The fluorescent intensity If is linearly proportional to the absorption coefficient μ. For studies of surface structures around an element (κ) the fluorescence detection is often the mode of choice. However, the measurement may suffer from a self-absorption (SA) effect which nonlinearly distorts the spectra. The effect is severe when κ is concentrated or the measurements are carried out in certain geometries. Here, the correlations among emission events in compounds are examined following resonance X-ray core-electron excitation within κ. Under conditions leading to SA, If emitted from κ apparently has a conjugated relationship with the fluorescent intensities simultaneously emitted from other elements (ξ). Normalizing the former (κ) by the latter (ξ) will largely remove SA effects and reduce this nonlinear problem to a tractable linear problem. This does result in a moderate reduction of the spectral amplitude due to the so-called secondary emission from ξ excited by the emission from κ. Nonetheless, the resulting spectra will allow one to accurately determine bond distances and disorder and, in some respects, can be superior to spectra obtained via the absorption channel. For μξ < μκ and grazing incidence geometry, the amplitude reduction can be small and simple normalization is sufficient to restore the spectral integrity with remarkable accuracy. This has been instrumental in unravelling the surface and subsurface structures around cations in amorphous Ga–In–O and Zn–Sn–O films which are otherwise inaccessible due to severe SA effects. This method has also been applied to several samples with μξ ≃ μκ to examine its applicability. For these samples, the amplitude reduction is 12 ± 4% versus their standards for the data measured with the classical 45°/45° geometry. This experimental method is easy to implement. Since If from κ and ξ are measured by the same detector system, it is also superior to other methods in removing systematic errors such as detector system nonlinearity, electronic noise, and some beam instabilities, and in removing spectral imperfections due to, for example, SA effects, diffraction effects and sample inhomogeneity. The distortions resulting from the latter can be severe in the spectra measured in transmission mode.

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

Lifetime measurements of states of S 35 , S 36 , S 37 , and S 38 using the AGATA γ -ray tracking spectrometer

Lifetimes or lifetime limits of a small number of excited states of the sulfur isotopes with mass numbers A = 35 , 36, 37, and 38 have been measured using the differential recoil-distance method. The isotopes of sulfur were populated in binary grazing reactions initiated by a beam of 36 S ions of energy 225 MeV incident on a thin 208 Pb target which was mounted in the Cologne plunger apparatus. The combination of the PRISMA magnetic spectrometer and an early implementation of the AGATA γ-ray tracking array was used to detect γ rays in coincidence with projectile-like nuclear species. Lifetime measurements of populated states were measured within the range from about 1 to 100 ps. The number of states for which lifetime measurements or lifetime limits were possible was limited by statistics. For 35 S, the lifetime was determined for the first 1/2 + state at 1572 keV; the result is compared with a previous published lifetime value. The lifetime of the 3 - state of 36 S at 4193 keV was determined and compared with earlier measurements. No previous lifetime information exists for the (6 + ) state at 6690 keV; a lifetime measurement with large associated error was made in the present work. For 37 S, the states for which lifetime limits were established were those at 646 keV with J π = 3/2 - and at 2776 keV with J π = 11/2 - ; there are no previously published lifetime values for excited states of 37 S. Finally, a lifetime limit was established for the J π =(6 + ) state of 38 S at 3675 keV; no lifetime information exists for this state in the literature. Measured lifetime values were compared with the results of state-of-the-art shell-model calculations based on the PSDPF, SDPF-U, and FSU effective interactions. In addition, nuclear magnetic-dipole and electric-quadrupole moments, branching ratios, mixing ratios, and electromagnetic transition rates, where available, have been compared with shell-model values. The current work suffers from poor statistics; nevertheless, lifetime values and limits have been possible, allowing a useful discussion of the ability of state-of-the-art shell-model calculations to reproduce the experimental results.

electromagnetic transitions↗
DOE OSTI2022

Field validation of data-driven BSDF and peak extraction models for light-scattering fabric shades

Shading and daylighting systems affect cooling, heating, and lighting energy use by modulating solar radiation through the building façade. Characterizing shading systems holistically and accurately helps designers and engineers evaluate shading systems to achieve energy and non-energy performance goals. These complex fenestration systems can be modeled using Bidirectional Scattering Distribution Functions (BSDF), which map incident radiation to hemispherical distributions of outgoing radiation. Data-driven, tabulated BSDFs are derived from interpolated goniophotometer measured data, then sampled during the raytracing calculation. A peak extraction (PE) algorithm was developed to circumvent limits in BSDF angular resolution, where the specular peak is extracted during simulation by evaluating the BSDF in the through direction and surrounding region. The objective of this study was to validate this measurement and modeling workflow using field monitored data from a full scale testbed with eleven installed fabrics of different weaves, openness factors, and colors and assess the accuracy of the workflow under different adaptation and contrast conditions. Test conditions were limited to clear sky conditions with the sun in the field of view. Results showed that, for tensor tree datasets, vertical illuminance, solar luminance (2.5° apex), and daylight glare probability (DGP) were predicted to within a mean bias error (MBE) error of -456 lx (-12.3%), -3.46e5 (-38.4%), and -0.042 (-7.8%) when full PE occurred. With a binary classification of glare/ no glare, DGP was predicted accurately with a true positive rate of 0.98 and true negative rate of 1.0 using tensor tree data and less accurately with Klems BSDF data, particularly for cases of no glare. The workflow may be of insufficient accuracy to distinguish borderline performance between fabrics using the four-point glare scale, particularly under low adaptation, high contrast daylit conditions. Errors were due to reductions in peak shape and intensity across the BSDF interpolation and data reduction workflow. Future work is needed to better preserve measurement fidelity during interpolation and sampling, which in turn will improve PE performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Measurement of laser absorption in underdense plasmas using near-field imaging of the incident and transmitted beams

Measurements of laser absorption in high-temperature, underdense plasmas produced at the Omega Laser Facility are made using two near-field imaging detectors that diagnose the spatial profile and energy of the port P9 beam before and after it transmits through the plasma. The incident beam is sampled using a partial reflection from a full-aperture, (30 cm-diam) uncoated wedge pickoff located before the target chamber vacuum window and final focus lens assembly. A concave mirror reduces the reflected beam size, allowing it to be recorded directly using a charged-coupled device (CCD) camera. The P9 transmitted beam diagnostic (P9TBD) characterizes the transmitted light by terminating the expanded beam on a semi-transparent diffuser and imaging the illuminated surface using a lens and CCD camera. The P9TBD samples a numerical aperture twice as large as the input beam, allowing the energy of transmitted beams with moderate levels of beam spray to be measured accurately. Calibration shots with no plasma provide a path to infer absorption without absolute photometric calibration of either detector. Lastly, the cross-calibration between the two detectors was measured to remain stable at ±200 ppm, enabling measurements of total beam absorption below 1% with ±0.07% error.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Optimal design of acoustic metamaterial cloaks under uncertainty

In this work, we consider the problem of optimal design of an acoustic cloak under uncertainty and develop scalable approximation and optimization methods to solve this problem. The design variable is taken as an infinite-dimensional spatially-varying field that represents the material property, while an additive infinite-dimensional random field represents the variability of the material property or the manufacturing error. Discretization of this optimal design problem results in high-dimensional design variables and uncertain parameters. To solve this problem, we develop a computational approach based on a Taylor approximation and an approximate Newton method for optimization, which is based on a Hessian derived at the mean of the random field. We show our approach is scalable with respect to the dimension of both the design variables and uncertain parameters, in the sense that the necessary number of acoustic wave propagations is essentially independent of these dimensions, for numerical experiments with up to one million design variables and half a million uncertain parameters. Additionally, we demonstrate that, using our computational approach, an optimal design of the acoustic cloak that is robust to material uncertainty is achieved in a tractable manner. The optimal design under uncertainty problem is posed and solved for the classical circular obstacle surrounded by a ring-shaped cloaking region, subjected to both a single-direction single-frequency incident wave and multiple-direction multiple-frequency incident waves. Finally, we apply the method to a deterministic large-scale optimal cloaking problem with complex geometry, to demonstrate that the approximate Newton method’s Hessian computation is viable for large, complex problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Exploring the dependence of gas cooling and heating functions on the incident radiation field with machine learning

ABSTRACT Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (shap) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behaviour of cooling and heating functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Emulation of radiation transport in 3D stochastic media using 1D planar Monte Carlo stochastic media radiation transport algorithms

A subset of stochastic media radiation transport problems involves those in which radiation is incident on a thin slab of stochastic material. Particle tracking in 3D for such problems is expensive, and 1D planar models lack accuracy because they only allow the material to change in one dimension. Therefore, we propose dimensional emulation, which through a slight modification allows existing 1D planar geometry stochastic media radiation transport models to reproduce results from the equivalent 3D models by allowing the material to change in all three dimensions, reproducing the fidelity of the 3D model for the low computational cost of the 1D planar model. In this work, we apply dimensional emulation to three Monte Carlo stochastic media radiation transport models: Chord Length Sampling (CLS), the Local Realization Preserving method (LRP), and a variant of Conditional Point Sampling (CoPS). For a common Markovian benchmark set, the 3D emulation variants of these algorithms are numerically verified to reproduce the results of the 3D variants within statistics while running 1.3 to 2 times faster in the implementation within Sandia National Laboratories open-source research code PlaybookMC. The 3D emulation variants are also shown to yield a 72%–92% reduction in error for the thin slab problems in comparison to the 1D benchmark. As a result, the 3D emulation variant of CLS and CoPS-1 are shown to reproduce 3D CLS results that were used to approximate results for a 3D spherical inclusion geometry benchmark set.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Shedding light on the 239Pu fission source term with new high-precision experiments and advanced fission modeling

In the last decade, there has been a renaissance of fission research resulting in new high-precision experiments and advanced fission modeling. For instance, the Chi-Nu and CEA teams supplied, for the first time, the 239 Pu prompt fission neutron spectrum (PFNS) for broad ranges of incident and outgoing neutron energies. The CEA team also measured 239 Pu average prompt neutron multiplicities, ν ̄ p , with lower statistical uncertainties and a technique significantly different than the one used in the past. The NIFFTE collaboration provided 239 Pu( n ,f)/ 235 U( n ,f) cross section shape ratios with uncertainties below 1% utilizing a novel detector type. Advanced fission event generators were developed, among them CGMF, FIFRELIN, FREYA, and GEF, which calculate post-scission fission observables in a correlated manner. These new experimental data and more consistent fission models change the evaluated PFNS, ν ̄ p , and ( n ,f) cross sections, some only modestly, compared to ENDF/B-VIII.0. In turn, the individual new nuclear data distinctly change simulated effective neutron multiplication factors of fast critical assemblies, but their combined impact is small, while affecting the prediction of LLNL pulsed sphere neutron leakage spectra and reaction rates only within experimental uncertainties. Also, the parameters obtained from fitting to ν ̄ p reproduce various post-scission fission observables within the uncertainties of experimental data. This indicates that new differential experiments and consistent fission modeling reduce compensating errors present in ENDF/B-VIII.0.

fission cross section↗