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 415 records · Page 23

Temperature dependence of Lorentz air-broadening and pressure-shift coefficients of (12)CH4 lines in the 2.3-micron spectral region

High-resolution (0.01/cm) absorption spectra of lean mixtures of CH4 in dry air were recorded with the McMath-Pierce Fourier transform spectrometer (FTS) of the National Solar Observatory on Kitt Peak at various temperatures between 24 and -61 C. The spectra have been analyzed to determine the values at room temperature of pressure-broadened widths and pressure-induced shifts of more than 740 transitions. The temperature dependence of air-broadened widths and pressure-induced shifts was deduced for approx. 370 transitions in the nu(sub 1) + nu(sub 4), nu(sub 3) + nu(sub 4), and nu(sub 2) + nu(sub 3) bands of (12)CH4 located between 4118 and 4615/cm. These results were obtained by analyzing a total of 29 spectra simultaneously using a multi-spectral non-linear least-squares fitting technique. This new technique allowed the determination of correlated spectral line parameters (e.g. intensity and broadening coefficient) better than the procedure of averaging values obtained by fitting the spectra individually. This method also provided a direct determination of the uncertainties in the retrieved parameters due to random errors. For each band analysed in this study the dependence of the various spectral line parameters upon the tetrahedral symmetry species and the rotational quantum numbers of the transitions is also presented.

Devi, V. Malathy↗

Equatorial spacecraft-plasma interaction phenomenon observed with DE 2

Data from the retarding potential analyzer and ion drift meter on Dynamics Explorer 2 indicate that an unusual spacecraft-plasma interaction phenomemon occurs at times when the spacecraft velocity vector becomes nearly aligned with the local geomagnetic field lines. The primary signature of the interaction is a transient increase in the ion collection currents obtained with these instruments, up to approximately 15% in magnitude and typically a few tens of seconds duration. This signature is indicative of an increase in the net ion and electron currents to the satellite and is accompanied by a small positive increase in the spacecraft potential relative to the plasma. We present here a case study covering six of the strongest such events observed with DE 2 and discuss a possible physical mechanism. We suggest in particular that what might be called a collisional electron 'snowplow' effect may be occurring, and we derive a simple numerical model based on this scenario. Least squares fitting is employed to test the model and to derive new estimates of the ambient ion concentration at times when the measurements are being perturbed by the interaction.

Cragin, B. L.↗

Magnetic flux rope versus the spheromak as models for interplanetary magnetic clouds

Magnetic clouds form a subset of interplanetary ejecta with well-defined magnetic and thermodynamic properties. Observationally, it is well established that magnetic clouds expand as they propagate antisunward. The aim of this paper is to compare and contrast two models which have been proposed for the global magnetic field line topology of magnetic clouds: a magnetic flux tube geometry, on the one hand, and a spheromak geometry (including possible higher multiples), on the other. Traditionally, the magnetic structure of magnetic clouds has been modeled by force-free configurations. In a first step, we therefore analyze the ability of static force-free models to account for the asymmetries observed in the magnetic field profiles of magnetic clouds. For a cylindrical flux tube the magnetic field remains symmetric about closest approach to the magnetic axis on all spacecraft orbits intersecting it, whereas in a spheromak geometry one can have asymmetries in the magnetic field signatures along some spacecraft trajectories. The duration of typical magnetic cloud encounters at 1 AU (1 to 2 days) is comparable to their travel time from the Sun to 1 AU and thus magnetic clouds should be treated as strongly nonstationary objects. In a second step, therefore, we abandon the static approach and model magnetic clouds as self-similarly evolving MHD configurations. In our theory, the interaction of the expanding magnetic cloud with the ambient plasma is taken into account by a drag force proportional to the density and the velocity of expansion. Solving rigorously the full set of MHD equations, we demonstrate that the asymmetry in the magnetic signature may arise solely as a result of expansion. Using asymptotic solutions of the MHD equations, we least squares fit both theoretical models to interplanetary data. We find that while the central part of the magnetic cloud is adequately described by both models, the 'edges' of the cloud data are modeled better by the magnetic flux tube. Further comparisons of the two models necessarily involve thermodynamic properties, since real magnetic configurations are never exactly force-free and gas pressure plays an essential role. We consider a polytropic gas. Our theoretical analysis shows that the self-similar expansion of a magnetic flux tube requires the polytropic index gamma to be less than unity. For the spheromak, however, self-similar, radially expanding solutions are known only for gamma equal to 4/3. This difference, therefore, yields a good way of distinguishing between the two geometries. It has been shown recently that the polytropic relationship is applicable to magnetic clouds and that the corresponding polytropic index is approximately 0.5. This observational result is consistent with the self-similar model of the magnetic flux rope but is in conflict with the self-similar spheromak model.

Farrugia, C. J.↗

Thermodynamic Data For 50 Chemical Elements

Report presents data on thermo-dynamic properties of 50 chemical elements, isotope (deuterium) of one of elements, and electron gas, all in standard reference states. Data tabulated as functions of temperature; also given in form of least-squares-fit coefficients of two functional forms for heat capacities at constant pressure in standard states, with constants of integration for enthalpy and entropy.

Mcbride, Bonnie J.↗

Calibration of the advanced microwave sounding unit-A for NOAA-K

The thermal-vacuum chamber calibration data from the Advanced Microwave Sounding Unit-A (AMSU-A) for NOAA-K, which will be launched in 1996, were analyzed to evaluate the instrument performance, including calibration accuracy, nonlinearity, and temperature sensitivity. The AMSU-A on NOAA-K consists of AMSU-A2 Protoflight Model and AMSU-A1 Flight Model 1. The results show that both models meet the instrument specifications, except the AMSU-A1 antenna beamwidths, which exceed the requirement of 3.3 +/- 10%. We also studied the instrument's radiometric characterizations which will be incorporated into the operational calibration algorithm for processing the in-orbit AMSU-A data from space. Particularly, the nonlinearity parameters which will be used for correcting the nonlinear contributions from an imperfect square-law detector were determined from this data analysis. It was found that the calibration accuracies (differences between the measured scene radiances and those calculated from a linear two-point calibration formula) are polarization-dependent. Channels with vertical polarizations show little cold biases at the lowest scene target temperature 84K, while those with horizontal polarizations all have appreciable cold biases, which can be up to 0.6K. It is unknown where these polarization-dependent cold biases originate, but it is suspected that some chamber contamination of hot radiances leaked into the cold scene target area. Further investigation in this matter is required. The existence and magnitude of nonlinearity in each channel were established and a quadratic formula for modeling these nonlinear contributions was developed. The model was characterized by a single parameter u, values of which were obtained for each channel via least-squares fit to the data. Using the best-fit u values, we performed a series of simulations of the quadratic corrections which would be expected from the space data after the launch of AMSU-A on NOAA-K. In these simulations, the cosmic background radiance corresponding to a cold space temperature 2.73K was adopted as one of the two reference points of calibration. The largest simulated nonlinear correction is about 0.3K, which occurs at channel 15 when the instrument temperature is at 38.09 deg C. Others are less than 0.2K in the remaining channels. Possible improvement for future instrument calibration is also discussed.

Mo, Tsan↗

REQUEST: A Recursive QUEST Algorithm for Sequential Attitude Determination

In order to find the attitude of a spacecraft with respect to a reference coordinate system, vector measurements are taken. The vectors are pairs of measurements of the same generalized vector, taken in the spacecraft body coordinates, as well as in the reference coordinate system. We are interested in finding the best estimate of the transformation between these coordinate system.s The algorithm called QUEST yields that estimate where attitude is expressed by a quarternion. Quest is an efficient algorithm which provides a least squares fit of the quaternion of rotation to the vector measurements. Quest however, is a single time point (single frame) batch algorithm, thus measurements that were taken at previous time points are discarded. The algorithm presented in this work provides a recursive routine which considers all past measurements. The algorithm is based on on the fact that the, so called, K matrix, one of whose eigenvectors is the sought quaternion, is linerly related to the measured pairs, and on the ability to propagate K. The extraction of the appropriate eigenvector is done according to the classical QUEST algorithm. This stage, however, can be eliminated, and the computation simplified, if a standard eigenvalue-eigenvector solver algorithm is used. The development of the recursive algorithm is presented and illustrated via a numerical example.

Bar-Itzhack, Itzhack Y.↗

Kalman Filtering USNO's GPS Observations for Improved Time Transfer Predictions

The Global Positioning System (GPS) Master Control Station (MCS) performs the Coordinated Universal Time (UTC) time transfer mission by uploading and broadcasting predictions of the GPS-UTC offset in subframe 4 of the GS navigation message. These predictions are based on only two successive daily data points obtained from the US Naval Observatory (USNO). USNO produces these daily smoothed data points by performing a least-squares fit on roughly 38 hours worth of data from roughly 160 successive 13-minute tracks of GPS satellites. Though sufficient for helping to maintain a time transfer error specification of 28 ns (1 Sigma), the MCS's prediction algorithm does not make the best use of the available data from from USNO, and produces data that can degrade quickly over extended prediction spans. This paper investigates how, by applying Kalman filtering to the same available tracking data, the MCS could improve its estimate of GPS-UTC, and in particular, the GPS-UTC A(sub 1) term. By refining the A(sub 1) (frequency) estimate for GPS-UTC predictions, error in GPS time transfer could drop significantly. Additional, the risk of future spikes in GPS's time transfer error could similarly be minimized, by employing robust Kalman filtering for GPS-UTC predictions.

Hutsell, Steven T.↗

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗

Polarization Signals of Common Spacecraft Materials

This is the final report documenting the results of the polarization testing of near-planar objects with various reflectance properties. The purpose of this investigation was to determine the portion of the reflected signal which is polarized for materials commonly used in space applications. Tests were conducted on several samples, with surface characteristics ranging from highly reflective to relatively dark. The measurements were obtained by suspending the test object in a beam of collimated light. The amount of light falling on the sample was controlled by a circular aperture placed in the light field. The polarized reflectance at various phase angles was then measured. A nonlinear least squares fitting program was used for analysis. For the specular test objects, the reflected signals were measured in one degree increments near the specular point. Otherwise, measurements were taken every five degrees in phase angle. Generally, the more diffuse surfaces had lower polarized reflectances than their more specular counterparts. The reflected signals for the more diffuse surfaces were spread over a larger phase angle range, while the signals from the more specular samples were reflected almost entirely within five degrees of angular deviation from the specular point. The method used to test all the surfaces is presented. The results of this study will be used to support the NASA Orbital Debris Optical Signature Tests. These tests are intended to help better understand the reflectance properties of materials often used in space applications. This data will then be used to improve the capabilities for identification and tracking of space debris.

Gravseth, Ian↗

Equatorial Spacecraft-Plasma Interaction Phenomenon Observed With DE 2

Data from the retarding potential analyzer and ion drift meter on Dynamics Explorer 2 indicate that an unusual spacecraft-plasma interaction phenomenon occurs at times when the spacecraft velocity vector becomes nearly aligned with the local geomagnetic field lines. The primary signature of the interaction is a transient increase in the ion collection currents obtained with these instruments, up to approximately 15% in magnitude and typically a few tens of seconds duration. This signature is indicative of an increase in the net ion and electron currents to the satellite and is accompanied by a small positive increase in the spacecraft potential relative to the plasma. We present here a case study covering six of the strongest such events observed with DE 2 and discuss a possible physical mechanism. We suggest in particular that what might be called a collisional electron snowplow effect may be occurring, and we derive a simple numerical model based on this scenario. Least squares fitting is employed to test the model and to derive new estimates of the ambient ion concentration at times when the measurements are being perturbed by the interaction.

Cragin, B. L.↗

Land Surface Temperature Measurements form EOS MODIS Data

We have developed a physics-based land-surface temperature (LST) algorithm for simultaneously retrieving surface band-averaged emissivities and temperatures from day/night pairs of MODIS (Moderate Resolution Imaging Spectroradiometer) data in seven thermal infrared bands. The set of 14 nonlinear equations in the algorithm is solved with the statistical regression method and the least-squares fit method. This new LST algorithm was tested with simulated MODIS data for 80 sets of band-averaged emissivities calculated from published spectral data of terrestrial materials in wide ranges of atmospheric and surface temperature conditions. Comprehensive sensitivity and error analysis has been made to evaluate the performance of the new LST algorithm and its dependence on variations in surface emissivity and temperature, upon atmospheric conditions, as well as the noise-equivalent temperature difference (NE(Delta)T) and calibration accuracy specifications of the MODIS instrument. In cases with a systematic calibration error of 0.5%, the standard deviations of errors in retrieved surface daytime and nighttime temperatures fall between 0.4-0.5 K over a wide range of surface temperatures for mid-latitude summer conditions. The standard deviations of errors in retrieved emissivities in bands 31 and 32 (in the 10-12.5 micrometer IR spectral window region) are 0.009, and the maximum error in retrieved LST values falls between 2-3 K. Several issues related to the day/night LST algorithm (uncertainties in the day/night registration and in surface emissivity changes caused by dew occurrence, and the cloud cover) have been investigated. The LST algorithms have been validated with MODIS Airborne Simulator (MAS) dada and ground-based measurement data in two field campaigns conducted in Railroad Valley playa, NV in 1995 and 1996. The MODIS LST version 1 software has been delivered.

Wan, Zhengming↗

Measurement of Charge Transfer Rate Coefficient Between Ground-State N(2+) Ion and He at Electron-Volt Energies

The charge transfer rate coefficient for the reaction N(2+)(2p(sup 2)P(sup 0)) + He yields products is measured by recording the time dependence of the N(2+) ions stored in an ion trap. A cylindrical radio-frequency ion trap was used to store N(2+) ions produced by laser ablation of a solid titanium nitride target. The decay of the ion signals was analyzed by single exponential least-squares fits to the data. The measured rate coefficient is 8.67(0.76) x 10(exp -11)sq cm/s. The N(2+) ions were at a mean energy of 2.7 eV while He gas was at room temperature, corresponding to an equivalent temperature of 3.9 x 10(exp 3) K. The measured value is in good agreement with a recent calculation.

Fang, Z.↗

Retrieval of Surface Snow Grain Size and Melt Water from AVIRIS Spectra

The Earth's energy balance and hydrology are affected by the distribution and characteristics of snow cover on the surface. Snow grain size and snow melt influence surface albedo and hydrology. A model of snow reflectance that depends on both grain size and surface melt water was developed to derive these parameters from remote spectral measurements. This reflectance model is based on a discrete ordinate radiative transfer approach that uses Mie calculations of snow optical properties, which are based on the complex refractive index of ice and water. This snow model was linked to an atmospheric radiative transfer code and a nonlinear least squares fitting algorithm. The resulting combined algorithm was applied to an AVIRIS snow data set acquired over Mammoth Mountain, California. Maps of grain size and surface snow melt were generated that are consistent with the expected ranges and distributions for conditions at the site.

Green Robert O.↗

Collision Integrals for a Modified Stockmayer Potential

Collision integrals were calculated for the modified Stockmayer potential E(r) = 4 epsilon [(sigma/r)(exp 12) - (sigma/r)(exp 6) - delta(sigma/r)(exp 3)], which may be applied to polar molecules. It was assumed that the colliding molecules maintain their same relative orientation during the encounter. Calculations of the integrals were made for a large reduced temperature range and for a range of delta from 0 to 10. The results agree with other work on non-polar interactions (delta = 0). However for polar interactions the only previously published calculations have been found to be in error and do not agree with this work. Assuming that the molecules interact as alined dipoles of maximum attraction, values for sigma, epsilon, and delta were determined for various polar molecules by a least squares fit of experimental viscosity data. Satisfactory results were obtained for slightly polar molecules, but not for more highly polar molecules such as NH3 or H2O. Therefore, it appears that the assumed model of molecules interacting at all times as alined dipoles of maximum attraction is not satisfactory for estimating trans- port properties of polar molecules.

Itean, Engene C.↗

Evidence for Excitation of Polar Motion by Fortnightly Ocean Tides

The second-degree zonal tide raising potential, which is responsible for tidal changes in the Earth's rotation rate and length-of-day, is symmetric about the polar axis and hence can excite the Earth's polar motion only through its action upon nonaxisymmetric features of the Earth such as the oceans. Ocean tidal excitation of polar motion in the diurnal and semidiurnal tidal bands has been previously detected and examined. Here, the detection of ocean tidal excitation of polar motion in the long-period tidal band, specifically at the Mf' (13.63-day) and Mf (13.66-day) tidal frequencies, is reported. Spectra of the SPACE94 polar motion excitation series exhibit peaks at the prograde and retrograde fortnightly tidal periods. After removing effects of atmospheric wind and pressure changes, an empirical model for the effect of the fortnightly ocean tides upon polar motion excitation is obtained by least-squares fitting periodic terms at the Mf and Mf' tidal frequencies to the residual polar motion excitation series. The resulting empirical model is then compared with the predictions of two hydrodynamic ocean tide models.

Gross, Richard S.↗

Self-Broadening and Self-Shift Coefficients in the Fundamental Band of 12C 16O

High quality and precise measurements of self-broadened and self-shift coefficients in the fundamental band of C-12O-16 were made using spectra recorded at room temperature with the high-resolution (0.0027 cm(exp -1)) McMath-Pierce Fourier transform spectrometer located at the National Solar Observatory on Kitt Peak, Arizona. The spectral region under investigation (2008-2247 cm(exp -1)) contains the P(31) to R(31) transitions. The data were obtained using a high-purity natural isotopic sample ofcarbon monoxide and two absorption cells with pathlengths of 4.08 and 9.98 cm, respectively. Various pressures of CO were used, ranging between 0.25 and 201.2 Torr. The results were obtained by analyzing five spectra simultaneously, using a multispectrum nonlinear least-squares fitting technique. The self-broadened coefficients ranged from 0.0426(2) cm(exp -1) atm(exp -1) at 296 K to 0.0924(2) cm(exp -1) atm(exp -1) at 296 K, while the pressure-induced shift coefficients varied between -0.0042(3) cm(exp -1) atm(exp -1) at 296 K and +0.0005(l) cm(exp -1) atm(exp -1) at 296 K. The value in parentheses is the estimated uncertainty in units of the last digit. The self-broadened coefficients of lines with same values of m in the P and R branches agree close to within experimental uncertainties while the self-shift coefficients showed considerable variation within and between the two branches. The mean value of the ratios of P branch to R branch self-broadened coefficients was found to be 1.01 with a standard deviation of + or - 0.01. Comparisons of the results with other published data were made.

Devi, Malathy V.↗

Time-Frequency Analysis of the Dispersion of Lamb Modes

Accurate knowledge of the velocity dispersion of Lamb modes is important for ultrasonic nondestructive evaluation methods used in detecting and locating flaws in thin plates and in determining their elastic stiffness coefficients. Lamb mode dispersion is also important in the acoustic emission technique for accurately triangulating the location of emissions in thin plates. In this research, the ability to characterize Lamb mode dispersion through a time-frequency analysis (the pseudo-Wigner-Ville distribution) was demonstrated. A major advantage of time-frequency methods is the ability to analyze acoustic signals containing multiple propagation modes, which overlap and superimpose in the time domain signal. By combining time-frequency analysis with a broadband acoustic excitation source, the dispersion of multiple Lamb modes over a wide frequency range can be determined from as little as a single measurement. In addition, the technique provides a direct measurement of the group velocity dispersion. The technique was first demonstrated in the analysis of a simulated waveform in an aluminum plate in which the Lamb mode dispersion was well known. Portions of the dispersion curves of the AO, A I , So, and S2 Lamb modes were obtained from this one waveform. The technique was also applied for the analysis of experimental waveforms from a unidirectional graphite/epoxy composite plate. Measurements were made both along and perpendicular to the fiber direction. In this case, the signals contained only the lowest order symmetric and antisymmetric modes. A least squares fit of the results from several source to detector distances was used. Theoretical dispersion curves were calculated and are shown to be in good agreement with experimental results.

Prosser, W. H.↗

Gaseous abundances and methane supersaturation in Titan's troposphere

Various properties of Titan's troposphere are inferred from an analysis of Voyager I infrared spectrometer (IRIS) data between 200 and 600/ cm. Two homogeneous spectral averages acquired at widely separated emission angles are chosen for the analysis. Both data sets are associated with northern low latitudes very close to that of the radio science ingress occultation point. Solutions require simultaneous nonlinear least-squares fits to the two IRIS data sets, coupled with iteration of the radio occultation refractivity data. Values and associated 1-sigma uncertainties of several parameters are inferred from our analysis. These include mole fractions for molecular hydrogen (approx. 0.0011), argon (small), and methane near the surface ( approx. 0.057). Solutions are also obtained for the hydrogen parafraction (close to equilibrium, with considerable uncertainty), air temperature near the surface ( approx. 93 K), surface surface temperature discontinuity (approx. 1 K), and maximum degree of methane supersaturation in the upper troposphere (approx. 1.5). Actual values for the above-mentioned parameters depend on the amount of ethane cloud near the tropopause. There is no evidence for methane clouds in the upper troposphere, nor is their presence compatible with large degrees of supersaturation. A wave number dependence for the stratospheric haze opacity is inferred similar to that found for a polymeric residue created in laboratory discharge experiments. This haze appears to be uniformly distributed with latitude between altitudes of 40 and 160 km, provided those nighttime data at southern high latitudes that are discounted.

Samuelson, Robert E.↗