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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quasielastic least-squares reverse time migration of PS reflections

Conventional analysis of amplitude variation with offset for elastic PS reflections is based on analytical reflection coefficients in a layered medium, and wave-equation-based PS migration is mainly used to produce a structural image. To overcome this problem, we have developed a least-squares reverse time migration (LSRTM) method for elastic PS reflections based on a quasielastic wave equation. Here, the quasielastic wave equation can accurately model PS reflections with elastic amplitudes under the first-order Born approximation. Our LSRTM method inverts for perturbations of the S-wave velocity and density by minimizing the L 2 norm of the difference between recorded and predicted PS reflections modeled using a quasielastic wave equation. We refer to our new method as quasielastic LSRTM of PS reflections. Numerical tests on synthetic and field data indicate that our method can properly handle the amplitudes of elastic PS reflections and provides an accurate estimate of the perturbations of S-wave velocity and density. Extending the method to the 3D case is not straightforward and might require incorporating certain PS data processing techniques into the inversion itself.

58 GEOSCIENCES↗
NASA NTRS2008

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)↗
NASA NTRS1992

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.↗
NASA NTRS1984

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.↗
NASA NTRS1984

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.↗
NASA NTRS1992

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.↗
NASA NTRS2010

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.↗
NASA NTRS2013

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.↗
NASA NTRS2014

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↗
NASA NTRS1993

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.↗
NASA NTRS2018

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↗
DOE OSTI2023

A leaf-level spectral library to support high-throughput plant phenotyping: predictive accuracy and model transfer

Abstract Leaf-level hyperspectral reflectance has become an effective tool for high-throughput phenotyping of plant leaf traits due to its rapid, low-cost, multi-sensing, and non-destructive nature. However, collecting samples for model calibration can still be expensive, and models show poor transferability among different datasets. This study had three specific objectives: first, to assemble a large library of leaf hyperspectral data (n=2460) from maize and sorghum; second, to evaluate two machine-learning approaches to estimate nine leaf properties (chlorophyll, thickness, water content, nitrogen, phosphorus, potassium, calcium, magnesium, and sulfur); and third, to investigate the usefulness of this spectral library for predicting external datasets (n=445) including soybean and camelina using extra-weighted spiking. Internal cross-validation showed satisfactory performance of the spectral library to estimate all nine traits (mean R2=0.688), with partial least-squares regression outperforming deep neural network models. Models calibrated solely using the spectral library showed degraded performance on external datasets (mean R2=0.159 for camelina, 0.337 for soybean). Models improved significantly when a small portion of external samples (n=20) was added to the library via extra-weighted spiking (mean R2=0.574 for camelina, 0.536 for soybean). The leaf-level spectral library greatly benefits plant physiological and biochemical phenotyping, whilst extra-weight spiking improves model transferability and extends its utility.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

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↗
NASA NTRS2000

Initial OSSE Sky Survey Observation: Deep Survey in the Virgo Region

The OSSE Virgo survey covered about 1% of the sky, using eight weeks of observing time on the Compton Gamma Ray Observatory, and was performed in conjunction with the EGRET and COMPTEL instruments on CGRO. A fundamental aspect of the survey was the development of data reduction techniques for the scanning data from the OSSE detectors. A direct algebraic approach for spatially deconvolving the data was developed. For the survey work, no background observations isolated from the survey region were acquired, unlike normal discrete source observations. The survey, reduction method covers the mapped region with a grid of pixels, and estimates the flux from each pixel by least-squares agreement with the data, assuming positivity constraints. This technique has led to the detection of compact sources at hard X-ray energies within the Virgo region, and the OSSE survey studies have been extended to the south galactic pole and north Ecliptic pole, with similar analysis techniques.

Cameron, Robert↗
DOE OSTI2025

Deciphering the Scattering of Mechanically Driven Polymers Using Deep Learning

Here, we present a deep learning approach for analyzing two-dimensional scattering data of semiflexible polymers under external forces. In our framework, scattering functions are compressed into a three-dimensional latent space using a Variational Autoencoder (VAE), and two converter networks establish a bidirectional mapping between the polymer parameters (bending modulus, stretching force, and steady shear) and the scattering functions. The training data are generated using off-lattice Monte Carlo simulations to avoid the orientational bias inherent in lattice models, ensuring robust sampling of polymer conformations. The feasibility of this bidirectional mapping is demonstrated by the organized distribution of polymer parameters in the latent space. By integrating the converter networks with the VAE, we obtain a generator that produces scattering functions from given polymer parameters and an inferrer that directly extracts polymer parameters from scattering data. While the generator can be utilized in a traditional least-squares fitting procedure, the inferrer produces comparable results in a single pass and operates 3 orders of magnitude faster. This approach offers a scalable automated tool for polymer scattering analysis and provides a promising foundation for extending the method to other scattering models, experimental validation, and the study of time-dependent scattering data.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
NASA NTRS1989

Sampling characteristics of satellite orbits

The irregular space-time sampling of any finite region by an orbiting satellite raises difficult questions as to which frequencies and wavenumbers can be determined and which will alias into others. Conventional sampling theorems must be extended to account for both irregular data distributions and observational noise - the sampling irregularity making the system much more susceptible to noise than in regularly sampled cases. The problem is formulated here in terms of least-squares and applied to spacecraft in 10-day and 17-day repeating orbits. The 'diamond-pattern' laid down spatially in such repeating orbits means that either repeat period adequately samples the spatial variables, but the slow overall temporal coverage in the 17-day pattern leads to much greater uncertainty than in the shorter repeat cycle. The result is not definitive and it is not concluded that a 10-day orbit repeat is the most appropriate one. A major conclusion, however, is that different orbital choices have potentially quite different sampling characteristics which need to be analyzed in terms of the spectral characteristics of the moving sea surface.

Wunsch, Carl↗
NASA NTRS1990

Optimal model reduction and frequency-weighted extension

In this paper the quadratically optimal model reduction problem for single-input, single-output systems is considered. The reduced order model is determined by minimizing the integral of the magnitude-squared of the transfer function error. It is shown that the numerator coefficients of the optimal approximant satisfy a weighted least squares problem and, on this basis, a two-step iterative algorithm is developed combining a least squares solver with a gradient minimizer. The existence of globally optimal stable solutions to the optimization problem is established, and convergence of the algorithm to stationary values of the cost function is proved. The formulation is extended to handle the frequency-weighted optimal model reduction problem. Three examples demonstrate the optimization algorithm.

Spanos, J. T.↗
NASA NTRS1986

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.↗