Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized least squared”

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 253 records · Page 14

Decomposition of Electron Ionization Mass Spectra for Space Application Using a Monte-Carlo approach

Rationale: Quadrupole mass spectrometers equipped with an electron ionization (EI) sources have been widely used in space exploration to investigate the composition of planetary surfaces and atmospheres. However, the complexity of the samples and the minimal calibration for the fragmentation of molecules in the ionization chambers have prevented the deconvolution of the majority of the mass spectra obtained at different targets, thus limiting the determination of the exact composition of the samples analyzed. We propose a Monte‐Carlo approach to solve this issue mathematically. Methods: We decomposed simulated mass spectra of mixtures acquired with unit resolving power mass spectrometers and EI sources into the sum of the single components fragmentation patterns weighted by their relative concentration using interior‐point least‐square fitting. To fit compounds with poorly known fragmentation patterns, we used a Monte‐Carlo method to vary the intensity of individual fragment ions. We then decomposed the spectrum thousands of times to obtain a statistical distribution. Results: By performing the deconvolution on a mixture of seven different molecules with interfering fragmentation patterns (H2O, O2, CH4, Ar, N2, C2H4, and C2H6) we show that this approach retrieves the mixing ratio of the individual components more accurately than regular mass spectra decomposition methods that rely on fragmentation patterns from general databases. It also provides the probability density function for each species's mixing ratio. Conclusions: By removing the solution degeneracy in the decomposition of mass spectra, the method described herein could significantly increase the scientific retrieval from archived space flight mass spectrometry data, where calibration of the ionization source is no longer an option.

Thomas Gautier↗

Eclipsing binary stars as tests of gravity theories - The apsidal motion of AS Camelopardalis

AS Camelopardalis is an 8th-magnitude eclipsing binary that consists of two main-sequence (B8 V and a B9.5 V) components in an eccentric orbit (e = 0.17) with an orbital period of 3.43 days. Like the eccentric eclipsing system DI Herculis, and a few other systems, AS Cam is an important test case for studying relativistic apsidal motion. In these systems, the theoretical general relativistic apsidal motion is comparable to that expected from classical effects arising from tidal and rotational deformation of the stellar components. Accurate determinations of the orbital and stellar properties of AS Cam have been made by Hilditch (1972) and Khalliulin and Kozyreva (1983) that permit the theoretical relativistic and classical contributions to the apsidal motion to be determined reasonably well. All the published timings of primary and secondary minima have been gathered and supplemented with eclipse timings from 1899 to 1920 obtained from the Harvard plate collection. Least-squares solutions of the eclipse timings extending over an 80 yr interval yield a smaller than expected apsidal motion, in agreement with that found by Khalliulin and Kozyreva from a smaller set of data. The observed apsidal motion for AS Cam is about one-third that expected from the combined relativistic and classical effects. Thus, AS Cam joins DI Her in having an observed apsidal motion significantly less than that predicted from theory.

Maloney, Frank P.↗

A dynamic melting model for the origin of Apollo 15 olivine-normative and quartz-normative mare basalts

Early studies of mare basalts from the Apollo 15 site established that two distinct groups are represented: the olivine-normative basalts (ONB) and the quartz-normative basalts (QNB). The ONB and QNB suites are distinguished petrographically by their phenocryst assemblages (the ONB's are olivine-phyric, the QNB's are generally pyroxene-phyric) and chemically by their major element compositions: the QNB's are higher in SiO2 and MgO/FeO, and lower in FeO and TiO2 than ONB's with similar MgO contents. Experimental data show that the QNB suite is derived from a more magnesian, olivine-normative parent magma, a conclusion which is supported by the recent discovery of high-SiO2 olivine-normative basalt clasts in breccia 15498. The high-SiO2 ONB's fall on olivine control lines with primitive QNB's, and least-squares mixing calculations are consistent with the high-SiO2 ONB's being parental to the more evolved QNB suite. These high-SiO2 ONB's are included as part of the 'QNB suite'. Our major element modeling results also are consistent with the conclusions of earlier studies which showed that the ONB and QNB suites cannot be related to one another by low pressure crystal fractionation. The combination of high Mg#, high SiO2, and low TiO2 in the QNB suite precludes a relationship to the ONB suite by simple removal of liquidus minerals (olivine and pigeonite). Despite these significant differences in petrography and major element composition, both groups have nearly identical trace element concentrations and chondrite-normalized abundance patterns. The major question to be addressed by any petrogenetic model for Apollo 15 mare basalts is how to form mare basalt suites with distinctly different major element characteristics but nearly identical trace element compositions. The similarity in trace element concentrations imply compositionally similar source regions and similar percent melting, but these conclusions are not easily reconciled with the observed differences in major element compositions, which require sources with distinct mineralogies or large differences in percent melt.

Vetter, Scott K.↗

Determination of Local Experimental Heat-Transfer Coefficients on Combustion Side of an Ammonia-Oxygen Rocket

Local experimental heat-transfer coefficients were measured in the chamber and throat of a 2400-pound-thrust ammonia-oxygen rocket engine with a nominal chamber pressure of 600 pounds per square inch absolute. Three injector configurations were used. The rocket engine was run over a range of oxidant-fuel ratio and chamber pressure. The injector that achieved the best performance also produced the highest rates of heat flux at design conditions. The heat-transfer data from the best-performing injector agreed well with the simplified equation developed by Bartz at the throat region. A large spread of data was observed for the chamber. This spread was attributed generally to the variations of combustion processes. The spread was least evident, however, with the best-performing injector.

Liebert, Curt H.↗

Comparison of some biased estimation methods (including ordinary subset regression) in the linear model

Ridge, Marquardt's generalized inverse, shrunken, and principal components estimators are discussed in terms of the objectives of point estimation of parameters, estimation of the predictive regression function, and hypothesis testing. It is found that as the normal equations approach singularity, more consideration must be given to estimable functions of the parameters as opposed to estimation of the full parameter vector; that biased estimators all introduce constraints on the parameter space; that adoption of mean squared error as a criterion of goodness should be independent of the degree of singularity; and that ordinary least-squares subset regression is the best overall method.

Sidik, S. M.↗

Solutions to Kuessner's integral equation in unsteady flow using local basis functions

The computational procedure and numerical results are presented for a new method to solve Kuessner's integral equation in the case of subsonic compressible flow about harmonically oscillating planar surfaces with controls. Kuessner's equation is a linear transformation from pressure to normalwash. The unknown pressure is expanded in terms of prescribed basis functions and the unknown basis function coefficients are determined in the usual manner by satisfying the given normalwash distribution either collocationally or in the complex least squares sense. The present method of solution differs from previous ones in that the basis functions are defined in a continuous fashion over a relatively small portion of the aerodynamic surface and are zero elsewhere. This method, termed the local basis function method, combines the smoothness and accuracy of distribution methods with the simplicity and versatility of panel methods. Predictions by the local basis function method for unsteady flow are shown to be in excellent agreement with other methods. Also, potential improvements to the present method and extensions to more general classes of solutions are discussed.

Fromme, J. A.↗

Crustal interpretation of the MAGSAT data in the continental United States

The processing of MAGSAT scalar data to construct a crustal magnetic anomaly map over the continental U.S. involves removal of the reference field model, a path-by-path subtraction of a low order polynomial through a least-squares fit to reduce orbital offset errors, and a two dimensional spectral filtering to mitigate the spectral bias induced by the path-by-path orbital correction scheme. The resultant anomaly map shows reasonably good correlations with an aeromagnetic map derived from the project MAGNET. Prominent satellite magnetic anomalies are identified in terms of geological provinces and age boundaries. An inversion method was applied to MAGSAT data which produces both the Curie depth topography and laterally varying magnetic susceptibility of the crust. A contoured Curie depth map thus derived shows general agreements with a crustal thickness map based on seismic data.

Won, I. J.↗

On the estimation of the correlation dimension and its application to radar reflector discrimination

Recently, system theorists have recognized that low order systems of nonlinear differential equations can give rise to solutions which are neither periodic, constant, nor predictable in steady state, but which are nonetheless bounded and deterministic. This behavior, which was first described in the study of weather systems, has been termed 'chaotic.' Much study of chaotic systems has concentrated on analysis of the systems' phase space attractors. It has been recognized that invariant measures of the attractor possess inherent information about the system. One such measure is the dimension of the attractors. The dimension of a chaotic attractor has been shown to be noninteger, leading to the term 'strange attractor;' the attractor is said to have a fractal structure. The correlation dimension has become one of the most popular measures of dimension. However, many problems have been identified in correlation dimension estimation from time sequences. The most common methods for obtaining the correlation dimension have been least squares curves fitting to find the slope of the correlation integral and the Takens Estimator. However, these estimates show unacceptable sensitivity to the upper limit on the distance chosen. Here, a new method is proposed which is shown to be rather insensitive to the upper limit and to perform in a very stable manner, at least in the absence of noise. The correlation dimension is also shown to be an effective discriminant in distinguishing between radar returns resulting from weather and those from the ground. The weather returns are shown to have a correlation dimension generally between 2.0 and 3.0, while ground returns have a correlation dimension exceeding 3.0.

Barnett, Kevin D.↗

Computer user's manual for a generalized curve fit and plotting program

A FORTRAN coded program has been developed for generating plotted output graphs on 8-1/2 by 11-inch paper. The program is designed to be used by engineers, scientists, and non-programming personnel on any IBM 1130 system that includes a 1627 plotter. The program has been written to provide a fast and efficient method of displaying plotted data without having to generate any additions. Various output options are available to the program user for displaying data in four different types of formatted plots. These options include discrete linear, continuous, and histogram graphical outputs. The manual contains information about the use and operation of this program. A mathematical description of the least squares goodness of fit test is presented. A program listing is also included.

Schlagheck, R. A.↗

An interactive graphics program to retrieve, display, compare, manipulate, curve fit, difference and cross plot wind tunnel data

The Aerodynamic Data Analysis and Integration System (ADAIS), developed as a highly interactive computer graphics program capable of manipulating large quantities of data such that addressable elements of a data base can be called up for graphic display, compared, curve fit, stored, retrieved, differenced, etc., was described. The general nature of the system is evidenced by the fact that limited usage has already occurred with data bases consisting of thermodynamic, basic loads, and flight dynamics data. Productivity using ADAIS of five times that for conventional manual methods of wind tunnel data analysis is routinely achieved. In wind tunnel data analysis, data from one or more runs of a particular test may be called up and displayed along with data from one or more runs of a different test. Curves may be faired through the data points by any of four methods, including cubic spline and least squares polynomial fit up to seventh order.

Elliott, R. D.↗

Data Processing of Miniaturized Laser Heterodyne Radiometer (mini-LHR) Ground Instrument Retrievals

The Miniaturized Laser Heterodyne Radiometer (mini-LHR) is a passive ground instrument that observes the mole fraction of carbon dioxide (CO) and methane (CH)in the atmospheric column by measuring their absorption of sunlight at 1.6 microns. A laser heterodyne radiometer is similar in design to the super heterodyne radio receiver that is well known by ham radio enthusiasts. While not previously a commercial technique, laser heterodyne radiometers have a history of measuring atmospheric trace gases that started in the 1960s. With the commercial availability of inexpensive, low-power, thumbnail-sized lasers developed for the telecommunications industry, it was possible to miniaturize this technique and ultimately commercialize it. The mini-LHR has been under development at NASA GSFC since 2009. During that time, in addition to signicant technical improvements, processing has also evolved and been streamlined. Here we present details of the processing approach for raw mini-LHR data to produce 30- and 60-minute data products of CH and CO column mole fractions. Processing occurs in two general stages: a python-based pre-processing of raw data, followed by ingestion into a Planetary Spectrum Generator (PSG) retrieval algorithm. Raw data processing involves removal of outliers, correcting for changes in air mass throughout the day, averaging scans, and ultimately converting averaged scans into transmittance vs. wavelength. The PSG retrieval simulates a spectra for the time/day/location of the scan with meteorological inputs from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data set and then perturbs the concentrations of CO and CH to obtain a t based on an iterative least-squares curve fitting procedure.

Giancarlo Roberto Zambrano↗

General purpose algorithms for characterization of slow and fast phase nystagmus

In the overall aim for a better understanding of the vestibular and optokinetic systems and their roles in space motion sickness, the eye movement responses to various dynamic stimuli are measured. The vestibulo-ocular reflex (VOR) and the optokinetic response, as the eye movement responses are known, consist of slow phase and fast phase nystagmus. The specific objective is to develop software programs necessary to characterize the vestibulo-ocular and optokinetic responses by distinguishing between the two phases of nystagmus. The overall program is to handle large volumes of highly variable data with minimum operator interaction. The programs include digital filters, differentiation, identification of fast phases, and reconstruction of the slow phase with a least squares fit such that sinusoidal or psuedorandom data may be processed with accurate results. The resultant waveform, slow phase velocity eye movements, serves as input data to the spectral analysis programs previously developed for NASA to analyze nystagmus responses to pseudorandom angular velocity inputs.

Lessard, Charles S.↗

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗

Improved solution accuracy for Landsat-4 (TDRSS-user) orbit determination

This paper presents the results of a study to compare the orbit determination accuracy for a Tracking and Data Relay Satellite System (TDRSS) user spacecraft, Landsat-4, obtained using a Prototype Filter Smoother (PFS), with the accuracy of an established batch-least-squares system, the Goddard Trajectory Determination System (GTDS). The results of Landsat-4 orbit determination will provide useful experience for the Earth Observing System (EOS) series of satellites. The Landsat-4 ephemerides were estimated for the January 17-23, 1991, timeframe, during which intensive TDRSS tracking data for Landsat-4 were available. Independent assessments were made of the consistencies (overlap comparisons for the batch case and convariances for the sequential case) of solutions produced by the batch and sequential methods. The filtered and smoothed PFS orbit solutions were compared with the definitive GTDS orbit solutions for Landsat-4; the solution differences were generally less than 15 meters.

Oza, D. H.↗

Geoid undulations and gravity anomalies over the Aral Sea, the Black Sea and the Caspian Sea from a combined GEOS-3/SEASAT/GEOSAT altimeter data set

Satellite-based altimetric data taken by GOES-3, SEASAT, and GEOSAT over the Aral Sea, the Black Sea, and the Caspian Sea are analyzed and a least squares collocation technique is used to predict the geoid undulations on a 0.25x0.25 deg. grid and to transform these geoid undulations to free air gravity anomalies. Rapp's 180x180 geopotential model is used as the reference surface for the collocation procedure. The result of geoid to gravity transformation is, however, sensitive to the information content of the reference geopotential model used. For example, considerable detailed surface gravity data were incorporated into the reference model over the Black Sea, resulting in a reference model with significant information content at short wavelengths. Thus, estimation of short wavelength gravity anomalies from gridded geoid heights is generally reliable over regions such as the Black Sea, using the conventional collocation technique with local empirical covariance functions. Over regions such as the Caspian Sea, where detailed surface data are generally not incorporated into the reference model, unconventional techniques are needed to obtain reliable gravity anomalies. Based on the predicted gravity anomalies over these inland seas, speculative tectonic structures are identified and geophysical processes are inferred.

Au, Andrew Y.↗

ELAS: A general-purpose computer program for the equilibrium problems of linear structures. Volume 2: Documentation of the program

A general purpose digital computer program for the in-core solution of linear equilibrium problems of structural mechanics is documented. The program requires minimum input for the description of the problem. The solution is obtained by means of the displacement method and the finite element technique. Almost any geometry and structure may be handled because of the availability of linear, triangular, quadrilateral, tetrahedral, hexahedral, conical, triangular torus, and quadrilateral torus elements. The assumption of piecewise linear deflection distribution insures monotonic convergence of the deflections from the stiffer side with decreasing mesh size. The stresses are provided by the best-fit strain tensors in the least squares at the mesh points where the deflections are given. The selection of local coordinate systems whenever necessary is automatic. The core memory is used by means of dynamic memory allocation, an optional mesh-point relabelling scheme and imposition of the boundary conditions during the assembly time.

Utku, S.↗

JPL pulsar timing observations. I - The Vela pulsar

It is found through an analysis of 12 years of arrival time data for PSR 0833-45, by fitting a model of the pulse period to the data in a least-squares sense, that there exist regions of exceptionally stable behavior in which rms residuals are typically about 40 microsec. In addition, each residual exhibits a significant period parameter second derivative which is found to change little within a given interjump era. Parameter changes in the period P, and the derivatives introduced by each of the four large observed jumps in period, are recomputed to demonstrate that the post-jump decay is nearly exponential in form. It is found that random fluctuations in P are an order of magnitude less than the systematic variations caused by the seond derivative period parameter, whose component associated with the general slowdown of the pulsar can be reliably measured when the observations span 55 years or more.

Downs, G. S.↗

Optimum data weighting and error calibration for estimation of gravitational parameters

A new technique was developed for the weighting of data from satellite tracking systems in order to obtain an optimum least-squares solution and an error calibration for the solution parameters. Data sets from optical, electronic, and laser systems on 17 satellites in GEM-T1 Goddard Earth Model-T1 (GEM-T1) were employed toward application of this technique for gravity field parameters. Also GEM-T2 (31 satellites) was recently computed as a direct application of the method and is summarized. The method employs subset solutions of the data associated with the complete solution to agree with their error estimates. With the adjusted weights the process provides for an automatic calibration of the error estimates for the solution parameters. The data weights derived are generally much smaller than corresponding weights obtained from nominal values of observation accuracy or residuals. Independent tests show significant improvement for solutions with optimal weighting. The technique is general and may be applied to orbit parameters, station coordinates, or other parameters than the gravity model.

Lerch, Francis J.↗