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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological Implications from two Decades of Spectroscopic Surveys at the Apache Point observatory

We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Ly$\alpha$ forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, $r_d$, from eight different samples and six measurements of the growth rate parameter, $f\sigma_8$, from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, $\Lambda$CDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization the BAO data provide nearly an order of magnitude improvement on curvature constraints. The RSD measurements indicate a growth rate that is consistent with predictions from Planck primary data and with General Relativity. When combining the results of SDSS BAO and RSD with external data, all multiple-parameter extensions remain consistent with a $\Lambda$CDM model. Regardless of cosmological model, the precision on $\Omega_\Lambda$, $H_0$, and $\sigma_8$, remains at roughly 1\%, showing changes of less than 0.6\% in the central values between models. The inverse distance ladder measurement under a o$w_0w_a$CDM yields $H_0= 68.20 \pm 0.81 \, \rm km\, s^{-1} Mpc^{-1}$, remaining in tension with several direct determination methods. (abridged)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Design Power Control Strategies of Grid-Forming Inverters for Microgrid Application: Preprint

This paper compares two control schemes in the application control layer of a non-PLL grid-forming (GFM) inverter to gain the insights and understanding of how the two schemes affect the dynamic responses of the GFM inverters and the transition operation of microgrids. Both schemes can achieve smooth microgrid transition operation and power tracking, however, the first scheme expects to exhibit more transients due to the feedforward of the inverse dynamic model of the output filter. Simulation results show that both control schemes achieve successful transition operation with smooth transients, and the second control scheme exhibits slightly better transients in the system frequency during islanding operation and the first one has slightly superior transients in the system frequency during reconnection operation. The power tracking performance in the grid-connected mode is evaluated, which shows that the first scheme has better active power tracking and the second one has better reactive power tracking. The analysis and results are useful to develop reliable control schemes for non-PLL GFM inverters because more and more inverters will work as non-PLL grid-forming sources in the future grids due to better stability and reliability.

grid-forming inverter↗
NASA NTRS2011

Nonlinear Dynamic Inversion Baseline Control Law: Architecture and Performance Predictions

A model reference dynamic inversion control law has been developed to provide a baseline control law for research into adaptive elements and other advanced flight control law components. This controller has been implemented and tested in a hardware-in-the-loop simulation; the simulation results show excellent handling qualities throughout the limited flight envelope. A simple angular momentum formulation was chosen because it can be included in the stability proofs for many basic adaptive theories, such as model reference adaptive control. Many design choices and implementation details reflect the requirements placed on the system by the nonlinear flight environment and the desire to keep the system as basic as possible to simplify the addition of the adaptive elements. Those design choices are explained, along with their predicted impact on the handling qualities.

Miller, Christopher J.↗
DOE OSTI2024

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗
NASA NTRS1996

SSME Condition Monitoring Using Neural Networks and Plume Spectral Signatures

For a variety of reasons, condition monitoring of the Space Shuttle Main Engine (SSME) has become an important concern for both ground tests and in-flight operation. The complexities of the SSME suggest that active, real-time condition monitoring should be performed to avoid large-scale or catastrophic failure of the engine. In 1986, the SSME became the subject of a plume emission spectroscopy project at NASA's Marshall Space Flight Center (MSFC). Since then, plume emission spectroscopy has recorded many nominal tests and the qualitative spectral features of the SSME plume are now well established. Significant discoveries made with both wide-band and narrow-band plume emission spectroscopy systems led MSFC to develop the Optical Plume Anomaly Detection (OPAD) system. The OPAD system is designed to provide condition monitoring of the SSME during ground-level testing. The operational health of the engine is achieved through the acquisition of spectrally resolved plume emissions and the subsequent identification of abnormal emission levels in the plume indicative of engine erosion or component failure. Eventually, OPAD, or a derivative of the technology, could find its way on to an actual space vehicle and provide in-flight engine condition monitoring. This technology step, however, will require miniaturized hardware capable of processing plume spectral data in real-time. An objective of OPAD condition monitoring is to determine how much of an element is present in the SSME plume. The basic premise is that by knowing the element and its concentration, this could be related back to the health of components within the engine. For example, an abnormal amount of silver in the plume might signify increased wear or deterioration of a particular bearing in the engine. Once an anomaly is identified, the engine could be shut down before catastrophic failure occurs. Currently, element concentrations in the plume are determined iteratively with the help of a non-linear computer code called SPECTRA, developed at the USAF Arnold Engineering Development Center. Ostensibly, the code produces intensity versus wavelength plots (i.e., spectra) when inputs such as element concentrations, reaction temperature, and reaction pressure are provided. However, in order to provide a higher-level analysis, element concentration is not specified explicitly as an input. Instead, two quantum variables, number density and broadening parameter, are used. Past experience with OPAD data analysis has revealed that the region of primary interest in any SSME plume spectrum lies in the wavelength band of 3300 A to 4330 A. Experience has also revealed that some elements, such as iron, cobalt and nickel, cause multiple peaks over the chosen wavelength range whereas other elements (magnesium, for example) have a few, relatively isolated peaks in the chosen wavelength range. Iteration with SPECTRA as a part of OPAD data analysis is an incredibly labor intensive task and not one to be performed by hand. What is really needed is the "inverse" of the computer code but the mathematical model for the inverse mapping is tenuous at best. However, building generalized models based upon known input/output mappings while ignoring details of the governing physical model is possible using neural networks. Thus the objective of the research project described herein was to quickly and accurately predict combustion temperature and element concentrations (i.e., number density and broadening parameter) from a given spectrum using a neural network. In other words, a neural network had to be developed that would provide a generalized "inverse" of the computer code SPECTRA.

Hopkins, Randall↗
DOE OSTI2023

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗
NASA NTRS2018

Three-Dimensional Sensor Forward Modeling of Clouds and Precipitation in the Multi-Instrument Inverse Solver Testbed (MIIST)

Sensor forward models are an essential tool for interpreting remote sensing observations and performing quantitative estimates of geophysical parameters. Our three-dimensional forward modeling and retrieval framework allows us to perform detailed analyses of NASA field campaign datasets for a deeper understanding of the remote sensing of clouds and precipitation. This presentation details the componenets of this radiative transfer model used to simulate active (radar) and passive (microwave radiometer) observations, and we give some relevant examples based on both model precipitation systems and actual observations.

Adams, Ian S.↗
NASA NTRS2008

Impact of Surface Emissions to the Zonal Variability of Tropical Tropospheric Ozone and Carbon Monoxide for November 2004

The chemical and dynamical processes governing the zonal variability of tropical tropospheric ozone and carbon monoxide are investigated for November 2004 using satellite observations, in-situ measurements, and chemical transport models in conjunction with inverse-estimated surface emissions. Vertical ozone profile estimates from the Tropospheric Emission Spectrometer (TES) and ozone sonde measurements from the Southern Hemisphere Additional Ozonesondes (SHADOZ) network show the so called zonal 'wave-one' pattern, which is characterized by peak ozone concentrations (70-80 ppb) centered over the Atlantic, as well as elevated concentrations of ozone over Indonesia and Australia (60-70 ppb) in the lower troposphere. Observational evidence from TES CO vertical profiles and Ozone Monitoring Instrument (OMI) NO2 columns point to regional surface emissions as an important contributor to the elevated ozone over Indonesia. This contribution is investigated with the GEOS-Chem chemistry and transport model using surface emission estimates derived from an optimal inverse model, which was constrained by TES and Measurements Of Pollution In The Troposphere (MOPITT) CO profiles (Jones et al., 2007). These a posteriori estimates, which were over a factor of 2 greater than climatological emissions, reduced differences between GEOS-Chem and TES ozone observations by 30-40% and led to changes in GEOS-Chem upper tropospheric ozone of up to 40% over Indonesia. The remaining residual differences can be explained in part by upper tropospheric ozone produced from lightning NOx in the South Atlantic. Furthermore, model simulations from GEOS-Chem indicate that ozone over Indonesian/Australian is more sensitive to changes in surface emissions of NOx than ozone over the tropical Atlantic.

tropical tropospheric ozone↗
DOE OSTI2025

Air temperature and precipitation constraining the modelled wetland methane emissions in a boreal region in northern Europe

Wetland methane responses to temperature and precipitation are studied in a boreal wetland-rich region in northern Europe using ecosystem process models. Six ecosystem models (JSBACH-HIMMELI, LPX-Bern, LPJ-GUESS, JULES, CLM4.5, and CLM5) are compared to multi-model means of ecosystem models and atmospheric inversions from the Global Carbon Project and upscaled eddy covariance flux results for their temperature and precipitation responses and seasonal cycles of the regional fluxes. Two models with contrasting response patterns, LPX-Bern and JSBACH-HIMMELI, are used as priors in atmospheric inversions with Carbon Tracker Europe–CH4 (CTE-CH4) in order to find out how the assimilation of atmospheric concentration data changes the flux estimates and how this alters the interpretation of the flux responses to temperature and precipitation. Inversion moves wetland emissions of both models towards co-limitation by temperature and precipitation. Between 2000 and 2018, periods of high temperature and/or high precipitation often resulted in increased emissions. However, the dry summer of 2018 did not result in increased emissions despite the high temperatures. The process models show strong temperature and strong precipitation responses for the region (51 %–91 % of the variance explained by both). The month with the highest emissions varies from May to September among the models. However, multi-model means, inversions, and upscaled eddy covariance flux observations agree on the month of maximum emissions and are co-limited by temperature and precipitation. The setup of different emission components (peatland emissions, mineral land fluxes) has an important role in building up the response patterns. Considering the significant differences among the models, it is essential to pay more attention to the regional representation of wet and dry mineral soils and periodic flooding which contribute to the seasonality and magnitude of methane fluxes. The realistic representation of temperature dependence of the peat soil fluxes is also important. Furthermore, it is important to use process-based descriptions for both mineral and peat soil fluxes to simulate the flux responses to climate drivers.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1974

Techniques in Doppler gravity inversion

The types of Doppler gravity data available for local as opposed to planetwide geophysical modeling are reviewed. Those gravity fields that are determined dynamically in orbit determination programs yield a smoothed representation of the local gravity field that may be used for quantitative modeling. An estimate of the difference between smoothed and true fields can be considered as a noise limitation in generating local gravity models. A nonlinear inversion for the geometry, depth, and density of the Mare Serenitatis mascon using an ellipsoidal model yielded a global least squares minimum in horizontal dimensions, depth, and thickness-density contrast product. It was subsequently found, by using a linear model, that there were an infinite number of solutions corresponding to various combinations of depth and lateral inhomogeneity. Linear modeling was performed by means of generalized inverse theory.

Phillips, R. J.↗
NASA NTRS1974

An inversion in the atmosphere of Titan

A very detailed greenhouse model derives a methane to hydrogen ratio of unity and a minimum surface pressure of 0.4 atm. Based on a surface gravity g = 140 cm sec/2, the minimum CH4 abundance is 30-40 km-A and the minimum H2 abundance varies from 15 to 85 km-A. A model of the atmosphere of Titan is proposed which seems to be consistent with observations and requires a much smaller CH4 abundance (of the order or 2 km-atm). Although no H2 is required, the presence of some H2 is readily accommodated. In this model, a temperature inversion exists in the atmosphere due to absorption of blue and ultraviolet solar radiation by small particles. The absorbed radiation is re-radiated by the dust and by molecules having long wavelength bands such as CH4 7.7 micrometer and ethane at 12.2 micrometer. The brightness temperature at 20 micrometer is primarily due to re-radiation by the dust.

Danielson, R. E.↗
NASA NTRS2013

Adaptive Sampling of Spatiotemporal Phenomena with Optimization Criteria

This work was designed to find a way to optimally (or near optimally) sample spatiotemporal phenomena based on limited sensing capability, and to create a model that can be run to estimate uncertainties, as well as to estimate covariances. The goal was to maximize (or minimize) some function of the overall uncertainty. The uncertainties and covariances were modeled presuming a parametric distribution, and then the model was used to approximate the overall information gain, and consequently, the objective function from each potential sense. These candidate sensings were then crosschecked against operation costs and feasibility. Consequently, an operations plan was derived that combined both operational constraints/costs and sensing gain. Probabilistic modeling was used to perform an approximate inversion of the model, which enabled calculation of sensing gains, and subsequent combination with operational costs. This incorporation of operations models to assess cost and feasibility for specific classes of vehicles is unique.

Chien, Steve A.↗
NASA NTRS1993

The North Atlantic circulation in the early 1980s - An estimate from inversion of a finite-difference model

A finite-difference model of the North Atlantic is constructed for the purpose of making an estimate of the circulation through an inverse calculation. The data base is eclectic, and includes hydrography, oxygen, nutrients, current meter and float records, atmospheric momentum, heat and water vapor transfers, as well as estimates of certain integral fluxes. Owing to the available hydrographic database, the model resolution is restricted to 1 deg at best, and is much coarser in many aspects. This limited resolution is a major obstacle to accurate estimates of climatological fluxes. In its final form, there are about 9000 constraints in 29,000 formal unknowns plus 9000 noise unknowns. The system is solved as a tapered least-squares system by a sparse conjugate gradient algorithm. With the exception of a few float velocities, all constraints are found to be consistent within error estimates. The model produces estimates of large-scale fluxes and flux divergences for all conventional properties including heat and nutrients as well as carbon dioxide and alkalinity. Meridional fluxes of carbon are found to be indistinguishable from zero, whereas the North Atlantic tends to export nutrients to the south, but carry heat to the north. Traditional oceanographic depictions of the circulation through combination of nonsynoptic data into steady models may have reached their useful limit in the present calculation, as the conflicts between the data and physical requirements become quantitatively apparent.

Martel, Francoise↗
NASA NTRS2018

Vegetation Height Estimate in Rice Fields Using Single Polarization TanDEM-X Science Phase Data

This study presents the retrieval of rice paddy height using single polarization (single-pol) interferometric SAR (InSAR) data by means of model-based height inversion without an external DEM. A total of eight TanDEM-X (TDX) scenes were used and the TDX images were; acquired using a large cross-track baseline configuration during the TDX Science Phase (June–August 2015). A single-pol inversion approach for a flooded rice field is proposed and evaluated over the Buan test site in South Korea. A novel approach is adopted for the estimation of the ground (i.e., water level) interferometric phase within a flooded rice paddy from TDX data acquired during a period of early rice growth. It is consequently possible to apply the model-based inversion algorithm to the single-pol InSAR data. Rice height maps during the rice growth cycle presented and validated by field measurements. The results demonstrated the high performance of the inversion with a correlation coefficient of 0.78 and an RMSE of 0.10 m. The proposed methodology will be useful to monitor rice plants and to predict a gross rice yield, along with dual and fully polarimetric interferometric SAR data.

Lee, Seung-Kuk↗
DOE OSTI2020

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Three-dimensional cooperative inversion of airborne magnetic and gravity gradient data using deep-learning techniques

Using multiple geophysical methods has become a prevailing approach in numerous geophysical applications to investigate subsurface structures and parameters. These multimethod-based exploration strategies have the potential to greatly diminish uncertainties and ambiguities encountered during geophysical data analysis and interpretation. One of the applications is the cooperative inversion of airborne magnetic and gravity gradient data for the interpretation of data obtained in mineral, oil and gas, and geothermal explorations. In this paper, a unified cooperative inversion framework is designed by combining the standard separate inversions with a deep neural network (DNN), which serves as the link between different types of data. A well-trained DNN takes the separately inverted susceptibility and density models as the inputs and provides improved models that will be used as the initial models of deterministic inversions. A two-round iteration strategy is adopted to guarantee the reasonability of the recovered models and overall efficiency of the inversion. In addition, this deep-learning (DL)-based framework demonstrates excellent generalization abilities when tested on models that are entirely distinct from the training data sets. The framework can easily incorporate multiphysics without necessitating any structural changes to the network. Synthetic experiments validate that our DL-based method outperforms conventional separate inversions and cross-gradient-based joint inversion in view of the accuracy of the recovered models and inversion efficiency. Successful application to field data further verifies the effectiveness of our DL-based method.

Geochemistry & Geophysics↗