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. These 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.
more »
« less
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.
more »
« less
- NSF-PAR ID:
- 10145806
- Date Published:
- Journal Name:
- Physics of plasmas
- Volume:
- 27
- Issue:
- 2
- ISSN:
- 1089-7674
- Page Range / eLocation ID:
- 023703
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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 production 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.more » « less
-
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.more » « less
-
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 bipolar 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.more » « less
-
The microgravity environment of the Plasmakristall-4 experiment on the International Space Station provides a laboratory for exploring plasma-mediated interactions among charged dust grains in fully three-dimensional space. Away from the strong influence of Earth's gravity, the dust grains can levitate in the bulk of the plasma, where they have been observed to form extended filamentary structures aligned with the discharge tube axis. These structures can be used as a macroscopic analogue for other self-organizing systems, including electrorheological fluids and liquid crystals, and the success of the analogy depends on a thorough understanding of the mechanisms guiding the dust interaction potential. Here we present the results from molecular dynamics simulations of the ion flow past isolated dust chains within the dust cloud and the dust cloud macrostructure. Although dust grains are known to respond on the millisecond timescale, analysis reveals that periodic variations of plasma conditions on the microsecond timescale significantly affect dust structure formation. In addition to the expected formation of filamentary dust chains in the dust cloud macrostructure, dust grains in a large cloud are also observed to organize into ordered positions on the surface of nested cylinders, in agreement with experimental observations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/https://doi.org/10.1063/1.5124246
