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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2007

The Local Stellar Velocity Field via Vector Spherical Harmonics

We analyze the local field of stellar tangential velocities for a sample of 42,339 nonbinary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism. We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships, we determine the solar velocity with respect to the local stars of (V(sub X), V(sub Y), V(sub Z)) (10.5, 18.5, 7.3) +/- 0.1 km s(exp -1) not corrected for the asymmetric drift with respect to the local standard of rest. If only stars more distant than 100 pc are considered, the peculiar solar motion is (V(sub X), V(sub Y), V(sub Z)) (9.9, 15.6, 6.9) +/- 0.2 km s(exp -1). The adverse effects of harmonic leakage, which occurs between the reflex solar motion represented by the three electric vector harmonics in the velocity space and higher degree harmonics in the proper-motion space, are eliminated in our analysis by direct subtraction of the reflex solar velocity in its tangential components for each star. The Oort parameters determined by a straightforward least-squares adjustment in vector spherical harmonics are A=14.0 +/- 1.4, B=13.1 +/- 1.2, K=1.1 +/- 1.8, and C=2.9 +/- 1.4 km s(exp -1) kpc(exp -1). The physical meaning and the implications of these parameters are discussed in the framework of a general linear model of the velocity field. We find a few statistically significant higher degree harmonic terms that do not correspond to any parameters in the classical linear model. One of them, a third-degree electric harmonic, is tentatively explained as the response to a negative linear gradient of rotation velocity with distance from the Galactic plane, which we estimate at approximately -20 km s(exp -1) kpc(exp -1). A similar vertical gradient of rotation velocity has been detected for more distant stars representing the thick disk (z greater than 1 kpc), but here we surmise its existence in the thin disk at z less than 200 pc. The most unexpected and unexplained term within the Ogorodnikov-Milne model is the first-degree magnetic harmonic, representing a rigid rotation of the stellar field about the axis -Y pointing opposite to the direction of rotation. This harmonic comes out with a statistically robust coefficient of 6.2 +/- 0.9 km s(exp -1) kpc(exp -1) and is also present in the velocity field of more distant stars. The ensuing upward vertical motion of stars in the general direction of the Galactic center and the downward motion in the anticenter direction are opposite to the vector field expected from the stationary Galactic warp model.

Galaxy↗

NASA NTRS1987

Potential surfaces for O atom-polymer reactions

Ab initio quantum chemistry methods are used to study the energetics of interactions of O atoms with organic compounds. Polyethylene (CH2)n has been chosen as the model system to study the interactions of O(3P) and O(1D) atoms with polymers. In particular, H abstraction is investigated and polyethylene is represented by a C3 (propane) oligomeric model. The gradient method, as implemented in the GRADSCF package of programs, is used to determine the geometries and energies of products and reactants. The saddle point, barrier geometry is determined by minimizing the squares of the gradients of the potential with respect to the internal coordinates. To correctly describe the change in bonding during the reaction at least a two configuration MCSCF (multiconfiguration self consistent field) or GVB (generalized valence bond) wave function has to be used. Basis sets include standard Pople and Dunning sets, however, increased with polarization functions and diffuse p functions on both the C and O atoms. The latter is important due to the O(-) character of the wave function at the saddle point and products. Normal modes and vibrational energy levels are given for the reactants, saddle points and products. Finally, quantitative energetics are obtained by implementing a small CAS (complete active space) approach followed by limited configuration interaction (CI) calculations. Comparisons are made with available experimental data.

Laskowski, B. C.↗

NASA NTRS1991

Linear and circular digital spectral analysis of serial data

Two methods of digital spectral analysis of unevenly sampled data are developed and illustrated here. One method uses a linear function of time (or space), the other uses circular functions. The circular method turns out to be essentially equivalent to a least-squares sine-wave analysis. The linear, anharmonic method uses only the field of real numbers and elementary algebraic operations, and hence it can be made computationally very fast and accurate. Both methods are very general, properly handling all kinds of time series ranging from simple series consisting only of the times of events to complicated series consisting of pulses with long duty cycles. The two methods are here applied to the analysis of annual mean relative sunspot numbers.

Stothers, Richard B.↗

NASA NTRS1996

Polarization Signals of Common Spacecraft Materials

This is the final report documenting the results of the polarization testing of near-planar objects with various reflectance properties. The purpose of this investigation was to determine the portion of the reflected signal which is polarized for materials commonly used in space applications. Tests were conducted on several samples, with surface characteristics ranging from highly reflective to relatively dark. The measurements were obtained by suspending the test object in a beam of collimated light. The amount of light falling on the sample was controlled by a circular aperture placed in the light field. The polarized reflectance at various phase angles was then measured. A nonlinear least squares fitting program was used for analysis. For the specular test objects, the reflected signals were measured in one degree increments near the specular point. Otherwise, measurements were taken every five degrees in phase angle. Generally, the more diffuse surfaces had lower polarized reflectances than their more specular counterparts. The reflected signals for the more diffuse surfaces were spread over a larger phase angle range, while the signals from the more specular samples were reflected almost entirely within five degrees of angular deviation from the specular point. The method used to test all the surfaces is presented. The results of this study will be used to support the NASA Orbital Debris Optical Signature Tests. These tests are intended to help better understand the reflectance properties of materials often used in space applications. This data will then be used to improve the capabilities for identification and tracking of space debris.

Gravseth, Ian↗

NASA NTRS1985

Results of the May 1982 ATMAP Campaign at the Meteor Wind Station Kuehlungsborn (54 Deg N, 12 Deg E), GDR

A sinusoisal least squares fit over the entire data of 1 to 6 May 1982 was computed. The results are shown. Wind results refere to an average height of 95 km (no height determination). The wind results correspond, in general, to the average behavior of early May. The zonal prevailing west wind, however, is stronger than the term mean value (-5 m/s). Also, the amplitude of the semidiurnal idal wind is larger (mean value 1977 to 1983: 8 m/s).

Greisiger, K. M.↗

NASA NTRS1979

State and model error estimation for distributed parameter systems

In-flight estimation of large structure model errors in order to detect inevitable deficiencies in large structure controller/estimator models is discussed. Such an estimation process is particularly applicable in the area of shape control system design required to maintain a prescribed static structural shape and, in addition, suppress dynamic disturbances due to the vehicle vibrational modes. The paper outlines a solution to the problem of static shape estimation where the vehicle shape must be reconstructed from a set of measurements discretely located throughout the structure. The estimation process is based on the principle of least-squares that inherently contains the definition and explicit computation of model error estimates that are optimal in some sense. Consequently, a solution is provided for the problem of estimation of static model errors (e.g., external loads). A generalized formulation applicable to distributed parameters systems is first worked out and then applied to a one-dimensional beam-like structural configuration.

Rodriguez, G.↗

NASA NTRS1992

Gross-fragmentation of meteoroids and bulk density of Geminids from photographic fireball records

The explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid moving in any actual atmosphere was published several years ago. The solution yields the theoretical relation of l, the distance flown by the meteoroid in its trajectory, as a function of time, t, assuming that the height, h, is a known function of l. The photographic records of meteors and fireballs are coded by time marks, using a rotating shutter or a similar device to break the moving image. Time is, thus, the independent variable and for each time mark on a meteoroid trajector, the observed distance along the trajectory, l sub obs, as well as the double- or multiple- station photographs of the same meteoroid. Applying this solution to all available Prairie Network (PN) fireball-records, we recognized that the majority of them gave good solutions with standard deviations somewhat bigger than the intrinsic geometrical precision of the data. We also noticed that, on an average, previous methods of evaluation of the meteoroid velocities (interpolation polynomials, numerical differenciation of the observed l sub obs) used up to only several tens of percent of the intrinsic precision of the PN observational data. When residuals of these solutions, i.e. l sub obs - l sub com, were represented as a function of time for about 75 percent of solutions. The remaining 25 percent of residuals showed systematic changes with time exceeding one standard deviation. We tried to explain these systematic time course of residuals by using different meteoroids first computed theoretically and then analyzed by the same model as the natural PN fireballs were. The conclusion of these model computations is that systematic time changes of residuals in the nonfragmenting model exceeding one standard of deviation are caused by sudden gross fragmentation at one or more trajectory points. Thus, we generalized the explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid by allowing for one or more points, where sudden gross fragmentation can occur. Using this generalized solution, the distances along the meteoroid trajectory can be computed for any choice of input parameters and compared with the observed distances flown by the meteoroid. For the most precise and long fireball trajectories, the least-squares solution can, thus, yield the initial velocities, the ablation coefficients, the dynamical masses, the positions of gross-fragmentation points, and the terminal mass. At a gross-fragmentation point, the ratio of the main mass to all the remaining fragments can be compared with the dynamic mass determined from our gross-fragmentation model and, thus, the meteoroid bulk density can be evaluated. We applied the gross-fragmentation model to sever PN fireballs showing time changes of residuals, and we recognized that, in all these cases, the new computed bulk densities of meteoroids resulted higher in comparison with the meteoroid densities determined with the non-gross-fragmentation model. Other aspects of the study are discussed.

Ceplecha, Zdenek↗

NASA NTRS2000

Stress Recovery and Error Estimation for 3-D Shell Structures

The C1-continuous stress fields obtained from finite element analyses are in general lower- order accurate than are the corresponding displacement fields. Much effort has focussed on increasing their accuracy and/or their continuity, both for improved stress prediction and especially error estimation. A previous project developed a penalized, discrete least squares variational procedure that increases the accuracy and continuity of the stress field. The variational problem is solved by a post-processing, 'finite-element-type' analysis to recover a smooth, more accurate, C1-continuous stress field given the 'raw' finite element stresses. This analysis has been named the SEA/PDLS. The recovered stress field can be used in a posteriori error estimators, such as the Zienkiewicz-Zhu error estimator or equilibrium error estimators. The procedure was well-developed for the two-dimensional (plane) case involving low-order finite elements. It has been demonstrated that, if optimal finite element stresses are used for the post-processing, the recovered stress field is globally superconvergent. Extension of this work to three dimensional solids is straightforward. Attachment: Stress recovery and error estimation for shell structure (abstract only). A 4-node, shear-deformable flat shell element developed via explicit Kirchhoff constraints (abstract only). A novel four-node quadrilateral smoothing element for stress enhancement and error estimation (abstract only).

Riggs, H. R.↗

NASA NTRSDate not supplied

Real-Time State Estimation of Structural Modes for an Aeroelastic Wind Tunnel Model

A method is presented for estimating displacements, velocities, and accelerations of structural modes in generalized coordinates from measured sensor data in real time. Specifically, strain data from conventional strain gauges and fiber optic strain sensors (FOSS) were combined with the strain modes (obtained from a finite element model) in a least-squares estimator to produce structural mode displacement estimates and uncertainties. Likewise, accelerometer data were combined with displacement mode shapes in a second least-squares estimator to produce structural mode acceleration estimates. A Kalman filter was then used to refine the displacement estimates and obtain velocity estimates. The method was applied to the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. The technique was found to be useful for real-time control and system identification applications.

Jared A. Grauer↗

NASA NTRS1979

Development of a winter wheat adjustable crop calendar model

The author has identified the following significant results. Least squares techniques were applied for parameter estimation of functions to predict winter wheat phenological stage with daily maximum temperature, minimum temperature, daylength, and precipitation as independent variables. After parameter estimation, tests were conducted using independent data. It may generally be concluded that exponential functions have little advantage over polynomials. Precipitation was not found to significantly affect the fits. The Robertson triquadratic form, in general use for spring wheat, yielded good results, but special techniques and care are required. In most instances, equations with nonlinear effects were found to yield erratic results when utilized with averaged daily environmental values as independent variables.

Baker, J. R.↗

NASA NTRS1984

Fitting aerodynamic forces in the Laplace domain: An application of a nonlinear nongradient technique to multilevel constrained optimization

A technique which employs both linear and nonlinear methods in a multilevel optimization structure to best approximate generalized unsteady aerodynamic forces for arbitrary motion is described. Optimum selection of free parameters is made in a rational function approximation of the aerodynamic forces in the Laplace domain such that a best fit is obtained, in a least squares sense, to tabular data for purely oscillatory motion. The multilevel structure and the corresponding formulation of the objective models are presented which separate the reduction of the fit error into linear and nonlinear problems, thus enabling the use of linear methods where practical. Certain equality and inequality constraints that may be imposed are identified; a brief description of the nongradient, nonlinear optimizer which is used is given; and results which illustrate application of the method are presented.

Tiffany, S. H.↗

NASA NTRS1988

Disruption of giant comets in the solar system and around other stars

In a standard cometary mass distribution (dN/dM) alpha M(-a), a = 1.5 to 2.0) most of the mass resides in the largest comets. The maximum mass M sub max for which this distribution holds uncertain but there are theoretical and observational indications that M sub max is at least approx. 10(23)g. Chiron, although formally classified as an asteroid, is most likely a giant comet in this mass range. Its present orbit is unstable and it is expected to evolve into a more typical short period comet orbit on a timescale of approx. 10(6) to 10(7)yr. The breakup of a chiron-like comet of mass approx. 10(23)g could in principle produce approx. 10(5) Halley-size comets, or a distribution with an even larger number. If a giant comet was in a typical short period comet orbit, such a breakup could result in a relatively brief comet shower (duration approx. less than 10(6)yr) with some associated terrestrial impacts. However, the most significant climatic effects may not in general be due to the impacts themselves but to the greatly enhanced zodiacal dust cloud in the inner Solar System. (Although this is probably not the case for the unique K-T impact). Researchers used a least Chi square program with error analysis to confirm that the 2 to 5 micrometer excess spectrum of Giclas 29 to 38 can be adequately fitted with either a disk of small inefficient (or efficient) grains or a single temperature black body. Further monitoring of this star may allow discrimination between these two models.

Whitmire, D. P.↗

NASA NTRS1993

Multiple degree of freedom force-state component identification

Force-state mapping is extended to allow the quasistatic characterization of realistic multiple degree of freedom systems whose constitutive relations are nonlinear, dissipative, coupled and depend on memory of past states. A general constitutive force-state model is formulated, which includes dynamic hysteresis phenomena typical of that found in deployable structures. A two step parameter identification approach is developed, which employs a linear extended Kalman filter and a nonlinear least squares fit. A six degree of freedom force-state testing device was constructed, and used to obtain data on one bay of a simple truss and on one bay of the MODE deployable truss.

Masters, Brett P.↗

NASA NTRS1994

Measurement uncertainty and feasibility study of a flush airdata system for a hypersonic flight experiment

Presented is a feasibility and error analysis for a hypersonic flush airdata system on a hypersonic flight experiment (HYFLITE). HYFLITE heating loads make intrusive airdata measurement impractical. Although this analysis is specifically for the HYFLITE vehicle and trajectory, the problems analyzed are generally applicable to hypersonic vehicles. A layout of the flush-port matrix is shown. Surface pressures are related airdata parameters using a simple aerodynamic model. The model is linearized using small perturbations and inverted using nonlinear least-squares. Effects of various error sources on the overall uncertainty are evaluated using an error simulation. Error sources modeled include boundarylayer/viscous interactions, pneumatic lag, thermal transpiration in the sensor pressure tubing, misalignment in the matrix layout, thermal warping of the vehicle nose, sampling resolution, and transducer error. Using simulated pressure data for input to the estimation algorithm, effects caused by various error sources are analyzed by comparing estimator outputs with the original trajectory. To obtain ensemble averages the simulation is run repeatedly and output statistics are compiled. Output errors resulting from the various error sources are presented as a function of Mach number. Final uncertainties with all modeled error sources included are presented as a function of Mach number.

Whitmore, Stephen A.↗

NASA NTRS1992

A variational method for finite element stress recovery: Applications in one-dimension

It is well-known that stresses (and strains) calculated by a displacement-based finite element analysis are generally not as accurate as the displacements. In addition, the calculated stress field is typically discontinuous at element interfaces. Because the stresses are typically of more interest than the displacements, several procedures have been proposed to obtain a smooth stress field, given the finite element stresses, and to improve the accuracy. Hinton and Irons introduced global least squares smoothing of discrete data defined on a plane using a finite element formulation. Tessler and co-workers recently developed a conceptually similar formulation for smoothing of two-dimensional data based on a discrete least square approximation with a penalty constraint. The penalty constraint results in a stress field which is C(exp 1)-continuous, a result not previously obtained. The approach requires additional, 'smoothing' finite element analysis and for their two-dimensional application, they used a conforming C(exp 0)-continuous triangular finite element based on a conforming plate element. This paper presents the results of a detailed investigation into the application of Tessler's smoothing procedure to the smoothing of finite element stresses from one-dimensional problems. Although the one-dimensional formulation has some practical applicability, such as in truss, beam, axisymmetric mechanics, and one-dimensional heat conduction, the primary motivation for developing the one-dimensional smoothing case is to explore the characteristics of the general smoothing strategy. In particular, it is used to describe the behavior of the method and to explore the suitability of criteria proposed for the smoothing analysis. Prior to presenting numerical results, the variational formulation of the smoothing strategy is presented and a criterion for the smoothing analysis is described.

Riggs, H. Ronald↗

NASA NTRS1987

Waveform analysis for Geosat day 96

A detailed model waveform function has been least-squares-fitted to 10-second averages of waveform sampler data from selected over-ocean portions of day 96 of the Geosat mission. The results confirm that the height corrections already applied in the routine Geosat data processing are generally good to well within 10 centimeters. Additional height corrections are provided that can refine height estimates to within the several-centimeter level.

Hayne, George S.↗

NASA NTRS1990

System identification and controller design using experimental frequency response data

Recent findings from modeling and controller design for the NASA-Marshall Single Structure Control Facility have raised questions regarding the ability of modern control design techniques and modern modeling techniques to deal effectively with the stringent modeling and control design requirements of Large Space Structure Control. A brief and general discussion is presented of the results of studies into the modeling and control issues performed under sponsorship of the NASA/ASEE Summer Faculty Fellowship Program. Several issues are addressed. The first is a study of a modeling technique based on least squares identification of individual transfer functions from measured frequency response data. The second is a study of multiobjective optimization techniques applied to the modeling, or system identification, problem. The third issue is a study into the question of whether multiobjective optimization approaches can be effectively used for control system design using only frequency response data, thereby bypassing the difficult modeling problem. The last issue studied involves the resolution of seeming discrepancies between predicted and measured control computer time delays in the Single Structure Control Facility.

Irwin, R. Dennis↗

NASA NTRS1992

Ballistic limit curve regression for Freedom Station orbital debris shields

A procedure utilized at Marshall Space Flight Center to formulate ballistic limit curves for the Space Station Freedom's manned module orbital debris shields is presented. A stepwise linear least squares regression method similar to that employed by Burch (1967) is used to relate a penetration parameter to various projectile and target descriptors. A stepwise regression was also conducted with the model reduced to lower forms, thus eliminating the effects of generalized assumptions.

Jolly, William H.↗