Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “least-squares analysis”

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 145 records · Page 8

Development of the NASA/FLAGRO computer program

A preliminary version of a new computer program has been developed for calculating crack propagation in cyclically loaded structures. The new program, titled NASA/FLAGRO, has numerous enhanced features over previous programs used for safe life analysis on space systems. The essential features of the new program consist of modular design, promptive input using a cathode ray tube terminal, improved crack growth rate analytical methods, many new crack case solutions, and incorporation of user-defined files for material properties and fatigue spectrums. Also, five options made available in the program are to compute: (1) safe life of a part with a preexisting crack; (2) critical crack size for a specified stress level; (3) stress-intensity factor numerical values for making comparisons or checks of solutions; (4) least-squares fit of growth rate equations to crack growth rate data; and (5) the conversion of a versus N data to Delta K versus da/dN.

Forman, Royce G.↗

Catalog-to-catalog reductions - Results for the FK, N30, and GC catalogs

It is shown that previous computations of catalog-to-catalog differences, whether by the bin-and-average method or by least-squares adjustment combined with an expansion in an incomplete set of basis functions, have all been performed at inappropriate epoch(s) of place. As a consequence, such results are frequently much more informative regarding the systematic differences between the proper motion systems than they are with respect to the systematic differences between the positional systems. Moreover, and for the same reason, some sets of proper motion differences really do not portray the true state of the evolution of the proper motion systems. Finally, this has practical effects in other astrometrical work, for instance, in the analysis of earth rotation time series. This error is corrected by arguing for a more meaningful epoch of place for catalog-to-catalog comparisons. This method of infinitely overlapping circles is applied at these other epochs to the FK sequence of catalogs, the N30, and the GC, and as a corollary the effectiveness of this new technique is shown.

Taff, L. G.↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

When the radio link between two satellites (GEOS-3/ATS-6) is intercepted by the earth's ionosphere and neutral atmosphere, a change in the Doppler frequency results. Travel through the atmosphere causes the Doppler phase to be advanced in the ionosphere's portion and retarded in the neutral portion of the atmosphere. Analysis of the shortening and lengthening of the phase of the Satellite-to-Satellite Tracking (SST) data that passed within 40-700 km above the earth's surface during its ATS-6 to GEOS-3 to ATS-6 path, caused by the atmosphere, results in refractivity versus height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbation from the Moon, Sun and a 15th order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occultation period and a model ionosphere was estimated by a least-square adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the 2 methods were caused by orbital errors, which propagated into the solution. The SST data yielded refractive index profiles in a novel economical manner because no additional or special on-board equipment were required.

Liu, A. S.↗

Doppler line profile analysis for a multichannel Fabry-Perot interferometer

A new method of instrument calibration and data analysis is presented for single-etalon interferometric measurements of winds, temperatures, and emission line intensities. The technique has been developed for the multichannel Fabry-Perot interferometer on the Dynamics Explorer spacecraft. A numerical representation of the instrumental transfer function is used based on a truncated Fourier series with empirically determined coefficients. The numerical form is compared with the conventional analytic form. The Fourier coefficients describing the instrument function are generated at the wavelength of a stable He-Ne laser and are translated to other wavelengths using an interpolation technique for both phase and power. A quasi-linear least-squares fitting process involving matrices provides for a rapid and accurate data reduction.

Killeen, T. L.↗

State-space formulation of multi-shaker modal analysis

As modal testing techniques have improved in recent years, the demands placed on the accuracy of modal test results have increased. This paper describes a new multi-input, multi-output modal parameter estimation algorithm suitable for the identification of reduced-order models of structures. The method described is a frequency-domain method based on excitation and response spectra. The principal features of the algorithm are: (1) applicability to general linear, time-invariant systems, (2) direct use of multiple input and multiple response data, (3) identification of reduced-order models with minimum user interaction through the use of two automatic model-reduction techniques, (4) identification of a consistent set of modal parameters through the use of numerically stable least-squares techniques, and (5) complete user documentation.

Craig, Roy R., Jr.↗

Meteoroid and debris special investigation group; status of 3-D crater analysis from binocular imagery

During the 3 month deintegration of the LDEF, the M&D SIG generated approximately 5000 digital color stereo image pairs of impact related features from all space exposed surfaces. Currently, these images are being processed at JSC to yield more accurate feature information. Work is currently underway to determine the minimum number of data points necessary to parametrically define impact crater morphologies in order to minimize the man-hour intensive task of tie point selection. Initial attempts at deriving accurate crater depth and diameter measurements from binocular imagery were based on the assumption that the crater geometries were best defined by paraboloid. We made no assumptions regarding the crater depth/diameter ratios but instead allowed each crater to define its own coefficients by performing a least-squares fit based on user-selected tiepoints. Initial test cases resulted in larger errors than desired, so it was decided to test our basic assumptions that the crater geometries could be parametrically defined as paraboloids. The method for testing this assumption was to carefully slice test craters (experimentally produced in an appropriate aluminum alloy) vertically through the center resulting in a readily visible cross-section of the crater geometry. Initially, five separate craters were cross-sectioned in this fashion. A digital image of each cross-section was then created, and the 2-D crater geometry was then hand-digitized to create a table of XY position for each crater. A 2nd order polynomial (parabolic) was fitted to the data using a least-squares approach. The differences between the fit equation and the actual data were fairly significant, and easily large enough to account for the errors found in the 3-D fits. The differences between the curve fit and the actual data were consistent between the caters. This consistency suggested that the differences were due to the fact that a parabola did not sufficiently define the generic crater geometry. Fourth and 6th order equations were then fitted to each crater cross-section, and significantly better estimates of the crater geometry were obtained with each fit. Work is presently underway to determine the best way to make use of this new parametric crater definition.

Sapp, Clyde A.↗

An improved plate theory of order (1,2) for thick composite laminates

A new (1,2)-order theory is proposed for the linear elasto-static analysis of laminated composite plates. The basic assumptions are those concerning the distribution through the laminate thickness of the displacements, transverse shear strains and the transverse normal stress, with these quantities regarded as some weighted averages of their exact elasticity theory representations. The displacement expansions are linear for the inplane components and quadratic for the transverse component, whereas the transverse shear strains and transverse normal stress are respectively quadratic and cubic through the thickness. The main distinguishing feature of the theory is that all strain and stress components are expressed in terms of the assumed displacements prior to the application of a variational principle. This is accomplished by an a priori least-square compatibility requirement for the transverse strains and by requiring exact stress boundary conditions at the top and bottom plate surfaces. Equations of equilibrium and associated Poisson boundary conditions are derived from the virtual work principle. It is shown that the theory is particularly suited for finite element discretization as it requires simple C(sup 0)- and C(sup -1)-continuous displacement interpolation fields. Analytic solutions for the problem of cylindrical bending are derived and compared with the exact elasticity solutions and those of our earlier (1,2)-order theory based on the assumed displacements and transverse strains.

Tessler, A.↗

Center of mass estimation in closed vortices - A verification in principle and practice

The problem of tracking closed mesoscale vortices using center of mass estimation techniques is studied. Three estimators are evaluated using data from a warm core Gulf Stream ring. The comparisons show that a method based on the intersection of perpendicular bisectors and one using a least-squares fit of a conic section perform comparably. The perpendicular bisector algorithm is used in conjunction with a Gaussian ring model and a star-shaped survey pattern to produce an expected error curve as a function of vortex translation, survey speed and vortex size. For typical ring parameters, center estimation is usually possible to within + or - 5 km. The feasibility of using differing data sets to construct a history of ring motion based on a coordinate system moving with the ring is also investigated. In this way, the validity of using satellite-derived data and drifter trajectories to estimate the center of mass of a mesoscale feature is assessed. The results of the analysis demonstrate that the location of the deeper structure of the ring and the surface expression are sufficiently well correlated to permit dynamically relevant calculations based on surface measurements. It is shown that satellite-derived data can be used to approximate the center of mass trajectory to within the error in the individual center estimates for the period analyzed. The Lagrangian-drifter-derived centers are offset from the center of mass trajectory in a manner consistent with kinematic arguments.

Hooker, S. B.↗

Line-frequency measurements and analysis of N2O between 900 and 4700/cm

The N2O infrared spectrum from 900 to 4700/cm has been studied with a high-resolution Fourier-transform spectrometer. Measurements were made of line frequencies of several N2O isotopes for several ground-state and hot bands. These bands include the 1530-0330, 2330-0330, 0441-0440, 0441-0330, 0531-0330, and 1331-0330 bands of (N-14)2(O-16). Also, perturbation-enhanced transitions in the 0710-0110 and 0730-0110 bands were assigned and measured. The frequencies have been analyzed to obtain a unique set of effective vibration-rotation parameters for each vibrational state. The results obtained from this research are of high accuracy for the ground-state, and the first-excited-state bands of (N-14)2(O-16), in which the computed frequencies derived from the least-squares fits, are known to an absolute uncertainty of only +/-4 x 10 exp -5 for nonperturbed transitions.

Toth, Robert A.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗

DFW microburst model based on AA-539 data

Analysis of the August 2, 1985 crash for an L-1011 jumbo jet (DL-191) on approach to the Dallas-Ft. Worth International Airport (DFW) in a thunderstorm indicates that the severe windshear microburst that caused the crash was composed not only of a strong downflow and outflow but also included several large-scale vortex rings entrained in the flowfield. This paper presents a detailed two-dimensional model of the DFW microburst based on data from the MD-80 (AA-539) that followed behind DL-191 and flew through the microburst about two minutes after the crash of DL-191. The model was developed using wind-vector and flight-path data reconstructed by NASA Ames Research Center and a combination of interactive graphics and least-squares error best fit between the modeled and measured wind vectors along the AA-539 flight path. The model indicates that the flowfield contains some significant elements and vortices not previously reported. The alternating direction of rotation of the vortices in the model suggests a microburst structure based on a von Karman vortex street rather than on a Kelvin-Helmholtz instability. The model also indicates that the reconstructed wind-vector data contain a time lag of at least one second in the horizontal winds.

Grantham, Walter J.↗

The relationship between IR, optical, and UV extinction

An analysis is presented for the variability of absolute IR, optical, and UV extinction, A(sub lambda), derived through the ratio of total-to-selective extinction, R, for 31 lines of sight for which reliable UV extinction parameters were derived. These data sample a wide range of environments and are characterized by 2.5 is less than or equal to R is less than or equal to 6.0. It was found that there is a strong linear dependence between extinction expressed as A(sub lambda)/A(sub V) and 1/R for 1.25 micron is less than or equal to lambda is less than or equal to 0.12 micron. Differences in the general shape of extinction curves are largely due to variations in shape of optical/near-UV extinction corresponding to changes in R, with A(sub lambda)/A(sub V) decreasing for increasing R. From a least-squares fit of the observed R-dependence as a function of wavelength for 0.8/micron is less than or greater than 1/lambda is less than or equal to 8.3/micron, an analytic expression was generated from which IR, optical, and UV extinction curves of the form A(sub lambda)/A(sub V) can be reproduced with reasonable accuracy from a knowledge of R. It was also found that the absolute bump strength normalized to A(sub V) shows a general decrease with increasing R, suggesting that some fraction of bump grains may be selectively incorporated into coagulated grains. Finally, it was found that absolute extinction normalized by suitably chosen color indices results in a minimization of the R-dependence of portions of the UV curve, allowing A(sub lambda) to be estimated for these wavelengths independent of R.

Cardelli, Jason A.↗

Extracting galactic structure parameters from multivariated density estimation

Multivariate statistical analysis, including includes cluster analysis (unsupervised classification), discriminant analysis (supervised classification) and principle component analysis (dimensionlity reduction method), and nonparameter density estimation have been successfully used to search for meaningful associations in the 5-dimensional space of observables between observed points and the sets of simulated points generated from a synthetic approach of galaxy modelling. These methodologies can be applied as the new tools to obtain information about hidden structure otherwise unrecognizable, and place important constraints on the space distribution of various stellar populations in the Milky Way. In this paper, we concentrate on illustrating how to use nonparameter density estimation to substitute for the true densities in both of the simulating sample and real sample in the five-dimensional space. In order to fit model predicted densities to reality, we derive a set of equations which include n lines (where n is the total number of observed points) and m (where m: the numbers of predefined groups) unknown parameters. A least-square estimation will allow us to determine the density law of different groups and components in the Galaxy. The output from our software, which can be used in many research fields, will also give out the systematic error between the model and the observation by a Bayes rule.

Chen, B.↗

An airborne study of microwave surface sensing and boundary layer heat and moisture fluxes for Project FIFE

This is the third semiannual status report describing the University of Wyoming participation in project FIFE, under NASA Grant NAG 5-913. The reader is referred to the first and second semiannual status reports for this grant for descriptions of University of Wyoming participation in the field operations of 1987. This report is concerned with continuing analysis of field data. Most of the instrument calibration effort has concentrated on interpretation and use of the high-rate water vapor measurements with the University of Wyoming Lyman-alpha device. In addition, fluxes of water vapor have been corrected for sensible heat flux, and all the high-rate (10 Hz for the Wyoming King Air) data have been filtered before the flux calculations. For each flight the Lyman-alpha data (10 Hz) are calibrated against the Cambridge dew-point hygrometer (1 Hz), by using data from the take-off sounding. The reasoning here is that the takeoff sounding (held to 500 fpm rate-of-climb specifically for the calibration) covers the widest range of humidity conditions of any sequence of the flight. For each takeoff sounding, then, the Lyman-alpha data, averaged to 1 Hz, are matched to the dewpoint data with a second-order least-squares fit. The resulting quadratic equation is then used in all subsequent analyses to convert the Lyman-alpha voltages to vapor desities for that specific flight. An alternative method is to use the range of fluctuation for the same two devices from horizontal passes to establish the calibration function. However, since the Lyman-alpha data are used only in the flux calculations, i.e., only the fluctuation and not the absolute values are used, this method is thought unnecessary.

Kelly, Robert D.↗

Least-squares sequential parameter and state estimation for large space structures

This paper presents the formulation of simultaneous state and parameter estimation problems for flexible structures in terms of least-squares minimization problems. The approach combines an on-line order determination algorithm, with least-squares algorithms for finding estimates of modal approximation functions, modal amplitudes, and modal parameters. The approach combines previous results on separable nonlinear least squares estimation with a regression analysis formulation of the state estimation problem. The technique makes use of sequential Householder transformations. This allows for sequential accumulation of matrices required during the identification process. The technique is used to identify the modal prameters of a flexible beam.

Thau, F. E.↗

Optimal focal-plane restoration

Image restoration can be implemented efficiently by calculating the convolution of the digital image and a small kernel during image acquisition. Processing the image in the focal-plane in this way requires less computation than traditional Fourier-transform-based techniques such as the Wiener filter and constrained least-squares filter. Here, the values of the convolution kernel that yield the restoration with minimum expected mean-square error are determined using a frequency analysis of the end-to-end imaging system. This development accounts for constraints on the size and shape of the spatial kernel and all the components of the imaging system. Simulation results indicate the technique is effective and efficient.

Reichenbach, Stephen E.↗