Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector fitting”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Geometric accuracy assessment of Landsat-4 Thematic Mapper p-tapes

An analysis was performed on the P-format computer compatible tapes for the Landsat-4 Thematic Mapper (TM) scene 40109-15140 acquired November 2, 1982 over Washington, D.C. and environs. Three tests of sensor geometry were undertaken: (1) band-to-band registration, (2) line-to-line registration between swaths, and (3) geometric correction of sensor and spacecraft/ephemeris characteristics. The band-to-band and line-to-line registration was measured at one hundred pixel spacings along a line using the phase correlation image alignment method developed by Kuglin and Hines and adapted to a one-dimensional FFT correlation technique. The sensor and spacecraft geometric calibration analysis was checked by identifying 75 ground control points in the TM scene and on 1:24,000 topographic maps. A least squares fit of the gaps was computed and the vector offsets of the residuals was tabulated and plotted.

Bryant, N. A.↗
NASA NTRS1993

Numerically Integrated Orbits of the Major Saturnian Satellites fit to Earthbased Observations

We have fit numerically integrated orbits of the eight major satellites of Saturn to all available astrometric and meridian circle observations for the period of 1971 to 1992. The integration was carried out in cartesian coordinates in the J2000 system. The force model included the gravitational effects of the oblate primary, the mutual perturbations of the satellites, and perturbations due to Jupiter and the Sun. Values of the gravitational parameters of the Saturnian system, e.g. planet and satellite masses, were taken from Campbell, et. al., 1989, only the epoch state vectors of the satellites were adjusted to obtain orbits which fit the observations. All astrometric data was processed in the form of satellite relative positions which were weighted according to observer and opposition to reflect the varying data quality...

Jacobson, R. A.↗
NASA NTRS1974

Proper motions, cluster membership and reddening in NGC 6611

Approaches used in reductions for proper motion are based on methods developed by Vasilevskis et al. (1965). Probabilities of cluster membership were determined by fitting the distribution of stars in the proper motion vector diagram with two bivariate Gaussian distributions for the field and the cluster stars. The spectroscopic data indicate that the stellar content of the cluster is not anomalous. The color-magnitude diagram is considered along with questions regarding member stars lying above the main sequence, the cluster age, and aspects of spatial structure.

Kamp, L. W.↗
NASA NTRS1980

Observations of interstellar helium with a gas absorption cell - Implications for the structure of the local interstellar medium

A photometer sensitive at the 584 A line of He 1, incorporating a helium gas resonance absorption cell, was flown on the Apollo-Soyuz Test Project in July 1975. The instrument observed much of the night-time sky, and returned 42 min of usable data. The data were analyzed by fitting to a model of resonant scattering of solar 584 A flux from nearby interstellar helium. Good model fits were obtained for an interstellar gas bulk velocity vector pointing toward alpha = 72 deg, delta = +15 deg, with speed 20 km/s, with interstellar medium temperatures from 5000 to 20,000 K and with neutral interstellar helium density (8.9 plus or minus 10 to the -3rd/cu cm). In the context of theoretical studies of the interstellar medium by McKee and Ostriker (1977), the results may indicate that the sun lies in the warm, partially ionized periphery of a cold interstellar cloud, surrounded by a high-temperature gas heated by old supernova remnants.

Freeman, J.↗
NASA NTRS2009

Hybrid NN/SVM Computational System for Optimizing Designs

A computational method and system based on a hybrid of an artificial neural network (NN) and a support vector machine (SVM) (see figure) has been conceived as a means of maximizing or minimizing an objective function, optionally subject to one or more constraints. Such maximization or minimization could be performed, for example, to optimize solve a data-regression or data-classification problem or to optimize a design associated with a response function. A response function can be considered as a subset of a response surface, which is a surface in a vector space of design and performance parameters. A typical example of a design problem that the method and system can be used to solve is that of an airfoil, for which a response function could be the spatial distribution of pressure over the airfoil. In this example, the response surface would describe the pressure distribution as a function of the operating conditions and the geometric parameters of the airfoil. The use of NNs to analyze physical objects in order to optimize their responses under specified physical conditions is well known. NN analysis is suitable for multidimensional interpolation of data that lack structure and enables the representation and optimization of a succession of numerical solutions of increasing complexity or increasing fidelity to the real world. NN analysis is especially useful in helping to satisfy multiple design objectives. Feedforward NNs can be used to make estimates based on nonlinear mathematical models. One difficulty associated with use of a feedforward NN arises from the need for nonlinear optimization to determine connection weights among input, intermediate, and output variables. It can be very expensive to train an NN in cases in which it is necessary to model large amounts of information. Less widely known (in comparison with NNs) are support vector machines (SVMs), which were originally applied in statistical learning theory. In terms that are necessarily oversimplified to fit the scope of this article, an SVM can be characterized as an algorithm that (1) effects a nonlinear mapping of input vectors into a higher-dimensional feature space and (2) involves a dual formulation of governing equations and constraints. One advantageous feature of the SVM approach is that an objective function (which one seeks to minimize to obtain coefficients that define an SVM mathematical model) is convex, so that unlike in the cases of many NN models, any local minimum of an SVM model is also a global minimum.

Rai, Man Mohan↗
NASA NTRS1993

Development of iterative techniques for the solution of unsteady compressible viscous flows

During the past two decades, there has been significant progress in the field of numerical simulation of unsteady compressible viscous flows. At present, a variety of solution techniques exist such as the transonic small disturbance analyses (TSD), transonic full potential equation-based methods, unsteady Euler solvers, and unsteady Navier-Stokes solvers. These advances have been made possible by developments in three areas: (1) improved numerical algorithms; (2) automation of body-fitted grid generation schemes; and (3) advanced computer architectures with vector processing and massively parallel processing features. In this work, the GMRES scheme has been considered as a candidate for acceleration of a Newton iteration time marching scheme for unsteady 2-D and 3-D compressible viscous flow calculation; from preliminary calculations, this will provide up to a 65 percent reduction in the computer time requirements over the existing class of explicit and implicit time marching schemes. The proposed method has ben tested on structured grids, but is flexible enough for extension to unstructured grids. The described scheme has been tested only on the current generation of vector processor architecture of the Cray Y/MP class, but should be suitable for adaptation to massively parallel machines.

Hixon, Duane↗
NASA NTRS1974

Mars physical parameters as determined from Mariner 9 observations of the natural satellites and Doppler tracking

Mariner 9 Doppler tracking data and television photographs of Deimos and Phobos were analyzed to determine the gravity field, mass, and spin-axis direction of Mars and the natural satellite orbits. The solutions agree with previously published results. Radio data consisted of an apoapsis state vector for each of revolutions 5-195 obtained from one-revolution fits of Doppler data. Optical data consisted of TV photographs of Phobos and Deimos taken between revolutions 25 and 221. A first-order analytical theory, extended to include dominant second-order resonance effects on the Mariner 9 orbit, was used to calculate the motion of the spacecraft, Deimos, and Phobos. The feasibility of combining radio and optical data in long-arc solutions for accurate determination of orbits and physical parameters is demonstrated. The analytical theory developed for the evolution of a highly eccentric orbit in shallow resonance is accurate to plus or minus 1 km in the apoapsis state vector of Mariner 9 over a period of 200 revolutions.

Born, G. H.↗
NASA NTRS1986

Computational methods for the identification of spatially varying stiffness and damping in beams

A numerical approximation scheme for the estimation of functional parameters in Euler-Bernoulli models for the transverse vibration of flexible beams with tip bodies is developed. The method permits the identification of spatially varying flexural stiffness and Voigt-Kelvin viscoelastic damping coefficients which appear in the hybrid system of ordinary and partial differential equations and boundary conditions describing the dynamics of such structures. An inverse problem is formulated as a least squares fit to data subject to constraints in the form of a vector system of abstract first order evolution equations. Spline-based finite element approximations are used to finite dimensionalize the problem. Theoretical convergence results are given and numerical studies carried out on both conventional (serial) and vector computers are discussed.

Banks, H. T.↗
NASA NTRS1987

Ginzburg-Landau theory for the solid-liquid interface of bcc elements

Consideration is given to a simple order-parameter theory for the interfacial tension of body-centered-cubic solids in which the principal order parameter is the amplitude of the density wave at the smallest nonzero reciprocal-lattice vector of the solid. The parameters included in the theory are fitted to the measured heat of fusion, melting temperature, and solid-liquid density difference, and to the liquid structure factor and its temperature derivative at freezing. Good agreement is found with experiment for Na and Fe and the calculated anisotropy of the surface tension among different crystal faces is of the order of 2 percent. On the basis of various assumptions about the universal behavior of bcc crystals at melting, the formalism predicts that the surface tension is proportional to the heat of fusion per surface atom.

Shih, W. H.↗
NASA NTRS1992

A model for the characterization of the spatial properties in vestibular neurons

Quantitative study of the static and dynamic response properties of some otolith-sensitive neurons has been difficult in the past partly because their responses to different linear acceleration vectors exhibited no "null" plane and a dependence of phase on stimulus orientation. The theoretical formulation of the response ellipse provides a quantitative way to estimate the spatio-temporal properties of such neurons. Its semi-major axis gives the direction of the polarization vector (i.e., direction of maximal sensitivity) and it estimates the neuronal response for stimulation along that direction. In addition, the semi-minor axis of the ellipse provides an estimate of the neuron's maximal sensitivity in the "null" plane. In this paper, extracellular recordings from otolith-sensitive vestibular nuclei neurons in decerebrate rats were used to demonstrate the practical application of the method. The experimentally observed gain and phase dependence on the orientation angle of the acceleration vector in a head-horizontal plane was described and satisfactorily fit by the response ellipse model. In addition, the model satisfactorily fits neuronal responses in three-dimensions and unequivocally demonstrates that the response ellipse formulation is the general approach to describe quantitatively the spatial properties of vestibular neurons.

Non-NASA Center↗
NASA NTRS1991

A Newton algorithm for complex curve fitting

The problem of synthesizing transfer functions from frequency response measurements is considered. Given a complex vector representing the measured frequency response of a physical system, a transfer function of specified order is determined that minimizes the sum of the magnitude-squared of the frequency response errors. This nonlinear least squares minimization problem is solved by an iterative global descent algorithm of the Newton type which converges quadratically near the minimum. The unknown transfer function is expressed as a sum of second order rational polynomials, a parameterization that facilitates a numerically robust computer implementation. The algorithm is developed for single-input, single-output, causal, stable transfer functions.

Spanos, J. T.↗
NASA NTRS1993

A model for the motion of the Philippine Sea plate consistent with NUVEL-1 and geological data

We investigate angular velocity vectors of the Philippine Sea (PH) plate relative to the adjacent major plates, Eurasia (EU) and Pacific (PA), and the smaller Caroline (CR) plate. Earthquake slip vector data along the Philippine Sea plate are inverted, subject to the constraint that EU-PA motion equals that predicted by the global relative plate model NUVEL-1. The resulting solution fails to satisfy geological constraints along the Caroline-Pacific boundary: convergence along the Mussau Trench and divergence along the Sorol Trough. We then seek solutions satisfying both the CR-PA boundary conditions and the Philippine Sea slip vector data, by adjusting the PA-PH and EU-PH best fitting poles within their error ellipses. We also consider northern Honshu to be part of the North American plate and impose the constraint that the Philippine Sea plate subducts beneath northern Honshu along the Sagmi Trough in a NNW-NW direction. Of the solutions satisfying these conditions, we select the best EU-PH as 48.2 deg N, 157.0 deg E, 1.09 deg/my, corresponding to a pole far from Japan and south of Kamchatka, and PA-PH, 1.2 deg N, 134.2 deg E, 1.00 deg/my. Predicted NA-PH and EU-PH convergence rates in central Honshu are consistent with estimated seismic slip rates. Previous estimates of the EU-PH pole close to central Honshu are inconsistent with extension within the Bonin backarc implied by earthquake slip vectors and NNW-NW convergence of the Bonin forearc at the Sagami Trough.

Seno, Tetsuzo↗
NASA NTRS1991

Algorithms for l2 and l-infinity transfer function curve fitting

In this paper algorithms for fitting transfer functions to frequency response data are developed. Given a complex vector representing the measured frequency response of a physical system, a transfer function of specified order is determined that minimizes either of the following criteria: (1) the sum of the magnitude-squared of the frequency response errors, and (2) the magnitude of the maximum error. Both of these criteria are nonlinear in the coefficients of the unknown transfer function, and iterative minimization algorithms are proposed. A numerical example demonstrates the effectiveness of the proposed algorithms.

Spanos, John T.↗
NASA NTRS1999

Backus Effect on a Perpendicular Errors in Harmonic Models of Real vs. Synthetic Data

Measurements of geomagnetic scalar intensity on a thin spherical shell alone are not enough to separate internal from external source fields; moreover, such scalar data are not enough for accurate modeling of the vector field from internal sources because of unmodeled fields and small data errors. Spherical harmonic models of the geomagnetic potential fitted to scalar data alone therefore suffer from well-understood Backus effect and perpendicular errors. Curiously, errors in some models of simulated 'data' are very much less than those in models of real data. We analyze select Magsat vector and scalar measurements separately to illustrate Backus effect and perpendicular errors in models of real scalar data. By using a model to synthesize 'data' at the observation points, and by adding various types of 'noise', we illustrate such errors in models of synthetic 'data'. Perpendicular errors prove quite sensitive to the maximum degree in the spherical harmonic expansion of the potential field model fitted to the scalar data. Small errors in models of synthetic 'data' are found to be an artifact of matched truncation levels. For example, consider scalar synthetic 'data' computed from a degree 14 model. A degree 14 model fitted to such synthetic 'data' yields negligible error, but amplifies 4 nT (rmss) added noise into a 60 nT error (rmss); however, a degree 12 model fitted to the noisy 'data' suffers a 492 nT error (rmms through degree 12). Geomagnetic measurements remain unaware of model truncation, so the small errors indicated by some simulations cannot be realized in practice. Errors in models fitted to scalar data alone approach 1000 nT (rmss) and several thousand nT (maximum).

Voorhies, C. V.↗
NASA NTRS1977

Computer solutions for studying correlations between solar magnetic fields and Skylab X-ray observations

A method is described which correlates the NASA-Marshall Space Flight Center (MSFC) Image Data Processing System (IDAPS) and MSFC magnetograph data to X-ray and H-alpha observations from the Skylab mission. Solutions of Laplace's equation in three dimensions, based on the magnetograph data, are convolved with observed X-ray and H-alpha regions. Matched filtering (template matching) provides a best fit of the observed X-ray regions to the computed total magnetic vector magnitude between 10,000 and 15,000 km above the photosphere.

Teuber, D.↗
NASA NTRS1993

Extension of state-vector splitting to the Navier-Stokes equations

A multi-dimensional upwind scheme for solving the Euler equations has been extended to solve the Navier-Stokes equations. These methods share the property of not requiring body-fitted computational grids. A time advance scheme which is second-order accurate in time and space has been developed for use in conjunction with state-vector splitting to accurately resolve regions with significant gradients in flow quantities. The scheme has been tested on several simple viscous flow problems for which analytic solutions are known; the agreement of computation with theory is excellent.

Gooch, Carl F.↗
NASA NTRS1992

Registering and resampling images in STSDAS

Registering different images can be difficult, especially if the images to be registered are images at different wavelengths, where features in one image may look entirely different or be absent from the second image. Using two new packages soon to be added to the STSDAS package, REGISTER and RESAMPLE, this job is done automatically. The REGISTER package allows the user to determine the amount of translation, rotation, and/or magnification needed to make two images, spectra, or time series congruent. The methods implemented to compute the registration parameters use: a set of the pixel coordinates of the same features identified in two files, or the FITS coordinate transformation parameters in the headers of two data files, or a single feature identified as the peak of a cross-correlation between two vectors. The coefficients describing the registration are defined by the equations (for a two-dimensional image): x = a + bx + cy, y = d + ex + fy, where (X,Y) are the pixel coordinates of a feature in the reference image, (x,y) are the pixel coordinates of a feature in the secondary image, and the computed coefficients are a, b, c, d, e, and f. Results may be produced by linking the output of REGISTER to RESAMPLE in a command language procedure. The output from REGISTER and the input to RESAMPLE consists of a matrix of coefficients (a through f above) fully specifying the registration. The RESAMPLE package resamples simple vector or image data for a given amount of translation, rotation (images only), and magnification, or reflection of the science data. Specific options included are: image rotation about the FITS reference pixel; scale changes, i.e. magnification or demagnification (for images, independently on both axes); simple translation; reflection (for images, about one or both axes); and resampling and registration to a reference data set. Output from the RESAMPLE task is the resampled image which may then be displayed and compared with the reference image.

Williamson, R. L., II↗
NASA NTRS1979

Comet Encke - Precession of the spin axis, nongravitational motion, and sublimation

From the observed light curve of P/Encke the jet force from sublimation is calculated both as a (precessing) torque and as a (perturbing) force transverse to the radius vector. An integral iteration is carried out over 59 perihelion passages, 1786-1977, to fit the previously determined nongravitational transverse force and to derive the precession of the spin axis. It is shown that the spin axis turned more than 100 degrees in longitude and almost 30 degrees in latitude from 1786 to 1977, but appears to have been almost fixed in direction for hundreds of revolutions before 1700. It is suggested that ejected meteoroidal debris accumulated on the currently less active hemisphere, insulating it to maintain a low activity level. A tentative rotation period of 6 h 33 min is derived, using Whipple's halo method. The suggested spinup rate is 21 min/century, while the current rate of relative mass loss by sublimation is 0.09% of the comet's mass per revolution. Moreover, the mass of the nucleus is estimated at less than 10 to the 16th grams, and its oblateness at less than 4%.

Whipple, F. L.↗