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 91 records · Page 5

Development of the Model of Galactic Interstellar Emission for Standard Point-Source Analysis of Fermi Large Area Telescope Data

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point-source and extended-source studies rely on the modeling of this diffuse emission for accurate characterization. Here, we describe the development of the Galactic Interstellar Emission Model (GIEM),which is the standard adopted by the LAT Collaboration and is publicly available. This model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse-Compton emission produced in the Galaxy. In the GIEM, we also include large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20deg and we observe an enhanced emission toward their base extending in the north and south Galactic directions and located within approximately 4deg of the Galactic Center.

gamma rays: diffuse background – gamma rays: gen↗

Global atomic hydrogen density derived from OGO-6 Lyman-alpha measurements

The paper analyzes a one-year set of Lyman-alpha airglow data measured in the local zenith at altitudes from 400 to 1100 km by a UV photometer aboard OGO-6. The zenith-intensity data are fitted to theoretical airglow calculations in four spherically symmetric models of the hydrogen geocorona to determine both the Ly-alpha solar flux at line center and the average atomic hydrogen column density. After correcting for a loss of instrument sensitivity, the Ly-alpha flux is found to be linearly correlated with daily Zurich sunspot number. It is also found that the hydrogen density is inversely correlated with Jacchia exospheric temperature, but the dependence is not that predicted by steady-state models with Jeans evaporative escape as the only loss mechanism. It is suggested that charge-exchange production of fast hydrogen atoms from 'hot' ionospheric protons might provide the additional loss this result requires.

Thomas, G. E.↗

Automatic Incorporation of Flaws into Cad and FEA Models Based On Measured NDE Data

Automatic incorporation of non-destructive evaluation (NDE) data into CAD and finite element (FE) models for structural analysis is essential for accurate residual-life and strength analyses of composite structures. Currently, two significant but independent challenges exist for incorporating NDE data in such analyses: interpretation of the NDE signals to locate and characterize the structural defects and the incorporation of that information into structural analysis tools. The first challenge will be addressed through linear inversion of calculated synthetic responses generated by a simulation engine to evaluate flaw size, location, and other characteristics. The second challenge will be addressed by building on prior and parallel defect modeling efforts to integrate suitable flaw models into the CAD/finite-element representation of a specimen. Ultrasonic NDE information is usually represented using C-Scans of time-of-flight or intensity information. Other modalities such as thermography and computed tomography (CT) could also be used to characterize the structure. This information will first be projected into material coordinates and then aligned with the specimen’s finite element mesh. The NDE data can then be incorporated by inserting flaw geometry into the CAD and FE analysis model and automatic modification of local material properties and boundary conditions. This process will enable the use of structural analysis tools even for practical modalities such as ultrasound and thermography which may not directly give accurate-enough geometric information to fully define the flaw. The full flaw characteristics will be inferred using prior expert and modeling knowledge to categorize flaw signatures. Outputs for this project will include demonstrated methods for creating structural analysis geometric models that incorporate flaw characterization information; these models can then be used directly in popular commercial FEA software. In addition, the primary deliverable for the project will include a software that will be able to automatically incorporate different kinds of flaws into an existing FE models. The models will be generated in standard file formats to enable direct import into the analysis software.

Adarsh Krishnamurthy↗

Fuzzy logic based robotic controller

Existing Proportional-Integral-Derivative (PID) robotic controllers rely on an inverse kinematic model to convert user-specified cartesian trajectory coordinates to joint variables. These joints experience friction, stiction, and gear backlash effects. Due to lack of proper linearization of these effects, modern control theory based on state space methods cannot provide adequate control for robotic systems. In the presence of loads, the dynamic behavior of robotic systems is complex and nonlinear, especially where mathematical modeling is evaluated for real-time operators. Fuzzy Logic Control is a fast emerging alternative to conventional control systems in situations where it may not be feasible to formulate an analytical model of the complex system. Fuzzy logic techniques track a user-defined trajectory without having the host computer to explicitly solve the nonlinear inverse kinematic equations. The goal is to provide a rule-based approach, which is closer to human reasoning. The approach used expresses end-point error, location of manipulator joints, and proximity to obstacles as fuzzy variables. The resulting decisions are based upon linguistic and non-numerical information. This paper presents a solution to the conventional robot controller which is independent of computationally intensive kinematic equations. Computer simulation results of this approach as obtained from software implementation are also discussed.

Attia, F.↗

Confidence set inference with a prior quadratic bound

In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z=(z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y (sup 0) = (y (sub 1) (sup 0), ..., y (sub D (sup 0)), using full or partial knowledge of the statistical distribution of the random errors in y (sup 0). The data space Y containing y(sup 0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x, Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. Confidence set inference is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface.

Backus, George E.↗

Uniform apparent contrast noise: A picture of the noise of the visual contrast detection system

A picture which is a sample of random contrast noise is generated. The noise amplitude spectrum in each region of the picture is inversely proportional to spatial frequency contrast sensitivity for that region, assuming the observer fixates the center of the picture and is the appropriate distance from it. In this case, the picture appears to have approximately the same contrast everywhere. To the extent that contrast detection thresholds are determined by visual system noise, this picture can be regarded as a picture of the noise of that system. There is evidence that, at different eccentricities, contrast sensitivity functions differ only by a magnification factor. The picture was generated by filtering a sample of white noise with a filter whose frequency response is inversely proportional to foveal contrast sensitivity. It was then stretched by a space-varying magnification function. The picture summmarizes a noise linear model of detection and discrimination of contrast signals by referring the model noise to the input picture domain.

Ahumada, A. J., Jr.↗

Surface and Atmospheric Parameter Retrieval From AVIRIS Data: The Importance of Non-Linear Effects

AVIRIS data represent a new and important approach for the retrieval of atmospheric and surface parameters from optical remote sensing data. Not only as a test for future space systems, but also as an operational airborne remote sensing system, the development of algorithms to retrieve information from AVIRIS data is an important step to these new approaches and capabilities. Many things have been learned since AVIRIS became operational, and the successive technical improvements in the hardware and the more sophisticated calibration techniques employed have increased the quality of the data to the point of almost meeting optimum user requirements. However, the potential capabilities of imaging spectrometry over the standard multispectral techniques have still not been fully demonstrated. Reasons for this are the technical difficulties in handling the data, the critical aspect of calibration for advanced retrieval methods, and the lack of proper models with which to invert the measured AVIRIS radiances in all the spectral channels. To achieve the potential of imaging spectrometry, these issues must be addressed. In this paper, an algorithm to retrieve information about both atmospheric and surface parameters from AVIRIS data, by using model inversion techniques, is described. Emphasis is put on the derivation of the model itself as well as proper inversion techniques, robust to noise in the data and an inadequate ability of the model to describe natural variability in the data. The problem of non-linear effects is addressed, as it has been demonstrated to be a major source of error in the numerical values retrieved by more simple, linear-based approaches. Non-linear effects are especially critical for the retrieval of surface parameters where both scattering and absorption effects are coupled, as well as in the cases of significant multiple-scattering contributions. However, sophisticated modeling approaches can handle such non-linear effects, which are especially important over vegetated surfaces. All the data used in this study were acquired during the 1991 Multisensor Airborne Campaign (MAC-Europe), as part of the European Field Experiment on a Desertification-threatened Area (EFEDA), carried out in Spain in June-July 1991.

Green Robert O.↗

Neural Flight Control System

The Neural Flight Control System (NFCS) was developed to address the need for control systems that can be produced and tested at lower cost, easily adapted to prototype vehicles and for flight systems that can accommodate damaged control surfaces or changes to aircraft stability and control characteristics resulting from failures or accidents. NFCS utilizes on a neural network-based flight control algorithm which automatically compensates for a broad spectrum of unanticipated damage or failures of an aircraft in flight. Pilot stick and rudder pedal inputs are fed into a reference model which produces pitch, roll and yaw rate commands. The reference model frequencies and gains can be set to provide handling quality characteristics suitable for the aircraft of interest. The rate commands are used in conjunction with estimates of the aircraft s stability and control (S&C) derivatives by a simplified Dynamic Inverse controller to produce virtual elevator, aileron and rudder commands. These virtual surface deflection commands are optimally distributed across the aircraft s available control surfaces using linear programming theory. Sensor data is compared with the reference model rate commands to produce an error signal. A Proportional/Integral (PI) error controller "winds up" on the error signal and adds an augmented command to the reference model output with the effect of zeroing the error signal. In order to provide more consistent handling qualities for the pilot, neural networks learn the behavior of the error controller and add in the augmented command before the integrator winds up. In the case of damage sufficient to affect the handling qualities of the aircraft, an Adaptive Critic is utilized to reduce the reference model frequencies and gains to stay within a flyable envelope of the aircraft.

Gundy-Burlet, Karen↗

Confidence set inference with a prior quadratic bound

In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z = (z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y(0)=(y sub 1(0),...,y sub D(0)) knowledge of the statistical distribution of the random errors in y(0). The data space Y containing y(0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x (e.g., energy or dissipation rate), Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. CSI is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface. Neither the heat flow nor the energy bound is strong enough to permit estimation of B(r) at single points on the CMB, but the heat flow bound permits estimation of uniform averages of B(r) over discs on the CMB, and both bounds permit weighted disc-averages with continous weighting kernels. Both bounds also permit estimation of low-degree Gauss coefficients at the CMB. The heat flow bound resolves them up to degree 8 if the crustal field at satellite altitudes must be treated as a systematic error, but can resolve to degree 11 under the most favorable statistical treatment of the crust. These two limits produce circles of confusion on the CMB with diameters of 25 deg and 19 deg respectively.

Backus, George E.↗

Coronal temperatures of selected active cool stars as derived from low resolution Einstein observations

Mean coronal temperatures of some active G-K stars were derived from Rev1-processed Einstein-observatory's IPC-spectra. The combined X-ray and transition region emission line data are in rough agreement with static coronal loop models. Although the sample is too small to derive any statistically significant conclusions, it suggests that the mean coronal temperature depends linearly on the inverse Rossby-number, with saturation at short rotation periods.

Vilhu, Osmi↗

Spatial Statistical Data Fusion (SSDF)

As remote sensing for scientific purposes has transitioned from an experimental technology to an operational one, the selection of instruments has become more coordinated, so that the scientific community can exploit complementary measurements. However, tech nological and scientific heterogeneity across devices means that the statistical characteristics of the data they collect are different. The challenge addressed here is how to combine heterogeneous remote sensing data sets in a way that yields optimal statistical estimates of the underlying geophysical field, and provides rigorous uncertainty measures for those estimates. Different remote sensing data sets may have different spatial resolutions, different measurement error biases and variances, and other disparate characteristics. A state-of-the-art spatial statistical model was used to relate the true, but not directly observed, geophysical field to noisy, spatial aggregates observed by remote sensing instruments. The spatial covariances of the true field and the covariances of the true field with the observations were modeled. The observations are spatial averages of the true field values, over pixels, with different measurement noise superimposed. A kriging framework is used to infer optimal (minimum mean squared error and unbiased) estimates of the true field at point locations from pixel-level, noisy observations. A key feature of the spatial statistical model is the spatial mixed effects model that underlies it. The approach models the spatial covariance function of the underlying field using linear combinations of basis functions of fixed size. Approaches based on kriging require the inversion of very large spatial covariance matrices, and this is usually done by making simplifying assumptions about spatial covariance structure that simply do not hold for geophysical variables. In contrast, this method does not require these assumptions, and is also computationally much faster. This method is fundamentally different than other approaches to data fusion for remote sensing data because it is inferential rather than merely descriptive. All approaches combine data in a way that minimizes some specified loss function. Most of these are more or less ad hoc criteria based on what looks good to the eye, or some criteria that relate only to the data at hand.

Braverman, Amy J.↗

TPSAS-NF1676L-12839-DND

Hyper-spectral remote sensors, such as Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer (IASI), provides top of atmospheric radiances with high information content on atmospheric and surface properties. In order to analyze these data in real time, fast, and accurate forward and inverse models are needed. 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 method has been successfully applied to AIRS, IASI, and NAST-I data. By applying the same method to the current and future hyperspectral sounders, the derived products will have less errors due to biases introduced by different retrieval methods. We have compared the retrieved products to radiosondes, aircraft measurements and other validation dataset. Simulations have done to quantify retrieval errors associated with the retrievals. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. 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. 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 subset of the channels. The PCRTM forward model has also been used to performance end-to-end sensor performance simulations of the Climate Absolute Radiance and Refractivity Observatory (CLARREO).

Xu Liu↗

Physical Retrievals of Over-Ocean Rain Rate from Multichannel Microwave Imagery: Theoretical Characteristics of Normalized Polarization and Scattering Indices - Part 1

Microwave rain rate retrieval algorithms have most often been formulated in terms of the raw brightness temperatures observed by one or more channels of a satellite radiometer. Taken individually, single-channel brightness temperatures generally represent a near-arbitrary combination of positive contributions due to liquid water emission and negative contributions due to scattering by ice and/or visibility of the radiometrically cold ocean surface. Unfortunately, for a given rain rate, emission by liquid water below the freezing level and scattering by ice particles above the freezing level are rather loosely coupled in both a physical and statistical sense. Furthermore, microwave brightness temperatures may vary significantly (approx. 30-70 K) in response to geophysical parameters other than liquid water and precipitation. Because of these complications, physical algorithms which attempt to directly invert observed brightness temperatures have typically relied on the iterative adjustment of detailed micro-physical profiles or cloud models, guided by explicit forward microwave radiative transfer calculations. In support of an effort to develop a significantly simpler and more efficient inversion-type rain rate algorithm, the physical information content of two linear transformations of single-frequency, dual-polarization brightness temperatures is studied: the normalized polarization difference P of Petty and Katsaros (1990, 1992), which is intended as a measure of footprint-averaged rain cloud transmittance for a given frequency; and a scattering index S (similar to the polarization corrected temperature of Spencer et al.,1989) which is sensitive almost exclusively to ice. A reverse Monte Carlo radiative transfer model is used to elucidate the qualitative response of these physically distinct single-frequency indices to idealized 3-dimensional rain clouds and to demonstrate their advantages over raw brightness temperatures both as stand-alone indices of precipitation activity and as primary variables in physical, multichannel rain rate retrieval schemes. As a byproduct of the present analysis, it is shown that conventional plane-parallel analyses of the well-known foot-print-filling problem for emission-based algorithms may in some cases give seriously misleading results.

Petty, G. W.↗

Iterative methods for mixed finite element equations

Iterative strategies for the solution of indefinite system of equations arising from the mixed finite element method are investigated in this paper with application to linear and nonlinear problems in solid and structural mechanics. The augmented Hu-Washizu form is derived, which is then utilized to construct a family of iterative algorithms using the displacement method as the preconditioner. Two types of iterative algorithms are implemented. Those are: constant metric iterations which does not involve the update of preconditioner; variable metric iterations, in which the inverse of the preconditioning matrix is updated. A series of numerical experiments is conducted to evaluate the numerical performance with application to linear and nonlinear model problems.

Nakazawa, S.↗

Variable structure control of globally feedback-decoupled deformable vehicle maneuvers

The use of Cayley-Rodrigues attitude parameters as kinematic variables is shown to yield a globally linearized and decoupled model of the equations of motion of a deformable body, where the structural deformation state appears only in the coefficients of the inverse transformation. It is shown how commanded multiaxial attitude maneuvers can be encoded as switching surfaces for a variable-structure control implementation of the corresponding computed slew torques, automatically modulated in respnse only to detected angular rate error signs, for accurate tracking in the presence of separately damped or even uncontrolled (but stable) structural deformations.

Dwyer, T. A. W., III↗

Biophysical and spectral modeling for crop identification and assessment

The development of a technique for estimating all canopy parameters occurring in a canopy reflectance model from the measured canopy reflectance data is summarized. The Suits and the SAIL model for a uniform and homogeneous crop canopy were used to determine if the leaf area index and the leaf angle distribution could be estimated. Optimal solar/view angles for measuring CR were also investigated. The use of CR in many wavelengths or spectral bands and of linear and nonlinear transforms of CRs for various solar/view angles and various spectral bands is discussed as well as the inversion of rediance data inside the canopy, angle transforms for filtering out terrain slope effects, and modification of one dimensional models.

Goel, N. S.↗

On dynamics and control of multi-link flexible space manipulators

In this paper dynamics, inverse dynamics, and control problems for multi-link flexible space manipulators are presented. In deriving the flexible manipulator dynamics the following are assumed: flexible deformations are relatively small; angular rates of the links are much smaller than their fundamental frequencies; nonlinear terms (centrifugal and Coriolis forces) in the flexible manipulator model are the same as those in the rigid body model. These assumptions are reasonable for large space manipulators, such as the space crane. Flexible displacements are measured with respect to the rigid body configuration, for which a linear time-varying system is obtained. The inverse dynamics problem consists of determination of joint torques, given tip trajectory, such that joint angles in flexible configuration are equal to the angles in the rigid body configuration. The manipulator control system consists of the feedforward compensation and feedback control loops. Simulation results of a two-link space crane with large payload show that the performance of this linearized dynamics and control approach is reasonable and robust subject to parameter variations during slew operations.

Gawronski, W.↗

Ultrasonic characterization of changes in viscoelastic properties of epoxy during cure

A technique for using the longitudinal velocity (LV) of an ultrasonic wave to monitor the extent of cross linking (CL) during the cure of thermosetting resins is described. The method was developed by monitoring the rate of change in LV during the cure of a bisphenol-A epoxy resin with an amine adduct. The experiment included variations in the temperature and stoichiometry in order to express the rate of change in terms of the reaction kinetics. The pulse-echo method was used with a single transducer operating at 20 MHz. Numerical models were defined to account for the acoustic response of a single layer, the attenuation and the reflection coefficient. A linear relationship was observed between the inverse of the temperature and the log of the rate of change in the velocity, supporting the theory that the velocity could be used to monitor the extent of the cross-linking reaction. An activation energy of 11.9 kcal/mole was calculated for the mixture being investigated.

Winfree, W. P.↗