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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Extrapolation of the dna fragment-size distribution after high-dose irradiation to predict effects at low doses

The patterns of DSBs induced in the genome are different for sparsely and densely ionizing radiations: In the former case, the patterns are well described by a random-breakage model; in the latter, a more sophisticated tool is needed. We used a Monte Carlo algorithm with a random-walk geometry of chromatin, and a track structure defined by the radial distribution of energy deposition from an incident ion, to fit the PFGE data for fragment-size distribution after high-dose irradiation. These fits determined the unknown parameters of the model, enabling the extrapolation of data for high-dose irradiation to the low doses that are relevant for NASA space radiation research. The randomly-located-clusters formalism was used to speed the simulations. It was shown that only one adjustable parameter, Q, the track efficiency parameter, was necessary to predict DNA fragment sizes for wide ranges of doses. This parameter was determined for a variety of radiations and LETs and was used to predict the DSB patterns at the HPRT locus of the human X chromosome after low-dose irradiation. It was found that high-LET radiation would be more likely than low-LET radiation to induce additional DSBs within the HPRT gene if this gene already contained one DSB.

NASA Discipline Radiation Health↗
NASA NTRS2012

Algorithm for Compressing Time-Series Data

An algorithm based on Chebyshev polynomials effects lossy compression of time-series data or other one-dimensional data streams (e.g., spectral data) that are arranged in blocks for sequential transmission. The algorithm was developed for use in transmitting data from spacecraft scientific instruments to Earth stations. In spite of its lossy nature, the algorithm preserves the information needed for scientific analysis. The algorithm is computationally simple, yet compresses data streams by factors much greater than two. The algorithm is not restricted to spacecraft or scientific uses: it is applicable to time-series data in general. The algorithm can also be applied to general multidimensional data that have been converted to time-series data, a typical example being image data acquired by raster scanning. However, unlike most prior image-data-compression algorithms, this algorithm neither depends on nor exploits the two-dimensional spatial correlations that are generally present in images. In order to understand the essence of this compression algorithm, it is necessary to understand that the net effect of this algorithm and the associated decompression algorithm is to approximate the original stream of data as a sequence of finite series of Chebyshev polynomials. For the purpose of this algorithm, a block of data or interval of time for which a Chebyshev polynomial series is fitted to the original data is denoted a fitting interval. Chebyshev approximation has two properties that make it particularly effective for compressing serial data streams with minimal loss of scientific information: The errors associated with a Chebyshev approximation are nearly uniformly distributed over the fitting interval (this is known in the art as the "equal error property"); and the maximum deviations of the fitted Chebyshev polynomial from the original data have the smallest possible values (this is known in the art as the "min-max property").

Hawkins, S. Edward, III↗
NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗
NASA NTRS1992

Application of fiber bridging models to fatigue crack growth in unidirectional titanium matrix composites

Several fiber bridging models were reviewed and applied to study the matrix fatigue crack growth behavior in center notched (0)(sub 8) SCS-6/Ti-15-3 and (0)(sub 4) SCS-6/Ti-6Al-4V laminates. Observations revealed that fatigue damage consisted primarily of matrix cracks and fiber matrix interfacial failure in the (0)(sub 8) SCS-6/Ti-15-3 laminates. Fiber-matrix interface failure included fracture of the brittle reaction zone and cracking between the two carbon rich fiber coatings. Intact fibers in the wake of the matrix cracks reduce the stress intensity factor range. Thus, an applied stress intensity factor range is inappropriate to characterize matrix crack growth behavior. Fiber bridging models were used to determine the matrix stress intensity factor range in titanium metal matrix composites. In these models, the fibers in the wake of the crack are idealized as a closure pressure. An unknown constant frictional shear stress is assumed to act along the debond or slip length of the bridging fibers. The frictional shear stress was used as a curve fitting parameter to available data (crack growth data, crack opening displacement data, and debond length data). Large variations in the frictional shear stress required to fit the experimental data indicate that the fiber bridging models in their present form lack predictive capabilities. However, these models provide an efficient and relatively simple engineering method for conducting parametric studies of the matrix growth behavior based on constituent properties.

Bakuckas, J. G., Jr.↗
NASA NTRS1991

Exocentric judgements in real environments and stereoscopic displays

Spatial direction errors during interpretation of perspective images, such as 3D map displays, may originate from misjudgment of the orientation of the viewing direction used to make the display. One source of these errors could be perceptual evidence of the display surface. Two experiments are reported in which the same judgment exocentric task was presented, but the cues to the picture surface were reduced or eliminated by presenting the task as a stereoscopic, virtual image or by a geometrically matched physical model. A theory developed to model exocentric direction errors on perspective displays has been fitted to the data from these two experiments. The parameters estimated from the fit in both experiments indicate that the subjects may be more correctly estimating the viewing direction than in ordinary perspective displays. Consequently, in some real world or stereo viewing conditions, errors in estimating the viewing direction are not likely to dominate exocentric direction errors.

Ellis, Stephen R.↗
NASA NTRS2007

Gamma-Ray Burst Afterglows as Probes of Environment and Blastwave Physics II: The Distribution of p and Structure of the Circumburst Medium

We constrain blastwave parameters and the circumburst media of a subsample of BeppoSAX Gamma-Ray Bursts. For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical and nIR afterglow data. The spectral fits have been done in count space and include the effects of metallicity, and are compared with the previously reported optical and X-ray temporal behaviour. Assuming the fireball model, we can find a mean value of p for the sample as a whole of 2.035. A statistical analysis Of the distribution demonstrates that the p values in this sample are inconsistent with a single universal value for p at the 3sigma level or greater. This approach provides us with a measured distribution of circumburst density structures rather than considering only the cases of k = 0 (homogeneous) and k = 2 (wind-like). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly wind-like. The fifth source has a value of 0 less than or equal to k less than or equal to 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗
NASA NTRS2008

Gamma-Ray Burst Afterglows as Probes of Environment and Blast Wave Physics. II. The Distribution of rho and Structure of the Circumburst Medium

We constrain blast wave parameters and the circumburst media ofa subsample of 10 BeppoSAX gamma-ray bursts (GRBs). For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical, and NIR afterglow data. The spectral fits have been done in count space and include the effects ofmetallicity, and are compared with the previously reported optical and X-ray temporal behavior. Using the blast wave model and some assumptions which include on-axis viewing and standard jet structure, constant blast wave energy, and no evolution of the microphysical parameters, we find a mean value ofp for the sample as a whole of 9.... oa -0.003.0" 2 a_ statistical analysis of the distribution demonstrates that the p-values in this sample are inconsistent with a single universal value forp at the 3 _ level or greater, which has significant implications for particle acceleration models. This approach provides us with a measured distribution ofcircumburst density structures rather than considering only the cases of k ----0 (homogeneous) and k - 2 (windlike). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly windlike. The fifth source has a value of 0 < k < 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗
NASA NTRS2011

Episodes of Cross-Polar Transport in the Arctic Troposphere During July 2008 as Seen from Models, Satellite, and Aircraft Observations

During the POLARCAT summer campaign in 2008, two episodes (2 5 July and 7 10 July 2008) occurred where low-pressure systems traveled from Siberia across the Arctic Ocean towards the North Pole. The two cyclones had extensive smoke plumes from Siberian forest fires and anthropogenic sources in East Asia embedded in their associated air masses, creating an excellent opportunity to use satellite and aircraft observations to validate the performance of atmospheric transport models in the Arctic, which is a challenging model domain due to numerical and other complications. Here we compare transport simulations of carbon monoxide (CO) from the Lagrangian transport model FLEXPART and the Eulerian chemical transport model TOMCAT with retrievals of total column CO from the IASI passive infrared sensor onboard the MetOp-A satellite. The main aspect of the comparison is how realistic horizontal and vertical structures are represented in the model simulations. Analysis of CALIPSO lidar curtains and in situ aircraft measurements provide further independent reference points to assess how reliable the model simulations are and what the main limitations are. The horizontal structure of mid-latitude pollution plumes agrees well between the IASI total column CO and the model simulations. However, finer-scale structures are too quickly diffused in the Eulerian model. Applying the IASI averaging kernels to the model data is essential for a meaningful comparison. Using aircraft data as a reference suggests that the satellite data are biased high, while TOMCAT is biased low. FLEXPART fits the aircraft data rather well, but due to added background concentrations the simulation is not independent from observations. The multi-data, multi-model approach allows separating the influences of meteorological fields, model realisation, and grid type on the plume structure. In addition to the very good agreement between simulated and observed total column CO fields, the results also highlight the difficulty to identify a data set that most realistically represents the actual pollution state of the Arctic atmosphere.

Sodemann, H.↗
NASA NTRS2015

Enhancements of Bayesian Blocks; Application to Large Light Curve Databases

Bayesian Blocks are optimal piecewise linear representations (step function fits) of light-curves. The simple algorithm implementing this idea, using dynamic programming, has been extended to include more data modes and fitness metrics, multivariate analysis, and data on the circle (Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations, Scargle, Norris, Jackson and Chiang 2013, ApJ, 764, 167), as well as new results on background subtraction and refinement of the procedure for precise timing of transient events in sparse data. Example demonstrations will include exploratory analysis of the Kepler light curve archive in a search for "star-tickling" signals from extraterrestrial civilizations. (The Cepheid Galactic Internet, Learned, Kudritzki, Pakvasa1, and Zee, 2008, arXiv: 0809.0339; Walkowicz et al., in progress).

Kepler light curve archive↗
NASA NTRS2020

Characteristics of DSD Bulk Parameters: Implication to Radar Rain Retrieval

With the use of 213,456 one-minute measured data of droplet-size distribution (DSD) of rain collected during several National Aeronautics and Space Administration (NASA)-sponsored field campaigns, the relationships between rainfall rate R, mass-weighted diameter D(m) and normalized intercept parameter N(w) of the gamma DSD are studied. It is found, based on the simulations of the gamma DSD model, that R, D(m) and N(w) are closely interrelated, and that the ratio of R to N(w) is solely a function of D(m), independent of the shape factor μ of the gamma distribution. Furthermore, the model-produced ratio agrees well with those from the DSD data. When a power-law equation is applied to fit the model data, we have: R=aN(w)D(sup b, sub m) , where a=1.588×10(exp −4) , b=4.706 . Analysis of two-parameter relationships such as R–D(m), N(w)–R and N(w)–D(m) reveals that R and D(m) are moderately correlated while N(w) and D(m) are negatively correlated. N(w) and R, however, are uncorrelated. The gamma DSD model also reveals that variation of R–D(m) relation is caused primarily by N(w). For the application of the Ku- and Ka-band dual-frequency radar for the retrieval of the DSD bulk parameters as well as the specific radar attenuations, the study is carried out to relate the dual-frequency radar reflectivity factors to the DSD and attenuation parameters.

Liang Liao↗
NASA NTRSDate not supplied

A Misperception in Reliability Growth Modelling

The Duane reliability growth model is n(t)/t = k t^-alpha (1) The reliability growth rate is alpha, the downward slope of n(t)/t versus t. It usually varies from 0.2 to 0.6. k is a constant. Crow used a 56-failure data set to illustrate reliability growth.1 A graphical Duane model fit to this data gives n(t)/t = 0.640 t^-0.283 (2) A problem in using the Duane-Crow reliability growth model is that it assumes that reliability growth continues and the failure rate decreases throughout the test period. It is more usual that reliability growth stops when the failure rated is low enough. Growth testing is often followed by testing with a low constant failure rate due to rare or uncorrectable failure modes. As more and more low constant rate acceptable failures accumulate after the period of reliability growth, the reliability growth time exponent alpha decreases toward zero. This occurs if constant rate failures are treated as occurring during the reliability growth period. It is more accurate to model a period of initial reliability growth followed by testing without repair to more accurately determine the final constant failure rate. This is done in the abcd model. n(t)/t = a t^-b + c from t = 0 to td (3) = c + d after td, where d = a td^-b (4) The term a t^-b describes the continuous reliability growth that continues out to time td and c is the constant uncorrected failure rate. The parameter d represents an additional constant failure rate due to correctable but uncorrected failure modes. After the reliability growth process is terminated, the failure rate n(t)/t = c + d.

Harry W Jones↗
NASA NTRS1982

Active region coronal evolution

Scaling relations between coronal base pressure and longitudinal photospheric magnetic field strength are tested for the case of a single active region observed for five solar rotations from Skylab. The evolution of measureable quantities, such as coronal thermal energy content, total longitudinal photospheric magnetic flux, region scale size, and peak energy density, is traced throughout the five rotations observed. The theoretically derived scaling law of Golub et al. (1980) is found to provide an acceptable fit to the data throughout the entire evolutionary history of the region from an age of about 3 days to the fully evolved state in which the mature active region merges into the general large-scale structure of the quiet corona. An alternative scaling law obtained by including the results of Galeev et al. (1981), however, is found to provide a somewhat better fit to the data. The study is seen as providing additional justification for the belief that magnetic field-related heating is the operative mechanism in the solar corona.

Golub, L.↗
NASA NTRS1994

ROSAT PSPC observations of two X-ray-faint early-type galaxies: NGC 4365 and NGC 4382

We present the results of ROSAT Positive Sensitive Proportional Counter (PSPC) observations of the two early-type galaxies NGC 4365 and NGC 4382. These galaxies are among those observed with Einstein to have the lowest X-ray to optical flux ratios of early-type galaxies. The PSCP data show that for radii r greater than 50 arcsec the radial distributions of the X-ray surface brightness are consistent with the optical distributions of King (1978). We also find that these galaxies have X-ray spectra significantly different from those observed in X-ray-bright ellipticals, with a relative excess of counts detected in the softest spectral channels. This confirms earlier Einstein results. The characteristics of the ROSAT PSPC do not allow us to discriminate between possible spectral models. If we adopt a two-component thermal model on the grounds of physical plausibility, we find that the spectral data can be fitted with a very soft optically thin component, with kT approximately 0.2 keV, and a hard component with kT greater than (1.0-1.5) keV. The hard component has a luminosity consistent with that expected from the integrated emission of a population of low mass-X-ray binaries in these galaxies; the nature of the very soft component is more speculative. Candidates include the coronal emission of late-type stars, supersoft X-ray sources, RS CVn, and perhaps a hot Interstellar Medium (ISM). Alternatively, the spectal data may be fitted with a 0.6-1 keV bremsstrahlung spectrum (expontential plus Gaunt), and may suggest the presence of a totally new population of X-ray sources.

Fabbiano, G.↗
NASA NTRS1973

Observations of noise bands associated with the upper hybrid resonance by the Imp 6 radio astronomy experiment.

The intense noise bands occurring near the upper hybrid resonance frequency have been observed with the Imp 6 GSFC radio astronomy experiment in the plasmasphere. The identification of the upper hybrid resonance provides an accurate measure of the local electron density and allows the observed noise data to be fit to the scale of characteristic frequencies in the plasma. The data are consistent with earlier theoretical interpretations in which noise is generated between the upper hybrid and plasma frequencies and propagates to region 4 of the CMA diagram, where it is reflected at the L = 0 plasma cutoff.

Mosier, S. R.↗
NASA NTRS1975

The correlation of Skylab L-band brightness temperatures with antecedent precipitation

The S194 L-band radiometer flown on the Skylab mission measured terrestrial radiation at the microwave wavelength of 21.4 cm. The terrain emissivity at this wavelength is strongly dependent on the soil moisture content, which can be inferred from antecedent precipitation. For the Skylab data acquisition pass from the Oklahoma panhandle to southeastern Texas on 11 June 1973, the S194 brightness temperatures are highly correlated with antecedent precipitation from the preceding eleven day period, but very little correlation was apparent for the preceding five day period. The correlation coefficient between the averaged antecedent precipitation index values and the corresponding S194 brightness temperatures between 230 K and 270 K, the region of apparent response to soil moisture in the data, was -0.97. The equation of the linear least squares line fitted to the data was: API (cm) = 31.99 -0.114 TB (K), where API is the antecedent precipitation index and TB is the S194 brightness temperature.

Mcfarland, M. J.↗
NASA NTRS1994

Expert system shell to reason on large amounts of data

The current data base management systems (DBMS's) do not provide a sophisticated environment to develop rule based expert systems applications. Some of the new DBMS's come with some sort of rule mechanism; these are active and deductive database systems. However, both of these are not featured enough to support full implementation based on rules. On the other hand, current expert system shells do not provide any link with external databases. That is, all the data are kept in the system working memory. Such working memory is maintained in main memory. For some applications the reduced size of the available working memory could represent a constraint for the development. Typically these are applications which require reasoning on huge amounts of data. All these data do not fit into the computer main memory. Moreover, in some cases these data can be already available in some database systems and continuously updated while the expert system is running. This paper proposes an architecture which employs knowledge discovering techniques to reduce the amount of data to be stored in the main memory; in this architecture a standard DBMS is coupled with a rule-based language. The data are stored into the DBMS. An interface between the two systems is responsible for inducing knowledge from the set of relations. Such induced knowledge is then transferred to the rule-based language working memory.

Giuffrida, Gionanni↗
NASA NTRS1998

Incorporating Endmember Variability into Spectral Mixture Analysis Through Endmember Bundles

Variation in canopy structure and biochemistry induces a concomitant variation in the top-of-canopy spectral reflectance of a vegetation type. Hence, the use of a single endmember spectrum to track the fractional abundance of a given vegetation cover in a hyperspectral image may result in fractions with considerable error. One solution to the problem of endmember variability is to increase the number of endmembers used in a spectral mixture analysis of the image. For example, there could be several tree endmembers in the analysis because of differences in leaf area index (LAI) and multiple scatterings between leaves and stems. However, it is often difficult in terms of computer or human interaction time to select more than six or seven endmembers and any non-removable noise, as well as the number of uncorrelated bands in the image, limits the number of endmembers that can be discriminated. Moreover, as endmembers proliferate, their interpretation becomes increasingly difficult and often applications simply need the aerial fractions of a few land cover components which comprise most of the scene. In order to incorporate endmember variability into spectral mixture analysis, we propose representing a landscape component type not with one endmember spectrum but with a set or bundle of spectra, each of which is feasible as the spectrum of an instance of the component (e.g., in the case of a tree component, each spectrum could reasonably be the spectral reflectance of a tree canopy). These endmember bundles can be used with nonlinear optimization algorithms to find upper and lower bounds on endmember fractions. This approach to endmember variability naturally evolved from previous work in deriving endmembers from the data itself by fitting a triangle, tetrahedron or, more generally, a simplex to the data cloud reduced in dimension by a principal component analysis. Conceptually, endmember variability could make it difficult to find a simplex that both surrounds the data cloud and has vertices that are realistic endmember spectra with reflectances between 0 and 1. In this paper, we create endmember bundles and bounding fraction images for an AVIRIS subscene simulated with a plant canopy radiative transfer model. The simulated subscene is spatially patterned after a subscene from the AVIRIS image acquired August, 1993 over La Copita, Texas. In addition, for comparison, we performed a traditional unmixing with image endmembers.

Bateson, C. Ann↗
NASA NTRS2004

Representation of the Geosynchronous Plasma Environment for Spacecraft Charging Calculations

Historically, our ability to predict and postdict surface charging has suffered from both a lack of reliable secondary emission and backscattered electron yields and poor characterization of the plasma environment. One difficulty lies in the common practice of fitting the plasma data to a Maxwellian or Double Maxwellian distribution function, which may not represent the data well for charging purposes. For 13 years Los Alamos National Laboratory (LANL) has been accumulating measurements of electron and proton spectra from Magnetospheric Plasma Analyzer (MPA) instruments aboard a series of geosynchronous satellites. These data provide both a plasma characterization and the potential of the instrument ground. We use electron and ion flux spectra measured by the LANL MPA to examine how the use of different spectral representations of the charged particle environment in computations of spacecraft potentials during magnetospheric substorms affects the accuracy of the results. We calculate the spacecraft potential using both the measured fluxes and several different fits to these fluxes. These flux measurements and fits have been corrected for the difference between the measured and calculated potential. The potentials computed using the measured fluxes, the best available material properties of graphite carbon, and a secondary electron escape fraction of 81%, are within a factor of three of the measured potential for nearly all the data. Using a Kappa fit to the electron distribution function and a Maxwellian fit to the ion distribution function gives agreement similar to the calculations using the actual data. Alternative spectral representations, including Maxwellian and double Maxwellian for both species, lead to less satisfactory agreement between predicted and measured potentials.

Davis, V. A.↗