Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1983

Charged particle distributions in Jupiter's magnetosphere

In situ data from the Pioneer and Voyager spacecraft, supplemented by earth-based observations and theoretical considerations, are used as the basis for the present quantitative, compact model of the 1 eV-several MeV charged particle distribution in the Jovian magnetosphere. The thermal plasma parameters of convection speed, number density, and characteristic energy, are specified as functions of position for electrons and for the ion species H(+), O(+), O(2+), S(+), S(2+), S(3+), and Na(+). Major features of the magnetic field, thermal plasma, and trapped particle distributions, are modeled and results for each plasma region are compared with observed spectra. Comparisons show that the model represents the data to within a factor of 2 + or - 1, except where time variations are significant. Practical applications of the model to spacecraft near Jupiter are given.

Divine, N.↗
NASA NTRS1980

The effect of injection location on the spectrum of energetic magnetospheric particles

The role of local particle injections and accelerations in determining energetic particle spectra used as an indication of the radial diffusion of trapped magnetospheric particles is assessed. An idealized one-dimensional steady-state model of magnetospheric radial transport in which diffusion is balanced against particle sources and sinks is used to illustrate the effects of particle injection at a point and over a band of radii in which the observation point is immersed in particle spectra. For an injection spectrum uniformly distributed in space and a step function in energy, it is shown that the energy dependence of the measured spectrum is determined not only by adiabatic energization of the input spectrum but by the spatial structure of the injection and the radial dependence of the diffusion coefficient as well. The relevance of the results for observations of particle spectra in the terrestrial and Jovian magnetospheres is also considered

Birmingham, T. J.↗
NASA NTRS2003

A New Look into the Effect of Large Drops on Radiative Transfer Process

Recent studies indicate that a cloudy atmosphere absorbs more solar radiation than any current 1D or 3D radiation model can predict. The excess absorption is not large, perhaps 10-15 W/sq m or less, but any such systematic bias is of concern since radiative transfer models are assumed to be sufficiently accurate for remote sensing applications and climate modeling. The most natural explanation would be that models do not capture real 3D cloud structure and, as a consequence, their photon path lengths are too short. However, extensive calculations, using increasingly realistic 3D cloud structures, failed to produce photon paths long enough to explain the excess absorption. Other possible explanations have also been unsuccessful so, at this point, conventional models seem to offer no solution to this puzzle. The weakest link in conventional models is the way a size distribution of cloud particles is mathematically handled. Basically, real particles are replaced with a single average particle. This "ensemble assumption" assumes that all particle sizes are well represented in any given elementary volume. But the concentration of larger particles can be so low that this assumption is significantly violated. We show how a different mathematical route, using the concept of a cumulative distribution, avoids the ensemble assumption. The cumulative distribution has jumps, or steps, corresponding to the rarer sizes. These jumps result in an additional term, a kind of Green's function, in the solution of the radiative transfer equation. Solving the cloud radiative transfer equation with the measured particle distributions, described in a cumulative rather than an ensemble fashion, may lead to increased cloud absorption of the magnitude observed.

Marshak, Alexander↗
DOE OSTI2019

The Fall Speed Variability of Similarly Sized Ice Particle Aggregates

The terminal velocity ( V t ) of ice hydrometeors is of high importance to atmospheric modeling. V t is governed by the physical characteristics of a hydrometeor, including mass and projected area, as well as environmental conditions. When liquid hydrometeors coalesce to form larger hydrometeors, the resulting hydrometeor can readily be characterized by its spherical or near-spherical shape. For ice hydrometeors, it is more complicated because of the variability of ice shapes possible in the atmosphere as well as the inherent randomness in the aggregation process, which leads to highly variable characteristics. The abundance of atmospheric processes affecting ice particle dimensional characteristics creates potential for highly variable V t for ice particles that are predicted or measured to be of the “same size.” In this article we explore the variability of ice hydrometeor V t both theoretically and through the use of experimental observations. Theoretically, the variability in V t is investigated by analyzing the microphysical characteristics of randomly aggregated hexagonal shapes. The modeled dimensional characteristics are then compared to aircraft probe measurements to constrain the variability in atmospheric ice hydrometeor V t . Results show that the spread in V t can be represented with Gaussian distributions relative to a mean. Variability expressed as the full width at half maximum of the normalized Gaussian probability distribution function is around 20%, with somewhat higher values associated with larger particle sizes and warmer temperatures. Field campaigns where mostly convective clouds were sampled displayed low variability, while Arctic and midlatitude winter campaigns showed broader V t spectra.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Rainbow Fourier Transform

We present a novel technique for remote sensing of cloud droplet size distributions. Polarized reflectances in the scattering angle range between 135deg and 165deg exhibit a sharply defined rainbow structure, the shape of which is determined mostly by single scattering properties of cloud particles, and therefore, can be modeled using the Mie theory. Fitting the observed rainbow with such a model (computed for a parameterized family of particle size distributions) has been used for cloud droplet size retrievals. We discovered that the relationship between the rainbow structures and the corresponding particle size distributions is deeper than it had been commonly understood. In fact, the Mie theory-derived polarized reflectance as a function of reduced scattering angle (in the rainbow angular range) and the (monodisperse) particle radius appears to be a proxy to a kernel of an integral transform (similar to the sine Fourier transform on the positive semi-axis). This approach, called the rainbow Fourier transform (RFT), allows us to accurately retrieve the shape of the droplet size distribution by the application of the corresponding inverse transform to the observed polarized rainbow. While the basis functions of the proxy-transform are not exactly orthogonal in the finite angular range, this procedure needs to be complemented by a simple regression technique, which removes the retrieval artifacts. This non-parametric approach does not require any a priori knowledge of the droplet size distribution functional shape and is computationally fast (no look-up tables, no fitting, computations are the same as for the forward modeling).

Alexandrov, Mikhail D.↗
DOE OSTI2019

Hybrid particle-in-cell simulations of laser-driven plasma interpenetration, heating, and entrainment

Kinetic-ion, quasineutral, fluid-electron particle-in-cell simulations of interpenetrating carbon–carbon plasma flows in 2D RZ cylindrical geometry are presented. The simulations are initialized with solid density targets that are subsequently irradiated by 10 14 W/cm 2 intensity lasers using a raytracing package. The ablation, interpenetration, heating, slowing, entrainment, and stagnation of the plasma flows evolve self-consistently within the code. The particle density, velocity phase space, and fits to the velocity distribution functions are used, along with analytical collisional stopping rates, to interpret the dynamics of the flow evolution. Comparisons to multifluid simulations are described and used to highlight ion-kinetic effects in the setup. Synthetic Thomson scattering diagnostic signals are generated using detailed knowledge of the plasma distribution functions. The large scale of the system, 1 × 1 mm for 2 ns, and the detailed dynamics extracted demonstrate that such hybrid codes are powerful tools for the design and evaluation of laboratory-scale high-energy-density plasma physics experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗