Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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 1,315 records · Page 73

Energization and transport of ions of ionospheric origin in the terrestrial magnetosphere

The work of this grant has been predominantly focused on ion outflows from two data sets: Prognoz 7 and Dynamics Explorer. The Prognoz analysis studied ion densities, temperatures, and flow velocities in the magnetotail. The work performed under this contract consisted of developing a program to load the raw data, computing the background subtraction of a strong sun pulse, and using the net count to calculate the low order moments of the distribution function. The study confirms the results of ISEE with regard to the supply of plasma from the cusp as a major source of plasmasheet plasma and goes beyond this to discuss the use of ion velocities as a way to examine the motions of the magnetotail. The abstract of the work to be reported is included as an appendix. The work on the DE/Retarding Ion Mass Spectrometer is separated into two categories: (1) classification of low-energy ion flows from high-latitudes, and (2) studies of the polar wind. Major publications resulting from this work are also included as an appendix to this report. The polar wind is in a category by itself as a result of the thermal escape of hydrogen and helium because of ambipolar diffusion through the heavier, oxygen-dominated topside ionosphere. The analysis of the polar wind reports the flux variability as a function of season, magnetic activity, etc. Much effort has been expended under this grant to complete a follow on study of the thermal structure of the polar wind. Extensive display tools and analysis software have been developed and used in an attempt to carry out this thermal analysis. The present work uses a constrained fit scheme that combines the ion densities and flow velocities derived from Chandler et al. and a spacecraft potential derived from an empirical relation to the total ion density to determine the remaining fit parameter, the ion temperature, via a least squares fit to the RIMS data.

Waite, J. Hunter JR.↗

Automatic Data Filter Customization Using a Genetic Algorithm

This work predicts whether a retrieval algorithm will usefully determine CO2 concentration from an input spectrum of GOSAT (Greenhouse Gases Observing Satellite). This was done to eliminate needless runtime on atmospheric soundings that would never yield useful results. A space of 50 dimensions was examined for predictive power on the final CO2 results. Retrieval algorithms are frequently expensive to run, and wasted effort defeats requirements and expends needless resources. This algorithm could be used to help predict and filter unneeded runs in any computationally expensive regime. Traditional methods such as the Fischer discriminant analysis and decision trees can attempt to predict whether a sounding will be properly processed. However, this work sought to detect a subsection of the dimensional space that can be simply filtered out to eliminate unwanted runs. LDAs (linear discriminant analyses) and other systems examine the entire data and judge a "best fit," giving equal weight to complex and problematic regions as well as simple, clear-cut regions. In this implementation, a genetic space of "left" and "right" thresholds outside of which all data are rejected was defined. These left/right pairs are created for each of the 50 input dimensions. A genetic algorithm then runs through countless potential filter settings using a JPL computer cluster, optimizing the tossed-out data s yield (proper vs. improper run removal) and number of points tossed. This solution is robust to an arbitrary decision boundary within the data and avoids the global optimization problem of whole-dataset fitting using LDA or decision trees. It filters out runs that would not have produced useful CO2 values to save needless computation. This would be an algorithmic preprocessing improvement to any computationally expensive system.

Mandrake, Lukas↗

A comparison of time domain boundary conditions for acoustic waves in wave guides

Researchers consider several types of boundary conditions in the context of time domain models for acoustic waves. Experiments with four different duct terminations (hard wall, free radiation, foam, and wedge) were carried out in a wave duct from which reflection coefficients over a wide frequency range were measured. These reflection coefficients were used to estimate parameters in the time domain boundary conditions. A comparison of the relative merits of the models in describing the data is presented. Boundary conditions which yield a good fit of the model to the experimental data were found for all duct terminations except the wedge.

Banks, H. T.↗

Absence of diurnal variation of C-reactive protein concentrations in healthy human subjects

BACKGROUND: The concentration of C-reactive protein (CRP) in otherwise healthy subjects has been shown to predict future risk of myocardial infarction and stroke. CRP is synthesized by the liver in response to interleukin-6, the serum concentration of which is subject to diurnal variation. METHODS: To examine the existence of a time-of-day effect for baseline CRP values, we determined CRP concentrations in hourly blood samples drawn from healthy subjects (10 males, 3 females; age range, 21-35 years) during a baseline day in a controlled environment (8 h of nighttime sleep). RESULTS: Overall CRP concentrations were low, with only three subjects having CRP concentrations >2 mg/L. Comparison of raw data showed stability of CRP concentrations throughout the 24 h studied. When compared with cutoff values of CRP quintile derived from population-based studies, misclassification of greater than one quintile did not occur as a result of diurnal variation in any of the subjects studied. Nonparametric ANOVA comparing different time points showed no significant differences for both raw and z-transformed data. Analysis for rhythmic diurnal variation using a method fitting a cosine curve to the group data was negative. CONCLUSIONS: Our data show that baseline CRP concentrations are not subject to time-of-day variation and thus help to explain why CRP concentrations are a better predictor of vascular risk than interleukin-6. Determination of CRP for cardiovascular risk prediction may be performed without concern for diurnal variation.

Non-NASA Center↗

Instrumental Biases in Ionospheric Measurement Derived from GPS Data

Line-of-sight ionosphere measurements derived from differencing dual frequency GPS range data are corrupted by instrumental biases in both receiver and GPS satellite transmitters due to hardware delays in the L1 and L2 signal paths. Using a global ionospheric shell model to fit GPS based ionospheric delay data from a world wide network of 30 to 40 receivers, we can estimate, with a single fit, satellite biases for the entire GPS constellation and receiver biases for all the uncalibrated receivers.

GPS↗

Non-linear Least Square Fitting Technique for the Determination of Field Line Resonance Frequency in Ground Magnetometer Data: Application to Remote Sensing of Plasmaspheric Mass Density

The accurate determination of the Field Line Resonance (FLR) frequency of a resonating geomagnetic field line is necessary to remotely monitor the plasmaspheric mass density during geomagnetic storms and quiet times alike. Under certain assumptions the plasmaspheric mass density at the equator is inversely proportional to the square of the FLR frequency. The most common techniques to determine the FLR frequency from ground magnetometer measurements are the amplitude ratio and phase difference techniques, both based on geomagnetic field observations at two latitudinally separated ground stations along the same magnetic meridian. Previously developed automated techniques have used statistical methods to pinpoint the FLR frequency using the amplitude ratio and phase difference calculations. We now introduce a physics-based automated technique, using non-linear least square fitting of the ground magnetometer data to the analytical resonant wave equations, that reproduces the wave characteristics on the ground, and from those determine the FLR frequency. One of the advantages of the new technique is the estimation of physics-based errors of the FLR frequency, and as a result of the equatorial plasmaspheric mass density. We present analytical results of the new technique, and test it using data from the Inner-Magnetospheric Array for Geospace Science (iMAGS) ground magnetometer chain along the coast of Chile and the east coast of the United States. We compare the results with the results of previously published statistical automated techniques.

A. Boudouridis↗

Fatigue behavior of axial and pressure cycled butt and girth welds containing defects

This paper presents the results of a study directed at developing a data base for establishing fitness-for-purpose defect acceptance criteria for welds with defects. The study focused on A106 Grade B steel pipe. Data are presented for flat plate, wall segment, and vessel specimens and actual pipe sections containing either artificial or natural planar or volumetric defects. Defect acceptance criteria developed from the test data are discussed.

Leis, Brian N.↗

Aeroelastic Model Structure Computation for Envelope Expansion

Structure detection is a procedure for selecting a subset of candidate terms, from a full model description, that best describes the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modelling may be of critical importance in the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion which may save significant development time and costs. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of nonlinear aeroelastic systems. The LASSO minimises the residual sum of squares by the addition of an l(sub 1) penalty term on the parameter vector of the traditional 2 minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudolinear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. Applicability of this technique for model structure computation for the F/A-18 Active Aeroelastic Wing using flight test data is shown for several flight conditions (Mach numbers) by identifying a parsimonious system description with a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

Aeroelastic Model Structure Computation for Envelope Expansion

Structure detection is a procedure for selecting a subset of candidate terms, from a full model description, that best describes the observed output. This is a necessary procedure to compute an efficient system description which may afford greater insight into the functionality of the system or a simpler controller design. Structure computation as a tool for black-box modeling may be of critical importance in the development of robust, parsimonious models for the flight-test community. Moreover, this approach may lead to efficient strategies for rapid envelope expansion that may save significant development time and costs. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of non-linear aeroelastic systems. The LASSO minimises the residual sum of squares with the addition of an l(Sub 1) penalty term on the parameter vector of the traditional l(sub 2) minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudo-linear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. Applicability of this technique for model structure computation for the F/A-18 (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) Active Aeroelastic Wing project using flight test data is shown for several flight conditions (Mach numbers) by identifying a parsimonious system description with a high percent fit for cross-validated data.

Kukreja, Sunil L.↗

A Data Exploration Tool for Large Sets of Spectra

We present an exploration tool for very large spectrum data sets such as the SDSS (Sloan Digital Sky Survey), LAMOST (Large Sky Area Multi-Object Fiber Spectroscopic Telescope), and 4MOST (4-meter Multi-Object Spectroscopic Telescope) data sets. The tool works in two stages: the first uses batch processing and the second runs interactively. The latter employs the NASA hyperwall, a configuration of 128 workstation displays (8 by 16 array) controlled by a parallelized software suite running on NASA's Pleiades supercomputer. The stellar subset of the Sloan Digital Sky Survey, DR10, was chosen to show how the our tool may be used. In stage one, SDSS files for 569,740 stars are processed through our data pipeline. The pipeline fits each spectrum using an iterative continuum algorithm, distinguishing emission from absorption and handling molecular absorption bands correctly. It then measures 1659 discrete atomic and molecular spectral features that were carefully preselected based on their likelihood of being visible at some spectral type. The depths relative to the local continuum at each feature wavelength are determined for each spectrum: these depths, the local S/N (signal to noise ratio) level, and DR10-supplied variables such as magnitudes, colors, positions, and radial velocities are the basic measured quantities used on the hyperwall. In stage two, each hyperwall panel is used to display a 2-D scatter plot showing the depth of feature A vs the depth of feature B for all of the stars. A and B change from panel to panel. The relationships between the various (A,B) strengths and any distinctive clustering are immediately apparent when examining and inter-comparing the different panels on the hyperwall. The interactive software allows the user to select the stars in any interesting region of any 2-D plot on the hyperwall, immediately rendering the same stars on all the other 2-D plots in a unique color. The process may be repeated multiple times, each selection displaying a distinctive color on all the plots. At any time, the spectra of the selected stars may be examined in detail on a connected workstation display. We illustrate how our approach allows us to quickly isolate and examine such interesting stellar subsets as EMP (Extremely Metal‐Poor) stars, CV (Cataclymic Variable) stars and C (Carbon)-rich stars.

Data Exploration↗

Compound equation developed for postnatal growth of birds and mammals

Compound growth equation was developed in which the rate of this linear growth process is regarded as proportional to the mass already attained at any instant by an underlying Gompertz process. This compound growth model was fitted to the growth data of a variety of birds and mammals of both sexes.

Laird, A. K.↗

Scattering of 42-MeV alpha particles from Cu-65

The extended particle-core coupling model was used to predict the properties of low-lying levels of Cu-65. A 42-MeV alpha particle cyclotron beam was used for the experiment. The experiment included magnetic analysis of the incident beam and particle detection by lithium-drifted silicon semiconductors. Angular distributions were measured for 10 to 50 degrees in the center of mass system. Data was reduced by fitting the peaks with a skewed Gaussian function using a least squares computer program with a linear background search. The energy calibration of each system was done by pulsar, and the excitation energies are accurate to + or - 25 keV. The simple weak coupling model cannot account for the experimentally observed quantities of the low-lying levels of Cu-65. The extended particle-core calculation showed that the coupling is not weak and that considerable configuration mixing of the low-lying states results.

Stewart, W. M.↗

Gadolinium-hydrogen ion exchange of zirconium phosphate

The Gd(+3)/H(+) ion exchange on a commercial zirconium phosphate ion exchanger was investigated in chloride, sulfate, and phosphate solutions of Gd(+3) at gadolinium concentrations of 0.001 to 1 millimole per cc and in the pH range of 0 to 3.5. Relatively low Gd(+3) capacities, in the range of 0.01 to 0.1 millimole per g of ion exchanger were found at room temperature. A significant difference in Gd(+3) sorption was observed, depending on whether the ion exchanger was converted from initial conditions of greater or lesser Gd(+3) sorption than the specific final conditions. Correlations were found between decrease in Gd(+3) capacity and loss of exchanger phosphate groups due to hydrolysis during washing and between increase in capacity and treatment with H3PO4. Fitting of the experimental data to ideal ion exchange equilibrium expressions indicated that each Gd(+3) ion is sorbed on only one site of the ion exchanger. The selectivity quotient was determined to be 2.5 + or - 0.4 at room temperature on gadolinium desorption in chloride solutions.

Liu, D. C.↗

Refractive index and birefringence of 2H silicon carbide.

Measurement of the refractive indices of 2H SiC over the wavelength range from 435.8 to 650.9 nm by the method of minimum deviation. A curve fit of the experimental data to the Cauchy dispersion equation yielded, for the ordinary index, n sub zero = 2.5513 + 25,850/lambda squared + 8.928 x 10 to the 8th power/lambda to the 4th power and, for the extraordinary index, n sub e = 2.6161 + 28,230/lambda squared + 11.490 x 10 to the 8th power/lambda to the 4th power when lambda is expressed in nm. The estimated error (standard deviation) in these values is plus or minus 0.0006 for n sub zero and plus or minus 0.0009 for n sub e. The birefringence calculated from these expressions is about 20% less than previously published values.

Powell, J. A.↗

The F and K components of the solar corona.

Numerical calculations of the F component are conducted for a large range of parameters, particle size distribution, minimum particle size, minimum approach to the sun, albedo, and spatial distribution. Various parametrized functional forms for the electron spatial distribution are assumed, and the parametrical elongation dependence of the K component is calculated. New and old data are analyzed by fitting the theoretical parametrized elongation dependence to observations of the total intensity of the solar corona. The orbital behavior of the dust is estimated together with the temperature and thermal emission.

Calbert, R.↗