Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle statistics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

High Energy Particle Excitations as a Measure of Electronic Crosstalk in MODIS and VIIRS VIS/NIR Bands

The MODIS instruments aboard the Terra and Aqua satellites and the VIIRS instruments aboard the SNPP and NOAA-20 satellites each contain several arrays of Si detectors that measure Earth-reflected radiance in the visible and near-infrared spectral range. Even in the absence of incident light, the Si detectors are occasionally excited by high energy charged particles that pass through the spacecraft. These particle radiation events are, fortunately, infrequent enough that they do not lead to significant degradation of the detectors and they do not have a significant impact on the Earth scene radiance images. On the other hand, they are frequent enough that the cumulative data from many years on orbit may provide valuable diagnostic information about the sensors. In this paper, we provide some basic statistics on the frequency and magnitude of the particle excitation events for MODIS and VIIRS and explore the usefulness of this data as a measure of electronic crosstalk. Large amounts of crosstalk can degrade the quality of the Earth images, so it is crucial to have methods to characterize and correct for it on-orbit, which has previously been done for MODIS using lunar image analysis. The particle excitations can manifest as single-pixel spikes in the otherwise dark space view background, which may be an ideal source for evaluating crosstalk. We derive crosstalk coefficients between the NIR band detectors of Terra MODIS, and compare them to coefficients previously derived from lunar observations. The same approach is applied to SNPP VIIRS, which does not show any significant electronic crosstalk. While the HgCdTe detectors used in the MODIS and VIIRS infrared bands can also be excited by particle radiation, the magnitudes of the excitations are much smaller compared to the Si detectors and in general are not large enough to be useful for examining crosstalk.

MODIS↗
DOE OSTI2024

coh3

CoH3 (CoH ver.3) is an optical model, exciton pre-equilibrium, and Hauser-Feshbach statistical model code, which calculates nuclear reaction cross sections for medium to heavy targets in the keV to MeV energy region. This program is written in standard C++, divided into approximately 200 source and header files. CoH solves the Schroedinger equation for optical potentials defined in the code, and calculates differential elastic scattering, reaction, and total cross sections, for neutron, proton, deuteron, triton, 3He, and alpha-particle. Deformed optical potentials are solved with the coupled-channels method, in which the ground state rotational band members, or vibrational phonon states are coupled. The optical model gives particle transmission coefficients that are fed into the statistical model calculations. CoH includes the pre-equilibrium model (exciton model), the direct/semidirect capture model, and the multi-stage Hauser-Feshbach statistical decay with width fluctuation correction based on the Gaussian orthogonal ensemble. For weakly coupled levels, the DWBA (distorted wave Born approximation) method is used to calculate the direct inelastic scattering process to the excited states.

Kawano, Toshihiko↗
DOE OSTI2021

Formation of Aerosol Nanoparticles by Gas-Phase Hydrolysis Reaction of Uranium Hexafluoride

The aerosol physics of uranyl particle formation has been addressed in this research using advanced aerosol instrumentation and an aerosol dynamics model. Based on the research works, we conclude that the formation and growth of aerosol particles by gas-phase UF6 hydrolysis strongly depends on the availability of water molecules in our reactor conditions. The total number concentration of the UO 2 F 2 particulate material that could be produced in the hydrolysis reaction is also regulated primarily by the availability of water molecule concentration. The higher the water molecule concentration, the higher the number and the larger the size of UO 2 F 2 aerosol particles that could be produced in a reactor custom-built at ORNL. Although the aerosol reactor was enabling the study of particle formation kinetics, the instrumentation was still insufficient in characterizing the chemical composition of the produced particles as well as the time-dependent evolution of the particulate species. The temporal evolution could impact the eventual fate of the particles upon release to the environment (i.e., the physio-chemical transformation, transport, and removal). On uranyl particle formation kinetics, we found that the growth rates of aerosol particles appeared to approach a single number in the range of 0.05 ± 0.03 - 0.08 ± 0.04 nm/s, statistically, as the ω value becomes smaller than 1. The size of primary particles from the UF6 hydrolysis at water-deprived condition was estimated to be 3.6 ± 0.4 nm; the higher the availability of water molecules, the larger the primary particles. The ability to precisely control the availability of water molecules in the reaction could lead to the production of nearly monodispersed aerosol particles. In other words, the result suggests that one can precisely manipulate the size of UO 2 F 2 aerosol particles by controlling the water vapor availability and interaction of water molecules with U F6 in the reaction. This finding has significant implications in the engineering manufacturing of fuel powder materials and possibly to future development and deployment of an environmental sampling apparatus.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1991

Preferential concentration of particles by turbulence

Direct numerical simulation of isotropic turbulence was used to investigate the effect of turbulence on the concentration fields of heavy particles. The hydrodynamic field was computed using 64-cubed points, and a statistically stationary flow was obtained by forcing the low-wave-number components of the velocity field. The particles used in the simulations were time advanced according to Stokes drag law and were also assumed to be much more dense than the fluid. Properties of the particle cloud were obtained by following the trajectories of 1,000,000 particles through the simulated flow fields. Three values of the ratio of the particle time constant to large-scale turbulence time scale were used in the simulations: 0.075, 0.15, and 0.52. The simulations show that the particles collect preferentially in regions of low vorticity and high strain rate. This preferential collection was most pronounced for particles with the instantaneous number density as much as 25 times the mean value. The fact that dense particles collect in regions of low vorticity and high strain in turn implies that turbulence may actually inhibit rather than enhance mixing of particles.

Squires, Kyle D.↗
DOE OSTI2022

Particle transport in Markov stochastic mixtures with spatial gradients

Markov geometries are a prototype class of stochastic media that are widely used to model complex disordered systems. In a series of recent works, we have shown that spatially homogeneous (but possibly non-isotropic) Poisson tessellations can be conveniently used in order to generate three-dimensional realizations of Markov media and thus characterize the features of particle transport e.g. in fragmented fuel elements following severe nuclear accidents or Rayleigh-Taylor turbulent mixing for inertial confinement fusion. In this paper we show that Poisson tessellations can be extended to take into account the presence of spatial gradients, which commonly occur in real-world applications due to material stratification. For this purpose, we will provide an explicit construction of spatially-non-homogeneous Poisson tessellations, fully characterize their statistical properties, and finally illustrate the behaviour of particle transport in such random media. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Reconstruction of effective potential from statistical analysis of dynamic trajectories

The broad incorporation of microscopic methods is yielding a wealth of information on the atomic and mesoscale dynamics of individual atoms, molecules, and particles on surfaces and in open volumes. Analysis of such data necessitates statistical frameworks to convert observed dynamic behaviors to effective properties of materials. Here, we develop a method for the stochastic reconstruction of effective local potentials solely from observed structural data collected from molecular dynamics simulations (i.e., data analogous to those obtained via atomically resolved microscopies). Using the silicon vacancy defect in graphene as a model, we apply the statistical framework presented herein to reconstruct the free energy landscape from the calculated atomic displacements. Evidence of consistency between the reconstructed local potential and the trajectory data from which it was produced is presented, along with a quantitative assessment of the uncertainty in the inferred parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Single electron in a storage ring: a probe into the fundamental properties of synchrotron radiation and a powerful diagnostic tool

An experimental study of undulator radiation generated by a single electron was carried out at the Integrable Optics Test Accelerator (IOTA) storage ring at Fermilab. Photons were detected with a single-photon avalanche diode (SPAD) at an average rate of 20 kHz, corresponding to one detection every few hundred revolutions. Events were continuously recorded by a picosecond event timer for as long as 5 minutes at a time. The collected data were used to directly test the predictions of classical Poisson photostatistics. The study was motivated in part by a previous observation of sub-Poissonian statistics in a similar experiment with synchrotron radiation. In addition, we show that the time series of recorded events can be used to study the synchrotron motion of a single electron, yielding dynamical parameters such as the phase jitter of the radiofrequency cavity and the dependence of the synchrotron period on amplitude. Research plans for further experiments are also discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Mechanisms of energetic-particle transport in magnetically confined plasmas

Super-thermal ions and electrons occur in both space and fusion plasmas. Because these energetic particles (EP) have large velocities, EP orbits necessarily deviate substantially from magnetic surfaces. Orbits are described by conserved constants of motion that define topological boundaries for different orbit types. Electric and magnetic field perturbations produced by instabilities can disrupt particle orbits, causing the constants of motion to change. The statistics of the “kicks” associated with these perturbations determines the resulting cross field transport. A unifying theme of this tutorial is the importance of the perturbation’s phase at the particle’s position Θ = k • r – ωt, where k and ω are the wavevector and frequency of the perturbation, r is the EP position, and t is the time. A distinction is made between field perturbations that resonate with an aspect of the orbital motion and those that do not. Resonance occurs when the wave phase returns to its initial value in an integer multiple of an orbital period. Convective transport occurs when resonant particles experience an unvarying wave phase. Alternatively, multiple wave-particle resonances usually decorrelate the phase, resulting in diffusive transport. Large orbits increase the number of important resonances and can cause chaotic orbits even for relatively small amplitude waves. In contrast, in the case of non-resonant perturbations, orbital phase averaging reduces transport. Here, large field perturbations introduce additional effects, including nonlinear resonances at fractional values of the orbital motion. In summary, large orbits are a blessing and a curse: For non-resonant modes, orbit-averaging reduces transport but, for resonant transport, large orbits facilitate jumps across topological boundaries and enhance the number of important resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dynamic Rubies

Algorithmic art selected as a winner of the contest Pixel by Pixel: The Art of CHANCE, to be published in the April 2024 issue of CHANCE, the magazine of the American Statistical Association.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Lagrangian and Eulerian drag models that are consistent between Euler-Lagrange and Euler-Euler (two-fluid) approaches for homogeneous systems

The undisturbed flow of a particle is of fundamental importance since it controls both the undisturbed flow force and the perturbation force (which includes quasisteady, added-mass, and history forces). Here we use the pairwise interaction extended point particle framework to evaluate the undisturbed flow of each particle through superposition of the perturbation flow induced by all its neighbors. This approach allows calculation of various statistics related to undisturbed fluid velocity under conditions of both stationary and nonstationary particles. In a random distribution of stationary particles, while the macroscale undisturbed flow is slowly varying, the microscale undisturbed flow that arises due to the perturbation flow of neighbors varies substantially from one particle to another and this in turn leads to large variation in the hydrodynamic force exerted on the particles. The effect of particle motion is generally to increase the particle-to-particle variation in the undisturbed fluid velocity of the particles. We observe that this increase is greater for the transverse component than for the streamwise component. As a result, with increasing random particle motion, the distribution of undisturbed fluid-velocity fluctuation becomes isotropic. Three different normalized forces are defined: $Φ^L$ is the Lagrangian normalized force on an individual particle suitable for application in a microscale-informed Euler-Lagrange simulation, $Φ^E$ is the Eulerian normalized average force suitable for application in an Euler-Euler simulation, and $Φ^{LE}$ is the Lagrangian normalized force on an individual particle suitable for application in the standard Euler-Lagrange simulation. In this work, we establish precise relations between these different definitions. The drag laws developed based on particle-resolved direct numerical simulation results and experiments are appropriate for application only as the Eulerian normalized average force. We introduce the force consistency relation and use it to obtain an expression for $Φ^L$, which when applied to each particle and averaged over all the particles equals $Φ^E$. The results are first obtained in the limit of stationary particles and then extended to the general case of nonstationary particles.

42 ENGINEERING↗
NASA NTRS1977

Computer simulation of the velocity diffusion of cosmic rays

Monte Carlo simulation experiments were performed in order to study the velocity diffusion of charged particles in a static turbulent magnetic field. By following orbits of particles moving in a large ensemble of random magnetic field realizations with suitable chosen statistical properties, a pitch-angle diffusion coefficient is derived. Results are presented for a variety of particle rigidities and rms random field strengths and compared with the predictions of standard quasi-linear theory and the nonlinear partially averaged field theory.

Kaiser, T. B.↗
NASA NTRS1978

Computer simulation of the velocity diffusion of cosmic rays

Monte Carlo simulation experiments have been performed in order to study the velocity diffusion of charged particles in a static turbulent magnetic field. By following orbits of particles moving in a large ensemble of random magnetic field realizations with suitably chosen statistical properties, a pitch-angle diffusion coefficient is derived. Results are presented for a variety of particle rigidities and rms random field strengths and compared with the predictions of standard quasi-linear theory and the nonlinear partially averaged field theory.

Kaiser, T. B.↗
NASA NTRS2020

Upstream Ultra-Low Frequency Waves Observed by MESSENGER's Magnetometer: Implications for Particle Acceleration at Mercury's Bow Shock

We perform the first statistical analysis of the main properties of waves observed in the 0.05–0.41 Hz frequency range in the Hermean foreshock by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) Magnetometer. Although we find similar polarization properties to the “30 second” waves observed at the Earth's foreshock, the normalized wave amplitude (δ B /|B0| ∼0.2) and occurrence rate (∼0.5%) are much smaller. This could be associated with relatively lower backstreaming proton fluxes, the smaller foreshock size and/or less stable solar wind (SW) conditions around Mercury. Furthermore, we estimate that the speed of resonant backstreaming protons in the SW reference frame (likely source for these waves) ranges between 0.95 and 2.6 times the SW speed. The closeness between this range and what is observed at other planetary foreshocks suggests that similar acceleration processes are responsible for this energetic population and might be present in the shocks of exoplanets.

N. Romanelli↗
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
NASA NTRS1994

Statistical characteristics of bursty bulk flow events

Using a common methodology to analyze data from the Active Magnetospheric Particle Tracer Explorer/Ion Release Module (AMPTE/IRM) and International Sun-Earth Explorer 2 (ISEE 2) satellites we report on the statistical properties of bursty bulk flow events (BBFs) in the inner plasma sheet (IPS). A positive correlation between BBFs and the AE index suggests that BBFs are predominantly geomagnetically active time phenomena. Earthward BBFs are more frequent close to midnight and away from Earth, up to a distance of approximately 19 R(sub E). Tailward BBFs are very infrequent in the IRM data set and somewhat less infrequent in the ISEE 2 data set in the region of the satellites' spatial overlap, possibly due to the more active conditions prevailing during the ISEE 2 mission in that region. However, in both data sets the ratio of tailward to earthward BBFs increases with distance from Earth; more than 20% of all BBFs are anti-sunward tailward of X = -19 R(sub E) in the ISEE 2 data set. BBFs are responsible for 60-100% of the measured earthward transport of mass, energy and magnetic flux past the satellite in the regions of maximum occurrence rate, even though they last approximately 10-15% of the IPS observation time there. Thus BBFs represent the primary transport mechanism at those regions. The one-to-one correspondence between BBFs and substorm phase, as well as the relative contribution of BBFs to the total transport observed during substorms are questions that await further investigation based on multi instrument studies of individual events.

Angelopoulos, V.↗