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 73 records · Page 4

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗

Statistical characteristics of laser interference and its effect on IR optoelectronic observation systems

We consider statistical characteristics of laser interference recorded in the focal region of an IR optoelectronic observation system under field and out-of-field illumination by laser radiation. The applicability of the analytical description of experimental histograms of the laser interference signal distribution at the photodetector array output using the probability density approximations in the form of gamma distribution and the Gaussian function is substantiated. It is shown that the specific form of these functions is determined both by the mean value (mathematical expectation) of the interference signal and by the characteristic parameter M, which depends on the statistical properties of the laser interference. (laser applications and other issues in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimation of mean rain rate - Application to satellite observations

A method for the estimation of the mean area average rain rate from dependent data is developed and applied to the GARP Atlantic Tropical Experiment data. The method consists of fitting a mixed distribution, containing an atom at zero, by minimum chi-square in combination with certain time-space sampling designs. In modeling the continuous component of the mixed distribution, it is shown that the lognormal distribution provides a very close fit for the nonzero area average rainrates. A comparison with the gamma distribution shows that the lognormal distribution is a better choice as expressed by the minimum chi-square criterion. Some of the time-space sampling designs correspond to satellite sampling. The results indicate that a satellite visiting an area of about 350 x 350 sq km in the tropics approximately every 10 hours over a period can provide a rather close estimate for the mean area average rain rate.

Kedem, Benjamin↗

A probability distribution model for rain rate

A systematic approach is suggested for modeling the probability distribution of rain rate. Rain rate, conditional on rain and averaged over a region, is modeled as a temporally homogeneous diffusion process with appropiate boundary conditions. The approach requires a drift coefficient-conditional average instantaneous rate of change of rain intensity-as well as a diffusion coefficient-the conditional average magnitude of the rate of growth and decay of rain rate about its drift. Under certain assumptions on the drift and diffusion coefficients compatible with rain rate, a new parametric family-containing the lognormal distribution-is obtained for the continuous part of the stationary limit probability distribution. The family is fitted to tropical rainfall from Darwin and Florida, and it is found that the lognormal distribution provides adequate fits as compared with other members of the family and also with the gamma distribution.

Kedem, Benjamin↗

Vector wind profile gust model

A methodology was developed for the derivation and analysis of small scale perturbations in Jimsphere wind profiles. Gusts in various wavelength bands have been derived from these perturbations; the probability distribution of gust components and associated gust length has been shown to be accurately represented by a gamma distribution. Theoretical and observed distributions of component gust vary significantly with season, altitude, and wavelength range. The results of this study provide the basis for a vector wind model for Cape Kennedy, Florida.

Adelfang, S. I.↗

Parameterizing Grid-Averaged Longwave Fluxes for Inhomogeneous Marine Boundary Layer Clouds

This paper examines the relative impacts on grid-averaged longwave flux transmittance (emittance) for Marine Boundary Layer (MBL) cloud fields arising from horizontal variability of optical depth tau and cloud sides, First, using fields of Landsat-inferred tau and a Monte Carlo photon transport algorithm, it is demonstrated that mean all-sky transmittances for 3D variable MBL clouds can be computed accurately by the conventional method of linearly weighting clear and cloudy transmittances by their respective sky fractions. Then, the approximations of decoupling cloud and radiative properties and assuming independent columns are shown to be adequate for computation of mean flux transmittance. Since real clouds have nonzero geometric thicknesses, cloud fractions A'(sub c) presented to isotropic beams usually exceed the more familiar vertically projected cloud fractions A(sub c). It is shown, however, that when A(sub c)less than or equal to 0.9, biases for all-sky transmittance stemming from use of A(sub c) as opposed to A'(sub c) are roughly 2-5 times smaller than, and opposite in sign to, biases due to neglect of horizontal variability of tau. By neglecting variable tau, all-sky transmittances are underestimated often by more than 0.1 for A(sub c) near 0.75 and this translates into relative errors that can exceed 40% (corresponding errors for all-sky emittance are about 20% for most values of A(sub c). Thus, priority should be given to development of General Circulation Model (GCM) parameterizations that account for the effects of horizontal variations in unresolved tau, effects of cloud sides are of secondary importance. On this note, an efficient stochastic model for computing grid-averaged cloudy-sky flux transmittances is furnished that assumes that distributions of tau, for regions comparable in size to GCM grid cells, can be described adequately by gamma distribution functions. While the plane-parallel, homogeneous model underestimates cloud transmittance by about an order of magnitude when 3D variable cloud transmittances are less than or equal to 0.2 and by approx. 20% to 100% otherwise, the stochastic model reduces these biases often by more than 80%.

Barker, Howard W.↗

To What Extent Can Raindrop Size Be Determined by a Multiple-Frequency Radar?

In this paper, an analytical treatment of the atmospheric remote sensing problem of determining the raindrop size distribution (DSD) with a spaceborne multifrequency microwave nadir-looking radar system is presented. It is typically assumed that with two radar measurements at different frequencies one ought to be able to calculate two state variables of the DSD: a bulk quantity, such as the rain rate, and a distribution shape parameter. To determine if this nonlinear problem can indeed be solved, the DSD is modeled as a Gamma distribution and quadratic approximations to the corresponding radar-rain relations are used to examine the invertibility of the resulting system of equations in the case of two as well as three radar frequencies. From the investigation, it is found that for regions of DSD state space multiple solutions exist for two or even three different frequency radar measurements. This should not be surprising given the nonlinear coupled nature of the problem.

rain characteristics↗

Rain Drop Size Distributions Estimated from NOAA Snow-Level Radar Data

Using NOAA’s S-band High-Power Snow-Level Radar (HPSLR), a technique for estimating the rain drop size distribution (DSD) above the radar is presented. This technique assumes the DSD can be described by a four parameter, generalized gamma distribution (GGD). Using the radar’s measured average Doppler velocity spectrum and a value (assumed, measured, or estimated) of the vertical air motion w, an estimate of the GGD is obtained. Four different methods can be used to obtain w. One method that estimates a mean mass-weighted raindrop diameter D m from the measured reflectivity Z produces realistic DSDs compared to prior literature examples. These estimated DSDs provide evidence that the radar can retrieve the smaller drop sizes constituting the “drizzle” mode part of the DSD. Here, this estimation technique was applied to 19 h of observations from Hankins, North Carolina. Results support the concept that DSDs can be modeled using GGDs with a limited range of parameters. Further work is needed to validate the described technique for estimating DSDs in more varied precipitation types and to verify the vertical air motion estimates.

54 ENVIRONMENTAL SCIENCES↗

Inferring the spatial and energy distribution of gamma-ray burst sources. 1: Methodology

We describe Bayesian methods for analyzing the distribution of gamma-ray burst peak photon fluxes and directions. These methods fit the differential distribution, and have the following advantages over rival methods: (1) they do not destroy information by binning or averaging the data (as do, say, chi squared, the averaged value of V/V(sub max), and angular moment analyses); (2) they straightforwardly handle uncertainties in the measured quantities; (3) they analyze the strength and direction information jointly; (4) they use information available about nondetections; and (5) they automatically identify and account for biases and selection effects given a precise description of the experiment. In these methods, the most important information needed about the instrument threshold is not its value at the times of burst triggers, as is used in the average value of V/V(sub max) analyses, but rather the value of the threshold at times when no trigger occurred. We show that this information can be summarized as an average detection efficiency that is similar to the product of the exposure and efficiency reported in the First Burst and Transient Source Experiment (BATSE) Burst (1B) Catalog, but significantly different from it at low fluxes. We also quantify an important bias that results from estimating the peak flux by scanning the burst to find the peak number of counts in a window of specified duration, as was done for the 1B Catalog. When the duration of the peak of the light curve is longer than the window duration, a simple flux estimate based on the peak counts significantly overestimates the peak flux in a nonlinear fashion that distorts the shape of the log(N)-log(P) distribution. This distortion also corrupts analyses of the V/V(sub max) distribution that use ratios of counts above background to estimate V/V(sub max). The Bayesian calculation specifies how to account for this bias. Implementation of the Bayesian approach requires some changes in the way burst data are reported that we describe in detail. Subsequent papers will report analyses of the 1B Catalog data using the methods described here.

Loredo, Thomas J.↗

Spatial modelling of total storm rainfall

The spatial structure of the depth of rainfall from a stationary storm event is investigated by using point process techniques. Cells are assumed to be stationary and to be distributed in space either independently according to a Poisson process, or with clustering according to a Neyman-Scott scheme. Total storm rainfall at the centre of each cell is a random variable and rainfall is distributed around the centre in a way specified by a spread function that may incorporate random parameters. The mean, variance and covariance structure of the precipitation depth at a point are obtained for different spread functions. For exponentially distributed centre depth and a spread function having quadratically exponential decay, the total storm depth at any point in the field is shown to have a gamma distribution. The probability of zero rainfall at a point is investigated, as is the stochastic variability of model parameters from storm to storm. Data from the Upper Rio Guaire basin in Venezuela are used in illustration.

Rodriguez-Iturbe, I.↗

Newton/Poisson-Distribution Program

NEWTPOIS, one of two computer programs making calculations involving cumulative Poisson distributions. NEWTPOIS (NPO-17715) and CUMPOIS (NPO-17714) used independently of one another. NEWTPOIS determines Poisson parameter for given cumulative probability, from which one obtains percentiles for gamma distributions with integer shape parameters and percentiles for X(sup2) distributions with even degrees of freedom. Used by statisticians and others concerned with probabilities of independent events occurring over specific units of time, area, or volume. Program written in C.

Bowerman, Paul N.↗

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↗

Minimizing Timing Jitter’s Impact on Ground-Penetrating Radar Array Coupling Signals

This article presents a novel investigation into the effect of receiver timing jitter on the quality of imaging of ground-penetrating radar systems. We explicitly show the impact of timing jitter throughout an entire example processing pipeline. The process of coupling removal in systems whose receivers have a random Gaussian-distributed jitter with a high standard deviation leaves randomly distributed residue in the radar image. This residue, called jitter-induced coupling noise, extends even into the expected target domain which may cause false alarms. This residue approximately follows a sum of Gaussian and Gamma distributions based on the higher order derivatives of the coupling signal. Once migrated to create an image, the coupling residue is best described by a Weibull distribution. We show that simple filtering does not remove enough residue and propose two preventative design specifications called the linear slope and migrated probability specifications. Antenna and system designers should use these specifications to assess the true severity of coupling signals in a given jitter environment.

47 OTHER INSTRUMENTATION↗

Towards a Global Water Scarcity Risk Assessment Framework: Incorporation of Probability Distributions and Hydro-Climatic Variability

Changing hydro-climatic and socioeconomic conditions increasingly put pressure on fresh water resources and are expected to aggravate water scarcity conditions towards the future. Despite numerous calls for risk-based water scarcity assessments, a global-scale framework that includes UNISDR's definition of risk does not yet exist. This study provides a first step towards such a risk based assessment, applying a Gamma distribution to estimate water scarcity conditions at the global scale under historic and future conditions, using multiple climate change and population growth scenarios. Our study highlights that water scarcity risk, expressed in terms of expected annual exposed population, increases given all future scenarios, up to greater than 56.2% of the global population in 2080. Looking at the drivers of risk, we find that population growth outweigh the impacts of climate change at global and regional scales. Using a risk-based method to assess water scarcity, we show the results to be less sensitive than traditional water scarcity assessments to the use of fixed threshold to represent different levels of water scarcity. This becomes especially important when moving from global to local scales, whereby deviations increase up to 50% of estimated risk levels.

Populations↗

Resolving Ice Cloud Optical Thickness Biases Between CALIOP and MODIS Using Infrared Retrievals

Despite its importance as one of the key radiative properties that determines the impact of upper tropospheric clouds on the radiation balance, ice cloud optical thickness (IOT) has proven to be one of the more challenging properties to retrieve from space-based remote sensing measurements. In particular, optically thin upper tropospheric ice clouds (cirrus) have been especially challenging due to their tenuous nature, extensive spatial scales, and complex particle shapes and light scattering characteristics. The lack of independent validation motivates the investigation presented in this paper, wherein systematic biases between MODIS Collection 5 (C5) and CALIOP Version 3 (V3) unconstrained retrievals of tenuous IOT (< 3) are examined using a month of collocated A-Train observations. An initial comparison revealed a factor of two bias between the MODIS and CALIOP IOT retrievals. This bias is investigated using an infrared (IR) radiative closure approach that compares both products with MODIS IR cirrus retrievals developed for this assessment. The analysis finds that both the MODIS C5 and the unconstrained CALIOP V3 retrievals are biased (high and low, respectively) relative to the IR IOT retrievals. Based on this finding, the MODIS and CALIOP algorithms are investigated with the goal of explaining and minimizing the biases relative to the IR. For MODIS we find that the assumed ice single scattering properties used for the C5 retrievals are not consistent with the mean IR COT distribution. The C5 ice scattering database results in the asymmetry parameter (g) varying as a function of effective radius with mean values that are too large. The MODIS retrievals have been brought into agreement with the IR by adopting a new ice scattering model for Collection 6 (C6) consisting of a modified gamma distribution comprised of a single habit (severely roughened aggregated columns); the C6 ice cloud optical property models have a constant g approx. = 0.75 in the mid-visible spectrum, 5-15% smaller than C5. For CALIOP, the assumed lidar ratio for unconstrained retrievals is fixed at 25 sr for the V3 data products.This value is found to be inconsistent with the constrained (predominantly nighttime) CALIOP retrievals. An experimental data set was produced using a modified lidar ratio of 32 sr for the unconstrained retrievals (an increase of 28%), selected to provide consistency with the constrained V3 results. These modifications greatly improve the agreement with the IR and provide consistency between the MODIS and CALIOP products. Based on these results the recently released MODIS C6 optical products use the single habit distribution given above, while the upcoming CALIOP V4 unconstrained algorithm will use higher lidar ratios for unconstrained retrievals.

Cloud↗

Simulation of radar reflectivity and surface measurements of rainfall

Raindrop size distributions (RSDs) are often estimated using surface raindrop sampling devices (e.g., disdrometers) or optical array (2D-PMS) probes. A number of authors have used these measured distributions to compute certain higher-order RSD moments that correspond to radar reflectivity, attenuation, optical extinction, etc. Scatter plots of these RSD moments versus disdrometer-measured rainrates are then used to deduce physical relationships between radar reflectivity, attenuation, etc., which are measured by independent instruments (e.g., radar), and rainrate. In this paper RSDs of the gamma form as well as radar reflectivity (via time series simulation) are simulated to study the correlation structure of radar estimates versus rainrate as opposed to RSD moment estimates versus rainrate. The parameters N0, D0 and m of a gamma distribution are varied over the range normally found in rainfall, as well as varying the device sampling volume. The simulations are used to explain some possible features related to discrepancies which can arise when radar rainfall measurements are compared with surface or aircraft-based sampling devices.

Chandrasekar, V.↗