Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code validation”

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 577 records · Page 32

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗

Engineering-Scale Integrated Energy System Data Projection Demonstration via the Dynamic Energy Transport and Integration Laboratory

The objective of this study is to demonstrate and validate the Dynamic Energy Transport and Integration Laboratory (DETAIL) preliminary scaling analysis using Modelica language system-code Dymola. The DETAIL preliminary scaling analysis includes a multisystem integral scaling package between thermal-storage and hydrogen-electrolysis systems. To construct the system of scaled equations, dynamical system scaling (DSS) was applied to all governing laws and closure relations associated with the selected integral system. The existing Dymola thermal-energy distribution system (TEDS) facility and high-temperature steam electrolysis (HTSE) facility models in the Idaho National Laboratory HYBRID repository were used to simulate a test case and a corresponding scaled case for integrated system HYBRID demonstration and validation. The DSS projected data based on the test-case simulations and determined scaling ratios were generated and compared with scaled case simulations. The preliminary scaling analysis performance was evaluated, and scaling distortions were investigated based on data magnitude, sequence, and similarity. The results indicated a necessity to change the normalization method for thermal storage generating optimal operating conditions of 261 kW power and mass flow rate of 6.42 kg/s and the possibility of reselecting governing laws for hydrogen electrolysis to improve scaling predictive properties. To enhance system-scaling similarity for TEDS and HTSE, the requirement for scaling validation via physical-facility demonstration was identified.

08 HYDROGEN↗

Testing the Activation Analysis for Fusion in OpenMC

OpenMC is a community-developed Monte Carlo neutron and photon transport simulation code. It can perform fission simulations such as fixed-source, k-eigenvalue, and subcritical multiplication calculations on models built using either a constructive solid geometry or CAD representation. To explore the use of OpenMC for fusion activation analysis, a detailed model of the Fusion Neutronics Science Facility (FNSF) was first developed for comparisons against an existing SERPENT model. A 90-degree model of FNSF in Standard-Triangle-Language (STL) CAD format was converted to Constructive Solid Geometry (CSG) using each code's built-in functions, and the geometries were validated by ensuring no cells overlapped and no particles were lost during simulations. The neutron fluxes were calculated and compared for multiple components close to the plasma. The results show differences mostly below 1% in fluxes and averaged 8% for activity and decay heat. Here, the work described in this study tests the CAD-based geometry using the DagMC toolkit in OpenMC and compares the activation analysis of OpenMC to SERPENT code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sierra/SD – Verification Test Manual – 5.22

Verification and validation (V&V) of scientific computing programs are important at Sandia National Labs due to the expanding role of computational simulation in managing the United States nuclear stockpile. The complexities of structural response calculations used to analyze physical problems, the varieties of codes applied to the calculations, and the importance of accurate predictions when assessing field conditions demand confidence in the consistency and accuracy of computer codes. Confidence in the accuracy of the predictions arising from computer simulations must ultimately be gained through verification and validation. The Sierra salinas structural dynamics analysis code, Sierra/SD, is used at the DOE Laboratories, and in several DOD projects. The roles of Sierra/SD in the qualification of weapon systems and components for normal and hostile environments throughout the Stockpile-to-Target Sequence include to, • Redesign weapon components. • Certify weapon components and systems for target environments such as hypersonic vehicles. • Certify that components will survive the thermal mechanical shock loads associated with hostile environments. • Evaluate current stockpile issues, including issues associated with uncertainty quantification. • Address many other problems that are encountered in stockpile management. The Sierra/SD verification plan is described, and an evolving set of key verification tests are described in detail. The verification tests ensure the correctness of the mathematics and numerical algorithms associated with functionality describing engineering phenomena. Development is in accordance with a set of tailored Software Quality Engineering (SQE) practices. SQE practices guide the overall verification and validation effort.

97 MATHEMATICS AND COMPUTING↗

Autodifferentiable Spectrum Model for High-dispersion Characterization of Exoplanets and Brown Dwarfs

We present an autodifferentiable spectral modeling of exoplanets and brown dwarfs. This model enables a fully Bayesian inference of the high-dispersion data to fit the ab initio line-by-line spectral computation to the observed spectrum by combining it with the Hamiltonian Monte Carlo in recent probabilistic programming languages. An open-source code, ExoJAX, developed in this study, was written in Python using the GPU/TPU compatible package for automatic differentiation and accelerated linear algebra, JAX. We validated the model by comparing it with existing opacity calculators and a radiative transfer code and found reasonable agreements for the output. As a demonstration, we analyzed the high-dispersion spectrum of a nearby brown dwarf, Luhman 16 A, and found that a model including water, carbon monoxide, and H 2 /He collision-induced absorption was well fitted to the observed spectrum (R = 10 5 and 2.28-2.30 μm). As a result, we found that T0 = 1295$^{-32}_{+35}$ K at 1 bar and C/O = 0.62 ± 0.03, which is slightly higher than the solar value. This work demonstrates the potential of a full Bayesian analysis of brown dwarfs and exoplanets as observed by high-dispersion spectrographs and also directly imaged exoplanets as observed by high-dispersion coronagraphy.

79 ASTRONOMY AND ASTROPHYSICS↗

Subchannel codes: CTF and VIPRE-01

In this chapter, two state-of-the-art and well-used subchannel codes are presented: CTF and VIPRE-01. As part of the Consortium for Advanced Simulation of Light Water Reactors, CTF is being used for single- and two-phase flow analysis under light-water reactor normal and anticipated transient operating conditions. Accurate determination of flow distribution, pressure drop, and void content is crucial for predicting margins to thermal crisis and ensuring more efficient nuclear power plant performance. In preparation for the intended applications, CTF is extensively validated against data from experimental facilities. Meanwhile, the licensed commercial code VIPRE-01 is used to generate a baseline set of simulations for the targeted tests and the solution parameters are compared to CTF results. The closure models in CTF that are in the greatest need of improvement are highlighted and discussed.

Zhao, Xingang↗

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Practical overview of the ICRP Publication 145 phantoms for use with MCNP

Abstract This paper provides a practical overview of the mesh-type reference phantoms (MRCPs) published by the ICRP by providing supplementary guidance on the proper implementation of using the MRCPs in computational codes, such as Monte Carlo N-Particle. Various validation scenarios for photon and neutron exposures were also examined where published values were reproduced using the MRCPs for both external and internal exposures. Generally, the MRCPs performed well in the validation scenarios with the values reproduced being within 10% of the published values.

Hetrick, Lucas (ORCID:0000000178978007)↗

HTS Accelerator Magnets Conceptual Design for Future Lepton Colliders

There is an interest in designing superconducting magnet systems for future lepton circular colliders. This application requires many low-field iron-dominated dipole and quadrupole magnets. Conventional room-temperature magnets are often used because of their low field, low total current, and low power losses. However, high electricity bills for large accelerators drive magnet design to superconductivity. High-temperature superconducting (HTS) magnets can substantially reduce energy losses in magnet systems. The present study investigated the conceptual design of HTS dipole and quadrupole magnets operating in persistent current mode. Energy is transferred into the magnet from an external detachable power source. A continuously circulating current generates a stable magnetic field. The iron-dominated magnet system concept was investigated using OPERA3D code, and the results confirmed the proposed approach’s validity.

43 PARTICLE ACCELERATORS↗

PyFloT v0.1

PyFlot is a combination of python and C++ which aims at tuning floating point precision inside scientific applications. PyFloT uses link time interposition for portability and generality. The set of functions of interest are wrapped with code that enables either profiling or their lowering based on the classification criteria: SLOC (Source code LOCation) or backtrace. PyFloT implementation has only be validated using the GNU compiler toolchain.

Brunie, Hugo↗

Models for synthetic data generation

The software includes a suite of probabilistic statistical/machine learning models that can generate discrete synthetic data. Each model is trained on a set of real (private) data and then it can be used to generate synthetic but statistically similar data. Once ready, the model can generate as many samples as we want. Finally, in addition to the actual models, the software includes code to process data, evaluate results (based on cross validation), and produce reports.

De Oliveira Sales, Ana Paula↗

Transfer Factors for the FRMAC Assessment Manual and Turbo FRMAC to Improve Radiological Dose Assessment

The Turbo FRMAC analysis tool is used to perform complex calculations to quickly evaluate radiological consequences and aid in decision making during an emergency response by assessing impacts to the public, workers, and the food supply. Turbo FRMAC calculations are based on methods established by the Federal Radiological Monitoring and Assessment Center (FRMAC). To be able to assess impacts, input data called transfer factors that describe radionuclide uptake by local plants and animals are required. During the code application exercises in late 2016, identifying, finding, and validating nonstandard transfer factors proved to be time consuming and diverted the teams’ activity away from other critical tasks. As a result, a task was undertaken to dramatically expand the list of available transfer factors (food and non-food) and incorporate these factors into Turbo FRMAC. This will ultimately result in improved efficiency of the assessment team to perform calculations during times when the FRMAC is activated and provide more defensible, vetted data from which to calculate results. As a result of these and subsequent exercises, transfer factors were needed for the following items: bell peppers, Christmas tree, deer, flowers, fresh cucumbers, tomatoes, grapefruit, lichens, mushrooms, oranges, snap beans, squash, strawberries, sugarcane for sugar and seed, sweet corn, tea, tobacco, tree bark, and watermelon. To expand the applicability of the information tables, generic transfer data were also provided for common categories and recommendations were made for expanding the list of chemicals based on chemical similarity. The data presented in this report were compiled from recent literature with most of the data encompassing the period from 2000 to 2018. The following radionuclides were targeted during the literature search: elements associated with reactor accidents or nuclear detonations (Sr 89/90, Cs 134/137, Ce 141/144, Ru 103/106, I 129/131/133, Pu 238/239, Am 241, Zr 95, and Nb 95) and elements associated with industrial accidents or dirty bombs (Ir 192 and Co 60). Data from other elements were evaluated if they were identified during the literature search. In addition, reports were evaluated that were recommended by the research consultants. For each plant or animal transfer factor, the goal was to determine the geometric mean, the geometric standard deviation, the minimum, the maximum, and the number of measurements used. If only one measurement was available that was presented as the mean. Concentration ratios with large geometric standard deviations (GSDs) were generally the result of a paucity of measurements or a few disparate measurements. Measurement disparity was observed for data from different soil types. The geometric mean (GM) is a good reference value for planning and responses purposes, but the location-specific concentrations are unlikely to be similar to the model results. Recommendations for further work include: segregating the data to reflect the influence of soil type, developing approaches to incorporate animal data based on aggregated transfer measurements, and including data for foliar deposition on plants.

61 RADIATION PROTECTION AND DOSIMETRY↗

Experimental Validation of Thermal Hydraulic Behavior in Sodium Fast Reactors (SFR) with the Thermal Hydraulic Experimental Test Article (THETA)

Thermal stratification and transition to natural circulation pose two of the largest sources of uncertainty in systems-level modeling of liquid metal-cooled fast reactors. As these phenomena typically develop during transient event sequences, licensing-basis events analyzed using systemslevel models may have considerable uncertainties associated with thermal-hydraulic parameters of the system to account for these phenomena. As a result, the validation basis for these phenomena for systems-level codes is insufficient to fully support the wide range of liquid metal fast reactors being developed in the US. Currently, the most viable path for licensing a design is to take significant conservatisms and maintain sufficiently large safety margins to account for this uncertainty.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY24: Progress Report on the MELCOR Modeling of the Liquid Salt Test Loop

This report outlines the activities conducted in FY24 focused on updating the liquid salt test loop (LSTL) model through the integration of a new test section featuring 16 radial tubes coupled to a filter section. Some comparative analyses of the updated model with existing experimental data for the LSTL was made and benchmarked against other computational tools, such as the ORNL code, SAM. These actions are part of a comprehensive validation effort and promote collaboration among laboratories participating in the Molten Salt Reactor (MSR) campaign.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation Gap Assessment of a TN B1 Package Containing Fresh, Metallic Sodium-Cooled Fuel

This report describes an assessment of sodium fast reactor assemblies loaded inside a Transnuclear B1 package. The assessment was performed to determine the state of the validation basis for sodium fast reactor fuel within a transportation package originally intended for fresh, light-water reactor fuel assemblies. The study employed the TSUNAMI-3D and TSUNAMI-IP sequences from the SCALE code system to determine correlation coefficients between various configurations and databases of validation assessments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effective Lanthanide Diffusivity through U-Zr Metallic Fuel

Fuel-cladding chemical interaction (FCCI) is one of the main concerns for the performance of U-Zr metallic fuels. Although FCCI results in complex phase transformations and morphological changes in both the fuel and cladding sides of the interface, the formation of a brittle layer in the cladding side of the interface (sometimes referred to as wastage) has one of the most significant impacts on the cladding mechanical integrity. This brittle layer is associated with the formation of intermetallic compounds between the cladding constituents Fe, Cr and the lanthanide fission product species. The formation of this cladding wastage layer has been the focus of recent efforts to develop a mechanistic modeling framework, including production of lanthanides in the fuel, transport through the fuel to the fuel-cladding interface, and phase transformation to form intermetallic phases. To calculate the rate of transport of lanthanides through the fuel, a multi-scale computational approach has been used. Atomistic calculations have been performed to determine the diffusivities of Nd, the most prevalent lanthanide species, through the solid fuel matrix and along pore surfaces. These diffusivities have been used as input parameters for a mesoscale model to determine an effective diffusion coefficient that accounts for the formation of porosity, its interconnection, and infiltration with bond sodium from the fuel-cladding gap. The effective diffusivity has been used in engineering-scale simulations using the BISON fuel performance code, and BISON simulations have been compared to representative validation cases from the EBR-II reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of Process Flowsheets for the Recovery of Metals by Solvent Extraction to Recycle Battery Materials: Evaluation of Application of AMUSE Code. Final CRADA Report

RSR Technologies (RSRT) wished to investigate the use of the AMUSE solvent extraction code, housed at Argonne National Laboratory (ANL) to evaluate the recycling of battery cells for the reclamation of metal values. The modeling studies were undertaken by the group at ANL to evaluate several separation schemes for battery metals using solvent extraction. These studies would develop a framework to be experimentally validated at RSRT at its facility in Dallas, TX. The AMUSE code modeling work would speed the development of recovery flowsheets for spent batteries.

25 ENERGY STORAGE↗