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 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS1981

The Effect of Sea-Surface Sun Glitter on Microwave Radiometer Measurements

A relatively simple model for the microwave brightness temperature of sea surface Sun glitter is presented. The model is an accurate closeform approximation for the fourfold Sun glitter integral. The model computations indicate that Sun glitter contamination of on orbit radiometer measurements is appreciable over a large swath area. For winds near 20 m/s, Sun glitter affects the retrieval of environmental parameters for Sun angles as large as 20 to 25 deg. The model predicted biases in retrieved wind speed and sea surface temperature due to neglecting Sun glitter are consistent with those experimentally observed in SEASAT SMMR retrievals. A least squares retrieval algorithm that uses a combined sea and Sun model function shows the potential of retrieving accurate environmental parameters in the presence of Sun glitter so long as the Sun angles and wind speed are above 5 deg and 2 m/s, respectively.

Wentz, F. J.↗

NASA NTRS1990

Loop current eddy shedding estimated using Geosat altimeter data

A method for extracting periodic signals from altimeter data is applied in the Gulf of Mexico using data from the first two years of the Geosat Exact Repeat Mission. The least squares fitting of selected frequencies of sine and cosine functions to the time series of sea surface height at each point of the Geosat ground track is used to estimate the dominant frequency in the gulf for this time period. Optimal interpolation of the sine and cosine coefficients at the dominant frequency to a uniformly spaced grid allows reconstruction of the time series in the gulf and shows eddies shedding from the Loop Current with a period of approximately 10.5 months, consistent with shedding periods determined from numerical models and analysis of satellite infrared data. This technique is a novel method for investigation of periodic signals in Geosat altimetry and is not severely affected by the data outages present in the time series.

Jacobs, Gregg A.↗

DOE OSTI2020

Algorithms in Diffraction Profile Analysis

This chapter introduces several approaches to crystallographic structure refinement from diffraction profile data including the Rietveld method, global optimization methods, and Bayesian modeling. First, the Rietveld method and related nonlinear least squares approaches to analyzing diffraction profiles are briefly discussed. We then review the most common algorithms within global optimization as alternative methods of structure solution: grid search, simulated annealing, and genetic algorithms. As the Rietveld method is the standard tool for analyzing diffraction profile data, we catalogue its challenges and limitations. A Bayesian approach to diffraction profile analysis is presented as an alternative, along with several sampling algorithms including Gibbs sampling, variants of Metropolis sampling, Hamiltonian Monte Carlo, and approximate Bayesian computing. These sampling algorithms are applied to the fitting of a neutron diffraction profile from a National Institute of Standards and Technology silicon standard reference material, and their results are compared.

Singh, Susheela↗

NASA NTRS1983

Rotational CARS application to simultaneous and multiple-point temperature and concentration determination in a turbulent flow

Coherent anti-Stokes Raman scattering (CARS) from the pure rotational Raman lines of N2 is employed to measure the instantaneous (approximately 10 ns) rotational temperature of N2 gas at room temperature and below with good spatial resolution (0.2 x 0.2 x 3.0 cu mm). A broad bandwidth dye laser is used to obtain the entire rotational spectrum from a single laser pulse; the CARS signal is then dispersed by a spectrograph and recorded on an optical multichannel analyzer. A best fit temperature is found in several seconds with the aid of a computer for each experimental spectrum by a least squares comparison with calculated spectra. The model used to calculate the theoretical spectra incorporates the temperature and pressure dependence of the pressure-broadened rotational Raman lines, includes the nonresonant background susceptibility, and assumes that the pump laser has a finite linewidth. Temperatures are fit to experimental spectra recorded over the temperature range of 135 to 296 K, and over the pressure range of .13 to 15.3 atm.

Snow, J. B.↗

NASA NTRS1990

A DFW microburst model based on DL-191 data

Analysis of the August 2, 1985 crash of an L-1011 jumbo jet (DL-191), on approach to the Dallas-Ft. Worth International Airport (DFW) in a thunderstorm indicates that the severe windshear microburst which caused the crash was composed not only of a strong downflow and outflow but also included several large-scale vortex rings entrained in the flowfield. A detailed two-dimensional model of the DFW microburst, based on data from DL-191, is presented. The model uses wind vector and flight path data reconstructed by NASA-Ames Research Center from the L-1011's digital flight data records and ground-based radar measurements. The model was developed using a combination of interactive graphics and a least-square-error best fit between the modeled and measured wind vectors along the DL-191 flight path. The model confirms a microburst structure based on a von Karman vortex street rather than on a Kelvin-Helmholtz instability and indicates that the L-1011 flew near or directly through several strong vortices. In addition, the results also confirm that the reconstructed wind vector data contain a time lag in the horizontal winds. For the DL-191 this lag is approximately three seconds.

Grantham, Walter J.↗

NASA NTRS2001

Properties of Stratospheric Aerosol Estimated from HALOE Data

Extinction coefficients for stratospheric aerosols at 8 HALOE (HALogen Occultation Experiment) wavelengths are determined by comparing transmittances data for two adjacent solar occultation measurements, where one limb path is loaded with aerosols but the other path is free of aerosols. These extinction coefficients are used to infer the aerosol properties such as composition and size distribution parameters. Mie theory has been used to calculate the extinction coefficients, and a nonlinear least square method is applied to determine the aerosol properties. Sixteen cases are selected for the retrieval in southern hemisphere at latitudes from 21 to 48 deg S for the period of 29 Mar - 31 May 1992. Retrieved size width ranges from 1.1 to 1.5 and radius ranges from 0.25 to 0.45 micrometers. These size parameters are within the ranges of in situ measurements at Laramie, Wyoming. Retrieved weight percent of H2SO4 is larger than the equilibrium value by about 5 to approximately 10 weight percent, similar to the results for northern hemisphere at latitudes 20 to 55 deg N for the period from Nov 1991 to Feb. 1992.

Lee, Kwang-Mog↗

NASA NTRS2001

A Combined HIPPARCOS and Multichannel Astrometric Photometer Study of the Proposed Planetary System of Rho Coronae Borealis

Hipparcos and Multichannel Astrometric Photometer (MAP) observations of rho Coronae Borealis independently display astrometric motion at the period of the proposed extrasolar planetary companion to the star. Individual least-squares fits to each astrometric data set yield independent estimates of the semimajor axis, inclination, and node angle that are in excellent agreement. A combined solution of the Hipparcos and MAP data yields an inclination of 0.5 deg, a node at 30.5 +/- 12.4, and a semimajor axis of 1.66 +/- 0.35 mas, indicating a companion mass of 0.14 +/- 0.05 solar masses over two orders of magnitude greater than the minimum mass for the companion as determined by radial velocity studies. This mass is approximately that of an M dwarf star, the companion cannot be a planetary object.

Gatewood, George↗

NASA NTRS2019

GOES-17 Magnetometer On-Orbit Calibration Design and Results

The Geostationary Operational Environmental Satellite-R (GOES-R) is a series of four NOAA weather satellites. Each spacecraft flies a Magnetometer instrument which includes two boom mounted fluxgate sensors located 6.3 and 8.5 meters from the spacecraft. Sensor readings are corrupted by both internal (bias and scale factor) and external (spacecraft emitted fields and misalignments) error sources. In order to estimate the bias, an on-orbit calibration is performed, slewing the spacecraft about multiple axes and then performing a fit of the data to determine zero offset and misalignment terms. The first satellite in the series, GOES-16 was launched in November 2016. Its magnetometer calibration maneuver was found to be inadequate providing insufficient rotation about orthogonal axes. The uncertainty in the estimated biases was found to be in the 3-5 nanotesla range, exceeding the magnetometer accuracy requirement of 1.0 nanotesla. Based on lessons learned, the maneuver was redesigned to provide multiple 360 degree rotations about two orthogonal axes, providing optimal visibility into bias and misalignment. Varying ambient fields over the duration of the calibration maneuvers, corrupts the ability to accurately determine parameters of interest. To address this issue, rather than assuming a constant field for the duration of the maneuver, a non-linear least squares filter was employed that uses a series of spline functions to better fit the ambient field. The number of splines can be varied to minimize the least squares filter residuals. Simulations were performed using 360 days of GOES-16 data to bound the uncertainty in the biases for given residuals of the fit to the ambient field. Following the launch of the second satellite in the series, GOES-17 in March 2018, the new calibration maneuver and enhanced bias estimation algorithm were exercised twice. Based on the residuals of the fit, the uncertainty of the estimated zero offset was reduced by an order of magnitude from 3-5 nanotesla to approximately 0.4 nanotesla per axis. This paper provides details of the maneuver design and the innovative algorithm used to estimate the GOES-17 biases and misalignments.

GOES-17↗

NASA NTRS1983

Spline approximation of quantile functions

The study reported here explored the development and utility of a spline representation of the sample quantile function of a continuous probability distribution in providing a functional description of a random sample and a method of generating random variables. With a spline representation, the random samples are generated by transforming a sample of uniform random variables to the interval of interest. This is useful, for example, in simulation studies in which a random sample represents the only known information about the distribution. The spline formulation considered here consists of a linear combination of cubic basis splines (B-splines) fit in a least squares sense to the sample quantile function using equally spaced knots. The following discussion is presented in five parts. The first section highlights major results realized from the study. The second section further details the results obtained. The methodology used is described in the third section, followed by a brief discussion of previous research on quantile functions. Finally, the results of the study are evaluated.

Schiess, J. R.↗

DOE OSTI2023

Probabilistic partition of unity networks for high–dimensional regression problems

We explore the probabilistic partition of unity network (PPOU-Net) model in the context of high-dimensional regression problems and propose a general framework focusing on adaptive dimensionality reduction. With the proposed framework, the target function is approximated by a mixture of experts model on a low-dimensional manifold, where each cluster is associated with a fixed-degree polynomial. We present a training strategy that leverages the expectation maximization (EM) algorithm. During the training, we alternate between (i) applying gradient descent to update the DNN coefficients; and (ii) using closed-form formulae derived from the EM algorithm to update the mixture of experts model parameters. Under the probabilistic formulation, step (ii) admits the form of embarrassingly paralleliazable weighted least-squares solves. The PPOU-Nets consistently outperform the baseline fully-connected neural networks of comparable sizes in numerical experiments of various data dimensions. Here, we also explore the proposed model in applications of quantum computing, where the PPOU-Nets act as surrogate models for cost landscapes associated with variational quantum circuits.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Enrichment of heavy nuclei in He-3-rich flares

IMP-7 observations of five solar-flare particle events characterized by He-3 enrichment are reported which show that such events are also enriched in nuclei with charges (Z) of at least 6. The ratio of these nuclei to H-1 at approximately 3 MeV/nucleon was found to be enriched by about 10 to 100 times, while the ratio with respect to He-4 was enriched by about 3 to 30 times. It is suggested that the simultaneous enhancement of He-3 and the heavier nuclei as well as the absence of H-2 and H-3 during the observed events may be partly due to a preferential acceleration process which depends on the ratio of the square of the charge to the atomic weight of the nuclei.

Hurford, G. J.↗

DOE OSTI2026

Validation and Verification of Python based Neutron Spectrum Unfolding Software

To validate and verify the python-based code (PySL), designed to replicate the programs used by STAYSL for Beam Correction Factor (BCF) and Self-Shielding Factor (SHIELD), a series of tests were performed. To test BCF a python script was written to generate a random flux history file and both versions of the code processed the data. The test verified matching values up to at least one decimal place, approximately 10,000 tests where run and each one passed. Isotopes began to fail the tests once neutron saturation was reached. To verify this the total time of exposure was varied the isotopes that failed were compared to a list of their half-lives. The test process for SHIELD was very similar but, in this case, the code began by producing an input file with varying thickness and device type/environment for the SHIELD input. The failure condition for this test was if any of the data points for an isotope had a difference above 3%. Approximately 40 of these tests were run and there were only 3 isotopes that had reoccurring failures but only 2% of their points were above the 3% difference. A visual comparison was conducted by plotting the results from both programs. Although the test failed, the differences between their values were minuscule, and the self-shielding factor’s shape was preserved when plotted. Next steps for this project will be validating and verifying the python-based SigPhi code and then reproducing and testing the least squares unfolding performed by STAYSL.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Network-Wide Traffic Signal Control Using Bilinear System Modeling and Adaptive Optimization

This study proposes a new multi-input multi-output optimal bilinear signal control method in which a bilinear dynamic model approximation is used to capture the nonlinear dynamics of the urban traffic networks. With signal green time splits as the control input and traffic delay changes as the output for each intersections in the network, a bilinear system model was developed, which, on the basis of linear system modeling, takes interactions among traffic delays and signal timing splits into consideration. Based on the bilinear system modeling framework, we conducted two steps in each time interval to derive traffic control strategies: (1) we used the normalized least-squared algorithm to estimate system parameters; and (2) we solved an online optimization problem to obtain the updated traffic control inputs for the signal timing that minimizes future traffic delays. We evaluated the proposed method in a microscopic traffic simulation environment (VISSIM) with a 35-intersection network of Bellevue city in Washington. Two different traffic demand patterns: (1) normal traffic demands; and (2) time-varying traffic demands were simulated to compare the performance of different control strategies. Experimental results show that (1) the proposed bilinear system model can better describe traffic system dynamics than linear-model based methods, such as our previously developed linear-quadratic regulator control; and (2) the proposed method outperforms the state-of-the-art signal control strategies, namely the max-pressure and the self-organizing traffic light control methods. We have also shown that the proposed method is applicable to all other possible network layouts and signal controller phasing structures.

42 ENGINEERING↗

NASA NTRS1978

The prediction and mapping of geoidal undulations from GEOS-3 altimetry

From the adjusted altimeter data an approximation to the geoid height in ocean areas is obtained. Methods are developed to produce geoid maps in these areas. Geoid heights are obtained for grid points in the region to be mapped, and two of the parameters critical to the production of an accurate map are investigated. These are the spacing of the grid, which must be related to the half-wavelength of the altimeter signal whose amplitude is the desired accuracy of the contour; and the method adopted to predict the grid values. Least squares collocation was used to find geoid undulations on a 1 deg grid in the mapping area. Twenty maps, with their associated precisions, were produced and are included. These maps cover the Indian Ocean, Southwestern and Northeastern portions of the Pacific Ocean, and Southwest Atlantic and the U.S. Calibration Area.

Kearsley, W.↗

NASA NTRS1974

Post-flight differential correction analysis using Vinti's spheroidal method for the small astronomy satellite orbit

The results of an intensive analysis of a differential orbit improvement method utilizing observational data for a 550-kilometer altitude, near-circular, near-equatorial satellite orbit are presented. Observations of the Small Astronomy Satellite (SAS-1) are in the form of direction cosines as measured at two ground interferometer tracking stations near the Equator during the first 22 orbital revolutions (approximately 37 hours) after launch of the spacecraft. Numerical results, in both tabular and graphical form, are displayed for numerous iterated fittings of various observational arcs by differential correction of the orbital elements. Parameters varied in these comparative cases include the time duration of the observational data block, the number of pairs of direction cosine data and the number of tracking station passes included in the solution, the distribution of such passes between the two available tracking stations, and the acceptance criterion for the observational residuals in the least squares fitting procedure. It was found that three observational pairs of direction cosine data, the minimum number possible for a uniquely determined solution in theory, are sufficient to promote covergence to an accurate solution, if properly selected.

Walden, H.↗

NASA NTRS1985

Effect of enclosure shape on natural convection velocities in microgravity

Natural convection velocities within two-dimensional enclosures of various shapes are estimated with reasonable accuracy, for Rayleigh numbers up to at least 5000, by considering the area to be equivalent to a circle of equal area. For the class of figures including the square, upright half-circle, and upright 2 x 1 aspect ratio rectangle, excellent agreement is obtained by considering the temperature difference across the maximum horizontal dimension to be equal to that across the equal area circle horizontal diameter. The agreement of two upright oblong shapes indicates probable agreement for any similar upright oblong shape of roughly the same aspect ratio. The horizontal oblong shapes have natural convection velocities approximately one-half that estimated based on the equal area circle, with the temperature difference across the maximum horizontal dimension taken to be equal to that across the equal area circle diameter.

Robertson, S. J.↗

NASA NTRS1972

A computer program for Moessbauer data processing

A computer program to analyze Mossbauer data is presented in detail. The least-squares curve fitting techniques described apply to single line spectra, single hyperfine spectra, or when the constituent spectra are separated well enough to let the individual absorption peaks stand alone. The present program is not adapted for complex spectra resulting from the existence of several local environments in the absorber iron alloy. Sample problems are presented to aid the user in setting up and running the program. The program is written in FORTRAN 4 language for the Control Data 6000 series digital computer with the SCOPE 3.0 operating system and requires approximately 115,000 octal locations of core storage. A typical case with one absorption peak runs in 20 seconds, and a typical problem with six absorption peaks requires 50 seconds.

Howser, L. M.↗

NASA NTRS1999

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↗