Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: ##EQU1## wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absoption information. An algorithm, which combines a two dimensional (2D) matrix inversion with a one-dimensional (1D) Fourier transform inversion is used to obtain images of three dimensional hidden objects in turbid scattering media.

Alfano, Robert R.↗
NASA NTRS2001

Computationally efficient control allocation

A computationally efficient method for calculating near-optimal solutions to the three-objective, linear control allocation problem is disclosed. The control allocation problem is that of distributing the effort of redundant control effectors to achieve some desired set of objectives. The problem is deemed linear if control effectiveness is affine with respect to the individual control effectors. The optimal solution is that which exploits the collective maximum capability of the effectors within their individual physical limits. Computational efficiency is measured by the number of floating-point operations required for solution. The method presented returned optimal solutions in more than 90% of the cases examined; non-optimal solutions returned by the method were typically much less than 1% different from optimal and the errors tended to become smaller than 0.01% as the number of controls was increased. The magnitude of the errors returned by the present method was much smaller than those that resulted from either pseudo inverse or cascaded generalized inverse solutions. The computational complexity of the method presented varied linearly with increasing numbers of controls; the number of required floating point operations increased from 5.5 i, to seven times faster than did the minimum-norm solution (the pseudoinverse), and at about the same rate as did the cascaded generalized inverse solution. The computational requirements of the method presented were much better than that of previously described facet-searching methods which increase in proportion to the square of the number of controls.

Durham, Wayne↗
NASA NTRS2008

1D-VAR Retrieval Using Superchannels

Since modern ultra-spectral remote sensors have thousands of channels, it is difficult to include all of them in a 1D-var retrieval system. We will describe a physical inversion algorithm, which includes all available channels for the atmospheric temperature, moisture, cloud, and surface parameter retrievals. Both the forward model and the inversion algorithm compress the channel radiances into super channels. These super channels are obtained by projecting the radiance spectra onto a set of pre-calculated eigenvectors. The forward model provides both super channel properties and jacobian in EOF space directly. For ultra-spectral sensors such as Infrared Atmospheric Sounding Interferometer (IASI) and the NPOESS Airborne Sounder Testbed Interferometer (NAST), a compression ratio of more than 80 can be achieved, leading to a significant reduction in computations involved in an inversion process. Results will be shown applying the algorithm to real IASI and NAST data.

Liu, Xu↗
NASA NTRS2007

Reconfigurable Control with Neural Network Augmentation for a Modified F-15 Aircraft

This paper describes the performance of a simplified dynamic inversion controller with neural network supplementation. This 6 DOF (Degree-of-Freedom) simulation study focuses on the results with and without adaptation of neural networks using a simulation of the NASA modified F-15 which has canards. One area of interest is the performance of a simulated surface failure while attempting to minimize the inertial cross coupling effect of a [B] matrix failure (a control derivative anomaly associated with a jammed or missing control surface). Another area of interest and presented is simulated aerodynamic failures ([A] matrix) such as a canard failure. The controller uses explicit models to produce desired angular rate commands. The dynamic inversion calculates the necessary surface commands to achieve the desired rates. The simplified dynamic inversion uses approximate short period and roll axis dynamics. Initial results indicated that the transient response for a [B] matrix failure using a Neural Network (NN) improved the control behavior when compared to not using a neural network for a given failure, However, further evaluation of the controller was comparable, with objections io the cross coupling effects (after changes were made to the controller). This paper describes the methods employed to reduce the cross coupling effect and maintain adequate tracking errors. The IA] matrix failure results show that control of the aircraft without adaptation is more difficult [leas damped) than with active neural networks, Simulation results show Neural Network augmentation of the controller improves performance in terms of backing error and cross coupling reduction and improved performance with aerodynamic-type failures.

Burken, John J.↗
NASA NTRS2008

Modeling-Error-Driven Performance-Seeking Direct Adaptive Control

This paper presents a stable discrete-time adaptive law that targets modeling errors in a direct adaptive control framework. The update law was developed in our previous work for the adaptive disturbance rejection application. The approach is based on the philosophy that without modeling errors, the original control design has been tuned to achieve the desired performance. The adaptive control should, therefore, work towards getting this performance even in the face of modeling uncertainties/errors. In this work, the baseline controller uses dynamic inversion with proportional-integral augmentation. Dynamic inversion is carried out using the assumed system model. On-line adaptation of this control law is achieved by providing a parameterized augmentation signal to the dynamic inversion block. The parameters of this augmentation signal are updated to achieve the nominal desired error dynamics. Contrary to the typical Lyapunov-based adaptive approaches that guarantee only stability, the current approach investigates conditions for stability as well as performance. A high-fidelity F-15 model is used to illustrate the overall approach.

Kulkarni, Nilesh V.↗
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 NTRS2010

Helicity and Filament Channels? The Straight Twist!

One of the most important and most puzzling features of the coronal magnetic field is that it appears to have smooth magnetic structure with little evidence for non-potentiality except at special locations, photospheric polarity inversions lines where the non-potentiality is observed as a filament channel. This characteristic feature of the closed-field corona is highly unexpected given that photospheric motions continuously tangle its magnetic field. Although reconnection can eliminate some of the injected structure, it cannot destroy the helicity, which should build up to produce observable complexity. We propose that an inverse cascade process transports the injected helicity from the interior of closed flux regions to their boundaries, polarity inversion lines, creating filament channels. We describe how the helicity is injected and transported and calculate the relevant rates. We argue that one process, helicity transport, can explain both the observed lack and presence of structure in the coronal magnetic field.

Antiochos, Spiro K.↗
NASA NTRS2010

The Origins of Magnetic Structure in the Corona and Wind

One of the most important and most puzzling features of the coronal magnetic field is that it appears to have smooth magnetic structure with little evidence for non-potentiality except at two special locations: photospheric polarity inversions lines. (non-potentiality observed as a filament channel) and coronal hole boundaries, (observed as the slow solar wind). This characteristic feature of the closed-field corona is highly unexpected given that its magnetic field is continuously tangled by photospheric motions. Although reconnection can eliminate some of the injected structure, it cannot destroy the helicity, which should build up to produce observable complexity. I propose that an inverse cascade process transports the injected helicity from the interior of closed flux regions to their boundaries inversion lines and coronal holes, creating both filament channels and the slow wind. We describe how the helicity is injected and transported and calculate the relevant rates. I argue that one process, helicity transport, can explain both the observed lack and presence of structure in the coronal magnetic field. This work has been supported by the NASA HTP, SR&T, and LWS programs.

Antiochos, Spiro K.↗
NASA NTRS2010

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗
NASA NTRS2005

Mapping From an Instrumented Glove to a Robot Hand

An algorithm has been developed to solve the problem of mapping from (1) a glove instrumented with joint-angle sensors to (2) an anthropomorphic robot hand. Such a mapping is needed to generate control signals to make the robot hand mimic the configuration of the hand of a human attempting to control the robot. The mapping problem is complicated by uncertainties in sensor locations caused by variations in sizes and shapes of hands and variations in the fit of the glove. The present mapping algorithm is robust in the face of these uncertainties, largely because it includes a calibration sub-algorithm that inherently adapts the mapping to the specific hand and glove, without need for measuring the hand and without regard for goodness of fit. The algorithm utilizes a forward-kinematics model of the glove derived from documentation provided by the manufacturer of the glove. In this case, forward-kinematics model signifies a mathematical model of the glove fingertip positions as functions of the sensor readings. More specifically, given the sensor readings, the forward-kinematics model calculates the glove fingertip positions in a Cartesian reference frame nominally attached to the palm. The algorithm also utilizes an inverse-kinematics model of the robot hand. In this case, inverse-kinematics model signifies a mathematical model of the robot finger-joint angles as functions of the robot fingertip positions. Again, more specifically, the inverse-kinematics model calculates the finger-joint commands needed to place the fingertips at specified positions in a Cartesian reference frame that is attached to the palm of the robot hand and that nominally corresponds to the Cartesian reference frame attached to the palm of the glove. Initially, because of the aforementioned uncertainties, the glove fingertip positions calculated by the forwardkinematics model in the glove Cartesian reference frame cannot be expected to match the robot fingertip positions in the robot-hand Cartesian reference frame. A calibration must be performed to make the glove and robot-hand fingertip positions correspond more precisely. The calibration procedure involves a few simple hand poses designed to provide well-defined fingertip positions. One of the poses is a fist. In each of the other poses, a finger touches the thumb. The calibration subalgorithm uses the sensor readings from these poses to modify the kinematical models to make the two sets of fingertip positions agree more closely.

Goza, Michael↗
NASA NTRS2011

Linear Approximation to Optimal Control Allocation for Rocket Nozzles with Elliptical Constraints

In this paper we present a straightforward technique for assessing and realizing the maximum control moment effectiveness for a launch vehicle with multiple constrained rocket nozzles, where elliptical deflection limits in gimbal axes are expressed as an ensemble of independent quadratic constraints. A direct method of determining an approximating ellipsoid that inscribes the set of attainable angular accelerations is derived. In the case of a parameterized linear generalized inverse, the geometry of the attainable set is computationally expensive to obtain but can be approximated to a high degree of accuracy with the proposed method. A linear inverse can then be optimized to maximize the volume of the true attainable set by maximizing the volume of the approximating ellipsoid. The use of a linear inverse does not preclude the use of linear methods for stability analysis and control design, preferred in practice for assessing the stability characteristics of the inertial and servoelastic coupling appearing in large boosters. The present techniques are demonstrated via application to the control allocation scheme for a concept heavy-lift launch vehicle.

Orr, Jeb S.↗
NASA NTRS2010

Simultaneous Retrieval of Aerosol and Surface Optical Properties from Combined Airborne- and Ground-Based Direct and Diffuse Radiometric Measurements

This paper presents a new method for simultaneously retrieving aerosol and surface reflectance properties from combined airborne and ground-based direct and diffuse radiometric measurements. The method is based on the standard Aerosol Robotic Network (AERONET) method for retrieving aerosol size distribution, complex index of refraction, and single scattering albedo, but modified to retrieve aerosol properties in two layers, below and above the aircraft, and parameters on surface optical properties from combined datasets (Cloud Absorption Radiometer (CAR) and AERONET data). A key advantage of this method is the inversion of all available spectral and angular data at the same time, while accounting for the influence of noise in the inversion procedure using statistical optimization. The wide spectral (0.34-2.30 m) and angular range (180 ) of the CAR instrument, combined with observations from an AERONET sunphotometer, provide sufficient measurement constraints for characterizing aerosol and surface properties with minimal assumptions. The robustness of the method was tested on observations made during four different field campaigns: (a) the Southern African Regional Science Initiative 2000 over Mongu, Zambia, (b) the Intercontinental Transport Experiment-Phase B over Mexico City, Mexico (c) Cloud and Land Surface Interaction Campaign over the Atmospheric Radiation Measurement (ARM) Central Facility, Oklahoma, USA, and (d) the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) over Elson Lagoon in Barrow, Alaska, USA. The four areas are dominated by different surface characteristics and aerosol types, and therefore provide good test cases for the new inversion method.

Gatebe, C. K.↗
NASA NTRS2011

Proximity Effects and Nonequilibrium Superconductivity in Transition-Edge Sensors

We have recently shown that normal-metal/superconductor (N/S) bilayer TESs (superconducting Transition-Edge Sensors) exhibit weak-link behavior.l Here we extend our understanding to include TESs with added noise-mitigating normal-metal structures (N structures). We find TESs with added Au structures also exhibit weak-link behavior as evidenced by exponential temperature dependence of the critical current and Josephson-like oscillations of the critical current with applied magnetic field. We explain our results in terms of an effect converse to the longitudinal proximity effect (LoPE) 1, the lateral inverse proximity effect (LaiPE), for which the order parameter in the N/S bilayer is reduced due to the neighboring N structures. Resistance and critical current measurements are presented as a function of temperature and magnetic field taken on square Mol Au bilayer TESs with lengths ranging from 8 to 130 {\mu}m with and without added N structures. We observe the inverse proximity effect on the bilayer over in-plane distances many tens of microns and find the transition shifts to lower temperatures scale approximately as the inverse square of the in- plane N-structure separation distance, without appreciable broadening of the transition width. We also present evidence for nonequilbrium superconductivity and estimate a quasiparticle lifetime of 1.8 \times 10-10 s for the bilayer. The LoPE model is also used to explain the increased conductivity at temperatures above the bilayer's steep resistive transition.

Sadleir, John E.↗
NASA NTRS2011

Applications of the BIOPHYS Algorithm for Physically-Based Retrieval of Biophysical, Structural and Forest Disturbance Information

Canopy reflectance model inversion using look-up table approaches provides powerful and flexible options for deriving improved forest biophysical structural information (BSI) compared with traditional statistical empirical methods. The BIOPHYS algorithm is an improved, physically-based inversion approach for deriving BSI for independent use and validation and for monitoring, inventory and quantifying forest disturbance as well as input to ecosystem, climate and carbon models. Based on the multiple-forward mode (MFM) inversion approach, BIOPHYS results were summarized from different studies (Minnesota/NASA COVER; Virginia/LEDAPS; Saskatchewan/BOREAS), sensors (airborne MMR; Landsat; MODIS) and models (GeoSail; GOMS). Applications output included forest density, height, crown dimension, branch and green leaf area, canopy cover, disturbance estimates based on multi-temporal chronosequences, and structural change following recovery from forest fires over the last century. Good correspondences with validation field data were obtained. Integrated analyses of multiple solar and view angle imagery further improved retrievals compared with single pass data. Quantifying ecosystem dynamics such as the area and percent of forest disturbance, early regrowth and succession provide essential inputs to process-driven models of carbon flux. BIOPHYS is well suited for large-area, multi-temporal applications involving multiple image sets and mosaics for assessing vegetation disturbance and quantifying biophysical structural dynamics and change. It is also suitable for integration with forest inventory, monitoring, updating, and other programs.

Peddle, Derek R.↗
NASA NTRS2012

Interpretation of the Total Magnetic Field Anomalies Measured by the CHAMP Satellite Over a Part of Europe and the Pannonian Basin

In this study we interpret the magnetic anomalies at satellite altitude over a part of Europe and the Pannonian Basin. These anomalies are derived from the total magnetic measurements from the CHAMP satellite. The anomalies reduced to an elevation of 324 km. An inversion method is used to interpret the total magnetic anomalies over the Pannonian Basin. A three dimensional triangular model is used in the inversion. Two parameter distributions: Laplacian and Gaussian are investigated. The regularized inversion is numerically calculated with the Simplex and Simulated Annealing methods and the anomalous source is located in the upper crust. A probable source of the magnetization is due to the exsolution of the hematite-ilmenite minerals.

Kis, K. I.↗
NASA NTRS2012

Forest Canopy Cover and Height from MISR in Topographically Complex Southwestern US Landscape Assessed with High Quality Reference Data

This study addresses the retrieval of spatially contiguous canopy cover and height estimates in southwestern USforests via inversion of a geometric-optical (GO) model against surface bidirectional reflectance factor (BRF) estimates from the Multi-angle Imaging SpectroRadiometer (MISR). Model inversion can provide such maps if good estimates of the background bidirectional reflectance distribution function (BRDF) are available. The study area is in the Sierra National Forest in the Sierra Nevada of California. Tree number density, mean crown radius, and fractional cover reference estimates were obtained via analysis of QuickBird 0.6 m spatial resolution panchromatic imagery usingthe CANopy Analysis with Panchromatic Imagery (CANAPI) algorithm, while RH50, RH75 and RH100 (50, 75, and 100 energy return) height data were obtained from the NASA Laser Vegetation Imaging Sensor (LVIS), a full waveform light detection and ranging (lidar) instrument. These canopy parameters were used to drive a modified version of the simple GO model (SGM), accurately reproducing patterns ofMISR 672 nm band surface reflectance (mean RMSE 0.011, mean R2 0.82, N 1048). Cover and height maps were obtained through model inversion against MISR 672 nm reflectance estimates on a 250 m grid.The free parameters were tree number density and mean crown radius. RMSE values with respect to reference data for the cover and height retrievals were 0.05 and 6.65 m, respectively, with of 0.54 and 0.49. MISR can thus provide maps of forest cover and height in areas of topographic variation although refinements are required to improve retrieval precision.

Chopping, Mark↗
NASA NTRS2012

CO2 Flux Estimation Errors Associated with Moist Atmospheric Processes

Vertical transport by moist sub-grid scale processes such as deep convection is a well-known source of uncertainty in CO2 source/sink inversion. However, a dynamical link between vertical transport, satellite based retrievals of column mole fractions of CO2, and source/sink inversion has not yet been established. By using the same offline transport model with meteorological fields from slightly different data assimilation systems, we examine sensitivity of frontal CO2 transport and retrieved fluxes to different parameterizations of sub-grid vertical transport. We find that frontal transport feeds off background vertical CO2 gradients, which are modulated by sub-grid vertical transport. The implication for source/sink estimation is two-fold. First, CO2 variations contained in moist poleward moving air masses are systematically different from variations in dry equatorward moving air. Moist poleward transport is hidden from orbital sensors on satellites, causing a sampling bias, which leads directly to small but systematic flux retrieval errors in northern mid-latitudes. Second, differences in the representation of moist sub-grid vertical transport in GEOS-4 and GEOS-5 meteorological fields cause differences in vertical gradients of CO2, which leads to systematic differences in moist poleward and dry equatorward CO2 transport and therefore the fraction of CO2 variations hidden in moist air from satellites. As a result, sampling biases are amplified and regional scale flux errors enhanced, most notably in Europe (0.43+/-0.35 PgC /yr). These results, cast from the perspective of moist frontal transport processes, support previous arguments that the vertical gradient of CO2 is a major source of uncertainty in source/sink inversion.

moist atmospheric processes↗
NASA NTRS2013

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 increase 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 therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) 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 rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗