Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “APPROXIMATIONS”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermophysical properties of FLiBe using moment tensor potentials

Fluoride salts are prospective materials for applications in some next-generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (66.6% LiF – 33.3% BeF2) predicted by the ab initio (DFT-D3) method. Here, we then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi/Fumi.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Structural coherence model for predicting molten salt thermal conductivity informed by the pair distribution function

To enable thermal behavior prediction and design optimization of molten salt reactors, thermal conductivity of molten salts must be characterized in terms of salt composition and temperature. Current theoretical models fail to provide consistent approximations for all halide mixtures, particularly actinide-bearing melts. This study aims to link the short-range order structure of molten salts to the mean free path of energy carriers through a simple structural coherence model informed by the partial pair distribution function. The proposed method is used to predict the thermal conductivity of 33 alkali and alkaline earth halide salts. Predictions approximate experimental measurements with a mean absolute error of 15.7% for dissociating, complexing, and actinide salts, including unary LiCl, NaCl, and MgCl 2 as well as mixtures LiF–NaF–KF (FLiNaK), LiF–BeF 2 (FLiBe), and NaCl–UCl 3 . The work provides evidence for the validity of energy carrier descriptions of molecular-level heat transfer in molten salts, with implications for improved theories of liquid energy transport in general.

Actinide mixtures↗
DOE OSTI2019

High temperature neutron diffraction investigation of U 3 Si 2

The U-Si system is actively undergoing studies due to its promise as a component of an accident tolerant nuclear fuel. In this work, the crystal structure of the U 3 Si 2 compound in the U-Si system was investigated as a function of temperature from room temperature to 1298 K for a sample of approximately 42.02 at.% Si composition and 1373 K for a sample with approximately 42.10 at.% Si using high temperature time-of-flight neutron diffraction on the High-Pressure Preferred Orientation (HIPPO) diffractometer at Los Alamos Neutron Science Center (LANSCE). The simultaneous Rietveld refinement of five histograms from the five HIPPO detector backs (40°, 60°, 90°, 120°, and 145°) provided datasets for the lattice parameters, anisotropic atomic displacement parameters, and atomic positions as a function of temperature. To explore the possibility of a homogeneity range, two sample compositions were analyzed, a 40.02 at.% Si (near stoichiometric) and a potentially hyperstoichiometric (40.10 at.% Si) sample. While minor differences in the anisotropic atomic displacement parameters between the two samples were observed, over the entire investigated temperature range no additional phases were found indicative of a single phase sample and the ability of the U 3 Si 2 compound to accommodate nonstoichiometry, suggesting that the U 3 Si 2 compound is not a line compound confirming previous determinations. However, the differences observed in the average coefficient of linear thermal expansion identified between the two compositions warrant future investigation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DeepONet-grid-UQ: A trustworthy deep operator framework for predicting the power grid’s post-fault trajectories

This paper proposes a novel data-driven method for the reliable prediction of the power grid’s post-fault trajectories, i.e., the power grid’s dynamic response after a disturbance or fault. Here, the proposed method is based on the recently proposed concept of Deep Operator Networks (DeepONets). Unlike traditional neural networks that learn to approximate functions, DeepONets are designed to approximate nonlinear operators, i.e., mappings between infinite-dimensional spaces. Under this operator framework, we design a novel and efficient DeepONet that (i) takes as inputs the trajectories collected before and during the fault and (ii) outputs the predicted post-fault trajectories. In addition, we endow our method with the much-needed ability to balance efficiency with reliable/trustworthy predictions via uncertainty quantification. To this end, we propose and compare two novel methods that enable quantifying the predictive uncertainty. First, we propose a Bayesian DeepONet (B-DeepONet) that uses stochastic gradient Hamiltonian Monte-Carlo to sample from the posterior distribution of the DeepONet trainable parameters. Then, we design a Probabilistic DeepONet (Prob-DeepONet) that uses a probabilistic training strategy to enable quantifying uncertainty at virtually no extra computational cost. Finally, we validate the proposed methods’ predictive power and uncertainty quantification capability using the New York-New England power grid model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Jensen–Shannon divergence based novel loss functions for Bayesian neural networks

Bayesian neural networks (BNNs) are state-of-the-art machine learning methods that can naturally regularize and systematically quantify uncertainties using their stochastic parameters. Kullback–Leibler (KL) divergence-based variational inference used in BNNs suffer from unstable optimization and challenges in approximating light-tailed posteriors due to the unbounded nature of the KL divergence. To resolve these issues, we formulate a novel loss function for BNNs based on a new modification to the generalized Jensen–Shannon (JS) divergence, which is bounded. In addition, we propose a Geometric JS divergence-based loss, which is computationally efficient since it can be evaluated analytically. We found that the JS divergence-based variational inference is intractable, and hence employed a constrained optimization framework to formulate these losses. Our theoretical analysis and empirical experiments on multiple regression and classification data sets suggest that the proposed losses perform better than the KL divergence-based loss, especially when the data sets are noisy or biased. Specifically, there are approximately 5% and 8% improvements in accuracy for a noise-added CIFAR-10 dataset and a regression dataset, respectively. There is about 13% reduction in false negative predictions of a biased histopathology dataset. Additionally, we quantify and compare the uncertainty metrics for the regression and classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantifying the generalization error in deep learning in terms of data distribution and neural network smoothness

We report the accuracy of deep learning, i.e., deep neural networks, can be characterized by dividing the total error into three main types: approximation error, optimization error, and generalization error. Whereas there are some satisfactory answers to the problems of approximation and optimization, much less is known about the theory of generalization. Most existing theoretical works for generalization fail to explain the performance of neural networks in practice. To derive a meaningful bound, we study the generalization error of neural networks for classification problems in terms of data distribution and neural network smoothness. We introduce the cover complexity (CC) to measure the difficulty of learning a data set and the inverse of the modulus of continuity to quantify neural network smoothness. A quantitative bound for expected accuracy/error is derived by considering both the CC and neural network smoothness. Although most of the analysis is general and not specific to neural networks, we validate our theoretical assumptions and results numerically for neural networks by several data sets of images. The numerical results confirm that the expected error of trained networks scaled with the square root of the number of classes has a linear relationship with respect to the CC. We also observe a clear consistency between test loss and neural network smoothness during the training process. In addition, we demonstrate empirically that the neural network smoothness decreases when the network size increases whereas the smoothness is insensitive to training dataset size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

SympNets: Intrinsic structure-preserving symplectic networks for identifying Hamiltonian systems

Here, we propose new symplectic networks (SympNets) for identifying Hamiltonian systems from data based on a composition of linear, activation and gradient modules. In particular, we define two classes of SympNets: the LA-SympNets composed of linear and activation modules, and the G-SympNets composed of gradient modules. Correspondingly, we prove two new universal approximation theorems that demonstrate that SympNets can approximate arbitrary symplectic maps based on appropriate activation functions. We then perform several experiments including the pendulum, double pendulum and three-body problems to investigate the expressivity and the generalization ability of SympNets. The simulation results show that even very small size SympNets can generalize well, and are able to handle both separable and non-separable Hamiltonian systems with data points resulting from short or long time steps. In all the test cases, SympNets outperform the baseline models, and are much faster in training and prediction. We also develop an extended version of SympNets to learn the dynamics from irregularly sampled data. This extended version of SympNets can be thought of as a universal model representing the solution to an arbitrary Hamiltonian system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

First 3D vector tracking of helium recoils for fast neutron measurements at Super$\mathrm{KEKB}$

Here we present results from the first deployment of novel, high definition, compact gas Time Projection Chambers (TPCs) with pixel chip readout as part of the BEAST II beam background measurement project at SuperKEKB. The TPCs provide detailed 3D imaging of ionization from neutron-induced nuclear recoils in a helium and carbon dioxide target gas mixture at standard temperature and pressure. We present the TPC performance and the neutron backgrounds observed during the initial stage of collider commissioning. We find excellent electron background rejection, leading to background-free nuclear recoil measurements above 50 keV ee , despite the extreme high-background environment. We measure an angular resolution better than 20° for recoil tracks longer than 1.7 mm, corresponding to an average ionization energy of approximately 100 keV ee . We also obtain the full 3D vector direction of helium recoils by utilizing charge profile measurements along the recoil axis, with a correct head/tail assignment efficiency of approximately 80%. With this performance, we present comparisons between measured and simulated event rates, recoil energy spectra, and directional distributions originating from beam–gas and Touschek beam losses at SuperKEKB. We utilize head/tail recognition to distinguish neutron components traveling with positive radial velocity in the Belle II coordinate system from those traveling in the opposite direction. Finally, we present a novel method of discriminating beam–gas interactions from Touschek beam losses that can eliminate the need for dedicated accelerator runs for background measurements. This method is still statistics-limited. However, future studies should be able to verify this method, which in turn could lead to neutron background analysis runs symbiotic with normal Belle II operation. The capabilities demonstrated here also suggest that high definition recoil imaging in gas TPCs is applicable to low energy, low-background experiments, such as directional dark matter searches.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

A flexible event reconstruction based on machine learning and likelihood principles

Event reconstruction is a central step in many particle physics experiments, turning detector observables into parameter estimates; for example estimating the energy of an interaction given the sensor readout of a detector. A corresponding likelihood function is often intractable, and approximations need to be constructed. Here, in our work, we first show how the full likelihood for a many-sensor detector can be broken apart into smaller terms, and secondly how we can train neural networks to approximate all terms solely based on forward simulation. Our technique results in a fast, flexible, and close-to-optimal surrogate model proportional to the likelihood and can be used in conjunction with standard inference techniques allowing for a consistent treatment of uncertainties. We illustrate our technique for parameter inference in neutrino telescopes based on maximum likelihood and Bayesian posterior sampling. Given its great flexibility, we also showcase our method for geometry optimization enabling to learn optimal detector designs. Lastly, we apply our method to realistic simulation of a ton-scale water-based liquid scintillator detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Absolute light yield of the EJ-204 plastic scintillator

The absolute light yield of a scintillator, defined as the number of scintillation photons produced per unit energy deposited, is a useful quantity for scintillator development, research, and applications. Yet, literature data on the absolute light yield of organic scintillators are limited. The goal of this work is to assess the suitability of the EJ-204 plastic scintillator from Eljen Technology to serve as a reference standard for measurements of the absolute light yield of organic scintillators. Four EJ-204 samples were examined: two manufactured approximately four months prior and stored in high-purity nitrogen, and two aged approximately eleven years and stored in ambient air. The scintillator response was measured using a large-area avalanche photodiode calibrated using low energy $γ$-ray and X-ray sources. The product of the quantum efficiency of the photodetector and light collection efficiency of the housing was characterized using an experimentally-benchmarked optical photon simulation. The average absolute light yield of the fresh samples, 9100 ± 400 photons per MeV, is lower than the manufacturer-reported value of 10400 photons per MeV. Moreover, the aged samples demonstrated significantly lower light yields, deviating from the manufacturer specification by as much as 26%. These results are consistent with recent work showcasing environmental aging in plastic scintillators and suggest that experimenters should use caution when deploying plastic scintillators in photon counting applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Measured and projected beam backgrounds in the Belle II experiment at the SuperKEKB collider

The Belle II experiment at the SuperKEKB electron–positron collider aims to collect an unprecedented data set of 50 ab -1 to study CP -violation in the B -meson system and to search for Physics beyond the Standard Model. SuperKEKB is already the world’s highest-luminosity collider. In order to collect the planned data set within approximately one decade, the target is to reach a peak luminosity of 6 x 10 35 cm −2 s −1 by further increasing the beam currents and reducing the beam size at the interaction point by squeezing the betatron function down to β $^*_y$ = 0.3 mm. To ensure detector longevity and maintain good reconstruction performance, beam backgrounds must remain well controlled. We report on current background rates in Belle II and compare these against simulation. We find that a number of recent refinements have significantly improved the background simulation accuracy. Finally, we estimate the safety margins going forward. We predict that backgrounds should remain high but acceptable until a luminosity of at least 2.8 x 10 35 cm −2 s −1 is reached for β $^*_y$ = 0.6 mm. At this point, the most vulnerable Belle II detectors, the Time-of-Propagation (TOP) particle identification system and the Central Drift Chamber (CDC), have predicted background hit rates from single-beam and luminosity backgrounds that add up to approximately half of the maximum acceptable rates.

Detector background↗
DOE OSTI2025

Charge collection efficiency of diamond and silicon sensors irradiated with alpha particles

To evaluate the viability of using semiconductors as sensor materials in a detector for the Associated Particle Imaging technique, the radiation hardness of silicon and diamond diodes to alpha particles has been assessed. Here, the detector lifetimes for both silicon and diamond sensors were measured under the prolonged exposure to alpha particles emitted by an 241 Am source. The silicon detector was exposed to alpha radiation for approximately two months, reaching an accumulated fluence of ~ 1.5 x 10 12 α cm –2 . Additionally, by using a high purity single-crystal diamond with coplanar electrodes operating with full charge collection, the diamond detector response was measured over approximately ten months reaching an accumulated fluence of over 6 x 10 12 α cm –2 cm.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Validation of SOLPS-ITER simulations with kinetic, fluid, and hybrid neutral models for JET-ILW low-confinement mode plasmas

For JET L-mode plasmas in low-recycling conditions (electron temperature at the outer strike point, T e,ot ≳ 30 eV ), SOLPS-ITER simulations agree within the error bars for the experimental profiles at the low-field side (LFS) divertor target. The peak Balmer-α (Dα ) emission in the LFS divertor agrees within the error bars of the KS3 filterscope diagnostic, but is approximately 30% lower than the peak value of the KT1 spectrometer. Simulations have been performed with fluid, kinetic, and hybrid models for the neutrals. The large fluid-kinetic discrepancies of more than a factor 2 are successfully corrected by using a hybrid fluid-kinetic approach, for which kinetic atoms are transferred to the fluid population when the local Knudsen number of the atom becomes smaller than a user-defined transition Knudsen number Kn t . The hybrid-kinetic discrepancies are limited to a few % for Kn t ≤ 100 . When increasing the upstream density to high-recycling conditions, at the onset of detachment ( T e,ot ≈ 5 eV ), the simulations predict more than a factor 2 lower peak ion saturation current to the LFS divertor than the experiments. Also the Dα emission is underpredicted with approximately a factor 2. For these high-recycling conditions, the fluid-kinetic discrepancies are limited to maximum 50%, which are again corrected by using the hybrid approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Catalyzed oxidation of IG-110 nuclear graphite by simulated fission products Ag and Pd nanoparticles

To evaluate the stability of nuclear materials in high temperature gas reactors under air ingress conditions, catalytic oxidation of IG-110 graphite by two simulated fission products, metallic Pd and Ag, was studied in oxidative atmosphere and at temperatures up to 1000°C using an integrated furnace, mass spectroscopy and infrared spectroscopy system. Transmission electron microscopy and X-ray diffraction studies show that Pd and Ag nanoparticles were successfully introduced onto powdery IG-110 graphite through an impregnation and subsequent heat-treatment process. The combined mass spectroscopy and infrared spectroscopy methods allow simultaneous analysis of two gaseous products, CO and CO 2 , and separate measurements of activation energy for their formation reactions. It was found that the introduction of Pd or Ag to IG-110 graphite substantially catalyzed the oxidation of graphite, characteristic of decreased onset temperatures for the oxidation of graphite. Moreover, the catalytic effects by Pd and Ag are considerably different based on measured concentration ratios of CO 2 to CO as a function of oxidation temperatures. Ag makes the graphite oxidation commence at approximately 400°C with CO 2 being the dominant product. In contrast, Pd significantly increases the concentration ratio of CO 2 to CO at temperatures higher than approximately 690°C, although it decreases the onset temperature for the oxidation reaction to around 525°C. To understand the catalytic difference, the mechanism of the graphite oxidation was discussed based on the changes of surface oxygen species on Ag and Pd.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Molecular dynamics study of interstitial He clusters in nickel

This study presents a molecular dynamics analysis focusing on the behavior of interstitial helium (He) clusters in nickel (Ni), examining their formation, stability, and migration energetics. Consistent with previous research, we found that the binding energies of interstitial helium within a helium cluster are positive and increase with the cluster size, indicating a preference for helium atoms to cluster together. However, our findings also reveal that while the formation energy increases monotonically with cluster size, the increase in binding energy is non-monotonic. Importantly, small He clusters were observed to be thermally unstable at reactor operational temperatures (approximately 600 K), with the He 2 cluster exhibiting instability even at room temperature. With a binding energy of 0.49 eV for a He 4 cluster, we hypothesize that for helium bubbles to form via homogeneous nucleation (i.e., through trap mutation) at reactor operating temperatures, the helium concentration must be high enough to facilitate the formation of helium clusters of at least size 4 or larger. As expected, interstitial helium and small helium clusters are highly mobile. This mobility was observed not only at room temperature but also at temperatures as low as approximately 200 K. Furthermore, the mean squared displacement method has been utilized to determine the migration barriers and the corresponding prefactors for clusters ranging from He 1 to He 6

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Development and assessment of a reactor system prognosis model with physics-guided machine learning

Autonomous control systems provide recommendations to help operators in decision-making during plant operations ranging from normal operation to accident management. An important step of autonomous control is prognosis. In nuclear engineering domain, prognosis is the process of predicting future conditions of a system or equipment based on present signs and symptoms of a fault. The prognosis model allows predicting future reactor states for possible candidate control strategies so that the outcomes can be evaluated to determine the best control strategy. The prognosis model requires representing direct relationships between the symptoms and the predictions. In nuclear engineering, computational simulations are approximate representations of the operation of the real system. However, prognosis with computational simulations requires high computation power and time due to possible large number of scenarios. Necessary computation resources can be reduced with machine learning (ML) approach for fast predictions by building a surrogate function using the simulation data. A critical issue is, ML models are ignorant of physical knowledge, and these models approximate statistical relationships between the system variables. This ignorance can produce results that are inconsistent with physical laws, even if an optimal result is achieved from a mathematical point of view. Physics-guided machine learning (PGML) is an approach to tackle this issue. Here, this work formulates and illustrates a framework to guide development and assessment of the ML-based prognosis model for autonomous control systems. The development of the prognosis model considers the training of a ML model which consists of optimizing many aspects of the ML approach. The assessment of the prognosis model considers training data limitations and uncertainties of the ML approach. Prognosis models with standalone ML and PGML are developed and assessed on the loss-of-flow scenario of Experimental Breeder Reactor II. The results indicate that PGML based prognosis model has the best performance compared to other prognosis models.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗
DOE OSTI2024

Wrought FeCrAl alloy (C26M) cladding behavior and burst under simulated loss-of-coolant accident conditions

Cladding burst experiments for FeCrAl cladding were performed in the Severe Accident Test Station facility at Oak Ridge National Laboratory. These experiments were simulated using the BISON fuel performance code to better understand the cladding plastic behavior and failure under simulated loss-of-coolant accident conditions. 3D cladding surface boundary conditions were generated using composite axial and azimuthal profiles from experiment thermocouple data. To improve the simulation analysis capabilities in BISON for cladding burst behavior, new thermal creep, plasticity, and failure stress models specific to C26M, a wrought FeCrAl alloy, were developed and implemented. Initial cladding burst results indicated a general underprediction in the failure temperature of the six cladding burst simulations versus the observed failure temperatures. Close investigation of the experiment timing versus the underlying tensile test data revealed that, compared with the tensile specimens, the cladding tubes did not experience the same long holding time at high temperatures. New tensile tests were performed at high temperatures using a temperature ramp similar to the simulated loss-of-coolant accident experiments. These new tensile curves showed an approximately 80% increase in the ultimate tensile strength of the C26M alloy, indicating that a holding time of 10 min at 700 °C and 800 °C allows annealing to change the material microstructure. Using the updated tensile properties, the burst temperatures and stresses from the simulations showed remarkable agreement with the experimental results. This study was then extended by varying the initial pressure to highlight the burst temperature difference between standard Zircaloy-4 and C26M cladding under equivalent conditions. The results show that C26M has a burst temperature that is approximately 70–130 K greater than that of Zircaloy-4. In conclusion, these modeling predictions can be further improved by collecting high-temperature tensile data for C26M beyond the temperature ranges used in this work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗