Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correlation histogram”

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 91 records · Page 5

Elemental composition of individual chondrules from ordinary chondrites.

Sequential nondestructive neutron activation analysis was used to determine the bulk abundance of Fe, Al, Na, Mn, Cr, Sc, Co, and Ir in approximately 300 individual chondrules from 16 chondrites representing the H(3-5), L4, and LL(3-6) compositional and petrologic classes. The histograms indicate that the most probable abundances for lithophilic elements, except Cr, are enriched in the chondrules, while the siderophilic elements are depleted in the chondrules compared to the whole chondrite. Some of the abundance populations, such as Al and Fe, appear to be multimodal. Systematic variations in the composition of the chondrules with increasing petrologic type were observed; most consistent are an increasing Na-Al and Cr-Al correlation, a decreasing Na-Mn correlation, increasing Na abundance, and decreasing Na and Mn dispersions among chondrules.

Osborn, T. W.↗

A comparative statistical study of long-term agroclimatic conditions affecting the growth of US winter wheat: Distributions of regional monthly average precipitation on the Great Plains and the state of Maryland and the effect of agroclimatic conditions on yield in the state of Kansas

A histogram analysis of average monthly precipitation over 30 and 84 year periods for both Maryland and Kansas was made and the results compared. A second analysis, a statistical assessment of the effect of average monthly precipitation on Kansas winter wheat yield was made. The data sets covered the three periods of 1941-1970, 1887-1970, and 1887-1921. Analyses of the limited data sets used (only the average monthly precipitation and temperature were correlated against yield) indicated that fall precipitation values, especially those of September and October, were more important to winter wheat yield than were spring values, particularly for the period 1941-1970.

Welker, J.↗

Generic Signatures of the Time Profiles of BATSE Cosmic Gamma-Ray Bursts

A new method is proposed, which allows the study of generic signatures of cosmic gamma-ray burst time histories. We average the 64 ms resolution time profiles of 275 bright bursts detected by BATSE. The profile of each burst is normalized by the maximum number of counts at the peak of the primary pulse, and individual pulses and interpulse valleys are selected from the normalized profiles by identical selection criteria. New generic temporal parameters are introduced, which characterize the duration and equivalent width of each pulse and the duration of each valley. The histograms of the total equivalent pulse width and summed pulse duration are bimodal. Bimodality is also seen in the histogram of the mean duration of individual pulses. Bursts from the short and long peaks of these distributions correspond to the two modes of the Third BATSE Burst Catalog T(sub 50) and T(sub 90) distributions. Therefore, these new burst parameters demonstrate that the observed bimodal temporal behavior results from properties of the pulsed emission of gamma-ray bursts. The long mode of the T(sub 90) histogram includes bursts with from one to approximately 20 pulses; the logarithmic mean pulse duration is 1.17 +/- 0.09 S; for the long events with more than one pulse, the logarithmic mean valley duration is 1.28 +/- 0.15 S. Bursts of the short mode of T(sub 90) are mainly single-pulse events, and the logarithmic mean pulse duration is much smaller, 0.20 +/- 0.01 s. For multipulse bursts of the T(sub 90) long mode, marginal correlations were found between the parameters of the pulses and valleys and the number of pulses. The basic signatures of the evolution of pulses and valleys along the time course of bursts are examined. Conclusions are drawn concerning the physics of gamma-ray emission by taking into account these signatures.

Mitrofanov, Igor G.↗

SHARPEN: A Scheme to Restore the Distribution of Averaged Precipitation Fields

A key strategy in obtaining complete global coverage of high-resolution precipitation is to combine observations from multiple fields, such as the intermittent passive microwave observations, precipitation propagated in time using motion vectors, and geosynchronous infrared observations. These separate precipitation fields can be combined through weighted averaging, which produces estimates that are generally superior to the individual parent fields. However, the process of averaging changes the distribution of the precipitation values, leading to an increase in precipitating area and a decrease in the values of high precipitation rates, a phenomenon observed in IMERG. To mitigate this issue, we introduce a new scheme called SHARPEN (Scheme for Histogram Adjustment with Ranked Precipitation Estimates in the Neighborhood), which recovers the distribution of the averaged precipitation field based on the idea of quantile mapping applied to the local environment. When implemented in IMERG, precipitation estimates from SHARPEN exhibit a distribution that resembles that of the original instantaneous observations, with matching precipitating area and peak precipitation rates. Case studies demonstrate its improved ability in bridging between the parent precipitation fields. Evaluation against ground observations reveals a distinct improvement in precipitation detection skill, but also a slightly reduced correlation likely because of a sharper precipitation field. The increased computational demand of SHARPEN can be mitigated by striding over multiple grid boxes, which has only marginal impacts on the accuracy of the estimates. SHARPEN can be applied to any precipitation algorithm that produces an average from multiple input precipitation fields and is being considered for implementation in IMERG V07.

Jackson Tan↗

Crab pulsar giant pulses: Simultaneous radio and GRO observations

Observations are reported of the Crab pulsar made at radio frequencies concurrent with Oriented Scintillation Spectrometer Experiment (OSSE) observations from 15 to 27 May 1991. Using the 43 m telescope at Green Bank at 0.8 and 1.4 GHz, samples were made continuously for 10 hrs/day at intervals of 100 to 300 microsecs. The analysis of the radio data includes calculation of histograms of pulse intensities, absolute timing to about 20 microsec precision, and characterization of intensity variations on time scales from the 33 ms spin period to days. The most detailed analysis is presented made of giant pulses. The ultimate goal is to bin the radio data into giant and nongiant pulses and to form average waveforms of OSSE data for the corresponding pulse periods. A test is done to see whether the violet radio fluctuations (which are not seen in other radio pulsars to the same degree) are correlated with low energy gamma rays, yielding constraints on the radio coherence mechanism and the steadiness of the electron-positron outflow in the magnetosphere. Timing analysis of the radio data provides a well defined ephemeris over the specified range of epochs. The gamma ray pulse phase was predicted with an error of less than 70 microsecs.

Lundgren, Scott C.↗

Fast correlation function calculator: A high-performance pair-counting toolkit

A novel high-performance exact pair-counting toolkit called fast correlation function calculator (FCFC) is presented. With the rapid growth of modern cosmological datasets, the evaluation of correlation functions with observational and simulation catalogues has become a challenge. High-efficiency pair-counting codes are thus in great demand. We introduce different data structures and algorithms that can be used for pair-counting problems, and perform comprehensive benchmarks to identify the most efficient algorithms for real-world cosmological applications. We then describe the three levels of parallelisms used by FCFC, SIMD, OpenMP, and MPI, and run extensive tests to investigate the scalabilities. Finally, we compare the efficiency of FCFC with alternative pair-counting codes. The data structures and histogram update algorithms implemented in FCFC are shown to outperform alternative methods. FCFC does not benefit greatly from SIMD because the bottleneck of our histogram update algorithm is mainly cache latency. Nevertheless, the efficiency of FCFC scales well with the numbers of OpenMP threads and MPI processes, even though speedups may be degraded with over a few thousand threads in total. FCFC is found to be faster than most (if not all) other public pair-counting codes for modern cosmological pair-counting applications.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation of the two-region Rossi-alpha model for reflected assemblies

Nuclear criticality safety, nonproliferation and safeguards, emergency response, and stockpile stewardship utilize estimates of the k eff multiplication factor. The value of k eff cannot be directly measured, but it can be inferred from the prompt neutron period. One modality of measuring the prompt neutron period is the Rossi-alpha method, which is predicated on fitting a histogram of detection time differences due to the nonrandom temporal distribution of same-fission-chain neutrons. Recent works have developed the motivation and theory to expand traditional one-region point kinetic models resulting in one-exponential histogram fits to two-region models that result in two-exponential fits. This paper validates the new two-region model using organic scintillator measurements of copper-reflected weapons-grade plutonium (0:83 ≤ k eff ≤ 0:94) and high-density- polyethylene-reflected highly enriched uranium (0:73 ≤ k eff ≤ 0:95). Furthermore, the results show that more thermal systems have shorter prompt neutron periods in the core region due to increased induced-fission probabilities for moderated neutrons. Here, a new parameter introduced by the two- region model is also shown to be correlated to the amount of reflection, and may be used to infer assembly properties such as type and amount of reflector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Einstein soft X-ray survey of the Pleiades

The results of a 0.1-4.5-keV X-ray survey of a 2 x 2-deg area centered on the Pleiades open cluster, performed using the imaging proportional counter and high-resolution imager of the Einstein Observatory on four days in 1980-1981, are presented in extensive tables, graphs, maps, histograms, and finding charts and characterized. A total of 61 sources are detected, and 44 of these are identified with cluster members of spectral types B-M. Findings discussed include Lx/Lbol of 10 the -7th for early-type stars; F-star mean Lx like that of F stars in the Hyades and in the field (denying time evolution of Lx); mean G-star Lx = 3.7 x 10 to the 29th erg/sec; G-star activity decay more gradual than 1/sq rt t, with sharp fall for t greater than 1 Gyr; and no evidence for X-ray-emissivity/rotational-velocity correlation in a homogeneous sample of K dwarfs or a sample of K and M stars with established V sin i data.

Caillault, J.-P.↗

Characterizing Oscillatory Bursts in Single-Trial EEG Data

Oscillatory bursts in numerous bands ranging from low (theta) to high frequencies (e.g., gamma) undoubtedly play an important role in cortical dynamics. Largely because of the inadequacy of existing analytic techniques. however, oscillatory bursts and their role in cortical processing remains poorly understood. To study oscillatory bursts effectively one must be able to isolate them and characterize them in the single trial. We describe a series of straightforward analysis techniques that produce useful indices of burst characteristics. First, stimulus-evoked responses are estimated using Differentially Variable Component Analysis (dVCA), and are subtracted from the single-trial. The single-trial characteristics of the evoked responses are stored to identify possible correlations with burst activity. Time-frequency (T-F), or wavelet, analyses are then applied to the single trial residuals. While T-F plots have been used in recent studies to identify and isolate bursts, we go further by fitting each burst in the T-F plot with a two-dimensional Gaussian. This provides a set of burst characteristics, such as, center time. burst duration, center frequency. frequency dispersion. and amplitude, all of which contribute to the accurate characterization of the individual burst. The burst phase can also be estimated. Burst characteristics can be quantified with several standard techniques (e.g.. histogramming and clustering), as well as Bayesian techniques (e.g., blocking) to allow a more parametric description analysis of the characteristics of oscillatory bursts, and the relationships of specific parameters to cortical excitability and stimulus integration.

Knuth, K. H.↗

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗

C 4 fluxes from the sun as a star and the correlation with magnetic flux

A total of 144 C 4 wavelength 1548 SMM-UVSP spectroheliograms of solar plages were analyzed, some of which are series of exposures of the same region on the same day. Also analyzed were C 4 wavelength 1551 rasters of plages and C 4 1548 rasters of the quiet sun. The sample contains data on 17 different plages, observed on 50 different days. The center-to-limb variations of the active regions show that the optical thickness effects in the C 4 wavelength 1548 line can be neglected in the conversion from intensity to flux density. As expected for the nearly optically thin situation, the C 4 1548 line is twice as bright as the C 4 wavelength 1551 line. The average C 4 wavelength 1548 flux density for a quiet is 2700 erg/cm/s and, with surprisingly little scatter, 18,000 erg/cm/s for plages. The intensity histograms of rasters obtained at disk centers can be separated into characteristic plage and quiet-sun contributions with variable relative filling factors. The disk-averaged flux density in the C 4 doublet and the disk-averaged magnitude of the magnetic flux density are related. The relationship between the C 4 and magnetic flux densities for spatially resolved data is inferred to be almost the same, with only an additional factor of order unity in the constant of proportionality.

Schrijver, C. J.↗

C IV fluxes from the Sun as a star, and the correlation with magnetic flux

A total of 144 C IV wavelength 1548 Solar Maximum Mission (SMM)-UVSP spectroheliograms of solar plages were analyzed, some of which are series of exposures of the same region on the same day. Also analyzed were the C IV wavelength 1551 rasters of plages and C IV wavelength 1548 rasters of the quiet sun. The sample contained data on 17 different plages, observed on 50 different days. The center-to-limb variations of the active regions show that the optical thickness effects in the C IV wavelength 1548 line can be neglected in the conversion from intensity to flux density. As expected for the nearly optically thin situation, the C IV wavelength 1548 line is twice as bright as the C IV 1551 line. The average C IV wavelength 1548 flux density for a quiet region is 2700 ergs/cm/s and, with surprisingly little scatter, 18,000 erg/cm/s for plages. The intensity histograms of rasters obtained at disk center can be separated into characteristic plage and quiet sun contributions with variable relative filling factors. The relationship between the C IV and magnetic flux densities for spatially resolved data is inferred to be almost the same, with only an additional factor of order unity in the constant of proportionality.

Schrijver, C. J.↗

Technical note: Recommendations for diagnosing cloud feedbacks and rapid cloud adjustments using cloud radiative kernels

Abstract. The cloud radiative kernel method is a popular approach to quantify cloud feedbacks and rapid cloud adjustments to increased CO2 concentrations and to partition contributions from changes in cloud amount, altitude, and optical depth. However, because this method relies on cloud property histograms derived from passive satellite sensors or produced by passive satellite simulators in models, changes in obscuration of lower-level clouds by upper-level clouds can cause apparent low-cloud feedbacks and adjustments, even in the absence of changes in lower-level cloud properties. Here, we provide a methodology for properly diagnosing the impact of changing obscuration on cloud feedbacks and adjustments and quantify these effects across climate models. Averaged globally and across global climate models, properly accounting for obscuration leads to weaker positive feedbacks from lower-level clouds and stronger positive feedbacks from upper-level clouds while simultaneously removing a mostly artificial anti-correlation between them. Given that the methodology for diagnosing cloud feedbacks and adjustments using cloud radiative kernels has evolved over several papers, and obscuration effects have only occasionally been considered in recent papers, this paper serves to establish recommended best practices and to provide a corresponding code base for community use.

54 ENVIRONMENTAL SCIENCES↗

Satellite Remote Sensing of Tropical Precipitation and Ice Clouds for GCM Verification

This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.

Evans, K. Franklin↗

The effects of cloud inhomogeneities upon radiative fluxes, and the supply of a cloud truth validation dataset

The ASTER polar cloud mask algorithm is currently under development. Several classification techniques have been developed and implemented. The merits and accuracy of each are being examined. The classification techniques under investigation include fuzzy logic, hierarchical neural network, and a pairwise histogram comparison scheme based on sample histograms called the Paired Histogram Method. Scene adaptive methods also are being investigated as a means to improve classifier performance. The feature, arctan of Band 4 and Band 5, and the Band 2 vs. Band 4 feature space are key to separating frozen water (e.g., ice/snow, slush/wet ice, etc.) from cloud over frozen water, and land from cloud over land, respectively. A total of 82 Landsat TM circumpolar scenes are being used as a basis for algorithm development and testing. Numerous spectral features are being tested and include the 7 basic Landsat TM bands, in addition to ratios, differences, arctans, and normalized differences of each combination of bands. A technique for deriving cloud base and top height is developed. It uses 2-D cross correlation between a cloud edge and its corresponding shadow to determine the displacement of the cloud from its shadow. The height is then determined from this displacement, the solar zenith angle, and the sensor viewing angle.

Welch, Ronald M.↗

Calculation of design load for the MOD-5A 7.3 mW wind turbine system

Design loads are presented for the General Electric MOD-SA wind turbine. The MOD-SA system consists of a 400 ft. diameter, upwind, two-bladed, teetered rotor connected to a 7.3 mW variable-speed generator. Fatigue loads are specified in the form of histograms for the 30 year life of the machine, while limit (or maximum) loads have been derived from transient dynamic analysis at critical operating conditions. Loads prediction was accomplished using state of the art aeroelastic analyses developed at General Electric. Features of the primary predictive tool - the Transient Rotor Analysis Code (TRAC) are described in the paper. Key to the load predictions are the following wind models: (1) yearly mean wind distribution; (2) mean wind variations during operation; (3) number of start/shutdown cycles; (4) spatially large gusts; and (5) spatially small gusts (local turbulence). The methods used to develop statistical distributions from load calculations represent an extension of procedures used in past wind programs and are believed to be a significant contribution to Wind Turbine Generator analysis. Test/theory correlations are presented to demonstrate code load predictive capability and to support the wind models used in the analysis. In addition MOD-5A loads are compared with those of existing machines. The MOD-5A design was performed by the General Electric Company, Advanced Energy Program Department, under Contract DEN3-153 with NASA Lewis Research Center and sponsored by the Department of Energy.

Mirandy, L.↗

Evidence for long-term variability in the ultra high energy photon flux from Cygnus X-3

A time-correlation analysis of atmospheric Cerenkov pulses by a wide-angle photomultiplier system was previously shown to have present in it a nonrandom component which seemed associated with the Right Ascension (RA) range approx. 20+or-04h. A recent examination of multi-muon events recorded by a photon-decay detector shows a similar time-dependent effect, closely matching the previous results, supporting the suggestion that the effect is of cosmic origin. However, even though Cyg. X-3 lies well inside the region of peak intensity, it does not seem possible to ascribe to it the whole effect, for the implied photon flux appears too large to be reconciled to various gamma-ray measurements of Cyg. X-3. The original data were subjected to a phase-histogram analysis and it as found that only 2.5% of overall recorded data are compatible with a phase-dependent emission from Cyg. X-3. Assuming these events to be gamma rays yields a detected flux of (2.6 + or - 0.3) x 10 to the minus 12th power gamma cm -2s-1 above 5 x 10 to the 14th power eV. Comparing this value with more recent ultra high energy (UHE) photon data from the same source, it is suggested that the available data generally favor a long-term reduction in the Cyg. X-3 inferred luminosity ( 10 to the 13th power eV) by a factor of (1.8 + or - 0.3) per year.

Bhat, C. L.↗

Non-Gaussian Ensemble Filtering and Adaptive Inflation for Soil Moisture Data Assimilation

The rank histogram filter (RHF) and the ensemble Kalman filter (EnKF) are assessed for soil moisture estimation using perfect model (identical twin) synthetic data assimilation experiments. The primary motivation is to gauge the impact on analysis quality attributable to the consideration of non-Gaussian forecast error distributions. Using the NASA Catchment land surface model, the two filters are compared at 18 globally distributed single-catchment locations for a 10-yr experiment period. It is shown that both filters yield adequate estimates of soil moisture, with the RHF having a small but significant performance advantage. Most notably, the RHF consistently increases the normalized information contribution (NIC) score of the mean absolute bias by 0.05 over that of the EnKF for surface, root-zone, and profile soil moisture. The RHF also increases the NIC score for the anomaly correlation of surface soil moisture by 0.02 over that of the EnKF (at a 5% significance level). Results additionally demonstrate that the performance of both filters is somewhat improved when the ensemble priors are adaptively inflated to offset the negative effects of systematic errors.

Rolf Reichle↗