Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,495 records · Page 83

The Preliminary Verification of the 2-D Transport Module of Radiation Shielding Code ARES - Paper 33

ARES is a multi-group SN particle transport code system with arbitrary order anisotropic scattering. It is primarily used for the nuclear systems' shielding calculation and analysis in the case of one dimensional, two dimensional and three dimensional systems. ARES2D is the 2-D transport module of ARES code system. It is used to solve the first-order form of the steady-state Boltzmann transport equation in the two dimensional geometry with uniform rectangular mesh. The equation's energy dependence is modeled using multiple energy groups. The angular dependence is modeled using discrete ordinates method. The spatial variable is modeled using finite differencing scheme. The overall solution technique is source iteration with Partial Current Rebalance acceleration of the scattering source. ARES2D is capable of calculating k eigenvalue and fixed source problem with arbitrary order anisotropic scattering in two dimensional geometry. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Consideration for the application of artificial neural networks to criticality code validation - 14270

The artificial neural network (ANN) is a consistent example for non-parametric parameter estimation and regression, with powerful universal approximation properties, that is distribution independent. ANN application to criticality code validation centers on aligning mathematical relationships between code validation parameters and neural network prediction for those parameter values. Two specific goals are: (1) an appropriate ANN architecture to determine validation parameter values, and (2) bootstrapping individual ANNs as a 'Prediction Interval' (PI) technique to estimate the Calculation Margin and Upper Subcritical Limit. The individual ANN workflow is segregated into Training, Validation (network), and Testing Sets, where the network predictive properties and performance are expressed in the latter set. The Bootstrap ANN replicates a large set of individual ANNs (taken from the original data distribution) where (1-α/2)100% PI estimates can be determined. Results from the ANN method are compared to a variety of validation examples: a 'toy example' from ANSI/ANS-8.24, Appendix D, a variety of methods noted in a Nuclear Criticality Safety Division 2009 paper, and several validation campaigns for highly enriched uranium and plutonium systems. The ANN method is robust and can be applied to single or multivariate inputs. A trained ANN can provide both an average estimate (e.g., single value) and an estimate over a range of independent variable values where the kUSL can be formed in a regression relationship suitable for extrapolation. The input dependence has been tested, ranging from 15 to over 400 sample size distributions. An expanded Shapiro-Wilk test (for sample sizes up to 5000) was implemented to confirm that the ANN method is transparent to distribution function. The Principle Component Analysis method was applied as an input preprocessing step embedded within the ANN in order to reduce dimensionality and correlation for multivariable inputs. This can potentially aid in determining accurate margins over applying arbitrary margins. The automatic optimization of the Bootstrap ANN PI is left for future consideration. (authors)

42 ENGINEERING↗

The international research reactor conversion effort and its contribution to the validation of reactor physics and thermal-hydraulics codes

Full text of publication follows. With about 200 reactors in operation, civilian Research and Test Reactors (RTRs) represent approximately a third of the global nuclear fleet. RTRs make use of core designs that are drastically different from commercial power plants to perform a wide variety of non-power applications that greatly benefit society. Because they often rely on high neutron fluxes, RTRs are designed with relatively compact cores and as a result, prior to the 1980's, were often deployed using Highly Enriched Uranium fuel (HEU, {sup 235}U/U = 20 wt. %). Due to proliferation concerns, the international community aims at eliminating the use of HEU in civilian facilities and favor instead the use of Low Enriched Uranium fuel (LEU, {sup 235}U/U < 20 wt. %). A program to support conversion of the world's RTRs to LEU fuel has been initiated in 1978 by the U.S. Department of Energy (DOE). This program is still alive today and has achieved more than 103 conversion metrics. Today, the program focuses heavily on the conversion of so-called high-performance RTRs, which are far more challenging than previous conversions as they require new fuel element designs and the use of new, higher density LEU fuel forms. Development and qualification of new LEU fuel elements is ongoing and requires extensive engineering analysis and testing. Both activities require the development and validation of codes and methods for reactor physics, thermal-hydraulics, and multi-physics, which in turn rely on experiments performed in RTRs or other experimental facilities. This talk will present the collection of RTR experimental data and benchmark analyses from the international conversion program that contribute to the validation of computer codes and methods.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Walnut v1.2 Code Abstract

Abstract for publication for copyright assertion for Walnut v1.2 software code to satisfy prime contract requirements. HQ approval to assert copyright on this code was granted to NNSS on 9/30/2024

47 OTHER INSTRUMENTATION↗

ASME Design Code Rule Changes for Nuclear Graphite

The American Society of Mechanical Engineers Boiler Pressure and Vessel Code (ASME BPVC) Section III, Division 5, Article HHA-3000 outlines graphite core component and graphite core assembly design guidelines. Graphite core components are defined as ?components manufactured from graphite that are installed to form a graphite core assembly within the reactor pressure vessel of a high temperature, graphite moderated fission reactor.? (p. 413) Graphites? inherent defect distributions do not allow for deterministic material reliability. Rather, graphite has variable strength distributions which change by grade. Article HHA-3000 outlines two semi-probabilistic methods, the full and simplified assessments, which set design load limit targets for each of three component structural reliability classes. The Design Task Group was officially recognized as a specialized task group within ASME November of 2023, though we?ve been collaborating since 2022. The purpose of the Design Task Group is to correct, clarify, and make HHA-3000 function as intended. The Design Task Group will sunset once we?ve achieved our objectives. The Design Task Group was specifically told to not write new Code. While there may be more precise and more accurate methods to determine reliability targets, the current methods are conservative, relatively simple to implement, and have thus far been considered satisfactory for setting design reliability targets. Much of the ground-work to write proposal files and background documents for records to make the changes needed to achieve our objective have been completed. The Design Task Group has documented much of their work through papers, presentations, and memorandums. Three memorandums in which INL team members had substantial contributions are found in the Appendices: FEA Modeling for the Baseline Program, Evaluating the Effects on Margin of Updating the Threshold and Shape Parameters in the Full Assessment, and Interpretations of the Full and Simplified Assessments in ASME BPVC. Most of the on-going work to achieve the Design Task Group?s objective will be addressing comments on existing records and moving records through the balloting process. The Design Task Group met bi-weekly mostly through the end of FY2023. Since February 2024, the Design Task Group has mostly been completed with solving and documenting the technical issues associated with the assessments. Unless new tasks are identified, the remaining work of the Design Task Group will be political and editorial.

97 MATHEMATICS AND COMPUTING↗

Reductive Analysis with Compiler-Guided Large Language Models for Input-Centric Code Optimizations

Input-centric program optimization aims to optimize code by considering the relations between program inputs and program behaviors. Despite its promise, a long-standing barrier for its adoption is the difficulty of automatically identifying critical features of complex inputs. This paper introduces a novel technique, reductive analysis through compiler-guided Large Language Models (LLMs), to solve the problem through a synergy between compilers and LLMs. It uses a reductive approach to overcome the scalability and other limitations of LLMs in program code analysis. The solution, for the first time, automates the identification of critical input features without heavy instrumentation or profiling, cutting the time needed for input identification by 44× (or 450× for local LLMs), reduced from 9.6 hours to 13 minutes (with remote LLMs) or 77 seconds (with local LLMs) on average, making input characterization possible to be integrated into the workflow of program compilations. Optimizations on those identified input features show similar or even better results than those identified by previous profiling-based methods, leading to optimizations that yield 92.6% accuracy in selecting the appropriate adaptive OpenMP parallelization decisions, and 20-30% performance improvement of serverless computing while reducing resource usage by 50-60%.

Input-Centric Optimization↗

Revisiting Source Convergence Diagnostics in the KENO Monte Carlo Neutron Transport Codes [Abstract]

Monte Carlo criticality transport codes, which rely on the power iteration procedure, are a fundamental tool for nuclear criticality safety practitioners in assessing the neutron multiplication factor (k eff ) for problems involving fissile material. In these calculations, ensuring the convergence of both the fission source distributions and the k eff estimate for accurate results is crucial. However, a converged k eff estimate does not necessarily mean the fission source distribution is also converged because the fission source and flux distribution may continue to evolve even after k eff convergence. Therefore, most Monte Carlo transport criticality codes now offer various diagnostic tests to assess fission source convergence in addition to the k eff convergence by analyzing the trends of these quantities over multiple generations.

AZURE↗

Multiphysics Demonstration of Temperature-Driven Assembly Bowing in SFRs using MOOSE-Based Codes

Core bowing is an important passive safety mechanism in liquid metal cooled fast reactors. When the core restraint system is properly designed, temperature and flux gradients influence assemblies in the core to bow into less reactive configurations during accident scenarios, resulting in negative reactivity feedback. Prediction of core bowing involves complex interplay of radiation transport, impacts of fluid flow and heat transfer on duct temperature, and mechanical responses to the induced temperature and flux gradients. Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program [1], an integrated multiphysics approach is being developed to model the core bowing phenomena in liquid metal-cooled fast reactors with the Multiphysics Object Oriented Simulation Environment (MOOSE) [2]. In this methodology, the MOOSE-based reactor physics code Griffin [3] will solve the neutron transport equation and determine the power distribution. With the detailed power distribution from Griffin, the subchannel analysis codes MOOSE-Subchannel [4] and Pronghorn [5] are utilized to calculate the assembly temperature distribution. MOOSE’s Solid Mechanics [6] and Contact [7] Modules are leveraged to calculate the thermal expansion and duct bowing displacement with the duct wall temperature from thermal hydraulics calculation. In this work, an initial one-way coupling demonstration of the integrated multiphysics approach has been performed on a seven-assembly problem based on the sodium-cooled fast reactor ABR-1000 design [8]. The neutronics calculation with Griffin is not yet involved in the current simulation. MOOSE-Subchannel and Pronghorn evaluate fluid and solid temperature based on a fixed power distribution. In addition, one-way coupling is utilized in this coupled calculation, via Pronghorn passing the duct temperature data to the MOOSE Solid Mechanics calculation. An assessment of the Solid Mechanics module was performed in parallel to verify duct bowing behavior with duct-to-duct contact phenomenon [9]. The displacement from MOOSE Solid Mechanics is not yet transferred back and utilized in the Pronghorn and MOOSE-Subchannel calculation. This model will be available on the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) repository [10]. Future stages of this work will involve solving problems of increasing complexity as well as adding more physics (e.g. reactor physics) to the integrated workflow to reach the end goal of modeling the core bowing phenomenon with an integrated multiphysics workflow.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Sockeye Code Validation Against UNIST Heat Pipe - Poster - Internship 2025

Sockeye is an engineering level code being developed under the MOOSE framework for modeling heat pipes. A heat pipe is a sealed tube with a wick structure filled with a working fluid that transfers heat efficiently and passively via phase change of the working fluid. Heat applied to the evaporator end creates a pressure gradient which causes vapor to migrate to the condenser end where the vapor deposits its energy and condenses. Capillary force generated by the wick structure then draws the liquid back to the evaporator. A validation study was performed to compare the Sockeye code against data produced at the Ulsan National Institute of Standards and Technology for a slightly overfilled sodium heat pipe. The heat pipe was operated under natural convection cooling and a decreasing inactive length was observed. The experiment was modeled in Sockeye using the condenser pool model, the front non-condensable gas (NCG) model, and the mixture NCG model separately to reproduce the inactive length phenomenon. Also, a new capability was implemented in MOOSE to allow for conjugate heat transfer from the condenser based on the Churchill-Chu correlation for natural convection. The condenser pool model showed an increasing inactive length, demonstrating that the behavior was likely not caused by a pool of excess liquid. The front NCG model was able to show good agreement with the experiment, but it suffered convergence issues at the front. Finally, the mixture NCG model gave good results when the axial mesh was sufficiently refined. This work culminated in additions to the Sockeye documentation and a contribution to a journal article that will be published later. This poster is a summary of my work which I can take back to my university for presentation.

42 - ENGINEERING↗

COR-0011 Rev 6 MARVEL Project Code of Record

This Code of Record identifies the codes, standards, and procedures necessary to design, develop, construct, and startup the Microreactor Applications Research Validation and Evaluation (MARVEL) Project at the Materials and Fuel Complex (MFC) Transient Reactor Test (TREAT) Facility and the Idaho National Laboratory (INL). The MARVEL Project is an INL test microreactor funded by the United States Department of Energy (DOE) via the Microreactor Program (MRP). The goal of the project is to establish an operational nuclear applications test bed that can generate combined heat and power to enable integration and R&D with end-user technologies, as well as allow microreactor technologists to test next-generation control systems. The microreactor is a thermal reactor utilizing Uranium Zirconium Hydride (UZrH) fuel with review and authorization by the Department of Energy Idaho Operations Office (DOE-ID) for National Environmental Policy Act (NEPA) compliance, safety review, and supplemental readiness assessments for startup and operation. To enable rapid deployment, the MARVEL reactor will reside in the Transient Reactor Test (TREAT) Facility and utilize the existing operating Category B reactor facility, approved facility safety basis, operating crews, and recent re-start experience.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Neural Network Code MATLAB

SAND2022-12606 O Neural Network Code MATLAB reads in an excel file, processes the data, trains a neural network, and tests the neural network. The code generates predictions of diffusion constants based on the data in the excel file. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Allers, Joshua↗

Study of components and statistical reaction mechanism in simulation of nuclear process for optimized production of {sup 64}Cu and {sup 67}Ga medical radioisotopes using TALYS, EMPIRE and LISE++ nuclear reaction and evaporation codes

Production of medical radioisotopes is one of the most important tasks in the field of nuclear technology. These radioactive isotopes are mainly produced through variety nuclear process. In this research, excitation functions and nuclear reaction mechanisms are studied for simulation of production of these radioisotopes in the TALYS, EMPIRE and LISE++ reaction codes, then parameters and different models of nuclear level density as one of the most important components in statistical reaction models are adjusted for optimum production of desired radioactive yields.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interfacing VPSC with finite element codes. Demonstration of irradiation growth simulation in a cladding tube

This Milestone report shows good progress in interfacing VPSC with the FE codes ABAQUS and MOOSE, to perform component-level simulations of irradiation-induced deformation in Zirconium alloys. In this preliminary application, we have performed an irradiation growth simulation in the quarter geometry of a cladding tube. We have benchmarked VPSC-ABAQUS and VPSC-MOOSE predictions with VPSC-SA predictions to verify the accuracy of the VPSCFE interface. Predictions from the FE simulations are in general agreement with VPSC-SA simulations and also with experimental trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Code Benchmark of Depressurized Conduction Cooldown Transient in the High Temperature Test Facility

This paper presents results from modeling of a depressurized conduction cooldown (DCC) transient at the High Temperature Test Facility (HTTF) as part of the OECD-NEA Thermal Hydraulics Code Validation Benchmark for High-Temperature Gas-Cooled Reactors using HTTF Data . This paper briefly describes the benchmark and the models being used. It then presents a comparison of steady state and transient results based on the Problem 2 Exercise 1A and 1B definitions. We compare block and helium temperature distributions, mass flow distribution, and energy balance in steady state. All models show comparable mass flow distributions and energy balances. The temperatures within the core and outer regions are comparable in all models too, but inner reflector temperatures can vary significantly. Despite that, we find that the models are in good agreement for the full-power steady state. In the DCC, we look at block temperature at the core midplane and RCCS water exit temperature. The INL and ANL models are found to be in excellent agreement with one another on block temperature over time, while the agreement when the KAERI and NRG models are added into consideration is good. Differences in the transient heat removal from the RCCS cause the differences in block temperature over time in these models. The CNL models show similar trends to the INL, ANL, KAERI, and NRG models, but the temperatures are high because the volumes used in calculating the average temperature include the heater rods in the CNL models only. The HUN-REN model shows results that suggest significantly lower heat removal in the RCCS which merit further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computer code for double beta decay QRPA based calculations

The computer code developed by our group some years ago for the evaluation of nuclear matrix elements, within the QRPA and PQRPA nuclear structure models, involved in neutrino-nucleus reactions, muon capture and β{sup ±} processes, is extended to include also the nuclear double beta decay.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Plasma transport in an Eulerian AMR code

A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions to flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unit of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Visualizing Hierarchical Performance Profiles of Parallel Codes Using CallFlow

Calling context trees (CCTs) couple performance metrics with call paths, helping understand the execution and performance of parallel programs. To identify performance bottlenecks, programmers and performance analysts visually explore CCTs to form and validate hypotheses regarding degraded performance. However, due to the complexity of parallel programs, existing visual representations do not scale to applications running on a large number of processors. In this work, we present CallFlow, an interactive visual analysis tool that provides a high-level overview of CCTs together with semantic refinement operations to progressively explore CCTs. Using a flow-based metaphor, we visualize a CCT by treating execution time as a resource spent during the call chain, and demonstrate the effectiveness of our design with case studies on large-scale, production simulation codes.

97 MATHEMATICS AND COMPUTING↗

Building America Guidance for Identifying and Overcoming Code, Standard, and Rating Method Barriers

This guidance document was prepared using the input from the meeting summarized in the draft CSI Roadmap to provide Building America research teams and partners with specific information and approaches to identifying and overcoming potential barriers to Building America innovations arising in and/or stemming from codes, standards, and rating methods.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗