SEARCH · Engineering Papers

Results for “linear 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗

NASA NTRS2004

Pilots Rate Augmented Generalized Predictive Control for Reconfiguration

The objective of this paper is to report the results from the research being conducted in reconfigurable fight controls at NASA Ames. A study was conducted with three NASA Dryden test pilots to evaluate two approaches of reconfiguring an aircraft's control system when failures occur in the control surfaces and engine. NASA Ames is investigating both a Neural Generalized Predictive Control scheme and a Neural Network based Dynamic Inverse controller. This paper highlights the Predictive Control scheme where a simple augmentation to reduce zero steady-state error led to the neural network predictor model becoming redundant for the task. Instead of using a neural network predictor model, a nominal single point linear model was used and then augmented with an error corrector. This paper shows that the Generalized Predictive Controller and the Dynamic Inverse Neural Network controller perform equally well at reconfiguration, but with less rate requirements from the actuators. Also presented are the pilot ratings for each controller for various failure scenarios and two samples of the required control actuation during reconfiguration. Finally, the paper concludes by stepping through the Generalized Predictive Control's reconfiguration process for an elevator failure.

Soloway, Don↗

NASA NTRS1999

Mars Crustal Remanent Magnetism: An Extinct Dynamo Leaves a Record of Field Reversals in the Heavily Cratered Highlands

The Mars Global Surveyor spacecraft, in a highly elliptical polar orbit about Mars, obtained vector magnetic field measurements just above the surface of Mars (altitudes > 100 kilometers). Crustal magnetization, largely confined to the most ancient, heavily cratered Mars highlands, is frequently organized in east-west trending linear features, the largest of which extends over 2000 km. A representative set of survey passes are modeled using uniformly magnetized thin plates and a generalized inverse methodology. Crustal remanent magnetization exceeds that deduced for the largest terrestrial magnetic anomalies by more than an order of magnitude. Groups of quasi-parallel linear features of alternating magnetic polarity are found. They are reminiscent of similar magnetic features associated with sea floor spreading and crustal genesis on Earth but with a much larger spatial scale.

Connerney, John E.↗

NASA NTRS1990

Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data

Inverse methods for estimating the surface ciculation of the equatorial Pacific by combining a linear reduced-gravity shallow-water model with the Tropical Ocean-Global Atmosphere ship-of-opportunity expendable bathythermograph (TOGA SOP XBT) observing program are examined. It is demonstrated that a simple linear model of the upper circulation of the equatorial Pacific can be successfully used as a weak constraint when smoothing the TOGA SOP XBT data. A circulation is sought as the weighted least squares fit to the dynamics and the data. The solution method is an expansion in representer functions, and the generalized inverse problem is thereby reduced from a functional problem to an algebraic problem for the coefficients of the representer. A specific inverse calculation using synthetic forcing and data is presented.

Bennett, Andrew F.↗

NASA NTRS1995

Analyzing rocket plume spectral data with neural networks

The Optical Plume Anomaly Detection (OPAD) system is under development to provide early-warning failure detection in support of ground-level testing of the Space Shuttle Main Engine (SSME). Failure detection is to be achieved through the acquisition of spectrally resolved plume emissions and subsequent identification of abnormal levels indicative of engine corrosion or component failure. Two computer codes (one linear and the other non-linear) are used by the OPAD system to iteratively determine specific element concentrations in the SSME plume, given emission intensity and wavelength information. Since this analysis is extremely labor intensive, a study was initiated to develop neural networks that would model the 'inverse' of these computer codes. Optimally connected feed-forward networks with imperceptible prediction error have been developed for each element modeled by the linear code, SPECTRA4. Radial basis function networks were developed for the non-linear code, SPECTRA5, and predict combustion temperature in addition to element concentrations.

Whitaker, Kevin W.↗

NASA NTRS1979

Total aircraft flight-control system - Balanced open- and closed-loop control with dynamic trim maps

The availability of the airborne digital computer has made possible a Total Aircraft Flight Control System (TAFCOS) that uses virtually the complete nonlinear propulsive and aerodynamic data for the aircraft to construct dynamic trim maps that represent an inversion of the aircraft model. The trim maps, in series with the aircraft, provide essentially a linear feed-forward path. Basically, open-loop trajectory control is employed with only a small perturbation feedback signal required to compensate for inaccuracy in the aircraft model and for external disturbances. Simulation results for application to an automatic carrier-landing system are presented. Flight-test results for a STOL aircraft operating automatically over a major portion of its flight regime are presented. The concept promises a more rapid and straightforward design from aerodynamic principles, particularly for highly nonlinear configurations, and requires substantially less digital computer capacity than conventional automatic flight-control system designs.

Smith, G. A.↗

NASA NTRS1994

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗

NASA NTRS1992

Relations between surface conductance and spectral vegetation indices at intermediate (100 m sq to 15 km sq) length scales

The relationship between surface conductance and spectral vegetation indices is investigated utilizing the FIFE data set, principally the surface flux station data and images from the TM instrument. It is found that the unstressed canopy conductance for a given site for a given day is near-linearly related to the incident PAR flux. Estimates of unstressed canopy conductance were acquired via a model inversion that separated the soil and vegetation contributions to evapotranspiration and made adjustments for the effects of vapor pressure deficit and soil moisture stress.

Sellers, Piers J.↗

NASA NTRS2012

Iterative Atmospheric Correction Scheme and the Polarization Color of Alpine Snow

Characterization of the Earth's surface is crucial to remote sensing, both to map geomorphological features and because subtracting this signal is essential during retrievals of the atmospheric constituents located between the surface and the sensor. Current operational algorithms model the surface total reflectance through a weighted linear combination of a few geometry-dependent kernels, each devised to describe a particular scattering mechanism. The information content of these measurements is overwhelmed by that of instruments with polarization capabilities: proposed models in this case are based on the Fresnel reflectance of an isotropic distribution of facets. Because of its remarkable lack of spectral contrast, the polarized reflectance of land surfaces in the shortwave infrared spectral region, where atmospheric scattering is minimal, can be used to model the surface also at shorter wavelengths, where aerosol retrievals are attempted based on well-established scattering theories. In radiative transfer simulations, straightforward separation of the surface and atmospheric contributions is not possible without approximations because of the coupling introduced by multiple reflections. Within a general inversion framework, the problem can be eliminated by linearizing the radiative transfer calculation, and making the Jacobian (i.e., the derivative expressing the sensitivity of the reflectance with respect to model parameters) available at output. We present a general methodology based on a Gauss-Newton iterative search, which automates this procedure and eliminates de facto the need of an ad hoc atmospheric correction. In this case study we analyze the color variations in the polarized reflectance measured by the NASA Goddard Institute of Space Studies Research Scanning Polarimeter during a survey of late-season snowfields in the High Sierra. This insofar unique dataset presents challenges linked to the rugged topography associated with the alpine environment and a likely high water content due to melting. The analysis benefits from ancillary information provided by the NASA Langley High Spectral Resolution Lidar deployed on the same aircraft. The results obtained from the iterative scheme are contrasted against the surface polarized reflectance obtained ignoring multiple reflections, via the simplistic subtraction of the atmospheric scattering contribution. Finally, the retrieved reflectance is modeled after the scattering properties of a dense collection of ice crystals at the surface. Confirming that the polarized reflectance of snow is spectrally flat would allow to extend the techniques already in use for polarimetric retrievals of aerosol properties over land to the large portion of snow-covered pixels plaguing orbital and suborbital observations.

Polarized BRDF↗

NASA NTRS2015

Retrieval of Aerosol Microphysical Properties from AERONET Photo-Polarimetric Measurements. 2: A New Research Algorithm and Case Demonstration

A new research algorithm is presented here as the second part of a two-part study to retrieve aerosol microphysical properties from the multispectral and multiangular photopolarimetric measurements taken by Aerosol Robotic Network's (AERONET's) new-generation Sun photometer. The algorithm uses an advanced UNified and Linearized Vector Radiative Transfer Model and incorporates a statistical optimization approach.While the new algorithmhas heritage from AERONET operational inversion algorithm in constraining a priori and retrieval smoothness, it has two new features. First, the new algorithmretrieves the effective radius, effective variance, and total volume of aerosols associated with a continuous bimodal particle size distribution (PSD) function, while the AERONET operational algorithm retrieves aerosol volume over 22 size bins. Second, our algorithm retrieves complex refractive indices for both fine and coarsemodes,while the AERONET operational algorithm assumes a size-independent aerosol refractive index. Mode-resolved refractive indices can improve the estimate of the single-scattering albedo (SSA) for each aerosol mode and thus facilitate the validation of satellite products and chemistry transport models. We applied the algorithm to a suite of real cases over Beijing_RADI site and found that our retrievals are overall consistent with AERONET operational inversions but can offer mode-resolved refractive index and SSA with acceptable accuracy for the aerosol composed by spherical particles. Along with the retrieval using both radiance and polarization, we also performed radiance-only retrieval to demonstrate the improvements by adding polarization in the inversion. Contrast analysis indicates that with polarization, retrieval error can be reduced by over 50% in PSD parameters, 10-30% in the refractive index, and 10-40% in SSA, which is consistent with theoretical analysis presented in the companion paper of this two-part study.

aerosol retrieval algorithm↗

NASA NTRS1996

High Resolution Doppler Lidar

This Grant supported the development of an incoherent lidar system to measure winds and aerosols in the lower atmosphere. During this period the following activities occurred: (1) an active feedback system was developed to improve the laser frequency stability; (2) a detailed forward model of the instrument was developed to take into account many subtle effects, such as detector non-linearity; (3) a non-linear least squares inversion method was developed to recover the Doppler shift and aerosol backscatter without requiring assumptions about the molecular component of the signal; (4) a study was done of the effects of systematic errors due to multiple etalon misalignment. It was discovered that even for small offsets and high aerosol loadings, the wind determination can be biased by as much as 1 m/s. The forward model and inversion process were modified to account for this effect; and (5) the lidar measurements were validated using rawinsonde balloon measurements. The measurements were found to be in agreement within 1-2 m/s.

Source record↗

NASA NTRS1976

The evolution of the clear air convective layer revealed by surface-based remote sensors

Results are reported for simultaneous observations of the growth and decay of the clear-air convective mixing layer near a coastline, which were made with an FM-CW radar, a high-power narrow-beam S-band radar, and an acoustic echo sounder. The main purpose of this study was to determine the relationship between the rise rate of the convective depth and the lapse rate of temperature, particularly in the morning hours. The results indicate that the three remote sensors can provide excellent mutually supporting data on the convective depth. It is found that this depth is well behaved during the day and that its rise rate varies roughly linearly with the inverse square root of the temperature lapse rate during the morning. The data suggest that some models concerning the rise rate require modification, since these models imply that the surface heat flux would have to be unreasonably large to produce the observed relationship.

Noonkester, V. R.↗

NASA NTRS1976

Alternative methods for the design of jet engine control systems

Various alternatives to linear quadratic design methods for jet engine control systems are discussed. The main alternatives are classified into two broad categories: nonlinear global mathematical programming methods and linear local multivariable frequency domain methods. Specific studies within these categories include model reduction, the eigenvalue locus method, the inverse Nyquist method, polynomial design, dynamic programming, and conjugate gradient approaches.

Sain, M. K.↗

NASA NTRS2019

Lidar-Polarimeter Retrieval OSSEs in Support of NASA's Aerosols and Clouds-Convection-Precipitation (ACCP) Study

The 2017 Decadal Survey (DS) highlighted Earth System Science themes, science and application questions, and several high priority objectives that have led to the inclusion of Aerosols (A) and Clouds-Convection-Precipitation (CCP) as Designated Observables (DOs). On June 1, 2018, several NASA centers (GSFC, LaRC, JPL, MSFC, GRC and ARC) submitted a joint Study Plan to the NASA Earth Science Division for the Aerosol (A) and Cloud, Convection, and Precipitation (CCP) Pre-formulation Study (A-CCP). The DS and the A-CCP team recognized the science merit in combining the A and CCP DOs for both enhancing the ability to address a number of science objectives and also to provide an expanded capability to address additional objectives beyond those addressed by individual DOs.A critical element of the A-CCP observing strategy is to make extensive use of new passive and active sensors as well as of the so-called Program-of-Record (PoR), complemented by a fully integrated sub-orbital component. Central to this observing system design is the adoption of a Value Framework in which quantitative assessment of the science benefits of space-and air-borne assets is a key element. Given pre-defined A-CCP science objectives and geophysical variables with desired accuracies, A-CCP relies on a spectrum of Observing System Simulation Experiments (OSSEs) aimed at addressing pixel level retrieval uncertainties and sampling trade-offs. In this talk we will discuss a subset of Retrieval OSSEs being considered for A-CCP, namely, synergistic lidar-polarimeter retrievals of particular relevance for the A-CCP aerosol science objectives. Starting with aerosol states from the GEOS-5 Nature Run (G5NR) sampled along specific satellite orbits, we simulate polarized radiances at the desired polarimeter wavelengths with the Vector Linearized Direct Ordinate Radiative Transfer (VLIDORT) model, alongside the lidar signal for the relevant lidars with realistic error characterization. Next, inversions are performed with the Generalized Retrieval of Aerosol and Surface Properties (GRASP) system and the accuracy of the retrieved geophysical variables are assessed. In this presentation we will highlight results for key architectures being considered for A-CCP.

da Silva, Arlindo↗

NASA NTRSDate not supplied

Lidar-Polarimeter Retrieval OSSEs using a Nature Run in Support of NASA's Aerosols, Clouds, Convection and Precipitation (ACCP) Study

The 2017 Decadal Survey (DS) highlighted Earth System Science themes, science and application questions, and several high priority objectives that have led to the inclusion of Aerosols (A) and Clouds-Convection-Precipitation (CCP) as Designated Observables (DOs). On June 1, 2018, several NASA centers (GSFC, LaRC, JPL, MSFC, GRC and ARC) submitted a joint Study Plan to the NASA Earth Science Division for the Aerosol (A) and Cloud, Convection, and Precipitation (CCP) Pre-formulation Study (ACCP). The DS and the ACCP team recognized the science merit in combining the A and CCP DOs for both enhancing the ability to address a number of science objectives and also to provide an expanded capability to address additional objectives beyond those addressed by individual DOs. A critical element of the ACCP observing strategy is to make extensive use of new passive and active sensors as well as of the so-called Program-of-Record (PoR), complemented by a fully integrated sub-orbital component. Central to this observing system design is the adoption of a Value Framework in which quantitative assessment of the science benefits of space- and air-borne assets is a key element. Given pre-defined ACCP science objectives and geophysical variables with desired accuracies, ACCP relies on a spectrum of Observing System Simulation Experiments (OSSEs) aimed at addressing pixel level retrieval uncertainties and sampling trade-offs. In this talk we will discuss a subset of Retrieval OSSEs being considered for ACCP, namely, synergistic lidar-polarimeter retrievals based on observation simulations from the GEOS-5 Nature Run (G5NR). Starting with aerosol states from the GEOS-5 Nature Run (G5NR) sampled along specific satellite orbits, we simulate polarized radiances at the desired polarimeter wavelengths with the Vector Linearized Direct Ordinate Radiative Transfer (VLIDORT) model, alongside the lidar signal for the relevant lidars with realistic error characterization. Next, inversions are performed with the Generalized Retrieval of Aerosol and Surface Properties (GRASP) system and the accuracy of the retrieved geophysical variables are assessed. In this presentation we will highlight results for key architectures being considered for ACCP with emphasis on geophysical variables of particular relevance to ACCP’s aerosol science objectives.

A Da Silva↗

NASA NTRS1992

Accelerated complete-linearization method for calculating NLTE model stellar atmospheres

Two approaches to accelerating the method of complete linearization for calculating NLTE model stellar atmospheres are suggested. The first one, the so-called Kantorovich variant of the Newton-Raphson method, consists of keeping the Jacobi matrix of the system fixed, which allows us to calculate the costly matrix inversions only a few times and then keep them fixed during the subsequent computations. The second method is an application of the Ng acceleration. Both methods are extremely easy to implement with any model atmosphere code based on complete linearization. It is demonstrated that both methods, and especially their combination, yield a rapidly and globally convergent algorithm, which takes 2 to 5 times less computer time, depending on the model at hand and the required accuracy, than the ordinary complete linearization. Generally, the time gain is more significant for more complicated models. The methods were tested for a broad range of atmospheric parameters, and in all cases they exhibited similar behavior. Ng acceleration applied on the Kantorovich variant thus offers a significant improvement of the standard complete-linearization method, and may now be used for calculating relatively involved NLTE model stellar atmospheres.

Hubeny, I.↗

NASA NTRS1973

Range of earth structure nonuniqueness implied by body wave observations.

The Herglotz-Wiechert integral for the direct inversion of ray parameter versus distance curves can be manipulated to find the envelope of all possible models consistent with geometrical body wave observations (travel time and ray parameter versus distance). Such an extremal inversion approach has been used to find the uncertainty bounds for the velocity structure in the mantle and core. It is found, for example, that there is an uncertainty of plus or minus 40 km in the radius of the inner core boundary, plus or minus 18 km at the core-mantle boundary, and plus or minus 35 km at the 435-km transition zone. The velocity uncertainty is about plus or minus 0.08 km/sec for P and S waves in the lower mantle and about plus or minus 0.20 km/sec in the core. Experiments with various combinations of ray types in the core indicate that rather crude observations of SKKS-SKS travel times confine the range of possible models far more dramatically than do the most precise estimates of PmKP travel times. Comparisons of results from extremal inversion and linearized perturbation inversions indicate that body wave behavior is too strongly nonlinear for linearized schemes to be effective for predicting uncertainty.

Wiggins, R. A.↗

NASA NTRS1986

Determining the general circulation of the ocean - A nonlinear inverse problem

The usual inversion procedures for determining the oceanic circulation take the density field as given. This does not explicitly account for the mismatch between the data (which are typically a single realization of the density field) and the circulation being estimated (which may be a time mean or some other smoothed quantity), but explicit consideration of the uncertainty of the density data renders the problem nonlinear, and new estimation methods have to be used. As an example, a nonlinear inversion is successfully applied here to the simple analytical model of Stommel and Veronis (1981). It is also shown that their linear procedure, if properly used, yields a much better result than they obtained.

Mercier, H.↗