Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle distributions”

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 487 records · Page 27

Corona splashing triggered by a loose monolayer of particles

In nature, high-speed raindrops often impact and spread on particulate surfaces (e.g., soil, plant leaves with spores or pollen). We study the dynamics of droplet impact on a loosely packed monolayer of particles by combining experimental and mathematical approaches. We find that the presence of mobile particles lowers the critical impact velocity at which the droplet exhibits corona splashing, as the particle area fraction is systematically increased. We rationalize this experimental observation by considering the jamming of frictional particles at the spreading rim. Elucidating the splashing transition of the drop on a particulate bed can lead to a better understanding of soil loss and erosion from falling raindrops

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A mathematical framework for ejecta cloud dynamics with application to source models and piezoelectric mass measurements

We present a mathematical framework for describing the dynamical evolution of an ejecta cloud generated by a generic ejecta source model. We consider a piezoelectric sensor fielded in the path of an ejecta cloud, for experimental configurations in which the ejecta are created at a singly shocked planar surface and fly ballistically through vacuum to the stationary sensor. To do so, we introduce the concept of a time- and velocity-dependent ejecta “areal mass function.” We derive expressions for the analytic (“true”) accumulated ejecta areal mass at the sensor and the measured (“inferred”) value obtained via the standard method for analyzing piezoelectric voltages. In this way, we derive an exact expression and upper bound for the error imposed upon a piezoelectric ejecta mass measurement (in a perfect system) by the assumption of instantaneous creation, which is commonly required for momentum diagnostic analyses. This error term is zero for truly instantaneous source models; otherwise, the standard piezoelectric analysis is guaranteed to overestimate the true mass. When combined with a piezoelectric dataset, this framework provides a unique solution for the ejecta particle velocity distribution, subject to the assumptions inherent in the data analysis. The framework also leads to strong boundary conditions that any ejecta source model must satisfy in order to be consistent with apparently global properties of piezoelectric measurements from a wide range of experiments. We demonstrate this methodology by applying it to the Richtmyer–Meshkov instability+self-similar velocity distribution ejecta source model currently under development at Los Alamos National Laboratory.

97 MATHEMATICS AND COMPUTING↗

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-yield and high-angular-fluence neutron generation from deuterons accelerated by laser-driven collisionless shock

Since a bright collimated neutron source is an essential tool for global security missions and fundamental scientific research. In this paper, we study a compact high-yield and high-angular-fluence neutron source particularly suitable for high-energy neutron applications utilizing the breakup reaction of laser-driven deuterons in a 9 Be converter. The neutron generation scaling from such a reaction is used to guide the choice and optimization of the acceleration process for bulk ions in a low density CD 2 foam. In particular, the collisionless shock acceleration mechanism is exploited with proper choice in the laser and target parameter space to accelerate these ions toward energies above the temperature of the distribution. Particle-in-cell and Monte Carlo simulations are coupled here to investigate this concept and possible adverse effects as well as the contribution from the surface ions accelerated and the optimal converter design. The simulation results indicated that our design can be a practical approach to increase both the neutron yield and angular fluence of laser-driven neutron sources, reaching >10 11 neutron/pulse (or >10 8 neutron/J) and >10 11 neutron/sr (or >10 8 neutron/sr/J) with present-day kJ-class high-power lasers. Such developments will advance fundamental neutron science, high precision radiography, and other global security applications with laser-driven sources.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An investigation of shock formation vs shock mitigation of colliding plasma jets

Here, this work studies the interaction between colliding plasma jets to understand regimes in which jet merging results in shock formation vs regimes in which the shock formation is mitigated due to the collisionless interpenetration of the jets. A kinetic model is required for this study because fluid models will always produce a shock upon the collision of plasma jets. The continuum-kinetic, Vlasov–Maxwell–Dougherty model with one velocity dimension is used to accurately capture shock heating, along with a novel coupling with a moment equation to evolve perpendicular temperature for computational efficiency. As a result, this relatively inexpensive simulation can be used for detailed scans of the parameter space toward predictions of shocked vs shock-mitigated regimes, which is of interest for several fusion concepts such as plasma-jet-driven magneto-inertial fusion, high-energy-density plasmas, astrophysical phenomena, and other laboratory plasmas. The initial results obtained using this approach are in agreement with the preliminary outcomes of the plasma liner experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of a central barrier potential on asymmetry-induced transport in a coaxial Malmberg–Penning trap

Experiments are presented measuring the radial particle flux produced when a central barrier potential is applied to the central ring of a coaxial Malmberg-Penning trap at the same time as the voltages producing the main asymmetry potential $\phi_1(r)\cos{(kz)}\cos{(l\theta-\omega t)}$. When a negative DC voltage is applied to the entire central ring, the peak flux produced by the main asymmetry is reduced by a factor $e^{(V_{eff}/1.2)}$ where $V_{eff}$ is the effective applied potential in volts. This barrier potential also increases the frequency of this peak flux. When asymmetric $\pm$ voltages are applied to the two azimuthally-divided halves of the central ring, DC voltages increase the radial flux while AC voltages decrease it. It is shown that these asymmetric barrier voltages produce their own transport and that the resulting modification of the plasma may be a factor in changing the observed flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mapping the structural trends in zinc aluminosilicate glasses

The structure of zinc aluminosilicate glasses with the composition (ZnO) x (Al 2 O 3 ) y (SiO 2 ) 1-x-y , where 0 ≤ x < 1, 0 ≤ y < 1, and x + y < 1, was investigated over a wide composition range by combining neutron and high-energy x-ray diffraction with 27 Al magic angle spinning nuclear magnetic resonance spectroscopy. The results were interpreted using an analytical model for the composition-dependent structure in which the zinc ions do not act as network formers. Four-coordinated aluminum atoms were found to be in the majority for all the investigated glasses, with five-coordinated aluminum atoms as the main minority species. Mean Al–O bond distances of 1.764(5) and 1.855(5) Å were obtained for the four- and five-coordinated aluminum atoms, respectively. The coordination environment of zinc was not observed to be invariant. Instead, it is dependent on whether zinc plays a predominantly network-modifying or charge-compensating role and, therefore, varies systematically with the glass composition. The Zn–O coordination number and bond distance were found to be 4.36(9) and 2.00(1) Å, respectively, for the network-modifying role vs 5.96(10) and 2.08(1) Å, respectively, for the charge-compensating role. The more open coordination environment of the charge-compensator is related to an enhanced probability of zinc finding bridging oxygen atoms as nearest-neighbors, reflecting a change in the connectivity of the glass network comprising four-coordinated silicon and aluminum atoms as the alumina content is increased.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Discovering invariant spatial features in electron energy loss spectroscopy images on the mesoscopic and atomic levels

Over the last two decades, Electron Energy Loss Spectroscopy (EELS) imaging with a scanning transmission electron microscope has emerged as a technique of choice for visualizing complex chemical, electronic, plasmonic, and phononic phenomena in complex materials and structures. The availability of the EELS data necessitates the development of methods to analyze multidimensional data sets with complex spatial and energy structures. Traditionally, the analysis of these data sets has been based on analysis of individual spectra, one at a time, whereas the spatial structure and correlations between individual spatial pixels containing the relevant information of the physics of underpinning processes have generally been ignored and analyzed only via the visualization as 2D maps. Here, we develop a machine learning-based approach and workflows for the analysis of spatial structures in 3D EELS data sets using a combination of dimensionality reduction and multichannel rotationally invariant variational autoencoders. This approach is illustrated for the analysis of both the plasmonic phenomena in a system of nanowires and in the core excitations in functional oxides using low loss and core-loss EELS, respectively. The code developed in this manuscript is open sourced and freely available and provided as a Jupyter notebook for the interested reader.

36 MATERIALS SCIENCE↗

Interactions of laser-driven tin ejecta microjets over phase transition boundaries

Ejecta microjets offer an experimental methodology to study high-speed particle laden-flow interactions, as microjets consist of millions of particulates traveling at velocities of several kilometers per second and are easily generated by most common shock drives. Previous experiments on the OMEGA Extended Performance laser found that collisions between two counter-propagating laser-driven tin ejecta microjets varied as a function of drive pressure; jets generated near shock pressures of 10 GPa passed through each other without interacting, whereas jets generated at shock pressures of over 100 GPa interacted strongly, forming a cloud around the center interaction point. In this paper, we present a more systematic scan of tin ejecta microjet collisions over intermediate pressure regimes to identify how and at what shock pressure interaction behavior onsets. Radiographs of interacting microjets at five different laser drive energies qualitatively demonstrate that interaction behavior onsets slowly as a function of laser drive energy. Quantitative mass and density metrics from each radiograph provide trends on jet characteristics and collisional mass dispersion. It is observed that jetting mass, jet densities, and mass dispersion increase with increasing drive pressures and that the increased jet density at the higher drive energies may account for the increased mass dispersion. This work provides an important step in the understanding of tin ejecta microjet collisions and paves the way for future studies on the physics dominating high-speed particle-laden flow interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An implicit solution for Asay foil trajectories generated by separable, sustained-production ejecta source models

We present a simple implicit solution for the time-dependent trajectory of a thin Asay foil ejecta diagnostic for the general case where the impinging ejecta cloud is generated by a source function characterized by an arbitrary (sustained) time dependence and a time-independent (stationary) particle velocity distribution. In the limit that the source function time dependence becomes a delta function, this solution—which is amenable to rapid numerical calculations of arbitrary accuracy—exactly recovers a previously published solution for the special case of instantaneous ejecta production. We also derive simple expressions for the free-surface arrival (catch-up) time as well as the true ejecta areal mass accumulation on the accelerating foil and place bounds on the level of error incurred when applying instant-production mass solutions to a sustained-production trajectory. We demonstrate these solutions with example calculations for hypothetical source functions spanning a wide range of ejecta production durations, velocity distributions, and temporal behaviors. These calculations demonstrate how the foil trajectory is often insensitive to the temporal dependence of the source function, instead being dominated by the velocity distribution. We quantify this insensitivity using a “compatibility score” metric. Under certain conditions, one may capitalize upon this insensitivity to obtain a good approximation of the second integral of the velocity distribution from the observed foil trajectory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Poor confinement in stellarators at high energy II

This work discusses particle resonances in stellarators. There are two important effects. First, resonances provide locations for Alfvén modes, which can cause significant particle loss. Second, if a resonance matches the period of the modulation of the equilibrium field magnitude it produces islands in particle orbits increasing in size with particle energy, also possibly causing loss. A previous publication suffered from a restriction to the consideration of deeply passing particle orbits, clearly showing the resonance islands. Large particle magnetic moment increases the effect of particle cross field drift, essential for this phenomenon, so passing particles with larger values of magnetic moment give larger islands. It is possible to avoid designs exhibiting such resonances.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium↗

A quantitative imaging framework for lithium morphology: Linking deposition uniformity to cycle stability in lithium metal batteries

Characterizing the morphology of lithium (Li) is crucial for developing long-lasting lithium metal batteries. It is well established that more uniform Li deposition correlates with better cell performance. Li morphology is often characterized through qualitative analysis of scanning electron microscopy (SEM) images; however, there are no widely accepted metrics to quantitatively describe deposition uniformity. Here, we propose a framework to quantify uniformity through SEM image analysis via the index of dispersion (ID) metric, which is defined and presented in the context of Li metal batteries. We also explore experimental impacts of sampling protocols onIDmeasurements. Our results demonstrate that theIDmetric is highly sensitive to variations in deposition uniformity, including the coexistence and uniformity of multiple morphologies, uniformity within a single morphology, and particle size distribution uniformity. Furthermore, it is demonstrated that uniformity, as measured by theID, can be related to the average potential of Li||Li symmetric cells over cycling. Higher capacity cycling leads to more pronounced changes in bothIDand average cell potential. Local minima/maxima are found consistently in bothIDand average cell potential immediately before cells short-circuit, which we suggest may indicate a collapse of the microstructure prior to failure. We put forward this framework as a more robust approach to quantify Li deposition uniformity, advancing the development of Li metal batteries that are safer and longer lasting.

Science & Technology - Other Topics↗

Preliminary Modeling of Chloride Deposition on Spent Nuclear Fuel Canisters in Dry Storage Relevant to Stress Corrosion Cracking

This study presents a method that can be used to gain information relevant to determining the corrosion risk for spent nuclear fuel (SNF) canisters during extended dry storage. Currently, it is known that stainless steel canisters are susceptible to chloride-induced stress corrosion cracking (CISCC). However, the rate of CISCC degradation and the likelihood that it could lead to a through-wall crack is unknown. This study uses well-developed computational fluid dynamics and particle-tracking tools and applies them to SNF storage to determine the rate of deposition on canisters. The deposition rate is determined for a vertical canister system and a horizontal canister system, at various decay heat rates with a uniform particle size distribution, ranging from 0.25 to 25 µm, used as an input. In all cases, most of the dust entering the overpack passed through without depositing. Most of what was retained in the overpack was deposited on overpack surfaces (e.g., inlet and outlet vents); only a small fraction was deposited on the canister itself. These results are provided for generalized canister systems with a generalized input; as such, this paper is intended to demonstrate the technique. This study is a part of an ongoing effort funded by the U.S. Department of Energy, Nuclear Energy Office of Spent Fuel Waste Science and Technology, which is tasked with doing research relevant to developing a sound technical basis for ensuring the safe extended storage and subsequent transport of SNF. This work is being presented to demonstrate a potentially useful technique for SNF canister vendors, utilities, regulators, and stakeholders to utilize and further develop for their own designs and site-specific studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Properties of different LLM-105 preparations

LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide) has been prepared at several different size quantities (~10–70 Kg) and not subjected to further purification, such as recrystallization. Chemical and physical properties were compared—small-scale safety testing (impact, friction, ESD), thermal (chemical reactivity, differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), one-dimensional time-to-explosion (ODTX)), morphological (particle-size distribution (PSD), scanning electron microscopy (SEM), powder x-ray diffraction (PXRD), skeletal density). Each preparation had similar properties to the other preparations except for PSD and PXRD. For PSD, all preparations exhibited a mono-modal distribution with a variation in median values from ~ 40 μm to over 80 μm. All PXRD values indicated the same morphology with minor variations in crystal orientations. Possible polymorphism was also observed in a few cases. SEM images and PXRD indicate all preparations to be diamond-type and X-type crystals with some twinning. Impact, friction, and ESD and thermal sensitivities values are consistent with these type of crystal habits as seen in vast literature studies.

Explosive characterization↗

Applications of emulation and Bayesian methods in heavy-ion physics

Abstract Heavy-ion collisions provide a window into the properties of many-body systems of deconfined quarks and gluons. Understanding the collective properties of quarks and gluons is possible by comparing models of heavy-ion collisions to measurements of the distribution of particles produced at the end of the collisions. These model-to-data comparisons are extremely challenging, however, because of the complexity of the models, the large amount of experimental data, and their uncertainties. Bayesian inference provides a rigorous statistical framework to constrain the properties of nuclear matter by systematically comparing models and measurements. This review covers model emulation and Bayesian methods as applied to model-to-data comparisons in heavy-ion collisions. Replacing the model outputs (observables) with Gaussian process emulators is key to the Bayesian approach currently used in the field, and both current uses of emulators and related recent developments are reviewed. The general principles of Bayesian inference are then discussed along with other Bayesian methods, followed by a systematic comparison of seven recent Bayesian analyses that studied quark-gluon plasma properties, such as the shear and bulk viscosities. The latter comparison is used to illustrate sources of differences in analyses, and what it can teach us for future studies.

Paquet, Jean-François (ORCID:0000000187368171)↗