Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Aircraft automatic flight control system with model inversion

A simulator study was conducted to verify the advantages of a Newton-Raphson model-inversion technique as a design basis for an automatic trajectory control system in an aircraft with highly nonlinear characteristics. The simulation employed a detailed mathematical model of the aerodynamic and propulsion system performance characteristics of a vertical-attitude takeoff and landing tactical aircraft. The results obtained confirm satisfactory control system performance over a large portion of the flight envelope. System response to wind gusts was satisfactory for various plausible combinations of wind magnitude and direction.

Smith, G. A.↗
DOE OSTI2023

Conditional Karhunen–Loève regression model with Basis Adaptation for high-dimensional problems: Uncertainty quantification and inverse modeling

Here, we propose a methodology for improving the accuracy of surrogate models of the observable response of physical systems as a function of the systems’ spatially heterogeneous parameter fields, with applications to uncertainty quantification and parameter estimation in high-dimensional problems. Practitioners often formulate finite-dimensional representations of spatially heterogeneous parameter fields using truncated unconditional Karhunen–Loève expansions (KLEs) for a certain choice of unconditional covariance kernel and construct surrogate models of the observable response with respect to the KLE coefficients. When direct measurements of the parameter fields are available, we propose improving the accuracy of these surrogate models by representing the parameter fields via conditional Karhunen-Loève expansions (CKLEs). CKLEs are constructed by conditioning the covariance kernel of the unconditional expansion on the direct measurements of the parameter field via Gaussian process regression, and then truncating the corresponding KLE. We apply the proposed methodology to constructing surrogate models via the Basis Adaptation (BA) method of the stationary hydraulic head response, measured at spatially discrete observation locations, of a groundwater flow model of the Hanford Site, as a function of the 1000-dimensional representation of the model’s log-transmissivity field. We find that BA surrogate models of the hydraulic head based on CKLEs are more accurate than BA surrogate models based on unconditional expansions for forward uncertainty quantification tasks. Furthermore, we find that inverse estimates of the hydraulic transmissivity field computed using CKLE-based BA surrogate models are more accurate than those computed using unconditional BA surrogate models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Dataset for SIAM MPI23 Project "Model inversion for complex physical systems using low-dimensional surrogates"

Training and testing data for the SIAM MPI23 workshop "Model inversion for complex physical systems using low-dimensional surrogates". The data consists of ensembles of input and output pairs corresponding to queries of a 2D saturated groundwater flow model. The inputs consist of random vectors which map to discretized model parameter fields via a Kosambi-Karhunen-Loève expansion. Output corresponds to discretized pressure fields.

basis adaptation, conditional Karhunen-Loève expan↗
NASA NTRS1989

Land-surface temperature measurement from space - Physical principles and inverse modeling

To apply the multiple-wavelength (split-window) method used for satellite measurement of sea-surface temperature from thermal-infrared data to land-surface temperatures, the authors statistically analyze simulations using an atmospheric radiative transfer model. The range of atmospheric conditions and surface temperatures simulated is wide enough to cover variations in clear atmospheric properties and surface temperatures, both of which are larger over land than over sea. Surface elevation is also included in the simulation as the most important topographic effect. Land covers characterized by measured or modeled spectral emissivities include snow, clay, sands, and tree leaf samples. The empirical inverse model can estimate the surface temperature with a standard deviation less than 0.3 K and a maximum error less than 1 K, for viewing angles up to 40 degrees from nadir under cloud-free conditions, given satellite measurements in three infrared channels. A band in the region from 10.2 to 11.0 microns will usually give the most reliable single-band estimate of surface temperature. In addition, a band in either the 3.5-4.0-micron region or in the 11.5-12.6-micron region must be included for accurate atmospheric correction, and a band below the ozone absorption feature at 9.6 microns (e.g., 8.2-8.8 microns) will increase the accuracy of the estimate of surface temperature.

Wan, Zhengming↗
NASA NTRSDate not supplied

Quantifying the benefits of improved satellite remote-sensing observations for inverse modeling of NOx and NMVOC emissions

This study aims to demonstrate the benefits of using novel high spatiotemporal retrieval products from newer satellites for top-down emission estimates of nitrogen oxides (NO x ) and non-methane volatile organic compounds (NMVOCs) for the summer of 2019 over the contiguous United States. Recent satellite retrievals have not only advanced spatiotemporal resolution but also greatly reduced error and uncertainty due to reduced noise in the retrievals compared to spaceborne sensors launched in the past. We applied inverse modeling techniques using tropospheric nitrogen dioxide (NO 2 ) and formaldehyde (H-CHO) column retrieval products from the Ozone Monitoring Instrument (OMI) and TROPOspheric Monitoring Instrument (TROPOMI) in conjunction with the Weather Research Forecast and Community Multiscale Air Quality Modeling system (WRF-CMAQ). In order to provide a better representation of background chemical composition and avoid misalignment of emission adjustment, we applied monthly scaling factors for ozone (O 3 ) and CO boundary concentrations in addition to the inclusion of lightning and aviation emissions. Satellite-constrained NO x and NMVOCs posterior emissions showed a mitigated discrepancy between observed and modeled columns. The improvement in the model performance was greater when using TROPOMI, primarily benefiting from reduced errors/biases of the satellite retrievals that enabled us to explore corresponding changes in O 3 concentrations and production sensitivity regimes using the ratio of H-CHO and NO 2 .

remote sensing↗
DOE OSTI2022

Intra-Specific Variability in Plant Hydraulic Parameters Inferred From Model Inversion of Sap Flux Data

Understanding plant hydraulic regulation is critical for predicting plant and ecosystem responses to projected increases in drought stress. Plant hydraulic regulation is controlled by observable, diverse plant hydraulic traits that can vary as much across individuals of the same species as they do across different species. Direct measurements of plant hydraulic traits from a range of ecosystems remain limited in comparison to other, more readily measured traits (e.g., specific leaf area). Furthermore, plant hydraulic trait measurements, often made at leaf or branch levels, are not easily scaled to whole-plant values that are typically used to predict plant and ecosystem fluxes. In this study, multiple whole-plant hydraulic parameters are inferred from observations of plant water use (i.e., sap flow), soil properties, and meteorological data. We use a Markov Chain Monte Carlo model inversion approach to obtain the best estimates and uncertainty of plant hydraulic parameters that capture whole-plant effective embolism resistance and stomatal sensitivity to decreasing plant water potential. We then use the inferred values in the model to estimate whole-tree water use and isohydricity. This approach reliably infers whole-plant parameter values with enough specificity to resolve inter- and intra-specific differences, and thus supplements time- and labor-intensive direct measurements of traits.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2013

Generalized Ocean Color Inversion Model for Retrieving Marine Inherent Optical Properties

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future ensemble applications.

Werdell, P. Jeremy↗
NASA NTRS2011

SO2 Emissions and Lifetimes: Estimates from Inverse Modeling Using In Situ and Global, Space-Based (SCIAMACHY and OMI) Observations

Top-down constraints on global sulfur dioxide (SO2) emissions are inferred through inverse modeling using SO2 column observations from two satellite instruments (SCIAMACHY and OMI). We first evaluated the S02 column observations with surface SO2 measurements by applying local scaling factors from a global chemical transport model (GEOS-Chem) to SO2 columns retrieved from the satellite instruments. The resulting annual mean surface SO2 mixing ratios for 2006 exhibit a significant spatial correlation (r=0.86, slope=0.91 for SCIAMACHY and r=0.80, slope = 0.79 for OMI) with coincident in situ measurements from monitoring networks throughout the United States and Canada. We evaluate the GEOS-Chem simulation of the SO2 lifetime with that inferred from in situ measurements to verity the applicability of GEOS-Chem for inversion of SO2 columns to emissions. The seasonal mean SO2 lifetime calculated with the GEOS-Chem model over the eastern United States is 13 h in summer and 48 h in winter, compared to lifetimes inferred from in situ measurements of 19 +/- 7 h in summer and 58 +/- 20 h in winter. We apply SO2 columns from SCIAMACHY and OMI to derive a top-down anthropogenic SO2 emission inventory over land by using the local GEOS-Chem relationship between SO2 columns and emissions. There is little seasonal variation in the top-down emissions (<15%) over most major industrial regions providing some confidence in the method. Our global estimate for annual land surface anthropogenic SO2 emissions (52.4 Tg S/yr from SCIAMACHY and 49.9 Tg S / yr from OMI) closely agrees with the bottom-up emissions (54.6 Tg S/yr) in the GEOS-Chem model and exhibits consistency in global distributions with the bottom-up emissions (r = 0.78 for SCIAMACHY, and r = 0.77 for OMI). However, there are significant regional differences.

Lee, Chulkyu↗
NASA NTRS2019

Developing a soil inversion model framework for regional permafrost monitoring

Currently, the community lacks capabilities to assess and monitor landscape scale permafrost active layer dynamics over large extents. To address this need, we developed a concept of a remote sensing based Soil Inversion Model for regional Permafrost (SIM-P) monitoring. The current SIM-P framework includes a satellite-based soil process model and a soil dielectric model. We are also working on incorporating a radar scattering model for Arctic tundra into the SIM-P framework. A unified soil parameterization scheme was developed to harmonize key soil thermal, hydraulic and dielectric parameters in the soil process and radar models that can be used in the joint soil-radar inversion framework. The soil parameter retrievals of the SIM-P framework include soil organic content (SOC) and active layer thickness (ALT). Initial tests of SIM-P using in-situ soil permittivity observations showed reasonable accuracy in predicting site-level SOC and soil temperature profiles at an Alaska tundra site and ALT in Arctic Alaska. SIM-P will be further tested using airborne P- and L-band radar data collected during NASA’s Arctic Boreal Vulnerability Experiment (ABoVE) to evaluate the sensitivity of longwave radar to active layer properties.

Miller, Charles E.↗
DOE OSTI2026

Data and Code for: Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits

This repository contains the simulation outputs and processing scripts associated with the study of winter wheat traits across the United States, utilizing the Ecosys agroecosystem model. The dataset includes model results for both rainfed and irrigated winter wheat systems, supporting the findings presented in the manuscript titled "Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits." Data includes the original Ecosys simulation outputs (archived in .db format within the compressed .zip files) and extracted analysis data (stored in .pkl files for efficient processing). Python code for data processing and figure generation is provided in a Jupyter notebook. External Observational Datasets should refer to the following official repositories for the input and validation data used in this study. The eddy covariance data from the AmeriFlux network (https://ameriflux.lbl.gov/). Climate-forcing data of NLDAS-2 from NASA LDAS (https://ldas.gsfc.nasa.gov/nldas/nldas-2-forcing-data). Soil data from the Gridded Soil Survey Geographic Database (gSSURGO), available at (https://www.nrcs.usda.gov/resources/data-and-reports/gridded-soil-survey-geographic-gssurgo-database). Crop yields, planting and harvest dates from the USDA public databases (https://quickstats.nass.usda.gov/; https://webapp.rma.usda.gov/apps/actuarialinformationbrowser/CropCriteria.aspx). Satellite-derived SLOPE GPP data from ORNL DAAC (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1786). Land use and crop progress information from the USDA Crop Data Layer and Crop Progress and Condition Gridded Layers (https://www.nass.usda.gov/Research_and_Science/). The Ecosys model code is available online at https://github.com/jinyun1tang/ECOSYS.

Wheat↗
NASA NTRS2024

Inverse Modeling of the Initial Stage of the 1991 Pinatubo Volcanic Cloud Accounting for Radiative Feedback of Volcanic Ash

The way volcanic clouds evolve is very sensitive to the initial spatial 3D distributions of volcanic materials, which are often unknown. In this study, we conducted inverse modeling of the Mt. Pinatubo cloud using TOMS 2D mapping of Aerosol Index and SO 2 loading during the first three post-eruption days to estimate the time-dependent emissions profiles and initial 3D spatial distributions of volcanic ash and SO 2 . We account for aerosol radiative feedback and dynamic lofting of volcanic ash in the inversion calculations for the first time. This resulted in a lower ash injection height (by 1.5 km for ash) than without ash radiative feedback. The Pinatubo eruption ejected ≈77% of fine ash at 12 to 23 km, ≈65% of SO 2 at 18–25 km. In contrast with previous studies, which suggested that all volcanic materials were emitted above the tropopause, a significant fraction of SO 2 (5.1 of 15.5 Mt) and fine ash (37.2 of 66.5 Mt) were ejected in the troposphere, where SO 2 quickly oxidized into sulfate aerosol that is short-lived in the troposphere. This explains the early presence of sulfate aerosols in the plume and why the models can reproduce the observed volcanic aerosols' optical depth (AOD), assuming lower-than-observed SO 2 emission in the stratosphere. Despite the quicker than in observations build-up of sulfate AOD, in a month after the eruption, the evolution of the Pinatubo AOD simulated using the obtained ash and SO 2 initial distributions converges with the available stratospheric aerosol and gas experiment (SAGE) observations.

Pinatubo volcanic ash and SO2 emissions↗
NASA NTRS2011

An Inverse Modeling Approach to Estimating Phytoplankton Pigment Concentrations from Phytoplankton Absorption Spectra

Phytoplankton absorption spectra and High-Performance Liquid Chromatography (HPLC) pigment observations from the Eastern U.S. and global observations from NASA's SeaBASS archive are used in a linear inverse calculation to extract pigment-specific absorption spectra. Using these pigment-specific absorption spectra to reconstruct the phytoplankton absorption spectra results in high correlations at all visible wavelengths (r(sup 2) from 0.83 to 0.98), and linear regressions (slopes ranging from 0.8 to 1.1). Higher correlations (r(sup 2) from 0.75 to 1.00) are obtained in the visible portion of the spectra when the total phytoplankton absorption spectra are unpackaged by multiplying the entire spectra by a factor that sets the total absorption at 675 nm to that expected from absorption spectra reconstruction using measured pigment concentrations and laboratory-derived pigment-specific absorption spectra. The derived pigment-specific absorption spectra were further used with the total phytoplankton absorption spectra in a second linear inverse calculation to estimate the various phytoplankton HPLC pigments. A comparison between the estimated and measured pigment concentrations for the 18 pigment fields showed good correlations (r(sup 2) greater than 0.5) for 7 pigments and very good correlations (r(sup 2) greater than 0.7) for chlorophyll a and fucoxanthin. Higher correlations result when the analysis is carried out at more local geographic scales. The ability to estimate phytoplankton pigments using pigment-specific absorption spectra is critical for using hyperspectral inverse models to retrieve phytoplankton pigment concentrations and other Inherent Optical Properties (IOPs) from passive remote sensing observations.

Moisan, John R.↗
DOE OSTI2025

Inverse model based error detection in beamline optics

Optics tuning in transfer lines and LINACs can be challenging due to the fact that multiple combinations of machine settings can lead to the same diagnostic output. Moreover, the lack of a periodic solution can limit the ability to infer optics in the same way as rings from BPM signals. Model based approaches are often used to assist with the optics tuning in combination with optimization or parameter estimation. Here we have developed a novel approach using machine learning inverse models trained on a known configuration to detect variations in quadrupole settings without explicitly including them in the model. This paper shows a comparison of neural network models and linear models on both a simulation based study and experimental studies conducted at the AGS to RHIC transfer line at Brookhaven National Lab.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Large Ensemble Particle Filter for Spatial Climate Reconstructions Using a Linear Inverse Model

Proxy records that document the last 2000 years of climate provide evidence for the wide range of the natural climate variability from inter-annual to secular timescales not captured by the short window of recent direct observations. Assessing climate models ability to reproduce such natural variations is crucial to understand climate sensitivity and impacts of future climate change. Paleoclimate data assimilation (PDA) offers a powerful way to extend the short instrumental period by optimally combining the physics described by General Circulation Climate Models (GCMs) with information from available proxy records while taking into account their uncertainties. Here we present a new PDA approach based on a sequential importance resampling (SIR) Particle filter (PF) that uses Linear Inverse Modeling (LIM) as an emulator of several CMIP-class GCMs. We examine in a perfect-model framework the skill of the various LIMs to forecast the dynamics of the surface temperatures and provide spatial field reconstructions over the last millennium in a SIR PF. Our results show that the LIMs allow for skillful ensemble forecasts at 1-year lead-time based on GCMs dynamical knowledge with best prediction in the tropics and the North Atlantic. The PDA further provides a set of physically consistent spatial fields allowing robust uncertainty quantification related to climate models biases and proxy spatial sampling. Our results indicate that the LIM yields dynamical memory improving climate variability reconstructions and support the use of the LIM as a GCM-emulator in real reconstruction to propagate large ensembles of particles at low cost in SIR PF.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

Inverse Modelling to Obtain Head Movement Controller Signal

Experimentally obtained dynamics of time-optimal, horizontal head rotations have previously been simulated by a sixth order, nonlinear model driven by rectangular control signals. Electromyography (EMG) recordings have spects which differ in detail from the theoretical rectangular pulsed control signal. Control signals for time-optimal as well as sub-optimal horizontal head rotations were obtained by means of an inverse modelling procedures. With experimentally measured dynamical data serving as the input, this procedure inverts the model to produce the neurological control signals driving muscles and plant. The relationships between these controller signals, and EMG records should contribute to the understanding of the neurological control of movements.

Kim, W. S.↗
DOE OSTI2022

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Forward and Inverse Models for Satellite Remote Sensors using Principal Component Analysis

Satellite remote sensors such as AIRS on Aqua, CrIS on S-NPP, NOAA20 and JPSS-2, IASI on Metop A, B, and C make millions of observations each day with thousands of spectral channels for each observation; this poses challenges for efficiently inversion of the inherently large dataset as needed to retrieve atmospheric and surface properties. This presentation will illustrate the use of Principal Component Analysis (PCA) to speed up radiative transfer forward model calculations and to stabilize the inversion algorithms. A Principal Component-based radiative transfer model (PCRTM) developed at NASA Langley Research Center can simulate top of atmosphere (TOA) radiance or reflectance spectra from 50 cm-1 to 50000 cm-1 (200 m to 0.20 m quickly and accurately. PCRTM demonstrated very high accuracy relative to reference line-by-line radiative transfer models and it saves orders of magnitude computational time. Examples of the PCRTM model developed for hyperspectral sensors such as AIRS, CrIS, IASI, NAST-I, SHIS, CPF, TEMPO, SBG, OMI, and SCIAMACHY will be presented. In addition to using the PCRTM as forward model, the NASA Langley developed inversion algorithm also uses PCA to compress the state vector into a compressed dimension to speed up and stabilize the inversion process. Examples of retrieved atmospheric temperature, water vapor, CO2, CO, CH4, N2O, and O3 profiles, cloud properties (optical depth, size, phase, and height), and surface properties (surface emissivity spectra and skin temperatures) will be presented. This algorithm is being transitioned to the NASA Sounder SIPS and NASA's Goddard Earth Sciences Data and Information Services Center (GES DISC).

forward model↗
NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗