SEARCH · Engineering Papers

Results for “LEAST SQUARE FIT”

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

All sources NASA NTRS DOE OSTI

NASA NTRS2019

Retrieval of Snow Freeboard of Antarctic Sea Ice Using Waveform Fitting of CryoSat-2 Returns

In this paper we develop a CryoSat-2 algorithm to retrieve the surface elevation of the air–snow interface over Antarctic sea ice. This algorithm utilizes a two-layer physical model that accounts for scattering from a snow layer atop sea ice as well as scattering from below the snow surface. The model produces waveforms that are fit to CryoSat-2 level 1B data through a bounded trust region least-squares fitting process. These fit waveforms are then used to track the air–snow interface and retrieve the surface elevation at each point along the CryoSat-2 ground track, from which the snow freeboard is computed. To validate this algorithm, we compare retrieved surface elevation measurements and snow surface radar return power levels with those from Operation IceBridge, which flew along a contemporaneous CryoSat-2 orbit in October 2011 and November 2012. Average elevation differences (standard deviations) along the flight lines (IceBridge Airborne Topographic Mapper, ATM – CryoSat-2) are found to be 0.016 cm (29.24 cm) in 2011 and 2.58 cm (26.65 cm) in 2012. The spatial distribution of monthly average pan-Antarctic snow freeboard found using this method is similar to what was observed from NASA's Ice, Cloud, and land Elevation Satellite (ICESat), where the difference (standard deviation) between October 2011–2017 CryoSat-2 mean snow freeboard and spring 2003–2007 mean freeboard from ICESat is 1.92 cm (9.23 cm). While our results suggest that this physical model and waveform fitting method can be used to retrieve snow freeboard from CryoSat-2, allowing for the potential to join laser and radar altimetry data records in the Antarctic, larger (∼30 cm) regional differences from ICESat and along-track differences from ATM do exist, suggesting the need for future improvements to the method. Snow–ice interface elevation retrieval is also explored as a potential to obtain snow depth measurements. However, it is found that this retrieval method often tracks a strong scattering layer within the snow layer instead of the actual snow–ice interface, leading to an overestimation of ice freeboard and an underestimation of snow depth in much of the Southern Ocean but with promising results in areas such as the East Antarctic sector.

Fons, Steven W.↗

NASA NTRS1997

Determination of the Charon/Pluto Mass Ratio from Center-of-Light Astrometry

The Charon/Pluto mass ratio is a fundamental but poorly known parameter of the two-body system. Previous values for the mass ratio have ranged from 0.0837 plus or minus 0.0147 (Null et al., 1993, Astron. J. 105, 2319-2335) to 0.1566 plus or minus 0.0035 (Young et al., 1994, Icarus 108,186-199). We report here a new determination of the Charon/Pluto mass ratio, using five sets of groundbased images taken at four sites in support of Pluto occultation predictions. Unlike the Null et al. and Young et A determinations, where the centers of light for Pluto and Charon could be determined separately, this technique examines the motion of the center of light of the blended Pluto-Charon image. We compute the offsets of the observed center-of-light position of Pluto-Charon from the ephemeris position of the system and fit these offsets to a model of the Pluto-Charon system. The least-squares fits to the five data sets agree within their errors, and the weighted mean mass ratio is 0.117 plus or minus 0.006. The effects of errors in the Charon light fraction, semimajor axis, and ephemeris have been examined and are equal to only a small fraction of the formal error from the fit. This result is intermediate between those of Null et al., and Young et al. and matches a new value of 0.124 plus or minus 0.008 by Null and Owen (1996, Astron. J. 111, 1368-1381). The mass ratio and resulting individual masses and densities of Pluto and Charon are consistent with a collisional origin for the Pluto-Charon system.

Foust, Jeffrey A.↗

NASA NTRS1975

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

NASA NTRS1974

The 6300 A O/1-D/ airglow and dissociative recombination

Measurements of night-time 6300 A airglow intensities at the Arecibo Observatory have been compared with dissociative recombination calculations based on electron densities derived from simultaneous incoherent backscatter measurements. The agreement indicates that the nightglow can be fully accounted for by dissociative recombination. The comparisons are examined to determine the importance of quenching, heavy ions, ionization above the F-layer peak, and the temperature parameter of the model atmosphere. Comparable fits between the observed and calculated intensities are found for several available model atmospheres. The least-squares fitting process, used to make the comparisons, produces comparable fits over a wide range of combinations of neutral densities and of reaction constants. Yet, the fitting places constraints upon the possible combinations; these constraints indicate that the latest laboratory chemical constants and densities extrapolated to a base altitude are mutually consistent.

Wickwar, V. B.↗

NASA NTRS1972

Zeeman effect of P II.

An investigation was conducted of the Zeeman effect of PII because it provides unique J-value assignments together with a sensitive check on the accuracy of wave functions obtained from least-squares fitting of energy levels. About forty new lines were added to the spectrum, and g values were obtained for 76 levels of PII. By incorporating this new experimental material in a systematic program of least-squares calculations, supplemented by ab initio calculations, it was possible to complete the description of the 3p5p and 3p5d configurations.

Li, H.↗

NASA NTRS1981

Simulation of the magnetic structure of the inner heliosphere by means of non-spherical source surface

A new method for mapping the Sun's magnetic field B from the photosphere through the corona and interplanetary space is presented. The method entails the derivation of B from a scalar potential within a current-free annular volume bounded inside by photosphere and outside by a prescribed nonspherical source surface to which B is made (as nearly as possible) perpendicular. As usual we obtain the potential for the part of B that arises from currents inside the Sun by fitting an expansion to the observed line-of-sight component of B at the photosphere. A second least-squares fit is introduced to obtain the part of B that arises from currents outside the source surface. Comparisons are made between this model and observed coronal and interplanetary structures. There is evidence that observation data underestimate the strength of photospheric polar magnetic fields.

Levine, R. H.↗

NASA NTRS1990

Numerical experiments in geomagnetic modeling

Numerical tests were made, using least squares fitting of a spherical harmonic model, to a selection of Magsat data to determine the practical limits of this technique with modern computers. The resulting (M102189) model, whose coefficients were adjusted up to n = 50, was compared with M07AV6, a previous model which used least squares (on vector data) for coefficients up to n = 29, and Gauss-Legendre quadrature (on Z residuals) to adjust the coefficients up to n = 63. For the new least squares adjustment to n = 50 a condition number of 115 was obtained for the solution matrix, with a resulting precision of 11 significant figures. The M102189 model shows a lower and more Gaussian residual distribution than did M07AV6, though the Gaussian envelope fits to the residual distributions, even for the scalar field, gives "standard deviations' never lower than 6 nT, a factor of three higher than the estimated Magsat observational errors. Ionospheric currents are noted to have a significant effect on the coefficients of the internal potential functions.

Cain, Joseph C.↗

NASA NTRS2007

Numerical Technique for Analyzing Rotating Rake Mode Measurements in a Duct With Passive Treatment and Shear Flow

A technique is presented for the analysis of measured data obtained from a rotating microphone rake system. The system is designed to measure the interaction modes of ducted fans. A Fourier analysis of the data from the rotating system results in a set of circumferential mode levels at each radial location of a microphone inside the duct. Radial basis functions are then least-squares fit to this data to obtain the radial mode amplitudes. For ducts with soft walls and mean flow, the radial basis functions must be numerically computed. The linear companion matrix method is used to obtain both the eigenvalues of interest, without an initial guess, and the radial basis functions. The governing equations allow for the mean flow to have a boundary layer at the wall. In addition, a nonlinear least-squares method is used to adjust the wall impedance to best fit the data in an attempt to use the rotating system as an in-duct wall impedance measurement tool. Simulated and measured data are used to show the effects of wall impedance and mean flow on the computed results.

Dahl, Milo D.↗

DOE OSTI2021

Parameter inference from event ensembles and the top-quark mass

One of the key tasks of any particle collider is measurement. In practice, this is often done by fitting data to a simulation, which depends on many parameters. Sometimes, when the effects of varying different parameters are highly correlated, a large ensemble of data may be needed to resolve parameter-space degeneracies. An important example is measuring the top-quark mass, where other physical and unphysical parameters in the simulation must be profiled when fitting the top-quark mass parameter. We compare four different methodologies for top-quark mass measurement: a classical histogram fit similar to one commonly used in experiment augmented by soft-drop jet grooming; a 2D profile likelihood fit with a nuisance parameter; a machine-learning method called DCTR; and a linear regression approach, either using a least-squares fit or with a dense linearly-activated neural network. Despite the fact that individual events are totally uncorrelated, we find that the linear regression methods work most effectively when we input an ensemble of events sorted by mass, rather than training them on individual events. Although all methods provide robust extraction of the top-quark mass parameter, the linear network does marginally best and is remarkably simple. For the top study, we conclude that the Monte-Carlo-based uncertainty on current extractions of the top-quark mass from LHC data can be reduced significantly (by perhaps a factor of 2) using networks trained on sorted event ensembles. More generally, machine learning from ensembles for parameter estimation has broad potential for collider physics measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1994

Vega: A rapidly rotating pole-on star

High-dispersion (2.4 A/mm), ultrahigh signal-to-noise ratio (3000:1) Reticon spectra of Vega revealed two distinct types of profiles. The strong lines exhibit classical rotational profiles with enhanced wings, but the weak lines have distinctly different, flat-bottomed profiles. Using ATLAS9 model atmopheres and SYNTHE synthetic spectra, Vega has been modeled as a rapidly rotating, pole-on star with a gradient in temperature and gravity over the photosphere. By fitting to the flat-bottomed line profiles of Fe 1 lambda 4528 and Ti 2 lambda 4529, we find least-squares fit values of V sin i = 21.8 plus or minus 0.2 km/sec polar T(sub eff) = 9695 plus or minus 25 K, polar log(base 10)g = 3.75 plus or minus 0.02 dex, V(sub eq) = 245 plus or minus 15 km/sec, and inclination 5 deg .1 plus or minus 0 deg .3. The variations in T(sub eff) and log(base 10)g over the photosphere total 390 K and 0.08 dex, respectively. Assuming V sin i = 21.8 km/sec, an independent fit to the observed continuous flux from 1200 to 10,500 A produced a similar set of values with polar T(sub eff) = 9595 plus or minus 20 K, polar log(base 10)g = 3.80 plus or minus 0.03 dex, and inclination 6 deg .0 plus or minus 0 deg .7.

Gulliver, Austin F.↗

NASA NTRS1967

Computer program for network synthesis by frequency response fit

Computer program synthesizes a passive network by minimizing the difference in desired and actual frequency response. The program solves for the critical points of the error function /weighted least squares fit between calculated and desired frequency response/ by the multivariable Newton-Raphson method with components constrained to an admissible region.

Green, S.↗

NASA NTRS1994

Year-round measurements of ozone at 66 deg S with a visible spectrometer

In March 1990, a zenith-sky UV-visible spectrometer of the design 'Systeme Automatique d'Obervation Zenithal' (SAOZ) was installed at Faraday in Antarctica (66.3 deg S, 64.3 deg W). SAOZ records spectra between 290 and 600 nm during daylight. Its analysis program fits laboratory spectra of constituents, at various wavelengths, to the differential of the ratio of the observed spectrum and a reference spectrum. The least-squares fitting procedure minimizes the sum-of-squares of residuals. Ozone is deduced from absorption in its visible bands between 500 and 560 nm. The fortunate colocation of this SAOZ with the well-calibrated Dobson at Faraday has allowed us to examine the calibration of the zero of the SAOZ, difficult at visible wavelengths because of the small depth of absorption. Here we describe recent improvements and limitations to this calibration, and discuss SAOZ measurements of ozone during winter in this important location at the edge of the Antarctic vortex.

Roscoe, Howard K.↗

NASA NTRS2014

An Improved Cryosat-2 Sea Ice Freeboard Retrieval Algorithm Through the Use of Waveform Fitting

We develop an empirical model capable of simulating the mean echo power cross product of CryoSat-2 SAR and SAR In mode waveforms over sea ice covered regions. The model simulations are used to show the importance of variations in the radar backscatter coefficient with incidence angle and surface roughness for the retrieval of surfaceelevation of both sea ice floes and leads. The numerical model is used to fit CryoSat-2 waveforms to enable retrieval of surface elevation through the use of look-up tables and a bounded trust region Newton least squares fitting approach. The use of a model to fit returns from sea ice regions offers advantages over currently used threshold retrackingmethods which are here shown to be sensitive to the combined effect of bandwidth limited range resolution and surface roughness variations. Laxon et al. (2013) have compared ice thickness results from CryoSat-2 and IceBridge, and found good agreement, however consistent assumptions about the snow depth and density of sea ice werenot used in the comparisons. To address this issue, we directly compare ice freeboard and thickness retrievals from the waveform fitting and threshold tracker methods of CryoSat-2 to Operation IceBridge data using a consistent set of parameterizations. For three IceBridge campaign periods from March 20112013, mean differences (CryoSat-2 IceBridge) of 0.144m and 1.351m are respectively found between the freeboard and thickness retrievals using a 50 sea ice floe threshold retracker, while mean differences of 0.019m and 0.182m are found when using the waveform fitting method. This suggests the waveform fitting technique is capable of better reconciling the seaice thickness data record from laser and radar altimetry data sets through the usage of consistent physical assumptions.

CryoSat-2↗

DOE OSTI2020

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2001

On-Orbit Performance of Autonomous Star Trackers

This paper presents the results of a performance study of the autonomous star trackers (ASTs) on the IMAGE and the EO-1 spacecraft. IMAGE is a spinning spacecraft without gyros or redundant precision attitude sensors, so the statistical properties of the AST are estimated simply by comparing the output observed quaternions with a rigid rotator model with constant angular momentum. The initial conditions are determined by a least-squares fit to minimize the AST residuals. An additional fit is used to remove the remaining systematic error and to obtain the inherent sensor noise. Gyro rate data are available for the EO-1 mission, so the AST noise statistics are obtained from the residuals after solving for an epoch attitude and gyro bias also using a least-squares method.

Airapetian, V.↗

NASA NTRS1999

Optimized System Identification

In system identification, one usually cares most about finding a model whose outputs are as close as possible to the true system outputs when the same input is applied to both. However, most system identification algorithms do not minimize this output error. Often they minimize model equation error instead, as in typical least-squares fits using a finite-difference model, and it is seen here that this distinction is significant. Here, we develop a set of system identification algorithms that minimize output error for multi-input/multi-output and multi-input/single-output systems. This is done with sequential quadratic programming iterations on the nonlinear least-squares problems, with an eigendecomposition to handle indefinite second partials. This optimization minimizes a nonlinear function of many variables, and hence can converge to local minima. To handle this problem, we start the iterations from the OKID (Observer/Kalman Identification) algorithm result. Not only has OKID proved very effective in practice, it minimizes an output error of an observer which has the property that as the data set gets large, it converges to minimizing the criterion of interest here. Hence, it is a particularly good starting point for the nonlinear iterations here. Examples show that the methods developed here eliminate the bias that is often observed using any system identification methods of either over-estimating or under-estimating the damping of vibration modes in lightly damped structures.

Juang, Jer-Nan↗

NASA NTRS1992

Polynomial coefficients for calculating O2 Schumann-Runge cross sections at 0.5/cm resolution

O2 cross sections from 49,000 to 57,000/cm have been fitted with temperature dependent polynomial expressions, providing an accurate and efficient means of determining Schumann-Runge band cross sections for temperatures between 130 and 500 K. The least squares fits were carried out on a 0.5/cm spectral grid, using cross sections obtained from a Schumann-Runge line-by-line model that incorporates the most recent spectroscopic data. The O2 cross sections do not include the underlying Herzberg continuum, but they do contain contributions from the temperature dependent Schumann-Runge continuum. The cross sections are suitable for use in UV transmission calculations at high spectral resolution. They should also prove useful for updating existing parameterizations of ultraviolet transmission and O2 photolysis.

Minschwaner, K.↗

NASA NTRS1982

Line positions and strengths of HDO in the 2400-3300 per cm region

The vibration-rotation spectrum of HDO in the 2400-3300 per cm region has been studied at 0.01 per cm resolution with a Fourier transform spectrometer. Values of the ground state levels and levels in the (100) and (020) states have been determined from the data. The ground state levels are fitted to a Watson-type Hamiltonian from which 17 rotational constants are obtained. Measurements of the strengths of over 1100 lines have been made using a nonlinear least-squares fitting technique, and 548 of these lines are analyzed to determine the band strengths, the coefficients of the F factors, and the Coriolis interaction parameters. The strengths of the A- and B-type (100) bands and the A- and B-type (020) bands are 53.2 + or - 0.9, 0.43 + or - 0.07, 6.00 + or - 0.54, and 2.75 + or - 0.12 per sq cm per atm at 296 K, respectively.

Toth, R. A.↗