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 613 records · Page 34

Simultaneous solution for core magnetic field and fluid flow beneath an electrically conducting mantle

The effects of laterally homogeneous mantle electrical conductivity were included in steady, frozen-flux core surface flow estimation along with refinements in method and weighting. The refined method allows simultaneous solution for both the initial radial geomagnetic field component at the core-mantle boundary (CMB) and the sub-adjacent fluid motion; it also features Gauss' method for solving the non-linear inverse problem associated with steady motional induction. The tradeoff between spatial complexity of the derived flows and misfit to the weighted Definitive Geomagnetic Reference Field models (DGRF's) is studied for various mantle conductivity profiles. For simple flow and a fixed initial geomagnetic condition, a fairly high deep-mantle conductivity performs better than either insulating or weakly conducting profiles; however, a thin, very high conductivity layer at the base of the mantle performs almost as well. Simultaneous solution for both initial geomagnetic field and flow reduces the misfit per degree of freedom even more than does changing the mantle conductivity profile. Moreover, when both core field and flow are estimated, the performance of the solutions and the derived flows become insensitive to the conductivity profile.

Voorhies, Coerte V.↗

Control theory based airfoil design using the Euler equations

This paper describes the implementation of optimization techniques based on control theory for airfoil design. In our previous work it was shown that control theory could be employed to devise effective optimization procedures for two-dimensional profiles by using the potential flow equation with either a conformal mapping or a general coordinate system. The goal of our present work is to extend the development to treat the Euler equations in two-dimensions by procedures that can readily be generalized to treat complex shapes in three-dimensions. Therefore, we have developed methods which can address airfoil design through either an analytic mapping or an arbitrary grid perturbation method applied to a finite volume discretization of the Euler equations. Here the control law serves to provide computationally inexpensive gradient information to a standard numerical optimization method. Results are presented for both the inverse problem and drag minimization problem.

Jameson, Antony↗

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗

Some experience using subcritical response methods in wind-tunnel flutter model studies

Experiences obtained with four methods of predicting flutter of wind-tunnel models from subcritical response data are described. The four methods are: co/quad, randomdec, power spectra density, and the peak-hold spectrum. Model excitation techniques included both forced (sinusoidal sweep) and random (tunnel turbulence). These methods were successfully used to measure the frequency and damping (or an inverse response amplitude proportional to the damping) in the predominant flutter modes. Implementation and application of each method are discussed. Some results and comparisons between methods are presented.

Foughner, J. T., Jr.↗

Investigating Sources of Variability and Error in Simulations of Carbon Dioxide in an Urban Region

Greenhouse gas (GHG) emissions estimation methods that use atmospheric trace gas observations, including inverse modeling techniques, perform better when carbon dioxide (CO2) fluxes are more accurately transported and dispersed in the atmosphere by a numerical model. In urban areas, transport and dispersion is particularly difficult to simulate using current mesoscale meteorological models due, in part, to added complexity from surface heterogeneity and fine spatial/temporal scales. It is generally assumed that the errors in GHG estimation methods in urban areas are dominated by errors in transport and dispersion. Other significant errors include, but are not limited to, those from assumed emissions magnitude and spatial distribution. To assess the predictability of simulated trace gas mole fractions in urban observing systems using a numerical weather prediction model, we employ an Eulerian model that combines traditional meteorological variables with multiple passive tracers of atmospheric CO2 from anthropogenic inventories and a biospheric model. The predictability of the Eulerian model is assessed by comparing simulated atmospheric CO2 mole fractions to observations from four in situ tower sites (three urban and one rural) in the Washington DC/Baltimore, MD area for February 2016. Four different gridded fossil fuel emissions inventories along with a biospheric flux model are used to create an ensemble of simulated atmospheric CO2 observations within the model. These ensembles help to evaluate whether the modeled observations are impacted more by the underlying emissions or transport. The spread of modeled observations using the four emission fields indicates the model's ability to distinguish between the different inventories under various meteorological conditions. Overall, the Eulerian model performs well; simulated and observed average CO2 mole fractions agree within 1% when averaged at the three urban sites across the month. However, there can be differences greater than 10% at any given hour, which are attributed to complex meteorological conditions rather than differences in the inventories themselves. On average, the mean absolute error of the simulated compared to actual observations is generally twice as large as the standard deviation of the modeled mole fractions across the four emission inventories. This result supports the assumption, in urban domains, that the predicted mole fraction error relative to observations is dominated by errors in model meteorology rather than errors in the underlying fluxes in winter months. As such, minimizing errors associated with atmospheric transport and dispersion may help improve the performance of GHG estimation models more so than improving flux priors in the winter months. We also find that the errors associated with atmospheric transport in urban domains are not restricted to certain times of day. This suggests that atmospheric inversions should use CO2 observations that have been filtered using meteorological observations rather than assuming that meteorological modeling is most accurate at certain times of day (such as using only mid-afternoon observations).

CO2↗

Inverse dynamics of a 3 degree of freedom spatial flexible manipulator

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

Bayo, Eduardo↗

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, H.↗

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid points are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, Hillel↗

Land-surface temperature measurement from space - Physical principles and inverse modeling

To apply the multiple-wavelength (split-window) method used for satellite measurement of sea-surface temperature from thermal-infrared data to land-surface temperatures, the authors statistically analyze simulations using an atmospheric radiative transfer model. The range of atmospheric conditions and surface temperatures simulated is wide enough to cover variations in clear atmospheric properties and surface temperatures, both of which are larger over land than over sea. Surface elevation is also included in the simulation as the most important topographic effect. Land covers characterized by measured or modeled spectral emissivities include snow, clay, sands, and tree leaf samples. The empirical inverse model can estimate the surface temperature with a standard deviation less than 0.3 K and a maximum error less than 1 K, for viewing angles up to 40 degrees from nadir under cloud-free conditions, given satellite measurements in three infrared channels. A band in the region from 10.2 to 11.0 microns will usually give the most reliable single-band estimate of surface temperature. In addition, a band in either the 3.5-4.0-micron region or in the 11.5-12.6-micron region must be included for accurate atmospheric correction, and a band below the ozone absorption feature at 9.6 microns (e.g., 8.2-8.8 microns) will increase the accuracy of the estimate of surface temperature.

Wan, Zhengming↗

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

Effect of Malmquist bias on correlation studies with IRAS data base

The relationships between galaxy properties in the sample of Trinchieri et al. (1989) are reexamined with corrections for Malmquist bias. The linear correlations are tested and linear regressions are fit for log-log plots of L(FIR), L(H-alpha), and L(B) as well as ratios of these quantities. The linear correlations for Malmquist bias are corrected using the method of Verter (1988), in which each galaxy observation is weighted by the inverse of its sampling volume. The linear regressions are corrected for Malmquist bias by a new method invented here in which each galaxy observation is weighted by its sampling volume. The results of correlation and regressions among the sample are significantly changed in the anticipated sense that the corrected correlation confidences are lower and the corrected slopes of the linear regressions are lower. The elimination of Malmquist bias eliminates the nonlinear rise in luminosity that has caused some authors to hypothesize additional components of FIR emission.

Verter, Frances↗

Test images for the maximum entropy image restoration method

One of the major activities of any experimentalist is data analysis and reduction. In solar physics, remote observations are made of the sun in a variety of wavelengths and circumstances. In no case is the data collected free from the influence of the design and operation of the data gathering instrument as well as the ever present problem of noise. The presence of significant noise invalidates the simple inversion procedure regardless of the range of known correlation functions. The Maximum Entropy Method (MEM) attempts to perform this inversion by making minimal assumptions about the data. To provide a means of testing the MEM and characterizing its sensitivity to noise, choice of point spread function, type of data, etc., one would like to have test images of known characteristics that can represent the type of data being analyzed. A means of reconstructing these images is presented.

Mackey, James E.↗

On strategies for inverting remote sensing data

Attention is given to a number of methods for inverting remote sensing data obtained in a variety of astronomical applications. Applications include image restoration, inversion of helioseismological data to obtain the internal rotation rate of stars such as the sun, fitting of spectra (especially thermal line spectra) to grating or other dispersed observed spectra, differential emission measure analysis, and reconstruction of images derived from interferometric observations. The results consider the tradeoff between resolution and variance and the stability properties for each method and propose an inversion stragegy using the available techniques.

Jeffrey, W.↗

Optimal application of Morrison's iterative noise removal for deconvolution. Appendices

Morrison's iterative method of noise removal, or Morrison's smoothing, is applied in a simulation to noise-added data sets of various noise levels to determine its optimum use. Morrison's smoothing is applied for noise removal alone, and for noise removal prior to deconvolution. For the latter, an accurate method is analyzed to provide confidence in the optimization. The method consists of convolving the data with an inverse filter calculated by taking the inverse discrete Fourier transform of the reciprocal of the transform of the response of the system. Various length filters are calculated for the narrow and wide Gaussian response functions used. Deconvolution of non-noisy data is performed, and the error in each deconvolution calculated. Plots are produced of error versus filter length; and from these plots the most accurate length filters determined. The statistical methodologies employed in the optimizations of Morrison's method are similar. A typical peak-type input is selected and convolved with the two response functions to produce the data sets to be analyzed. Both constant and ordinate-dependent Gaussian distributed noise is added to the data, where the noise levels of the data are characterized by their signal-to-noise ratios. The error measures employed in the optimizations are the L1 and L2 norms. Results of the optimizations for both Gaussians, both noise types, and both norms include figures of optimum iteration number and error improvement versus signal-to-noise ratio, and tables of results. The statistical variation of all quantities considered is also given.

Ioup, George E.↗

Feasibility study for automatic reduction of phase change imagery

The feasibility of automatically reducing a form of pictorial aerodynamic heating data is discussed. The imagery, depicting the melting history of a thin coat of fusible temperature indicator painted on an aerodynamically heated model, was previously reduced by manual methods. Careful examination of various lighting theories and approaches led to an experimentally verified illumination concept capable of yielding high-quality imagery. Both digital and video image processing techniques were applied to reduction of the data, and it was demonstrated that either method can be used to develop superimposed contours. Mathematical techniques were developed to find the model-to-image and the inverse image-to-model transformation using six conjugate points, and methods were developed using these transformations to determine heating rates on the model surface. A video system was designed which is able to reduce the imagery rapidly, economically and accurately. Costs for this system were estimated. A study plan was outlined whereby the mathematical transformation techniques developed to produce model coordinate heating data could be applied to operational software, and methods were discussed and costs estimated for obtaining the digital information necessary for this software.

Nossaman, G. O.↗

On Characterizing Particle Shape

It is well known that particle shape affects flow characteristics of granular materials, as well as a variety of other solids processing issues such as compaction, rheology, filtration and other two-phase flow problems. The impact of shape crosses many diverse and commercially important applications, including pharmaceuticals, civil engineering, metallurgy, health, and food processing. Two applications studied here include the dry solids flow of lunar simulants (e.g. JSC-1, NU-LHT-2M, OB-1), and the flow properties of wet concrete, including final compressive strength. A multi-dimensional generalized, engineering method to quantitatively characterize particle shapes has been developed, applicable to both single particle orientation and multi-particle assemblies. The two-dimension, three dimension inversion problem is also treated, and the application of these methods to DEM model particles will be discussed. In the case of lunar simulants, flow properties of six lunar simulants have been measured, and the impact of particle shape on flowability - as characterized by the shape method developed here -- is discussed, especially in the context of three simulants of similar size range. In the context of concrete processing, concrete construction is a major contributor to greenhouse gas production, of which the major contributor is cement binding loading. Any optimization in concrete rheology and packing that can reduce cement loading and improve strength loading can also reduce currently required construction safety factors. The characterization approach here is also demonstrated for the impact of rock aggregate shape on concrete slump rheology and dry compressive strength.

Ennis, Bryan J.↗

Interpretation of solar extinction data for stratospheric aerosols

This paper discusses the inversion problem for aerosols using the solar extinction method. A series of numerical experiments is described in which solar extinction measurement systems are modeled. A numerical model of a solar extinction measurement system has been coupled with model atmospheres that exhibit fine scale structures to produce numerically generated data signals. These signals were then inverted to study the effect that measurement errors and desired vertical resolution produce in the inverted results. Knowledge o2 the trade off between vertical resolution and the accuracy of inversion aid in the interpretation of the inverted results.

Pepin, T. J.↗

Soot Volume Fraction Imaging

A new technique is described for the full-field determination of soot volume fractions via laser extinction measurements. This technique differs from previously reported point-wise methods in that a two-dimensional array (i.e., image) of data is acquired simultaneously. In this fashion, the net data rate is increased, allowing the study of time-dependent phenomena and the investigation of spatial and temporal correlations. A telecentric imaging configuration is employed to provide depth-invariant magnification and to permit the specification of the collection angle for scattered light. To improve the threshold measurement sensitivity, a method is employed to suppress undesirable coherent imaging effects. A discussion of the tomographic inversion process is provided, including the results obtained from numerical simulation. Results obtained with this method from an ethylene diffusion flame are shown to be in close agreement with those previously obtained by sequential point-wise interrogation.

Greenberg, Paul S.↗