SEARCH · Engineering Papers

Results for “surface codes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A Python Script to Read MCNP6.3 Surface-Source Files

This report provides a Python script to read an MCNP ® surface source file created with the SSW card with SYM = 0 (i.e., the default symmetry treatment). For background: the general format of an MCNP surface-source file is described in; however, that document did not provide coding and/or a tool to interrogate such files. The current format will not be given in this document other than through the record-read statements necessary for the script to function. The reader capability in this report is augmented with the ability to directly write a couple demonstrative outputs: 1. A comma-separated value (CSV) file containing particle phase-space state information and 2. A Matplotlib histogram of the energy distribution of the particles. This report also describes accompanying verification work that shows the script performing as required with MCNP6.2, MCNP6.3, and (expected) MCNP6.4 surface-source files. However, users of the enclosed script must still verify that the script is behaving correctly for their own work.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Spatial Study 2021: Sample-Based Surface Water Chemistry and Organic Matter Characterization across Watersheds in the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin. The dataset provides geochemistry and organic matter characterization data generated from samples collected during the same two-week period at 47 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. Related sensor data are published at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1892054. This data package was originally published September 2022. It was updated May 2023 (modified files) and November 2024 (new and modified files). See the change history section in the readme for more details. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) dissolved inorganic carbon (DIC), dissolved organic carbon (DOC; reported as non-purgeable organic carbon; NPOC), total nitrogen (TN), total suspended solids (TSS), ions, and benzene polycarboxylic acid (BPCA) concentration and stable isotope data; (5) averaged values from water chemistry data; (6) surface water sampling protocol; (7) sensor protocol (8) readme; (9) methods codes; (10) international generic-sample number (IGSN) mapping file; and (11) folder of high resolution characterization of organic matter via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). This folder contains two subfolders, one containing the .xml data files and the other containing instructions for using Formultitude (https://github.com/PNNL-Comp-Mass-Spec/Formultitude) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .R, .ref, or .xml. We thank the United States Forest Service, Washington Department of Natural Resources, Washington Department of Fish and Wildlife, Washington State Parks, Confederated Tribes and Bands of the Yakama Nation, and Cowiche Canyon Conservancy for access to field locations where these samples were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Carbon Organisms Rhizosphere and Protection in Soil Environment model script and input data for soil moisture-respiration responses in tropical forests

Objectives: Climatic drying is predicted for many tropical forests, yet models remain poorly parameterized for tropical forests, hampering predictions of forest-climate feedbacks. We applied an integrated model–experiment approach, parameterizing an ecosystem model Carbon Organisms Rhizosphere and Protection in the Soil Environment (CORPSE) with tropical forest observational data, and comparing model predictions with a field drying manipulation. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We used the field data to parameterize and run tests in the model.Results: Measured CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. We used this data to parameterize the model, which then predicted increased soil CO2 fluxes in wetter and fertile forests with drying, and decreased fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment in the forests initially suppressed soil CO2 fluxes across forests, with sustained suppression after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season), as predicted by the model. The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Code files:CORPSE_array.py: Defines the equations of the CORPSE modelCORPSE_solvers: Functions for running the CORPSE model using either iterative or ordinary differential equation (ODE) solversrun_Panama_sims.py: Read in datasets and run the model simulations for this studyInput data:PanamaGradientEcosystemChem_BT_CPools_20152016CO2_DC_20190615.xlsx: Plot characteristics used in running model simulationsLiCor compiled surface flux only to 2020_03 DC_20200825.xlsx: Surface gas exchange fluxes used in model-data comparisonsPARCHED litterfall data for Ben Sulman LD 20200902.xlsx: Litterfall data used to drive model simulationsInitialization data:state_500y_20190823.csv: Initial state of model pools based on previous spinup runsOutput data:Outputs/prev_moisture_response.csv: Simulations of multiple sites using original model moisture response function.Outputs/updated_moisture_response.csv: Simulations of multiple sites using updated model moisture response function.Outputs/dry15_prev_moisture_response.csv: Simulations with soil moisture reduced by 15%, using original moisture response function.Outputs/dry15_updated_moisture_response.csv: Simulations with soil moisture reduced by 15%, using updated moisture response function.Outputs/dry30_prev_moisture_response.csv: Simulations with soil moisture reduced by 30%, using original moisture response function.Outputs/dry30_updated_moisture_response.csv: Simulations with soil moisture reduced by 30%, using updated moisture response function.Outputs/latestart_prev_moisture_response.csv: Simulations with extended dry season, using original moisture response function.Outputs/latestart_updated_moisture_response.csv: Simulations with extended dry season, using updated moisture response function.Outputs/[site name]_oneyear.csv: One-year simulation for each site in expanded site list using original moisture response function.Outputs/[site name]_oneyear_dried.csv: One-year simulation for each site in expanded site list using original moisture response function, with soil moisture reduced by 25%.Outputs/[site name]_oneyear_updated_moisture_response.csv: One-year simulation for each site in expanded site list using updated moisture response function.Outputs/[site name]_oneyear_updated_moisture_response_dried.csv: One-year simulation for each site in expanded site list using updated moisture response function, with soil moisture reduced by 25%.Field plot location data:There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Space weather observations, modeling, and alerts in support of human exploration of Mars

Space weather observations and modeling at Mars have begun but they must be significantly increased to support the future of Human Exploration on the Red Planet. A comprehensive space weather understanding of a planet without a global magnetosphere and a thin atmosphere is very different from our situation at Earth so there is substantial fundamental research remaining. It is expected that the development of suitable models will lead to a comprehensive operational Mars space weather alert (MSWA) system that would provide rapid dissemination of information to Earth controllers, astronauts in transit, and those in the exploration zone (EZ) on the surface by producing alerts that are delivered rapidly and are actionable. To illustrate the importance of such a system, we use a magnetohydrodynamic code to model an extreme Carrington-type coronal mass ejection (CME) event at Mars. The results show a significant induced surface field of nearly 3,000 nT on the dayside that could radically affect unprotected electrical systems that would dramatically impact human survival on Mars. Other associated problems include coronal mass ejection (CME) shock-driven acceleration of solar energetic particles producing large doses of ionizing radiation at the Martian surface. In summary, along with working more closely with international partners, the next Heliophysics Decadal Survey must include a new initiative to meet expected demands for space weather forecasting in support of humans living and working on the surface of Mars. It will require significant effort to coordinate NASA and the international community contributions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Virtual prototyping of liquid metal blanket performance in fusion pilot plant

Liquid metal blanket is a dominant design option for the next step fusion devices responsible for harvesting energy from fusion reaction, and simultaneously producing fuel for the same reaction through tritium breeding. Liquid metal blankets introduce additional complexity to the design due to fluid motion, fluid structure interaction, and magnetohydrodynamic (MHD) effects arising from the motion of the conducting fluid through the magnetic field. They are also directly affected by the plasma heat flux and neutronic fluence. PPPL is currently developing a virtual prototyping system for numerical analysis of the liquid metal blankets for future fusion devices. The system has a customized 3D computational fluid dynamics (CFD) code in its core, allowing MHD flow and conjugate heat transfer analysis in blankets fluids and solids. The code was successfully used before for dual coolant blanket analysis [A. Khodak et al., Fusion Eng. and Des. 137 (2018)]. Recently the same code was modified to allow verified simulation of MHD flows at high Hartmann numbers of several thousand typical for blanket applications. CFD code receives volumetric heat source distribution from the neutronic analysis based on MCNP code. In addition, direct tritium breeding simulation will be performed allowing optimization of the blanket performance. 2D axisymmetric version of neutronics code will be used for rapid optimization, with 3D version employed for detailed analysis. The surface heat distribution on the plasma facing wall will be defined by the software HEAT allowing 3D modeling of the heat flux based on the magnetic field distribution including gyro-orbit effects. Results of thermal analysis are imported into structural analysis code also included in the system. Finally, direct import of CAD geometry will be used for analyzing all components and as a result design option can be efficiently optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗

DOE OSTI2022

WHONDRS Surface Water and Sediment Non-Purgeable Organic Carbon and FTICR-MS across Glacial Features in Svalbard 2021 (v3)

This dataset supports a broader study examining land-atmosphere gas exchange associated with permafrost and glaciers in Svalbard. The dataset provides geochemistry and organic matter characterization data generated from surface water and sediment. Samples were collected from 8 spatially diverse sites across central Svalbard and also from a more comprehensive temporal and spatial study within one glacial catchment area throughout one summer melt season. Related data were collected and will be published separately in collaboration with Yde and Kleber.This data package was originally published in September 2022. It was updated in March 2023 (v2; new and modified files), and July 2025 (v3; modified files). See the change history section in the readme for details.For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) dissolved organic carbon (DOC; measure as non-purgeable organic carbon; NPOC); (5) surface water sampling protocol; (6) sediment extraction protocol; (7) readme; (8) methods codes; (9) international generic sample number (IGSN) mapping file; and (10) folder of high resolution characterization of organic matter via 7 and 12 Tesla (7T and 12T) Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) generated through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). The FTICR folder contains two subfolders containing the 7T and 12T .xml data files and another subfolder containing instructions for using Formularity (https://omics.pnl.gov/software/formularity) and an R script to process the data based on the user's specific needs. All files are .csv, .pdf, .R, .ref, or .xml.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

ParFlow Sand Tank: A tool for groundwater exploration

The ParFlow Sand Tank model is an open source application designed to allow users to interactively simulate and visualize groundwater movement through the subsurface. The app is designed for both research and education; teaching hydrogeology concepts and making it easy explore and run sophisticated groundwater simulations. Our goal is to support increased accessibility and usability of research grade hydrology tools for research and teaching. The Sand Tank application simulates groundwater and surface water fluxes as well as contaminant transport in real time using the integrated physical hydrology model ParFlow (Kollet & Maxwell, 2006; Maxwell & Miller, 2005; Osei-Kuffuor et al., 2014) and the particle tracking code EcoSlim (Maxwell et al., 2019). ParFlow is a numerical hydrology model that simulates spatially distributed groundwater and surface water flow. It is a well established research tool with more than 90 publications documenting its development use to advance our understanding of groundwater dynamics and groundwater surface water interactions from the hillslope to the continental scale e.g. (Condon et al., 2020; Condon & Maxwell, 2019; Maxwell & Condon, 2016). It is designed for efficient parallel computation and has been run on many platforms spanning from laptops to supercomputers. However, one of the challenges of ParFlow is that it requires significant training and hydrologic expertise to develop simulations. The Sand Tank application makes this model accessible to anyone for education and exploration. Our application uses ParFlow for its simulation backend and ParaView for the data loading and processing. The communication infrastructure relies on the ParaViewWeb framework. We use model templates deployed in Docker images to setup the Sand Tank framework. Users can build the application locally or interact with it through our web deployment. When interacting with a template users can interactively change model parameters like subsurface processes or pump/inject water into the subsurface and watch the system respond to their changes in real time as the simulation runs. Additionally, our template setup will allow more advanced users to build custom templates of increasing complexity for both research and educational purposes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Cubit for MCNP Unstructured Mesh Analysis of Oktavian Benchmarks

The Monte Carlo N-Particle (MCNP) transport code developed by Los Alamos National Laboratory (LANL) can be used to transport various particles across user defined three dimensional (3D) geometries. Traditionally, these geometries are created as constructive solid geometry (CSG), involving the use of Boolean operators on defined surfaces to create 3D regions known as cells. A newer method of geometry definition in the MCNP code is unstructured mesh (UM) embedded within a CSG cell using the universe and fill repeated-structure features. An MCNP UM calculation requires an MCNP input file and mesh geometry file. The MCNP code cannot be used to generate UM geometry models. A computer-aided design (CAD) software is typically used to construct a solid geometry model, and a CAD file is then imported into a mesh generation software to prepare and generate a mesh model. Some mesh generation software packages have the ability to create a solid geometry and thus a separate CAD software for creating a CAD model is not needed. In this work, Cubit, a geometry creation and meshing software developed by Sandia National Laboratories, is used to construct UM models for MCNP simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Assessment of RELAP5-3D Code with Molten Salt Heat Transfer Experiments

The accurate thermal-hydraulic assessment of thermal storage systems is crucial for addressing the integrity and performance of the storage system design. This is particularly important for thermal storage systems using molten salts as a heat transport and storage medium, as the chemical corrosion and erosion by the high-temperature molten salts can cause serious damage to the structural components of the storage system. To understand these effects for TerraPower’s Natrium® Demonstration Reactor design, the system thermal-hydraulic analysis code, RELAP5-3D, is utilized to analyze the thermal storage system. This study evaluates two heat transfer correlations implemented in the RELAP5-3D code—the Dittus-Boelter and Gnielinski correlations—using selected benchmark cases from two molten salt heat transfer experiments conducted at Xi’an Jiao Tong University and the German Aerospace Center. The Nusselt numbers calculated using the Gnielinski correlation agree with the experimental data within ±5% of the relative deviations at 300°C and 400°C. With the increasing Reynolds numbers, the Dittus-Boelter correlation underestimates the Nusselt numbers by up to 22% compared to the Gnielinski correlation. Comparison of the RELAP5-3D calculation results with experimental data at a high temperature of 550°C revealed that the thermo-chemical characteristics of solar salt at high temperatures surpassing the decomposition temperatures of nitrate salts reduce the accuracy of heat transfer model of the RELAP5-3D code. Specifically, it has been observed that the temperatures of the wall surface and liquid near the wall can surpass the decomposition temperatures, even though the bulk temperature of molten salts remains significantly lower than these decomposition temperatures. In summary, this study recommends the use of the Gnielinski correlation for the heat transfer analysis of molten salt, rather than the Dittus-Boelter correlation.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2024

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Interactions involved in the adsorption of ethylene glycol and 2-hydroxyethoxide on the Au(111) surface: a Density Functional Theory study

Abstract Context:The monolayers of ethylene glycol and 2-hydroxyethoxide on gold surfaces have been used in hybrid materials as biosensors. In this article, the adsorption of ethylene glycol and 2-hydroxyethoxide on the Au(111) surface was analyzed. For the first system, ethylene glycol on Au(111), there are Au$$\cdot \cdot \cdot $$ · · · O and Au$$\cdot \cdot \cdot $$ · · · H interactions. To the best of our knowledge, the Au$$\cdot \cdot \cdot $$ · · · H interaction has been overlooked until now. However, in this work, there is strong evidence that this interaction is important to stabilize the system. For the second system, the atomic interactions mentioned previously are also predicted, although there is an additional interaction between 2-hydroxyethoxide molecules. Such an interaction induces the link -O-H-O-, with high values of the electron density at the critical points of the corresponding bond path of the O-H interaction. These links suggest the forming of ethylene glycol chains.Methods:The calculations were performed using two exchange-correlation functionals: BEEF-vdW and C09$$_{x}$$ x -vdW; both functionals incorporate dispersion effects within the Kohn-Sham approach in Density Functional Theory as implemented in GPAW code and ASE computational packages. The contacts between the molecules considered in this article and the Au(111) surface were analyzed through the Quantum Theory of Atoms in Molecules implemented in GPUAM code.

Biochemistry & Molecular Biology↗

DOE OSTI2022

Gamma-Ray Imaging to Monitor Accumulating Holdup

Over time, nuclear materials can accumulate in process equipment and keeping track of such material is of interest to many stakeholder communities including safeguards, process monitoring, criticality safety, and radiological safety. Passive gamma-ray imaging provides a visual and quantitative tool that can be applied to this problem, allowing one to determine the amount of material present from the images. In this work we emulate the accumulation of highly enriched uranium (HEU) in a large rectangular stainless-steel duct by stacking 5, 27 × 50 cm2 cards, each of which has approximately 11 g of 235U in oxide form distributed across its surface, behind each other. A series of images were collected using a coded-aperture gamma-ray imager, starting with a single card, adding a second card, etc. Long integrations were collected to allow evaluating the performance of the approach as a function of integration time. The images are analyzed using first principles to determine the amount of material in each image. The results are compared to the known amount of material in each configuration, and as a function of time. The analysis includes a first-principles error propagation and to further evaluate the uncertainty, a boot-strap sampling approach is also applied to the data. The work offers a proof of concept for unattended monitoring of uranium holdup with a quantitative gamma-ray imaging technique.

Ziock, Klaus-Peter↗

DOE OSTI2022

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

hippynn Python Package

hippynn is a python package for defining, training, and applying neural networks to atomistic systems. In particular, it focuses on Hierarchical Interacting Particle Neural Networks (HIP-NNs), a deep learning architecture for atomistic systems. HIP-NNs take input data describing the properties of atomistic systems (often generated using ab-initio quantum mechanics) to learn fast and accurate models. A trained HIP-NN can predict potential energy surfaces, atomic charges, and more. hippynn uses PyTorch for portable high-performance code, including both CPU and GPU support. hippynn allows for extensive customization, including user-defined models and loss functions, to facilitate future research into extensions of HIP-NN as well as other atomistic deep learning models.

Lubbers, Nicolas↗

DOE OSTI2022

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Simulation of plasma and neutral transport in PISCES-RF using SOLPS-ITER

In this research, the fluid plasma transport code SOLPS-ITER is applied and validated against experimental data from the plasma interaction surface component experimental station (PISCES)-RF linear plasma device to establish a physics basis for plasma and neutral transport in its two magnetic field (B-field) geometry setups-(1)the cusp and (2) non-cusp or linear B-field. The main focus of this study is to understand (1) radial plasma transport (2) heat and particle loads on the upstream dump and downstream target plate, and (3) the physics of plasma-neutral interactions in PISCES-RF. The simulation setup adheres to typical PISCES-RF experimental conditions, with a 2D helicon power deposition profile as an input heating source. SOLPS-ITER simulations reproduce experimental conditions with the Bohm diffusion model for both B-field configurations of the PISCES-RF experiment. Major energy loss channels include neutral radiation and power deposited on the wall and dump plate, with only 1% of the input power reaching the target. The ionization front is well confined near the dump plate due to the heating and puffing regions. Additionally, SOLPS-ITER simulation results are also found to be in very good agreement with the particle-in-cell calculations using the code-PICOS++ which supports the validity of SOLPS in low collisionality regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗