Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse prediction”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Machine learning based approach to predict ductile damage model parameters for polycrystalline metals

Damage models for ductile materials typically need to be parameterized, often with the appropriate parameters changing for a given material depending on the loading conditions. This can make parameterizing these models computationally expensive, since an inverse problem must be solved for each loading condition. Using standard inverse modeling techniques typically requires hundreds or thousands of high-fidelity computer simulations to estimate the optimal parameters. Additionally, the time of a human expert is required to set up the inverse model. Machine learning has recently emerged as an alternative approach to inverse modeling in these settings, where the machine learning model is trained in an offline manner and new parameters can be quickly generated on the fly, after training is complete. Here, this work utilizes such a workflow to enable the rapid parameterization of a ductile damage model called TEPLA with a machine learning inverse model. The machine learning model can efficiently estimate the model parameters much faster, as compared to previously employed methods, such as Bayesian calibration. The results demonstrate good accuracy on a synthetic test dataset and is validated against experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient mapping between void shapes and stress fields using Deep Convolutional Neural Networks with sparse data

Establishing fast and accurate structure-to-property relationships is an important component in the design and discovery of advanced materials. Physics-based simulation models like the finite element method (FEM) are often used to predict deformation, stress, and strain fields as a function of material microstructure in material and structural systems. Such models may be computationally expensive and time intensive if the underlying physics of the system is complex. This limits their application to solve inverse design problems and identify structures that maximize performance. In such scenarios, surrogate models are employed to make the forward mapping computationally efficient to evaluate. However, the high dimensionality of the input microstructure and the output field of interest often renders such surrogate models inefficient, especially when dealing with sparse data. Deep convolutional neural network (CNN) based surrogate models have shown great promise in handling such high-dimensional problems. In this paper, a single ellipsoidal void structure under a uniaxial tensile load represented by a linear elastic, high-dimensional and expensive-to-query, FEM model. We consider two deep CNN architectures, a modified convolutional autoencoder framework with a fully connected bottleneck and a UNet CNN, and compare their accuracy in predicting the von Mises stress field for any given input void shape in the FEM model. Additionally, a sensitivity analysis study is performed using the two approaches, where the variation in the prediction accuracy on unseen test data is studied through numerical experiments by varying the number of training samples from 20 to 100.

surrogate modeling; convolutional neural networks;↗
DOE OSTI2020

Fine Root and Soil Organic Carbon Depth Distributions are Inversely Related Across Fertility and Rainfall Gradients in Lowland Tropical Forests

Humid tropical forests contain some of the largest soil organic carbon (SOC) stocks on Earth. Much of this SOC is in subsoil, yet variation in the distribution of SOC through the soil profile remains poorly characterized across tropical forests. We used a correlative approach to quantify relationships among depth distributions of SOC, fine root biomass, nutrients and texture to 1 m depths across 43 lowland tropical forests in Panama. The sites span rainfall and soil fertility gradients, and these are largely uncorrelated for these sites. We used fitted b parameters to characterize depth distributions, where b is a numerical index based on an asymptotic relationship, such that larger b values indicate greater concentrations of root biomass or SOC at depth in the profile. Root b values ranged from 0.82 to 0.95 and were best predicted by soil pH and extractable potassium (K) stocks. For example, the three most acidic (pH < 4) and Kpoor (< 20 g K m -2 ) soils contained 76 ± 5% of fine root biomass from 0 to 10 cm depth, while the three least acidic (pH > 6.0) and most K-rich (> 50 g K m -2 ) soils contained only 41 ± 9% of fine root biomass at this depth. Root b and SOC b values were inversely related, such that a large fine root biomass in surface soils corresponded to large SOC stocks in subsoils (50–100 cm). SOC b values were best predicted by soil pH and base cation stocks, with the three most base-poor soils containing 34 ± 8% of SOC from 50 to 100 cm depth, and the three most base-rich soils containing just 9 ± 2% of SOC at this depth. Nutrient depth distributions were not related to Root b or SOC b values. Here, these data show that large surface root biomass stocks are associated with large subsoil C stocks in strongly weathered tropical soils. Further studies are required to evaluate why this occurs, and whether changes in surface root biomass, as may occur with global change, could in turn influence SOC storage in tropical forest subsoils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Autonomous Synthesis and Inverse Design of Electrochromic Polymers with High Efficiency and Accuracy

Here, the design and synthesis of functional polymers, aimed at targeted properties through specific structures, have long been challenged by their complex and often nonlinear structure–property relationships. Key processes, including knowledge accumulation for predictive design and experimental refinement and validation, are traditionally labor-insensitive and time-consuming, making it difficult to balance accuracy and efficiency. Here, we introduce an accelerated, autonomous system for the on-demand synthesis of electronic polymers that achieves the desired electrochromic functionality with high accuracy and efficiency. Our approach leverages large language model-assisted data mining, a physics-informed copolymer machine learning model, and an AI-driven autonomous robotic workflow in the Polybot lab. Within 72 h, Polybot autonomously synthesized electrochromic polymers (ECPs) with targeted, previously-unreported color values, including green polymers with specific absorption profiles, precisely fine-tuning copolymer structures with a 5% step size in comonomer composition within a three-monomer system. A publicly accessible ECP informatics database has also been created to foster knowledge exchange.

AI-driven Robotic Lab↗
DOE OSTI2025

Artificial intelligence driven laser parameter search: Inverse design of photonic surfaces using greedy surrogate-based optimization

Photonic surfaces designed with specific optical characteristics are becoming increasingly crucial for novel energy harvesting and storage systems. The design of these surfaces can be achieved by texturing materials using lasers. The optimal adjustment of laser fabrication parameters to achieve target surface optical properties is an open challenge. Thus, we develop a surrogate-based optimization approach. Our framework employs the Random Forest algorithm to model the forward relationship between the laser fabrication parameters and the resulting optical characteristics. During the optimization process, we use a greedy, prediction-based exploration strategy that iteratively selects batches of laser parameters to be used in experimentation by minimizing the predicted discrepancy between the surrogate model’s outputs and the user-defined target optical characteristics. This strategy allows for efficient identification of optimal fabrication parameters without the need to model the error landscape directly. We demonstrate the efficiency and effectiveness of our approach on two synthetic benchmarks and two specific experimental applications of photonic surface inverse design targets. By calculating the average performance of our algorithm compared to other state of the art optimization methods, we show that our algorithm performs, on average, twice as well across all benchmarks. Additionally, a warm starting inverse design technique for changed target optical characteristics enhances the performance of the introduced approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Updates to the Regional Seismic Travel Time (RSTT) Model: 1. Tomography

Abstract A function of global monitoring of nuclear explosions is the development of Earth models for predicting seismic travel times for more accurate calculation of event locations. Most monitoring agencies rely on fast, distance-dependent one-dimensional (1D) Earth models to calculate seismic event locations quickly and in near real-time. RSTT (Regional Seismic Travel Time) is a seismic velocity model and computer software package that captures the major effects of three-dimensional crust and upper mantle structure on regional seismic travel times, while still allowing for fast prediction speed (milliseconds). We describe updates to the RSTT model using a refined data set of regional phases (i.e., Pn, Pg, Sn, Lg) using the Bayesloc relative relocation algorithm. The tomographic inversion shown here acts to refine the previous RSTT public model ( rstt201404um ) and displays significant features related to areas of global tectonic complexity as well as further reduction in arrival residual values. Validation of the updated RSTT model demonstrates significant reduction in median epicenter mislocation (15.3 km) using all regional phases compared to the iasp91 1D model (22.1 km) as well as to the current station correction approach used at the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (18.9 km).

58 GEOSCIENCES↗
DOE OSTI2020

Deep learning of interface structures from simulated 4D STEM data: cation intermixing vs. roughening ∗

Abstract Interface structures in complex oxides remain an active area of condensed matter physics research, largely enabled by recent advances in scanning transmission electron microscopy (STEM). Yet the nature of the STEM contrast in which the structure is projected along the given direction precludes separation of possible structural models. Here, we utilize deep convolutional neural networks (DCNN) trained on simulated 4D STEM datasets to predict structural descriptors of interfaces. We focus on the widely studied interface between LaAlO 3 and SrTiO 3 , using dynamical diffraction theory and leveraging high performance computing to simulate thousands of possible 4D STEM datasets to train the DCNN to learn properties of the underlying structures on which the simulations are based. We test the DCNN on simulated data and show that it is possible (with >95% accuracy) to identify a physically rough from a chemically diffuse interface and create a DCNN regression model to predict step positions. We quantify the applicability of the model to different thicknesses and the transferability of the approach. The method shown here is general and can be applied for any inverse imaging problem where forward models are present.

42 ENGINEERING↗
NASA NTRS2009

A Simple and Efficient Computational Approach to Chafed Cable Time-Domain Reflectometry Signature Prediction

A method for the prediction of time-domain signatures of chafed coaxial cables is presented. The method is quasi-static in nature, and is thus efficient enough to be included in inference and inversion routines. Unlike previous models proposed, no restriction on the geometry or size of the chafe is required in the present approach. The model is validated and its speed is illustrated via comparison to simulations from a commercial, three-dimensional electromagnetic simulator.

Kowalski, Marc Edward↗
DOE OSTI2026

Nuclear Data Adjustment for Nonlinear Applications in the OECD/NEA WPNCS SG14 Benchmark—A Bayesian Inverse UQ-Based Approach for Data Assimilation

The Organisation for Economic Co-operation and Development Working Party on Nuclear Criticality Safety has proposed a benchmark exercise to assess the performance of current nuclear data adjustment techniques applied to nonlinear applications and experiments with low correlation to applications. This work introduces Bayesian inverse uncertainty quantification (IUQ) employing scientific machine learning surrogate models as a method for nuclear data adjustments in this benchmark, and compares IUQ to the more traditional methods of generalized linear least squares (GLLS) and Monte Carlo Bayes (MOCABA). Posterior predictions from IUQ showed agreement with GLLS and MOCABA for linear applications. Here, when comparing GLLS, MOCABA, and IUQ posterior predictions to computed model responses using adjusted parameters, we observe that the GLLS predictions failed to replicate the computed response distributions for nonlinear applications, while MOCABA showed near agreement, and IUQ used the computed model responses directly. We also discuss observations on why experiments with low correlation to applications can be informative to nuclear data adjustments and identify some properties useful in selecting experiments for inclusion in nuclear data adjustment. Performance in this benchmark indicates potential for Bayesian IUQ in nuclear data adjustments.

Bayesian calibration↗
DOE OSTI2021

B-PINNs: Bayesian Physics-informal Neural Networks for Forward and Inverse PDE Problems with Noisy Data

We propose a Bayesian physics-informed neural network (B-PINN) to solve both forward and inverse nonlinear problems described by partial differential equations (PDEs) and noisy data. In this Bayesian framework, the Bayesian neural network (BNN) combined with a PINN for PDEs serves as the prior while the Hamiltonian Monte Carlo (HMC) or the variational inference (VI) could serve as an estimator of the posterior. B-PINNs make use of both physical laws and scattered noisy measurements to provide predictions and quantify the aleatoric uncertainty arising from the noisy data in the Bayesian framework. Compared with PINNs, in addition to uncertainty quantification, B-PINNs obtain more accurate predictions in scenarios with large noise due to their capability of avoiding overfitting. We conduct a systematic comparison between the two different approaches for the B-PINNs posterior estimation (i.e., HMC or VI), along with dropout used for quantifying uncertainty in deep neural networks. Our experiments show that HMC is more suitable than VI with mean field Gaussian approximation for the B-PINNs posterior estimation, while dropout employed in PINNs can hardly provide accurate predictions with reasonable uncertainty. Finally, we replace the BNN in the prior with a truncated Karhunen-Loève (KL) expansion combined with HMC or a deep normalizing flow (DNF) model as posterior estimators. The KL is as accurate as BNN and much faster but this framework cannot be easily extended to high-dimensional problems unlike the BNN based framework.

Non linear PDEs, Noisy data, Bysian physics inform↗
DOE OSTI2021

The Highest Energy HAWC Sources are Likely Leptonic and Powered by Pulsars

The HAWC Collaboration has observed gamma rays at energies above 56 TeV from a collection of nine sources. It has been suggested that this emission could be hadronic in nature, requiring that these systems accelerate cosmic-ray protons or nuclei up to PeV-scale energies. In this paper, we instead show that the spectra of these objects favor a leptonic (inverse Compton) origin for their emission. More specifically, the gamma-ray emission from these objects can be straightforwardly accommodated within a model in which $\sim \mathcal{O}(10\%)$ of the host pulsar's spindown power is transferred into the acceleration of electrons and positrons with a power-law spectrum that extends to several hundred TeV or higher. The spectral break that is observed among these sources is naturally explained within the context of this simple model, and occurs at the energy where the timescale for energy losses matches the age of the pulsar. In contrast, this spectral feature cannot be straightforwardly accommodated in hadronic scenarios. Furthermore, hadronic models predict that these sources should produce more emission at GeV-scale energies than is observed. In light of these considerations, we conclude that HAWC's highest energy sources should be interpreted as TeV halos or pulsar wind nebulae, which produce their emission through inverse Compton scattering, and are powered by the rotational kinetic energy of their host pulsar.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Model data for numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

This model-data archive pertains to a study aimed at the numerical evaluation of the photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively). Separating effects of SSZ and HSZ are important for accurately representing stream function as HSZs and SSZs expose solutes to significantly different biogeochemical conditions. We perform numerical experiments using a multiscale reactive transport model for stream corridors implemented in ATS code, which allows for representing multiple storage zones with their respective travel time distributions and biogeochemistry. For each of the numerical experiment, we provide in this model-data archive python wrapper to drive ATS (forward_model.py), synthetic observation (synthetic_btc.py, BTC_observed.csv), forward model files including ATS input (multiscale_transport.tpl), PFLOTRAN inputs (reactions_channel.tpl, reactions_hz.tpl, reactions_sz.tpl), MCMC (mcmc_run.py), predictive uncertainty (pred_uncert.py, BTCs_simulated_day.csv, BTCs_simulated_night.csv), MCMC outputs (tracer_test_GR.npy, tracer_test_logps.npy, tracer_test_parameters.npy) and Jupyter notebook for post-processing and visualization (post-processing.ipynb). For the denitrification application, ATS input file (denitrification_multisubgrid.tpl), PFLOTRAN input files (denitrification_channel.in, denitrification_hz.in), predictive uncertainty (pred_uncert_denitrification.py, BTCs_simulated_DO.csv, BTCs_simulated_DOC.csv, BTCs_simulated_Nitrate.csv).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Highly tunable band inversion in AB 2 X 4 (A=Ge, Sn, Pb; B=As, Sb, Bi; X=Se, Te) compounds

Topological materials have been discovered so far largely by searching for existing compounds in crystallographic databases, but there are potentially new topological materials with desirable features that have not been synthesized. One of the desirable features is high tunability resulting from the band inversion with a very small direct band gap, which can be tuned by changes in pressure or strain to induce a topological phase transition. Here, using density-functional theory (DFT) calculations, we have studied the septuple layered AB 2 X 4 series compounds, where A=(Ge, Sn and Pb), B=(As, Sb and Bi), and X=(Se and Te). With the DFT thermodynamic stability validated by the already-reported compounds in these series, we predict stable Se compounds, which are not found in crystallographic database. Among them, we find that GeBi 2 Se 4 and GeSb 2 Se 4 having a small direct band gap at the Z point are very close to a strong topological insulator, which can be tuned by a moderate pressure to induce the band inversion. Importantly, the topological features with the small direct band gap are well isolated in both momentum and energy windows, which offers high tunability for studying the topological phase transition.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An initial exploration of Bayesian model calibration for estimating the composition of rocks and soils on Mars

The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Visible-Light-Driven Triplet Sensitization of Polycyclic Aromatic Hydrocarbons Using Thionated Perinones

Metal-free chromophores that efficiently generate triplet excited states represent enticing alternatives with respect to transition metal-containing photosensitizers, such as those featuring metal-to-ligand charge transfer excited states. Yet, such molecular constructs have remained underexplored due to the unclear relationship(s) between molecular structure and efficient/rapid intersystem crossing. In this regard, we present a series of three thionated perinone chromophores serving as a newly conceived class of heavy metal-free triplet photosensitizers. We demonstrate that thionation of the lone C=O substituent in each highly fluorescent perinone imparts red-shifted absorbance bands that maintain intense extinction coefficients across the visible spectrum, as well as unusually efficient triplet excited state formation as inferred from the measured singlet O 2 quantum yields at 1270 nm (Φ Δ = 0.78–1.0). Electronic structure calculations revealed the emergence of a low energy S 1 (n → π*) excited state in the proximity of a slightly higher energy S 2 (π → π*) excited state. The distinct character in each of the two lowest-lying singlet state manifolds resulted in the energetic inversion of the corresponding triplet excited states due to differences in electron exchange interactions. Rapid S 1 → T 1 intersystem crossing was thereby facilitated in this manner through spin–orbit coupling as predicted by the El Sayed rules. The lifetimes of the resultant triplet excited states persisted into the microsecond time regime, as measured by transient absorbance spectroscopy, enabling effective bimolecular triplet sensitization of some common polycyclic aromatic hydrocarbons. The synthetically facile interchange of a single O atom to an S atom in the investigated perinones resulted in marked changes to their photophysical properties, namely, conversion of dominant singlet state fluorescence in the former to long-lived triplet excited states in the latter. The combined results suggest a general strategy for accessing long-lived triplet excited states in organic chromophores featuring a lone C=O moiety residing within its structure, valuable for the design of metal-free triplet photosensitizers.

14 SOLAR ENERGY↗
DOE OSTI2020

Majorana-based quantum computing in nanowire devices

The boundary of one-dimensional topological superconductors might lead to the appearance of Majorana zero modes, whose nontrivial exchange statistics can be used for quantum computing. In branched nanowire networks, one can exchange Majorana modes by time dependently tuning topologically nontrivial parameter regions. In this paper, we simulate the exchange of four Majorana modes in T-shaped junctions made out of $\textit{p}$-wave superconducting Rashba wires. Specifically, we derive concrete experimental predictions for (quasi-)adiabatic braiding times and determine geometric conditions for successful Majorana exchange processes. Moreover, we prove that, in the adiabatic limit, the gating time needs to be smaller than the inverse of the squared superconducting order parameter and scales linearly with the gating potential. Furthermore, we show how to circumvent the formation of additional Majorana modes in branched nanowire systems, arising at wire intersection points of narrow junctions. Finally, we propose a multiqubit setup, allowing for universal quantum computing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1971

Pollution detection by digital correlation of multispectral, stero-image pairs.

Remote detection of air pollution circulation patterns is proposed to eventually predict the accumulation of hazardous surface concentrations in time for preventive emission control operations. Earth observations from space platforms will contain information on the height, mean velocity and lateral mixing scales of inversion layers and pollution plumes. Although this information is often not visible on photographs, it could conceivably be retrieved through a digital cross-correlation of multispectral stereo image pairs. Laboratory and field test results are used to illustrate the detection of non-visual inversion layers, the reduction of dominant signal interference, and the spectroscopic identification of combustion products.

Krause, F. R.↗
NASA NTRS1988

The fluorescence of cometary OH

A new solar spectral atlas, improved oscillator strengths, and additional transitions have been incorporated in calculations of the fluorescent equilibrium of the OH molecule in comets as a function of heliocentric distance and radial velocity. Comparisons of detailed, ultraviolet band profiles have been made with the observed, high-resolution spectra of eight comets. While a purely fluorescent equilibrium calculation is sufficient to reproduce the major variations seen in the band structure from one comet to another, some collisional quenching of the inversion of the X2II3/2 Lambda-doublet ground state is required to match the detailed features. The relative population inversion of the upper and lower levels of the ground state is also calculated as a function of both distance and velocity, with the results being consistent with 18 cm radio observations. Predictions of the infrared OH vibrational bands are made; the 1-0 band should be observable if a sufficiently high spectral resolution is employed.

Schleicher, David G.↗