SEARCH · Engineering Papers

Results for “least square approximations”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Comparative Studies for the Sodium and Potassium Atmospheres of the Moon and Mercury

A summary discussion of recent sodium and potassium observations for the atmospheres of the Moon and Mercury is presented with primary emphasis on new full-disk images that have become available for sodium. For the sodium atmosphere, image observations for both the Moon and Mercury are fitted with model calculations (1) that have the same source speed distribution, one recently measured for electron-stimulated desorption and thought to apply equally well to photon-stimulated desorption, (2) that have similar average surface sodium fluxes, about 2.8 x 10(exp 5) to 8.9 x 10(exp 5) atoms cm(exp -2)s(exp -1) for the Moon and approximately 3.5 x 10(exp 5) to 1.4 x 10(exp 6) atoms cm(exp -2)s(exp -1) for Mercury, but (3) that have very different distributions for the source surface area. For the Moon, a sunlit hemispherical surface source of between approximately 5.3 x 10(exp 22) to 1.2 x 10(exp 23) atoms/s is required with a spatial dependence at least as sharp as the square of the cosine of the solar zenith angle. For Mercury, a time dependent source that varies from 1.5 x 10(exp 22) to 5.8 x l0(exp 22) atoms/s is required which is confined to a small surface area located at, but asymmetrically distributed about, the subsolar point. The nature of the Mercury source suggest that the planetary magnetopause near the subsolar point acts as a time varying and partially protective shield through which charged particles may pass to interact with and liberate gas from the planetary surface. Suggested directions for future research activities are discussed.

Smyth, William H.↗

NASA NTRS1999

Gaussian Decomposition of Laser Altimeter Waveforms

We develop a method to decompose a laser altimeter return waveform into its Gaussian components assuming that the position of each Gaussian within the waveform can be used to calculate the mean elevation of a specific reflecting surface within the laser footprint. We estimate the number of Gaussian components from the number of inflection points of a smoothed copy of the laser waveform, and obtain initial estimates of the Gaussian half-widths and positions from the positions of its consecutive inflection points. Initial amplitude estimates are obtained using a non-negative least-squares method. To reduce the likelihood of fitting the background noise within the waveform and to minimize the number of Gaussians needed in the approximation, we rank the "importance" of each Gaussian in the decomposition using its initial half-width and amplitude estimates. The initial parameter estimates of all Gaussians ranked "important" are optimized using the Levenburg-Marquardt method. If the sum of the Gaussians does not approximate the return waveform to a prescribed accuracy, then additional Gaussians are included in the optimization procedure. The Gaussian decomposition method is demonstrated on data collected by the airborne Laser Vegetation Imaging Sensor (LVIS) in October 1997 over the Sequoia National Forest, California.

Hofton, Michelle A.↗

NASA NTRS2016

Unsteady Aerodynamic Force Sensing from Measured Strain

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm. A cantilevered rectangular wing built and tested at the NASA Langley Research Center (Hampton, Virginia, USA) in 1959 is used to validate the simple approach. Unsteady aerodynamic forces as well as wing deflections, velocities, accelerations, and strains are computed using the CFL3D computational fluid dynamics (CFD) code and an MSC/NASTRAN code (MSC Software Corporation, Newport Beach, California, USA), and these CFL3D-based results are assumed as measured quantities. Based on the measured strains, wing deflections, velocities, accelerations, and aerodynamic forces are computed using the proposed approach. These computed deflections, velocities, accelerations, and unsteady aerodynamic forces are compared with the CFL3D/NASTRAN-based results. In general, computed aerodynamic forces based on the lifting surface theory in subsonic speeds are in good agreement with the target aerodynamic forces generated using CFL3D code with the Euler equation. Excellent aeroelastic responses are obtained even with unsteady strain data under the signal to noise ratio of -9.8dB. The deflections, velocities, and accelerations at each sensor location are independent of structural and aerodynamic models. Therefore, the distributed strain data together with the current proposed approaches can be used as distributed deflection, velocity, and acceleration sensors. This research demonstrates the feasibility of obtaining induced drag and lift forces through the use of distributed sensor technology with measured strain data. An active induced drag control system thus can be designed using the two computed aerodynamic forces, induced drag and lift, to improve the fuel efficiency of an aircraft. Interpolation elements between structural finite element grids and the CFD grids and centroids are successfully incorporated with the unsteady aeroelastic computation scheme. The most critical technology for the success of the proposed approach is the robust on-line parameter estimator, since the least-squares curve fitting method depends heavily on aeroelastic system frequencies and damping factors.

shape sensing↗

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 NTRS1979

M2, S2, K1 models of the global ocean tide

Ocean tidal signals appear in many geophysical measurements. Geophysicists need realistic tidal models to aid in interpretation of their data. Because of the closeness to resonance of dissipationless ocean tides, it is difficult for numerical models to correctly represent the actual open ocean tide. As an approximate solution to this problem, test functions derived by solving Laplace's Tidal Equations with ocean loading and self gravitation are used as a basis for least squares dynamic interpolation of coastal and island tidal data for the constituents M2, S2, and Kl. The resulting representations of the global tide are stable over at least a ?5% variation in the mean depth of the model basin, and they conserve mass. Maps of the geocentric tide, the induced free space potential, the induced vertical component of the solid earth tide, and the induced vertical component of the gravitational field for each contituent are presented.

Parke, M. E.↗

NASA NTRS1980

Retrieval of aerosol size distributions from scattering and extinction measurements in the presence of multiple scattering

A fast technique, based on approximation to atmospheric radiative transfer in the solar aureole, is presented for the retrieval of aerosol size distribution from multispectral measurements of both scattered sky radiance and extinction by aerosol media. In this approximation, it is assumed that for a relatively clear day, the sky radiance is due to single scattering by molecules and aerosols, and to multiple scattering by molecules alone. The accuracy of this approximation is discussed in connection with the inversion of scattered radiation data. Results for aerosol size distributions retrieved from multispectral solar extinction measurements by the nonlinear least squares method are discussed, and the retrievals are found to be accurate when checked against ground truth measurements.

Deepak, A.↗

NASA NTRS1993

Space Shuttle Main Engine performance analysis

For a number of years, NASA has relied primarily upon periodically updated versions of Rocketdyne's power balance model (PBM) to provide space shuttle main engine (SSME) steady-state performance prediction. A recent computational study indicated that PBM predictions do not satisfy fundamental energy conservation principles. More recently, SSME test results provided by the Technology Test Bed (TTB) program have indicated significant discrepancies between PBM flow and temperature predictions and TTB observations. Results of these investigations have diminished confidence in the predictions provided by PBM, and motivated the development of new computational tools for supporting SSME performance analysis. A multivariate least squares regression algorithm was developed and implemented during this effort in order to efficiently characterize TTB data. This procedure, called the 'gains model,' was used to approximate the variation of SSME performance parameters such as flow rate, pressure, temperature, speed, and assorted hardware characteristics in terms of six assumed independent influences. These six influences were engine power level, mixture ratio, fuel inlet pressure and temperature, and oxidizer inlet pressure and temperature. A BFGS optimization algorithm provided the base procedure for determining regression coefficients for both linear and full quadratic approximations of parameter variation. Statistical information relative to data deviation from regression derived relations was also computed. A new strategy for integrating test data with theoretical performance prediction was also investigated. The current integration procedure employed by PBM treats test data as pristine and adjusts hardware characteristics in a heuristic manner to achieve engine balance. Within PBM, this integration procedure is called 'data reduction.' By contrast, the new data integration procedure, termed 'reconciliation,' uses mathematical optimization techniques, and requires both measurement and balance uncertainty estimates. The reconciler attempts to select operational parameters that minimize the difference between theoretical prediction and observation. Selected values are further constrained to fall within measurement uncertainty limits and to satisfy fundamental physical relations (mass conservation, energy conservation, pressure drop relations, etc.) within uncertainty estimates for all SSME subsystems. The parameter selection problem described above is a traditional nonlinear programming problem. The reconciler employs a mixed penalty method to determine optimum values of SSME operating parameters associated with this problem formulation.

Santi, L. Michael↗

DOE OSTI2022

maestro

MÆSTRO stands for Multi-fidelity Adaptive Ensemble Stochastic Trust Region Optimization and it is a plug n play derivate fee stochastic optimization solver. The problem being considered in MÆSTRO involves fitting Monte Carlo simulations that describe complex phenomena to experiments. This is done by finding parameters of the resource intensive and noisy simulation that yield the least squares objective function value to the noisy experimental data. This problem is solved using a stochastic trust-region optimization algorithm where in each iteration, a local approximation of the simulation signal and of the simulation noise is constructed over data, which is obtained by running the simulation at strategically placed design points within the trust-region around the current iterate. Then the simulation components of the objective are replaced by their approximations and this analytical and closed-form optimization problem is solved to find the next iterate within the trust-region. Then the trust region is moved and the iterations continue until a satisfactory convergence criteria is met.

KRISHNAMOORTHY, MOHAN↗

NASA NTRS1988

Recursive inversion of externally defined linear systems

The approximate inversion of an internally unknown linear system, given by its impulse response sequence, by an inverse system having a finite impulse response, is considered. The recursive least squares procedure is shown to have an exact initialization, based on the triangular Toeplitz structure of the matrix involved. The proposed approach also suggests solutions to the problems of system identification and compensation.

Bach, Ralph E., Jr.↗

NASA NTRS2002

Locally-Based Kernal PLS Smoothing to Non-Parametric Regression Curve Fitting

We present a novel smoothing approach to non-parametric regression curve fitting. This is based on kernel partial least squares (PLS) regression in reproducing kernel Hilbert space. It is our concern to apply the methodology for smoothing experimental data where some level of knowledge about the approximate shape, local inhomogeneities or points where the desired function changes its curvature is known a priori or can be derived based on the observed noisy data. We propose locally-based kernel PLS regression that extends the previous kernel PLS methodology by incorporating this knowledge. We compare our approach with existing smoothing splines, hybrid adaptive splines and wavelet shrinkage techniques on two generated data sets.

Rosipal, Roman↗

NASA NTRS1978

The determination of gravity anomalies from geoid heights using the inverse Stokes' formula, Fourier transforms, and least squares collocation

A numerical method for the determination of gravity anomalies from geoid heights is described using the inverse Stokes formula. This discrete form of the inverse Stokes formula applies a numerical integration over the azimuth and an integration over a cubic interpolatory spline function which approximates the step function obtained from the numerical integration. The main disadvantage of the procedure is the lack of a reliable error measure. The method was applied on geoid heights derived from GEOS-3 altimeter measurements in the calibration area of the GEOS-3 satellite.

Rummel, R.↗

NASA NTRS1972

Zr-92(d,p)Zr-93 and Zr-92(d,t)Zr-91

The structures of Zr-93 and Zr-91 were studied by the stripping reaction Zr-92(d,p)Zr-93 and the pick-up reaction Zr-92(d,t)Zr-91 using 13 MeV incident deuterons. The reaction product particles were detected by counter telescope. Typical spectra from the reactions were analyzed by a nonlinear least squares peak fitting program which included a background search. Spin and parity assignments to observed excited levels were made by comparing experimental angular distributions with distorted wave Born approximation calculations.

Baron, N.↗

NASA NTRS1972

Polynomial approximations of thermodynamic properties of arbitrary gas mixtures over wide pressure and density ranges

Computer programs for flow fields around planetary entry vehicles require real-gas equilibrium thermodynamic properties in a simple form which can be evaluated quickly. To fill this need, polynomial approximations were found for thermodynamic properties of air and model planetary atmospheres. A coefficient-averaging technique was used for curve fitting in lieu of the usual least-squares method. The polynomials consist of terms up to the ninth degree in each of two variables (essentially pressure and density) including all cross terms. Four of these polynomials can be joined to cover, for example, a range of about 1000 to 11000 K and 0.00001 to 1 atmosphere (1 atm = 1.0133 x 100,000 N/m sq) for a given thermodynamic property. Relative errors of less than 1 percent are found over most of the applicable range.

Allison, D. O.↗

NASA NTRS1985

The determination of physical and dynamical parameters of Pluto/Charon and binary asteroids by least-square formation of a matched filter for a time series of images. I - The operational theory

A theory is derived for the determination of the masses, radii, and orbital elements of the Pluto/Charon, or similar, system based on the prediction of an image distribution over space and time and its comparison with observation. The comparison may be ultimately through the theory of least squares or the application of a matched filter to the observations as a three-dimensional signal stream at an initial or intermediate state. The theory is an approximation correct to fifth order in the diameters of celestial bodies. The theory of astronomical seeing that is used is based on Kolmogorov turbulence in the long-exposure limit. The images must be photometric. Linear tracking errors that can be removed are preferable to either automatic or manual guiding, in the collection of candidate observations.

Wildey, R. L.↗

NASA NTRS1989

Recursive inversion of externally defined linear systems by FIR filters

The approximate inversion of an internally unknown linear system, given by its impulse response sequence, by an inverse system having a finite impulse response, is considered. The recursive least-squares procedure is shown to have an exact initialization, based on the triangular Toeplitz structure of the matrix involved. The proposed approach also suggests solutions to the problem of system identification and compensation.

Bach, Ralph E., Jr.↗

NASA NTRS1989

Tau approximation techniques for identification of coefficients in parabolic PDE

A variant of the Tau method, called the weak Tau method, is developed on the basis of the weak form of the PDE for use in least-squares parameter estimation; also presented is a suitable abstract convergence framework. The emphasis is on the theoretical framework that allows treatment of the weak Tau method when it is applied to a wide class of inverse problems, including those for diffusion-advection equations, the Fokker-Planck model for population dynamics, and damped beam equations. Extensive numerical testing of the weak Tau method has demonstrated that it compares quite favorably with existing methods.

Banks, H. T.↗

NASA NTRS2017

Unsteady Aerodynamic Force Sensing from Strain Data

aerodynamic load sensing↗

NASA NTRS1960

Determination of Satellite Orbits from Radar Data

An optimum method for determining satellite orbits from radar data is presented in this report. Offering a good combination of speed and accuracy, the method makes use of orbit inclination and orbit elements in the plane, and proceeds with a differential correction of the orbit elements. Rapid, accurate methods of computing orbit elements are required to predict satellite positions for acquisition by other radars at points later along the orbit. In some cases the data are limited to a single pass over the observing station. The dynamical method is described in detail, and its accuracy is compared with those of two other methods: the purely geometrical, and the least-squares geometrical. By this optimum method the computing time, including the differential correction time, is 1 minute. Without differential correction, the rough determination takes from 10 to 20 seconds with approximately 5 miles positional uncertainty.

Cahill, W. F.↗