Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Testing Nesting DXTRAN Spheres for Neutron Transport in Air

The deterministic transport sphere (DXTRAN sphere) is a technology that can be utilized in MCNP for improving sampling in the vicinity of a particular geometry, when the solid angle, between the source and particular geometry, is small. At a collision event, a DXTRAN particle is created and deterministically transported to the DXTRAN sphere, while the colliding particle (NONDXTRAN particle) continues transporting as if no DXTRAN event had occurred. Weight is balanced by killing the NONDXTRAN particle when traversing the DXTRAN sphere. If set up correctly, the DXTRAN particle placed on the DXTRAN sphere has a high probability of interacting with the geometry within the DXTRAN sphere. The Central Limit Theorem (CLT) states that the estimated mean approaches a normal distribution with known variance as the number of random samples grows sufficiently large. In reality the variance is seldom known and therefore must be estimated. In a Monte Carlo calculation, one must determine when the number of random samples is sufficiently large for the CLT to be satisfied. For the CLT to be satisfied, the first two moments (mean and variance) must exist. The second moment, the variance of a normal distribution is estimated. The slope test, which uses the largest 201 history scores of each tally fluctuation chart (TFC) bin fit to a Pareto function, is used to determine if the largest history scores decrease faster than 1/x{sup 3}. If the slope is less than 3, then enough histories have not yet been sampled such to meet the second moment existence requirement of the CLT. Usually, a slope of less than 3 indicates that the effects of high weight scores have not been satisfactorily sampled; such a consequence can also be manifested in examining trends in the 4. moment, the variance of the variance (VOV). For example, when examining the TFC, if the VOV is monotonically decreasing for a few sets of histories and then drastically increases for the next set of histories, chances are a rare high weight score, significantly greater than the average score, was encountered. As a result of the 1/R{sup 2} term, a point detector can achieve an extremely high weight, much greater than the average, as the distance from collision to point detector approaches zero. To remedy this difficulty, a 'radius of exclusion' is applied, where within the radius of exclusion 'the point detector estimate is assumed to be the average flux uniformly distributed'. A DXTRAN sphere can suffer from the same weight fluctuation as collisions in the vicinity of the DXTRAN sphere produce much larger weights than collisions happening far from the DXTRAN sphere. If the rare large weight put on the DXTRAN sphere then scores in a tally, there will be a significant increase in the VOV and the slope may drop below 3. To compliment that analysis, here we utilize an air composition and density that varies as a function of altitude, detailing the trials and tribulations of attempting to pass the 10 statistical checks for a variety of single and nested DXTRAN sphere approaches. The objective here is to illuminate important considerations for radiation transport problems with small solid angles in a scatter medium of few mean free paths (MFPs)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Formulation of autoconversion and drop spectra shape in shallow cumulus clouds

Two-moment autoconversion parameterizations as compared to accretion parameterizations exhibit significant errors suggesting that additional moments are needed to increase their accuracy. We develop a three-moment autoconversion parameterization using output from an LES model with size-resolved microphysics. Adding the third moment decreases the errors of parameterization and improves precipitation prediction. However, the errors are still significantly larger than errors of accretion rate. An analysis of the cloud drop size distributions (DSDs) in the simulated tropical convective cloud system reveals that most of DSDs have a significant fraction of cloud liquid water content (qc) in the mid-size droplet range (radii from 20 to 40 microns). Our data indicates that more than 30% of DSDs have over half of qc contained in the mid-size range and about 60% of spectra have, at least, one third of qc in this range. Even when the rain/drizzle mode is small (radar reflectivity Z < -10 dBZ), there is a significant number of spectra in which fraction of qc in the mid-size range is as large as 60%. These DSDs are more complex than the frequently used single-mode Gamma or Log-normal distributions, which usually do not have a considerable tail extending to the mid-size range and can be defined by three microphysical moments. The need to define DSDs by more than three moments explains the large errors in the three-moment autoconversion parameterization. The limitation of three-parameter Gamma or Log-normal distributions should be kept in mind when applying them in precipitating shallow Cu clouds.

Kogan, Yefim↗
DOE OSTI2015

SU-E-J-176: Characterization of Inter-Fraction Breast Variability and the Implications On Delivered Dose

Purpose: To systematically characterize inter-fraction breast variability and determine implications on delivered dose. Methods: Weekly port films were used to characterize breast setup variability. Five evenly spaced representative positions across the contour of each breast were chosen on the electronic port film in reference to graticule, and window and level was set such that the skin surface of the breast was visible. Measurements from the skin surface to treatment field edge were taken on each port film at each position and compared to the planning DRR, quantifying the variability. The systematic measurement technique was repeated for all port films for 20 recently treated breast cancer patients. Measured setup variability for each patient was modeled as a normal distribution. The distribution was randomly sampled from the model and applied as isocentric shifts in the treatment planning computer, representing setup variability for each fraction. Dose was calculated for each shifted fraction and summed to obtain DVHs and BEDs that modeled the dose with daily setup variability. Patients were categorized in to relevant groupings that were chosen to investigate the rigorousness of immobilization types, treatment techniques, and inherent anatomical difficulties. Mean position differences and dosimetric differences were evaluated between planned and delivered doses. Results: The setup variability was found to follow a normal distribution with mean position differences between the DRR and port film between − 8.6–3.5 mm with sigma range of 5.3–9.8 mm. Setup position was not found to be significantly different than zero. The mean seroma or whole breast PTV dosimetric difference, calculated as BED, ranged from a −0.23 to +1.13Gy. Conclusion: A systematic technique to quantify and model setup variability was used to calculate the dose in 20 breast cancer patients including variable setup. No statistically significant PTV or OAR BED differences were found between immobilization, treatment or anatomical differences investigated.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Identification of Accretion as Grain Growth Mechanism in Astrophysically Relevant Water–Ice Dusty Plasma Experiment

The grain growth process in the Caltech water–ice dusty plasma experiment has been studied using a high-speed camera and a long-distance microscope lens. It is observed that (i) the ice grain number density decreases fourfold as the average grain major axis increases from 20 to 80 μ m, (ii) the major axis length has a log-normal distribution rather than a power-law dependence, and (iii) no collisions between ice grains are apparent. The grains have a large negative charge resulting in strong mutual repulsion and this, combined with the fractal character of the ice grains, prevents them from agglomerating. In order for the grain kinetic energy to be sufficiently small to prevent collisions between ice grains, the volumetric packing factor (i.e., ratio of the actual volume to the volume of a circumscribing ellipsoid) of the ice grains must be less than ∼0.1 depending on the exact relative velocity of the grains in question. Thus, it is concluded that direct accretion of water molecules is very likely to dominate the observed ice grain growth.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Constraining baryon annihilation in the hadronic phase of heavy-ion collisions via event-by-event fluctuations

We point out that the variance of net-baryon distribution normalized by the Skellam distribution baseline, κ 2 [$B-\overline{B}$]/$\langle B+\overline{B} \rangle$, is sensitive to the possible modification of (anti)baryon yields due to $B\overline{B}$ annihilation in the hadronic phase. The corresponding measurements can thus place stringent limits on the magnitude of the $B\overline{B}$ annihilation and its inverse reaction. We perform Monte Carlo simulations of the hadronic phase in Pb-Pb collisions at the LHC via the recently developed subensemble sampler + UrQMD afterburner and show that the effect survives in net-proton fluctuations, which are directly accessible experimentally. The available experimental data of the ALICE Collaboration on net-proton fluctuations disfavors a notable suppression of (anti)baryon yields in $B\overline{B}$ annihilations predicted by the present version of UrQMD if only global baryon conservation is incorporated. On the other hand, the annihilations improve the data description when local baryon conservation is imposed. The two effects can be disentangled by measuring κ 2 [$B-\overline{B}$]/$\langle B+\overline{B} \rangle$, which at the LHC is notably suppressed by annihilations but virtually unaffected by baryon number conservation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quantifying the unknown impact of segmentation uncertainty on image-based simulations

Abstract Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.

42 ENGINEERING↗
DOE OSTI2015

Ion heating, burnout of the high-frequency field, and ion sound generation under the development of a modulation instability of an intense Langmuir wave in a plasma

The development of one-dimensional parametric instabilities of intense long plasma waves is considered in terms of the so-called hybrid models, with electrons being treated as a fluid and ions being regarded as particles. The analysis is performed for both cases when the average plasma field energy is lower (Zakharov's hybrid model—ZHM) or greater (Silin's hybrid model—SHM) than the plasma thermal energy. The efficiency of energy transfer to ions and to ion perturbations under the development of the instability is considered for various values of electron-to-ion mass ratios. The energy of low-frequency oscillations (ion-sound waves) is found to be much lower than the final ion kinetic energy. We also discuss the influence of the changes in the damping rate of the high-frequency (HF) field on the instability development. The decrease of the absorption of the HF field inhibits the HF field burnout within plasma density cavities and gives rise to the broadening of the HF spectrum. At the same time, the ion velocity distribution tends to the normal distribution in both ZHM and SHM.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Linewidth broadening and tunneling of excitons bound to N pairs in dilute GaAs:N

The exciton bound to a pair of nitrogen atoms situated at nearby lattice sites in dilute GaAs:N provides an energetically uniform electronic system, spectrally distinct from pairs with larger or smaller separations, and can even be grown with a uniform pair orientation in the crystal. We use photoluminescence excitation spectroscopy on an ensemble of N pairs to study the narrow continuous energy distribution within two of the individual exchange- and symmetry-split exciton states. Inhomogeneous linewidths of 50-60 ueV vary across the crystal on a mesoscopic scale and can be 30 ueV at microscopic locations indicating that the homogeneous linewidth inferred from previous time-domain measurements is still considerably broadened. While excitation and emission linewidths are similar, results show a small energy shift between them indicative of exciton transfer via phonon-assisted tunneling between spatially separated N pairs. We numerically simulate the tunneling in a spatial network of randomly distributed pairs having a normal distribution of bound exciton energies. Comparing the ensemble excitation-emission energy shift with the measured results shows that the transfer probability is higher than expected from the dilute pair concentration and what is known of the exciton wavefunction spatial extent. Both the broadening and the exciton transfer have implications for the goal of pair-bound excitons as a single- or multi-qubit system.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TH-CD-BRA-02: 3D Remote Dosimetry for MRI-Guided Radiation Therapy: A Hybrid Approach

Purpose: To validate the dosimetric accuracy of a commercially available MR-IGRT system using a combination of 3D dosimetry measurements (with PRESAGE(R) radiochromic plastic and optical-CT readout) and an in-house developed GPU-accelerated PENELOPE Monte-Carlo dose calculation system. Methods: {sup 60}Co IMRT subject to a 0.35T lateral magnetic field has recently been commissioned in our institution following AAPM’s TG-119 recommendations. We performed PRESAGE(R) sensitivity studies in 4ml cuvettes to verify linearity, MR-compatibility, and energy-independence. Using 10cm diameter PRESAGE(R), we delivered an open calibration field to examine the percent depth dose and a symmetrical 3-field plan with three adjacent regions of varying dose to determine uniformity within the dosimeter under a magnetic field. After initial testing, TG-119 plans were created in the TPS and then delivered to 14.5cm 2kg PRESAGE(R) dosimeters. Dose readout was performed via optical-CT at a second institution specializing in remote 3D dosimetry. Absolute dose was measured using an IBA CC01 ion chamber and the institution standard patient-specific QA methods were used to validate plan delivery. Calculated TG-119 plans were then compared with an independent Monte Carlo dose calculation (gPENELOPE). Results: PRESAGE(R) responds linearly (R{sup 2}=0.9996) to {sup 60}Co irradiation, in the presence of a 0.35T magnetic field, with a sensitivity of 0.0305(±0.003)cm{sup −1}Gy{sup −1}, within 1% of a 6MV non-MR linac irradiation (R{sup 2}=0.9991) with a sensitivity of 0.0302(±0.003)cm{sup −1}Gy{sup −1}. Analysis of TG-119 clinical plans using 3D-gamma (3%/3mm, 10% threshold) give passing rates of: HN 99.1%, prostate 98.0%, C-shape 90.8%, and multi-target 98.5%. The TPS agreed with gPENELOPE with a mean gamma passing rate of 98.4±1.5% (2%/2mm) with the z-score distributions following a standard normal distribution. Conclusion: We demonstrate for the first time that 3D remote dosimetry using both experimental and computational methods is a feasible and reliable approach to commissioning MR-IMRT, which is particularly useful for less specialized clinics in adopting this new treatment modality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Chandra measurements of gas homogeneity and turbulence at intermediate radii in the Perseus Cluster

ABSTRACT We present a Chandra study of surface brightness fluctuations in the diffuse intracluster medium of the Perseus Cluster. Our study utilizes deep, archival imaging of the cluster core as well as a new mosaic of 29 short $\, 5\,$-ks observations extending in eight different directions out to radii of $r_{500} \sim 2.2\, r_{2500}$. Under the assumption that the distribution of densities at a given radius is log-normally distributed, two important quantities can be derived from the width of the log-normal density distribution on a given spatial scale: the density bias, that is equal to the square root of the clumping factor C; and the one-component turbulent velocity, vk,1D. We forward-model all contributions to the measured surface brightness, including astrophysical and particle background components, and account for the Poisson nature of the measured signal. Measuring the distribution of surface brightness fluctuations in 1 arcmin2 regions, spanning the radial range $0.3\rm {-}2.2\, r_{2500}$ ($7.8\rm {-}57.3\,$ arcmin), we find a small to moderate average density bias of around $3{{\ \rm per\ cent}}$ at radii below $1.6\, r_{2500}$. We also infer an average turbulent velocity at these radii of V1D < 400 km s−1. Direct confirmation of our results on turbulent velocities inferred from surface brightness fluctuations should be possible using the X-ray calorimeter spectrometers to be flown aboard the XRISM and Athena observatories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Studies of the 3D surface roughness height

Nowadays nano-coatings occupy more and more significant place in technology. Innovative, functional coatings acquire new aspects from the point of view of modern technologies, considering the aggregate of physical properties that can be achieved manipulating in the production process with the properties of coatings’ surfaces on micro- and nano-level. Nano-coatings are applied on machine parts, friction surfaces, contacting parts, corrosion surfaces, transparent conducting films (TCF), etc. The equipment available at present for the production of transparent conducting oxide (TCO) coatings with highest quality is based on expensive indium tin oxide (ITO) material; therefore cheaper alternatives are being searched for. One such offered alternative is zink oxide (ZnO) nano-coatings. Evaluating the TCF physical and mechanical properties and in view of the new ISO standard (EN ISO 25178) on the introduction of surface texture (3D surface roughness) in the engineering calculations, it is necessary to examine the height of 3D surface roughness, which is one of the most significant roughness parameters. The given paper studies the average values of 3D surface roughness height and the most often applied distribution laws are as follows: the normal distribution and Rayleigh distribution. The 3D surface is simulated by a normal random field.

36 MATERIALS SCIENCE↗
DOE OSTI2014

On bias of kinetic temperature measurements in complex plasmas

The kinetic temperature in complex plasmas is often measured using particle tracking velocimetry. Here, we introduce a criterion which minimizes the probability of faulty tracking of particles with normally distributed random displacements in consecutive frames. Faulty particle tracking results in a measurement bias of the deduced velocity distribution function and hence the deduced kinetic temperature. For particles with a normal velocity distribution function, mistracking biases the obtained velocity distribution function towards small velocities at the expense of large velocities, i.e., the inferred velocity distribution is more peaked and its tail is less pronounced. The kinetic temperature is therefore systematically underestimated in measurements. We give a prescription to mitigate this type of error.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Composited structure of non-precipitating shallow cumulus clouds

The normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large-eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water ( q l ) decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (θ v ), and potential temperature (θ) perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, θ' v becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power-law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud-top downdrafts cannot be captured well.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

SMALL-SCALE AND GLOBAL DYNAMOS AND THE AREA AND FLUX DISTRIBUTIONS OF ACTIVE REGIONS, SUNSPOT GROUPS, AND SUNSPOTS: A MULTI-DATABASE STUDY

In this work, we take advantage of 11 different sunspot group, sunspot, and active region databases to characterize the area and flux distributions of photospheric magnetic structures. We find that, when taken separately, different databases are better fitted by different distributions (as has been reported previously in the literature). However, we find that all our databases can be reconciled by the simple application of a proportionality constant, and that, in reality, different databases are sampling different parts of a composite distribution. This composite distribution is made up by linear combination of Weibull and log-normal distributions—where a pure Weibull (log-normal) characterizes the distribution of structures with fluxes below (above) 10{sup 21}Mx (10{sup 22}Mx). Additionally, we demonstrate that the Weibull distribution shows the expected linear behavior of a power-law distribution (when extended to smaller fluxes), making our results compatible with the results of Parnell et al. We propose that this is evidence of two separate mechanisms giving rise to visible structures on the photosphere: one directly connected to the global component of the dynamo (and the generation of bipolar active regions), and the other with the small-scale component of the dynamo (and the fragmentation of magnetic structures due to their interaction with turbulent convection)

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Analytical Methods of Approximating the Fisher Information Matrix for the Lognormal Distribution

Modelling the matter density in the universe through time has become one of the most important open problems in cosmology. One distribution that works particularly well is the lognormal distribution, where if Y follows a lognormal distribution, ln Y follows a normal distribution. In many cases, we want to know how much information a variable carries about some unknown parameter; this can be accomplished by calculating the Fisher information matrix. However, for the lognormal distribution, the Fisher matrix is much more complicated to evaluate in closed form. In this paper, we will look at some analytical methods to evaluate the Fisher matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Effect Of Ionization Variance On Nuclear-recoil Dark Matter Searches

Direct detection dark matter experiments are increasingly interested in the low-mass dark matter parameter space, but zero-background low-mass searches require event separation between the electron and nuclear recoil bands, which requires a proper understanding of detector energy reconstruction. Previous simulations have shown that we do not entirely understand the ionization efficiency (yield) for electron and nuclear recoils, as the assumption that the distribution for the yield is normally distributed for a true recoil energy is violated. Since the yield distribution may directly affect dark matter low-mass limits, it is crucial we understand how the yield is distributed. A component to understanding the yield distribution is the variance in the number of electron-hole pairs produced or ionization variance. This effect has been studied relatively infrequently as experiments have been interested in large energy deposits (10 - 100 keV) and could accurately separate electron and nuclear recoil events. For electron recoils, the ionization variance is described by a “Fano” factor. For nuclear recoils the effect can be parameterized by an “effective” Fano factor, which has similar definition but a different physical origin. The nuclear recoil “effective” Fano factor is shown to be much larger than the electron-recoil version above around 10 keV deposited energy.

Matheny, Mitchell Douglas↗
DOE OSTI2015

STOCHASTIC MODEL OF THE SPIN DISTRIBUTION OF DARK MATTER HALOS

We employ a stochastic approach to probing the origin of the log-normal distributions of halo spin in N-body simulations. After analyzing spin evolution in halo merging trees, it was found that a spin change can be characterized by a stochastic random walk of angular momentum. Also, spin distributions generated by random walks are fairly consistent with those directly obtained from N-body simulations. We derived a stochastic differential equation from a widely used spin definition and measured the probability distributions of the derived angular momentum change from a massive set of halo merging trees. The roles of major merging and accretion are also statistically analyzed in evolving spin distributions. Several factors (local environment, halo mass, merging mass ratio, and redshift) are found to influence the angular momentum change. The spin distributions generated in the mean-field or void regions tend to shift slightly to a higher spin value compared with simulated spin distributions, which seems to be caused by the correlated random walks. We verified the assumption of randomness in the angular momentum change observed in the N-body simulation and detected several degrees of correlation between walks, which may provide a clue for the discrepancies between the simulated and generated spin distributions in the voids. However, the generated spin distributions in the group and cluster regions successfully match the simulated spin distribution. We also demonstrated that the log-normality of the spin distribution is a natural consequence of the stochastic differential equation of the halo spin, which is well described by the Geometric Brownian Motion model.

79 ASTRONOMY AND ASTROPHYSICS↗