Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverse optimization”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Uncertainty Assessment of Space-Borne Passive Soil Moisture Retrievals

The uncertainty associated with passive soil moisture retrieval is hard to quantify, and known to be underlain by various, diverse, and complex causes. Factors affecting space-borne retrieved soil moisture estimation include: (i) the optimization or inversion method applied to the radiative transfer model (RTM), such as e.g. the Single Channel Algorithm (SCA), or the Land Parameter Retrieval Model (LPRM), (ii) the selection of the observed brightness temperatures (Tbs), e.g. polarization and incidence angle, (iii) the definition of the cost function and the impact of prior information in it, and (iv) the RTM parameterization (e.g. parameterizations officially used by the SMOS L2 and SMAP L2 retrieval products, ECMWF-based SMOS assimilation product, SMAP L4 assimilation product, and perturbations from those configurations). This study aims at disentangling the relative importance of the above-mentioned sources of uncertainty, by carrying out soil moisture retrieval experiments, using SMOS Tb observations in different settings, of which some are mentioned above. The ensemble uncertainties are evaluated at 11 reference CalVal sites, over a time period of more than 5 years. These experimental retrievals were inter-compared, and further confronted with in situ soil moisture measurements and operational SMOS L2 retrievals, using commonly used skill metrics to quantify the temporal uncertainty in the retrievals.

Quets, Jan↗
NASA NTRS2017

Hyperspectral Radiative Transfer Modeling to Explore the Combined Retrieval of Biophysical Parameters and Canopy Fluorescence from FLEX - Sentinel-3 Tandem Mission Multi-Sensor Data

The FLuorescence EXplorer (FLEX) satellite mission, selected as ESA's 8th Earth Explorer, has been designed forthe measurement of sun-induced fluorescence (F) spectra emitted by plants. This will be accomplished through amulti-sensor approach by placing it in a common orbit in tandem with the Sentinel-3 (S3) mission, which willhave two optical sensors on board, OLCI (Ocean and Land Colour Instrument) and SLSTR (Sea and Land SurfaceTemperature Radiometer) to complement FLEX. These S3 instruments will be used in combination with theimaging spectrometers on board FLEX to provide data useful for atmospheric correction of FLEX data. However,a fully synergetic approach, i.e. by exploiting the spectral and directional information from all tandem missioninstruments together, is an attractive alternative which is explored in this paper. By employing all combined topof-atmosphere (TOA) spectral radiance data, one can (i) characterize the relevant optical properties of the atmosphere,(ii) retrieve biophysical canopy properties including the associated reflectance anisotropy, and (iii)retrieve a more accurate and consistent canopy F.Regarding retrieval methods, Fraunhofer Line Depth (FLD) and Spectral Fitting (SF) are well-known techniquesapplied to hyperspectral data. Both methods depend on a high spectral resolution and assume aLambertian (isotropic) canopy reflectance. However, most vegetation canopies are non-Lambertian. This impliesthat, in particular when ignoring the anisotropic surface reflection, substantial retrieval errors can occur due tothe interaction between atmospheric absorption bands and surface reflectance anisotropy. In this paper, a novelmethod based on spectral radiative transfer (RT) modeling is proposed, in which coupled RT models are used tosimulate TOA radiance spectra. These are then matched with ‘measured' spectra in order to retrieve surfacefluorescence, along with a suite of biophysical parameters, by model inversion through optimization. By applyingcoupled RT models of the soil-leaf-canopy and the surface-atmosphere systems, TOA radiance spectra canbe simulated for all optical sensors of this tandem mission. In this way, complex effects due to surface reflectanceanisotropy and the spectral sampling by the various instruments, which are difficult to compensate for in the endproducts, are properly taken into account by their incorporation in the forward modeling. Next, by model inversionof TOA radiance data via optimization, the most accurate F retrievals can be achieved in a consistentmanner, along with important canopy level biophysical parameters that may help interpret the F spectrum, suchas chlorophyll content and leaf area index (LAI). The potential of this approach has been explored in a numericalexperiment, and the results are presented in this paper. We find that, with the assumed well-characterized andplausible FLEX/S3 instrument performances, the simultaneous retrieval of biophysical canopy parameters and Fspectra would be possible with a remarkable accuracy, provided the correct atmospheric characterization isavailable.

Verhoef, Wouter↗
NASA NTRSDate not supplied

Optimizing spectroheliogram deconvolution methods : MaGIXS - A case study

Over the past five years, new methods to reconstruct spectrally pure maps of the Sun from spectroheliogram data have emerged. In particular, with the rebirth of long-abandoned slitless imaging spectroscopy, one may obtain both spatial and spectral information over a large field of view simultaneously. Yet, depending on the size of the extended source combined with the extent of spectral dispersion, there will be locations in the focal plane where spectral lines from different spatial locations overlap and must be deconvolved. An unfolding method has been successfully developed and demonstrated on the recent rocket flight MaGIXS, which observed several strong emission lines (9 to 30$\AA$) from different portions of two active regions. In the work we are going to present, we conduct a systematic investigation of the parameters that controls and optimizes the inversion method to unfold overlappogram data. We also demonstrate a derived method of the inversion that does not require previous assumptions on the thermal and ionization equilibrium and abundance state of the plasma.

Spectroheliogram↗
NASA NTRSDate not supplied

Homogenization of Satellite Based Hyperspectral Infrared Sounder Data To Build Long Term Climate Record

Building long term climate record using data from multiple hyperspectral Infrared (IR)sounders requires the homogenization of different data record to ensure the consistency andcontinuity. Such a requirement comes from two perspectives: 1) the need to adjust theoverlapping measurements of different sounders to ensure the radiometric consistency in thespectral radiance domain; 2) the need for a rigorously defined scheme to ensure the radiometricconsistency being transferred to the essential climate variables derived from the radiance record.We develop a solution that uses a spectral fingerprinting scheme to derive anomalies of keyclimate variables from long term spectral radiance data record constructed using both AIRS andCrIS observations aboard AQUA, SNPP and JPSS satellites. The fingerprinting scheme usescommon radiative kernels for all sounder measurements and therefore effectively avoids thealgorithm introduced inconsistency in retrieved geophysical variables. Our approach uses aunified sampling scheme to match AIRS and CrIS observations in both spectral and spatial-temporal domain, facilitating the intercomparison of spectral radiances from different sensors(platforms). The optimized liner inversion scheme allows the direct quantification and thereforethe adjustment for the impact on the derived climate anomalies imposed by potential radiometricinconsistency between the overlapping measurements. Such a scheme also enables the low-latency data processing of long term hyperspectral sounder data records. This paper provides a detailed introduction of the spectral fingerprinting methodology.Also introduced here is the climate fingerprinting Sounder Product (ClimFiSP) developed basedon the fingerprinting methodology. ClimFiSP products include the space-time averagedproperties of key climate variables that are derived from the long-term, space-time averagedradiances from AQUA-AIRS, SNPP-CrIS, and JPSS1-CrIS. ClimFiSP will be available to usersthrough NASA'sGoddard Earth Sciences Data and Information Services Center (GES DISC).

Wan Wu↗
NASA NTRS2019

Is There an Optimal CO2 Partial Column for Flux Inversions?

The fidelity of flux estimates from an atmospheric inversion depends on the ability of atmospheric transport models to simulate the measured quantity. For species such as CO2, with surface fluxes and a large network of surface measurements, this means correctly simulating the dynamics of the planetary boundary layer (PBL), which is one of the most uncertain aspects of atmospheric transport modeling. In contrast, the simulated total column average mole fraction of CO2 (XCO2) is largely insensitive to simulated PBL dynamics. Therefore, measurements of XCO2 provided by current and future short wave infrared (SWIR) greenhouse gas (GHG) satellites such as GOSAT and the OCO family would seem to be more appropriate to flux inversions, as far as minimizing transport model errors (the "noise") is concerned. Unfortunately, the flux-induced variation of CO2 (the "signal") is the largest within the PBL and smallest in XCO2. Therefore, assimilating XCO2 as opposed to PBL CO2 need not give us the strongest "signal to noise" in flux inversions. Recent work on GOSAT and OCO2 retrievals suggest that SWIR satellite spectra may be used to estimate a lower partial column CO2, which could be assimilated in a flux inversion, instead of XCO2.Here we report on a study to assess whether there is an optimal partial column average CO2, intermediate between PBL CO2 and XCO2, whose assimilation might yield the best signal to noise in flux inversions, where (as before) "signal" is the flux-induced variation and "noise" is the error in transport modeling. We simulate atmospheric CO2 with five different global transport models and a common surface CO2 flux over ten years. We consider the spread across the five models to be a proxy for transport model error (the "noise"), and the common variation of CO2 in all five models to be a proxy for the "signal". We compare these signals and noises at different spatiotemporal scales for different partial column specifications to investigate whether there exists an optimal partial column that has large surface flux-driven variations and yet is relatively insensitive to errors in transport models. Finally, we comment on the feasibility of estimating such a partial column from current and future SWIR GHG satellites in the light of recent work on vertically resolved CO2 from current SWIR GHG satellites.

Basu, Sourish↗
NASA NTRS2003

Modeling the Absorbing Aerosol Index

We propose a scheme to model the absorbing aerosol index and improve the biomass carbon inventories by optimizing the difference between TOMS aerosol index (AI) and modeled AI with an inverse model. Two absorbing aerosol types are considered, including biomass carbon and mineral dust. A priori biomass carbon source was generated by Liousse et al [1996]. Mineral dust emission is parameterized according to surface wind and soil moisture using the method developed by Ginoux [2000]. In this initial study, the coupled CCM1 and GRANTOUR model was used to determine the aerosol spatial and temporal distribution. With modeled aerosol concentrations and optical properties, we calculate the radiance at the top of the atmosphere at 340 nm and 380 nm with a radiative transfer model. The contrast of radiance at these two wavelengths will be used to calculate AI. Then we compare the modeled AI with TOMS AI. This paper reports our initial modeling for AI and its comparison with TOMS Nimbus 7 AI. For our follow-on project we will model the global AI with aerosol spatial and temporal distribution recomputed from the IMPACT model and DAO GEOS-1 meteorology fields. Then we will build an inverse model, which applies a Bayesian inverse technique to optimize the agreement of between model and observational data. The inverse model will tune the biomass burning source strength to reduce the difference between modelled AI and TOMS AI. Further simulations with a posteriori biomass carbon sources from the inverse model will be carried out. Results will be compared to available observations such as surface concentration and aerosol optical depth.

Penner, Joyce↗
NASA NTRS1992

Trajectory optimization for the National Aerospace Plane

The primary objective of this research is to develop an efficient and robust trajectory optimization tool for the optimal ascent problem of the National Aerospace Plane (NASP). This report is organized in the following order to summarize the complete work: Section two states the formulation and models of the trajectory optimization problem. An inverse dynamics approach to the problem is introduced in Section three. Optimal trajectories corresponding to various conditions and performance parameters are presented in Section four. A midcourse nonlinear feedback controller is developed in Section five. Section six demonstrates the performance of the inverse dynamics approach and midcourse controller during disturbances. Section seven discusses rocket assisted ascent which may be beneficial when orbital altitude is high. Finally, Section eight recommends areas of future research.

Lu, Ping↗
NASA NTRS1996

A PDE Sensitivity Equation Method for Optimal Aerodynamic Design

The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Borggaard, Jeff↗
NASA NTRS1987

Optimal control of large space structures via generalized inverse matrix

Independent Modal Space Control (IMSC) is a control scheme that decouples the space structure into n independent second-order subsystems according to n controlled modes and controls each mode independently. It is well-known that the IMSC eliminates control and observation spillover caused when the conventional coupled modal control scheme is employed. The independent control of each mode requires that the number of actuators be equal to the number of modelled modes, which is very high for a faithful modeling of large space structures. A control scheme is proposed that allows one to use a reduced number of actuators to control all modeled modes suboptimally. In particular, the method of generalized inverse matrices is employed to implement the actuators such that the eigenvalues of the closed-loop system are as closed as possible to those specified by the optimal IMSC. Computer simulation of the proposed control scheme on a simply supported beam is given.

Nguyen, Charles C.↗
NASA NTRSDate not supplied

Asteroseismic Inversions of Mixed Acoustic-Gravity Modes to Probe the Stellar Core Structure

The discovery of mixed acoustic-gravity modes of oscillations of moderate mass stars opens a unique opportunity to infer the structure of the inner energy-generating cores and thus test the stellar evolution theory. The mixed modes have properties of internal gravity waves (g-modes) in the convectively stable helium core and properties of acoustic modes outside the core. We select several sets of the oscillation mode frequencies in the mass range from about 1.3 to 1.6 solar masses from the Kepler Legacy database, and apply the optimally localized averaging inversion technique previously developed for low-degree helioseismology. The inversion technique takes into account the uncertainties in the determination of the mass and radius of the stars, as well as the surface effects. The methodology provides sensitivity kernels for various structure properties, including the sound speed, density, and Ledoux parameter of convective stability, and, thus, the direct relationship between the stellar properties and the deviation of observed frequencies from the reference models. The inversion results reveal significant deviations in the core structure from the reference models calculated using the MESA evolutionary code for the stellar parameters obtained by the asteroseismic model grid fitting. Our analysis shows that the best resolution of the inner helium core and surrounding shell is achieved in inversions for the Ledoux parameter.

SMD↗
NASA NTRS2017

TPSAS-NF1676L-28149-DND

An Optimal Estimation (OE) inversion method has been developed to retrieve vertical profiles of aerosol size distributions, volume/surface concentrations, complex refractive indices, and single scattering albedos using the High Spectral Resolution Lidar (HSRL) and polarimeter measurements. By combining the active and passive remote sensing data, we can take advantages of the high vertical resolution of the lidar measurements and the high information content of the combined observations. There are three components to the new retrieval system: 1) simultaneous aerosol vertical profile retrieval using HSRL data, 2) standalone aerosol retrieval using Research Scanning Polarimeter (RSP), and 3) combined HSRL and RSP optimal estimation retrieval. For the LIDAR retrieval, we solve for number concentrations, mode fractions, aerosol size distributions, and real and imaginary part of refractive indices simultaneously. Based on this information, other parameters such as effective particle radius, surface and volume concentration, single scattering albedo can be derived. The combined HSRL-polarimeter algorithm will provide improved single scattering albedo and refractive index retrievals. The newly developed retrieval algorithm is very fast, and it provides an optimal estimation solution with associated error estimates. By solving for the aerosol vertical profiles using HSRL extinctions and backscattering coefficients from different layers, the new OE method can effectively reduce the impact of measurement noise on the retrieved aerosol microphysical properties. It also enables us to effectively combing HSRL and RSP data in the single OE retrieval system since the passive polarimeter instruments measure path-integrated intensity and polarizations. We have performed sensitivity studies using simulated data and applied the OE algorithms to NASA airborne field data measured by the HSRL-2 and the RSP instruments.

Xu Liu↗
NASA NTRS2021

Computational Inference of Vibratory System with Incomplete Modal Information Using Parallel, Interactive and Adaptive Markov Chains

Inverse analysis of vibratory system is an important subject in fault identification, model updating, and robust design and control. It is challenging subject because 1) the problem is oftentimes underdetermined while the measurements are limited and/or incomplete; 2) many combinations of parameters may yield results that are similar with respect to actual response measurements; and 3) uncertainties inevitably exist. The aim of this research is to leverage upon computational intelligence through statistical inference to facilitate an enhanced, probabilistic framework using incomplete modal response measurement. This new framework is built upon efficient inverse identification through optimization, whereas Bayesian inference is employed to account for the effect of uncertainties. To overcome the computational cost barrier, we adopt Markov chain Monte Carlo (MCMC) to characterize the target function/distribution. Instead of using single Markov chain in conventional Bayesian approach, we develop a new sampling theory with multiple parallel, interactive and adaptive Markov chains and incorporate it into Bayesian inference. This can harness the collective power of these Markov chains to realize the concurrent search of multiple local optima. The number of required Markov chains and their respective initial model parameters are automatically determined via Monte Carlo simulation-based sample pre-screening followed by K-means clustering analysis. These enhancements can effectively address the aforementioned challenges in finite element inverse analysis. The validity of this framework is systematically demonstrated through case studies.

K Zhou↗
NASA NTRS1972

Application of a stochastic inverse to the geophysical inverse problem

The inverse problem for gross earth data can be reduced to an undertermined linear system of integral equations of the first kind. A theory is discussed for computing particular solutions to this linear system based on the stochastic inverse theory presented by Franklin. The stochastic inverse is derived and related to the generalized inverse of Penrose and Moore. A Backus-Gilbert type tradeoff curve is constructed for the problem of estimating the solution to the linear system in the presence of noise. It is shown that the stochastic inverse represents an optimal point on this tradeoff curve. A useful form of the solution autocorrelation operator as a member of a one-parameter family of smoothing operators is derived.

Jordan, T. H.↗
NASA NTRS1989

Theory of microemulsions in a gravitational field

A theory of microemulsions developed previously is extended to include the effect of a gravitational field. It predicts variation with position of drop size, drop volume fraction, and area per molecule in the surfactant films within a microemulsion phase. Variation in volume fraction is greatest and occurs in such a way that oil content increases with increasing elevation, as has been found experimentally. Large composition variations are predicted within a middle phase microemulsion near optimal conditions because inversion from the water-continuous to the oil-continuous arrangement occurs with increasing elevation. Generally speaking, gravity reduces solubilization within microemulsions and promotes separation of excess phases.

Jeng, J. F.↗
NASA NTRS1991

Hypersonic waverider analysis - A comparison of numerical and experimental results

In this study two waverider configurations are analyzed numerically and experimentally. The geometries are designed using an inverse design and optimization code with the optimization procedure driven by maximized L/D and limited by geometric constraints. Euler and Navier-Stokes numerical simulations are performed with a partially flux-split finite-difference code developed at NASA Ames, F3D. Experimental data is obtained at NASA Langley in the Unitary Plan Wind Tunnel (UPWT) and the 20-Inch Mach 6 wind tunnel for flow visualization, force/moment, and surface pressure measurements at both on- and off-design conditions. Euler simulations are compared with analytic results obtained from conical flow theory, and Navier-Stokes simulations are compared with experimental results both qualitatively and quantitatively for on- and off-design conditions.

Jones, K. D.↗
NASA NTRS1998

Deterministic Reconfigurable Control Design for the X-33 Vehicle

In the event of a control surface failure, the purpose of a reconfigurable control system is to redistribute the control effort among the remaining working surfaces such that satisfactory stability and performance are retained. Four reconfigurable control design methods were investigated for the X-33 vehicle: Redistributed Pseudo-Inverse, General Constrained Optimization, Automated Failure Dependent Gain Schedule, and an Off-line Nonlinear General Constrained Optimization. The Off-line Nonlinear General Constrained Optimization approach was chosen for implementation on the X-33. Two example failures are shown, a right outboard elevon jam at 25 deg. at a Mach 3 entry condition, and a left rudder jam at 30 degrees. Note however, that reconfigurable control laws have been designed for the entire flight envelope. Comparisons between responses with the nominal controller and reconfigurable controllers show the benefits of reconfiguration. Single jam aerosurface failures were considered, and failure detection and identification is considered accomplished in the actuator controller. The X-33 flight control system will incorporate reconfigurable flight control in the baseline system.

Wagner, Elaine A.↗
NASA NTRS2004

Applications of Response Surface-Based Methods to Noise Analysis in the Conceptual Design of Revolutionary Aircraft

Due to the growing problem of noise in today's air transportation system, there have arisen needs to incorporate noise considerations in the conceptual design of revolutionary aircraft. Through the use of response surfaces, complex noise models may be converted into polynomial equations for rapid and simplified evaluation. This conversion allows many of the commonly used response surface-based trade space exploration methods to be applied to noise analysis. This methodology is demonstrated using a noise model of a notional 300 passenger Blended-Wing-Body (BWB) transport. Response surfaces are created relating source noise levels of the BWB vehicle to its corresponding FAR-36 certification noise levels and the resulting trade space is explored. Methods demonstrated include: single point analysis, parametric study, an optimization technique for inverse analysis, sensitivity studies, and probabilistic analysis. Extended applications of response surface-based methods in noise analysis are also discussed.

Hill, Geoffrey A.↗