Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle distribution functions”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A predicted energetic ion signature of a neutral line in the geomagnetic tail

Under certain assumptions one can predict the formation of a magnetic neutral line on theoretical grounds, generally as a result of an unstable current sheet. Neutral points also seen to form spontaneously in a wide variety of plasma simulatioons, from MHD to fully kinetic. Hence, such magnetic null regions are central to some models of magnetotail energization and magnetic substorms. Yet the observational evidence for such regions remains controversial. This paper presents a new signature of a neutral line: the field-aligned ion beam characteristic of current sheet acceleration becomes dispersed into a 'ridge' in the velocity space distribution function as a neutral line is approached. The results were obtained using single-particle simulations in a fixed magnetic field model. The origin of the ridge is shown to be due to initially earthward going ions which pass through the neutral line region. These ions come from a smaller part of the initially antiearthward flowing distribution, generating the depletion of the distribution function below the ridge.

Martin, R. F., Jr.↗
NASA NTRS1981

Analytical description of charged particles transport along arbitrary guiding-field configurations

A new description of focused transport was developed which is valid for arbitrary spatial dependences of both the scattering mean free path and the focusing length. In particular, it describes the supercoherent transition from the coherent modes that occurs when focusing is strong compared to scattering to the diffusive mode that occurs when focusing is weak. Because it specifies the pitch angle distribution in terms of relatively simple functions, the new formulation facilitates not only the detailed interpretation of solar particle events, but also the qualitative understanding of focused transport. In the steady state, it predicts that the angular distribution consists of a highly collimated 'strahl' which is superimposed on an anisotropic 'halo'. This picture provides new insight into observed configurations of steady-state particle fluxes.

Earl, J. A.↗
DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗
DOE OSTI2020

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

A Radiation Chemistry Code Based on the Greens Functions of the Diffusion Equation

Ionizing radiation produces several radiolytic species such as.OH, e-aq, and H. when interacting with biological matter. Following their creation, radiolytic species diffuse and chemically react with biological molecules such as DNA. Despite years of research, many questions on the DNA damage by ionizing radiation remains, notably on the indirect effect, i.e. the damage resulting from the reactions of the radiolytic species with DNA. To simulate DNA damage by ionizing radiation, we are developing a step-by-step radiation chemistry code that is based on the Green's functions of the diffusion equation (GFDE), which is able to follow the trajectories of all particles and their reactions with time. In the recent years, simulations based on the GFDE have been used extensively in biochemistry, notably to simulate biochemical networks in time and space and are often used as the "gold standard" to validate diffusion-reaction theories. The exact GFDE for partially diffusion-controlled reactions is difficult to use because of its complex form. Therefore, the radial Green's function, which is much simpler, is often used. Hence, much effort has been devoted to the sampling of the radial Green's functions, for which we have developed a sampling algorithm This algorithm only yields the inter-particle distance vector length after a time step; the sampling of the deviation angle of the inter-particle vector is not taken into consideration. In this work, we show that the radial distribution is predicted by the exact radial Green's function. We also use a technique developed by Clifford et al. to generate the inter-particle vector deviation angles, knowing the inter-particle vector length before and after a time step. The results are compared with those predicted by the exact GFDE and by the analytical angular functions for free diffusion. This first step in the creation of the radiation chemistry code should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.

Plante, Ianik↗
DOE OSTI2021

Surface Sampling Techniques for the Canister Deposition Field Demonstration

This report describes plans for dust sampling and analysis for the multi-year Canister Deposition Field Demonstration. The demonstration will use three commercial 32PTH2 NUHOMS welded stainless steel storage canisters, which will be stored at an ISFSI site in Advanced Horizontal Storage Modules. One canister will be unheated; the other two will have heaters to achieve canister surface temperatures that match, to the degree possible, spent nuclear fuel (SNF) loaded canisters with heat loads of 10 kW and 40 kW. Surface sampling campaigns will take place on a yearly or bi-yearly basis. The goal of the planned dust sampling and analysis is to determine important environmental parameters that impact the potential occurrence of stress corrosion cracking on SNF dry storage canisters. Specifically, the size, morphology, and composition of the deposited dust and salt particles will be quantified, as well as the soluble salt load per unit area and the rate of deposition, as a function of canister surface temperature, location, time, and orientation. Sampling locations on the canister surface will nominally include 25 locations, corresponding to 5 circumferential locations at each of the 5 longitudinal locations. At each sampling location, a 2x2 sampling grid (containing 4 sample cells) will be painted onto the metal surface. During each sampling campaign, two samples at each sampling location will be collected, in a specific routine to measure both periodic (yearly or bi-yearly) and cumulative deposition rates. For each sample, a wet and a dry sample will be collected. Wet samples will be analyzed to determine the composition of the soluble salt fraction and to estimate salt loading per unit area. Dry samples will be analyzed to assess particle size, morphology, mineralogy, and identity (e.g. for floral/faunal fragments). The data generated by this proposed sampling plan will provide detailed information on dust and salt aerosol deposits on spent nuclear fuel canister surfaces. The anticipated results include information regarding particle compositions, size distributions, and morphologies, in addition to particle deposition rates as a function of canister surface location, orientation, time, and temperature. The information gathered during the Canister Deposition Field Demonstration is critical for ongoing efforts to develop a detailed understanding of the potential for stress corrosion cracking on SNF dry storage canisters

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗
NASA NTRS2018

Microphysical Properties of Tropical Tropopause Layer Cirrus

Pervasive cirrus clouds in the tropical tropopause layer (TTL) play an important role in determining the composition of stratospheric air through dehydration of tropospheric air entering the stratosphere. This dehydration affects Earth's energy budget and climate, yet uncertainties remain regarding the microphysical processes that govern TTL cirrus. TTL cirrus were sampled with the NASA Global Hawk UAV for over 30 hr in the Western Pacific in 2014 during the Airborne Tropical TRopopause EXperiment. In situ measurements by a Fast Cloud Droplet Probe and Hawkeye probe (combination Fast Cloud Droplet Probe, Two‐Dimensional Stereo optical array probe, and Cloud Particle Imager) provided particle concentrations and sizing between 1‐ and 1,280‐μm diameter and high resolution images for habit identification. We present the variability in ice concentrations, size distributions, and habits as functions of temperature, altitude, and time since convective influence. Observed ice particles were predominantly small and quasi‐spheroidal in shape, with the percentage of quasi‐spheroids increasing with decreasing temperature. In comparison to the large fraction of the population consisting of quasi‐spheroids, faceted habits (columns, plates, rosettes, and budding rosettes) constituted a smaller percentage of the overall population and exhibited the opposite correlation with temperature. The trend of higher percentages of faceted crystals occurring at warmer temperatures may be due to diffusional growth or aggregation as particles descend through cloud, and/or the more rapid diffusional growth rate at warmer temperatures. Sampling was typically well away from deep convection, however, and very few aggregates were observed, so the trend of higher percentages of faceted habits is likely attributable to diffusional growth.

Woods, Sarah↗
DOE OSTI2023

Breakups are complicated: an efficient representation of collisional breakup in the superdroplet method

Abstract. A key constraint of particle-based methods for modeling cloud microphysics is the conservation of total particle number, which is required for computational tractability. The process of collisional breakup poses a particular challenge to this framework, as breakup events often produce many droplet fragments of varying sizes, which would require creating new particles in the system. This work introduces a representation of collisional breakup in the so-called “superdroplet” method which conserves the total number of superdroplets in the system. This representation extends an existing stochastic collisional-coalescence scheme and samples from a fragment size distribution in an additional Monte Carlo step. This method is demonstrated in a set of idealized box model and single-column warm-rain simulations. We further discuss the effects of the breakup dynamic and fragment size distribution on the particle size distribution, hydrometeor population, and microphysical process rates. Box model experiments serve to characterize the impacts of properties such as coalescence efficiency and fragmentation function on the relative roles of collisional breakup and coalescence. The results demonstrate that this representation of collisional breakup can produce a stationary particle size distribution, in which breakup and coalescence rates are approximately equal, and that it recovers expected behavior such as a reduction in precipitate-sized particles in the column model. The breakup algorithm presented here contributes to an open-source pythonic implementation of the superdroplet method, PySDM, which will facilitate future research using particle-based microphysics.

de Jong, Emily (ORCID:0000000253104554)↗
NASA NTRS1984

Characteristics of suprathermal H(+) and He(++) in plasmoids in the distant magnetotail

Rest frame distribution functions were obtained for suprathermal H(+) and He(2+) ions during 20 particle intensity increases that were observed in the distant magnetotail by ISEE-3. All of these increases were characterized by isotropic, approximately 100-keV electrons, preceded by beaming electrons and velocity-dispersed protons distributions whose signatures are typically associated with plasmoids. In each of the plasmoids, both H(+) and He(2+) are found to be convected tailward, with speeds that range from 200 to 1000 km/sec. The acceleration mechanism is found to depend on particle velocity, and it is established that the solar wind is a dominant plasma source. The plasmoids appear to be expanding at about 100 km/sec during their tailward propagation.

Gloeckler, G.↗
NASA NTRS1993

The large-scale morphology of IRAS galaxies

At present, visual inspection is the only method for comparing the large-scale morphologies in the distribution of galaxies to those in model universes generated by N-body simulations. To remedy the situation, we have developed a set of three structure functions (S1, S2, S3) that quantify the degree of large-scale prolateness, oblateness, and sphericity/uniformity of a 3-D particle distribution and have applied them to a volume-limited (less than = 4000 km/s) sample of 699 IRAS galaxies with f sub 60 greater than 1.2 Jy. To determine the structure functions, we randomly select 500 galaxies as origins of spherical windows of radius R sub w, locate the centroid of the galaxies in the window (assuming all galaxies have equal mass) and then, compute the principal moments of inertia (I sub 1, I sub 2, I sub 3) about the centroid. Each S sub i is a function of (I sub 2)/(I sub 1) and (I sub 3)/I sub 1). S1, S2, and S3 tend to unity for highly prolate, oblate, and uniform distributions, respectively and tend to zero otherwise. The resulting 500 values of S sub i at each scale R sub w are used to construct a histogram.

Babul, Arif↗
NASA NTRS1983

Modelling atmospheric aerosol backscatter at CO2 laser wavelengths. I - Aerosol properties, modeling techniques, and associated problems. II - Modeled values in the atmosphere. III - Effects of changes in wavelength and ambient conditions

The various methods of calculating the atmospheric aerosol backscattering function, beta(CO2), both from measured aerosol characteristics and from optical measurements made at other wavelengths, are discussed in detail, with limits placed on their accuracy. The most significant factor in determining beta(CO2) is found to be the aerosol size distribution and concentration; this should be known accurately for particle radii up to at least 1 micron for stratospheric particles and 5 microns for tropospheric particles. Results are then presented from the modeling of the aerosol backscattering function at a wavelength of 10.6 microns in the lowest 20 km of the atmosphere. It is found that beta(CO2) varies from 10 to the -6th per m per sr in the planetary boundary layer to less than 10 to the -11th per m per sr in the stratosphere. It is next shown that, with the exception of (NH4)2SO4-containing aerosols, whose size distributions have relatively large numbers of small particles, the variation of backscattering with CO2 wavelength is less than a factor of approximately 3. For such (NH4)2SO4 aerosol distributions, however, the variation of backscatter function with CO2 wavelengths between 9.1 and 11.1 microns may reach one order of magnitude.

Kent, G.S.↗
DOE OSTI2023

Bi-Kappa Proton Mirror and Cyclotron Instabilities in the Solar Wind

Abstract The charged particles in the solar wind are often observed to possess a nonthermal tail in the velocity distribution function, a feature that can be fitted with the Kappa model. The anisotropic, or bi-Kappa, model of protons, electrons, and other charged particles is thus adopted in the literature for interpreting the data as well as in the context of the analysis of wave–particle interactions. The present paper develops an approximate but efficient theory of the mirror and cyclotron instabilities excited by the bi-Kappa protons in the solar wind. A velocity moment-based quasi-linear theory of these instabilities is also formulated in order to investigate the saturation behavior. Applications of the formalism are made for instabilities close to the marginally unstable state, which is typical of the solar wind near 1 au.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1980

Simulation of spacecraft charging environments by monoenergetic beams

To simulate the distributed spectra of plasmas in space, the differential current density spectrum of the plasma was divided into a number of energy bands and the beam energy and current were calculated for each band to provide a piecewise reproduction of the distributed spectrum. Beam energies and current densities were chosen to match the velocity moments of the plasma distribution function. The velocity moments are averages related to physical quantities such as particle density, flux, pressure, and energy flux, and have been used extensively to characterize the measured properties of plasmas in space. Combinations of one, two, and three beams were found to match two to six velocity moments of Maxwellian distributions. A computational model was used to compare the charging of a spacecraft by plasmas with distributed spectra and by monoenergetic beams.

Halverson, W.↗
NASA NTRSDate not supplied

Tracing the Cosmic Energy Flow with Parker Solar Probe: Universal Processes

Heliophysics is inherently interdisciplinary, born from multiple subjects to understand how the Sun interacts with its surrounding environment. Heliophysicists and Astrophysicists seldom exchange ideas since the scientific infrastructure partitions research areas by spatial location in the cosmos. Parker Solar Probe (PSP) is sending back puzzling information in regions we have never been before, presenting an opportunity to search for answers in non-traditional ways that are inherently cross-disciplinary. The fundamental processes that govern the physics of the universe, such as wave-particle interactions, turbulence, magnetic reconnection, and shocks, are mechanisms to study cosmic energy flow. We present PSP observations that may enable us to assess the universality of these fundamental processes in – and across boundaries of –various space and astrophysical regimes. Specifically, we focus on the exotic shapes of proton velocity distribution functions and subsequent evolution leading to multi-scale wave generation and particle acceleration. We welcome ideas, tools, and methodologies from all disciplines to understand the systems of systems governing our universe.

J L Verniero↗
DOE OSTI2024

Bottomonium suppression in 5.02 and 8.16 TeV 𝑝-Pb collisions

Here, we compute the suppression of ϒ⁡(1⁢𝑆), ϒ⁡(2⁢𝑆), and ϒ⁡(3⁢𝑆) states in 𝑝-Pb collisions relative to 𝑝⁢𝑝 collisions, including nuclear parton distribution function (nPDF) effects, coherent energy loss, momentum broadening, and final-state interactions in the quark-gluon plasma. We employ the EPPS21 nPDFs and calculate the uncertainty resulting from variation over the associated error sets. To compute coherent energy loss and momentum broadening, we follow the approach of Arleo, Peigne, and collaborators. The 3+1⁢D viscous hydrodynamical background evolution of the quark-gluon plasma is generated by anisotropic hydrodynamics. The in-medium suppression of bottomonium in the quark-gluon plasma is computed using a next-to-leading-order open quantum system framework formulated within potential nonrelativistic quantum chromodynamics. We find that inclusion of all these effects provides a reasonable description of experimental data from the ALICE, ATLAS, CMS, and LHCb Collaborations for the suppression of ϒ⁡(1⁢𝑆), ϒ⁡(2⁢𝑆), and ϒ⁡(3⁢𝑆) as a function of both transverse momentum and rapidity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Parallel transport modeling of linear divertor simulators with fundamental ion cyclotron heating *

Abstract The Material Plasma Exposure eXperiment (MPEX) is a steady state linear device with the goal to perform plasma material interaction studies at future fusion reactor relevant conditions. A prototype of MPEX referred as ‘Proto-MPEX’ is designed to carry out research and development related to source, heating and transport concepts on the planned full MPEX device. The auxiliary heating schemes in MPEX are based on cyclotron resonance heating with radio frequency (RF) waves. Ion cyclotron heating (ICH) and electron cyclotron heating in MPEX are used to independently heat the ions and electrons and provide fusion divertor conditions ranging from sheath-limited to fully detached divertor regimes at a material target. A hybrid particle-in-cell code- PICOS++ is developed and applied to understand the plasma parallel transport during ICH in MPEX/Proto-MPEX to the target. With this tool, evolution of the distribution function of MPEX/Proto-MPEX ions is modeled in the presence of (a) Coulomb collisions, (b) volumetric particle sources and (c) quasi-linear RF-based ICH. The code is benchmarked against experimental data from Proto-MPEX and simulation data from B2.5 EIRENE. The experimental observation of ‘density-drop’ near the target in Proto-MPEX and MPEX during ICH is demonstrated and explained via physics-based arguments using PICOS++ modeling. In fact, the density drops at the target during ICH in Proto-MPEX/MPEX to conserve the flux and to compensate for the increased flow during ICH. Furthermore, sensitivity scans of various plasma parameters with respect to ICH power are performed for MPEX to investigate its role on plasma transport and particle and energy fluxes at the target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗