Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “design source code”

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 253 records · Page 14

Investigating resource-efficient neutron/gamma classification ML models targeting eFPGAs

There has been considerable interest and resulting progress in implementing machine learning (ML) models in hardware over the last several years from the particle and nuclear physics communities. A big driver has been the release of the Python package, hls4ml, which has enabled porting models specified and trained using Python ML libraries to register transfer level (RTL) code. So far, the primary end targets have been commercial field-programmable gate arrays (FPGAs) or synthesized custom blocks on application specific integrated circuits (ASICs). However, recent developments in open-source embedded FPGA (eFPGA) frameworks now provide an alternate, more flexible pathway for implementing ML models in hardware. These customized eFPGA fabrics can be integrated as part of an overall chip design. In general, the decision between a fully custom, eFPGA, or commercial FPGA ML implementation will depend on the details of the end-use application. In this work, we explored the parameter space for eFPGA implementations of fully-connected neural network (fcNN) and boosted decision tree (BDT) models using the task of neutron/gamma classification with a specific focus on resource efficiency. We used data collected using an AmBe sealed source incident on Stilbene, which was optically coupled to an OnSemi J-series silicon photomultiplier (SiPM) to generate training and test data for this study. We investigated relevant input features and the effects of bit-resolution and sampling rate as well as trade-offs in hyperparameters for both ML architectures while tracking total resource usage. The performance metric used to track model performance was the calculated neutron efficiency at a gamma leakage of 10 -3 . The results of the study will be used to aid the specification of an eFPGA fabric, which will be integrated as part of a test chip.

47 OTHER INSTRUMENTATION↗

Layered CAD/CSG geometry for spatially complex radiation transport scenarios

Many spatially complex fission, fusion, and national security Monte Carlo (MC) radiation transport scenarios involve combining computer-aided design (CAD) models with constructive solid geometry (CSG) models. A layered geometry method has been implemented in the Shift MC code to address this need. With layered geometry, multiple CAD and/or CSG models can be clipped, translated, rotated, and placed in overlapping layers to form transport-ready geometries. Here, the utility of this method is demonstrated with two problems: (1) a fixed-source simulation with a layered geometry consisting of a LiDAR-generated CAD model of the Combined Arms Collective Training Facility urban environment overlaid with CSG models of a mock hotel and a detector apparatus, and (2) a k-eigenvalue calculation using a layered geometry model of the Transformational Challenge Reactor consisting of CAD fuel elements placed in a CSG core. Tallied particle flux distributions match expectations, but tracking robustness must be improved prior to general-purpose use.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DRiFT - Release 2.1.0: Organic Scintillators and Gas Detectors

DRiFT (a Detector Response Function Toolkit) is LANL-developed software that post-processes output from the extensively validated radiation transport code, MCNP, and generates realistic nuclear instrumentation response. DRiFT is designed to be flexible, enabling users to specify detector type and many experimental settings, as well as accommodating the addition of their own desired features. The focus of this release is on organic scintillator, gas detector, and associated capabilities, while semiconductor features are still under development. Organic scintillators are widely used in the areas of nuclear safeguards and nuclear non-proliferation efforts. DRiFT has several diagnostic and detector physics features relevant to detailed scintillator simulations including: tracking source particle information, scintillation light production, the effects of PMT quantum efficiency and gain, and digitizer settings. Users can select responses from many scintillator and PMT types supported natively by DRiFT, or add their own by following the instructions in this document. DRiFT also has several diagnostic and detector physics features relevant to neutron gas detector simulations including: gas detector wall effects, inactive areas at the end of detector tubes, and effects of a preamplifier. We acknowledge that DRiFT is under active development, bug reports and general questions and comments should be directed to Madison Andrews, madison@lanl.gov. This manual is divided into four parts: I) An overview of DRiFT, including how to obtain and install the executable, II) A description of the detector physics related to scintillators available, III) a description of more general DRiFT features the user may find useful, and IV) a description of the test suite and examples made available with the code release.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Block2 : A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond

block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.

Chemistry↗

Trilinos: Enabling Scientific Computing across Diverse Hardware Architectures at Scale

Trilinos is a community-developed, open-source software framework that facilitates building large-scale, complex, multiscale, multiphysics simulation code bases for scientific and engineering problems. Since the Trilinos framework has undergone substantial changes to support new applications and new hardware architectures, this document is an update to “An Overview of the Trilinos project” by Heroux et al. (ACM Transactions on Mathematical Software, 31(3):397–423, 2005). It describes the design of Trilinos, introduces its new organization in product areas, and highlights established and new features available in Trilinos. Particular focus is put on the modernized software stack based on the Kokkos ecosystem to deliver performance portability across heterogeneous hardware architectures. This article also outlines the organization of the Trilinos community and the contribution model to help onboard interested users and contributors.

Heterogeneous Hardware Architectures↗

Integrating ytopt and libEnsemble to autotune OpenMC

Ytopt is a Python machine-learning-based autotuning software package developed within the ECP PROTEAS-TUNE project. The ytopt software adopts an asynchronous search framework that consists of sampling a small number of input parameter configurations and progressively fitting a surrogate model over the input-output space until exhausting the user-defined maximum number of evaluations or the wall-clock time. libEnsemble is a Python toolkit for coordinating workflows of asynchronous and dynamic ensembles of calculations across massively parallel resources developed within the ECP PETSc/TAO project. libEnsemble helps users take advantage of massively parallel resources to solve design, decision, and inference problems and expands the class of problems that can benefit from increased parallelism. In this paper we present our methodology and framework to integrate ytopt and libEnsemble to take advantage of massively parallel resources to accelerate the autotuning process. Specifically, we focus on using the proposed framework to autotune the ECP ExaSMR application OpenMC, an open source Monte Carlo particle transport code. OpenMC has seven tunable parameters some of which have large ranges such as the number of particles in-flight, which is in the range of 100,000 to 8 million, with its default setting of 1 million. Setting the proper combination of these parameter values to achieve the best performance is extremely time-consuming. Therefore, we apply the proposed framework to autotune the MPI/OpenMP offload version of OpenMC based on a user-defined metric such as the figure of merit (FoM) (particles/s) or energy efficiency energy-delay product (EDP) on Crusher at Oak Ridge Leadership Computing Facility. In conclusion, the experimental results show that we achieve the improvement up to 29.49% in FoM and up to 30.44% in EDP.

Autotuning↗

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - I : Theory

This paper discusses two cross section interpolation methods commonly found in popular nodal codes; the partial derivatives and multiple tables models. The motivation for choosing a model, and thus a case matrix structure, is a trade off between accuracy and computational cost. Due to decades of experience, there are default structures that are sufficient for current light water reactor analysis. However, this is not necessarily the case for advanced reactor designs. Therefore, it is advantageous to understand the sources of error in cross section interpolation models so that the quality of a case matrix may be improved. A mathematical framework for these models is presented in this work that provides a more rigorous connection between the nuclear engineering field's cross section interpolation methods and the broader mathematical field of function approximation. The two cross section models examined in this paper were found to utilize Lagrange interpolation and are a subset of Lagrange tensor products. Classical results of Lagrange polynomial error analysis were then applied to the partial derivative and multiple tables models to derive expressions for the total point-wise error. The analytical results classify the total error into two parts: the model form error and interpolation error. Finally, based on our observations, a better foundation for improving the quality of a case matrix is proposed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance Assurance Planning Guide for Utility Energy Service Contracts: 2025 Edition

Administered by the U.S. Department of Energy's (DOE) Federal Energy Management Program (FEMP), the Utility Program has fostered collaboration among federal agencies and their serving utilities for more than 25 years. The Utility Program supports agencies using Utility Energy Service Contracts (UESCs), a well-developed, effective contracting vehicle that enable the latest approaches to cost-effective energy management at federal sites. Federal agencies have successfully used UESCs to award over 2,000 energy and water efficiency and renewable energy projects, investing approximately $\$$2.8 billion in furthering the Federal Government's efforts to reduce energy intensity. Authorized by 42 U.S. Code section 8256 (10 U.S. Code section 2913 for the Department of Defense), a UESC is a limited-source acquisition between a federal agency and an eligible serving utility for energy management services that generate savings from the implementation of energy- and water -conservation measures (collectively referred to as ECMs), with 42 U.S. Code section 8287 (Defense Federal Acquisition Regulation Supplement, Part 241), providing the term of a UESC, which may extend up to 25 years. Through a UESC, the utility partner assesses designs, and implements the desired ECMs - which can range from lighting retrofits and renewable energy systems, to combined heat and power plants or other technologies and strategies, and may provide financing for the project. The agency may use any combination of appropriations and third-party financing to pay for the project, providing useful flexibility. There is no limit to the project size, big or small, that can be implemented using a UESC. To assist agencies implementing a UESC, FEMP has developed a Utility Energy Service Contract Guide and this companion guidance document to help agencies and their utility partners better understand the best practices for to ensure UESCs continue to perform and generate savings throughout their performance period. These best practices utilize a combination of effective project management, communication, documentation, and a detailed Performance Assurance Plan. This plan is a project specific set of actionable protocols that define important tasks and responsibilities throughout the contract term and reflects the site conditions, complexities, agency capabilities, and operating and maintaining planned ECMs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DRiFT - Release 1.0.0 Organic Scintillators

DRiFT (a Detector Response Function Toolkit) is LANL-developed software that postprocesses output from the extensively validated radiation transport code, MCNP [1], and generates realistic nuclear instrumentation response. DRiFT is designed to be flexible, enabling users to specify detector type and many experimental settings, as well as accommodating the addition of their own desired features. Although DRiFT development has included scintillator [2], gas [3], and semiconductor features [4], the focus of this release is on organic scintillator and associated capabilities. Organic scintillators are widely used in the areas of nuclear safeguards and nuclear non-proliferation efforts [5, 6]. DRiFT has several diagnostic and detector physics features relevant to detailed scintillator simulations including: tracking source particle information, scintillation light production, the effects of PMT quantum efficiency and gain, and digitizer settings. Users can select responses from many scintillator and PMT types supported natively by DRiFT, or add their own by following the instructions in this document. We acknowledge that DRiFT is under active development, bug reports and general questions and comments should be directed to Madison Andrews, madison@lanl.gov. This manual is divided into four parts: I) An overview of DRiFT, including how to obtain and install the executable, II) A description of the detector physics related to scintillators available, III) a description of more general DRiFT features the user may find useful, and IV) a description of the test suite and examples made available with the code release.

61 RADIATION PROTECTION AND DOSIMETRY↗

A nodal transport method for coupled fast-thermal reactor analysis

We report that for efficient core analyses of coupled fast-thermal reactors, a new deterministic method has been developed based on the variational nodal transport method of VARIANT. A new multigroup cross section generation procedure was devised by combining Monte Carlo lattice calculations for thermal assemblies and the two-step procedure of the MC 2 -3 code for fast assemblies. To reduce assembly homogenization errors in nodal transport calculations, a new nodal equivalence method was developed based on the partial current discontinuity factor (PCDF) and incorporated in VARIANT. A practical procedure to calculate PCDFs with fixed source supercell calculations was also devised. The performance of the proposed method was investigated using a test problem derived from the versatile coupled test reactor (VCTR) design. The new procedure produced multigroup cross sections accurately for all assemblies. VARIANT transport calculations with PCDFs produced accurate multiplication factor and power distribution compared to reference Serpent-2 Monte Carlo solutions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Code-verification techniques for hypersonic reacting flows in thermochemical nonequilibrium

The study of hypersonic flows and their underlying aerothermochemical reactions is particularly important in the design and analysis of vehicles exiting and reentering Earth's atmosphere. Computational physics codes can be employed to simulate these phenomena; however, code verification of these codes is necessary to certify their credibility. To date, few approaches have been presented for verifying codes that simulate hypersonic flows, especially flows reacting in thermochemical nonequilibrium. In this work, we present our code-verification techniques for verifying the spatial accuracy and thermochemical source term in hypersonic reacting flows in thermochemical nonequilibrium. Additionally, we demonstrate the effectiveness of these techniques on the Sandia Parallel Aerodynamics and Reentry Code (SPARC).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improvements to PVWatts for Fixed and One-Axis Tracking Systems

This work presents improvements to the widely used NREL PVWatts photovoltaic system energy model to improve modeling accuracy for typical fixed and one axis system designs. The aim is to calculate losses in the PV system assuming typical modern system design practices, while maintaining simplicity by keeping the required set of input parameters small. These improvements allow users to more credibly and quickly evaluate competing system designs in early stage feasibility. Common submodels for module cover, spectral, snow, tracker, transformer, plant controller, and self-shading losses, in addition to a bifacial gain option, are incorporated into the PVWatts model, and are shown to improve PVWatts' system performance prediction capabilities without major impact to ease of use. We anticipate including these improvements in a future release of NREL's open source PVWatts code, and some of the features may become available in the System Advisor Model (SAM) desktop software as well as the popular PVWatts web application.

14 SOLAR ENERGY↗

Virtual Tokamak for Test and Development of Plasma Control Applied to NSTX-U

Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural Network Analysis of Nuclear Magnetic Resonance and Infrared Spectra

Nuclear magnetic resonance (NMR) spectroscopy and infrared (IR) spectroscopy are powerful chemical characterization techniques with broad general usage. However, the manual evaluation of the resulting spectra is time-consuming and requires significant expertise, preventing insights from being used in real-time applications. With recent advances in computation and artificial intelligence (AI), new tools are available for automating spectral interpretation. In this work, machine learning (ML) algorithms using 1-dimensional convolutional neural networks (CNNs) were applied to identify common functional groups from spectral information. Raw spectra were collected virtually from the Human Metabolome Database (HMDB) and National Institute of Standards and Technology (NIST) Chemistry WebBook and processed into a suitable standard. Algorithm design was tailored to best fit the nature of the problem, with built-in flexibility to accommodate relevant parameters beyond the raw spectral input, specifically solvent identity and magnetic frequency for NMR. The predictive capability of the algorithm in identifying functional groups is displayed in several examples. This methodology has been compiled into a code repository and could easily be modified to adapt alternative data sources, including other spectrum types. To mitigate overfitting, a common problem in mathematical modeling where overfamiliarity with training data produces trends that are not representative of the general data, a novel metric was developed, referred to as Accufit. Accufit includes a parameter that penalizes substantial differences in the training accuracy and the accuracy of an independent validation set. Examples are presented showing the effectiveness of Accufit in maintaining the model’s predictive capability while controlling the overfitting when used as a custom metric for hyperparameter tuning.

Sturgill, James↗

Application of Monte Carlo code SHIFT for MSBR dose rate calculations

This paper presents a study that applied the Monte Carlo code SHIFT to calculate the radiological environment within a detailed Molten Salt Breeder Reactor (MSBR) model. It represents one of the first applications of the SHIFT code in stand-alone simulations for non-light water reactors to help demonstrate the code's potential uses in the design, licensing, and operation of advanced reactors. The radiological conditions of the MSBR were modeled when the reactor is at two different operation modes: normal full power and a drained state. The Forward Weighted-Consistent Adjoint Driven Importance Sampling (FW-CADIS) hybrid method in SHIFT was applied successfully to calculate the ex-core neutron and gamma dose rates for the MSBR at full power. Neutron and gamma dose rates within the drain cell were also calculated for the MSBR at the drained state by integrating the source terms obtained from an ORIGEN-S depletion calculation into the SHIFT simulation. The results indicate that in the drained state, the delayed gammas from the depleted fuel salt are the main contributors to the dose rates. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Density reconstruction in convergent high-energy-density systems using x-ray radiography and Bayesian inference

X-ray radiography is a technique frequently used to diagnose convergent high-energy-density (HED) systems, such as inertial confinement fusion implosion, and provides unique information that is not available through self-emission measurements. In this work, we investigate the scope and limits of that information using a radiography simulation combined with the Bayesian inference workflow. The accuracy of density reconstruction from simulated radiographs of spherical implosions driven with 27 kJ laser energy is assessed, including the increase or decrease in accuracy due to the addition of Lagrangian marker layers, Poisson noise, and improved prior information. This work is the first to present the full uncertainty distributions inferred from radiography analysis in HED systems and demonstrates the importance of constructing the full posterior probability density, as opposed to a point estimate, due to the modal structure of the likelihood surface introduced by typical experimental noise sources. This general methodology can be used both for robust analysis of radiographic data and for an improved design of radiography experiments by modeling the full experimental system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design

GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗