Search NASA⌕ Search

SEARCH · Search NASA

Results for “least squares solution”

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 235 records · Page 13

In-plane, flexural, twisting and thickness-shear coefficients for stiffness and damping of a monolayer filamentary composite, part 1

Elastic and damping analyses resulting in determinations of the various stiffnesses and associated loss tangents for the complete characterization of the elastic and damping behavior of a monofilament composite layer are presented. For the determination of the various stiffnesses, either an elementary mechanics-of-materials formulation or a more rigorous mixed-boundary-value elasticity formulation is used. The solution for the latter formulation is obtained by means of the boundary-point least-square error technique. Kimball-Lovell type damping is assumed for each of the constituent materials. For determining the loss tangents associated with the various stiffnesses, either the viscoelastic correspondence principle or an energy analysis based on the appropriate elastic stress distribution is used.

Bert, C. W.↗

Application of advanced filtering methods to the determination of the interplanetary orbit of Mariner '71.

This paper presents the results of the applications of advanced filtering methods to the determination of the interplanetary orbit of the Mariner '71 spacecraft. The advanced techniques are specific extensions of the Kalman filter. The special problems associated with applying these techniques are discussed and the particular algorithmic implementations are outlined. The advanced methods are compared against the weighted least squares filters of conventional application. The results reveal that relatively simple advanced filter configurations yield solutions superior to those of the conventional methods when applied to the Mariner '71 radio measurements.

Rourke, K. H.↗

Aerodynamic influence coefficient method using singularity splines

A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.

Mercer, J. E.↗

Geodetic linear estimation technique and the norm choice problem

In this work the mathematical and probabilistic background of standard linear estimation techniques used in geodesy is clarified, and their interrelationship is revealed with the help of best approximation theory and the normal equations. Emphasis is given to the separation of the deterministic solution to the approximation problem from the probabilistic justification of the metric of the approximation. Least squares prediction has been related to deterministic (exact) collocation, and minimum error bound has been identified as a prediction optimality criterion in the latter. Criteria for the optimal choice of norm in Hilbert space collocation are proposed for gravimetric geodesy problems.

Dermanis, A.↗

Spherical harmonic expansion of the Levitus Sea surface topography

Prior information for the stationary sea surface topography (SST) may be needed in altimetric solutions that intend to simultaneously improve the gravity field and determine the SST. For this purpose the oceanographically derived SST estimates are represented by a spherical harmonic expansion. The spherical harmonic coefficients are computed from a least squares adjustment of the data covering the majority of the oceanic regions of the world. Several tests are made to determine the optimum maximum degree of solution and the best configuration of the geometry of the data in order to obtain a solution that fits the data and also provides a good spectral representation of the SST.

Engelis, Theodossios↗

Observations of the 8 December 1987 occultation of AG+40 deg 0783 by 324 Bamberga

The occultation of AG+40 deg 0783 by 324 Bamberga on 8 December 1987 was observed at 13 sites in the United States, Japan, and China. At four sites the event was observed photoelectrically; the other observations were visual. A least-squares fit of a circular limb profile to the data gives a diameter of 227.6 + or - 1.9 km. However, this solution is inconsistent with a negative visual observation near the northern edge of the ground track. The inconsistency cannot be removed by assuming an elliptical profile. The data suggest that Bamberga, despite its low-amplitude lightcurve, may depart significantly from a spherical or ellipsoidal shape. The asteroid also appears to be at least 10 percent smaller than indicated by infrared radiometry.

Millis, R. L.↗

Observations of the 8 December 1987 occultation of AG +40 deg 0783 by 324 Bamberga

The occultation of AG+40 deg 0783 by 324 Bamberga on 8 December 1987 was observed at 13 sites in the United States, Japan, and China. At four sites the event was observed photoelectrically; the other observations were visual. A least-squares fit of a circular limb profile to the data gives a diameter of 227.6 + or - 1.9 km. However, this solution is inconsistent with a negative visual observation near the northern edge of the ground track. The inconsistency cannot be removed by assuming an elliptical profile. The data suggest that Bamberga, despite its low-amplitude lightcurve, may depart significantly from a spherical or ellipsoidal shape. The asteroid also appears to be at least 10 percent smaller than indicated by infrared radiometry.

Millis, R. L.↗

The L sub 1 finite element method for pure convection problems

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Jiang, Bo-Nan↗

The L sub 1 finite element method for pure convection problems

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Jiang, Bo-Nan↗

Detailed gravity anomalies from Geos 3 satellite altimetry data

Detailed gravity anomalies are computed from a combination of Geos 3 satellite altimeter and terrestrial gravity data using least-squares principles. The mathematical model used is based on the Stokes' equation modified for a nonglobal solution. Using Geos 3 data in the calibration area, the effects of several anomaly parameter configurations and data densities/distributions on the anomalies and their accuracy estimates are studied. The accuracy estimates for 1 deg x 1 deg mean anomalies from low density altimetry data are of the order of 4 mgal. Comparison of these anomalies with the terrestrial data and also with Rapp's data derived using collocation techniques show rms differences of 7.2 and 4.9 mgal, respectively. Indications are that the anomaly accuracies can be improved to about 2 mgal with high density data. Estimation of 30 in. x 30 in. mean anomalies indicates accuracies of the order of 5 mgal. Proper verification of these results will be possible only when accurate ground truth data become available.

Gopalapillai, G. S.↗

On correcting radial orbit errors for altimetric satellites using crossover analysis

A method is proposed for correcting radial orbit error in measurements of sea surface height using a satellite altimeter. Traditionally, the orbit error is modeled in terms of a Fourier series with the Fourier coefficients determined by minimizing the residual crossover difference in a least-squares sense. In this method an a priori constraint must be imposed to obtain a unique solution. It is shown that by using singular value decomposition, no such constraint is needed. This proposed method leaves the geographically dependent errors unchanged and makes only those corrections warranted by the information contained in crossover differences. Thus, the resultant ocean topography is free from any undue distortion which might be incurred by an a priori constraint. It is suggested that the method is useful for application to high-accuracy altimetric mission, because the orbit error can be reduced without compromising the accuracy of the measured mean ocean topography.

Fu, Lee-Lueng↗

On the sea-state bias of the Geosat altimeter

The sea-state bias in a satellite altimeter's range measurement is caused by the influence of ocean waves on the radar return pulse; it results in an estimate of sea level that is too low according to some function of the wave height. This bias is here estimated for Geosat by correlating collinear differences of altimetric sea-surface heights with collinear differences of significant wave heights (H1/3). Corrections for satellite orbit error are estimated simultaneously with the sea-state bias. Based on twenty 17-day repeat cycles of the Geosat Exact Repeat Mission, the solution for the sea-state bias is 2.6 + or - 0.2 percent of H1/3. The least-squares residuals, however, show a correlation with wind speed U, so the traditional model of the bias has been supplemented with a second term: H1/3 + alpha-2H1/3U. This second term produces a small, but statistically significant, reduction in variance of the residuals. Both systematic and random errors in H1/3 and U tend to bias the estimates of alpha-1 and alpha-2, which complicates comparisons of the results with ground-based measurements of the sea-state bias.

Ray, Richard D.↗

CAN-DO, CFD-based Aerodynamic Nozzle Design and Optimization program for supersonic/hypersonic wind tunnels

A design program is developed which incorporates a modern approach to the design of supersonic/hypersonic wind-tunnel nozzles. The approach is obtained by the coupling of computational fluid dynamics (CFD) with design optimization. The program can be used to design a 2D or axisymmetric, supersonic or hypersonic, wind-tunnel nozzles that can be modeled with a calorically perfect gas. The nozzle design is obtained by solving a nonlinear least-squares optimization problem (LSOP). The LSOP is solved using an iterative procedure which requires intermediate flowfield solutions. The nozzle flowfield is simulated by solving the Navier-Stokes equations for the subsonic and transonic flow regions and the parabolized Navier-Stokes equations for the supersonic flow regions. The advantages of this method are that the design is based on the solution of the viscous equations eliminating the need to make separate corrections to a design contour, and the flexibility of applying the procedure to different types of nozzle design problems.

Korte, John J.↗

Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes

The accuracy of the least-squares technique for gradient reconstruction on unstructured meshes is examined. While least-squares techniques produce accurate results on arbitrary isotropic unstructured meshes, serious difficulties exist for highly stretched meshes in the presence of surface curvature. In these situations, gradients are typically under-estimated by up to an order of magnitude. For vertex-based discretizations on triangular and quadrilateral meshes, and cell-centered discretizations on quadrilateral meshes, accuracy can be recovered using an inverse distance weighting in the least-squares construction. For cell-centered discretizations on triangles, both the unweighted and weighted least-squares constructions fail to provide suitable gradient estimates for highly stretched curved meshes. Good overall flow solution accuracy can be retained in spite of poor gradient estimates, due to the presence of flow alignment in exactly the same regions where the poor gradient accuracy is observed. However, the use of entropy fixes has the potential for generating large but subtle discretization errors.

Mavriplis, Dimitri J.↗

Aerodynamic influence coefficient method using singularity splines.

A new numerical formulation with computed results, is presented. This formulation combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of the loading function methods. The formulation employs a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfies all of the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise (termed 'spline'). Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral.

Mercer, J. E.↗

Sequential Least-Squares Using Orthogonal Transformations

Square root information estimation, starting from its beginnings in least-squares parameter estimation, is considered. Special attention is devoted to discussions of sensitivity and perturbation matrices, computed solutions and their formal statistics, consider-parameters and consider-covariances, and the effects of a priori statistics. The constant-parameter model is extended to include time-varying parameters and process noise, and the error analysis capabilities are generalized. Efficient and elegant smoothing results are obtained as easy consequences of the filter formulation. The value of the techniques is demonstrated by the navigation results that were obtained for the Mariner Venus-Mercury (Mariner 10) multiple-planetary space probe and for the Viking Mars space mission.

Bierman, G. J.↗

Aerodynamic coefficient identification package dynamic data accuracy determinations: Lessons learned

The errors in the dynamic data output from the Aerodynamic Coefficient Identification Packages (ACIP) flown on Shuttle flights 1, 3, 4, and 5 were determined using the output from the Inertial Measurement Units (IMU). A weighted least-squares batch algorithm was empolyed. Using an averaging technique, signal detection was enhanced; this allowed improved calibration solutions. Global errors as large as 0.04 deg/sec for the ACIP gyros, 30 mg for linear accelerometers, and 0.5 deg/sec squared in the angular accelerometer channels were detected and removed with a combination is bias, scale factor, misalignment, and g-sensitive calibration constants. No attempt was made to minimize local ACIP dynamic data deviations representing sensed high-frequency vibration or instrument noise. Resulting 1sigma calibrated ACIP global accuracies were within 0.003 eg/sec, 1.0 mg, and 0.05 deg/sec squared for the gyros, linear accelerometers, and angular accelerometers, respectively.

Heck, M. L.↗

Polar motion from laser range measurements of GEOS-3

Using two-day arcs of GEOS-3 laser data, simultaneous solutions for pole position components, x sub p and y sub p, and orbit elements have been obtained for the period spanning 3 February to 6 March 1976 using three NASA Goddard Space Flight Center laser stations located near Washington, D.C. (STALAS) and on the islands of Bermuda and Grand Turk. The results are in general agreement with the BIH results. However, because of the locations of the laser sites, the x sub p solution is weaker than the y sub p solution. The x sub p and y sub p estimates were smoothed with a straight line by weighted least squares using the variance associated with the pole estimates as weights in order to reflect the effect of widely different data distributions. The smoothed y sub p differs by one meter with respect to the BIH smoothed values and the smoothed x sub p differs by about two meters. Spectral analysis of the results has identified frequencies associated with the orbital motion indicating the need for further improvements in the model of the physical system.

Schutz, B. E.↗