Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gamma distribution”

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 163 records · Page 9

Reconciling Contrasting Relationships Between Relative Dispersion and Volume-Mean Radius of Cloud Droplet Size Distributions

Cloud droplet spectral relative dispersion is critical to parameterizations of cloud radiative properties, warm-rain initiation, and aerosol-cloud interactions in models; however, there is no consistent relationship between relative dispersion and volume-mean radius in literature, which hinders improving relative dispersion parameterization and calls for physical explanation. Here we show, by analyzing aircraft observations of cumulus clouds during Routine AAF [Atmospheric Radiation Measurement (ARM) Aerial Facility] Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations, that the correlation between relative dispersion and volume-mean radius changes from positive to negative as volume-mean radius increases. With the new observation, we postulate that the sign of the correlation is determined by whether or not condensation (evaporation) occurs simultaneously with significant new activation (deactivation). The hypothesis is validated by simulations of both an adiabatic cloud parcel model and a parcel model accounting for entrainment-mixing. A new quantity, first bin strength, is introduced to quantify this new observation. Theoretical analysis of truncated gamma and modified gamma size distributions further supports the hypothesis and reconciles the contrasting relationships between relative dispersion and volume-mean radius, including the results in polluted fog observations. The results could shed new light on the so-called “twilight zone” between cloudy and cloud-free air, which in turn affects evaluation of aerosol-cloud interactions and retrieval of aerosol optical depth.

Lu, Chunsong↗

Gamma-Ray Burst Intensity Distributions

We use the lag-luminosity relation to calculate self-consistently the redshifts, apparent peak bolometric luminosities L(sub B1), and isotropic energies E(sub iso) for a large sample of BATSE bursts. We consider two different forms of the lag-luminosity relation; for both forms the median redshift, for our burst database is 1.6. We model the resulting sample of burst energies with power law and Gaussian dis- tributions, both of which are reasonable models. The power law model has an index of a = 1.76 plus or minus 0.05 (95% confidence) as opposed to the index of a = 2 predicted by the simple universal jet profile model; however, reasonable refinements to this model permit much greater flexibility in reconciling predicted and observed energy distributions.

Band, David L.↗

Gamma-ray spectroscopy using angular distribution of Compton scattering

We present a novel method to perform gamma-ray spectroscopy suitable for extreme pileup situations, wherein traditional pulse-height spectroscopy is infeasible or impractical. One example of such application is spectral characterization of intense laser-based radiation sources that employ laser wakefield acceleration, in which a large number of photons can be produced in a well-collimated beam on a picosecond or shorter timescale. The method relies on the angularly-resolved calorimetric measurement of Compton scattering. The probability for a photon to undergo Compton scattering into a given solid angle is a function of its energy and is described by the Klein–Nishina formula. By first Compton scattering the incident beam and then measuring the total energies deposited in multiple detectors placed at various scattering angles, the incident beam energy spectrum can be reconstructed without resolving the individual events in detectors. Here, we compare the performance of an maximum likelihood expectation maximization algorithm with that of an artificial neural network to measure the energy spectrum of a beam of monoenergetic gamma rays from 137Cs decay. We demonstrate that the method can be used to reconstruct the characteristic inverse Compton scattering spectra.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Improvement of Stratospheric Aerosol Extinction Retrieval from OMPS/LP Using a New Aerosol Model

The Ozone Mapping and Profiler Suite Limb Profiler (OMPS/LP) has been flying on the Suomi National Polar-orbiting Partnership (S-NPP) satellite since October 2011. It is designed to produce ozone and aerosol vertical profiles at 2-kilometer vertical resolution over the entire sunlit globe. Aerosol extinction profiles are computed with Mie theory using radiances measured at 675 nanometers. The operational Version 1.0 (V1.0) aerosol extinction retrieval algorithm assumes a bimodal lognormal aerosol size distribution (ASD) whose parameters were derived by combining an in situ measurement of aerosol microphysics with the Stratospheric Aerosol and Gas Experiment (SAGE II) aerosol extinction climatology. Internal analysis indicates that this bimodal lognormal ASD does not sufficiently explain the spectral dependence of LP-measured radiances. In this paper we describe the derivation of an improved aerosol size distribution, designated Version 1.5 (V1.5), for the LP retrieval algorithm. The new ASD uses a gamma function distribution that is derived from Community Aerosol and Radiation Model for Atmospheres (CARMA)-calculated results. A cumulative distribution fit derived from the gamma function ASD gives better agreement with CARMA results at small particle radii than bimodal or unimodal functions. The new ASD also explains the spectral dependence of LP-measured radiances better than the V1.0 ASD. We find that the impact of our choice of ASD on the retrieved extinctions varies strongly with the underlying reflectivity of the scene. Initial comparisons with collocated extinction profiles retrieved at 676 nanometers from the SAGE III instrument on the International Space Station (ISS) show a significant improvement in agreement for the LP V1.5 retrievals. Zonal mean extinction profiles agree to within 10 percent between 19 and 29 kilometers, and regression fits of collocated samples show improved correlation and reduced scatter compared to the V1.0 product. This improved agreement will motivate development of more sophisticated ASDs from CARMA results that incorporate latitude, altitude and seasonal variations in aerosol properties.

Ozone Mapping Profiler Suite Limb Profiler (OMPS/L↗

Simple model for scanning-angle distribution of planetary albedo gamma-rays.

The planetary atmosphere is scanned from a satellite located at a point at an altitude of 500 km from the subsatellite point on the surface (the altitude of the satellite OSO-3). Results point out the possibility that by scanning a planet from the subsatellite point to the limb with a high resolution gamma ray detector on an orbiting satellite, the mean total baryonic scale height may be unfolded from the angular distribution of the gamma-ray flux without having to determine chemical composition and temperature.

Stecker, F. W.↗

An Exact and Near-Exact Distribution Approach to the Behrens–Fisher Problem

The Behrens–Fisher problem occurs when testing the equality of means of two normal distributions without the assumption that the two variances are equal. This paper presents approaches based on the exact and near-exact distributions for the test statistic of the Behrens–Fisher problem, depending on different combinations of even or odd sample sizes. We present the exact distribution when both sample sizes are odd and the near-exact distribution when one or both sample sizes are even. The near-exact distributions are based on a finite mixture of generalized integer gamma (GIG) distributions, used as an approximation to the exact distribution, which consists of an infinite series. The proposed tests, based on the exact and the near-exact distributions, are compared with Welch’s t-test through Monte Carlo simulations, in particular for small and unbalanced sample sizes. The results show that the proposed approaches are competent solutions to the Behrens–Fisher problem, exhibiting precise sizes and better powers than Welch’s approach for those cases. Numerical studies show that the Welch’s t-test tends to be a bit more conservative than the test statistics based on the exact or near-exact distribution, in particular when sample sizes are small and unbalanced, situations in which the proposed exact or near-exact distributions obtain higher powers than Welch’s t-test.

Hong, Serim↗

Radial distribution of cosmic ray intensity in the galaxy from gamma-ray data

The radial distributions of Galactic gamma-ray emissivity and cosmic ray intensity are derived by unfolding the SAS2 line flux data. The results show a marked gradient, with the cosmic ray intensity in the region 3,R,7 kpc being 3 to 4 times as much as locally. There is also an indication for a higher gradient for cosmic rays of lower energy.

Goned, A.↗

Assessment of Ku- and Ka-Band Dual-Frequency Radar for Retrieval of Snow Properties

Current scattering look-up tables for snow assume a constant mass density along with a Gamma particle size distribution (PSD). The first assumption is tested against scattering parameters from simulated particle models generated by Dr. Kwo-Sen Kuo at GSFC and Dr. Guosheng Liu at FSU. Good agreement of the scattering parameters is found with the FSU results if the mass is taken to be the same as the mass of the simulated particle and the mass density is taken to be 0.2 g/cm cu. For the GSFC data base, good agreement is found if the mass density is taken to be between 0.1 and 0.2 g/cm cu. The second assumption of a Gamma PSD is tested against measured PSD's along with a m-d (mass-dimension) relationship. The degree of agreement depends on the value of, the 'shape' parameter in the Gamma distribution but to a lesser degree on the m-d relationship (of the three that were examined). A shortcoming of the simulated snow particle data bases is the lack of large particles. As a consequence, larger values of DFR (dual-frequency ratio) that are commonly seen in airborne and spaceborne measurements cannot be reproduced from the tables. This situation is expected to improve as scattering parameters from larger particles are included in the databases.

Liang Liao↗

Inference of Ice Cloud Properties from High-spectral Resolution Infrared Observations

The theoretical basis is explored for inferring the microphysical properties of ice crystal from high-spectral resolution infrared observations. A radiative transfer model is employed to simulate spectral radiances to address relevant issues. The extinction and absorption efficiencies of individual ice crystals, assumed as hexagonal columns for large particles and droxtals for small particles, are computed from a combination of the finite- difference time-domain (FDTD) technique and a composite method. The corresponding phase functions are computed from a combination of FDTD and an improved geometric optics method (IGOM). Bulk scattering properties are derived by averaging the single- scattering properties of individual particles for 30 particle size distributions developed from in situ measurements and for additional four analytical Gamma size distributions for small particles. The non-sphericity of ice crystals is shown to have a significant impact on the radiative signatures in the infrared (IR) spectrum; the spherical particle approximation for inferring ice cloud properties may result in an overest&ation of the optical thickness and an inaccurate retrieval of effective particle size. Furthermore, we show that the error associated with the use of the Henyey-Greenstein phase function can be as larger as 1 K in terms of brightness temperature for larger particle effective size at some strong scattering wavenumbers. For small particles, the difference between the two phase functions is much less, with brightness temperatures generally differing by less than 0.4 K. The simulations undertaken in this study show that the slope of the IR brightness temperature spectrum between 790-960/cm is sensitive to the effective particle size. Furthermore, a strong sensitivity of IR brightness temperature to cloud optical thickness is noted within the l050-1250/cm region. Based on this spectral feature, a technique is presented for the simultaneous retrieval of the visible optical thickness and effective particle size from high spectral resolution infrared data under ice cloudy con&tion. The error analysis shows that the uncertainty of the retrieved optical thickness and effective particle size has a small range of variation. The error for retrieving particle size in conjunction with an uncertainty of 5 K in cloud'temperature, or a surface temperature uncertainty of 2.5 K, is less than 15%. The corresponding e m r in the uncertainty of optical thickness is within 5-2096, depending on the value of cloud optical thickness. The applicability of the technique is demonstrated using the aircraft-based High- resolution Interferometer Sounder (HIS) data from the Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) in 1996 and the First ISCCP Regional Experiment - Arctic Clouds Experiment (FIRE-ACE) in 1998.

Huang, Hung-Lung↗

Atmospheric gamma ray angle and energy distributions from sea level to 3.5 g/sq cm and 2 to 25 MeV

Differential fluxes of gamma rays were calculated for energies of 2-25 MeV, zenith angles of 0-50 deg and 180-130 deg, and atmospheric depths from nominal sea level, 1000 g/sq cm, to float altitude, 3.5 g/sq cm residual atmosphere. Above 100 g/sq cm growth curves were constructed to estimate the contribution of the extraterrestrial gamma ray flux to the total downward-moving flux, while the upward-moving gamma rays were taken to be strictly of atmospheric origin. Below 100 g/sq cm, all gamma rays originate in the atmosphere. The downward atmospheric flux increases by almost two orders of magnitude between float altitude and the Pfotzer maximum, while the extraterrestrial flux is attenuated exponentially. Gamma rays produced by neutron interactions with the carbon in the scintillator liquid are eliminated by constructing growth curves for downward-moving gamma rays at high altitudes and are negligible compared with downward-moving gamma rays at lower altitudes and upward-moving gamma rays at all altitudes.

Ryan, J. M.↗

Comparison of high energy gamma rays from absolute value of b greater than 30 deg with the galactic neutral hydrogen distribution

High-energy gamma-ray (energy above 35 MeV) data from the SAS 2 satellite have been used to compare the intensity distribution of gamma rays with that of neutral hydrogen (H I) density along the line of sight, at high galactic latitudes (absolute values greater than 30 deg). A model has been constructed for the case where the observed gamma-ray intensity has been assumed to be the sum of a galactic component proportional to the H I distribution plus an isotropic extragalactic emission. A chi-squared test of the model parameters indicates that about 30% of the total high-latitude emission may originate within the Galaxy.

Ozel, M. E.↗

Parameter Estimation in Astronomy with Poisson-Distributed Data: The (CHI)2(gamma) Statistic - 1

Applying the standard weighted mean formula, [Sigma (sub i)n(sub i)ssigma(sub i, sup -2)], to determine the weighted mean of data, n(sub i), drawn from a Poisson distribution, will, on average, underestimate the true mean by approx. 1 for all true mean values larger than approx.3 when the common assumption is made that the error of the i th observation is sigma(sub i) = max square root of n(sub i), 1).This small, but statistically significant offset, explains the long-known observation that chi-square minimization techniques which use the modified Neyman'chi(sub 2) statistic, chi(sup 2, sub N) equivalent Sigma(sub i)((n(sub i) - y(sub i)(exp 2)) / max(n(sub i), 1), to compare Poisson - distributed data with model values, y(sub i), will typically predict a total number of counts that underestimates the true total by about 1 count per bin. Based on my finding that weighted mean of data drawn from a Poisson distribution can be determined using the formula [Sigma(sub i)[n(sub i) + min(n(sub i), 1)](n(sub i) + 1)(exp -1)] / [Sigma(sub i)(n(sub i) + 1)(exp -1))], I propose that a new chi(sub 2) statistic, chi(sup 2, sub gamma) equivalent, should always be used to analyze Poisson- distributed data in preference to the modified Neyman's chi(exp 2) statistic. I demonstrated the power and usefulness of,chi(sub gamma, sup 2) minimization by using two statistical fitting techniques and five chi(exp 2) statistics to analyze simulated X-ray power - low 15 - channel spectra with large and small counts per bin. I show that chi(sub gamma, sup 2) minimization with the Levenberg - Marquardt or Powell's method can produce excellent results (mean slope errors approx. less than 3%) with spectra having as few as 25 total counts.

Mighell, Kenneth J.↗

Implications of the BATSE data for a helocentric origin of gamma-ray bursts

Several possible models have been suggested to explain the observed distribution of gamma-ray bursts: heliocentric distributions such as the Oort cloud, large galactic halos, and cosmological models. We report here on an investigation into the implications of the Burst and Transient Source Experiment (BATSE) gamma-ray burst distribution (Meegan et al. 1992a) data on the possible helocentric origin of gamma-ray bursts. We find no statistically significant anisotropy in the angular distribution of the bursts in a Sun-referenced coordinate system; there is no dipole moment in the direction of the Sun, and no quardrupole moment associated with the ecliptic plane. We have employed direct analytic calculations and Monte Carlo simulations of sources in the Oort cloud to constrain possible helicentric burst distributions. These can produce distributions consistent with the observed angular isotropy, the meal value of V/V(sub max), and the observed C/C(sub min) distribution of BATSE, and provide limits to burst energy of a few times approximately 10(exp 27) ergs. However, the agreement of the heliocentric C/C(sub min) distributions with the BATSE data is attributable to the relatively limited sampling of strong, nearby bursts. These bursts are known from observation to be homogeneously distributed, yet the density of sources in the Oort cloud is not constant in this region. Integral number-intensity distributions from the Oort cloud for larger numbers of bursts cannot reproduce the known homogeneity of the strong bursts without modification to the computed cometary number density and are therefore unlikely explanations of the gamma-ray burst distribution.

Horack, J. M.↗

Compton Gamma Ray Observatory Guest Investigator Program

This paper presents a final report for the Compton Gamma Ray Observatory Guest Investigator Program from 06/01/91-07/31/97. The topics include: 1) Solar Flare Neutron Spectra and Accelerated Ions; 2) Gamma Ray Lines From The Orion Complex; 3) Implications of Nuclear Line Emission From The Orion Complex; 4) Possible Sites of Nuclear Line Emission From Massive OB Associations; 5) Gamma-Ray Burst Repitition and BATSE Position Uncertainties; 6) Effects of Compton Scattering on BATSE Gamma-Ray Burst Spectra; and 7) Selection Biases on the Spectral and Temporal Distribution of Gamma Ray Bursts.

Lingenfelter, Richard E.↗

Effective Ice Particle Densities for Cold Anvil Cirrus

This study derives effective ice particle densities from data collected from the NASA WB-57F aircraft near the tops of anvils during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers (CRYSTAL) Florida Area Cirrus Experiment (FACE) in southern Florida in July 2002. The effective density, defined as the ice particle mass divided by the volume of an equivalent diameter liquid sphere, is obtained for particle populations and single sizes containing mixed particle habits using measurements of condensed water content and particle size distributions. The mean effective densities for populations decrease with increasing slopes of the gamma size distributions fitted to the size distributions. The population-mean densities range from near 0.91 g/cu m to 0.15 g/cu m. Effective densities for single sizes obey a power-law with an exponent of about -0.55, somewhat less steep than found from earlier studies. Our interpretations apply to samples where particle sizes are generally below 200-300 microns in maximum dimension because of probe limitations.

Heymsfield, Andrew J.↗

Surrogate Distributed Radiological Sources—Part III: Quantitative Distributed Source Reconstructions

In this third part of a multi-paper series, we present quantitative image reconstruction results from aerial measurements of eight different surrogate distributed gamma-ray sources on flat terrain. Here, we show that our quantitative imaging methods can accurately reconstruct the expected shapes, and, after appropriate calibration, the absolute activity of the distributed sources. We conduct several studies of imaging performance versus various measurement and reconstruction parameters, including detector altitude and raster pass spacing, data and modeling fidelity, and regularization type and strength. The imaging quality performance is quantified using various quantitative image quality metrics. Our results confirm the utility of point source arrays as surrogates for truly distributed radiological sources, and advance the quantitative capabilities of Scene Data Fusion gamma-ray imaging methods.

Airborne survey↗

Primary $\gamma$-ray intensities and $\gamma$-strength functions from discrete two-step $\gamma$-ray cascades in radiative proton-capture experiments

Background: Reaction rates of radiative capture reactions can play a crucial role in the nucleosynthesis of heavy nuclei in explosive stellar environments. These reaction rates depend strongly on $\gamma$-ray decay widths in the reaction products, which are, for nonresonant capture reactions at high excitation energies, derived from the $\gamma$-ray strength function and the nuclear level density. Recently, the ratio method was applied to primary $\gamma$ rays observed from ($\textit{d, p}$) reactions and nuclear resonance fluorescence measurements to extract the dipole strength in atomic nuclei and to test the generalized Brink-Axel hypothesis. Purpose: The purpose of this work is to apply the ratio method to primary $\gamma$-ray intensities of the 63,65 Cu($\textit{p},\gamma$) reactions to extract $\gamma$-ray strength information on the nuclei 64,66 Zn. Here, the impact of spin distribution, total $\gamma$-ray decay widths, level densities, and width fluctuations on the application of the ratio method will be discussed. Additionally, by comparing the relative $\gamma$-ray strength at different excitation energies, conclusions on the validity of the generalized Brink-Axel hypothesis can be made. Method: The radiative proton capture reaction measurements have been performed at the HORUS $\gamma$-ray spectrometer of the University of Cologne at one excitation energy for each reaction. Primary $\gamma$-ray intensities have been determined by normalizing secondary $\gamma$-ray transitions in two-step cascades using their absolute branching ratio. The ratio method was applied to the measured primary $\gamma$-ray intensities as well as to previous measurements by Erlandsson et al. at different excitation energies. Results: The relative strength function curve for 64 Zn from our measurement shows no significant deviation from the previous measurement at a different excitation energy. The same is true for 66 Zn where both measurements were at almost the same excitation energy. Absolute $\gamma$-strength function values have been obtained by normalizing the relative curves to quasiparticle random phase approximation calculations because of the absence of experimental data in the respective energy region. Conclusion: The generalized Brink-Axel hypothesis, i.e., the independence of the strength function on the excitation energy, seems to hold in the studied energy region and nuclei. The method to obtain primary $\gamma$-ray intensities from two-step cascade spectra was shown to be a valuable and sensitive tool although its uncertainties are connected to the knowledge of the low-energy level scheme of the investigated nucleus. The scaling in the ratio method should be taken with care, because the relative strength is not a simple sum of $f_{E1}$ and $f_{M1}$ but a somewhat complex linear combination dependent on the excitation energy of the nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗