Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model interpretability”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Reduced-Order Modeling of Hydrogen Releases from Vent Stacks and with Wind Effects

Here, the physical release behavior of hydrogen is important to understand from a safety and design perspective. The consequences of unignited pressurized gaseous hydrogen plumes exiting vent stacks were considered by extending and modifying existing general hydrogen plume models. Entrainment, vent stack backpressure, and the flow regime of hydrogen exiting the vent were found to be significant factors affecting plume shape and size, but further investigation and validation with unchoked, low-Froude-number flows is recommended to improve the model’s robustness. Additionally, models for the effects of wind on unignited plume momentum and entrainment were added to explore this behavior. Wind was assumed to increase mixing of hydrogen with the ambient air, and to affect the momentum of the released jet. Introducing wind into the plume model led to a shorter plume for all wind and jet directions. A high counter-flowing wind led to non-physical results and challenges in interpreting the visualization. The proposed jet plume wind sub-models (specifically entrainment coefficients) were fit and compared to experimental data of different releases of hydrogen into a wind tunnel, but the quantity of data available and experimental conditions were limited. Thus, collection of more empirical data and for a wider range of conditions is recommended for improvement of the proposed computational models. Developing reduced-order models for these physical phenomena can improve accessibility to predicted physical behavior and the rate at which hydrogen systems can be safely designed and deployed.

entrainment↗
DOE OSTI2023

A New Modeling Approach for a Priori Uncertainties of Laser Tracker Angle Measurements

Methods for modeling the uncertainty in laser tracker angle measurements vary within the metrology industry, leading to confusion and questionable stochastic modeling for survey network adjustments and error propagation analysis. Interpreting the published laser tracker manufacturer performance specifications to determine an a priori sigma value for weighting azimuth and zenith angle measurements can be confusing and has led to differing implementations. Furthermore, this paper proposes a unique way to model survey network a priori laser tracker angular uncertainties based on laser tracker manufacturers’ published maximum permissible error (MPE) values referenced to current standards for weighting survey network measurements. This paper’s proposed model takes into account the disparate effects that pointing errors, target centering errors, and leveling errors have on azimuth and zenith angular uncertainties for measurements with steep sightings and at near ranges.

42 ENGINEERING↗
DOE OSTI2021

YAHFC: A Code Framework to Model Nuclear Reactions and Estimate Correlated Uncertainties

Reaction modeling is a key ingredient in designing experiments and interpreting their results, and is an essential component in the process of evaluating nuclear data and assembling nuclear data libraries used in nuclear technology applications. Typically, experimental data are available only for a handful of reaction channels and theory models are used to fill in the gaps. In addition, theory is often called upon as the arbitrator between discrepant data. Most importantly, theory and modeling are required for an accurate determination of uncertainties in the evaluated data and the correlations between the multiple channels. A fast, accurate, and flexible modeling capability has been developed at LLNL with the code system YAHFC (Yet Another Hauser-Feshbach Code). YAHFC is a Monte Carlo, Hauser-Feshbach code framework, making full use of dynamic memory allocation, derived types, and parallel computing. YAHFC can generate events to simulate experiments and is guiding experiments designed to measure inelastic neutron scattering from actinide targets. YAHFC is also being used to analyze decays from surrogate experiments, thereby enabling the inference of reaction cross sections inaccessible by direct measurement. Finally, by modeling nuclear reactions with constraints from experimental data, YAHFC can deliver complete nuclear data libraries, with evaluated uncertainties, using the modernized Generalized Nuclear Data Structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Generation of random geological models using multi-randomization for machine learning

Generating high-fidelity geological models is essential for advancing machine learning (ML) methods in automated seismic interpretation. For instance, seismic images paired with corresponding fault labels are foundational for ML-based fault detection from seismic migration sections. While several open-access datasets of random geological models exist, open-source tools specifically designed to produce large volumes of such models for ML applications remain scarce. To address this gap, we present RGM (Random Geological Model), an open-source software package for efficiently generating 2D and 3D synthetic geological models tailored for ML workflows. RGM supports the creation of diverse model components, including medium property distributions (P-/S-wave velocities and density), seismic reflectivity images (i.e., synthetic migration sections), relative geological time, and discrete fault attributes such as probability, dip, strike, rake, and displacement. It also accommodates the creation of complex geological features such as salt bodies and unconformities. The model generation algorithm employs a multi-randomization strategy, yielding an effectively infinite-dimensional model space that encompasses a wide range of geological scenarios and associated seismic features. Furthermore, RGM incorporates a method to generate synthetic elastic migration images using analytical elastic reflection coefficients combined with frequency-dependent scaling. This functionality enables the creation of training datasets for ML models that leverage elastic seismic images. RGM is implemented in modern object-oriented Fortran, allowing users to flexibly control statistical parameters governing model variability. We demonstrate the capability, performance, and geological realism of the package through comprehensive 2D and 3D examples.

58 GEOSCIENCES↗
DOE OSTI2023

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗
DOE OSTI2021

Hamiltonian reduction of Vlasov–Maxwell to a dark slow manifold

In this work, we show that non-relativistic scaling of the collisionless Vlasov–Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov–Maxwell phase space. Vlasov–Maxwell dynamics restricted to the slow manifold recovers the Vlasov–Poisson and Vlasov–Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov–Darwin model more generally. The slow manifold may be interpreted to all orders in perturbation theory as a collection of formal Vlasov–Maxwell solutions that do not excite light waves, and are therefore ‘dark’. We provide a heuristic lower bound for the time interval over which Vlasov–Maxwell solutions initialized optimally near the slow manifold remain dark. We also show how the dynamics on the slow manifold naturally inherits a Hamiltonian structure from the underlying system. After expressing this structure in a simple form, we use it to identify a manifestly Hamiltonian correction to the Vlasov–Darwin model. The derivation of higher-order terms is reduced to computing the corrections of the system Hamiltonian restricted to the slow manifold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Prediction of Complex Observed Shear Wave Splitting Patterns at Ryukyu Subduction Zone Using a Strong Intra-Slab Anisotropy Model

Complex shear wave splitting (SWS) patterns in subduction zones are often interpreted geodynamically as resulting from complex mantle flow; however, this may not always be necessary. We analyzed 7,093 high-quality SWS measurements from teleseismic S waves recorded by Hi-net stations across the Ryukyu arc in Japan. Our findings show a systematic rotation of the fast S polarization from trench-parallel to trench-perpendicular depending on the earthquake backazimuth. For the same earthquake, the measured splitting patterns also vary spatially across the southwest Japan. Using full-wave seismic modeling, we showed that a dipping slab with ~30% shear anisotropy of the tilted transverse isotropy (TTI) type, with a symmetry axis perpendicular to the slab interface, can predict the observed delay times and polarization rotation. Our results highlight the importance of considering dipping anisotropic slabs in interpreting SWS at subduction zones.

58 GEOSCIENCES↗
DOE OSTI2023

Ultrafast radiographic imaging and tracking: An overview of instruments, methods, data, and applications

Ultrafast radiographic imaging and tracking (U-RadIT) use state-of-the-art ionizing particle and light sources to experimentally study sub-nanosecond transients or dynamic processes in physics, chemistry, biology, geology, materials science and other fields. These processes are fundamental to modern technologies and applications, such as nuclear fusion energy, advanced manufacturing, communication, and green transportation, which often involve one mole or more atoms and elementary particles, and thus are challenging to compute by using the first principles of quantum physics or other forward models. One of the central problems in U-RadIT is to optimize information yield through, e.g. high-luminosity X-ray and particle sources, efficient imaging and tracking detectors, novel methods to collect data, and large-bandwidth online and offline data processing, regulated by the underlying physics, statistics, and computing power. We review and highlight recent progress in: (a.) Detectors such as high-speed complementary metal-oxide semiconductor (CMOS) cameras, hybrid pixelated array detectors integrated with Timepix4 and other application-specific integrated circuits (ASICs), and digital photon detectors; (b.) U-RadIT modalities such as dynamic phase contrast imaging, dynamic diffractive imaging, and four-dimensional (4D) particle tracking; (c.) U-RadIT data and algorithms such as neural networks and machine learning, and (d.) Applications in ultrafast dynamic material science using XFELs, synchrotrons and laser-driven sources. Hardware-centric approaches to U-RadIT optimization are constrained by detector material properties, low signal-to-noise ratio, high cost and long development cycles of critical hardware components such as ASICs. Interpretation of experimental data, including comparisons with forward models, is frequently hindered by sparse measurements, model and measurement uncertainties, and noise. Alternatively, U-RadIT make increasing use of data science and machine learning algorithms, including experimental implementations of compressed sensing. Machine learning and artificial intelligence approaches, refined by physics and materials information, may also contribute significantly to data interpretation, uncertainty quantification and U-RadIT optimization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Accurate contrast determination for X-ray speckle visibility spectroscopy

X-ray speckle visibility spectroscopy using X-ray free-electron lasers has long been proposed as a probe of fast dynamics in noncrystalline materials. In this paper, numerical modeling is presented to show how the data interpretation of visibility spectroscopy can be impacted by the nonidealities of real-life X-ray detectors. Using simulated detector data, this work provides a detailed analysis of the systematic errors of several contrast extraction algorithms in the context of low-count-rate X-ray speckle visibility spectroscopy and their origins are discussed. Here, it was found that the finite detector charge cloud and pixel size lead to an unavoidable `degeneracy' in photon position determination, and that the contrasts extracted using different algorithms can all be corrected by a simple linear model. The results suggest that experimental calibration of the correction coefficient at the count rate of interest is possible and essential. This allows computationally lightweight algorithms to be implemented for on-the-fly analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

PQ axiverse

We show that the strong CP problem is solved in a large class of compactifications of string theory. The Peccei-Quinn mechanism solves the strong CP problem if the CP-breaking effects of the ultraviolet completion of gravity and of QCD are small compared to the CP-preserving axion potential generated by low-energy QCD instantons. We characterize both classes of effects. To understand quantum gravitational effects, we consider an ensemble of flux compactifications of type IIB string theory on orientifolds of Calabi-Yau hypersurfaces in the geometric regime, taking a simple model of QCD on D7-branes. We show that the D-brane instanton contribution to the neutron electric dipole moment falls exponentially in N 4 , with N the number of axions. In particular, this contribution is negligible in all models in our ensemble with N > 17. We interpret this result as a consequence of large N effects in the geometry that create hierarchies in instanton actions and also suppress the ultraviolet cutoff. We also compute the CP breaking due to high-energy instantons in QCD. In the absence of vectorlike pairs, we find contributions to the neutron electric dipole moment that are not excluded, but that could be accessible to future experiments if the scale of supersymmetry breaking is sufficiently low. The existence of vectorlike pairs can lead to a larger dipole moment. Finally, we show that a significant fraction of models are allowed by standard cosmological and astrophysical constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Search for a CP-odd Higgs boson decaying into a heavy CP-even Higgs boson and a Z boson in the ${\ell}^{+}{\ell}^{-}t\overline{t}$ and $\nu\overline{\nu}b\overline{b}$ final states using 140 fb –1 of data collected with the ATLAS detector

A search for a heavy CP-odd Higgs boson, A, decaying into a Z boson and a heavy CP-even Higgs boson, H, is presented. It uses the full LHC Run 2 dataset of pp collisions at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector, corresponding to an integrated luminosity of 140 fb –1 . The search for A → ZH is performed in the $ℓ^+ℓ^–t\bar{t}$ and $ν\bar{ν}b\bar{b}$ final states and surpasses the reach of previous searches in different final states in the region with mH > 350 GeV and mA > 800 GeV. No significant deviation from the Standard Model expectation is found. Upper limits are placed on the production cross-section times the decay branching ratios. Limits with less model dependence are also presented as functions of the reconstructed $m(t\bar{t}$) and $m(b\bar{b}$) distributions in the $ℓ^+ℓ^–t\bar{t}$ and $ν\bar{ν}b\bar{b}$ channels, respectively. In addition, the results are interpreted in the context of two-Higgs-doublet models

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Observation of electroweak production of two jets in association with an isolated photon and missing transverse momentum, and search for a Higgs boson decaying into invisible particles at 13 TeV with the ATLAS detector

This paper presents a measurement of the electroweak production of two jets in association with a Zγ pair, with the Z boson decaying into two neutrinos. It also presents a search for invisible or partially invisible decays of a Higgs boson with a mass of 125 GeV produced through vector-boson fusion with a photon in the final state. These results use data from LHC proton–proton collisions at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector and corresponding to an integrated luminosity of 139 fb -1 . The event signature, shared by all benchmark processes considered for the measurements and searches, is characterized by a significant amount of unbalanced transverse momentum and a photon in the final state, in addition to a pair of forward jets. Electroweak $Zγ$ production in association with two jets is observed in this final state with a significance of 5.2 (5.1 expected) standard deviations. The measured fiducial cross-section for this process is 1.31 ± 0.29 fb. An observed (expected) upper limit of 0.37 (0.34$^{+0.15}_{-0.10}$) at 95% confidence level is set on the branching ratio of a 125 GeV Higgs boson to invisible particles, assuming the Standard Model production cross-section. The signature is also interpreted in the context of decays of a Higgs boson into a photon and a dark photon. An observed (expected) 95% CL upper limit on the branching ratio for this decay is set at 0.018 (0.017$^{+0.007}_{-0.005}$), assuming the Standard Model production cross-section for a 125 GeV Higgs boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

KGML-ag: a modeling framework of knowledge-guided machine learning to simulate agroecosystems: a case study of estimating N<sub>2</sub>O emission using data from mesocosm experiments

Abstract. Agricultural nitrous oxide (N2O) emission accounts for a non-trivial fraction of global greenhouse gas (GHG) budget. To date, estimating N2O fluxes from cropland remains a challenging task because the related microbial processes (e.g., nitrification and denitrification) are controlled by complex interactions among climate, soil, plant and human activities. Existing approaches such as process-based (PB) models have well-known limitations due to insufficient representations of the processes or uncertainties of model parameters, and due to leverage recent advances in machine learning (ML) a new method is needed to unlock the “black box” to overcome its limitations such as low interpretability, out-of-sample failure and massive data demand. In this study, we developed a first-of-its-kind knowledge-guided machine learning model for agroecosystems (KGML-ag) by incorporating biogeophysical and chemical domain knowledge from an advanced PB model, ecosys, and tested it by comparing simulating daily N2O fluxes with real observed data from mesocosm experiments. The gated recurrent unit (GRU) was used as the basis to build the model structure. To optimize the model performance, we have investigated a range of ideas, including (1) using initial values of intermediate variables (IMVs) instead of time series as model input to reduce data demand; (2) building hierarchical structures to explicitly estimate IMVs for further N2O prediction; (3) using multi-task learning to balance the simultaneous training on multiple variables; and (4) pre-training with millions of synthetic data generated from ecosys and fine-tuning with mesocosm observations. Six other pure ML models were developed using the same mesocosm data to serve as the benchmark for the KGML-ag model. Results show that KGML-ag did an excellent job in reproducing the mesocosm N2O fluxes (overall r2=0.81, and RMSE=3.6 mgNm-2d-1 from cross validation). Importantly, KGML-ag always outperforms the PB model and ML models in predicting N2O fluxes, especially for complex temporal dynamics and emission peaks. Besides, KGML-ag goes beyond the pure ML models by providing more interpretable predictions as well as pinpointing desired new knowledge and data to further empower the current KGML-ag. We believe the KGML-ag development in this study will stimulate a new body of research on interpretable ML for biogeochemistry and other related geoscience processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Self-Supervised and Interpretable Anomaly Detection Using Network Transformers

Machine learning and deep neural networks (DNNs) have been proposed as a tool to identify anomalies in computer network communications. However, due the obfuscated nature of off-the-shelf machine learning models, their output often does not provide enough information to isolate the source of the anomaly to take corrective measures. In this article, we introduce the network transformer (NeT), a DNN model for anomaly detection that incorporates the graph structure of the communication network in order to improve interpretability. Further, the presented approach has the following advantages: first, enhanced interpretability by incorporating the graph structure of computer networks; second, provides a hierarchical set of features that enables analysis at different levels of granularity; second, self-supervised training that does not require labeled data. The NeT model was evaluated on a set of anomalous scenarios executed in a real industrial control system. The presented approach successfully identified the anomalies, the devices affected, and the specific connections causing the anomalies, providing a data-driven hierarchical approach to analyze the behavior of a cyber network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Permeation of CO 2 and N 2 through glassy poly(dimethyl phenylene) oxide under steady- and presteady-state conditions

Glassy polymers are often used for gas separations because of their high selectivity. Although the dual-mode permeation model correctly fits their sorption and permeation isotherms, its physical interpretation is disputed, and it does not describe permeation far from steady state, a condition expected when separations involve intermittent renewable energy sources. To develop a more comprehensive permeation model, we combine experiment, molecular dynamics, and multiscale reaction–diffusion modeling to characterize the time-dependent permeation of N 2 and CO 2 through a glassy poly(dimethyl phenylene oxide) membrane, a model system. Simulations of experimental time-dependent permeation data for both gases in the presteady-state and steady-state regimes show that both single- and dual-mode reaction–diffusion models reproduce the experimental observations, and that sorbed gas concentrations lag the external pressure rise. The results point to environment-sensitive diffusion coefficients as a vital characteristic of transport in glassy polymers.

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

Measurement of high-mass $t\bar{t}\ell ^{+}\ell ^{-}$ production and lepton flavour universality-inspired effective field theory interpretations at $\sqrt{s}=13$ $\text {T}\text {e}\hspace{-1.00006pt}\text {V}$ with the ATLAS detector

Measurements of $t\bar{t}\ell ^{+}\ell ^{-}$ production in the region of high dilepton invariant mass with effective field theory (EFT) interpretations are presented. They are performed using final states with three isolated leptons (electrons or muons) and are based on $\sqrt{s} = 13$ TeV proton–proton collision data with an integrated luminosity of $140\,\textrm{fb}^{-1}$, recorded from 2015 to 2018 with the ATLAS detector at the Large Hadron Collider. Measurements of the $t\bar{t}\ell ^{+}\ell ^{-}$ signal strength and cross-section upper-limits are performed inclusively in lepton flavour and separately for electrons and muons. The study also aims to probe anomalous four-fermion interactions including to test for possible lepton flavor universality violation. No significant deviations from the Standard Model predictions are observed and the measurements are interpreted through the EFT formalism to provide new constraints on the relevant operators.

Aad, G. [CNRS/IN2P3] (ORCID:0000000266654934)↗
DOE OSTI2024

Elucidating Photoinduced Processes of Photosystem I Via Multidimensional Electronic and Vibrational Spectroscopies

This project was motivated by an overarching goal to elucidate the mechanism of energy and electron transfer that governs the efficient charge separation in photosystem I (PSI) complexes. PSI is a natural light harvesting complex that drives oxygenic photosynthesis in plants, algae, and cyanobacteria. It uses ~300 tightly packed chlorophylls (Chls) to absorb photons, transfer the excitation energy to the reaction center (RC), and generate a charge separated state with near unity quantum efficiency (QE). A better understanding of the mechanism of energy transfer and charge separation in PSI is required for understanding the high QE of natural light harvesting complexes, and it could lead to the further development of artificial photosynthetic systems for solar energy conversion and modification of light harvesting complexes to improve crop yields. We applied two-dimensional optical spectroscopies to different cyanobacterial photosystem I complexes, including PSI complexes that contain Chl f molecules, to map energy transfer pathways and gain insight into the efficient light harvesting of PSI. We used two-dimensional electronic spectroscopies (2DES) to map energy transfer in Chl a and Chl f containing PSI complexes. To investigate the Chl f PSI complexes, we modified our spectrometer to probe the lower energy states associated with Chl f molecules. We interpreted the 2DES spectra through global analysis procedures to generate maps of energy transfer. We also constructed a two-dimensional electronic vibrational (2DEV) spectrometer that will be used to investigate charge transfer transitions and dynamics within PSI complexes. Measurements were performed on model systems to establish general data analysis procedures for interpreting 2D spectra and gain insight into protein cofactor interactions.

14 SOLAR ENERGY↗
DOE OSTI2020

Kinetic Modeling and Experimental Studies of Capacitively Coupled Macro- and Microplasmas at Microwave Frequencies

The effect of driving frequency in the range of 13.56 MHz to 73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct temperatures for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2 eV to 10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25 % of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically-derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗