Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

Abridged spectral matrix inversion: parametric fitting of X-ray fluorescence spectra following integrative data reduction

Recent improvements in both X-ray detectors and readout speeds have led to a substantial increase in the volume of X-ray fluorescence data being produced at synchrotron facilities. This in turn results in increased challenges associated with processing and fitting such data, both temporally and computationally. Herein an abridging approach is described that both reduces and partially integrates X-ray fluorescence (XRF) data sets to obtain a fivefold total improvement in processing time with negligible decrease in quality of fitting. The approach is demonstrated using linear least-squares matrix inversion on XRF data with strongly overlapping fluorescent peaks. This approach is applicable to any type of linear algebra based fitting algorithm to fit spectra containing overlapping signals wherein the spectra also contain unimportant (non-characteristic) regions which add little (or no) weight to fitted values, e.g. energy regions in XRF spectra that contain little or no peak information.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

On-line identification and attitude control for SCOLE

This paper documents on-line linear least-square identification and attitude control of SCOLE, a laboratory apparatus representing an offset-feed antenna attached to the Space Shuttle. Identification is done autonomously by starting a slew maneuver in pitch or roll with reaction jets and observing the time history data of associated Euler angles when the jets are quiescent. Linear least-square analysis is used to select the parameters that best fit the output of an Autoregressive (AR) model to the data. The control effectiveness of the jets is determined in a subsequent test, again using linear least squares. The parameters so derived are used to design switching lines for time-optimal attitude control. This report describes the identification and control algorithms and the experimental apparatus and procedures used. Also, experimental data are presented that reflect the performance of the identification algorithms and the attitude control system.

Montgomery, R. C.↗

Temporal Ni K-Edge X-ray Absorption Spectroscopy Study Reveals the Kinetics of the Ni Redox Behavior of the Iron-Nickel Oxide Bimetallic OER Catalyst

Operando X-ray absorption spectroscopy (XAS) can be utilized to probe the phase and structural changes of FeNiOx and similar transition metal oxide electrocatalysts during electrocatalytic reactions. However, capturing the temporal changes occurring in the operando chemistry of electrocatalysts has been little studied. Here, in this work, we successfully capture the time-resolved changes at the Ni K-edge for both the X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) regions of an aqueous phase synthesized FeNiO x bimetallic nanoparticle electrocatalyst. During a stepped voltage experiment, the temporal change in the Ni K-edge was observed as the applied voltage was stepped from 0.7 to 0.8 V versus the silver/silver chloride reference electrode, and the change observed is associated with the Ni redox transition from the 2+ to 3+/4+ oxidation state. Individual XAS spectra were obtained in 90 s, and a total of 9 scans post voltage step showed a unique transition from a hydroxide to an oxyhydroxide phase. The shift in absorption edge energy position and the changes in spectral shape for individual scans explain the typically observed broadening of the Ni K-edge XANES spectrum during time-averaged operando XAS, and a kinetic analysis revealed a first-order observed rate constant, |k obs |, of 0.00426 s –1 and a half-life (t 1/2 ) of 163 s. Multivariate curve resolution-alternating least squares analysis followed by linear combination fitting analysis for time-averaged versus time-resolved Ni K-edge changes upon voltage step show distinct differences in estimated contributions from Ni 2+ oxide/hydroxide phases. Detailed EXAFS modeling shows the phase transition from hydroxide to oxyhydroxide on top of the unchanged metallic core.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The complex magnetic field of Jupiter

An analysis of the characteristics of the magnetic field of the planet Jupiter is presented. The data were obtained during the flight of Pioneer 11 space probe, using a high field triaxial fluxgate magnetometer. The data are analyzed in terms of traditional Schmitt normalized spherical harmonic expansion fitted to the observations in a least squares sense. Tables of data and graphs are provided to summarize the findings.

Acuna, M. H.↗

Computer modeling of multiple-channel input signals and intermodulation losses caused by nonlinear traveling wave tube amplifiers

The multiple channel input signal to a soft limiter amplifier as a traveling wave tube is represented as a finite, linear sum of Gaussian functions in the frequency domain. Linear regression is used to fit the channel shapes to a least squares residual error. Distortions in output signal, namely intermodulation products, are produced by the nonlinear gain characteristic of the amplifier and constitute the principal noise analyzed in this study. The signal to noise ratios are calculated for various input powers from saturation to 10 dB below saturation for two specific distributions of channels. A criterion for the truncation of the series expansion of the nonlinear transfer characteristic is given. It is found that he signal to noise ratios are very sensitive to the coefficients used in this expansion. Improper or incorrect truncation of the series leads to ambiguous results in the signal to noise ratios.

Stankiewicz, N.↗

Strength, Width, and Pressure Shift Measurements of 54 Lines in the Oxygen A-Band

The absorption band of molecular oxygen, centered at 760]en1] nm, is the atmospheric absorber for the Differential Absorption Lidar (DIAL) systems used to measure atmospheric temperature, pressure, and density. To provide accurate line parameters for such systems, a careful spectroscopic study was made of the A-band, with measurements of line strengths, widths, pressure-induced frequency shifts, and collisional narrowing effects. The width and shift parameters were measured over a temperature range of -20 to 100 C so that the temperature dependence of these parameters can also be determined. To analyze the results, a least-squares fiting routine was written to fit standard line profiles to the observed profiles. These measurements, which include the first observations of pressure shifts and collisional narrowing in the band, are an important contribution to lidar system utilizing the A-band.

Ritter, K. J.↗

A method of calculating the total flow from a given sea surface topography

Using a simple dynamical model of a wind-driven ocean circulation of the Stommel type, and an analytical basis developed to objectively analyze the sea surface height residuals from an altimeter and, in the process, to determine the total flow instead of just the near surface geostrophic component associated with the given sea surface topography. The method is based on first deriving the solution to the forced problem for a given wind stress required to develop a hypothetical true or perfect data field and to establishing the basis for the objective analysis. The stream function and the surface height field for the forced problem are developed in terms of certain characteristic functions with the same expansion coefficients for both fields. These characteristic functions are simply the solutions for a homogeneous elliptic equation for the stream function and the solutions of an inhomogeneous balance equation for the height field. For the objective analysis, using a sample of randomly selected height values from the true data field, the height field characteristic functions are used to fit the given topography in a least squares sense. The resulting expansion coefficients then permit the synthesis of the total flow field via the stream function characteristic modes and the solution is perfectly well behaved even along the equator. The method of solution is easily adaptable to realistic ocean basis by straight forward numerical methods. The analytical basis of the theory and the results for an ideal rectangular basin on a beta plane are described.

Rao, Desiraju B.↗

Deconvolution of spectra for intimate mixtures

Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.

Mustard, John F.↗

Multiple boundary condition testing error analysis

Techniques for interpreting data from multiple-boundary-condition (MBC) ground tests of large space structures are developed analytically and demonstrated. The use of MBC testing to validate structures too large to stand alone on the ground is explained; the generalized least-squares mass and stiffness curve-fitting methods typically applied to MBC test data are reviewed; and a detailed error analysis is performed. Consideration is given to sensitivity coefficients, covariance-matrix theory, the correspondence between test and analysis modes, constraints and step sizes, convergence criteria, and factor-analysis theory. Numerical results for a simple beam problem are presented in tables and briefly characterized. The improved error-updating capabilities of MBC testing are confirmed, and it is concluded that reasonably accurate results can be obtained using a diagonal covariance matrix.

Glaser, R. J.↗

Magnetic field structure of interplanetary magnetic clouds at 1 AU

Interplanetary magnetic clouds emerge as a feature of the solar wind at 1 AU, exhibiting enhanced field strength and lower plasma temperature and density than the surrounding plasma. A least-squares program has been developed which fits magnetic field data within a cloud, while estimating such cloud properties as its size, maximum field strength, and axis inclination. The results obtained from a study of 12 clouds observed at 1 AU point to a probable cloud axis direction within 15 deg of the ecliptic plane and about 100 deg from the sun's direction, when projected into the ecliptic plane. A wide variety of orientations is observed; some extend to 80 deg from the ecliptic.

Lepping, R. P.↗

Are existing magnetospheric models excessively stretched?

The distribution of the average north-south magnetic field component Bz in the vicinity of the neutral sheet has been investigated. This component is crucial for mappings between the nightside polar ionosphere and the equatorial magnetosphere. Data sets consisting of about 0.5 R(E) averages of magnetic field observations by the IMP/HEOS and ISEE spacecraft have been compared to the field predicted by the Tsyganenko (1987, 1989) models T87 and T89. It was found that both T87 and T89 underestimate Bz in the near tail region by as much as a factor of 2. Modified versions of the T87 model, incorporating plasma sheet warping, were obtained by fitting the model parameters via nonlinear least squares to the ISEE data set and yielded Bz values in agreement with the ISEE data. The study reveals an enormous scatter among the observed baseline values of Bz (on a time scale of 10-20 min), as well as intrinsic biases imposed by the mathematical structure of tail models, and these two factors (especially the first one) greatly limit the accuracy of model predictions of tail Bz.

Peredo, Mauricio↗

Robust pupil center detection using a curvature algorithm

Determining the pupil center is fundamental for calculating eye orientation in video-based systems. Existing techniques are error prone and not robust because eyelids, eyelashes, corneal reflections or shadows in many instances occlude the pupil. We have developed a new algorithm which utilizes curvature characteristics of the pupil boundary to eliminate these artifacts. Pupil center is computed based solely on points related to the pupil boundary. For each boundary point, a curvature value is computed. Occlusion of the boundary induces characteristic peaks in the curvature function. Curvature values for normal pupil sizes were determined and a threshold was found which together with heuristics discriminated normal from abnormal curvature. Remaining boundary points were fit with an ellipse using a least squares error criterion. The center of the ellipse is an estimate of the pupil center. This technique is robust and accurately estimates pupil center with less than 40% of the pupil boundary points visible.

Non-NASA Center↗

Steady surficial core motions - An alternate method

A new method for deriving steady surficial core motions from geomagnetic field models is presented and applied. The original method determines the steady velocity field at the top of a frozen-flux core which best fits, in both the spatial and temporal linear least squares sense, a model of the large scale geomagnetic secular variation (SV) at the base of a source-free mantle (Voorhies, 1986). New method is based upon the same physical assumptions, but fits SV at earth's surface instead of the core-mantle boundary (CMB). The newly derived flow differs somewhat from prior solutions, and shows but slight westward drift; however, other key global properties are virtually unchanged. Over a 15 year interval the new method provides a better fit to the time-varying geomagnetic field at earth's surface than does the old.

Voorhies, Coerte V.↗

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↗

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↗

Frequency domain state-space system identification

An algorithm for identifying state-space models from frequency response data of linear systems is presented. A matrix-fraction description of the transfer function is employed to curve-fit the frequency response data, using the least-squares method. The parameters of the matrix-fraction representation are then used to construct the Markov parameters of the system. Finally, state-space models are obtained through the Eigensystem Realization Algorithm using Markov parameters. The main advantage of this approach is that the curve-fitting and the Markov parameter construction are linear problems which avoid the difficulties of nonlinear optimization of other approaches. Another advantage is that it avoids windowing distortions associated with other frequency domain methods.

Chen, Chung-Wen↗

Parameter estimation of the vibrational model for the SCOLE experimental facility

The objective of this study is to experimentally determine an empirical model of the vibrational dynamics of the Spacecraft COntrol Laboratory Experiment (SCOLE) facility. The first two flexible modes of this test article are identified using a linear least-square identification procedure and the data utilized for this procedure are obtained by exciting the structure from a quiescent state with torque wheels. The time history data of rate gyro sensors and accelerometers due to excitation and after excitation in terms of free-decay are used in the parameter estimation of the vibrational model. The free-decay portion of the data is analyzed using the Discrete Fourier transform to determine the optimal model order to use in modelling the response. Linear least-square analysis is then used to select the parameters that best fit the output of an Autoregressive (AR) model to the data. The control effectiveness of the torque wheels is then determined using the excitation portion of the test data, again using linear least squares.

Crotts, B. D.↗

Empirical Model of Puncture Energy for Metals

The purpose of this work is to fit a previously developed empirical equation for puncture energy to simulation data. The conservative puncture energy equation could be used to expedite the process of performing calculations in the development of safety measures, avoiding the need to create complex finite element models for specific puncture scenarios. A total of 108 simulations are developed by varying coupon thickness, coupon material, probe shape, and probe diameter. The simulations are comprised of a low-velocity probe puncturing the coupons, from which the probe kinetic energy change is calculated. The empirical equation is fit to the dimensions, material properties, and energy results using a non-linear least-squares regression method within Python, which determines the two constant parameters for each fit. More statistically significant fit results are achieved by separating the data by probe shape and coupon material.

42 ENGINEERING↗