Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Observation and Simulation of Methane Plumes During the Morning Boundary Layer Transition

Abstract Methane (CH 4 ) contributes significantly to global warming. However, accurate identification of CH 4 sources for reducing CH 4 emissions is often hampered by inadequate accuracy and spatiotemporal coverage of CH 4 detection, and lack of accurate CH 4 forward modeling used in top‐down inversion systems. In this study, a field experiment was conducted in Pampa, Texas using two CH 4 sensors (LI‐COR and OGI camera) to detect CH 4 releases. We investigated whether high‐resolution simulations using the Weather Research and Forecasting (WRF) model with greenhouse gases (WRF‐GHG) could accurately simulate the CH 4 plumes in the presence of evolving atmospheric boundary layer from sunrise to noon. CH 4 plumes showed substantial variation in time. At a release rate of ∼17.5 kg hr −1 , the maximum enhancement of CH 4 measured by LI‐COR was 2.6 ppm at sunrise (7:36 a.m.), 250 m from the release location. Within half an hour after sunrise, this enhancement decreased to 0.3–0.4 ppm. The enhancement was 0.2 ppm by 10:00 a.m. and further dropped to less than 0.1 ppm after 11:30 a.m. Due to the low temperature at sunrise, the OGI camera failed to detect the CH 4 plume. The WRF‐GHG large‐eddy simulation (LES) with 32 m grid spacing successfully reproduced these CH 4 enhancements. In situ measurements together with numerical simulations illustrate the impact of the transition from a stable boundary layer in the early morning to a convective boundary layer at noon on the dispersion of CH 4 plumes. Additionally, CH 4 plumes from a cattle farm in Oklahoma are briefly examined using the same modeling approach.

Hu, Xiao‐Ming [Center for Analysis and Prediction ↗
DOE OSTI2026

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Recent advances and applications of deep learning methods in materials science

Deep learning (DL) is one of the fastest-growing topics in materials data science, with rapidly emerging applications spanning atomistic, image-based, spectral, and textual data modalities. DL allows analysis of unstructured data and automated identification of features. The recent development of large materials databases has fueled the application of DL methods in atomistic prediction in particular. In contrast, advances in image and spectral data have largely leveraged synthetic data enabled by high-quality forward models as well as by generative unsupervised DL methods. In this article, we present a high-level overview of deep learning methods followed by a detailed discussion of recent developments of deep learning in atomistic simulation, materials imaging, spectral analysis, and natural language processing. For each modality we discuss applications involving both theoretical and experimental data, typical modeling approaches with their strengths and limitations, and relevant publicly available software and datasets. We conclude the review with a discussion of recent cross-cutting work related to uncertainty quantification in this field and a brief perspective on limitations, challenges, and potential growth areas for DL methods in materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2024

AI-enabled Lorentz microscopy for quantitative imaging of nanoscale magnetic spin textures

The manipulation and control of nanoscale magnetic spin textures are of rising interest as they are potential foundational units in next-generation computing paradigms. Achieving this requires a quantitative understanding of the spin texture behavior under external stimuli using in situ experiments. Lorentz transmission electron microscopy (LTEM) enables real-space imaging of spin textures at the nanoscale, but quantitative characterization of in situ data is extremely challenging. Here, we present an AI-enabled phase-retrieval method based on integrating a generative deep image prior with an image formation forward model for LTEM. Our approach uses a single out-of-focus image for phase retrieval and achieves significantly higher accuracy and robustness to noise compared to existing methods. Furthermore, our method is capable of isolating sample heterogeneities from magnetic contrast, as shown by application to simulated and experimental data. This approach allows quantitative phase reconstruction of in situ data and can also enable near real-time quantitative magnetic imaging.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantitative imaging and automated fuel pin identification for passive gamma emission tomography

Compliance of member States to the Treaty on the Non-Proliferation of Nuclear Weapons is monitored through nuclear safeguards. The Passive Gamma Emission Tomography (PGET) system is a novel instrument developed within the framework of the International Atomic Energy Agency (IAEA) project JNT 1510, which included the European Commission, Finland, Hungary and Sweden. The PGET is used for the verification of spent nuclear fuel stored in water pools. Advanced image reconstruction techniques are crucial for obtaining high-quality cross-sectional images of the spent-fuel bundle to allow inspectors of the IAEA to monitor nuclear material and promptly identify its diversion. In this work, we have developed a software suite to accurately reconstruct the spent-fuel cross sectional image, automatically identify present fuel rods, and estimate their activity. Unique image reconstruction challenges are posed by the measurement of spent fuel, due to its high activity and the self-attenuation. While the former is mitigated by detector physical collimation, we implemented a linear forward model to model the detector responses to the fuel rods inside the PGET, to account for the latter. The image reconstruction is performed by solving a regularized linear inverse problem using the fast-iterative shrinkage-thresholding algorithm. We have also implemented the traditional filtered back projection (FBP) method based on the inverse Radon transform for comparison and applied both methods to reconstruct images of simulated mockup fuel assemblies. Higher image resolution and fewer reconstruction artifacts were obtained with the inverse-problem approach, with the mean-square-error reduced by 50%, and the structural-similarity improved by 200%. We then used a convolutional neural network (CNN) to automatically identify the bundle type and extract the pin locations from the images; the estimated activity levels finally being compared with the ground truth. The proposed computational methods accurately estimated the activity levels of the present pins, with an associated uncertainty of approximately 5%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Electron energy loss spectroscopy database synthesis and automation of core-loss edge recognition by deep-learning neural networks

Abstract The ionization edges encoded in the electron energy loss spectroscopy (EELS) spectra enable advanced material analysis including composition analyses and elemental quantifications. The development of the parallel EELS instrument and fast, sensitive detectors have greatly improved the acquisition speed of EELS spectra. However, the traditional way of core-loss edge recognition is experience based and human labor dependent, which limits the processing speed. So far, the low signal–noise ratio and the low jump ratio of the core-loss edges on the raw EELS spectra have been challenging for the automation of edge recognition. In this work, a convolutional-bidirectional long short-term memory neural network (CNN-BiLSTM) is proposed to automate the detection and elemental identification of core-loss edges from raw spectra. An EELS spectral database is synthesized by using our forward model to assist in the training and validation of the neural network. To make the synthesized spectra resemble the real spectra, we collected a large library of experimentally acquired EELS core edges. In synthesize the training library, the edges are modeled by fitting the multi-Gaussian model to the real edges from experiments, and the noise and instrumental imperfectness are simulated and added. The well-trained CNN-BiLSTM network is tested against both the simulated spectra and real spectra collected from experiments. The high accuracy of the network, 94.9%, proves that, without complicated preprocessing of the raw spectra, the proposed CNN-BiLSTM network achieves the automation of core-loss edge recognition for EELS spectra with high accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding Peelle’s Pertinent Puzzle bias in generalized least squares regression through eigenspectrum analysis

Certain correlation structures in the data covariance matrix (DCM) used for generalized least squares (GLS) regression can result in biased estimates, commonly known in the field of nuclear data evaluation as Peele’s Pertinent Puzzle (PPP). This article introduces a generative, forward modeling framework within which the PPP bias is characterized through an eigenspectrum analysis of the DCM. This analysis highlights the root cause of the bias, generalizes the problem beyond the nuclear data field, and provides insight to the problem regimes where it can occur. What follows is an understanding that the bias can show up for any experimental neutron time-of-flight data for which systematic uncertainties have been quantified. Lastly, a discussion of the adaptation of cross validation approaches that require pre-whitening to incorporate the known ‘fix’ to the PPP bias in the GLS estimator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Machine learning inversion of interatomic force constants from single-crystal inelastic neutron scattering

Atomic vibrations govern many macroscopic properties of materials, but experiments to comprehensively probe them remain challenging. Inelastic neutron scattering (INS) is a powerful technique to map phonon dispersions in crystals, especially when leveraging modern time-of-flight (ToF) spectrometers with large detectors. However, efficiently and robustly extracting interatomic force constants (FCs) parameterizing phonon dynamics from experimental spectra remains a bottleneck due to the complexity and high dimensionality of ToF INS datasets. Here, we present a machine learning approach for the direct inversion of FCs from single-crystal INS measurements. The framework leverages synthetic training data generated using universal machine-learned force fields and an efficient physics-based forward model. We benchmark two neural architectures–one emphasizing structured latent representation learning and the other direct, supervised spectral regression–across simulated datasets for two materials under idealized and noisy conditions. The latent-representation model is subsequently applied to experimental single-crystal INS data on germanium. The model is shown to reproduce FCs derived from both first-principles simulations and from iterative optimization, and furthermore achieves reliable inference even from sparse, single-orientation measurements representing short data acquisitions. Analysis of the learned latent space reveals semantically continuous and physically interpretable encodings that support strong cross-domain generalization. By bridging theoretical and experimental domains, we establish a path toward rapid inversion of experimental spectra and data-driven interpretation of temperature-dependent lattice dynamics.

42 ENGINEERING↗
DOE OSTI2021

Simultaneous inference of the compressibility and inelastic response of tantalum under extreme loading

We study the deformation of tantalum under extreme loading conditions. Experimental velocity data are drawn from both ramp loading experiments on Sandia’s Z-machine and gas gun compression experiments. The drive conditions enable the study of materials under pressures greater than 100 GPa. We provide a detailed forward model of the experiments including a model of the magnetic drive for the Z-machine. Utilizing these experiments, we simultaneously infer several different types of physically motivated parameters describing equation of state, plasticity, and anelasticity via the computational device of Bayesian model calibration. Characteristics of the resulting calculated posterior distributions illustrate relationships among the parameters of interest via the degree of cross correlation. The calibrated velocity traces display good agreement with the experiments up to experimental uncertainty as well as improvement over previous calibrations. Examining the Z-shots and gun-shots together and separately reveals a trade-off between accuracy and transferability across different experimental conditions. Implications for model calibration, limitations from model form, and suggestions for improvements are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Diagnosing the pedestal magnetic field and magnetohydrodynamics radial structure with pedestal–scrape of layer electron cyclotron emission radiation inversion in H-mode plasma (invited)

Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design study of an edge current density diagnostic using new high-performance single-channel beam emission spectrometers at DIII-D

A novel Motional Stark Effect spectroscopy system has been designed for application at the DIII-D tokamak. The system is optimized for studies of the poloidal and toroidal magnetic field in the plasma pedestal region with frame rates of up to 10 kHz. Further, light from an existing high-photon-throughput collection lens is analyzed using four single-channel f/2.8 Czerny–Turner spectrometers that use custom-made lens systems instead of mirrors. Each spectrometer has two separate outgoing legs and is operated in a positive grating order, which allows for simultaneous observations of D-alpha and D-beta spectra. Forward modeling using the code FIDASIM shows a radial resolution of the system close to 0.6 cm and sufficiently good spectral resolution when masking the high throughput light collection lens in the horizontal direction to avoid overly strong Doppler broadening of beam emission lines. Moreover, a detailed sensitivity study considering realistic levels of readout and photon noise shows that the poloidal and toroidal magnetic field strengths can be inferred with an uncertainty of less than 1%, which will allow the inference of changes of the plasma current during transient events.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Pedestal fluctuation measurements with charge exchange imaging at the DIII-D tokamak

A new high radial resolution 2D multichannel Charge eXchange Imaging (CXI) diagnostic is under development for deployment at DIII-D. Here, the diagnostic system will measure low-to-intermediate radial wavenumber carbon density fluctuations by observing the n = 8 – 7 (λ = 529.06 nm) C-VI emission line, resulting from charge exchange collisions between heating neutral beam atoms and the intrinsic carbon ion density. The new CXI diagnostic will provide measurements with ΔR ~ 0.4 cm to access higher k r instabilities (k r < 8 cm –1 ) predicted to arise in the steep-gradient region of the H-mode pedestal. The CXI system will feature 60 fiber bundles in a 12 × 5 arrangement, with each bundle consisting of four 1 mm fibers. A custom optical system has been designed to filter and image incoming signals onto an 8 × 8 avalanche photodiode array. Additionally, a novel electronics suite has been designed and commissioned to amplify and digitize the relatively low-intensity carbon signal at a 2 MHz bandwidth. Forward modeling results of the active C-VI emission suggest sufficient signal to noise ratios to resolve turbulent fluctuations. Prototype measurements demonstrate the ability to perform high frequency pedestal measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Toward streaked collective Thomson scattering measurements on an extreme ultraviolet plasma light source

Here we show through forward modeling calculations that streaked collective Thomson scattering measurements are feasible on laser-produced tin plasmas generated under conditions relevant for extreme ultraviolet lithography. Using a 532 nm probe laser beam, the feasibility of simultaneous measurements of electron plasma wave (EPW) and ion acoustic wave (IAW) spectra is investigated. Absolute photon counts for laser scattering off both waves are calculated. Probe laser electron heating and bremsstrahlung background radiation effects are accounted for. While a large spatiotemporal region can be successfully probed based on the IAW feature, only one measurement location can be accessed through the EPW as a result of the low signal to noise ratio. A portable/traveling tabletop system is proposed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Recovering non-Maxwellian particle velocity distribution functions from collective Thomson-scattered spectra

Collective optical Thomson scattering (TS) is a diagnostic commonly used to characterize plasma parameters. These parameters are typically extracted by a fitting algorithm that minimizes the difference between a measured scattered spectrum and an analytic spectrum calculated from the velocity distribution function (VDF) of the plasma. However, most existing TS analysis algorithms assume that the VDFs are Maxwellian, and applying an algorithm that makes this assumption does not accurately extract the plasma parameters of a non-Maxwellian plasma due to the effect of non-Maxwellian deviations on the TS spectra. We present new open-source numerical tools for forward modeling analytic spectra from arbitrary VDFs and show that these tools are able to more accurately extract plasma parameters from synthetic TS spectra generated by non-Maxwellian VDFs compared to standard TS algorithms. Estimated posterior probability distributions of fits to synthetic spectra for a variety of example non-Maxwellian VDFs are used to determine uncertainties in the extracted plasma parameters and show that correlations between parameters can significantly affect the accuracy of fits in plasmas with non-Maxwellian VDFs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Optimization of the gamma reaction history diagnostic for double-shell pusher areal density and reaction history measurements on the National Ignition Facility

The double-shell inertial confinement fusion campaign, which consists of an aluminum ablator, a foam cushion, a high-Z pusher (tungsten or molybdenum), and liquid deuterium–tritium (DT) fuel, aims for its first DT filled implosions on the National Ignition Facility (NIF) in 2024. The high-Z, high density pusher does not allow x-rays to escape the double-shell capsule. Therefore, nuclear diagnostics such as the Gamma Reaction History (GRH) diagnostic on the NIF are crucial for understanding high-Z implosion performance. To optimize the GRH measurement of fusion reaction history and the pusher’s areal density, the MCNP6.3-based forward model of the detector was built. When calculating the neutron-induced inelastic gamma ray production, the interaction of neutrons with the compressed fuel was additionally included. By folding the calculated gamma ray spectrum output and the previously calibrated GRH detector responses, the optimum set of GRH energy thresholds for measuring the pusher areal density is determined to be 2.9 and 6.3 MeV for DT double-shell experiments. In addition, the effect of the down-scattering of neutrons on the gamma ray spectrum, the minimum required yield for measurements, and the attenuation of the gamma rays through the pusher are analyzed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Effects of mosaic crystal instrument functions on x-ray Thomson scattering diagnostics

Mosaic crystals, with their high integrated reflectivities, are widely employed in spectrometers used to diagnose high energy density systems. X-ray Thomson scattering (XRTS) has emerged as a powerful diagnostic tool of these systems, providing in principle direct access to important properties such as the temperature via detailed balance. However, the measured XRTS spectrum is broadened by the spectrometer instrument function (IF), and without careful consideration of the IF one risks misdiagnosing system conditions. Here, we consider in detail the IF of 40 and 100 μm mosaic Highly Annealed Pyrolytic Graphite crystals, and how the broadening varies across the spectrometer in an energy range of 6.7–8.6 keV. Notably, we find a strong asymmetry in the shape of the IF toward higher energies. As an example, we consider the effect of the asymmetry in the IF on the temperature inferred via XRTS for simulated 80 eV CH plasmas and find that the temperature can be overestimated if an approximate symmetric IF is used. We, therefore, expect a detailed consideration of the full IF will have an important impact on system properties inferred via XRTS in both forward modeling and model-free approaches.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A compact x-ray spectrometer for measurements of electron temperature distributions in inertial confinement fusion implosions at OMEGA

The Wedge Range Filter (WRF), commonly used for proton spectroscopy at the OMEGA Laser Facility and National Ignition Facility, is adapted to measure the x-ray continuum spectrum through transmission measurement using a continuous-gradient filter. Continuum x rays emitted from the hotspot of an implosion contain information about the plasma composition and electron temperature. The WRF data are leveraged to probe this distribution, specifically the electron temperature distribution. In this work, the data recorded with the WRF are forward modeled using a temperature distribution model folded with the WRF response function. An uncertainty analysis is conducted through a Bayesian regression algorithm using a Hamiltonian Monte Carlo sampler. This analysis enables the uncertainties in the instrument response to be folded into the uncertainty estimation of the electron temperature and absolute x-ray emission. Data analysis for a series of OMEGA implosions is presented and compared with radiation hydrodynamic simulations.

Lasers↗
DOE OSTI2025

Reconstruction of 2D line-integrated electron density using angular filter refractometry and a fast marching Eikonal solver

Refraction of an optical probe beam by a plasma can be measured with angular filter refractometry (AFR), which produces an image of the beam’s 2D spatial profile that contains intensity contours corresponding to curves of constant refraction angle. Further analysis is required to reconstruct the underlying line-integrated electron density. Most prior efforts to calculate density from AFR data have been limited to 1D analysis or forward-fitting techniques. Here, in this paper, we detail the use of a fast-marching Eikonal solver to directly invert AFR data and obtain the full 2D line-integrated electron density. The analysis method is first verified with synthetic data and then applied to experimental measurements of single and colliding plasma plumes collected at the OMEGA EP Laser Facility. The calculated densities agree with 1D results and are shown to be consistent with the original AFR measurements via forward modeling. We also discuss ways to improve the precision of this technique.

McCluskey, B. [Princeton Univ., NJ (United States)↗