Engineering Papers⌕ Search

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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Multiwavelength and Statistical Research in Space Astrophysics

The accomplishments in the following three research areas are summarized: multiwavelength study of active galactic nuclei; magnetic activity of young stellar objects; and statistical methodology for astronomical data analysis. The research is largely based on observations of the ROSAT and ASCA X-ray observatories, complemented by ground-based optical and radio studies. Major findings include: discovery of inverse Compton X-ray emission from radio galaxy lobes; creation of the largest and least biased available sample of BL Lac objects; characterization of X-ray and nonthermal radio emission from T Tauri stars; obtaining an improved census of young stars in a star forming region and modeling the star formation history and kinematics; discovery of X-ray emission from protostars; development of linear regression methods and codes for interpreting astronomical data; and organization of the first cross-disciplinary conferences for astronomers and statisticians.

Feigelson, Eric D.↗
NASA NTRS2011

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu↗
NASA NTRS1992

High-resolution gravity model of Venus

The anomalous gravity field of Venus shows high correlation with surface features revealed by radar. We extract gravity models from the Doppler tracking data from the Pioneer Venus Orbiter by means of a two-step process. In the first step, we solve the nonlinear spacecraft state estimation problem using a Kalman filter-smoother. The Kalman filter has been evaluated through simulations. This evaluation and some unusual features of the filter are discussed. In the second step, we perform a geophysical inversion using a linear Bayesian estimator. To allow an unbiased comparison between gravity and topography, we use a simulation technique to smooth and distort the radar topographic data so as to yield maps having the same characteristics as our gravity maps. The maps presented cover 2/3 of the surface of Venus and display the strong topography-gravity correlation previously reported. The topography-gravity scatter plots show two distinct trends.

Reasenberg, R. D.↗
NASA NTRS1993

Venus gravity

The anomalous gravity field of Venus shows high correlation with surface features revealed by radar. We extract gravity models from the Doppler tracking data from the Pioneer Venus Orbiter (PVO) by means of a two-step process. In the first step, we solve the nonlinear spacecraft state estimation problem using a Kalman filter-smoother. The Kalman filter was evaluated through simulations. This evaluation and some unusual features of the filter are discussed. In the second step, we perform a geophysical inversion using a linear Bayesian estimator. To allow an unbiased comparison between gravity and topography, we use a simulation technique to smooth and distort the radar topographic data so as to yield maps having the same characteristics as our gravity maps. The maps presented cover 2/3 of the surface of Venus and display the strong topography-gravity correlation previously reported. The topography-gravity scatter plots show two distinct trends.

Reasenberg, Robert D.↗
NASA NTRS2013

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps to reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co‐located at one corner. Due to this limited geographical extent, near‐surface ray coverage from moonquakes is low, but increases with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefore well suited to joint inversion, since the complementary information reduces inherent model ambiguity. Previous joint inversions of the Apollo seismic data (seismic phase arrival times) and Clementine‐ or Lunar Prospector‐derived gravity data (mass and moment of inertia) attempted to recover the subsurface structure of the Moon by focusing on hypothetical lunar compositions that explore the density/velocity relationship. These efforts typically search for the best fitting thermodynamically calculated velocity/density model, allowing variables like core size, velocity, and/or composition to vary freely. Seismic velocity profiles previously derived from the Apollo seismic data through inversion of travel times vary both in the depth of the crust and mantle layers, and the seismic velocities and densities assigned to those layers. The lunar mass and moment of inertia likewise only constrain gross variations in the density profile beyond that of a uniform density sphere. As a result, composition and structure models previously obtained by jointly inverting these data retain the original uncertainties inherent in the input data sets. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by the GRAIL lunar gravity mission, in order to recover seismic velocities and density as a function of latitude, longitude, and depth within the Moon. We will relate density to seismic velocity using a linear relationship that is allowed to be depth‐dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of density, velocity, and B‐coefficient perturbations that minimize (in a least‐squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗
NASA NTRS2008

Radiative Transfer and Retrievals in EOF Domain

The Infrared Atmospheric Sounding Interferometer (IASI) is a hyperspectral sensor with 8461 spectral channels and a nominal spectral resolution of 0.25 cm(sup -1). It is computationally intensive to perform radiative transfer calculations and inversions using all these channels. We will present a Principal Component-based Radiative Transfer Model (PCRTM) and a retrieval algorithm which perform all the necessary calculations in EOF domain. Since the EOFs are orthogonal to each other, only about 100 principal components are needed to represent the information content of the 8461 channels. The PCRTM provides the EOF coefficients and associated derivatives with respect to atmospheric and surface parameters needed by the inversion algorithm. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance and a priori information as constraints. The retrieved parameters include atmospheric temperature, moisture and ozone profiles, cloud parameters, surface skin temperature, and surface emissivities. To make the retrieval system even more compact and stable. The atmospheric vertical profiles are compressed into the EOF space as well. The surface emissivities are also compressed into EOF space.

Liu, Xu↗
NASA NTRS1989

Absorption line profiles in a moving atmosphere - A single scattering linear perturbation theory

An integral equation is derived which linearly relates Doppler perturbations in the spectrum of atmospheric absorption features to the wind system which creates them. The perturbation theory is developed using a single scattering model, which is validated against a multiple scattering calculation. The nature and basic properties of the kernels in the integral equation are examined. It is concluded that the kernels are well behaved and that wind velocity profiles can be recovered using standard inversion techniques.

Hays, P. B.↗
NASA NTRS1990

Scattering from coated structures and antenna pattern control using impedance surfaces, part A/B

The scattering from coated, conducting structures, specifically the coated dihedral corner reflector configuration and the coated strip/plate configuration is examined. The formulation uses impedance-wedge Uniform Theory of Diffraction scattering coefficients to calculate the diffracted fields. A finite-thickness coating is approximated using the impedance boundary condition to arrive at an equivalent impedance for the coating. The formulation of the impedance wedge coefficients is outlined. Far-field, perfectly conducting approximations are discussed. Problems with the present dihedral corner reflector model for certain angles of incidence and observation are discussed along with a potentially rectifying modification. Also, the capacity to measure the electromagnetic properties of lossy materials was developed. The effects of using multiple material coatings on the radiation pattern of the horn antenna were studied. Numerous computations were devoted toward the inverse problem of synthesizing desired radiation patterns using the impedance surfaces. Stabilizing the equivalent sheet impedance using the linear control condition was attempted, and it was found to be a very difficult task.

Balanis, Constantine A.↗
NASA NTRS2009

2-Micron Laser Transmitter for Coherent CO2 DIAL Measurement

Carbon dioxide (CO2) has been recognized as one of the most important greenhouse gases. It is essential for the study of global warming to accurately measure the CO2 concentration in the atmosphere and continuously record its variation. A high repetition rate, highly efficient, Q-switched 2-micron laser system as the transmitter of a coherent differential absorption lidar for CO2 measurement has been developed in NASA Langley Research Center. This laser system is capable of making a vertical profiling of CO2 from ground and column measurement of CO2 from air and space-borne platform. The transmitter is a master-slave laser system. The master laser operates in a single frequency, either on-line or off-line of a selected CO2 absorption line. The slave laser is a Q-switched ring-cavity Ho:YLF laser which is pumped by a Tm:fiber laser. The repetition rate can be adjusted from a few hundred Hz to 10 kHz. The injection seeding success rate is from 99.4% to 99.95%. For 1 kHz operation, the output pulse energy is 5.5mJ with the pulse length of ~50 ns. The optical-to-optical efficiency is 39% when the pump power is 14.5W. A Ho:YLF laser operating in the range of 2.05 micrometers can be tuned over several characteristic lines of CO2 absorption. Experimentally, a diode pumped Ho:Tm:YLF laser has been successfully used as the transmitter of coherent differential absorption lidar for the measurement of CO2 with a repetition rate of 5 Hz and pulse energy of 75 mJ. For coherent detection, high repetition rate is required for speckle averaging to obtain highly precise measurements. However, a diode pumped Ho:Tm:YLF laser can not operate in high repetition rate due to the large heat loading and up-conversion. A Tm:fiber laser pumped Ho:YLF laser with low heat loading can operate in high repetition rate. A theoretical model has been established to simulate the performance of Tm:fiber laser pumped Ho:YLF lasers. For continuous wave (CW) operation, high pump intensity with small beam size is suitable for high efficiency. For Q-switched operation, the optimal energy extraction relies on the pump intensity, pump volume, and pump duration which is inversely proportion to the repetition rate. CW and Q-switched Ho:YLF lasers with different linear cavity configurations have been designed and demonstrated for a 30 W Tm:fiber pump laser. The CW Ho laser slope efficiency and optical-to-optical efficiencies reach 65% and 55%, respectively. The pulsed laser efficiency depends on the repetition rate. For 1 kHz operation, the optical-to-optical efficiency is 39% when the pump power is 14.5W. Currently, the injection seeding success rate is between 99.4% and 99.95%. After a ten thousand pulses, the standard deviation of the laser frequency jitter is about 3 MHz. It meets the requirements of highly precise CO2 concentration measurement. In conclusion, an injection seeded, high repetition rate, Q-switched Ho:YLF laser has been developed for a coherent CO2 differential absorption lidar. This master-slave laser system has high optical-to-optical efficiency and seeding success rate. It can potentially meet the requirements of the coherent detection of CO2 concentration by a differential absorption lidar technique.

Singh, Upendra N.↗
NASA NTRS2014

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increases with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefore well suited to joint inversion, since the complementary information reduces inherent model ambiguity. Previous joint inversions of the Apollo seismic data (seismic phase arrival times) and Clementine- or Lunar Prospector-derived gravity data (mass and moment of inertia) attempted to recover the subsurface structure of the Moon by focusing on hypothetical lunar compositions that explored the density/velocity relationship. These efforts typically searched for the best fitting thermodynamically calculated velocity/density model, and allowed variables like core size, velocity, and/or composition to vary freely. Seismic velocity profiles derived from the Apollo seismic data through travel time inversion vary both in the depth of the crust and mantle layers, and the seismic velocities and densities assigned to those layers. The lunar mass and moment of inertia likewise only constrain gross variations in the density profile beyond that of a uniform density sphere. As a result, composition and structure models previously obtained by jointly inverting these data retain the original uncertainties inherent in the input data sets. We perform a joint inversion of Apollo seismic delay times and gravity data collected by the GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude, and depth within the Moon. We relate density (p) to seismic velocity (v) using a depth-dependent linear relationship. The corresponding coefficient (B) can reflect a variety of material properties, including temperature and composition. The inversion seeks to recover the set of p, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗
NASA NTRS1986

A simple algorithm for computing canonical forms

It is well known that all linear time-invariant controllable systems can be transformed to Brunovsky canonical form by a transformation consisting only of coordinate changes and linear feedback. However, the actual procedures for doing this have tended to be overly complex. The technique introduced here is envisioned as an on-line procedure and is inspired by George Meyer's tangent model for nonlinear systems. The process utilizes Meyer's block triangular form as an intermedicate step in going to Brunovsky form. The method also involves orthogonal matrices, thus eliminating the need for the computation of matrix inverses. In addition, the Kronecker indices can be computed as a by-product of this transformation so it is necessary to know them in advance.

Ford, H.↗
NASA NTRS2000

A Reconfiguration Scheme for Accommodating Actuator Failures in Multi-Input, Multi-Output Flight Control Systems

A multi-input, multi-output reconfigurable flight control system design utilizing a robust controller and an adaptive filter is presented. The robust control design consists of a reduced-order, linear dynamic inversion controller with an outer-loop compensation matrix derived from Quantitative Feedback Theory (QFT). A principle feature of the scheme is placement of the adaptive filter in series with the QFT compensator thus exploiting the inherent robustness of the nominal flight control system in the presence of plant uncertainties. An example of the scheme is presented in a pilot-in-the-loop computer simulation using a simplified model of the lateral-directional dynamics of the NASA F18 High Angle of Attack Research Vehicle (HARV) that included nonlinear anti-wind up logic and actuator limitations. Prediction of handling qualities and pilot-induced oscillation tendencies in the presence of these nonlinearities is included in the example.

Siwakosit, W.↗
NASA NTRS2016

A Method for Improving Hotspot Directional Signatures in BRDF Models Used for MODIS

The semi-empirical, kernel-driven, linear RossThick-LiSparseReciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is used to generate the routine MODIS BRDFAlbedo product due to its global applicability and the underlying physics. A challenge of this model in regard to surface reflectance anisotropy effects comes from its underestimation of the directional reflectance signatures near the Sun illumination direction; also known as the hotspot effect. In this study, a method has been developed for improving the ability of the RTLSR model to simulate the magnitude and width of the hotspot effect. The method corrects the volumetric scattering component of the RTLSR model using an exponential approximation of a physical hotspot kernel, which recreates the hotspot magnitude and width using two free parameters (C(sub 1) and C(sub 2), respectively). The approach allows one to reconstruct, with reasonable accuracy, the hotspot effect by adjusting or using the prior values of these two hotspot variables. Our results demonstrate that: (1) significant improvements in capturing hotspot effect can be made to this method by using the inverted hotspot parameters; (2) the reciprocal nature allow this method to be more adaptive for simulating the hotspot height and width with high accuracy, especially in cases where hotspot signatures are available; and (3) while the new approach is consistent with the heritage RTLSR model inversion used to estimate intrinsic narrowband and broadband albedos, it presents some differences for vegetation clumping index (CI) retrievals. With the hotspot-related model parameters determined a priori, this method offers improved performance for various ecological remote sensing applications; including the estimation of canopy structure parameters.

airborne measurements↗
NASA NTRS2004

Optimal Aerosol Parameterization for Remote Sensing Retrievals

We have developed a new algorithm for the retrieval of aerosol and gases from SAGE It1 solar transmission measurements. This algorithm improves upon the NASA operational algorithm in several key aspects, including solving the problem non-linearly and incorporating a new methodology for separating the contribution of aerosols and gases. In order to extract aerosol information we have built a huge database of aerosol models for both stratospheric and tropospheric aerosols, and polar stratospheric cloud particles. This set of models allows us to calculate a vast range of possible extinction spectra for aerosols. and from these, derive a set of eigenvectors which then provide the basis set used in our inversion algorithm. Our aerosol algorithm and retrievals are described in several articles (listed in References Section) published under this grant. In particular they allow us to analyze the spectral properties of aerosols and PSCs and ultimately derive their microphysical properties. We have found some considerable differences between our spectra and the ones derived from the SAGE III operational algorithm. These are interesting as they provide an independent check on the validity of published aerosol data and, in particular, on their associated uncertainties. In order to understand these differences, we are assembling independent aerosol data from other sources with which to make comparisons. We have carried out extensive comparisons of our ozone retrievals with both SAGE III and independent lidar, ozonesonde, and satellite measurements (Polyakov et al., 2004). These show very good agreement throughout the stratosphere and help to quantify differences which can be attributed to natural variation in ozone versus that produced by algorithmic differences. In the mid - upper stratosphere, agreement with independent data was generally within 5 - 20%. but in the lower stratosphere the differences were considerably larger. We believe that a large proportion of this discrepancy in the lower stratosphere is attributable to natural variation, and is also seen in comparisons between lidar and ozonesonde measurements. NO2 profiles obtained with our algorithm were compared to those obtained through the SAGE III operational algorithm and exhibited differences of 20 - 40%. Our retrieved profiles agree with the HALOE NO2 measurements significantly better than those of the operational retrieval. In other work (described below), we are extending our aerosol retrievals into the infrared regime and plan to perform retrievals from combined uv-visible-infrared spectra. This work will allow us to use the spectra to derive the size and composition of aerosols, and we plan to employ our algorithms in the analysis of PSC spectra. We are presently also developing a limb-scattering algorithm to retrieve aerosol data from limb measurements of solar scattered radiation.

Newchurch, Michael J.↗
NASA NTRS1984

Algorithms for solving large sparse systems of simultaneous linear equations on vector processors

Very efficient algorithms for solving large sparse systems of simultaneous linear equations have been developed for serial processing computers. These involve a reordering of matrix rows and columns in order to obtain a near triangular pattern of nonzero elements. Then an LU factorization is developed to represent the matrix inverse in terms of a sequence of elementary Gaussian eliminations, or pivots. In this paper it is shown how these algorithms are adapted for efficient implementation on vector processors. Results obtained on the CYBER 200 Model 205 are presented for a series of large test problems which show the comparative advantages of the triangularization and vector processing algorithms.

David, R. E.↗
NASA NTRS1991

Application of nonlinear feedback control theory to supermaneuverable aircraft

Controlled flight at extremely high angles of attack, far exceeding the stall angle, and/or at high angular rates is sometimes referred to as supermaneuvering flight. The objective was to examine methods for design of control laws for aircraft performing supermaneuvers. Since the equations which govern the motion of aircraft during supermaneuvers are nonlinear, this study concentrated on nonlinear control law design procedures. The two nonlinear techniques considered were Nonlinear Quadratic Regulator (NLQR) theory and nonlinear dynamic inversion. A conventional gain scheduled proportional plus integral (P + I) controller was also developed to serve as a baseline design typical of current control laws used in aircraft. A mathematical model of a generic supermaneuverable aircraft was developed from data obtained from the literature. A detailed computer simulation of the aircraft was also developed. This simulation allowed the flying of proposed supermaneuvers and was used to evaluate the performance of the control law designs and to generate linearized models of the aircraft at different flight conditions.

Garrard, William L.↗
NASA NTRS1992

Commentary 2 to Cox and Little: radiation-induced oncogenic transformation: the interplay between dose, dose protraction, and radiation quality

There is now a substantial body of evidence for end points such as oncogenic transformation in vitro, and carcinogenesis and life shortening in vivo, suggesting that dose protraction leads to an increase in effectiveness relative to a single, acute exposure--at least for radiations of medium linear energy transfer (LET) such as neutrons. Table I contains a summary of the pertinent data from studies in which the effect is seen. [table: see text] This phenomenon has come to be known as the "inverse dose rate effect," because it is in marked contrast to the situation at low LET, where protraction in delivery of a dose of radiation, either by fractionation or low dose rate, results in a decreased biological effect; additionally, at medium and high LET, for radiobiological end points such as clonogenic survival, the biological effectiveness is independent of protraction. The quantity and quality of the published reports on the "inverse dose rate effect" leaves little doubt that the effect is real, but the available evidence indicates that the magnitude of the effect is due to a complex interplay between dose, dose rate, and radiation quality. Here, we first summarize the available data on the inverse dose rate effect and suggest that it follows a consistent pattern in regard to dose, dose rate, and radiation quality; second, we describe a model that predicts these features; and, finally, we describe the significance of the effect for radiation protection.

Review, Tutorial↗
NASA NTRS1989

Inverse dynamics of a 3 degree of freedom spatial flexible manipulator

A technique is presented for solving the inverse dynamics and kinematics of 3 degree of freedom spatial flexible manipulator. The proposed method finds the joint torques necessary to produce a specified end effector motion. Since the inverse dynamic problem in elastic manipulators is closely coupled to the inverse kinematic problem, the solution of the first also renders the displacements and rotations at any point of the manipulator, including the joints. Furthermore the formulation is complete in the sense that it includes all the nonlinear terms due to the large rotation of the links. The Timoshenko beam theory is used to model the elastic characteristics, and the resulting equations of motion are discretized using the finite element method. An iterative solution scheme is proposed that relies on local linearization of the problem. The solution of each linearization is carried out in the frequency domain. The performance and capabilities of this technique are tested through simulation analysis. Results show the potential use of this method for the smooth motion control of space telerobots.

Bayo, Eduardo↗