skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Self-consistent calculations of the electric charge, ion drag force, and the drift velocity of spherical grains using Langevin dynamics and comparisons against canonical experiments
We present trajectory simulation-based modeling to capture the interactions between ions and charged grains in dusty or complex plasmas. Our study is motivated by the need for a self-consistent and experimentally validated approach for accurately calculating the ion drag force and grain charge that determine grain collective behavior in plasmas. We implement Langevin dynamics in a computationally efficient predictor–corrector approach to capture multiscale ion and grain dynamics. Predictions of grain velocity, grain charge, and ion drag force are compared with prior measurements to assess our approach. The comparisons reveal excellent agreement to within ±20% between predicted and measured grain velocities [Yaroshenko et al., Phys. Plasmas 12, 093503 (2005) and Khrapak et al., Europhys. Lett. 97, 35001 (2012)] for 0.64, 1.25 μm grains at ∼20−500 Pa. Comparisons with the measured grain charge [Khrapak et al., Phys. Rev. E 72, 016406 (2005)] under similar conditions reveal agreement to within ∼20% as well. Measurements of the ion drag force [Hirt et al., Phys. Plasmas 11, 5690 (2004); IEEE Trans. Plasma Sci. 32, 582 (2004)] are used to assess the viability of the presented approach to calculate the ion drag force experienced by grains exposed to ion beams of well-defined energy. Excellent agreement between calculations and measurements is obtained for beam energies >10 eV, and the overprediction below 10 eV is attributed to the neglect of charge exchange collisions in our modeling. Along with critical assessments of our approach, suggestions for future experimental design to probe charging of and momentum transfer onto grains that capture the effect of space charge concentration and external fields are outlined.  more » « less
Award ID(s):
1903432
PAR ID:
10508794
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
AIP
Date Published:
Journal Name:
Physics of Plasmas
Volume:
30
Issue:
12
ISSN:
1070-664X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. The comparisons reveal that the ion/electron concentration and gas temperature in the afterglow critically influence the particle charge and the predictions are generally in qualitative agreement with the measurements. Along with critical assessment of the modeling assumptions, several recommendations are presented for future experimental design to probe charging in afterglows. 
    more » « less
  2. Based on the prior work of Chahl and Gopalakrishnan (2019) to infer particle-ion collision time distributions using a Langevin Dynamics (LD) approach, we develop a model for the non-dimensional diffusion charging collision kernel β_i or H that is applicable for 0≤Ψ_E≤60,0≤Ψ_I/Ψ_E ≤1,Kn_D≤2000 (defined in the main text). The developed model for β_i for attractive Coulomb and image potential interactions, along with the model for β_i for repulsive Coulomb and image potential interactions from Gopalakrishnan et al. (2013b), is tested against published diffusion charging experimental data. Current state of the art charging models, Fuchs (1963) and Wiedensohler (1988) regression for bipolar charging, are also evaluated and discussed. Comparisons reveal that the LD-based model accurately describes unipolar fractions for 10 – 100 nm particles measured in air (Adachi et al., 1985), nitrogen and argon but not in helium (Adachi et al., 1987). Fuchs model and the LD-based model yield similar predictions in the experimental conditions considered, except in helium. In the case of bipolar charging, the LD-based model captures the experimental trends quantitatively (within ±20%) across the entire size range of 4 – 40 nm producing superior agreement than Wiedensohler’s regression. The latter systematically underpredicts charge fraction below ~20 nm in air (by up to 40%) for the data presented in Adachi et al. (1985). Comparison with the data of Gopalakrishnan et al. (2015), obtained in UHP air along with measurements of the entire ion mass-mobility distribution, shows excellent agreement with the predictions of the LD-based model. This demonstrates the capability to accommodate arbitrary ion populations in any background gas, when such data is available. Wiedensohler’s regression, derived for bipolar charging in air using average ion mass-mobility, also describes the data reasonably well in the conditions examined. However, both models failed to capture the fraction of singly and doubly charged particles in carbon dioxide warranting further investigation. 
    more » « less
  3. Abstract Air‐sea drag governs the momentum transfer between the atmosphere and the ocean and remains largely unknown in hurricane winds. We revisit the momentum budget and eddy covariance methods to estimate the surface drag coefficient in the laboratory. Our drag estimates agree with field measurements in low‐to‐moderate winds and previous laboratory measurements in hurricane‐force winds. The drag coefficient saturates at 2.6×10−3andU10≈25 m s−1, in agreement with previous laboratory results by Takagaki et al. (2012,). During our analysis, we discovered an error in the original source code used by Donelan et al. (2004,). We present the corrected data and describe the correction procedure. Although the correction to the data does not change the key finding of drag saturation in strong winds, its magnitude and wind speed threshold are significantly changed. Our findings emphasize the need for an updated and unified drag parameterization based on field and laboratory data. 
    more » « less
  4. A molecular dynamics simulation of ion flow past dust grains is used to investigate the interaction between a pair of charged dust particles and streaming ions. The charging and dynamics of the grains are coupled and derived from the ion–dust interactions, allowing for detailed analysis of the ion wakefield structure and wakefield-mediated interaction as the dust particles change position. When a downstream grain oscillates vertically within the wake, it decharges by up to 30% as it approaches the upstream grain and then recharges as it recedes. There is an apparent hysteresis in charging depending on whether the grain is approaching or receding from a region of higher ion density. Maps of the ion-mediated dust–dust interaction force show that the radial extent of the wake region, which provides an attractive restoring force on the downstream particle, increases as the ion flow velocity decreases, though the restoring effect becomes weaker. As also shown in recent numerical results, there is no net attractive vertical force between the two grains. Instead, the reduced ion drag on the downstream particle allows it to “draft” in the wakefield of the upstream particle. 
    more » « less
  5. Abstract Architecting grain crystallographic orientation can modulate charge distribution and chemomechanical properties for enhancing the performance of polycrystalline battery materials. However, probing the interplay between charge distribution, grain crystallographic orientation, and performance remains a daunting challenge. Herein, we elucidate the spatially resolved charge distribution in lithium layered oxides with different grain crystallographic arrangements and establish a model to quantify their charge distributions. While the holistic “surface-to-bulk” charge distribution prevails in polycrystalline particles, the crystallographic orientation-guided redox reaction governs the charge distribution in the local charged nanodomains. Compared to the randomly oriented grains, the radially aligned grains exhibit a lower cell polarization and higher capacity retention upon battery cycling. The radially aligned grains create less tortuous lithium ion pathways, thus improving the charge homogeneity as statistically quantified from over 20 million nanodomains in polycrystalline particles. This study provides an improved understanding of the charge distribution and chemomechanical properties of polycrystalline battery materials. 
    more » « less