Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion methods”

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,063 records · Page 59

Low-Cost Evaluation of EO-1 Hyperion and ALI for Detection and Biophysical Characterization of Forest Logging in Amazonia (NCC5-481)

Major uncertainties exist regarding the rate and intensity of logging in tropical forests worldwide: these uncertainties severely limit economic, ecological, and biogeochemical analyses of these regions. Recent sawmill surveys in the Amazon region of Brazil show that the area logged is nearly equal to total area deforested annually, but conversion of survey data to forest area, forest structural damage, and biomass estimates requires multiple assumptions about logging practices. Remote sensing could provide an independent means to monitor logging activity and to estimate the biophysical consequences of this land use. Previous studies have demonstrated that the detection of logging in Amazon forests is difficult and no studies have developed either the quantitative physical basis or remote sensing approaches needed to estimate the effects of various logging regimes on forest structure. A major reason for these limitations has been a lack of sufficient, well-calibrated optical satellite data, which in turn, has impeded the development and use of physically-based, quantitative approaches for detection and structural characterization of forest logging regimes. We propose to use data from the EO-1 Hyperion imaging spectrometer to greatly increase our ability to estimate the presence and structural attributes of selective logging in the Amazon Basin. Our approach is based on four "biogeophysical indicators" not yet derived simultaneously from any satellite sensor: 1) green canopy leaf area index; 2) degree of shadowing; 3) presence of exposed soil and; 4) non-photosynthetic vegetation material. Airborne, field and modeling studies have shown that the optical reflectance continuum (400-2500 nm) contains sufficient information to derive estimates of each of these indicators. Our ongoing studies in the eastern Amazon basin also suggest that these four indicators are sensitive to logging intensity. Satellite-based estimates of these indicators should provide a means to quantify both the presence and degree of structural disturbance caused by various logging regimes. Our quantitative assessment of Hyperion hyperspectral and ALI multi-spectral data for the detection and structural characterization of selective logging in Amazonia will benefit from data collected through an ongoing project run by the Tropical Forest Foundation, within which we have developed a study of the canopy and landscape biophysics of conventional and reduced-impact logging. We will add to our base of forest structural information in concert with an EO-1 overpass. Using a photon transport model inversion technique that accounts for non-linear mixing of the four biogeophysical indicators, we will estimate these parameters across a gradient of selective logging intensity provided by conventional and reduced impact logging sites. We will also compare our physical ly-based approach to both conventional (e.g., NDVI) and novel (e.g., SWIR-channel) vegetation indices as well as to linear mixture modeling methods. We will cross-compare these approaches using Hyperion and ALI imagers to determine the strengths and limitations of these two sensors for applications of forest biophysics. This effort will yield the first physical ly-based, quantitative analysis of the detection and intensity of selective logging in Amazonia, comparing hyperspectral and improved multi-spectral approaches as well as inverse modeling, linear mixture modeling, and vegetation index techniques.

Asner, Gregory P.↗

Combined linear theory/impact theory method for analysis and design of high speed configurations

Pressure distributions on a wing body at Mach 4.63 are calculated. The combined theory is shown to give improved predictions over either linear theory or impact theory alone. The combined theory is also applied in the inverse design mode to calculate optimum camber slopes at Mach 4.63. Comparisons with optimum camber slopes obtained from unmodified linear theory show large differences. Analysis of the results indicate that the combined theory correctly predicts the effect of thickness on the loading distributions at high Mach numbers, and that finite thickness wings optimized at high Mach numbers using unmodified linear theory will not achieve the minimum drag characteristics for which they are designed.

Brooke, D.↗

Irrigation Requirement Estimation using MODIS Vegetation Indices and Inverse Biophysical Modeling; A Case Study for Oran, Algeria

Human demand for food influences the water cycle through diversion and extraction of fresh water needed to support agriculture. Future population growth and economic development alone will substantially increase water demand and much of it for agricultural uses. For many semi-arid lands, socio-economic shifts are likely to exacerbate changes in climate as a driver of future water supply and demand. For these areas in particular, where the balance between water supply and demand is fragile, variations in regional climate can have potentially predictable effect on agricultural production. Satellite data and biophysically-based models provide a powerful method to quantify the interactions between local climate, plant growth and water resource requirements. In irrigated agricultural lands, satellite observations indicate high vegetation density while the precipitation amount indicates otherwise. This inconsistency between the observed precipitation and the observed canopy leaf density triggers the possibility that the observed high leaf density is due to an alternate source of water, irrigation. We explore an inverse process approach using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS), climatological data, and the NASA's Simple Biosphere model, SiB2, to quantitatively assess water demand in a semi-arid agricultural land by constraining the carbon and water cycles modeled under both equilibrium (balance between vegetation and prevailing local climate) and nonequilibrium (water added through irrigation) conditions. We postulate that the degree to which irrigated lands vary from equilibrium conditions is related to the amount of irrigation water used. We added water using two distribution methods: The first method adds water on top of the canopy and is a proxy for the traditional spray irrigation. The second method allows water to be applied directly into the soil layer and serves as proxy for drip irrigation. Our approach indicates that over the study site, for the month of July, spray irrigation resulted in an irrigation amount of about 1.4 mm per occurrence with an average frequency of occurrence of 24.6 hours. The simulated total monthly irrigation for July was 34.85 mm. In contrast, the drip irrigation resulted in less frequent irrigation events with an average water requirement about 57% less than that simulated during the spray irrigation case. The efficiency of the drip irrigation method rests on its reduction of the canopy interception loss compared to the spray irrigation method. When compared to a country-wide average estimate of irrigation water use, our numbers are quite low. We would have to revise the reported country level estimates downward to 17% or less

Bounoua, L.↗

Diagonal dominance for the multivariable Nyquist array using function minimization

A new technique for the design of multivariable control systems using the multivariable Nyquist array method was developed. A conjugate direction function minimization algorithm is utilized to achieve a diagonal dominant condition over the extended frequency range of the control system. The minimization is performed on the ratio of the moduli of the off-diagonal terms to the moduli of the diagonal terms of either the inverse or direct open loop transfer function matrix. Several new feedback design concepts were also developed, including: (1) dominance control parameters for each control loop; (2) compensator normalization to evaluate open loop conditions for alternative design configurations; and (3) an interaction index to determine the degree and type of system interaction when all feedback loops are closed simultaneously. This new design capability was implemented on an IBM 360/75 in a batch mode but can be easily adapted to an interactive computer facility. The method was applied to the Pratt and Whitney F100 turbofan engine.

Leininger, G. G.↗

Hot-gas cold-dust pumping for water masers associated with H II regions

A collisional pump with an internal sink is proposed for the water masers associated with H II regions, where the population inversion occurs due to the absorption by cold ice-mantle grains in a highly dusty cloud of the far-infrared line radiation of hot water vapor. A new escape probability method is developed to calculate the transfer of line radiation in dusty medium. The pump mechanism explains the power of usual maser sources associated with H II regions and the enormous power of the sources associated with W49 N and external galaxies. Models of maser clouds have a radius of 5 x 10 to the 15th-10 to the 16th cm, an H2 number density of 4 x 10 to the 9th/cu cm, an expansion velocity of 10-30 km/s, a kinetic temperature of 350 K, and a grain temperature of 100 K. Giant maser sources require grains of the size about 1 micron. The apparent size of the emission spots (approximately 10 to the 13th cm) observed by VLBI is interpreted as due to a fluctuation in the cloud, and the assembly of the spots is spread within a size of 10 to the 16th cm. The temperature difference between the dust and gas is due to a relaxation process after an infrared burst accompanying protostar formation.

Deguchi, S.↗

A Strassen-Newton algorithm for high-speed parallelizable matrix inversion

Techniques are described for computing matrix inverses by algorithms that are highly suited to massively parallel computation. The techniques are based on an algorithm suggested by Strassen (1969). Variations of this scheme use matrix Newton iterations and other methods to improve the numerical stability while at the same time preserving a very high level of parallelism. One-processor Cray-2 implementations of these schemes range from one that is up to 55 percent faster than a conventional library routine to one that is slower than a library routine but achieves excellent numerical stability. The problem of computing the solution to a single set of linear equations is discussed, and it is shown that this problem can also be solved efficiently using these techniques.

Bailey, David H.↗

Active control of sound radiation from a vibrating rectangular panel by sound sources and vibration inputs - An experimental comparison

Active control of sound radiation from a rectangular panel by two different methods has been experimentally studied and compared. In the first method a single control force applied directly to the structure is used with a single error microphone located in the radiated acoustic field. Global attenuation of radiated sound was observed to occur by two main mechanisms. For 'on-resonance' excitation, the control force had the effect of increasing the total panel input impedance presented to the nosie source, thus reducing all radiated sound. For 'off-resonance' excitation, the control force tends not significantly to modify the panel total response amplitude but rather to restructure the relative phases of the modes leading to a more complex vibration pattern and a decrease in radiation efficiency. For acoustic control, the second method, the number of acoustic sources required for global reduction was seen to increase with panel modal order. The mechanism in this case was that the acoustic sources tended to create an inverse pressure distribution at the panel surface and thus 'unload' the panel by reducing the panel radiation impedance. In general, control by structural inputs appears more effective than control by acoustic sources for structurally radiated noise.

Fuller, C. R.↗

Solar Confocal Interferometers for Sub-Picometer-Resolution Spectral Filters

The confocal Fabry-Perot interferometer allows sub-picometer spectral resolution of Fraunhofer line profiles. Such high spectral resolution is needed to keep pace with the higher spatial resolution of the new set of large-aperture solar telescopes. The line-of-sight spatial resolution derived for line profile inversions would then track the improvements of the transverse spatial scale provided by the larger apertures. The confocal interferometer's unique properties allow a simultaneous increase in both etendue and spectral power. Methods: We have constructed and tested two confocal interferometers. Conclusions: In this paper we compare the confocal interferometer with other spectral imaging filters, provide initial design parameters, show construction details for two designs, and report on the laboratory test results for these interferometers, and propose a multiple etalon system for future testing of these units and to obtain sub-picometer spectral resolution information on the photosphere in both the visible and near-infrared.

Gary, G. Allen↗

A Method for Localizing Energy Dissipation in Blazars Using Fermi Variability

The distance of the Fermi-detected blazar gamma-ray emission site from the supermassive black hole is a matter of active debate. Here we present a method for testing if the GeV emission of powerful blazars is produced within the sub-pc scale broad line region (BLR) or farther out in the pc-scale molecular torus (MT) environment. If the GeV emission takes place within the BLR, the inverse Compton (IC) scattering of the BLR ultraviolet (UV) seed photons that produces the gamma-rays takes place at the onset of the Klein-Nishina regime. This causes the electron cooling time to become practically energy independent and the variation of the gamma-ray emission to be almost achromatic. If on the other hand the -ray emission is produced farther out in the pc-scale MT, the IC scattering of the infrared (IR) MT seed photons that produces the gamma-rays takes place in the Thomson regime, resulting to energy-dependent electron cooling times, manifested as faster cooling times for higher Fermi energies. We demonstrate these characteristics and discuss the applicability and limitations of our method.

Dotson, Amanda↗

The inference of atmospheric ozone using satellite nadir measurements in the 1042/cm band

A description and detailed analysis of a technique for inferring atmospheric ozone information from satellite nadir measurements in the 1042 cm band are presented. A method is formulated for computing the emission from the lower boundary under the satellite which circumvents the difficult analytical problems caused by the presence of atmospheric clouds and the watervapor continuum absorption. The inversion equations are expanded in terms of the eigenvectors and eigenvalues of a least-squares-solution matrix, and an analysis is performed to determine the information content of the radiance measurements. Under favorable conditions there are only two pieces of independent information available from the measurements: (1) the total ozone and (2) the altitude of the primary maximum in the ozone profile.

Russell, J. M., III↗

An efficient numerical technique for calculating thermal spreading resistance

An efficient numerical technique for solving the equations resulting from finite difference analyses of fields governed by Poisson's equation is presented. The method is direct (noniterative)and the computer work required varies with the square of the order of the coefficient matrix. The computational work required varies with the cube of this order for standard inversion techniques, e.g., Gaussian elimination, Jordan, Doolittle, etc.

Gale, E. H., Jr.↗

MAGSAT science investigations

A preliminary magnetization model for the U.S. was produced from MAGSAT data. The double grid processing method resulted in a 100 km source spacing. While some spurious features are present in the map, much of the detail appears to be real. The determination of the optional source spacing for inversion of MAGSAT delta B data is in progress. A fine attitude vector data set was obtained for comparison with the theoretical vector field associated with the magnetization model.

Source record↗

Computational methods for estimation of parameters in hyperbolic systems

Approximation techniques for estimating spatially varying coefficients and unknown boundary parameters in second order hyperbolic systems are discussed. Methods for state approximation (cubic splines, tau-Legendre) and approximation of function space parameters (interpolatory splines) are outlined and numerical findings for use of the resulting schemes in model "one dimensional seismic inversion' problems are summarized.

Banks, H. T.↗

Group-kinetic theory and modeling of atmospheric turbulence

A group kinetic method is developed for analyzing eddy transport properties and relaxation to equilibrium. The purpose is to derive the spectral structure of turbulence in incompressible and compressible media. Of particular interest are: direct and inverse cascade, boundary layer turbulence, Rossby wave turbulence, two phase turbulence; compressible turbulence, and soliton turbulence. Soliton turbulence can be found in large scale turbulence, turbulence connected with surface gravity waves and nonlinear propagation of acoustical and optical waves. By letting the pressure gradient represent the elementary interaction among fluid elements and by raising the Navier-Stokes equation to higher dimensionality, the master equation was obtained for the description of the microdynamical state of turbulence.

Tchen, C. M.↗

A feasibility study for long-path multiple detection using a neural network

Least-squares inverse filters have found widespread use in the deconvolution of seismograms and the removal of multiples. The use of least-squares prediction filters with prediction distances greater than unity leads to the method of predictive deconvolution which can be used for the removal of long path multiples. The predictive technique allows one to control the length of the desired output wavelet by control of the predictive distance, and hence to specify the desired degree of resolution. Events which are periodic within given repetition ranges can be attenuated selectively. The method is thus effective in the suppression of rather complex reverberation patterns. A back propagation(BP) neural network is constructed to perform the detection of first arrivals of the multiples and therefore aid in the more accurate determination of the predictive distance of the multiples. The neural detector is applied to synthetic reflection coefficients and synthetic seismic traces. The processing results show that the neural detector is accurate and should lead to an automated fast method for determining predictive distances across vast amounts of data such as seismic field records. The neural network system used in this study was the NASA Software Technology Branch's NETS system.

Feuerbacher, G. A.↗

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↗

On-line parameter estimation using a high sensitivity estimator

An on-line parameter identification method is presented. The method is based on a recursive formulation of the maximum likelihood method, with a significant modification on the gains of the state estimator. In the conventional maximum likelihood method, the Kalman gains are used in the state estimator. This produces unbiased, minimum variance parameter estimates in the presence of process noise and measurement noise, but it also slows the convergence rate when the identification is done on-line. Here we suggest choosing the gains to maximize a measure of the sensitivities of the state estimates to parameter variations. One such criterion is to minimize the trace of the inverse information matrix. This increases the convergence rate significantly. After one or two time constants, the gains are switched to the Kalman values to assure unbiased, minimum-variance estimates. The state estimate will initially be nonoptimal, and may not be adequate for control purposes. In this case, a parallel Kalman filter which uses the identifier's parameter estimates can be used. This method is applied here for the identification of a simple first-order system, and for the identification of short-period stability derivatives of an F-8 aircraft from simulated data.

Mishne, D.↗

Remote monitoring of environmental particulate pollution - A problem in inversion of first-kind integral equations

The determination of the microstructure, chemical nature, and dynamical evolution of scattering particulates in the atmosphere is considered. A description is given of indirect sampling techniques which can circumvent most of the difficulties associated with direct sampling techniques, taking into account methods based on scattering, extinction, and diffraction of an incident light beam. Approaches for reconstructing the particulate size distribution from the direct and the scattered radiation are discussed. A new method is proposed for determining the chemical composition of the particulates and attention is given to the relevance of methods of solution involving first kind Fredholm integral equations.

Fymat, A. L.↗