Search NASA⌕ Search

SEARCH · Search NASA

Results for “extended least squares”

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 73 records · Page 4

Toward autonomous spacecraft

Ways in which autonomous behavior of spacecraft can be extended to treat situations wherein a closed loop control by a human may not be appropriate or even possible are explored. Predictive models that minimize mean least squared error and arbitrary cost functions are discussed. A methodology for extracting cyclic components for an arbitrary environment with respect to usual and arbitrary criteria is developed. An approach to prediction and control based on evolutionary programming is outlined. A computer program capable of predicting time series is presented. A design of a control system for a robotic dense with partially unknown physical properties is presented.

Fogel, L. J.↗

Time delay occultation data of the Helios spacecrafts and preliminary analysis for probing the solar corona

S-band time delay measurements were collected from the spacecrafts Helios A and B during three solar occultations in 1975/76 within heliocentric distances of about 3 and 215 solar radii in terms of range, Doppler frequency shift, and electron content. A description is given concerning some characteristic features of the methods of measurement and data processing. Typical data sets are discussed to probe the electron density distribution near the sun (west and east limb as well) including the outer and extended corona. Steady-state and dynamical aspects of the solar corona are presented and compared with earth-bound K-coronagraph measurements. Using a weighted least squares estimation 3 parameters of an average coronal electron density profile are derived in a preliminary analysis to yield electron densities of about 130 billion; 100 million; 7 million/cu m at r?3; 65; 215 solar radii. Transient phenomena are discussed and a velocity of propagation v approximately 900 km/s for plasma ejecta from a solar flare is determined from an extraordinary set of Helios B electron content measurements on April 30/May 1, 1976.

Edenhofer, P.↗

Least-squares finite element solution of 3D incompressible Navier-Stokes problems

Although significant progress has been made in the finite element solution of incompressible viscous flow problems. Development of more efficient methods is still needed before large-scale computation of 3D problems becomes feasible. This paper presents such a development. The most popular finite element method for the solution of incompressible Navier-Stokes equations is the classic Galerkin mixed method based on the velocity-pressure formulation. The mixed method requires the use of different elements to interpolate the velocity and the pressure in order to satisfy the Ladyzhenskaya-Babuska-Brezzi (LBB) condition for the existence of the solution. On the other hand, due to the lack of symmetry and positive definiteness of the linear equations arising from the mixed method, iterative methods for the solution of linear systems have been hard to come by. Therefore, direct Gaussian elimination has been considered the only viable method for solving the systems. But, for three-dimensional problems, the computer resources required by a direct method become prohibitively large. In order to overcome these difficulties, a least-squares finite element method (LSFEM) has been developed. This method is based on the first-order velocity-pressure-vorticity formulation. In this paper the LSFEM is extended for the solution of three-dimensional incompressible Navier-Stokes equations written in the following first-order quasi-linear velocity-pressure-vorticity formulation.

Jiang, Bo-Nan↗

Development of an unsteady aerodynamic analysis for finite-deflection subsonic cascades

An unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, was developed to predict aerodynamic forces and moments associated with free vibration or flutter phenomena in the fan, compressor, or turbine stages of modern jet engines. Based on the assumption of small amplitude blade motions, the unsteady flow is governed by linear equations with variable coefficients which depend on the underlying steady low. These equations were approximated using difference expressions determined from an implicit least squares development and applicable on arbitrary grids. The resulting linear system of algebraic equations is block tridiagonal, which permits an efficient, direct (i.e., noniterative) solution. The solution procedure was extended to treat blades with rounded or blunt edges at incidence relative to the inlet flow.

Verdon, J. M.↗

Scattering of 42-MeV alpha particles from Cu-65

The extended particle-core coupling model was used to predict the properties of low-lying levels of Cu-65. A 42-MeV alpha particle cyclotron beam was used for the experiment. The experiment included magnetic analysis of the incident beam and particle detection by lithium-drifted silicon semiconductors. Angular distributions were measured for 10 to 50 degrees in the center of mass system. Data was reduced by fitting the peaks with a skewed Gaussian function using a least squares computer program with a linear background search. The energy calibration of each system was done by pulsar, and the excitation energies are accurate to + or - 25 keV. The simple weak coupling model cannot account for the experimentally observed quantities of the low-lying levels of Cu-65. The extended particle-core calculation showed that the coupling is not weak and that considerable configuration mixing of the low-lying states results.

Stewart, W. M.↗

The Least-Squares Calibration on the Micro-Arcsecond Metrology Test Bed

The Space Interferometry Mission (S1M) will measure optical path differences (OPDs) with an accuracy of tens of picometers, requiring precise calibration of the instrument. In this article, we present a calibration approach based on fitting star light interference fringes in the interferometer using a least-squares algorithm. The algorithm is first analyzed for the case of a monochromatic light source with a monochromatic fringe model. Using fringe data measured on the Micro-Arcsecond Metrology (MAM) testbed with a laser source, the error in the determination of the wavelength is shown to be less than 10pm. By using a quasi-monochromatic fringe model, the algorithm can be extended to the case of a white light source with a narrow detection bandwidth. In SIM, because of the finite bandwidth of each CCD pixel, the effect of the fringe envelope can not be neglected, especially for the larger optical path difference range favored for the wavelength calibration.

interferometer calibrations↗

Multiwavelength pyrometry for nongray surfaces in the presence of interfering radiation

A NASA developed multiwavelength pyrometry technique for nongray surfaces was extended to also measure surface temperature in the presence of interfering radiation. This radiation is produced by heat lamps used to raise the temperature of the surface. The necessary instruments are a spectral radiometer, an auxiliary radiation source, and a computer. Four radiation spectra are recorded: (1) the unobstructed spectrum characterizing an auxiliary radiation source; (2) the unobstructed spectrum characterizing the interfering radiation; (3) the radiation spectrum consisting of surface emission plus the interfering radiation; and (4) a spectrum consisting of the radiations of (3) plus the reflected radiation due to the incidence of the auxiliary radiation source on this surface. With these spectra, application of two variable, nonlinear, least squares, curve fitting computer software determines the surface temperature and the spectral emissivity. Use of the method to measure the surface temperature of silicon carbide under a simulated interference condition is shown at a low temperature just above ambient. The instrumentation necessary to extend the method to elevated temperatures is discussed.

Ng, Daniel↗

Sparse Regression as a Sparse Eigenvalue Problem

We extend the l0-norm "subspectral" algorithms for sparse-LDA [5] and sparse-PCA [6] to general quadratic costs such as MSE in linear (kernel) regression. The resulting "Sparse Least Squares" (SLS) problem is also NP-hard, by way of its equivalence to a rank-1 sparse eigenvalue problem (e.g., binary sparse-LDA [7]). Specifically, for a general quadratic cost we use a highly-efficient technique for direct eigenvalue computation using partitioned matrix inverses which leads to dramatic x103 speed-ups over standard eigenvalue decomposition. This increased efficiency mitigates the O(n4) scaling behaviour that up to now has limited the previous algorithms' utility for high-dimensional learning problems. Moreover, the new computation prioritizes the role of the less-myopic backward elimination stage which becomes more efficient than forward selection. Similarly, branch-and-bound search for Exact Sparse Least Squares (ESLS) also benefits from partitioned matrix inverse techniques. Our Greedy Sparse Least Squares (GSLS) generalizes Natarajan's algorithm [9] also known as Order-Recursive Matching Pursuit (ORMP). Specifically, the forward half of GSLS is exactly equivalent to ORMP but more efficient. By including the backward pass, which only doubles the computation, we can achieve lower MSE than ORMP. Experimental comparisons to the state-of-the-art LARS algorithm [3] show forward-GSLS is faster, more accurate and more flexible in terms of choice of regularization

Exact Sparse Least Squares (ESLS)↗

Coded aperture devices for viewing extended objects from space

Coded aperture sensors for photons or energetic neutral atoms (ENAs), which incorporate FOV limiters and subdivide the object field into a number of elements which is smaller than the number of detector pixels, are described. A least squares fit to the data is made in reconstructing the object field. To evaluate the optics and reconstruction algorithms, two 'breadboard' sensors have been constructed, one based on a film camera and the other upon a UV-light sensitive microchannel plate detector system. Results obtained thus far show that the concept is viable, and no special difficulties should be encountered in adapting the detector geometries to neutral particle imaging systems. Charged particle deflection plates could be incorporated into the region between the FOV limiter and the aperture, or installed ahead of the limiter.

Curtis, C. C.↗

Light Scattered from Polished Optical Surfaces: Wings of the Point Spread Function

Random figure errors from the polishing process plus particles on the main mirrors in a telescope cause an extended point spread function (PSF) declining approximately as the inverse square of the sine of the angle from a star from about 100 micro-rad to a right angle. The decline in at least one case, and probably in general, proceeds as the inverse cube at smaller angles where the usual focal plane aperture radius is chosen. The photometric error due to misalignment by one Airy ring spacing with an aperture of n rings depends on the net variance in the figure. It is approximately 60/(n+1)(3) when using the data of Kormendy (1973). A typical value is 6 x 10 to the -5th power per ring of misalignment with n = 100 rings. The encircled power may be modulated on a time scale of hours by parts per thousand in a wavelength dependent manner due to relative humidity effects on mirror dust. The scattering according to an inverse power law is due to a random walk in aberration height caused by a multitude of facets and slope errors left by the polishing process. A deviation from such a law at grazing emergence may permit monitoring the dust effects.

Kenknight, C. E.↗

A linear aerodynamic analysis for unsteady transonic cascades

A potential flow analysis to predict unsteady airloads produced by the vibrations of turbomachinery blades operating at transonic Mach numbers is presented. The unsteady aerodynamic model includes the effects of blade geometry, finite mean pressure variation across the blade row, high frequency blade motion, and shock motion within the framework of a linearized, frequency domain formulation. The unsteady equations are solved implicit, least squares, finite difference approximation which is applicable on arbitrary grids. A numerical solution for the entire unsteady field is determined by matching a solution determined on a rectilinear type cascade mesh, which covers an extended blade passage region, to a solution determined on a detailed polar type local mesh, which covers and extends well beyond the supersonic region(s) adjacent to a blade surface. Cascades of double circular arc and flat plate blades demonstrate the unsteady analysis, and partially illustrate the effects of blade geometry, inlet Mach number, blade vibration frequency and shock motion on unsteady response.

Verdon, J. M.↗

Boundary layer relaminarization device

Relamination of a boundary layer formed in supersonic flow over the leading edge of a swept airfoil is accomplished by means of at least one band, especially a quadrangular band, and most preferably a square band. Each band conforms to the leading edge and the upper and lower surfaces of the airfoil as an integral part thereof and extends perpendicularly from the leading edge. Each band has a height of about two times the thickness of the maximum expected boundary layer.

Creel, Theodore R.↗

Evaluating Descent and Ascent Trajectories Near Non-Spherical Bodies

Spacecraft landing on small bodies pass through regions where conventional gravitation formulations using exterior spherical harmonics are inaccurate. An investigation shows that a formulation using interior solid spherical harmonics might be satisfactory. Interior spherical harmonic expansions are usable inside an imaginary, empty sphere. For this application, such a sphere could be positioned in empty space above the intended landing site and rotating with the body. When the spacecraft is inside this sphere, the interior harmonic expansion would be used instead of the conventional, exterior harmonic expansion. Coefficients can be determined by a least-squares fit to gravitation measurements synthesized from conventional formulations. Due to their unfamiliarity, recurrences for interior, as well as exterior, expansions are derived. Hotine's technique for partial derivatives of exterior spherical harmonics is extended to interior harmonics.

Werner, Robert A.↗

Further Development of Rotating Rake Mode Measurement Data Analysis

The Rotating Rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. After analysis of the measured data, the mode amplitudes and phases were quantified. For low-speed fans within axisymmetric ducts, mode power levels computed from rotating rake measured data would agree with the far-field power levels on a tone by tone basis. However, this agreement required that the sound from the noise sources within the duct propagated outward from the duct exit without reflection at the exit and previous studies suggested conditions could exist where significant reflections could occur. To directly measure the modes propagating in both directions within a duct, a second rake was mounted to the rotating system with an offset in both the axial and the azimuthal directions. The rotating rake data analysis technique was extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode levels at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode amplitudes for the modes propagating in both directions within the duct. The fit equations were also modified to allow evanescent mode amplitudes to be computed. This extension of the rotating rake data analysis technique was tested using simulated data, numerical code produced data, and preliminary in-duct measured data.

Dahl, Milo D.↗

Analysis of Dual Rotating Rake Data from the NASA Glenn Advanced Noise Control Fan Duct with Artificial Sources

The Rotating Rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. Initially, the mode amplitudes and phases were quantified from a single rake measurement at one axial location. To directly measure the modes propagating in both directions within a duct, a second rake was mounted to the rotating system with an offset in both the axial and the azimuthal directions. The rotating rake data analysis technique was then extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode levels at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode amplitudes for the modes propagating in both directions within the duct. Validation experiments have been conducted using artificial acoustic sources. Results are shown for the measurement of the standing waves in the duct from sound generated by one and two acoustic sources that are separated into the component modes propagating in both directions within the duct. Measured reflection coefficients from the open end of the duct are compared to analytical predictions.

Noise↗

Boundary Layer Relaminarization Device

Relamination of a boundary layer formed in supersonic flow over the leading edge of a swept airfoil is accomplished using at least one band, especially a quadrangular band, and most preferably a square band. Each band conforms to the leading edge and the upper and lower surfaces of the airfoil as an integral part thereof and extends perpendicularly from the leading edge. Each band has a height of about two times the thickness of the maximum expected boundary layer.

Creel, Theodore R.↗

Dual Rotating Rake Measurements of Higher-Order Duct Modes: Validation Using Experimental and Numerical Data

A rotating rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. After analysis of the measured data, the mode coefficient amplitudes and phases were quantified. Early studies using this system found that mode power levels computed from rotating rake measured data would agree with the far-field power levels. However, this agreement required that the sound from the noise sources within the duct propagated outward from the duct exit without reflection and previous studies suggested conditions could exist where significant reflections could occur. This paper shows that mounting a second rake to the rotating system, with an offset in both the axial and the azimuthal directions, measures the data necessary to determine the modes propagating in both directions within a duct. The rotating rake data analysis technique was extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode coefficients at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode coefficients for the modes propagating in both directions within the duct while accounting for the presence of evanescent modes. The validation of the dual-rotating-rake measurements was conducted using data from a combination of experiments and numerical calculations to compute reflection coefficients and other mode coefficient ratios. Compared to results from analytical and numerical computations, the results from dual-rotating-rake measured data followed the expected trends when frequency, mode number, and duct termination geometry were changed.

Validation↗

Multiwavelength pyrometry for nongray bodies

A multiwavelength technique was developed and applied to measure the temperatures of nongray surfaces. The instruments required are a spectral radiometer, a dedicated auxiliary radiation source, and a computer. In general, three radiation spectra are recorded: (1) spectrum S sub 0 of the auxiliary radiation source; (2) spectrum S sub 1 of the surface-emitted radiation; and (3) spectrum S sub 2, the sum of the radiation of S sub 1 plus the reflected radiation due to the incidence of the auxiliary radiation source on the surface. Subtracting spectrum S sub 1 from spectrum Sub 2 yields the reflection spectrum resulting from the incidence radiation. From these spectra, a quantity z(lambda) is derived and is related to the reflectivity r(lambda) by r(lambda) = z(lambda)/f, where f is a constant. Spectrum S sub 1 is represented mathematically as the product of a wavelength-dependent emissivity obtained from Kirchhoff's law and a Planck function of temperature T. Application of two-variable (lambda and z), nonlinear, least-squares curve-fitting computer software to fit spectrum S sub 1 to this mathematical expression yielded the surface temperature. This technique also measured the spectral reflectivity and emissivity of the surface. Instrumentation necessary to extend measurement to elevated temperatures and in the presence of reflective interference is discussed.

Ng, Daniel↗