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

Investigation of Point Doppler Velocimetry (PDV) for Transition Detection in Boundary Layers

A two-component Point Doppler Velocimeter (PDV) system has been improved through the use of vapor-limited iodine cells that have responses that are insensitive to temperature variations. Two-component PDV velocity measurements have been obtained for a 1 inch diameter uniform circular jet flow at a nominal exit velocity of 60 m/sec, corresponding to a Reynolds number of 100,000. Similar data have also been obtained for an annular jet and a swirling jet, These PDV data runs have been duplicated to judge the repeatability of these measurements, and also have been compared with hot wire anemometer data for the same flow conditions. PDV mean velocity results are repeatable to within approximately 1-2 meters per second; the PDV RMS velocity results are also quite repeatable. Exit profiles of PDV mean axial velocity data generally agree with hot wire anemometer results to within about 2 meters per second as well. However, the PDV RMS velocity results are consistently lower than the hot wire results everywhere but at the exit of the standard jet, where they are too high relative to the hot wire data. This is believed to at least be partially due to the method used to compute the RMS (Root Mean Square).

Kuhlman, John↗

Least-work analysis of the problem of shear lag in box beams

The distribution of stress in the cover sheets of thin-wall box beams is analyzed, with regard to the effect of shear deformation in the cover sheets, by the method of least work. Explicit results are obtained for a number of representative cases that show the influence of the following factors on the stress patterns. (1) Variation of stress in spanwise direction as given by elementary beam theory. (2) Value of a parameter called shear-lag aspect ratio which designates the product of span-width ratio of the beam and of the square root of the ratio of effective shear modulus and tensing modulus of the cover sheets. (3) Value of ratio of cover-sheet stiffness to side-web stiffness. (4) Variation of beam height in span direction. (5) Variation of beam width in span direction. (6) Variation of cover-sheet thickness in span direction. General conclusions are drawn from the results obtained. Among them the most important one appears to be the fact that the shear-lag effect depends primarily on the following tow quantities: (1) the value of the shear-lag aspect ratio. (2) the shape of the curve representing the product of the stress of elementary beam theory and of the cover-sheet thickness.

STRESS CALCULATIONS, LEAST-WORK↗

Colored noise effects on batch attitude accuracy estimates

The effects of colored noise on the accuracy of batch least squares parameter estimates with applications to attitude determination cases are investigated. The standard approaches used for estimating the accuracy of a computed attitude commonly assume uncorrelated (white) measurement noise, while in actual flight experience measurement noise often contains significant time correlations and thus is colored. For example, horizon scanner measurements from low Earth orbit were observed to show correlations over many minutes in response to large scale atmospheric phenomena. A general approach to the analysis of the effects of colored noise is investigated, and interpretation of the resulting equations provides insight into the effects of any particular noise color and the worst case noise coloring for any particular parameter estimate. It is shown that for certain cases, the effects of relatively short term correlations can be accommodated by a simple correction factor. The errors in the predicted accuracy assuming white noise and the reduced accuracy due to the suboptimal nature of estimators that do not take into account the noise color characteristics are discussed. The appearance of a variety of sample noise color characteristics are demonstrated through simulation, and their effects are discussed for sample estimation cases. Based on the analysis, options for dealing with the effects of colored noise are discussed.

Bilanow, Stephen↗

An Optimal Orthogonal Decomposition Method for Kalman Filter-Based Turbofan Engine Thrust Estimation

A new linear point design technique is presented for the determination of tuning parameters that enable the optimal estimation of unmeasured engine outputs such as thrust. The engine s performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters related to each major engine component. Accurate thrust reconstruction depends upon knowledge of these health parameters, but there are usually too few sensors to be able to estimate their values. In this new technique, a set of tuning parameters is determined which accounts for degradation by representing the overall effect of the larger set of health parameters as closely as possible in a least squares sense. The technique takes advantage of the properties of the singular value decomposition of a matrix to generate a tuning parameter vector of low enough dimension that it can be estimated by a Kalman filter. A concise design procedure to generate a tuning vector that specifically takes into account the variables of interest is presented. An example demonstrates the tuning parameters ability to facilitate matching of both measured and unmeasured engine outputs, as well as state variables. Additional properties of the formulation are shown to lend themselves well to diagnostics.

Litt, Jonathan S.↗

An Optimal Orthogonal Decomposition Method for Kalman Filter-Based Turbofan Engine Thrust Estimation

A new linear point design technique is presented for the determination of tuning parameters that enable the optimal estimation of unmeasured engine outputs, such as thrust. The engine's performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters related to each major engine component. Accurate thrust reconstruction depends on knowledge of these health parameters, but there are usually too few sensors to be able to estimate their values. In this new technique, a set of tuning parameters is determined that accounts for degradation by representing the overall effect of the larger set of health parameters as closely as possible in a least squares sense. The technique takes advantage of the properties of the singular value decomposition of a matrix to generate a tuning parameter vector of low enough dimension that it can be estimated by a Kalman filter. A concise design procedure to generate a tuning vector that specifically takes into account the variables of interest is presented. An example demonstrates the tuning parameters ability to facilitate matching of both measured and unmeasured engine outputs, as well as state variables. Additional properties of the formulation are shown to lend themselves well to diagnostics.

Litt, Jonathan S.↗

Lunar occultation of Saturn. I - The diameters of Tethys, Dione, Rhea, Titan, and Iapetus

The diameters of Tethys, Dione, Rhea, Titan and Iapetus were determined from observations of their March 30, 1974, lunar occultations, made with the Mauna Kea 224 and 61 cm telescopes. Light curves were obtained simultaneously in four colors, and the difference between the time of occultation at the two telescopes provided a direct measurement of the slope of the lunar limb, found to be small in all cases. The satellite diameters were determined by least-squares fits of model occultation light curves to the data. Titan's diameter must be at least 5800 km, much larger than the currently accepted value of 5000 km, making it the largest satellite in the solar system. This larger diameter implies a low mean density. For the other four satellites arguments are presented in favor of accepting the occultation diameters corresponding to limb darkened disks. Except for Titan, the lunar occultation diameters generally agree with previous diskmeter and radiometric determinations.

Elliot, J. L.↗

A decentralized square root information filter/smoother

A number of developments has recently led to a considerable interest in the decentralization of linear least squares estimators. The developments are partly related to the impending emergence of VLSI technology, the realization of parallel processing, and the need for algorithmic ways to speed the solution of dynamically decoupled, high dimensional estimation problems. A new method is presented for combining Square Root Information Filters (SRIF) estimates obtained from independent data sets. The new method involves an orthogonal transformation, and an information matrix filter 'homework' problem discussed by Schweppe (1973) is generalized. The employed SRIF orthogonal transformation methodology has been described by Bierman (1977).

Bierman, G. J.↗

Numerical algorithms for steady and unsteady incompressible Navier-Stokes equations

The numerical analysis of the incompressible Navier-Stokes equations are becoming important tools in the understanding of some fluid flow problems which are encountered in research as well as in industry. With the advent of the supercomputers, more realistic problems can be studied with a wider choice of numerical algorithms. An alternative formulation is presented for viscous incompressible flows. The incompressible Navier-Stokes equations are cast in a velocity/vorticity formulation. This formulation consists of solving the Poisson equations for the velocity components and the vorticity transport equation. Two numerical algorithms for the steady two-dimensional laminar flows are presented. The first method is based on the actual partial differential equations. This uses a finite-difference approximation of the governing equations on a staggered grid. The second method uses a finite element discretization with the vorticity transport equation approximated using a Galerkin approximation and the Poisson equations are obtained using a least squares method. The equations are solved efficiently using Newton's method and a banded direct matrix solver (LINPACK). The method is extended to steady three-dimensional laminar flows and applied to a cubic driven cavity using finite difference schemes and a staggered grid arrangement on a Cartesian mesh. The equations are solved iteratively using a plane zebra relaxation scheme. Currently, a two-dimensional, unsteady algorithm is being developed using a generalized coordinate system. The equations are discretized using a finite-volume approach. This work will then be extended to three-dimensional flows.

Hafez, Mohammed↗

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗

An Optimal Orthogonal Decomposition Method for Kalman Filter-Based Turbofan Engine Thrust Estimation

A new linear point design technique is presented for the determination of tuning parameters that enable the optimal estimation of unmeasured engine outputs, such as thrust. The engine s performance is affected by its level of degradation, generally described in terms of unmeasurable health parameters related to each major engine component. Accurate thrust reconstruction depends on knowledge of these health parameters, but there are usually too few sensors to be able to estimate their values. In this new technique, a set of tuning parameters is determined that accounts for degradation by representing the overall effect of the larger set of health parameters as closely as possible in a least-squares sense. The technique takes advantage of the properties of the singular value decomposition of a matrix to generate a tuning parameter vector of low enough dimension that it can be estimated by a Kalman filter. A concise design procedure to generate a tuning vector that specifically takes into account the variables of interest is presented. An example demonstrates the tuning parameters ability to facilitate matching of both measured and unmeasured engine outputs, as well as state variables. Additional properties of the formulation are shown to lend themselves well to diagnostics.

Litt, Jonathan S.↗

The single scattering phase functions of Jupiter's clouds

The determination of the single scattering phase functions of Jupiter's clouds and a thin upper haze by Tomasko et al. was refined and extended to seven latitudes in blue and red light. The phase function is well-constrained by the Pioneer 10 and 11 photometric data sets. Multiple scattering models were computed to match the limb darkening at each latitude at up to 15 phase angles from 12 deg to 151 deg. Ground-based observations were used for absolute calibration and to extend the data to lower phase angles. The phase functions were parameterized using the double Henyey-Greenstein function. The three Henyey-Greenstein parameters and the single scattering albedo were determined using a non-linear least squares method for the haze and the clouds below. The phase functions derived for the northen zone and belt are remarkably similar to the phase functions of the corresponding regions in the south, with most of the differences in brightness of the northern and southern features resulting from minor differences in single scattering albedo. Analysis of the Equatorial Region is complicated by the presence of numerous small features, but the phase function required is generally similar to that seen in the more homogeneous regions. Details of the phase functions of the haze and clouds are presented, and the differences between the cloud phase functions at low and high latitudes in red and blue light are discussed.

Doose, L. R.↗

Parameter identification in distributed systems

This paper describes a method for the identification of the parameters entering into the equations of motion of distributed systems. Because the motion of distributed systems is described in terms of partial differential equations, these parameters are in general continuous functions of the spatial variables. For vibrating systems, these parameters ordinarily represent the mass, stiffness and damping distributions. In this paper, these distributions are expanded in terms of finite series of known functions of the spatial variables multiplied by undetermined coefficients. It is assumed that the nature of the equations of motion is known and that a limited number of eigenvalues and eigenfunctions is identified in advance. Use is then made of the least squares method, in conjunction with the eigenfunctions' orthogonality, to compute the undetermined coefficients, thus identifying the system distributed parameters. A method for the identification of the eigensolution is also presented. The procedure for the identification of the eigensolution and of the system parameters is demonstrated via a numerical example.

Baruh, H.↗

Variability simulations with a steady, linearized primitive equations model

Solutions of the steady, primitive equations on a sphere, linearized about a zonally symmetric basic state are computed for the purpose of simulating monthly mean variability in the troposphere. The basic states are observed, winter monthly mean, zonal means of zontal and meridional velocities, temperatures and surface pressures computed from the 15 year NMC time series. A least squares fit to a series of Legendre polynomials is used to compute the basic states between 20 H and the equator, and the hemispheres are assumed symmetric. The model is spectral in the zonal direction, and centered differences are employed in the meridional and vertical directions. Since the model is steady and linear, the solution is obtained by inversion of a block, pente-diagonal matrix. The model simulates the climatology of the GFDL nine level, spectral general circulation model quite closely, particularly in middle latitudes above the boundary layer. This experiment is an extension of that simulation to examine variability of the steady, linear solution.

Kinter, J. L., III↗

Acoustooptic linear algebra processors - Architectures, algorithms, and applications

Architectures, algorithms, and applications for systolic processors are described with attention to the realization of parallel algorithms on various optical systolic array processors. Systolic processors for matrices with special structure and matrices of general structure, and the realization of matrix-vector, matrix-matrix, and triple-matrix products and such architectures are described. Parallel algorithms for direct and indirect solutions to systems of linear algebraic equations and their implementation on optical systolic processors are detailed with attention to the pipelining and flow of data and operations. Parallel algorithms and their optical realization for LU and QR matrix decomposition are specifically detailed. These represent the fundamental operations necessary in the implementation of least squares, eigenvalue, and SVD solutions. Specific applications (e.g., the solution of partial differential equations, adaptive noise cancellation, and optimal control) are described to typify the use of matrix processors in modern advanced signal processing.

Casasent, D.↗

Estimating errors in least-squares fitting

While least-squares fitting procedures are commonly used in data analysis and are extensively discussed in the literature devoted to this subject, the proper assessment of errors resulting from such fits has received relatively little attention. The present work considers statistical errors in the fitted parameters, as well as in the values of the fitted function itself, resulting from random errors in the data. Expressions are derived for the standard error of the fit, as a function of the independent variable, for the general nonlinear and linear fitting problems. Additionally, closed-form expressions are derived for some examples commonly encountered in the scientific and engineering fields, namely ordinary polynomial and Gaussian fitting functions. These results have direct application to the assessment of the antenna gain and system temperature characteristics, in addition to a broad range of problems in data analysis. The effects of the nature of the data and the choice of fitting function on the ability to accurately model the system under study are discussed, and some general rules are deduced to assist workers intent on maximizing the amount of information obtained form a given set of measurements.

Richter, P. H.↗

Efficient identification of critical stresses in structures subject to dynamic loads

Optimum structural design problems generally employ constraints which are parametric in terms of space and time variables. A parametric constraint may be replaced by equivalent critical point constraints at its local minima for optimization applications. In complex structures, accurate identification of such critical points is computationally expensive due to the cost of finite element analyses. Three techniques are described for efficiently and accurately identifying critical points for space- and time-dependent parametric constraints. An adaptive search technique and a spline interpolation technique are developed for exactly known response. A least squares spline approximation is suggested for noisy behavior. A helicopter tail-boom structure subjected to transient loading is used as an example to demonstrate the techniques described. All three techniques are shown to be computationally efficient for critical point identification and the least squares approximation also removes noise from the data. The case of multiple constraints per element is shown to be particularly suited to the use of spline techniques.

Haftka, R. T.↗

Shade images of forested areas obtained from LANDSAT MSS data

The pixel size in the present day Remote Sensing systems is large enough to include different types of land cover. Depending upon the target area, several components may be present within the pixel. In forested areas, generally, three main components are present: tree canopy, soil (understory), and shadow. The objective is to generate a shade (shadow) image of forested areas from multispectral measurements of LANDSAT MSS (Multispectral Scanner) data by implementing a linear mixing model, where shadow is considered as one of the primary components in a pixel. The shade images are related to the observed variation in forest structure, i.e., the proportion of inferred shadow in a pixel is related to different forest ages, forest types, and tree crown cover. The Constrained Least Squares (CLS) method is used to generate shade images for forest of eucalyptus and vegetation of cerrado using LANDSAT MSS imagery over Itapeva study area in Brazil. The resulted shade images may explain the difference on ages for forest of eucalyptus and the difference on three crown cover for vegetation of cerrado.

Shimabukuro, Yosio Edemir↗

Uncertainty Analysis of Slug Calorimeters in the HyMETS Arc-Jet Facility

The objective of this work is to perform an uncertainty analysis of the deduced stagnation heat flux environment on a slug calorimeter for conditions that span the performance envelope of the Hypersonic Materials Environmental Test System arc-jet facility located at NASA Langley Research Center. Analytical solutions are developed for boundary-value problems on the slug element accounting for non-ideal effects, including spatial variation in the slug heat flux, multi-dimensional thermal conduction, and back-face losses, which departs from the state-of-the-art method derived from the American Society of Testing and Materials. Boundary-value problem definitions are informed by preliminary finite element thermal analysis of the slug calorimeter assembly (including both slug and housing) and just the slug element. The analytical solutions are presented in a general sense and in a truncated form from error analysis. Results are shown in optimizing and validating the analytical models against available slug back-face thermal data. The optimization results indicate that the appropriate epistemic uncertainty of the deduced stagnation heat flux on the slug calorimeter is at most±2.5% for both a high-and low-enthalpy test condition. In addition, a numerical approach is used to determine the aleatory (probabilistic) uncertainty component in the slug stagnation heat flux by applying a marching least-squares slope routine through the steady-state portion of the slug back-face thermal response. Results indicate a compromise between the number of samples and the filter frequency of slug back-face thermal data points when evaluating the standard deviation of the deduced stagnation heat flux statistics. When combining the mixed uncertainty, both aleatory and epistemic, the interval of uncertainty in the deduced stagnation heat flux is determined to be up to ±4%, which is at least a 60% reduction from the standard uncertainty used in the state-of-the-art method.

uncertainty↗