Search NASA⌕ Search

SEARCH · Search NASA

Results for “normal 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 613 records · Page 34

Working with Bézier Curves as bases for Functional Expansion Tallies

Functional expansion tallies (FETs) are powerful tools for getting more information per history from Monte Carlo simulations, but in the past they have been constrained to orthogonal bases. Bézier curves are used widely in computer aided design (CAD) geometry kernels and could be well-suited for FETs due to their ability to assume many arbitrary shapes, but they use nonorthogonal bases. Recent developments in 2021 have made nonorthogonal FETs possible. The convergence of Bézier curve FETs in both polynomial order and number of samples is explored in this work. It is shown that these bases are well-suited for representing normal distributions and this opens the door to the possibility of other CAD-derived FET bases.

97 - MATHEMATICS AND COMPUTING↗

FREQUENCY DISTRIBUTIONS OF 90SR AND 137CS CONCENTRATIONS IN AN ECOSYSTEM OF THE 'RED FOREST' AREA IN THE CHERNOBYL EXCLUSION ZONE

In the most highly contaminated region of the Chernobyl Exclusion Zone: the 'Red Forest' site, the accumulation of the major dose-affecting radionuclides ({sup 90}Sr and {sup 137}Cs) within the components of an ecological system encompassing 3,000 m{sup 2} were characterized. The sampled components included soils (top 0-10 cm depth), Molina caerulea (blue moor grass), Camponotus vagus (carpenter ants) and Pelobates fuscus (spade-footed toad). In a comparison among the components of this ecosystem, the {sup 90}Sr and {sup 137}Cs concentrations measured in 40 separate grids exhibited significant differences, while the frequency distribution of the values were close to a logarithmically normal leptokurtic distribution with a significant right-side skew. While it is important to identify localized areas of high contamination or 'hot spots,' including these values in the arithmetic mean may overestimate the exposure risk. In component sample sets that exhibited logarithmically normal distribution, the geometrical mean more accurately characterizes a site. Ideally, risk assessment is most confidently achieved when the arithmetic and geometrical means are most similar, meaning the distribution approaches normal. Through bioaccumulation, the highest concentrations of {sup 90}Sr and {sup 137}Cs were measured in the blue moor grass and spade-footed toad. These components also possessed distribution parameters that shifted toward a normal distribution.

54 ENVIRONMENTAL SCIENCES↗

Retrieval of Normalized Gamma Size Distribution Parameters Using Precipitation Imaging Package (Pip) Snowfall Observations During Ice-Pop 2018

Parameters of the normalized gamma particle size distribution (PSD) have been retrieved from the Precipitation Image Package (PIP) snowfall observations collected during the International Collaborative Experiment - PyeongChang Olympics and Paralympic (ICE-POP 2018). Two of the gamma PSD parameters, the mass weighted particle diameter (Dmass) and the normalized intercept parameter NW, have median values of 1.15-1.31 mm and 2.84-3.04 log(mm-1 m-3), respectively. This range arises from the choice of the relationship between the maximum versus equivalent diameter, Dmx-Deq, and the relationship between the Reynolds and Best numbers, Re-X. Normalization of snow water equivalent rate (SWER) and ice water content (W) by NW reduces the range in NW resulting in well fitted power law relationship, between SWER/NW and Dmass and between W/NW and Dmass. The bulk descriptors of snowfall are calculated from PIP observations and from the gamma PSD with values of the shape parameter (μ) ranging from -2 to 10. NASA’s Global Precipitation Measurement (GPM) mission, which adopted the normalized gamma PSD, assumes μ = 2 and μ = 3 in its two separate algorithms. The mean fractional bias (MFB) of the snowfall parameters changes with μ, where the functional dependence on μ depends on the specific snowfall parameter of interest. The MFB of the total concentration was underestimated by 0.23-0.34 when μ = 2 and by 0.29-0.40 when μ = 3, while the MFB of SWER had a much narrower range (-0.03 to 0.04) for the same μ values.

Snowfall↗

Study of the normalized transverse momentum distribution of W bosons produced in $p\bar{p}$ collisions at $\sqrt{s}$ = 1.96 TeV

We present a study of the normalized transverse momentum distribution of W bosons produced in $p\bar{p}$ collisions, using data corresponding to an integrated luminosity of 4.35 fb -1 collected with the D0 detector at the Fermilab Tevatron collider at $\sqrt{s}$ = 1.96 TeV. The measurement focuses on the transverse momentum region below 15 GeV, which is of special interest for electroweak precision measurements; it relies on the same detector calibration methods which were used for the precision measurement of the W boson mass. The measured distribution is compared to different QCD predictions and a procedure is given to allow the comparison of any further theoretical models to the D0 data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

XPCA - R Implementation

Principal component analysis (PCA) is arguably the most popular tool in multivariate exploratory data analysis. We consider the question of how to handle heterogeneous variables that include continuous, binary, and ordinal. In the probabilistic interpretation of low-rank PCA, the data has a normal multivariate distribution and, therefore, normal marginal distributions for each column. If some marginals are continuous but not normal, the semiparametric copula-based principal component analysis (COCA) method is an alternative to PCA that combines a Gaussian copula with nonparametric marginals. If some marginals are discrete or semi-continuous, we present a new extended PCA (XPCA) method that also uses a Gaussian copula and nonparametric marginals and accounts for discrete variables in the likelihood calculation by integrating over appropriate intervals. Like PCA, the factors produced by XPCA can be used to find latent structure in data, build predictive models, and perform dimensionality reduction. We implement the new model, XPCA, and its ancestors PCA and COCA in both R and Python. This software provides implementations of XPCA, COCA, and probabilistic PCA in both R and Python. XPCA is written in arXiv report https://arxiv.org/abs/1808.07510. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525. SAND2019-10591 M

Kincher-Winoto, Kina↗

XPCA -Python Implementation

Principal component analysis (PCA) is arguably the most popular tool in multivariate exploratory data analysis. We consider the question of how to handle heterogeneous variables that include continuous, binary, and ordinal. In the probabilistic interpretation of low-rank PCA, the data has a normal multivariate distribution and, therefore, normal marginal distributions for each column. If some marginals are continuous but not normal, the semiparametric copula-based principal component analysis (COCA) method is an alternative to PCA that combines a Gaussian copula with nonparametric marginals. If some marginals are discrete or semi-continuous, we present a new extended PCA (XPCA) method that also uses a Gaussian copula and nonparametric marginals and accounts for discrete variables in the likelihood calculation by integrating over appropriate intervals. Like PCA, the factors produced by XPCA can be used to find latent structure in data, build predictive models, and perform dimensionality reduction. We implement the new model, XPCA, and its ancestors PCA and COCA in both R and Python. This software provides implementations of XPCA, COCA, and probabilistic PCA in both R and Python. XPCA is written in arXiv report https://arxiv.org/abs/1808.07510. SAND2019-10591 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Kincher-Winoto, Kina↗

MO-G-304-02: Knowledge Based DVH Prediction Using a Geometric Dose Transform

Purpose: To demonstrate a novel method for predicting patient dose-volume histograms (DVHs) using a prior database of optimized radiotherapy treatment plans. Such predicted DVHs could be useful for automating treatment planning. Methods: Our initial demonstration utilized a database of 100 prostate intensity-modulated radiotherapy (IMRT) data-sets. Each data-set contained a CT image with contours of the planning target volume (PTV), rectum, and bladder, the parameters of a clinically approved IMRT plan, and a corresponding simulated dose distribution. We applied a novel geometric transformation to remove the influence of the PTV size, shape, and location on the dose distribution. We termed the transformed distribution the geometrically normalized dose distribution (GNDD). This normalization transform was applied to 80 data-sets randomly selected from the database, and a population GNDD was computed as the average. Next, the population GNDD was mapped onto each of the remaining 20 patient datasets using the reverse of the geometric normalization transform, and predicted DVHs were calculated from the reverse transformed dose distributions (GNDD-DVHs). In addition, a state of the art machine learning based method from the literature was tested for comparison. Results: DVH prediction accuracy was quantified by calculating the relative root mean squared error (rRMSE) on predicted DVHs for the 20 test patients using their known DVHs. For bladder, rectum, and PTV average rRMSEs for the GNDD method were 9.7 ± 4.2%, 13.9 ± 6.0%, and 2.3 ± 0.5% respectively. Prediction results using GNDD were roughly equivalent to that from the machine learning method. Conclusion: We developed a new method for predicting DVH curves from a database of prior patient plans. We demonstrated that our simple approach achieves accuracy comparable to a method using a complicated machine learning based approach.

60 APPLIED LIFE SCIENCES↗

A General N -Moment Normalization Method for Deriving Raindrop Size Distribution Scaling Relationships

A general drop size distribution (DSD) normalization method is formulated in terms of generalized power series relating any DSD moment to any number and combination of reference moments. This provides a consistent framework for comparing the variability of normalized DSD moments using different sets of reference moments, with no explicit assumptions about the DSD functional form (e.g., gamma). It also provides a method to derive any unknown moment plus an estimate of its uncertainty from one or more known moments, which is relevant to remote sensing retrievals and bulk microphysics schemes in weather and climate models. Herein, the approach is applied to a large dataset of disdrometer-observed and bin microphysics-modeled DSDs. As expected, the spread of normalized moments decreases as the number of reference moments is increased, quantified by the logarithmic standard deviation of the normalized moments, σ. Averaging σ for all combinations of reference moments and normalized moments of integer order 0–10, 42.9%, 81.3%, 93.7%, and 96.9% of spread are accounted for applying one-, two-, three-, and four-moment normalizations, respectively. Thus, DSDs can be well characterized overall using three reference moments, whereas adding a fourth reference moment contributes little independent information. The spread of disdrometer-observed DSD moments from uncertainty associated with drop count statistics generally lies between values of σ using two- and three-moment normalizations. However, this uncertainty has little impact on the derived DSD scaling relationships or σ when considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the distribution of sunspot cycle periods

A comparison is made between the observed distribution of sunspot cycle periods and distributions based on uniform, normal, and bimodal distributions. The bimodal distribution, composed of short-period and long-period cycles, is found to best describe the observed distribution. Compared to the normal distribution for the most reliably determined cycles (cycles 8-20), the bimodal distribution has a residual (sum of squares of differences) that is about 86 percent smaller. Means for short-period and long-period cycles are estimated to be 122 + or - 4 months and 140 + or - 5 months, respectively.

Wilson, Robert M.↗

A New Bond Albedo for Performing Orbital Debris Brightness to Size Transformations

We have developed a technique for estimating the intrinsic size distribution of orbital debris objects via optical measurements alone. The process is predicated on the empirically observed power-law size distribution of debris (as indicated by radar RCS measurements) and the log-normal probability distribution of optical albedos as ascertained from phase (Lambertian) and range-corrected telescopic brightness measurements. Since the observed distribution of optical brightness is the product integral of the size distribution of the parent [debris] population with the albedo probability distribution, it is a straightforward matter to transform a given distribution of optical brightness back to a size distribution by the appropriate choice of a single albedo value. This is true because the integration of a powerlaw with a log-normal distribution (Fredholm Integral of the First Kind) yields a Gaussian-blurred power-law distribution with identical power-law exponent. Application of a single albedo to this distribution recovers a simple power-law [in size] which is linearly offset from the original distribution by a constant whose value depends on the choice of the albedo. Significantly, there exists a unique Bond albedo which, when applied to an observed brightness distribution, yields zero offset and therefore recovers the original size distribution. For physically realistic powerlaws of negative slope, the proper choice of albedo recovers the parent size distribution by compensating for the observational bias caused by the large number of small objects that appear anomalously large (bright) - and thereby skew the small population upward by rising above the detection threshold - and the lower number of large objects that appear anomalously small (dim). Based on this comprehensive analysis, a value of 0.13 should be applied to all orbital debris albedo-based brightness-to-size transformations regardless of data source. Its prima fascia genesis, derived and constructed from the current RCS to size conversion methodology (SiBAM Size-Based Estimation Model) and optical data reduction standards, assures consistency in application with the prior canonical value of 0.1. Herein we present the empirical and mathematical arguments for this approach and by example apply it to a comprehensive set of photometric data acquired via NASA's Liquid Mirror Telescopes during the 2000-2001 observing season.

Mulrooney, Mark K.↗

Rapid flattening of butterfly pitch angle distributions of radiation belt electrons by whistler-mode chorus

Van Allen radiation belt electrons exhibit complex dynamics during geomagnetically active periods. Investigation of electron pitch angle distributions (PADs) can provide important information on the dominant physical mechanisms controlling radiation belt behaviors. In this paper, we report a storm time radiation belt event where energetic electron PADs changed from butterfly distributions to normal or flattop distributions within several hours. Van Allen Probes observations showed that the flattening of butterfly PADs was closely related to the occurrence of whistler-mode chorus waves. Two-dimensional quasi-linear STEERB simulations demonstrate that the observed chorus can resonantly accelerate the near-equatorially trapped electrons and rapidly flatten the corresponding electron butterfly PADs. Finally, these results provide a new insight on how chorus waves affect the dynamic evolution of radiation belt electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

Gauges, loops, and polynomials for partition functions of graphical models

Graphical models represent multivariate and generally not normalized probability distributions. Computing the normalization factor, called the partition function, is the main inference challenge relevant to multiple statistical and optimization applications. The problem is #P-hard that is of an exponential complexity with respect to the number of variables. Here, aimed at approximating the partition function, we consider multi-graph models where binary variables and multivariable factors are associated with edges and nodes, respectively, of an undirected multi-graph. We suggest a new methodology for analysis and computations that combines the Gauge function technique from Chertkov and Chernyak with the technique developed in Anari and Oveis Gharan 2017 arXiv:1702.02937; Gurvits 2011 arXiv:1106.2844; Straszak and Vishnoi 2017 55th Annual Allerton Conf. on Communication, Control, and Computing, based on the recent progress in the field of real stable polynomials. We show that the Gauge function, representing a single-out term in a finite sum expression for the partition function which achieves extremum at the so-called belief-propagation gauge, has a natural polynomial representation in terms of gauges/variables associated with edges of the multi-graph. Moreover, Gauge function can be used to recover the partition function through a sequence of transformations allowing appealing algebraic and graphical interpretations. Algebraically, one step in the sequence consists of the application of a differential operator over gauges associated with an edge. Graphically, the sequence is interpreted as a repetitive elimination/contraction of edges resulting in multi-graph models on decreasing in size (number of edges) graphs with the same partition function as in the original multi-graph model. Even though the complexity of computing factors in the sequence of the derived multi-graph models and respective Gauge functions grow exponentially with the number of eliminated edges, polynomials associated with the new factors remain bi-stable if the original factors have this property. Moreover, we show that BP estimations in the sequence do not decrease, each low-bounding the partition function.

97 MATHEMATICS AND COMPUTING↗

Quantitative model of magnetic coupling between solar wind and magnetosphere

Preliminary results are presented of a quantitative three-dimensional model of an open steady-state magnetosphere configuration incorporating a normal-component distribution corresponding to the subsolar merging-line hypothesis. The distribution of the normal magnetic-field component at the magnetopause is used as input and is used to calculate an interconnection magnetic field that links the internal and external fields. The interconnected field is then used to map the solar-wind electric field onto the polar cap. The resulting polar-cap flow patterns are found to be in agreement with observations.

Toffoletto, F. R.↗

Automated quantification of mitral valve regurgitation based on normalized centerline velocity distribution

Previous echocardiographic techniques for quantifying valvular regurgitation are limited by factors including uncertainties for orifice location and a hemispheric convergence assumption that often results in over- and underestimation of flow rate and regurgitant orifice area. Using computational fluid dynamics simulations, these factors were eliminated, allowing a more accurate assessment of regurgitation. A model was developed to allow automated quantification of regurgitant orifice diameter based on the centerline velocity data available from color M-mode echocardiography. The model, validated using in vitro unsteady flow data, demonstrated improved accuracy for orifice diameter (y=0.95x + 0.38, r=0.96) and volume (y=1.18x - 4.72, r=0.93).

NASA Program Biomedical Research and Countermeasur↗

A simplified method for the calculation of airfoil pressure distribution

A method is presented for the rapid calculation of the pressure distribution over an airfoil section when the normal-force distribution and the pressure distribution over the "base profile" (i.e., the profile of the same airfoil were the camber line straight and the resulting airfoil at zero angle of attack) are known. This note is intended as a supplement to N.A.C.A. Report Nos. 631 and 634 wherein methods are presented for the calculation of the normal-force distribution over plain and flapped airfoils, respectively, but not of the pressures on the individual surfaces. Base-profile pressure-coefficient distributions for the usual N.A.C.A. family of airfoils, which are also suitable for several other commonly employed airfoils, are included in tabular form. With these tabulated base-profile pressures and the computed normal-force distributions, pressure distributions adequate for most engineering purposes can be obtained.

Allen, H Julian↗