Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse 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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Controlling inversion disorder in a stoichiometric spinel magnet

In the study of frustrated quantum magnets, it is essential to be able to control the nature and degree of site disorder during the growth process, as many measurement techniques are incapable of distinguishing between site disorder and frustration-induced spin disorder. Pyrochlore-structured spinel oxides can serve as model systems of geometrically frustrated three-dimensional quantum magnets; however, the nature of the magnetism in one well-studied spinel, ZnFe 2 O 4 , remains unclear. Here, we demonstrate simultaneous control of both stoichiometry and inversion disorder in the growth of ZnFe 2 O 4 single crystals, directly yielding a revised understanding of both the collective spin behavior and lattice symmetry. Crystals grown in the stoichiometric limit with minimal site inversion disorder contravene all the previously suggested exotic spin phases in ZnFe 2 O 4 . Furthermore, the structure is confirmed on the F$\bar{4}$3m space group with broken inversion symmetry that induces antiferroelectricity. The effective tuning of magnetic behavior by site disorder in the presence of robust antiferroelectricity makes ZnFe 2 O 4 of special interest to multiferroic devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Bilevel optimization, deep learning and fractional Laplacian regularization with applications in tomography

Here we consider a generalized bilevel optimization framework for solving inverse problems. We introduce fractional Laplacian as a regularizer to improve the reconstruction quality, and compare it with the total variation regularization. We emphasize that the key advantage of using fractional Laplacian as a regularizer is that it leads to a linear operator, as opposed to the total variation regularization which results in a nonlinear degenerate operator. Inspired by residual neural networks, to learn the optimal strength of regularization and the exponent of fractional Laplacian, we develop a dedicated bilevel optimization neural network with a variable depth for a general regularized inverse problem. We illustrate how to incorporate various regularizer choices into our proposed network. As an example, we consider tomographic reconstruction as a model problem and show an improvement in reconstruction quality, especially for limited data, via fractional Laplacian regularization. We successfully learn the regularization strength and the fractional exponent via our proposed bilevel optimization neural network. We observe that the fractional Laplacian regularization outperforms total variation regularization. This is specially encouraging, and important, in the case of limited and noisy data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

DEEPEN: Newberry Volcano MT and Gravity Data - 2022 and 2023 Acquisition and Processing

As part of DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments), a 3D play fairway analysis (PFA) was conducted at Newberry Volcano in Central Oregon for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). For use in this PFA, combined full tensor broadband magnetotelluric (MT) and gravity data were acquired, processed, and inverted by Enthalpion Energy LLC (Enthalpion) with support from NREL staff. The data collection efforts took place from June 19th to July 24. Data were collected with the goal of gaining an improved understanding of the South Flank and the extent of the magma chamber. This GDR submission includes the raw data, single inversions, joint inversions, resolution matrices, and a report on these processes. More detailed information about the folder structure of the datasets included in this submission may be found on pages 94-97 of Magnetotelluric and Gravity Survey Report.pdf below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Forward sensitivity analysis and mode dependent control for closure modeling of Galerkin systems

Model reduction by projection-based approaches is often associated with losing some of the important features that contribute towards the dynamics of the retained scales. As a result, a mismatch occurs between the predicted trajectories of the original system and the truncated one. We put forth a framework to apply a continuous time control signal in the latent space of the reduced order model (ROM) to account for the effect of truncation. We set the control input using parameterized models by following energy transfer principles. Our methodology relies on observing the system behavior in the physical space and using the projection operator to restrict the feedback signal into the latent space. Then, we leverage the forward sensitivity method (FSM) to derive relationships between the feedback and the desired mode-dependent control. We test the performance of the proposed approach using two test cases, corresponding to viscous Burgers and vortex merger problems at high Reynolds number. Results show that the ROM trajectory with the applied FSM control closely matches its target values in both the data-dense and data-sparse regimes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Probabilistic constrained Bayesian inversion for transpiration cooling

Abstract To enable safe operations in applications such as rocket combustion chambers, the materials require cooling to avoid material damage. Here, transpiration cooling is a promising cooling technique. Numerous studies investigate possibilities to simulate and evaluate the complex cooling mechanism. One naturally arising question is the amount of coolant required to ensure a safe operation. To study this, we introduce an approach that determines the posterior probability distribution of the Reynolds number using an inverse problem and constraining the maximum temperature of the system under parameter uncertainties. Mathematically, this chance inequality constraint is dealt with by a generalized polynomial chaos expansion of the system. The posterior distribution will be evaluated by different Markov chain Monte Carlo based methods. A novel method for the constrained case is proposed and tested among others on two‐dimensional transpiration cooling models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A fully kinetic model for orphan gamma-ray flares in blazars

ABSTRACT Blazars emit a highly variable non-thermal spectrum. It is usually assumed that the same non-thermal electrons are responsible for the IR-optical-UV emission (via synchrotron) and the gamma-ray emission (via inverse Compton). Hence, the light curves in the two bands should be correlated. Orphan gamma-ray flares (i.e. lacking a luminous low-frequency counterpart) challenge our theoretical understanding of blazars. By means of large-scale two-dimensional radiative particle-in-cell simulations, we show that orphan gamma-ray flares may be a self-consistent by-product of particle energization in turbulent magnetically dominated pair plasmas. The energized particles produce the gamma-ray flare by inverse Compton scattering an external radiation field, while the synchrotron luminosity is heavily suppressed since the particles are accelerated nearly along the direction of the local magnetic field. The ratio of inverse Compton to synchrotron luminosity is sensitive to the initial strength of turbulent fluctuations (a larger degree of turbulent fluctuations weakens the anisotropy of the energized particles, thus increasing the synchrotron luminosity). Our results show that the anisotropy of the non-thermal particle population is key to modelling the blazar emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Quantitative and three-dimensional assessment of holdup material

Nuclear material deposited in equipment, transfer lines, and ventilation systems of a processing facility is usually referred to as holdup. Operators need to know the location and amount of holdup in their facility for radiation safety and material accountability. In this work, we propose to use an array of detectors co-axial to the inspected pipe to measure the holdup material. This method is implementable into an automated system capable of crawling on surfaces and pipes of various curvatures, which would enable faster, easier, and more accurate holdup safeguards measurements. We first demonstrated that the current holdup assay procedure could lead to a non-negligible bias in the estimate of special nuclear material mass, due to the simplified assumption of deposited geometry introduced by the Generalized Geometry Holdup (GGH) model. The new approach consists of imaging the inner holdup material by characterizing the detector array’s response and unfolding it from the measured light output. Our experimental proof of principle consists of three NaI(Tl) detectors surrounding an aluminum pipe containing two cesium-137 ( 137 Cs) sources. We derived the source distribution inside the pipe by first calculating the detector response matrix using a method adaptive to the surface geometry of the object containing the measured holdup material. Creating a matrix of the detector array’s measured counts, we then proceed to solve an inverse problem, resulting in an accurately located source position and activity distribution within the response matrix’s spatial resolution. Furthermore, we then developed a simulated model of the envisioned experimental setup, which accurately described both the activity and position of the source in 2D. Finally, we extended our model onto a discretized three-dimensional model of the system, encompassing 36 detectors. The model was simulated to measure the source distribution and activity of multiple sources in pipes of varied geometries, which could then be conveniently translated to a real-time holdup measurement protocol. For the 3D simulation of four different source geometries, the model accurately localized the source position in 3D, while the activity retained a maximum relative error of ±5.32%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Dynamic and static corrosion of chromium-alumina refractory in simulated nuclear waste glass

Monofrax® K-3, a chromium–alumina refractory, is widely used in nuclear waste glass melters due to its high durability. However, aggressive vitrification conditions still lead to corrosion that limits melter lifespan. Here, this study investigates two key corrosion modes, subsurface and melt-line corrosion, under static and dynamic conditions, using a custom setup enabling precise control of melt velocity and temperature. Experimental results, interpreted using diffusion-limited dissolution models, show that (i) subsurface corrosion increases by approximately 10% per mm/s increase in melt velocity, while melt-line corrosion remains virtually unaffected, (ii) corrosion rate dependence on temperature is inversely proportional to melt viscosity, and (iii) glass composition affects both corrosion modes similarly, based on Cr 2 O 3 , Al 2 O 3 , and alkali contents. The experimental data also allowed the determination of diffusion coefficients of major dissolving species, which were found to be consistent between the melt-line and subsurface corrosion models, indicating their suitability for future CFD studies.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2023

Relating storm-snow avalanche instabilities to data collected from the Differential Emissivity Imaging Disdrometer (DEID)

Storm-snow avalanches are challenging to forecast due to complex alpine terrain and during rapidly changing weather conditions. They can result in loss of lives and significant economic impact. We describe how a new device that continuously measures with high-frequency snowflake mass, size, density, and type, the Differential Emissivity Imaging Disdrometer (DEID), and show how the DEID can be used to aid avalanche forecasting when coupled with a storm-snow stability model. DEID measurements of snow accumulation, snow water equivalent (SWE), and snow density obtained during seventeen storms taken at the mid-Collins Snow-Study Plot at Alta Ski Area in Utah's Central Wasatch mountain range during winter 2020–2021 show excellent agreement with infrequent manual measurements. Additionally, two new variables, the Shape Density Index (SDI) and Complexity, are proposed and used to classify snowflake habit and estimate storm-snow shear strength. We illustrate how these DEID-derived data can be used to identify layers of concern in the storm snow such as density inversions, in real-time without digging snow pits. Furthermore, the DEID-data are used to run four variations of the SNOw Slope Stability model (SNOSS) for the storms investigated. The results are evaluated with data collected from tilt-board tests, infrasound measurements, and visual observations of avalanches. For a total fourteen storms analyzed, the DEID-driven SNOSS-modeled minimum stability index predicts the general stability of the storm-snow as indicated by observed avalanches, both natural and of unknown cause. Finally, the results provide a promising approach for nowcasting instabilities within storm-snow layers with a single instrument.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

In situ aerobic methane oxidation rates in a stratified lake

Abstract Microbial aerobic methane oxidation is an important sink for aquatic methane worldwide. Despite its importance to global methane fluxes, few aerobic methane oxidation rates have been obtained in freshwater or marine environments without imposing changes to the microbial community through use of ex situ methods. A novel in situ incubation method for continuous time‐series measurements was used in Jordan Lake, North Carolina, during 2020–2021, to determine reaction kinetics for aerobic methane oxidation rates across a wide range of naturally varying methane (55–1833 nM) and dissolved oxygen (DO; 28–366 μM) concentrations and temperatures (17–30°C). Methane oxidation began immediately at the start of each of 21 incubations and methane oxidation rates were 1 st order with respect to methane. The data density allowed for accurate calculation of 1 st ‐order rate constants,k, that ranged from 0.018 to 0.462 h −1 (R 2 > 0.967). Addition of ammonium (20–45 μM) to natural concentrations ranging from 0.057 to 2.4 μM did not change aerobic methane oxidation rate kinetics, suggesting that the natural population of aerobic methane oxidizers in this eutrophic lake was not nitrogen limited. Values ofkinversely correlated most strongly with initial DO concentrations (R 2 = 0.82) rather than temperature. Values forkincreased with Julian day throughout our sampling period, suggesting seasonal influences on methane oxidation via responses to geochemical changes or shifts in microbial community abundance and composition. These experiments demonstrate a high variability in the enzymatic capacity for 1 st ‐order methane oxidation rates in this eutrophic lake that is tightly and inversely coupled to oxygen concentrations. Measurements of in situ aerobic methane oxidation rate constants allow for the direct quantification and modeling of the microbial community's capacity for methane oxidation over a wide range of natural methane concentrations.

Marine & Freshwater Biology↗
DOE OSTI2025

The Efficiency of Water Vapor on Top‐of‐Atmosphere Radiation

Earth's climate sensitivity is greatly affected by the compensation between temperature feedback and water vapor (WV) feedback. Using abrupt 4xCO 2 experiments, we show that the global-mean WV feedback is nearly a linear function of the temperature feedback, the slope of which is explained by the longwave radiative efficiency of WV (ϵ). Although ϵ remains constant across models in the global mean, it exhibits substantial spatial variations and is particularly weak in Antarctica, where near-surface inversions decouple the surface from the free troposphere. We introduce a surface–free troposphere temperature difference (SFTD) metric, showing that positive SFTD (e.g., high lifting condensation level) amplifies ϵ, while negative SFTD (e.g., strong surface inversion) suppresses it. These findings provide a clear explanation of how local climate conditions modulate the radiative compensation between temperature and WV feedbacks.

climate feedback↗
DOE OSTI2024

Source size of x rays from self-modulated laser wakefield accelerators

A comparative study of x-ray sources generated with different mechanisms from self-modulated laser wakefield acceleration (SM-LWFA) electrons was performed to compare the source size or spatial resolution for use in high energy density science applications. We examine the source size of betatron, inverse Compton scattering, and bremsstrahlung radiation with a Fresnel diffraction based formalism and a modified x-ray ray tracing model. We observe the dependence of source size on the radiation generation process, laser parameters, and compare to what is possible in other regimes of LWFA, as well as current methods. This information is significant as we begin to explore the use of light sources driven by SM-LWFA for use as a diagnostic at large-scale laser facilities where blowout regime LWFA is not possible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time-resolved site-selective imaging of predissociation and charge transfer dynamics: the CH 3 I B-band

The predissociation dynamics of the 6s (B 2 E) Rydberg state of gas-phase CH 3 I were investigated by time-resolved Coulomb-explosion imaging using extreme ultraviolet (XUV) free-electron laser pulses. Inner-shell ionization at the iodine 4d edge was utilized to provide a site-specific probe of the ensuing dynamics. The combination of a velocity-map imaging (VMI) spectrometer coupled with the pixel imaging mass spectrometry (PImMS) camera permitted three-dimensional ionic fragment momenta to be recorded simultaneously for a wide range of iodine charge states. In accord with previous studies, initial excitation at 201.2 nm results in internal conversion and subsequent dissociation on the lower-lying A-state surface on a picosecond time scale. Examination of the time-dependent yield of low kinetic energy iodine fragments yields mechanistic insights into the predissociation and subsequent charge transfer following multiple ionization of the iodine products. The effect of charge transfer was observed through differing delay-dependencies of the various iodine charge states, from which critical internuclear distances for charge transfer could be inferred and compared to a classical over-the-barrier model. Time-dependent photofragment angular anisotropy parameters were extracted from the central slice of the Newton sphere, without Abel inversion, and highlight the effect of rotation of the parent molecule before dissociation, as observed in previous works. Our results demonstrate the ability to perform three-dimensional ion imaging at high event rates and showcase the potential benefits of this approach, particularly in relation to further time-resolved studies at free-electron laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗
DOE OSTI2021

Validation of the HyRAM+ physics models for use with propane

The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses

08 HYDROGEN↗
DOE OSTI2020

NuSTAR Perspective on High-redshift MeV Blazars

With bolometric luminosities exceeding 10 48 erg s -1 , powerful jets, and supermassive black holes at their center, MeV blazars are some of the most extreme sources in the universe. Recently, the Fermi-Large Area Telescope detected five new γ-ray emitting MeV blazars beyond redshift z = 3.1. With the goal of precisely characterizing the jet properties of these extreme sources, we started a multiwavelength campaign to follow them up with joint Nuclear Spectroscopic Telescope Array, Swift, and the Southeastern Association for Research in Astronomy’s optical telescopes. We observe six high-redshift quasars, four of them belonging to the new γ-ray emitting MeV blazars. Thorough X-ray analysis reveals spectral flattening at soft X-ray for three of these objects. The source NVSS J151002+570243 also shows a peculiar rehardening of the X-ray spectrum at energies E > 6 keV. Adopting a one-zone leptonic emission model, this combination of hard X-rays and γ-rays enables us to determine the location of the Inverse Compton peak and to accurately constrain the jet characteristics. In the context of the jet-accretion disk connection, we find that all six sources have jet powers exceeding accretion disk luminosity, seemingly validating this positive correlation even beyond z > 3. Our six sources are found to have 10 9 M ⊙ black holes, further raising the space density of supermassive black holes in the redshift bin z = [3, 4].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Sequential spectral line analysis for accurate density and temperature diagnosis of laboratory opacity measurements

The accuracy of iron opacity calculated in stellar interiors has been questioned since the discovery of the “solar problem” and the discrepancies between the measured and modeled iron opacity reported in 2015. Experimental opacity benchmarks require accurate temperature and density measurements, which were inferred by analyzing tracer magnesium spectra in those experiments. Could the observed discrepancy be explained by insufficient accuracy in the inferred temperature, density, and their uncertainties? Previous analyses may have yielded biased results due to three limitations: (1) simultaneous multi-line fitting, (2) approximations in line-shape models, and (3) exclusion of certain spectral lines due to insufficient background characterization. Notably, the first issue is a common concern for many inversion methods, including Bayesian inferences. We present a refined analysis method that overcomes these limitations, applied to three categories of iron opacity experiments (Anchor 1, 2, and 3). In particular, the sequential fitting method yields unbiased results with more realistic uncertainties by accounting for line inconsistencies in the parameter uncertainties. The average electron temperature and density values are 162 ± 6 eV and (7.0 ± 1.9) × 10 21 cm −3 for six Anchor 1 experiments, 189 ± 7 eV and (3.4 ± 0.3) × 10 22 cm −3 for 21 Anchor 2 experiments, and 201 ± 6 eV and (4.8 ± 1.1) × 10 22 cm −3 for nine Anchor 3 experiments. These results show ∼4% temperature and ∼20% density reproducibility over a decade, which also aligns with the inferred parameter uncertainties. In conclusion, the resulting temperature and density uncertainties lead to a quasi-continuum iron opacity variation of ±4%–7% for wavelengths below 9.5 Å, which is insufficient to explain the significant model-data discrepancies reported in 2015.

Absorption spectroscopy↗
DOE OSTI2023

Machine learning from RANS and LES to inform coarse grid simulations

Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗