Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral Benchmarks”

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 649 records · Page 36

Domain Decomposition for Integer Optimal Control with Total Variation Regularization

Total variation integer optimal control problems admit solutions and necessary optimality conditions via geometric variational analysis. In spite of the existence of said solutions, algorithms which solve the discretized objective suffer from high numerical cost associated with the combinatorial nature of integer programming. Hence, such methods are often limited to small and medium-sized problems. We propose a globally convergent, coordinate descent–inspired algorithm that allows tractable subproblem solutions restricted to a partition of the domain. Our decomposition method solves relatively small trust-region subproblems that modify the control variable on a subdomain only. Given nontrivial subdomain overlap, we prove that a global first-order necessary optimality condition is equivalent to a first-order necessary optimality condition per subdomain. We additionally show that a sufficient decrease is achieved on a single subdomain by way of a trust-region subproblem solver using geometric measure–theoretic arguments, which we integrate with a greedy patch selection to prove convergence of our algorithm. In conclusion, we demonstrate the practicality of our algorithm on a benchmark large-scale, PDE-constrained integer optimal control problem and find that our method is faster than the state of the art.

domain decomposition↗

Laser Driven Hydrothermal Processing (LDHP) for Rapid Dissolution of Glassy Nuclear Debris

Discovered in 2009 on the SLAB laser system at LLNL, Laser-Driven Hydrothermal Processing (LDHP) both spalls and transiently dissolves a wetted surface (e.g., rock, concrete, paint) when that surface is struck by a laser beam pulse under specific conditions. Previous projects have explored the application of LDHP to expediting sample dissolution for post-detonation nuclear forensics. Early results were very promising, indicating that this method could decrease material dissolution time by a factor of 10 or more. However, the large and complicated lasers previously used are impractical for this application. Under the R^3 venture we aimed to benchmark the LDHP capabilities of turn-key commercial off-the-shelf (COTS) lasers, and to establish how this process could be integrated into an analytical workflow for nuclear forensic analysis. We demonstrated that geological rock standards could be quickly processed into fine powder by LDHP followed by complete acid digestion in 15 minutes, without elemental fractionation that would skew analytical results. We determined that the COTS laser can be used for LDHP that could potentially reduce dissolution time for solid nuclear debris samples in post detonation nuclear forensics. One downside to this approach, however, is that the laser interactions are different with every material, so extensive method validation would be required for all likely material compositions. The LDHP processing rates of several materials were tested and determined to be favorable in the beginning with significant slowdown due to particulate scattering. This could be mitigated either by stopping and centrifuging the sample at specific intervals, or through design of a continuous solvent flow through mechanism.

36 MATERIALS SCIENCE↗

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

TEX-HEU: Integral Experiment Execution with Polyethylene at Very Low Temperatures

The low-temperature variant of the TEX HEU (called Low-Temperature TEX or sometimes LT TEX) campaign is a highly anticipated and necessary experimental series by the greater nuclear science community. Fundamentally, the need for low-temperature integral experiments is required to perform validation of cross sections below room temperature. There has been substantial international interest in low-temperature benchmarks to validate below room temperature cross sections, namely talks given at the 2019 International Conference on Nuclear Criticality (ICNC): UK (Watson, 2019), France (Milin, 2019), and UK (Gan & Wilson, 2019). Additionally, NCSP funded thermal scattering laws (TSLs) were produced by North Carolina State University and require low-temperature benchmarks to validate them. Validation of low-temperature cross sections is also necessary for criticality safety applications. One particularly important application is to ensure that during transportation, fissile materials must remain subcritical under normal ambient conditions which is defined as temperatures down to -40°C/°F by the United States 10 CRF 71 as well as a regulation put forward by the International Atomic Energy Agency (IAEA).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spatio-Temporal Anomaly Detection with Graph Networks for Data Quality Monitoring of the Hadron Calorimeter

The Compact Muon Solenoid (CMS) experiment is a general-purpose detector for high-energy collision at the Large Hadron Collider (LHC) at CERN. It employs an online data quality monitoring (DQM) system to promptly spot and diagnose particle data acquisition problems to avoid data quality loss. In this study, we present a semi-supervised spatio-temporal anomaly detection (AD) monitoring system for the physics particle reading channels of the Hadron Calorimeter (HCAL) of the CMS using three-dimensional digi-occupancy map data of the DQM. We propose the GraphSTAD system, which employs convolutional and graph neural networks to learn local spatial characteristics induced by particles traversing the detector and the global behavior owing to shared backend circuit connections and housing boxes of the channels, respectively. Recurrent neural networks capture the temporal evolution of the extracted spatial features. We validate the accuracy of the proposed AD system in capturing diverse channel fault types using the LHC collision data sets. The GraphSTAD system achieves production-level accuracy and is being integrated into the CMS core production system for real-time monitoring of the HCAL. We provide a quantitative performance comparison with alternative benchmark models to demonstrate the promising leverage of the presented system.

43 PARTICLE ACCELERATORS↗

High-fidelity Multidisciplinary Sensitivity Analysis and Design Optimization for Rotorcraft Applications

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. A sensitivity-enabled fluid dynamics solver and a nonlinear flexible multibody dynamics solver are coupled to predict aerodynamic loads and structural responses of helicopter rotor blades. A discretely consistent adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system is validated by conducting simulations for a benchmark rotorcraft model and comparing solutions with established analyses and experimental data. Sensitivities of lift computed by the multidisciplinary sensitivity analysis are verified by comparison with the sensitivities obtained by complex-variable simulations. Finally the multidisciplinary sensitivity analysis is applied to a constrained gradient-based design optimization for a HART-II rotorcraft configuration.

Wang, Li↗

Accelerating the Inference of the Exa.TrkX Pipeline

Recently, graph neural networks (GNNs) have been successfully used for a variety of particle reconstruction problems in high energy physics, including particle tracking. The Exa.TrkX pipeline based on GNNs demonstrated promising performance in reconstructing particle tracks in dense environments. It includes five discrete steps: data encoding, graph building, edge filtering, GNN, and track labeling. All steps were written in Python and run on both GPUs and CPUs. In this work, we accelerate the Python implementation of the pipeline through customized and commercial GPU-enabled software libraries, and develop a C++ implementation for inferencing the pipeline. The implementation features an improved, CUDA-enabled fixed-radius nearest neighbor search for graph building and a weakly connected component graph algorithm for track labeling. GNNs and other trained deep learning models are converted to ONNX and inferenced via the ONNX Runtime C++ API. The complete C++ implementation of the pipeline allows integration with existing tracking software. We report the memory usage and average event latency tracking performance of our implementation applied to the TrackML benchmark dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗

Nuclear three-body problem in the complex energy plane: Complex-scaling Slater method

The physics of open quantum systems is an interdisciplinary area of research. The nuclear “openness” manifests itself through the presence of the many-body continuum representing various decay, scattering, and reaction channels. As the radioactive nuclear beam experimentation extends the known nuclear landscape toward the particle drip lines, the coupling to the continuum space becomes exceedingly more important. Of particular interest are weakly bound and unbound nuclear states appearing around particle thresholds. Theories of such nuclei must take into account their open quantum nature. To describe open quantum systems, we introduce a complex-scaling (CS) approach in the Slater basis. We benchmark it with the complex-energy Gamow shell model (GSM) by studying energies and wave functions of the bound and unbound states of the two-neutron halo nucleus 6 He viewed as an $α + n + n$ cluster system. Both CS and GSM approaches are applied to a translationally invariant Hamiltonian with the two-body interaction approximated by the finite-range central Minnesota force. In the CS approach, we use the Slater basis, which exhibits the correct asymptotic behavior at large distances. To extract particle densities from the back-rotated CS solutions, we apply the Tikhonov regularization procedure, which minimizes the ultraviolet numerical noise. In this paper, we show that the CS-Slater method is both accurate and efficient. Its equivalence to the GSM approach has been demonstrated numerically for both energies and wave functions of 6 He. One important technical aspect of our calculation was to fully retrieve the correct asymptotic behavior of a resonance state from the complex-scaled (square-integrable) wave function. While standard applications of the inverse complex transformation to the complex-rotated solution provide unstable results, the stabilization method fully reproduces the GSM benchmark. We also propose a method to determine the smoothing parameter of the Tikhonov regularization. The combined suite of CS-Slater and GSM techniques has many attractive features when applied to nuclear problems involving weakly bound and unbound states. While both methods can describe energies, total widths, and wave functions of nuclear states, the CS-Slater method—if it can be applied—can provide additional information about partial energy widths associated with individual thresholds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

OECD RIA Benchmark Phase II - Towards a better understanding of the RIA fuel modelling

This paper deals with the second phase of the RIA (reactivity initiated accident) fuel rod code benchmark for water-cooled reactors that began in 2014. This second phase has been organized around 2 complementary activities: the first activity is to compare the results of different simulations on simplified cases in order to provide additional bases for understanding the differences in modelling of the concerned phenomena, while the second activity is focused on the assessment of the uncertainty of the results. The detailed specification and comparison of the results from the first activity is presented in this paper. Based on the Phase-I and Phase-II conclusions, some generic recommendations can be made. First, fuel and clad thermomechanical models (with the associated material properties) should be further improved and validated more extensively against a sound RIA database. Secondly, build-up of a comprehensive and robust database consisting of both separate-effect tests and integral tests should be pursued in the short term to reinforce the basis for assessment of both individual model and model integration into codes. Thirdly, an assessment of the uncertainty of fuel thermo-mechanics is of high interest (which is consistent with the second activity of this RIA benchmark Phase-II). Finally, as RIA fuel codes are more and more likely to be used for reactor accident studies, particularly for those involving safety analyses, the fuel rod failure criteria (generally used in such analyses) will have to be carefully justified and validated. The current RIA fuel failure criteria are mainly based on the fuel thermal variables and the verification is based on 'conservative' assumptions for the heat transfer conditions. As all codes give rather consistent evaluations of such variables, it appears possible, taking into account adequate provisions, to derive criteria based on thermal variables from experimental values or from an analytical approach. However, if in the future more mechanistic modelling is ever to be used to establish fuel-failure criteria based on mechanical variables, weaknesses in the validation of the codes identified above will have to be accounted for.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NCSP Integral Experiments at Sandia in FY21 [Slides]

This presentation provides details regarding integral experiments at Sandia National Laboratory for fiscal year 2021. The experiments discussed are as follows: IER 230: Characterize the Thermal Capabilities of the 7uPCX; IER 304: Temperature Dependent Critical Benchmarks; IER 305: Critical Experiments with UO 2 Rods and Molybdenum Foils; IER 306: Critical Experiments with UO 2 Rods and Rhodium Foils ; IER 441: Epithermal HEX Lattices with SNL 7uPCX Fuel for Testing Nuclear Data; IER 452: Inversion Point of the Isothermal Reactivity Coefficient; and IER 523: Critical Experiments with ACRR UO 2 -BeO Fuel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluating Cloud Computing in the Proposed NASA DESDynI Ground Data System

The proposed NASA Deformation, Ecosystem Structure and Dynamics of Ice (DESDynI) mission would be a first-of-breed endeavor that would fundamentally change the paradigm by which Earth Science data systems at NASA are built. DESDynI is evaluating a distributed architecture where expert science nodes around the country all engage in some form of mission processing and data archiving. This is compared to the traditional NASA Earth Science missions where the science processing is typically centralized. What's more, DESDynI is poised to profoundly increase the amount of data collection and processing well into the 5 terabyte/day and tens of thousands of job range, both of which comprise a tremendous challenge to DESDynI's proposed distributed data system architecture. In this paper, we report on a set of architectural trade studies and benchmarks meant to inform the DESDynI mission and the broader community of the impacts of these unprecedented requirements. In particular, we evaluate the benefits of cloud computing and its integration with our existing NASA ground data system software called Apache Object Oriented Data Technology (OODT). The preliminary conclusions of our study suggest that the use of the cloud and OODT together synergistically form an effective, efficient and extensible combination that could meet the challenges of NASA science missions requiring DESDynI-like data collection and processing volumes at reduced costs.

Benchmark↗

Fracture Capabilities in Grizzly with the extended Finite Element Method (X-FEM)

Efforts are underway to develop fracture mechanics capabilities in the Grizzly code to enable it to be used to perform deterministic fracture assessments of degraded reactor pressure vessels (RPVs). A capability was previously developed to calculate three-dimensional interaction- integrals to extract mixed-mode stress-intensity factors. This capability requires the use of a finite element mesh that conforms to the crack geometry. The eXtended Finite Element Method (X-FEM) provides a means to represent a crack geometry without explicitly fitting the finite element mesh to it. This is effected by enhancing the element kinematics to represent jump discontinuities at arbitrary locations inside of the element, as well as the incorporation of asymptotic near-tip fields to better capture crack singularities. In this work, use of only the discontinuous enrichment functions was examined to see how accurate stress intensity factors could still be calculated. This report documents the following work to enhance Grizzly’s engineering fracture capabilities by introducing arbitrary jump discontinuities for prescribed crack geometries; X-FEM Mesh Cutting in 3D: to enhance the kinematics of elements that are intersected by arbitrary crack geometries, a mesh cutting algorithm was implemented in Grizzly. The algorithm introduces new virtual nodes and creates partial elements, and then creates a new mesh connectivity; Interaction Integral Modifications: the existing code for evaluating the interaction integral in Grizzly was based on the assumption of a mesh that was fitted to the crack geometry. Modifications were made to allow for the possibility of a crack front that passes arbitrarily through the mesh; and Benchmarking for 3D Fracture: the new capabilities were benchmarked against mixed-mode three-dimensional fracture problems with known analytical solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiscale modeling and characterization for performance and safety of lithium-ion batteries

Lithium-ion batteries are highly complex electrochemical systems whose performance and safety are governed by coupled nonlinear electrochemical-electrical-thermal-mechanical processes over a range of spatiotemporal scales. Gaining an understanding of the role of these processes as well as development of predictive capabilities for design of better performing batteries requires synergy between theory, modeling, and simulation, and fundamental experimental work to support the models. This paper presents the overview of the work performed by the authors aligned with both experimental and computational efforts. In this paper, we describe a new, open source computational environment for battery simulations with an initial focus on lithium-ion systems but designed to support a variety of model types and formulations. This system has been used to create a three-dimensional cell and battery pack models that explicitly simulate all the battery components (current collectors, electrodes, and separator). The models are used to predict battery performance under normal operations and to study thermal and mechanical safety aspects under adverse conditions. This paper also provides an overview of the experimental techniques to obtain crucial validation data to benchmark the simulations at various scales for performance as well as abuse. We detail some initial validation using characterization experiments such as infrared and neutron imaging and micro-Raman mapping. In addition, we identify opportunities for future integration of theory, modeling, and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An Orbital Basis Set for Double Photoionization of Atoms and Molecules

The ab initio theoretical treatment of one-photon double photoionization processes has been limited to atoms and diatomic molecules by the challenges posed by large grid-based representations of the double ionized continuum wave function. To provide a path for extensions to polyatomics, an energy-adapted orbital basis approach is demonstrated that reduces the dimensions of such representations and simultaneously allows larger time steps in time-dependent computational descriptions of double ionization. Additionally, an algorithm that exploits the diagonal nature of the two-electron integrals in the grid basis and dramatically accelerates the transformation between grid and orbital representations is presented. Excellent agreement between the present results and benchmark theoretical calculations is found for H – and Be atoms, as well as the hydrogen molecule, including for the triply differential cross sections that relate the angular distribution and energy sharing of all of the particles in the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward the performance assessment of advanced nuclear waste forms: temperature dependence of lanthanide borosilicate glass dissolution

Lanthanide borosilicate (LaBS) glasses are among the most promising waste forms for the immobilization of high-level radioactive waste generated from advanced nuclear fuel cycles. However, the temperature dependence of their dissolution kinetics remains poorly understood and constrained, limiting the integration of these materials into established performance assessment models. Here, we investigate the dissolution behavior of the legacy AmCm2-19 LaBS glass and the benchmark alkali aluminoborosilicate ISG-1 in deionized water between 50 °C and 250 °C using ASTM C1285 (Product Consistency Test-B) protocols. For AmCm2-19 LaBS glass, normalized elemental release rates for boron and silicon increase with temperature before plateauing near 150 °C, consistent with solubility-limited behavior. From data obtained at 50 °C and 100 °C, Arrhenius analysis yields activation energies of E a (B) = 24.8 ± 0.3 kJ mol⁻¹ and E a (Si) = 14.4 ± 0.2 kJ mol⁻¹, similar or slightly lower than those previously reported for two other compositions of LaBS glasses. No secondary phases or alteration layers were detected by SEM-EDX or pXRD. These results establish one of the first temperature-dependent kinetic datasets for LaBS glass dissolution, providing quantitative parameters to inform mechanistic corrosion models and predictive simulations of glass degradation in geological disposal environments.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measurement of L-shell emission from mid-Z targets under non-LTE conditions using Transmission Grating Spectrometer and DANTE power diagnostics

In this work, we present the measurement of L-band emission from buried Sc/V targets in experiments performed at the OMEGA laser facility. The goal of these experiments was to study non-local thermodynamic equilibrium plasmas and benchmark atomic physics codes. The L-band emission was measured simultaneously by the time resolved DANTE power diagnostic and the recently fielded time integrated Soreq-Transmission Grating Spectrometer (TGS) diagnostic. The TGS measurement was used to support the spectral reconstruction process needed for the unfolding of the DANTE data. The Soreq-TGS diagnostic allows for broadband spectral measurement in the 120 eV–2000 eV spectral band, covering L- and M-shell emission of mid- and high-Z elements, with spectral resolution λ/Δλ = 8–30 and accuracy better than 25%. The Soreq-TGS diagnostic is compatible with ten-inch-manipulator platforms and can be used for a wide variety of high energy density physics, laboratory astrophysics, and inertial confinement fusion experiments.

Fraenkel, M. (ORCID:0000000291939544)↗