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 793 records · Page 44

Limb-brightening observations from the OSO-7 satellite. II - Comparison of Abel-inverted intensities of Fe XIV and Fe XIII EUV emission lines with predictions

Intensities of Fe XIV and Fe XIII EUV emission lines obtained at coronal locations beyond the limb by the Goddard spectroheliograph on the OSO 7 satellite have been corrected for the wavelength dependence of the instrument's sensitivity and have been Abel-inverted to provide a valid comparison with theoretical predictions for each ion. Details of the Abel-inversion procedure are given, including explicit formulas for application of Bracewell's (1956) method. The intensity ratios of pairs of lines originating from a common level are compared with expected theoretical transition probability ratios over a range of heliocentric distance; deviations in some cases yield information about adjacent unclassified lines. Comparison of the observations with predictions for Fe XIV and Fe XIII shows generally good agreement, with a few interesting discrepancies that may imply a corresponding need for more accurate collisional excitation cross sections. The same comparison yields the variation of electron density with heliocentric radius for each ion separately; the two density functions are found to agree within a factor of three.

Kastner, S. O.↗

Numerical analysis of free vibrations of damped rotating structures

This paper is concerned with the efficient numerical solution of damped and undamped free vibration problems of rotating structures. While structural discretization is achieved by the finite element method, the associated eigenproblem solution is effected by a combined Sturm sequence and inverse iteration technique that enables the computation of a few required roots only without having to compute any other. For structures of complex configurations, a modal synthesis technique is also presented, which is based on appropriate combinations of eigenproblem solution of various structural components. Such numerical procedures are general in nature, which fully exploit matrix sparsity inherent in finite element discretizations, and prove to be most efficient for the vibration analysis of any damped rotating structure, such as rotating machineries, helicopter and turbine blades, spinning space stations, among others.

Gupta, K. K.↗

Stochastic inverse problem in the radiation of noise

The reported investigation is concerned with a stochastic inverse radiation problem in a uniform medium. The problem is illustrated with the aid of a simple model consisting of an array of point sources. The entropy functional is chosen to be the structural functional in determining the source distribution. A general theory for the stochastic inverse problem is introduced. It is shown that the general procedure yields the methods of the Lagrangian multiplier, when the structural and residual functionals are specialized. Tihonov's regularization and a method related to generalized or pseudoinverses are also obtained. Examples considered for purposes of illustration are related to a continuous source with the least noise intensity, a continuous source with a potential, and an axisymmetric line source.

Chow, P. L.↗

Computation of nonequilibrium radiating shock layers

A computational technique of coupling radiative transfer to fluid motion is developed for axisymmetric blunt body shock layer flows in a thermochemical nonequilibrium environment. The coupled formulation of radiation and flowfield leads to a governing set of integro-differential equations. This equation set is solved using a modified Gauss-Seidel line relaxation techniques which incorporates the inversion of full block matrix associated with radiative transfer using a block iteration method. The thermodynamic state of the gas is described by three temperatures: translational, rotational, and vibrational-electronic. Radiation phenomenon is assumed to be governed by the vibrational-electronic temperature. The radiative properties are described by a spectrally detailed model. The computations are presented for two cases, including the Fire II flight experiment. It is shown that the method converges and the calculated spectra qualitatively agree with the experimental data for the two test cases. The calculated total radiative flux underestimates the measured values owing to the low vibrational-electronic temperature predicted in the flowfield calculation.

Gokcen, Tahir↗

Boundary layer stratus clouds - Inferred from satellite infrared spectral measurements over oceans

A method to remotely sense marine stratus clouds, that are underneath the boundary layer inversion produced by strong subsidence, is developed with the help of the measurements made by the Infrared Interferometer Spectrometer (IRIS) on board the Nimbus 4 satellite. These measurements, which had a spectral resolution of 2.8/cm and a field view of about 100 km in the region 870-980/cm, allow one to compare the weak absorption due to water vapor lines with that of the weak nu3-nu1 band of CO2. From such a comparison it is possible to sense the presence of these boundary layer stratus clouds. Maps of these clouds over subsidence regions of the tropical and subtropical oceans deduced from IRIS data show the usefulness of the technique to climatological studies.

Prabhakara, C.↗

A New Neural Network Approach Including First-Guess for Retrieval of Atmospheric Water Vapor, Cloud Liquid Water Path, Surface Temperature and Emissivities Over Land From Satellite Microwave Observations

The analysis of microwave observations over land to determine atmospheric and surface parameters is still limited due to the complexity of the inverse problem. Neural network techniques have already proved successful as the basis of efficient retrieval methods for non-linear cases, however, first-guess estimates, which are used in variational methods to avoid problems of solution non-uniqueness or other forms of solution irregularity, have up to now not been used with neural network methods. In this study, a neural network approach is developed that uses a first-guess. Conceptual bridges are established between the neural network and variational methods. The new neural method retrieves the surface skin temperature, the integrated water vapor content, the cloud liquid water path and the microwave surface emissivities between 19 and 85 GHz over land from SSM/I observations. The retrieval, in parallel, of all these quantities improves the results for consistency reasons. A data base to train the neural network is calculated with a radiative transfer model and a a global collection of coincident surface and atmospheric parameters extracted from the National Center for Environmental Prediction reanalysis, from the International Satellite Cloud Climatology Project data and from microwave emissivity atlases previously calculated. The results of the neural network inversion are very encouraging. The r.m.s. error of the surface temperature retrieval over the globe is 1.3 K in clear sky conditions and 1.6 K in cloudy scenes. Water vapor is retrieved with a r.m.s. error of 3.8 kg/sq m in clear conditions and 4.9 kg/sq m in cloudy situations. The r.m.s. error in cloud liquid water path is 0.08 kg/sq m . The surface emissivities are retrieved with an accuracy of better than 0.008 in clear conditions and 0.010 in cloudy conditions. Microwave land surface temperature retrieval presents a very attractive complement to the infrared estimates in cloudy areas: time record of land surface temperature will be produced.

Aires, F.↗

A Nonlinear Reduced Order Method for Prediction of Acoustic Fatigue

The goal of this investigation is to assess the quality of high-cycle-fatigue life estimation via a reduced order method, for structures undergoing geometrically nonlinear random vibrations. Modal reduction is performed with several different suites of basis functions. After numerically solving the reduced order system equations of motion, the physical displacement time history is obtained by an inverse transformation and stresses are recovered. Stress ranges obtained through the rainflow counting procedure are used in a linear damage accumulation method to yield fatigue estimates. Fatigue life estimates obtained using various basis functions in the reduced order method are compared with those obtained from numerical simulation in physical degrees-of-freedom.

Przekop, Adam↗

A New Method to Cross Calibrate and Validate TOMS, SBUV/2, and SCIAMACHY Measurements

A unique method to validate back scattered ultraviolet (buv) type satellite data that complements the measurements from existing ground networks is proposed. The method involves comparing the zenith sky radiance measurements from the ground to the nadir radiance measurements taken from space. Since the measurements are compared directly, the proposed method is superior to any other method that involves comparing derived products (for example, ozone), because comparison of derived products involve inversion algorithms which are susceptible to several type of errors. Forward radiative transfer (RT) calculations show that for an aerosol free atmosphere, the ground-based zenith sky radiance measurement and the satellite nadir radiance measurements can be predicted with an accuracy of better than 1 percent. The RT computations also show that for certain values of the solar zenith angles, the radiance comparisons could be better than half a percent. This accuracy is practically independent of ozone amount and aerosols in the atmosphere. Experiences with the Shuttle Solar Backscatter Ultraviolet (SSBUV) program show that the accuracy of the ground-based zenith sky radiance measuring instrument can be maintained at a level of a few tenth of a percent. This implies that the zenith sky radiance measurements can be used to validate Total Ozone Mapping Spectrometer (TOMS), Solar Backscatter Ultraviolet (SBUV/2), and The SCanning Imaging Absorption SpectroMeter for Atmospheric CHartographY (SCIAMACHY) radiance data. Also, this method will help improve the long term precision of the measurements for better trend detection and the accuracy of other BUV products such as tropospheric ozone and aerosols. Finally, in the long term, this method is a good candidate to inter-calibrate and validate long term observations of upcoming operational instruments such as Global Ozone Monitoring Experiment (GOME-2), Ozone Mapping Instrument (OMI), Ozone Dynamics Ultraviolet Spectrometer (ODUS), and Ozone Mapping and Profiler Suite (OMPS).

Ahmad, Ziauddin↗

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We describe a variation of the Optimal Estimation (OE) method for greenhouse gas remote sensing retrievals using a singular value decomposition (SVD) and an uninformative prior. The SVD method is capable of discerning vertical information in column integrated absorption measurements. While traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We discuss the method, illustrate its use on an integrated path differential absorption CO2 lidar measurement model, and compare it to traditional optimal estimation using numerical simulations. In the absence of forward model error, the mathematics behind the SVD method guarantee it to be bias-free, which is confirmed by the numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. While the SVD approach can be used for most trace gas retrievals, it is particularly useful for situations where prior knowledge of the trace gas profile is poor. The SVD analysis is also useful in illustrating how vertical information is treated by the traditional OE approach.

Ramanathan, Anand K.↗

Investigation Of Prandtl-Ishlinskii Hysteresis Compensation for Deep Space Optical Communications Pointing Control

This paper studies the simulated performance of the Deep Space Optical Communication (DSOC) laser downlink pointing control system that has a pointing requirement at the sub-microradian scale. To achieve the required pointing precision, the DSOC Flight Laser Transceiver (FLT) point-ahead mirror (PAM) employs piezoelectric actuators, which are widely used in precision applications due to their high motion resolution, rapid response, and high output forces. However, a piezoelectric actuator's performance is often limited by nonlinearities in their response, most significantly hysteresis. To study the effect of hysteresis on our simulated pointing performance, an invertible hysteresis model is developed using the Classical Prantdl-Ishlinskii (CPI) method. The simulation includes models of DSOC hardware and electronic system components based on key parameters (bias, frequency response, noise), with the hysteresis modeled in the PAM piezo actuator dynamics. A Proportional-Integral controller provides feedback with measurements from two potential sources: strain gauge sensors or a photon-counting camera. To compensate for hysteresis behavior, a feedforward inverse CPI hysteresis model is implemented in the control system. Two CPI models are built using different sampling methods (uniform sampling and cosine spacing) and controller performance is compared. For the maximum command range of ±400 μrad at 1 mHz, results show that the feedforward combined with feedback control reduces output tracking error by 50% or 58% (to 0.20 or 0.38 μrad) compared to feedback control alone (from 0.40 or 0.90 μrad), corresponding to the uniform-sampled and cosine-sampled CPI models, respectively.

Marinello, Lindsey A.↗

Parallel Preconditioning for CFD Problems on the CM-5

Up to today, preconditioning methods on massively parallel systems have faced a major difficulty. The most successful preconditioning methods in terms of accelerating the convergence of the iterative solver such as incomplete LU factorizations are notoriously difficult to implement on parallel machines for two reasons: (1) the actual computation of the preconditioner is not very floating-point intensive, but requires a large amount of unstructured communication, and (2) the application of the preconditioning matrix in the iteration phase (i.e. triangular solves) are difficult to parallelize because of the recursive nature of the computation. Here we present a new approach to preconditioning for very large, sparse, unsymmetric, linear systems, which avoids both difficulties. We explicitly compute an approximate inverse to our original matrix. This new preconditioning matrix can be applied most efficiently for iterative methods on massively parallel machines, since the preconditioning phase involves only a matrix-vector multiplication, with possibly a dense matrix. Furthermore the actual computation of the preconditioning matrix has natural parallelism. For a problem of size n, the preconditioning matrix can be computed by solving n independent small least squares problems. The algorithm and its implementation on the Connection Machine CM-5 are discussed in detail and supported by extensive timings obtained from real problem data.

Simon, Horst D.↗

Replicated Composite Optics Development

Advanced optical systems for applications such as grazing incidence Wolter I x-ray mirror assemblies require extraordinary mirror surfaces in ten-ns of fine surface finish and figure. The impeccable mirror surface is on the inside of the rotational mirror form. One practical method of producing devices with these requirements is to first fabricate an exterior surface for the optical device then replicate that surface to have the inverse component with lightweight characteristics. The replicate optic is not better than the master or mandrel from which it is made. This task is a continuance of previous studies to identify methods and materials for forming these extremely low roughness optical components.

Engelhaupt, Darell↗

Optimal inputs for system identification

Identification criteria are presented for linear dynamic systems with and without process noise. With process noise, the state equations are replaced by the Kalman filter equations. If the identification performance index is expanded in a Taylor's series with respect to the parameters to be identified, then maximizing the weighting factor of the quadratic term with respect to the inputs will insure that an identification algorithm will converge more rapidly and to a more accurate result than with nonoptimal inputs. The expectation of this weighting factor is the Fisher information matrix, and its inverse is a lower bound for the covariance of the parameters. Direct and indirect methods of calculating the information matrix are presented for systems with and without process noise.

Reid, D. B.↗

Rotational excitation of symmetric top molecules by collisions with atoms: Close coupling, coupled states, and effective potential calculations for NH3-He

The formalism for describing rotational excitation in collisions between symmetric top rigid rotors and spherical atoms is presented both within the accurate quantum close coupling framework and also the coupled states approximation of McGuire and Kouri and the effective potential approximation of Rabitz. Calculations are reported for thermal energy NH3-He collisions, treating NH3 as a rigid rotor and employing a uniform electron gas (Gordon-Kim) approximation for the intermolecular potential. Coupled states are found to be in nearly quantitative agreement with close coupling results while the effective potential method is found to be at least qualitatively correct. Modifications necessary to treat the inversion motion in NH3 are discussed.

Green, S.↗

Spline methods for approximating quantile functions and generating random samples

Two cubic spline formulations are presented for representing the quantile function (inverse cumulative distribution function) of a random sample of data. Both B-spline and rational spline approximations are compared with analytic representations of the quantile function. It is also shown how these representations can be used to generate random samples for use in simulation studies. Comparisons are made on samples generated from known distributions and a sample of experimental data. The spline representations are more accurate for multimodal and skewed samples and to require much less time to generate samples than the analytic representation.

Schiess, J. R.↗

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.↗