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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Some functional limit theorems for compound Cox processes

An improved version of the functional limit theorem is proved establishing weak convergence of random walks generated by compound doubly stochastic Poisson processes (compound Cox processes) to Lévy processes in the Skorokhod space under more realistic moment conditions. As corollaries, theorems are proved on convergence of random walks with jumps having finite variances to Lévy processes with variance-mean mixed normal distributions, in particular, to stable Lévy processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Electric field stimulated growth of Zn whiskers

We have investigated the impact of strong (∼10{sup 4} V/cm) electric fields on the development of Zn whiskers. The original samples, with considerable whisker infestation were cut from Zn-coated steel floors and then exposed to electric fields stresses for 10-20 hours at room temperature. We used various electric field sources, from charges accumulated in samples irradiated by: (1) the electron beam of a scanning electron microscope (SEM), (2) the electron beam of a medical linear accelerator, and (3) the ion beam of a linear accelerator; we also used (4) the electric field produced by a Van der Graaf generator. In all cases, the exposed samples exhibited a considerable (tens of percent) increase in whiskers concentration compared to the control sample. The acceleration factor defined as the ratio of the measured whisker growth rate over that in zero field, was estimated to approach several hundred. The statistics of lengths of e-beam induced whiskers was found to follow the log-normal distribution known previously for metal whiskers. The observed accelerated whisker growth is attributed to electrostatic effects. These results offer promise for establishing whisker-related accelerated life testing protocols.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Detection of vapor nanobubbles by small angle neutron scattering (SANS)

Experiments using boiling water on untreated (roughness 100–300 nm) metal surfaces using small-angle neutron scattering (SANS) show the appearance of structures that are 50–70 nm in size when boiling is present. The scattering signal disappears when the boiling ceases, and no change in the signal is detected at any surface temperature condition below saturation. This confirms that the signal is caused by vapor nanobubbles. Two boiling regimes are evaluated herein that differ by the degree of subcooling (3–10 °C). A polydisperse spherical model with a log-normal distribution fits the SANS data well. In conclusion, the size distribution indicates that a large number of nanobubbles exist on the surface during boiling, and some of them grow into large bubbles.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro–Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density–temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Chemistry↗
DOE OSTI2015

Spectral dependences of extrinsic optical absorption in sillenite crystals

The influence of laser irradiation at wavelengths of 532 and 655 nm and annealing in air at temperatures from 200 to 370 °C on optical absorption spectra of undoped bismuth silicon oxide and bismuth germanium oxide and aluminium-doped bismuth titanium oxide crystals has been studied experimentally. The experimental data have been interpreted in terms of a model for extrinsic absorption that takes into account not only the contribution of the photoexcitation of electrons from deep donor centres with a normal distribution of their concentration with respect to ionisation energy but also that of intracentre transitions. (laser applications and other topics in quantum electronics)

36 MATERIALS SCIENCE↗
DOE OSTI2015

Enhanced enstrophy generation for turbulent convection in low-Prandtl-number fluids

Turbulent convection is often present in liquids with a kinematic viscosity much smaller than the diffusivity of the temperature. Here we reveal why these convection flows obey a much stronger level of fluid turbulence than those in which kinematic viscosity and thermal diffusivity are the same; i.e., the Prandtl number Pr is unity. We compare turbulent convection in air at Pr = 0.7 and in liquid mercury at Pr = 0.021. In this comparison the Prandtl number at constant Grashof number Gr is varied, rather than at constant Rayleigh number Ra as usually done. Our simulations demonstrate that the turbulent Kolmogorov-like cascade is extended both at the large- and small-scale ends with decreasing Pr. The kinetic energy injection into the flow takes place over the whole cascade range. In contrast to convection in air, the kinetic energy injection rate is particularly enhanced for liquid mercury for all scales larger than the characteristic width of thermal plumes. As a consequence, mean values and fluctuations of the local strain rates are increased, which in turn results in significantly enhanced enstrophy production by vortex stretching. The normalized distributions of enstrophy production in the bulk and the ratio of the principal strain rates are found to agree for both Prs. Finally, despite the different energy injection mechanisms, the principal strain rates also agree with those in homogeneous isotropic turbulence conducted at the same Reynolds numbers as for the convection flows. Thus, our results have interesting implications for small-scale turbulence modeling of liquid metal convection in astrophysical and technological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Welch Method and Bootstrapping Applied to Subcritical Gamma Noise

We measured the prompt neutron decay constant 𝛼 of the CROCUS zero-power reactor at the Swiss Federal Institute of Technology Lausanne using cross-power spectral density (CPSD) analysis of gamma-gamma correlations from two trans-stilbene organic scintillators positioned near the reactor core. We measured critical and subcritical states, with water levels ranging from 960 mm (critical) to 800 mm (𝜌=−1.4 $ subcritical). Our analysis used the Welch method, dividing signal segments for fast Fourier transform (FFT) frequency analysis and applying bootstrapping uncertainty quantification that uses Welch-defined segments. Results demonstrated a clear increase in the measured 𝛼 as reactor reactivity decreased, distinguishing critical from subcritical conditions. At the 960-mm critical level, 𝛼 was estimated at 155.9 ± 0.7 s −1 , and for the 800-mm subcritical level, 𝛼 increased significantly to 367.3 ± 6.9 s –1 . A linear regression of subcritical states yielded a critical estimate of 154.0 ± 3.1 s –1 , aligning with the static 𝛼 estimate at critical. The bootstrapping method produced normally distributed 𝛼 estimates, confirming data consistency. The gamma CPSD 𝛼 estimates clearly distinguish reactor states and improve monitoring of zero-power reactors. The future deployment of modular and microreactors as potential candidates for noise analysis is demonstrated in CROCUS, particularly zero-power mock-ups of new designs. The improvement of noise analysis in the subcritical domain from this work will support experimental data for reactor deployment and procedure.

CROCUS↗
DOE OSTI2016

Spatiotemporal modeling of node temperatures in supercomputers

Los Alamos National Laboratory (LANL) is home to many large supercomputing clusters. These clusters require an enormous amount of power (~500-2000 kW each), and most of this energy is converted into heat. Thus, cooling the components of the supercomputer becomes a critical and expensive endeavor. Recently a project was initiated to investigate the effect that changes to the cooling system in a machine room had on three large machines that were housed there. Coupled with this goal was the aim to develop a general good-practice for characterizing the effect of cooling changes and monitoring machine node temperatures in this and other machine rooms. This paper focuses on the statistical approach used to quantify the effect that several cooling changes to the room had on the temperatures of the individual nodes of the computers. The largest cluster in the room has 1,600 nodes that run a variety of jobs during general use. Since extremes temperatures are important, a Normal distribution plus generalized Pareto distribution for the upper tail is used to model the marginal distribution, along with a Gaussian process copula to account for spatio-temporal dependence. A Gaussian Markov random field (GMRF) model is used to model the spatial effects on the node temperatures as the cooling changes take place. This model is then used to assess the condition of the node temperatures after each change to the room. The analysis approach was used to uncover the cause of a problematic episode of overheating nodes on one of the supercomputing clusters. Lastly, this same approach can easily be applied to monitor and investigate cooling systems at other data centers, as well.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2010

ON THE ESTIMATE OF MAGNETIC NON-POTENTIALITY OF SUNSPOTS DERIVED USING HINODE SOT/SP OBSERVATIONS: EFFECT OF POLARIMETRIC NOISE

The accuracy of Milne-Eddington (ME) inversions, used to retrieve the magnetic field vector, depends upon the signal-to-noise ratio (S/N) of the spectro-polarimetric observations. The S/N in real observations varies from pixel to pixel; therefore the accuracy of the field vector also varies over the map. The aim of this work is to study the effect of polarimetric noise on the inference of the magnetic field vector and the magnetic non-potentiality of a real sunspot. To this end, we use the Hinode SOT/SP vector magnetogram of a real sunspot NOAA 10933 as an input to generate synthetic Stokes profiles under ME model assumptions. We then add normally distributed polarimetric noise of the level 0.5% of continuum intensity to these synthetic profiles and invert them again using the ME code. This process is repeated 100 times with different realizations of noise. It is found that within most of the sunspot areas (>90% area) the spread in the (1) field strength is less than 8 G, (2) field inclination is less than 1{sup 0}, and (3) field azimuth is less than 5{sup 0}. Further, we determine the uncertainty in the magnetic non-potentiality of a sunspot as determined by the force-free parameter {alpha}{sub g} and spatially averaged signed shear angle (SASSA). It is found that for the sunspot studied here these parameters are {alpha}{sub g} = -3.5 {+-} 0.37(x10{sup -9} m{sup -1}) and SASSA = -1.68 {+-} 0.{sup 0}014. This suggests that the SASSA is a less dispersed non-potentiality parameter as compared to {alpha}{sub g}. Further, we examine the effect of increasing noise levels, viz. 0.01%, 0.1%, 0.5%, and 1% of continuum intensity, and find that SASSA is less vulnerable to noise as compared to the {alpha}{sub g} parameter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

MAXIMUM LIKELIHOOD FITTING OF X-RAY POWER DENSITY SPECTRA: APPLICATION TO HIGH-FREQUENCY QUASI-PERIODIC OSCILLATIONS FROM THE NEUTRON STAR X-RAY BINARY 4U1608-522

High-frequency quasi-periodic oscillations (QPOs) from weakly magnetized neutron stars display rapid frequency variability (second timescales) and high coherence with quality factors up to at least 200 at frequencies about 800-850 Hz. Their parameters have been estimated so far from standard min({chi}{sup 2}) fitting techniques, after combining a large number of power density spectra (PDS), to have the powers normally distributed (the so-called Gaussian regime). Before combining PDS, different methods to minimize the effects of the frequency drift to the estimates of the QPO parameters have been proposed, but none of them relied on fitting the individual PDS. Accounting for the statistical properties of PDS, we apply a maximum likelihood method to derive the QPO parameters in the non-Gaussian regime. The method presented is general, easy to implement, and can be applied to fitting individual PDS, several PDS simultaneously, or their average, and is obviously not specific to the analysis of kHz QPO data. It applies to the analysis of any PDS optimized in frequency resolution and for low-frequency variability or PDS containing features whose parameters vary on short timescales, as is the case for kHz QPOs. It is equivalent to the standard {chi}{sup 2} minimization fitting when the number of PDS fitted is large. The accuracy, reliability, and superiority of the method is demonstrated with simulations of synthetic PDS, containing Lorentzian QPOs of known parameters. Accounting for the broadening of the QPO profile, due to the leakage of power inherent to windowed Fourier transforms, the maximum likelihood estimates of the QPO parameters are asymptotically unbiased and have negligible bias when the QPO is reasonably well detected. By contrast, we show that the standard min({chi}{sup 2}) fitting method gives biased parameters with larger uncertainties. The maximum likelihood fitting method is applied to a subset of archival Rossi X-ray Timing Explorer data of the neutron star X-ray binary 4U1608-522, for which we show that the lower kHz QPO parameters can be measured on timescales as short as 8 s. To demonstrate the potential use of the results of the maximum likelihood method, we show that in the observation analyzed the time evolution of the frequency is consistent with a random walk. We then show that the broadening of the QPO due to the frequency drift scales as {radical}T, as expected from a random walk (T is the integration time of the PDS). This enables us to estimate the intrinsic quality factor of the QPO to be {approx}260, whereas previous analysis indicated a maximum value around 200.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

MASS GROWTH AND MERGERS: DIRECT OBSERVATIONS OF THE LUMINOSITY FUNCTION OF LRG SATELLITE GALAXIES OUT TO z = 0.7 FROM SDSS AND BOSS IMAGES

We present a statistical study of the luminosity functions of galaxies surrounding luminous red galaxies (LRGs) at average redshifts (z) = 0.34 and (z) = 0.65. The luminosity functions are derived by extracting source photometry around more than 40,000 LRGs and subtracting foreground and background contamination using randomly selected control fields. We show that at both studied redshifts the average luminosity functions of the LRGs and their satellite galaxies are poorly fitted by a Schechter function due to a luminosity gap between the centrals and their most luminous satellites. We utilize a two-component fit of a Schechter function plus a log-normal distribution to demonstrate that LRGs are typically brighter than their most luminous satellite by roughly 1.3 mag. This luminosity gap implies that interactions within LRG environments are typically restricted to minor mergers with mass ratios of 1:4 or lower. The luminosity functions further imply that roughly 35% of the mass in the environment is locked in the LRG itself, supporting the idea that mass growth through major mergers within the environment is unlikely. Lastly, we show that the luminosity gap may be at least partially explained by the selection of LRGs as the gap can be reproduced by sparsely sampling a Schechter function. In that case LRGs may represent only a small fraction of central galaxies in similar mass halos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE STAR FORMATION RATE OF TURBULENT MAGNETIZED CLOUDS: COMPARING THEORY, SIMULATIONS, AND OBSERVATIONS

The role of turbulence and magnetic fields is studied for star formation in molecular clouds. We derive and compare six theoretical models for the star formation rate (SFR)-the Krumholz and McKee (KM), Padoan and Nordlund (PN), and Hennebelle and Chabrier (HC) models, and three multi-freefall versions of these, suggested by HC-all based on integrals over the log-normal distribution of turbulent gas. We extend all theories to include magnetic fields and show that the SFR depends on four basic parameters: (1) virial parameter {alpha}{sub vir}; (2) sonic Mach number M; (3) turbulent forcing parameter b, which is a measure for the fraction of energy driven in compressive modes; and (4) plasma {beta}=2M{sub A}{sup 2}/M{sup 2} with the Alfven Mach number M{sub A}. We compare all six theories with MHD simulations, covering cloud masses of 300 to 4 Multiplication-Sign 10{sup 6} M{sub Sun} and Mach numbers M=3-50 and M{sub A}=1-{infinity}, with solenoidal (b = 1/3), mixed (b = 0.4), and compressive turbulent (b = 1) forcings. We find that the SFR increases by a factor of four between M=5 and 50 for compressive turbulent forcing and {alpha}{sub vir} {approx} 1. Comparing forcing parameters, we see that the SFR is more than 10 times higher with compressive than solenoidal forcing for M=10 simulations. The SFR and fragmentation are both reduced by a factor of two in strongly magnetized, trans-Alfvenic turbulence compared to hydrodynamic turbulence. All simulations are fit simultaneously by the multi-freefall KM and multi-freefall PN theories within a factor of two over two orders of magnitude in SFR. The simulated SFRs cover the range and correlation of SFR column density with gas column density observed in Galactic clouds, and agree well for star formation efficiencies SFE = 1%-10% and local efficiencies {epsilon} = 0.3-0.7 due to feedback. We conclude that the SFR is primarily controlled by interstellar turbulence, with a secondary effect coming from magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE COLLISIONAL EVOLUTION OF DEBRIS DISKS

We explore the collisional decay of disk mass and infrared emission in debris disks. With models, we show that the rate of the decay varies throughout the evolution of the disks, increasing its rate up to a certain point, which is followed by a leveling off to a slower value. The total disk mass falls off {proportional_to}t {sup -0.35} at its fastest point (where t is time) for our reference model, while the dust mass and its proxy-the infrared excess emission-fades significantly faster ({proportional_to}t {sup -0.8}). These later level off to a decay rate of M{sub tot}(t){proportional_to}t {sup -0.08} and M{sub dust}(t) or L{sub ir}(t){proportional_to}t {sup -0.6}. This is slower than the {proportional_to}t {sup -1} decay given for all three system parameters by traditional analytic models. We also compile an extensive catalog of Spitzer and Herschel 24, 70, and 100 {mu}m observations. Assuming a log-normal distribution of initial disk masses, we generate model population decay curves for the fraction of stars harboring debris disks detected at 24 {mu}m. We also model the distribution of measured excesses at the far-IR wavelengths (70-100 {mu}m) at certain age regimes. We show general agreement at 24 {mu}m between the decay of our numerical collisional population synthesis model and observations up to a Gyr. We associate offsets above a Gyr to stochastic events in a few select systems. We cannot fit the decay in the far-infrared convincingly with grain strength properties appropriate for silicates, but those of water ice give fits more consistent with the observations (other relatively weak grain materials would presumably also be successful). The oldest disks have a higher incidence of large excesses than predicted by the model; again, a plausible explanation is very late phases of high dynamical activity around a small number of stars. Finally, we constrain the variables of our numerical model by comparing the evolutionary trends generated from the exploration of the full parameter space to observations. Amongst other results, we show that erosive collisions are dominant in setting the timescale of the evolution and that planetesimals on the order of 100 km in diameter are necessary in the cascades for our population synthesis models to reproduce the observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Stellar velocity dispersion and anisotropy of the Milky Way inner halo

We measure the three components of velocity dispersion, σ{sub R}, σ{sub θ}, σ{sub ϕ}, for stars within 6 < R < 30 kpc of the Milky Way using a new radial velocity sample from the MMT telescope. We combine our measurements with previously published data so that we can more finely sample the stellar halo. We use a maximum likelihood statistical method for estimating mean velocities, dispersions, and covariances assuming only that velocities are normally distributed. The alignment of the velocity ellipsoid is consistent with a spherically symmetric gravitational potential. From the spherical Jeans equation, the mass of the Milky Way is M(R⩽12 kpc)=1.3×10{sup 11} M{sub ⊙} with an uncertainty of 40%. We also find a region of discontinuity, 15 ≲ R ≲ 25 kpc, where the estimated velocity dispersions and anisotropies diverge from their anticipated values, confirming the break observed by others. We argue that this break in anisotropy is physically explained by coherent stellar velocity structure in the halo, such as the Sgr stream. To significantly improve our understanding of halo kinematics will require combining radial velocities with future Gaia proper motions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Statistics of general functions of a Gaussian field-application to non-Gaussianity from preheating

We provide a general formula for calculating correlators of arbitrary function of a Gaussian field. This work extends the standard leading-order approximation based on the δN formalism to the case where truncation of the δN at some low order does not yield the correct answer. As an application of this formula, we investigate 2, 3 and 4-point functions of the primordial curvature perturbation generated in the massless preheating model by approximating the mapping between the curvature perturbation and the Gaussian field as a sum of the many spiky normal distribution functions as suggested by lattice calculations. We also discuss observational consequences of this case and show that trispectrum would be a key observable to search signature of preheating in the CMB map. It is found the forms of the curvature correlation functions for any δN, at the leading order in the correlator of the Gaussian field, coincide with the standard local type ones. Within this approximation, it is also found that the standard formula for the non-linearity parameters given by the product of the derivatives of the e-folding number still holds after we replace the bare e-folding number appearing in the original δN expansion with the one smoothed in the field space with a Gaussian window function.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Detecting anomalies in CMB maps: a new method

Ever since WMAP announced its first results, different analyses have shown that there is weak evidence for several large-scale anomalies in the CMB data. While the evidence for each anomaly appears to be weak, the fact that there are multiple seemingly unrelated anomalies makes it difficult to account for them via a single statistical fluke. So, one is led to considering a combination of these anomalies. But, if we ''hand-pick'' the anomalies (test statistics) to consider, we are making an a posteriori choice. In this article, we propose two statistics that do not suffer from this problem. The statistics are linear and quadratic combinations of the a{sub ℓ m}'s with random co-efficients, and they test the null hypothesis that the a{sub ℓ m}'s are independent, normally-distributed, zero-mean random variables with an m-independent variance. The motivation for considering multiple modes is this: because most physical models that lead to large-scale anomalies result in coupling multiple ℓ and m modes, the ''coherence'' of this coupling should get enhanced if a combination of different modes is considered. In this sense, the statistics are thus much more generic than those that have been hitherto considered in literature. Using fiducial data, we demonstrate that the method works and discuss how it can be used with actual CMB data to make quite general statements about the incompatibility of the data with the null hypothesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Solving puzzles of GW150914 by primordial black holes

The black hole binary properties inferred from the LIGO gravitational wave signal GW150914 posed several serious problems. The high masses and low effective spin of black hole binary can be explained if they are primordial (PBH) rather than the products of the stellar binary evolution. Such PBH properties are postulated ad hoc but not derived from fundamental theory. We show that the necessary features of PBHs naturally follow from the slightly modified Affleck-Dine (AD) mechanism of baryogenesis. The log-normal distribution of PBHs, predicted within the AD paradigm, is adjusted to provide an abundant population of low-spin stellar mass black holes. The same distribution gives a sufficient number of quickly growing seeds of supermassive black holes observed at high redshifts and may comprise an appreciable fraction of Dark Matter which does not contradict any existing observational limits. Testable predictions of this scenario are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Stochastic particle production in a de Sitter background

Here, we explore non-adiabatic particle production in a de Sitter universe for a scalar spectator field, by allowing the effective mass m 2 (t) of this field and the cosmic time interval between non-adiabatic events to vary stochastically. Two main scenarios are considered depending on the (non-stochastic) mass M of the spectator field: the conformal case with M 2 =2H 2 , and the case of a massless field. We make use of the transfer matrix formalism to parametrize the evolution of the system in terms of the "occupation number", and two phases associated with the transfer matrix; these are used to construct the evolution of the spectator field. Assuming short-time interactions approximated by Dirac-delta functions, we numerically track the change of these parameters and the field in all regimes: sub- and super-horizon with weak and strong scattering. In all cases a log-normally distributed field amplitude is observed, and the logarithm of the field amplitude approximately satisfies the properties of a Wiener process outside the horizon. We derive a Fokker-Planck equation for the evolution of the transfer matrix parameters, which allows us to calculate analytically non-trivial distributions and moments in the weak-scattering limit.

79 ASTRONOMY AND ASTROPHYSICS↗