Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Benchmark testing”

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 613 records · Page 34

Verification and Performance Impact of the New Parallel MCNP6.3 Particle Track Output Capability for Subcritical Multiplication Simulations [Slides]

A separate MCNP6.3 V&V document reports on all the default calculations for all test suites. This report does not include the subcritical multiplication benchmark suite. After some additional clean-up and finalizing the post-processing and documentation steps, the subcritical multiplication benchmark suite will be released in the next version of our vnvstats repository. We tested the new HDF5 PTRAC feature in MCNP6.3 and found encouraging outcomes. Identical results coming out of the simulation with respect to the legacy PTRAC results. The overall runtime for all simulations is reduced by ~20% with the new HDF5 PTRAC capability. We consider giving the new HDF5 PTRAC features a try and using it for all subcritical multiplication and any other relevant (PTRAC) calculations.

97 MATHEMATICS AND COMPUTING↗

RELAP5-3D in the OECD-NEA HTGR Thermal Hydraulics Benchmark

The OECD-NEA HTGR Thermal Hydraulics benchmark uses data from the High Temperature Test Facility for a series of code-to-code and code-to-data exercises aimed at improving the state of knowledge on existing thermal hydraulics modeling and simulation tools for prismatic HTGR applications. This benchmark includes problems representing hot gas mixing in the lower plenum, the depressurized conduction cooldown accident and the pressurized conduction cooldown accident. This presentation discusses the benchmark and RELAP5-3D's role in that benchmark, including as a tool for predicting behavior in the core and as a tool for providing boundary conditions to computational fluid dynamics analysis in the lower plenum.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Evaluation of 3D Simulation Capabilities within xRAGE using a Suite of Test Problems

This presentation describes a 3D verification and validation suite for xRAGE. Three verification test problems for hydrodynamics are presented: the Kidder ball problem, the Verney shell problem, and a 5-material compression problem. There is excellent agreement between 2D and 3D XRAGE simulation results for these problems, and between the xRAGE results and the benchmark solutions. In addition, two 3D ICF test problems are presented, based on an OMEGA direct drive capsule experiment, and on a NIF indirect drive capsule experiment. These two ICF test problems are used to evaluate two hydrodynamic methods within xRAGE. In addition to hydrodynamics the simulations also include three temperature plasma physics, thermal conduction and radiation diffusion. When simulating either capsule in 3D the newer unsplit hydrodynamic method in xRAGE produces more vorticity relative to the older default method. For the indirect drive capsule the 3D simulations are in reasonable agreement with the experimental values of ion temperature and neutron production. Unstable 3D vortex rings are present near the DT ice/carbon interface which enhance the confinement of the DT ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reducing the Cost of CCSD Basis Set Extrapolation in Ab Initio Computational Thermochemistry

Here, a series of approximations to CCSD contributions in computational model chemistries is presented in the context of kcal mol –1 , kJ mol –1 , and 20 cm –1 theoretical predictions of total atomization energies, benchmarked within the HEAT+CH 4 test suite. A specific set of circumstances where MP2, without empirical scaling, may be used as an effective intermediate in the first two of these accuracy ranges was determined. However, SDQ-MP4, a method long used in pursuit of kcal mol –1 accuracy but relatively unstudied in the subchemical accuracy community, offers significant improvement over the quality of MP2 as a basis-set intermediate at significantly reduced cost compared to CCSD. Given this, we argue for SDQ-MP4 as the de facto CCSD basis-set intermediate in sub-chemical accuracy calculations when CCSD in a desired basis set becomes unaffordable. We additionally report on a “CBS-like” scheme, where MP2 and SDQ-MP4 are used in conjunction to create a “cheap” three-part approximation of large CCSD basis set limits. The data for the CCSD approximation schemes are organized in such a way that model chemistry developers can locate an analog of their current approach for the CCSD basis set limit and explore alternative intermediates that either decrease computational cost or increase computational accuracy. We also show, for a handful of molecules, that SDQ-MP4 shows promise as an effective basis-set intermediate for harmonic and fundamental frequency computations, allowing for zero-point corrections of nearly CCSD(T)/ANO1 quality using simple composite methods that only require CCSD(T)/ANO0.

Thorpe, James H. [Argonne National Laboratory (ANL↗

Hotspot Coevolution Is a Key Identifier of Near-Native Protein Complexes

Protein–protein interactions play a key role in mediating numerous biological functions, with more than half the proteins in living organisms existing as either homo- or hetero-oligomeric assemblies. Protein subunits that form oligomers minimize the free energy of the complex, but exhaustive computational search-based docking methods have not comprehensively addressed the challenge of distinguishing a natively bound complex from non-native forms. Current protein docking approaches address this problem by sampling multiple binding modes in proteins and scoring each mode, with the lowest-energy (or highest scoring) binding mode being regarded as a near-native complex. However, high-scoring modes often match poorly with the true bound form, suggesting a need for improvement of the scoring function. Here in this study, we propose a scoring function, KFC-E, that accounts for both conservation and coevolution of putative binding hotspot residues at protein–protein interfaces. We tested KFC-E on four benchmark sets of unbound examples and two benchmark sets of bound examples, with the results demonstrating a clear improvement over scores that examine conservation and coevolution across the entire interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-assembly of aramid amphiphiles into ultra-stable nanoribbons and aligned nanoribbon threads

Small molecule self-assembly is an established route for producing high surface area nanostructures with readily customizable chemistries and precise molecular organization. However, these structures are fragile, exhibiting molecular exchange, migration, and rearrangement - among other dynamic instabilities - and ultimately disassociate upon drying. Here we show a small molecule platform, the aramid amphiphile (AA), that overcomes these dynamic instabilities by incorporating a Kevlar-inspired domain into the molecular structure. Strong, anisotropic interactions between AAs suppress molecular exchange and elicit spontaneous self-assembly in water to form nanoribbons with lengths of up to 20 microns. The nanoribbons have a Young’s modulus of 1.7 GPa and tensile strength of 1.9 GPa. Here, we exploit this stability to extend small molecule self-assembly to hierarchically ordered macroscopic materials outside of solvated environments. Through an aqueous shear alignment process, we organize AA nanoribbons into arbitrarily long flexible threads that support 200 times their weight when dried. Tensile tests of the threads provide a benchmark for Young’s moduli (between ~400 and 600 MPa) and extensibilities (between ~0.6 and 1.1%) that depend on the counterion chemistry in solution. This bottom-up approach to macroscopic materials could benefit solid-state applications, historically inaccessible by self-assembled nanomaterials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlation consistent effective core potentials for late 3 d transition metals adapted for plane wave calculations

We construct a new modification of correlation consistent effective core potentials (ccECPs) for late 3d elements Cr–Zn with Ne-core that are adapted for efficiency and low energy cut-offs in plane wave calculations. The decrease in accuracy is rather minor, so that the constructions are in the same overall accuracy class as the original ccECPs. The resulting new constructions work with energy cut-offs at or below ≈400 Ry and, thus, make calculations of large systems with transition metals feasible for plane wave codes. We also provide the basic benchmarks for atomic spectra and molecular tests of this modified option that we denote as ccECP-soft.

Chemistry↗

Statistical inference of anomalous thermal transport with uncertainty quantification for interpretive 2D SOL models

The critical task of inferring anomalous cross-field transport coefficients is addressed in simulations of boundary plasmas with fluid models. A workflow for parameter inference in the UEDGE fluid code is developed using Bayesian optimization with parallelized sampling and integrated uncertainty quantification. In this workflow, transport coefficients are inferred by maximizing their posterior probability distribution, which is generally multidimensional and non-Gaussian. Uncertainty quantification is integrated throughout the optimization within the Bayesian framework that combines diagnostic uncertainties and model limitations. As a concrete example, we infer the anomalous electron thermal diffusivity $\chi_\perp$ from an interpretive 2D model describing electron heat transport in the conduction-limited region with radiative power loss. The workflow is first benchmarked against synthetic data and then tested on H-, L-, and I-mode discharges to match their midplane temperature and divertor heat flux profiles. We demonstrate that the workflow efficiently infers diffusivity and its associated uncertainty, generating 2D profiles that match 1D measurements. Future efforts will focus on incorporating more complicated fluid models and analyzing transport coefficients inferred from a large database of experimental results.

Bayesian optimization↗

Carbon-Carbon Aeroshell Chemical Reaction Model in an Open RPS in Martian Environment

A new RPS is being designed to be open to the vacuum of space to improve its thermal efficiency during deep space missions. This design, however, could lead to some potential complications if it is ever used in the low-pressure gas environment of mission destinations like Mars. This report documents the development of a chemical kinetics reaction model of the FWPF carbon-carbon aeroshell in the Martian gaseous environment. This model uses an “infinite sink” assumption where the RPS is open to the Martian environment with open exchange of flowing gases, so the gas chemistry over the aeroshell does not change with the reaction and represents the worst case scenario of the chemical reaction of the FWPF in the Martian environment. The model combines the kinetics equation of both oxidation reactions that would occur on Mars: the Langmuir-Hinshelwood kinetics of the Boudouard reaction of CO 2 with the chemical kinetics of O 2 oxidation of the carbon-carbon. The FWPF chemical reaction model was benchmarked using the Opila’s experimental chemical reaction test data (1) in a simulated Martian environment to increase the validity of the model’s predicted mass loss rate. The output of the benchmarked chemical reaction model showed a trade-off between the maximum temperature of the aeroshell to have less than 10% reaction at the mission life. For example, an RPS with a five-year intended mission life on Mars has a maximum aeroshell temperature of 835 K / 562 °C, increasing the mission life to ten years reduces the maximum operating temperature to 810 K / 527 °C. If the aeroshell operating temperature in the new RPS on the surface of Mars is 960 Kelvin, the model calculates that the life is limited to 0.22 years / 81 days. Increasing the aeroshell temperature to 1343 Kelvin, further reduces the life to 0.375 days / 9 hrs. This illustrates the thermally activated nature of the chemical kinetics model in limiting the life of the aeroshell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BISON TRISO Modeling Advancements and Validation to AGR-1 Data

BISON is a finite element-based nuclear fuel performance code. Among its unique characteristics are its ability to model 1D, 2D, and 3D geometries and its applicability to a wide variety of nuclear fuels. For the last eight years, BISON has included a beginning capability to model tri-structural isotropic (TRISO) fuel. Recently, interest in TRISO fuel has grown, and a significant effort has been made to improve BISON’s capabilities in this area. Capability development has occurred for each material present in TRISO fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. New elastic, creep, swelling, thermal expansion, thermal conductivity, and fission gas release (FGR) models are available. New models for the graphite matrix are also now available. Another important addition is the ability to perform statistical failure analysis of large samples of fuel particles. This new capability, which continues to grow, enables evaluation of failure due to pressure or crack formation by analyzing many thousands of particles. This enables realistic calculations of fission product release from the many particles in a TRISO-fueled reactor. These capabilities were checked via regression and verification tests. A large number of code benchmarking problems were also run, showing that BISON’s results closely match those of other software tools. Finally, a significant validation effort was completed in which fission product release, measured as part of the AGR-1 capsule experiments, was compared to BISON outputs. BISON outputs compared very well to the experimental data and to PARFUME results. Interest in BISON’s TRISO capabilities is growing, with the U.S. Nuclear Regulatory Commission (NRC) and Westinghouse Electric Company receiving training during the past year. Multiple other entities have expressed interest in or are actively using BISON. Kairos Power, LLC, has a strong partnership with Idaho National Laboratory (INL) regarding the use of BISON for TRISO analysis. While its capabilities still continue to grow, BISON has already become a powerful tool for TRISO analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

INDEN Fissile Actinide Issues [Slides]

This presentation addresses the outlook for the International Nuclear Data Evaluation Network (INDEN), specifically regarding general issues that have been solved, including Thermal Neutron Constants (TNC) and Prompt Fission Neutron Spectra (PFNS). Other topics discussed include outstanding issues in fissile evaluations, comments on the RPI quasi-differential experiment, and reaction rate testing of trial evaluations in fast assemblies. The presentation concludes by discussing how INDEN interactions strongly helped to find deficiencies in existing evaluations and highlighted potential solutions to existing challenges, that updated trial evaluations for U-235, U-233, Pu-239 are available for testing, that further work is expected on U-235, Pu-239, and U-233, and that performance was tested on the ICSBEP and/or SINBAD benchmarks, during which significant improvement was demonstrated. The INDEN scheme of international collaboration on nuclear data evaluation is working well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Zirconium Scattering Sensitivity in Neutron Transport Calculations of Multiplying Systems [Slides]

This presentation touches on motivation and applications of Zirconium. Performance of evaluated libraries on double and quasi differential experiments is considered. Additionally, the performance of evaluated libraries in select ICSBEP benchmarks is considered. Provided is a detailed look in the performance of evaluated libraries in preliminary benchmark configurations for IER-516: Zirconium Test Assembly (ZTA). Conclusions and recommended experiments to satisfy Zirconium scattering nuclear data needs are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation of ENDF/B-VIII.1 Nuclear Data Files [Slides]

ENDF-6 formatted files were processed into A Compact ENDF (ACE) files using NJOY2016. Several validation tests were performed: (1) LANL Legacy Benchmark Suite (2) “Modern” Benchmark Suite (3) HEU Benchmark Suite (4) LEU Benchmark Suite (5) Mixed (U+Pu) Benchmark Suite (6) Pu Benchmark Suite (7) 233 U Benchmark Suite.

97 MATHEMATICS AND COMPUTING↗

Finite-Geometry Corrections for Neutron Scattering Cross Section Measurements

In this report, we describe the corrections necessary for neutron scattering cross section measurements in experiments which use finite-geometry scattering samples. A new Geant4 based framework was developed to calculate these corrections which has been benchmarked against a well-measured test case. This framework also calculates the outgoing neutron energy dependence of the correction factors, which is crucial when resolving individual final states in the scatterer is not possible.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗