Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear analysis”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

EPR of Photoexcited Triplet-State Acceptor Porphyrins

The photoexcited triplet states of porphyrin architectures are of significant interest in a wide range of fields including molecular wires, nonlinear optics, and molecular spintronics. Electron paramagnetic resonance (EPR) is a key spectroscopic tool in the characterization of these transient paramagnetic states singularly well suited to quantify spin delocalization. Previous work proposed a means of extracting the absolute signs of the zero-field splitting (ZFS) parameters, D and E, and triplet sublevel populations by transient continuous wave, hyperfine measurements, and magnetophotoselection. Here, we present challenges of this methodology for a series of mesoperfluoroalkyl-substituted zinc porphyrin monomers with orthorhombic symmetries, where interpretation of experimental data must proceed with caution and the validity of the assumptions used in the analysis must be scrutinized. The EPR data are discussed alongside quantum chemical calculations, employing both DFT and CASSCF methodologies. Despite some success of the latter in quantifying the magnitude of the ZFS interaction, the results clearly provide motivation to develop improved methods for ZFS calculations of highly delocalized organic triplet states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Porting the Nonlinear Optimization Library HiOp to Accelerator-Based Hardware Architectures

While interior point method has been the centerpiece of nonlinear programming tools used in science and engineering, its reliance on linear solvers that can tackle sparse symmetric indefinite and highly ill-conditioned problems made it difficult to implement it effectively on hardware accelerators. HiOp optimization package attempts to provide an implementation of the interior point method suitable for hardware accelerators by compressing the original sparse problem to produce an underlying linear problem that is dense and of manageable size. Implementations of dense linear solvers are more mature and utilize hardware accelerators better than their sparse counterparts. There is a number of important domain problems, such as optimal power flow analysis for power grids, where the sparse problem can be effectively compressed and deploying dense linear solver within the interior point method can improve performance. Here we describe a portable implementation of HiOp optimization engine, which uses a linear solver from Magma library and runs entirely on hardware accelerators. To compress the problem, HiOp uses customized mixed dense-sparse linear algebra. All HiOp kernels are implemented using Umpire and RAJA portability libraries. We describe details of the implementation and discuss trade-offs between performance, portability and development cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Understanding latent timescales in neural ordinary differential equation models of advection-dominated dynamical systems

The neural ordinary differential equation (ODE) framework has shown considerable promise in recent years in developing highly accelerated surrogate models for complex physical systems characterized by partial differential equations (PDEs). For PDE-based systems, state-of-the-art neural ODE strategies leverage a two-step procedure to achieve this acceleration: a nonlinear dimensionality reduction step provided by an autoencoder, and a time integration step provided by a neural-network based model for the resultant latent space dynamics (the neural ODE). This work explores the applicability of such autoencoder-based neural ODE strategies for PDEs in which advection terms play a critical role. More specifically, alongside predictive demonstrations, physical insight into the sources of model acceleration (i.e., how the neural ODE achieves its acceleration) is the scope of the current study. Such investigations are performed by quantifying the effects of both autoencoder and neural ODE components on latent system time-scales using eigenvalue analysis of dynamical system Jacobians. To this end, the sensitivity of various critical training parameters – de-coupled versus end-to-end training, latent space dimensionality, and the role of training trajectory length, for example – to both model accuracy and the discovered latent system timescales is quantified. Furthermore, this work specifically uncovers the key role played by the training trajectory length (the number of rollout steps in the loss function during training) on the latent system timescales: larger trajectory lengths correlate with an increase in limiting neural ODE time-scales, and optimal neural ODEs are found to recover the largest time-scales of the full-order (ground-truth) system. Demonstrations are performed across fundamentally different unsteady fluid dynamics configurations influenced by advection: (1) the Kuramoto–Sivashinsky equations (2) Hydrogen-Air channel detonations (the compressible reacting Navier–Stokes equations with detailed chemistry), and (3) 2D Atmospheric flow.

Advection-dominated dynamical systems↗
DOE OSTI2021

Structural color from a coupled nanowire pair beyond the bonding and antibonding model

Optical resonances in nanostructures can be harnessed to produce a wide range of structural colors. Conversely, the analysis of structural colors has been used to clarify the nature of optical resonances. Here, we show that silicon nanowire (NW) pairs can display a wide range of structural colors by controlling their radiative coupling. This is accomplished by exciting a series of Fabry–Pérot-like modes where light is repeatedly scattered between two NWs. These modes are beyond the expectation from the conventional chemical bonding model under a quasi-electrostatic approximation, in which only bonding and antibonding modes can be formed in a pair system through modal hybridization. The additional eigenmodes found in a two-resonator system originate from the nonlinear, frequency-dependent coupling strength derived from the radiative nature of low- Q resonators. The Fabry–Pérot modes can be tuned across the entire visible frequency range by varying the distance between two NWs, leading to what we believe is a new type of universal building blocks that can provide structural color within a subwavelength footprint. The presented results pave the way toward the design and usage of highly tunable resonances that exploit the radiative coupling of high-index nanostructures.

Li, Qitong↗
DOE OSTI2021

GDSA PFLOTRAN Development (FY2021)

The Spent Fuel & Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2021 advances of the PFLOTRAN Development group of the SFWST Campaign. The mission of this group is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of generic disposal concepts. In FY 2021, development proceeded along three main thrusts: software infrastructure, code performance, and process model advancement. Software infrastructure improvements included implementing an Agile software development framework and making improvements to the QA Test Suite. Code performance improvements included development of advanced linear and nonlinear solvers as well as design of flexible smoothing algorithms for capillary pressure functions. Process modeling advancements included the addition of flexible thermal conductivity function definitions and refinement of multi-continuum reactive transport to support Sandia’s participation in DECOVALEX

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

PFLOTRAN Development FY2021

The Spent Fuel & Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is to develop a disposal system modeling and analysis capability for valuating disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2021 advances of the PFLOTRAN Development group of the SFWST Campaign. The mission of this group is to develop a geologic disposal system modeling capability for nuclear waste that can be used to probabilistically assess the performance of generic disposal concepts. In FY 2021, development proceeded along three main thrusts: software infrastructure, code performance, and process model advancement. Software infrastructure improvements included implementing an Agile software development framework and making improvements to the QA Test Suite. Code performance improvements included development of advanced linear and nonlinear solvers as well as design of flexible smoothing algorithms for capillary pressure functions. Process modeling advancements included the addition of flexible thermal conductivity function definitions and refinement of multi-continuum reactive transport to support Sandia’s participation in DECOVALEX. This report fulfills the GDSA PFLOTRAN Development Work Package Level 3 Milestone – PFLOTRAN Development, FY2021, M3SF-21SN010304072.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Mechanical Loads on Spent Nuclear Fuel in the General 30 cm Package Drop Scenario

The U.S. Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) research program is guided by the high-level goal of closing prioritized knowledge gaps related to spent nuclear fuel (SNF) storage and transportation, which are summarized by Saltzstein et al. (2020). One of the high-priority knowledge gaps is the identification and quantification of mechanical loads that are expected to affect SNF during normal conditions of transportation and storage to inform the range of physical SNF test programs. This report uses modeling and analysis methods to estimate the mechanical loads on spent nuclear fuel (SNF) in the general 30 cm package drop scenario. The drop scenario assumes impact limiters are in place in the transportation configuration and the impact surface is perfectly rigid. The goal of this analysis is to consider the universe of potential mechanical loading conditions that can happen to SNF and present the results in a manner that is useful for materials testing, decision making, and regulatory rule making purposes. This study uses validated finite element models and methods to perform a broad parametric study of key variables that can affect the mechanical loads on SNF during a hypothetical package free drop scenario. Physical drop test data from a cask and fuel assembly drop test campaign is the basis for model validation. Additionally, the results of the parametric study are used to inform a damage model, which uses multiple nonlinear regression to estimate the relationships between input variables and output response. The parametric finite element analyses consider thousands of input variable combinations, while the damage model estimates millions of combinations. The breadth of this study provides confidence that the potential range of mechanical loads that SNF might experience during the general 30 cm package drop scenario are characterized well enough to consider this knowledge gap closed. While this report documents the overall peak values calculated in this study, the 95 th percentile values, the histograms, and the observed trends are equally important. This study covered a large range of SNF temperatures, room temperature to 300°C, and burnups, 10 GWd/MTU to 62 GWd/MTU. Each temperature and burnup combination has a different cladding yield strain, so it is more meaningful to summarize the calculated cladding strain response as its factor of safety, which is defined relative to the yield strain. The factor of safety is calculated as the yield strain divided by the peak cladding strain. A factor of safety greater than unity indicates that the cladding remains below yield, whereas a value less than unity is indicative of plastic deformation. In all cases of this study a safety factor over 1.0 was calculated, although in the most limiting case at 300°C the safety factor was only 1.01, which suggests that yielding could occur when additional loads like rod internal pressure are included. When the temperature is restricted to 200°C the limiting safety factor increases to 1.28, which has significant margin to accommodate internal pressure and potential local cladding defects that could cause a local stress concentration. An important trend in the calculated fuel rod mechanical loads is that the 2 nd highest loaded fuel rod in an assembly tends to be significantly lower than the highest loaded rod. The implication is that even if one rod in an assembly experiences a failure the loads would have to be significantly higher to cause two or more rods to fail.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Effects of Accident Thermal Loading on In-Plane Shear Behavior of Steel-Plate Composite Walls

Structural walls in safety-related nuclear facilities are required to be designed for seismic and accident thermal (due to postulated high-energy pipe break events) loading combination. Current U.S. and international codes provide limited guidance for analysis and design of walls for this loading combination. This paper describes the experimental results and observations from tests conducted on a laboratory-scale (1:4 to 1:5) test unit representing steel-plate composite (SC) walls subjected to combined in-plane (seismic) and accident thermal loading. The test unit was subjected to surface temperatures of up to 450°F in combination with cyclic in-plane loading. Results of similar experiments recently conducted in Japan are also summarized (with surface temperatures up to 570°F). Surface heating combined with the low thermal conductivity and high specific heat of concrete resulted in nonlinear thermal gradients through the thickness of the specimens. These nonlinear thermal gradients and the associated self- or internal restraint led to extensive concrete cracking. This concrete cracking reduced the initial and secant stiffness of the specimens. The initial stiffness of the heated specimens reduced to 30 to 40% of the initial stiffness of the control (unheated) specimen. The secant stiffness of the heated specimens reduced up to 50% of the secant stiffness of the control (unheated) specimen. However, the in-plane shear strength of the heated SC specimens was still approximately 10 to 30% greater than the nominal in-plane shear strength calculated, for the limit state of steel plate von Mises yielding, using AISC N690 equations and measured material properties. Evaluation of the experimental results and observations suggests that the in-plane shear strength of SC walls subjected to typical accident surface temperatures (up to 570°F) can be estimated conservatively using the current provisions of AISC N690. The stiffness for accident thermal loading combinations can be considered to reduce from cracked composite stiffness at ambient temperature to fully cracked—that is, steel-only stiffness as the surface temperature increases up to typical accident value.

Construction & Building Technology↗
DOE OSTI2021

Theory of the charge density wave in A V 3 Sb 5 kagome metals

The family of metallic kagome compounds AV 3 Sb 5 (A = K, Rb, Cs) was recently discovered to exhibit both superconductivity and charge order. The nature of the charge density wave (CDW) phase is presently unsettled, which complicates the interpretation of the superconducting ground state. In this paper, we use group theory and density functional theory (DFT) to derive and solve a phenomenological Landau model for this CDW state. The DFT results reveal three unstable phonon modes with the same in-plane momentum but different out-of-plane momenta, whose frequencies depend strongly on the electronic temperature. This is indicative of an electronically driven CDW, stabilized by features of the in-plane electronic dispersion. Motivated by the DFT analysis, we construct a Landau free-energy expansion for coupled CDW order parameters with wave vectors at the M and L points of the hexagonal Brillouin zone. We find an unusual trilinear term coupling these different order parameters, which can promote the simultaneous condensation of both CDWs even if the two modes are not nearly degenerate. We classify the different types of coupled multi-Q CDW orders, focusing on those that break the sixfold rotational symmetry and lead to a unit-cell doubling along all three crystallographic directions, as suggested by experiments. Here, we determine a region in parameter space, characterized by large nonlinear Landau coefficients, where these phases—dubbed staggered trihexagonal and staggered Star of David—are the leading instabilities of the system. Finally, we discuss the implications of our results for the kagome metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Higher-order symmetry plane correlations in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The correlations between event-by-event fluctuations of symmetry planes are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ALICE detector at the Large Hadron Collider. This analysis is conducted using the Gaussian estimator technique, which is insensitive to biases from correlations between different flow amplitudes. The study presents, for the first time, the centrality dependence of correlations involving up to five different symmetry planes. The correlation strength varies depending on the harmonic order of the symmetry plane and the collision centrality. Comparisons with measurements from lower energies indicate no significant differences within uncertainties. Additionally, the results are compared with hydrodynamic model calculations. Although the model predictions provide a qualitative explanation of the experimental results, they overestimate the data for some observables. This is particularly true for correlators that are sensitive to the nonlinear response of the medium to initial-state anisotropies in the collision system. As these new correlators provide unique information—independent of flow amplitudes—their usage in future model developments can further constrain the properties of the strongly interacting matter created in ultrarelativistic heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Expanded Capabilities and Opportunities for Virtual Beamline Code

In the 25 years since LLNL broke ground on the National Ignition Facility (NIF) - a cornerstone of the National Nuclear Security Administration’s Stockpile Stewardship Program - the NIF and Photon Science (NIF&PS) Directorate has steadily pushed the boundaries of laser physics, nonlinear optics, and photonics in service of inertial confinement fusion experiments, and advanced photon sources development. Beyond the 192 beamlines of the primary laser facility with its world-record energy and fusion yields, NIF&PS capabilities have expanded to include chirped-pulse amplification, kilojoule petawatt-class short-pulse systems to generate hard x-ray radiographic probes, high-average power lasers as secondary sources for generation or wake field particle acceleration along directed energy applications, and other scientific pursuits. The Computing Directorate partners with NIF&PS on numerous projects and technologies from diagnostic measurements and high performance control systems to data analysis, information technology infrastructure, and scientific simulation code development. “This strong collaboration enables cutting-edge science, particularly with regard to physics-based modeling and simulation,” says physicist Jean-Michel Di Nicola. One crucial capability for NIF&PS is the Virtual Beamline (VBL) laser simulation code. It can model all the major laser physics and technology involved in the design optimization, commissioning, and operations of advanced laser architectures from tabletop to NIF-scale. VBL provides researchers with high-fidelity models and high-resolution calculations of laser performance predictions—including for the entire NIF laser system, the Advanced Radiographic Capability, parts of the High Repetition-Rate Advanced Petawatt Laser System laser delivered to the Czech Republic, and the Optical Science Laser. After more than two decades of experimentally verified physics and computing enhancements, this unique workhorse code is wrapping up another major milestone: migration from Java to C++ with a flurry of user interface and optimization features, as well as laser physics and high-resolution enhancements thanks to parallel execution on Livermore Computing platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Use of Sobol’ Variance-Based Global Sensitivity Analysis and Multidimensional Legendre Polynomial Fitting for Reduced Order Modeling

Sandia National Laboratories (SNL) has developed a novel reduced order modeling approach. Prioritization of inputs is accomplished using Sobo' indices obtained through a more efficient variance-based global sensitivity analysis. To determine the Sobo' functions, simulated input values are aligned to collocation points to permit the use of Gauss-Lobatto integration, thereby reducing the number of simulation trials needed by more than an order of magnitude compared to standard Monte Carlo approaches. Furthermore, by leveraging the orthogonality of Legendre polynomials in conjunction with those same simulations at the collocation nodes, an efficient fitting method is developed to represent the Sobo' functions from which a reduced order model (ROM) is constructed. The developed method is both more efficient computationally, and the resulting ROM is more accurate. The efficacy of this technique is demonstrated on a nonlinear polynomial test function as well as the nonlinear Ishigami and Sobo' g functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Wide-blocky veins explained by dependency of crystal growth rate on fracture surface type: Insights from phase-field modeling

Vein microstructures contain a wealth of information on coupled chemical and mechanical processes of fracturing, fluid transport, and crystal growth. Numerical simulations have been used for exploring the factors controlling the development of vein microstructures; however, they have not been quantitatively validated against natural veins. Here we combined phase-field modeling with microtextural analysis of previously unexplained wide-blocky calcite veins in natural limestone and of the fresh fracture surface in this limestone. Results show that the wide-blocky vein textures can only be reproduced if ~10%–20% of crystals grow faster than the rest. This fraction corresponds to the amount of transgranularly broken grains that were observed on the experimental fracture surfaces, which are dominantly intergranular. We hypothesize that transgranular fractures allow faster growth of vein minerals due to the lack of clay coatings and other nucleation discontinuities that are common along intergranular cracks. Our simulation results show remarkable similarity to the natural veins and reproduce the nonlinear relationship between vein crystal width and vein aperture. This allows accurate simulations of crystal growth processes and related permeability evolution in fractured rocks.

Geology↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

A comprehensive and fair comparison of two neural operators (with practical extensions) based on $\mathrm{FAIR}$ data

Neural operators can learn nonlinear mappings between function spaces and offer a new simulation paradigm for real-time prediction of complex dynamics for realistic diverse applications as well as for system identification in science and engineering. Herein, we investigate the performance of two neural operators, which have shown promising results so far, and we develop new practical extensions that will make them more accurate and robust and importantly more suitable for industrial-complexity applications. The first neural operator, DeepONet, was published in 2019 (Lu et al., 2019), and its original architecture was based on the universal approximation theorem of Chen & Chen (1995). The second one, named Fourier Neural Operator or FNO, was published in 2020, and it is based on parameterizing the integral kernel in the Fourier space. DeepONet is represented by a summation of products of neural networks (NNs), corresponding to the branch NN for the input function and the trunk NN for the output function; both NNs are general architectures, e.g., the branch NN can be replaced with a CNN or a ResNet. According to Kovachki et al. (2021), FNO in its continuous form can be viewed conceptually as a DeepONet with a specific architecture of the branch NN and a trunk NN represented by a trigonometric basis. In order to compare FNO with DeepONet computationally for realistic setups, we develop several extensions of FNO that can deal with complex geometric domains as well as mappings where the input and output function spaces are of different dimensions. We also develop an extended DeepONet with special features that provide inductive bias and accelerate training, and we present a faster implementation of DeepONet with cost comparable to the computational cost of FNO, which is based on the Fast Fourier Transform. Here we consider 16 different benchmarks to demonstrate the relative performance of the two neural operators, including instability wave analysis in hypersonic boundary layers, prediction of the vorticity field of a flapping airfoil, porous media simulations in complex-geometry domains, etc. We follow the guiding principles of FAIR (Findability, Accessibility, Interoperability, and Reusability) for scientific data management and stewardship. The performance of DeepONet and FNO is comparable for relatively simple settings, but for complex geometries the performance of FNO deteriorates greatly. We also compare theoretically the two neural operators and obtain similar error estimates for DeepONet and FNO under the same regularity assumptions.

42 ENGINEERING↗
DOE OSTI2025

An entropy-based debiasing approach to quantifying experimental coverage for novel applications of interest in the nuclear community

This manuscript proposes a novel information-theoretic approach to the quantification of experimental relevance, i.e., coverage, to achieve optimal data assimilation results for nuclear engineering applications. Specifically, this work posits the need for a new metric, called coverage (q C ) of an application’s quantity of interest, i.e., eigenvalue or power peaking for an advanced reactor concept, defined herein as the theoretically maximum achievable reduction in the quantity’s uncertainty given measurements from a pool of experiments in a manner that is independent of the data assimilation procedure employed. Currently, reduction in a quantity’s uncertainty is strongly biased by the underlying assumptions of the assimilation procedure to account for the under-determined nature of such problems and the similarity criterion employed to identify relevant experiments. To address this challenge, this work has developed a coverage metric, q C , based on mutual information, which establishes a new conceptual framework for assessing coverage, one that is independent of the model parameters and responses degree of variations in both the experimental and application domains, i.e., linear vs non-linear, and their prior uncertainty distributions, i.e., Gaussian vs. non-Gaussian. The q C is an entropic measure capable of addressing coverage for general nonlinear problems with non-Gaussian uncertainties and inclusive of the measurement uncertainties from multiple experiments. Numerical experiments from manufactured analytical problems as well as a set of benchmarks from the ICSBEP handbook are employed to demonstrate its theoretical and practical performance as compared to the c k -based experiment selection methodology, commonly employed in the neutronic community. The manuscript then employs other well-known adaptations to existing data assimilation methodologies for nonlinear and non-Gaussian problems capable of achieving the coverage posited by q C .

Bayesian data assimilation↗
DOE OSTI2020

Probabilistic Groundwater Modeling of the RDX Plume at Los Alamos National Laboratory to Support Risk Assessment - 20359

A plume of the contaminant hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with concentrations greater than the New Mexico tap water drinking standard (9.66 ppb) is present in the regional aquifer near the southwestern boundary of Los Alamos National Laboratory (LANL). A risk assessment is performed for exposure to regional aquifer groundwater, with long-term predictions of RDX concentrations provided by a calibrated, probabilistic, numerical fate and transport model. The structure of the model is hierarchical, with the RDX Regional Aquifer (RA) groundwater model acting as the primary tool for analysis of downgradient RDX concentrations. The RA model is deeply informed by the conceptual site model (CSM) and is calibrated using site RDX concentration data and hydraulic head measurements, along with other analyses. The model is calibrated using data through December 2019. Where data are scarce other lines of evidence are used to inform inputs, including the multiphase RDX Vadose Zone (VZ) model and Pipe and Disk analytical screening tool. Model inputs are described with informative prior distributions using a robust approach to development that incorporates all available lines of evidence for every parameter used as an input to the RA model. Calibration is performed using a classical nonlinear optimization routine, which is then used to initialize a Bayesian calibration. The Bayesian calibration constrains the uncertainty in the classical calibration, ultimately providing posterior distributions for all model parameters. The challenges of the calibration include high-dimensional parameter space, including spatially heterogeneous hydraulic conductivities, comparatively sparse data, and low RDX concentrations. Posterior parameter distributions developed in the calibration process are then used for stochastic predictive model runs into the future. The result of the forward runs is spatially and temporally explicit estimates of head and concentration with uncertainty at all points in space and time. The probabilistic modeling approach presented here includes innovative computational and statistical methods that leverage high-performance computing (HPC) resources. It makes use of prior modeling work performed at the chromium plume site in the central LANL area, with extensive updates. The risk assessment will ultimately be used to support decision-making at the site, using multiple appropriately weighted sources of information, as well as uncertainty, sensitivity, and value-of-information analyses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗