skip to main content

Title: Dust charging in dynamic ion wakes
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.
; ; ; ; ;
Award ID(s):
1740203 1707215
Publication Date:
Journal Name:
Physics of plasmas
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Particle acceleration behind a shock wave due to interactions between magnetic islands in the heliosphere has attracted attention in recent years. The downstream acceleration may yield a continuous increase of particle flux downstream of the shock wave. Although it is not obvious how the downstream magnetic islands are produced, it has been suggested that current sheets are involved in the generation of magnetic islands due to their interaction with a shock wave. We perform 2D hybrid kinetic simulations to investigate the interaction between multiple current sheets and a shock wave. In the simulation, current sheets are compressed by the shock wave and a tearing instability develops at the compressed current sheets downstream of the shock. As the result of this instability, the electromagnetic fields become turbulent and magnetic islands form well downstream of the shock wave. We find a “post-cursor” region in which the downstream flow speed normal to the shock wave in the downstream rest frame is decelerated to ∼ 1 V A immediately behind the shock wave, where V A is the upstream Alfvén speed. The flow speed then gradually decelerates to 0 accompanied by the development of the tearing instability. We also observe an efficient productionmore »of energetic particles above 100 E 0 during the development of the instability some distance downstream of the shock wave, where E 0 = m p V A 2 and m p is the proton mass. This feature corresponds to Voyager observations showing that the anomalous cosmic-ray intensity increase begins some distance downstream of the heliospheric termination shock.« less
  2. We develop a physical and computational model for performing fully coupled, grain-resolved direct numerical simulations of cohesive sediment, based on the immersed boundary method. The model distributes the cohesive forces over a thin shell surrounding each particle, thereby allowing for the spatial and temporal resolution of the cohesive forces during particle–particle interactions. The influence of the cohesive forces is captured by a single dimensionless parameter in the form of a cohesion number, which represents the ratio of cohesive and gravitational forces acting on a particle. We test and validate the cohesive force model for binary particle interactions in the drafting–kissing–tumbling (DKT) configuration. Cohesive sediment grains can remain attached to each other during the tumbling phase following the initial collision, thereby giving rise to the formation of flocs. The DKT simulations demonstrate that cohesive particle pairs settle in a preferred orientation, with particles of very different sizes preferentially aligning themselves in the vertical direction, so that the smaller particle is drafted in the wake of the larger one. This preferred orientation of cohesive particle pairs is found to remain influential for systems of higher complexity. To this end, we perform large simulations of 1261 polydisperse settling particles starting from rest. Thesemore »simulations reproduce several earlier experimental observations by other authors, such as the accelerated settling of sand and silt particles due to particle bonding, the stratification of cohesive sediment deposits, and the consolidation process of the deposit. They identify three characteristic phases of the polydisperse settling process, viz. (i) initial stir-up phase with limited flocculation, (ii) enhanced settling phase characterized by increased flocculation, and (iii) consolidation phase. The simulations demonstrate that cohesive forces accelerate the overall settling process primarily because smaller grains attach to larger ones and settle in their wakes. For the present cohesive number values, we observe that settling can be accelerated by up to 29 %. We propose physically based parametrization of classical hindered settling functions introduced by earlier authors, in order to account for cohesive forces. An investigation of the energy budget shows that, even though the work of the collision forces is much smaller than that of the hydrodynamic drag forces, it can substantially modify the relevant energy conversion processes.« less
  3. Particle shape strongly influences the diffusion charging of aerosol particles exposed to bipolar/unipolar ions and accurate modeling is needed to predict the charge distribution of non-spherical particles. A prior particle-ion collision kernel β_i model including Coulombic and image potential interactions for spherical particles is generalized for arbitrary shapes following a scaling approach that uses a continuum and free molecular particle length scale and Langevin dynamics simulations of non-spherical particle-ion collisions for attractive Coulomb-image potential interactions. This extended β_i model for collisions between unlike charged particle-ion (bipolar charging) and like charged particle-ion (unipolar charging) is validated by comparing against published experimental data of bipolar charge distributions for diverse shapes. Comparison to the bipolar charging data for spherical particles shows good agreement in air, argon, and nitrogen, while also demonstrating high accuracy in predicting charge states up to ±6. Comparisons to the data for fractal aggregates reveal that the LD-based β_i model predicts within overall ±30% without any systematic bias. The mean charge on linear chain aggregates and charge fractions on cylindrical particles is found to be in good agreement with the measurements (~±20% overall). The comparison with experimental results supports the use of LD-based diffusion charging models to predict the bipolarmore »and unipolar charge distribution of arbitrary shaped aerosol particles for a wide range of particle size, and gas temperature, pressure. The presented β_i model is valid for perfectly conducting particles and in the absence of external electric fields; these simplifications need to be addressed in future work on particle charging.« less
  4. Abstract Nonlinear frequency response analysis is a widely used method for determining system dynamics in the presence of nonlinearities. In dusty plasmas, the plasma–grain interaction (e.g. grain charging fluctuations) can be characterized by a single-particle non-linear response analysis, while grain–grain non-linear interactions can be determined by a multi-particle non-linear response analysis. Here a machine learning-based method to determine the equation of motion in the non-linear response analysis for dust particles in plasmas is presented. Searching the parameter space in a Bayesian manner allows an efficient optimization of the parameters needed to match simulated non-linear response curves to experimentally measured non-linear response curves.
  5. The margination and adhesion of micro-particles (MPs) have been extensively investigated separately, due to their important applications in the biomedical field. However, the cascade process from margination to adhesion should play an important role in the transport of MPs in blood flow. To the best of our knowledge, this has not been explored in the past. Here we numerically study the margination behaviour of elastic MPs to blood vessel walls under the interplay of their deformability and adhesion to the vessel wall. We use the lattice Boltzmann method and molecular dynamics to solve the fluid dynamics and particle dynamics (including red blood cells (RBCs) and elastic MPs) in blood flow, respectively. Additionally, a stochastic ligand–receptor binding model is employed to capture the adhesion behaviours of elastic MPs on the vessel wall. Margination probability is used to quantify the localization of elastic MPs at the wall. Two dimensionless numbers are considered to govern the whole process: the capillary number $Ca$ , denoting the ratio of viscous force of fluid flow to elastic interfacial force of MP, and the adhesion number $Ad$ , representing the ratio of adhesion strength to viscous force of fluid flow. We systematically vary them numerically and amore »margination probability contour is obtained. We find that there exist two optimal regimes favouring high margination probability on the plane $Ca{-}Ad$ . The first regime, namely region I, is that with high adhesion strength and moderate particle stiffness; the other one, region II, has moderate adhesion strength and large particle stiffness. We conclude that the existence of optimal regimes is governed by the interplay of particle deformability and adhesion strength. The corresponding underlying mechanism is also discussed in detail. There are three major factors that contribute to the localization of MPs: (i) near-wall hydrodynamic collision between RBCs and MPs; (ii) deformation-induced migration due to the presence of the wall; and (iii) adhesive interaction between MPs and the wall. Mechanisms (i) and (iii) promote margination, while (ii) hampers margination. These three factors perform different roles and compete against each other when MPs are located in different regions of the flow channel, i.e. near-wall region. In optimal region I, adhesion outperforms deformation-induced migration; and in region II, the deformation-induced migration is small compared to the coupling of near-wall hydrodynamic collision and adhesion. The finding of optimal regimes can help the understanding of localization of elastic MPs at the wall under the adhesion effect in blood flow. More importantly, our results suggest that softer MP or stronger adhesion is not always the best choice for the localization of MPs.« less