MCNP6.3 Electron Energy Deposition Validation with the Lockwood Experiments
Explore the source record for details and available documents.
Kulesza, Joel A.↗
SEARCH · Engineering Papers
Results for “Lockwood”
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.
Feasibility of a multigroup Boltzmann–Fokker–Planck solution for electron beam dose calculations
Legacy nuclear-reactor Boltzmann solvers start clinical deployment as an alternative to Monte Carlo (MC) codes and Fermi–Eyges semiemprical models in radiation oncology treatment planning. Today’s certified clinical solvers are limited to photon beams. In this paper, ELECTR, a state-of-the-art multigroup electron cross sections generation module in NJOY is presented and validated against Lockwood’s calorimetric measurements, EGS-nrc and GEANT-4 for 1–20 MeV unidirectional electron beams. The nuclear-reactor DRAGON-5 solver is upgraded to access the library and solve the Boltzmann–Fokker–Planck (BFP) equation. A variety of heterogeneous radiotherapy and radiosurgery phantom configurations were used for validation purpose. Case studies include a thorax benchmark, that of a typical breast Intra-Operative Radiotherapy and a high-heterogeneity patient-like benchmark. For all beams, 100% of the water voxels satisfied the American Association of Physicists in Medicine accuracy criterion for a BFP-MC dose error below 2%. At least, 97.0% of adipose, muscle, bone, lung, tumor and breast voxels satisfied the 2% criterion. The average BFP-MC relative error was about 0.56% for all voxels, beams and materials combined. By irradiating homogeneous slabs from Z = 1 (hydrogen) to Z = 99 (einsteinium), we reported performance and defects of the CEPXS mode [US. Sandia National Lab., SAND-89-1685] in ELECTR for the entire periodic table. For all Lockwood’s benchmarks, NJOY-DRAGON dose predictions are within the experimental data precision for 98% of voxels.
42 ENGINEERING↗
A Summary of Validation Studies for the Integrated TIGER Series Performed on ACORN Plus-up 218468/99
The Integrated TIGER Series (ITS) transport code is a valuable tool for photon-electron transport. A seven-problem validation suite exists to make sure that the ITS transport code works as intended. It is important to ensure that data from benchmark problems is correctly compared to simulated data. Additionally, the validation suite did not previously make use of a consistent quantitative metric for comparing experimental and simulated datasets. To this end, the goal of this long-term project was to expand the validation suite both in problem type and in the quality of the error assessment. To accomplish that, the seven validation problems in the suite were examined for potential drawbacks. When a drawback was identified, the problems were ranked based on severity of the drawback and approachability of a solution. We determined that meaningful improvements could be made to the validation suite by improving the analysis for the Lockwood Albedo problem and by introducing the Ross dataset as an eighth problem to the suite. The Lockwood error analysis has been completed and will be integrated in the future. The Ross data is unfinished, but significant progress has been made towards analysis.
97 MATHEMATICS AND COMPUTING↗
Experimental and kinetic modeling study of the low- temperature and high-pressure combustion chemistry of straight chain pentanol isomers: 1-, 2- and 3-Pentanol
Pentanols have received significant attention as a potential alternative fuel or fuel additive owing to their high energy densities and low vapor pressure. The development of robust chemical kinetic models for alternative fuels which can provide accurate and efficient predictions of combustion performance across a wide range of engine relevant conditions is important in developing cleaner, more efficient combustors. Although the high temperature oxidation kinetics of pentanol isomers has been researched considerably, their low temperature combustion chemistry needs further investigation. While previously proposed low temperature mechanisms for 1-pentanol based on analogy and rate rules need further refinement, the low temperature oxidation kinetics of 2-pentanol and 3-pentanol has not been studied previously by any means, experimentally or theoretically. A newly developed kinetic mechanism is presented in this work for the three straight chain pentanol isomers: 1-, 2- and 3-pentanol. Low temperature kinetics is based on a recent study by Lockwood et al., 2022 [20] involving theoretical calculations at the CCSD(T)/cc-pV∞Z level of theory for the oxidation pathways involving alcohol peroxy radicals. Rate of production analyses performed in this study highlight the importance of the newly added pressure-dependent reactions of the α-alcohol peroxy radical forming an RO 2 adduct. While the α-alcohol fuel radical reacts with O 2 to directly decompose via a chemically activated pathway at low pressures, the formation of the RO 2 adduct is favored at high pressures. Furthermore, the detailed model is comprehensively validated against new ignition experiments at low temperature and high pressure, together with the wide range of data available in the literature. Both qualitative and quantitative predictions of the experimental data using the proposed kinetic model are satisfactory for all three pentanol isomers studied here.
1-pentanol↗
Verification and validation testing and tools: comparison between MCNP code versions and nuclear data libraries [Slides]
This presentation discusses the primary goal of software testing which is to test the code for correctness. It also discusses the results for individual suites and the role of validation and verification also referred to in the presentation as V&V. In summation, the V&V framework enables easy comparison between calculations performed with different code versions and/or nuclear data libraries. This entire framework will be distributed with the upcoming MCNP6.3 release. V&V test suites shown and several that were not (Criticality, LAQGSM, Lockwood) will be distributed in the new framework.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Validation of electron backscattering in the MCNP6.3 code
A new electron backscattering validation suite has been developed for the MCNP® code. The calculations in this suite cover both condensed history and single-event transport methods, seven different elemental materials ranging from beryllium to uranium, and incident electron energies from 200 eV to 14.1 MeV. In general, the accuracy of the MCNP code for electron backscattering calculations is energy-dependent. For incident electron energies below 2 keV, the condensed history method is not applied while the single-event method shows poor agreement with experimental measurements due to poor accuracy of the underlying atomic data. For incident energies between 2 keV and 256 keV, the single-event method typically shows better agreement with experiments, while for incident energies above 256 keV the condensed history method gives better accuracy. Additionally, agreement with experimental data is generally worse for low- Z materials such as beryllium and carbon, which may reflect deficiencies in the data or physics underlying electron transport through light elements in the MCNP code. The new electron backscattering validation suite, along with the previously introduced electron stopping power and Lockwood energy deposition validation suites, provide MCNP users with a quantified understanding of the uncertainties in electron transport calculations. Future improvements to electron transport in the MCNP code are necessary for low energies (≤ 2 keV) and for low- Z materials.
36 MATERIALS SCIENCE↗
Architecture and performance of Perlmutter's 35 PB ClusterStor E1000 all-flash file system
NERSC's newest system, Perlmutter, features a 35 PB all-flash Lustre file system built on HPE Cray ClusterStor E1000. Here, we present its architecture, early performance figures, and performance considerations unique to this architecture. We demonstrate the performance of E1000 OSSes through low-level Lustre tests that achieve over 90% of the theoretical bandwidth of the SSDs at the OST and LNet levels. We also show end-to-end performance for both traditional dimensions of I/O performance (peak bulk-synchronous bandwidth) and nonoptimal workloads endemic to production computing (small, incoherent I/Os at random offsets) and compare them to NERSC's previous system, Cori, to illustrate that Perlmutter achieves the performance of a burst buffer and the resilience of a scratch file system. Finally, we discuss performance considerations unique to all-flash Lustre and present ways in which users and HPC facilities can adjust their I/O patterns and operations to make optimal use of such architectures.
97 MATHEMATICS AND COMPUTING↗
Understanding fundamental effects of biofuel structure on ignition and physical fuel properties
The development of biochemicals and biofuels with advantaged properties of higher combustion efficiency at lower emission levels can be aided by a priori prediction of critical global characteristics based on molecular structure. Existing predictive techniques are incomplete for comprehensive ignition and physical fuel property prediction. Here this investigation focuses on filling this knowledge gap by developing a priori prediction techniques for predicting important ignition as well as physical fuel properties of promising advantaged biofuels. We demonstrate the utility of this a priori prediction technique on saturated (≤C 5 ) alcohols, both linear and branched in nature by correlating their molecular structure with both important ignition and physical properties of the fuels. This work utilizes a semi-automated method of predicting ignition characteristics via accurate fundamental-based calculations of the fuel radical cascades that govern low-temperature combustion. This numerical methodology is supported by the relevant fuel ignition delay time (IDT) and research octane number (RON) data experimentally measured by the Advanced Fuel Ignition Delay Analyzer (AFIDA). Phyiscal fuel property trends are also analyzed by correlating the molecular structure of the fuel to its physical characteristics. Based on their molecular structure, recommendations are made on an a priori method for selecting biofuels that have advantaged properties for the selected fuel application. Lastly, the utility of this investigation is extended beyond the subset of alcohols considered by analyzing the IDT trends of i-Pentanol and n-Hexanol, revealing the fidelity of the novel kinetic modelling approach and the accuracy of the predicted trends in this study.
09 BIOMASS FUELS↗
Correlation of optical properties with particle size, morphology, and polymorph of fine- and nano-particle formulations of titanium dioxide powders
Titanium dioxide (TiO 2 ) particulates are known to exhibit different visible and infrared optical properties compared to the bulk material, showing strong dependence on particle size, crystal structure, and morphology. In this study, the optical properties, sizes, and morphologies of TiO 2 particles from two different sources (nano and fine powders) having a) nominally different particle sizes and b) various crystal polymorph mixture fractions are compared using a combination of single particle mass spectrometry, optical spectroscopies, and aerosol characterization methods. The nano sample was found to be largely particles of the anatase polymorph (88% by mass), while the fine sample was found to consist largely of rutile particles (95% by mass). Two distinct particle morphologies (fractal and compact) were found in each powder sample and could be identified and separated in-situ based on particle aerodynamic properties. The attenuation of near-infrared, visible and ultraviolet light by TiO 2 particles shows strong dependence on particle morphology. Furthermore, while the fine particles were found to have larger near-infrared (675–800 nm) extinction coefficients by mass than the nanoparticles, the reverse was true in the ultraviolet and visible regions (370–675 nm). However, for polydisperse particles with different sizes and shapes, the optical behaviors are not straightforward to directly correlate to a combination of physical parameters.
Lockwood, Schuyler P. [Pacific Northwest National ↗
Development and testing of a fission-like neutron field for dosimetry and instrument calibrations based on deuterium-tritium generator
For nearly five decades 252 Cf sources have been used for testing and calibration over a wide range of neutron detection devices used for nuclear safety and radiation protection. However, the vastly increased cost of 252 Cf sources, its short half-life, and concerns about future shortage have prompted nuclear facilities and radiological calibration laboratories that rely on high-intensity neutron sources to seek alternatives. There appear to be no other radionuclide neutron sources that would match “like-to-like” 252 Cf performance characteristics. Alternatives of producing a fission-like neutron spectrum based on commercially available deuterium-tritium (D-T) generators, and without any fissionable materials or beryllium, have been explored. This paper describes efforts of designing the spectrum shaping assembly, its construction and installation. The produced neutron field was evaluated both from the perspective of the neutron fluence spectrum and the resulting dosimetric outcome. Finally, it was tested for practical application by examining response of various health physics instruments and personal monitoring devices in comparison to the response of those devices when exposed to 252 Cf. In conclusion, the response of survey-type instruments tested within the prototype, surrogate fission, reference field are acceptably close to their response in the targeted 252 Cf field.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Advances in predictive chemistry enable a multi-scale rational design approach for biofuels with advantaged properties
Process of building predictive gas-phase chemical mechanisms faster and less computationally intensive using state of the art automated chemical tools and linking the predicted chemical reactivity to underlying chemical structure information.
09 BIOMASS FUELS↗
Photofragment imaging differentiates between one- and two-photon dissociation pathways in MgI+
The bond strength and photodissociation dynamics of MgI+ are determined by a combination of theory, photodissociation spectroscopy, and photofragment velocity map imaging. From 17 000 to 21 500 cm−1, the photodissociation spectrum of MgI+ is broad and unstructured; photofragment images in this region show perpendicular anisotropy, which is consistent with absorption to the repulsive wall of the (1) Ω = 1 or (2) Ω = 1 states followed by direct dissociation to ground state products Mg+ (2S) + I (2P3/2). Analysis of photofragment images taken at photon energies near the threshold gives a bond dissociation energy D0(Mg+-I) = 203.0 ± 1.8 kJ/mol (2.10 ± 0.02 eV; 17 000 ± 150 cm−1). At photon energies of 33 000–41 000 cm−1, exclusively I+ fragments are formed. Over most of this region, the formation of I+ is not energetically allowed via one-photon absorption from the ground state of MgI+. Images show the observed product is due to resonance enhanced two-photon dissociation. The photodissociation spectrum from 33 000 to 38 500 cm−1 shows vibrational structure, giving an average excited state vibrational spacing of 227 cm−1. This is consistent with absorption to the (3) Ω = 0+ state from ν = 0, 1 of the (1) Ω = 0+ ground state; from the (3) Ω = 0+ state, absorption of a second photon results in dissociation to Mg* (3P°J) + I+ (3PJ). From 38 500 to 41 000 cm−1, the spectrum is broad and unstructured. We attribute this region of the spectrum to one-photon dissociation of vibrationally hot MgI+ at low energy and ground state MgI+ at higher energy to form Mg (1S) + I+ (3PJ) products.
Chemistry↗
Design and Implementation of Burst Buffer Over-Subscription Scheme for HPC Storage Systems
Not Available
97 MATHEMATICS AND COMPUTING↗
Optical constants measurements of liquids for modeling aerosolized threats
Spectroscopic identification of aerosolized chemical threats is challenging due to the complex nature of the photon/particle interaction, as well as the diversity of possible particle sizes, morphologies, and compositions. Constructing a database of laboratory - measured transmittance spectra that covers each permutation is not practical. However, calculation of the spectra using the measured optical constants as a function of wavenumber, n(?) and k(?), for each of the liquids and/or solids composing the aerosol particles provides a viable alternative. These synthetic spectra (vis-à-vis laboratory-measured spectra) can be used to identify chemicals of interest and subtract out background interferents in field measurements. Using a well-established multiple pathlength approach, we measure the optical constants n/k for liquids that may be found as chemical species of interest or common background interferents aerosolized in plumes. These measurements are also used to generate several test case aerosol synthetic spectra.
Lockwood, Schuyler↗
Deriving the Infrared Complex Refractive Indices of Organic Powders for Optical Modeling: Comparison of Methods
Identification of solids via infrared reflection spectroscopy requires a spectral library of all solids likely to be encountered. A confounding factor in populating such a spectral library is that the reflectance spectra of solids vary with their form, including particle size, film thickness, and substrate. To reduce the efforts of experimentally constructing such a library, an alternate strategy is to use the wavelength-dependent optical constants, n and k, of a solid to calculate a series of reflectance spectra corresponding to each scenario or morphology. Because most n/k measurements are best performed on mm-sized crystals, however, the challenge of determining the optical constants increases when a solid is only readily available as a powder, as is often the case. Some organic solids, such as caffeine, are both unavailable in large crystals and difficult to press into pellets. In this study, the infrared optical constants, or complex refractive indices, of caffeine were determined using three different methods: single-angle reflectance, infrared spectroscopic ellipsometry, and quantitative absorbance measurements of KBr pellets. The n and k values derived through each method were used to model the hyperspectral imaging reflectance spectrum of a caffeine film on a steel planchet. Over 1,110 – 870 cm-1, the single-angle reflectance-derived n and k had the best correlation with the experimental spectrum. These results suggest different organic solids may require different methods to determine the most accurate infrared complex refractive indices for synthetic spectral libraries.
FTIR, hyperspectral imaging, infrared, infrared (I↗
Extension of Complex Refractive Index Measurements to the Near-Infrared for Liquids: Methodology and Uncertainty Analysis
Optical identification of liquid droplets, aerosols, or thin films is important for many applications. While reference spectra are sometimes available for such measurements, they are not always applicable to the observed spectrum or the given sample morphology. Reference spectra for many forms can be modeled, however, if the n/k vectors (real and imaginary refractive indices) are available. In previous work we have reported protocols to determine the n/k vectors for dozens of liquids, primarily in the mid-infrared (MIR) spectral range from 7500 to 400 cm –1 . In this work we extend the spectral range into the near-infrared (NIR) region, demonstrating a method to measure and merge the data sets to create composite n/k data ranging from 10 000 to 400 cm –1 (1.0 to 25 µm) with absorbance fidelity spanning over four orders of magnitude, and vastly improved signal-to-noise in the NIR. The precision of the composite data is evaluated for three different liquids, focusing primarily on the steps for converting the raw absorbance spectra to k values. The variability in both MIR and NIR data as well as in the final n/k vectors is also investigated for several liquids. For typical liquids, the overall variability (reported as 2σ) in the final n and k-vectors is determined to be ∼0.4% and 3%, respectively. Finally, the derived n/k data are used to calculate absorbance spectra for aerosol droplets, showing marginal variability due to the typical measurement errors in the final n/k vectors.
47 OTHER INSTRUMENTATION↗
Method to derive the infrared complex refractive indices n(λ) and k(λ) for organic solids from KBr pellet absorption measurements
Obtaining the complex refractive index vectors $n($$\tilde{v})$ and $k($$\tilde{v})$ allows calculation of the (infrared) reflectance spectrum that is obtained from a solid in any of its many morphological forms. Here, we report an adaptation to the KBr pellet technique using two gravimetric dilutions to derive quantitative $n($$\tilde{v})$/$k($$\tilde{v})$ for dozens of powders with greater repeatability. The optical constants of bisphenol A and sucrose are compared to those derived by other methods, particularly for powdered materials. The variability of the $k$ values for bisphenol A was examined by 10 individual measurements, showing an average coefficient of variation for $k$ peak heights of 5.6%. Though no established standards exist, the pellet-derived $k$ peak values of bisphenol A differ by 11% and 31% from their single-angle- and ellipsometry-derived values, respectively. These values provide an initial estimate of the precision and accuracy of complex refractive indices that can be derived using this method. Limitations and advantages of the method are discussed, the salient advantage being a more rapid method to derive $n$/$k$ for those species that do not readily form crystals or specular pellets.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Modeling aerosol transmission spectra from n(λ) and k(λ) infrared optical constants measurements of organic liquids and solids
The effects of light scattering and refraction play significantly different roles for aerosols than for bulk materials, making it challenging to identify aerosolized chemicals using traditional spectral methods or spectral reference libraries. Due to a potentially infinite number of particle morphologies, sizes, and compositions, constructing a database of laboratory-measured aerosol spectra is not a practical solution. Here, as an alternative approach, the measured n / k optical vectors of two example organic materials (diethyl phthalate and D-mannitol) are used in combination with particle absorption / scattering theory (Mie theory and FDTD) and the Beer-Lambert law to generate a series of synthetic infrared transmission / scattered light spectra. The synthetic spectra show significant differences versus simple slab transmission spectra, even for small changes in particle size (e.g., 5 vs. 10 µm) for both single particles and ensembles, potentially serving as useful reference data for aerosol sensing. For spherical single particles with diameters of 1 to 10 µm, FDTD simulations predict changes in the magnitudes of spectral shifts and the shapes of the peaks vs. particle size with only small deviations from Mie theory predictions, yet reliably capture the direction of the shifts. Typical spectral peak shifts in the longwave infrared correspond to Δλ ∼0.20 µm (∼34 cm -1 ) when compared to corresponding slab transmission spectra. Additionally, synthetic spectra generated from the n / k values derived using two different methods (KBr pellet transmission and single-angle reflectance) are compared using the Mie theory model.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗