An official website of the United States government
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
Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling
The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Potential applications and caveats to the usage of LD are included as a summary.
more »
« less
- Award ID(s):
- 1903432
- PAR ID:
- 10221107
- Date Published:
- Journal Name:
- Journal of aerosol science
- Volume:
- 155
- Issue:
- June
- ISSN:
- 0021-8502
- Page Range / eLocation ID:
- 105746
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Droplet size distributions in turbulent clouds: experimental evaluation of theoretical distributionsAbstract Precipitation efficiency and optical properties of clouds, both central to determining Earth's weather and climate, depend on the size distribution of cloud particles. In this work theoretical expressions for cloud droplet size distribution shape are evaluated using measurements from controlled experiments in a convective‐cloud chamber. The experiments are a unique opportunity to constrain theory because they are in steady‐state and because the initial and boundary conditions are well characterized compared to typical atmospheric measurements. Three theoretical distributions obtained from a Langevin drift‐diffusion approach to cloud formation via stochastic condensation are tested: (a) stochastic condensation with a constant removal time‐scale; (b) stochastic condensation with a size‐dependent removal time‐scale; (c) droplet growth in a fixed supersaturation condition and with size‐dependent removal. In addition, a similar Weibull distribution that can be obtained from the drift‐diffusion approach, as well as from mechanism‐independent probabilistic arguments (e.g., maximum entropy), is tested as a fourth hypothesis. Statistical techniques such as theχ2test, sum of squared errors of prediction, and residual analysis are employed to judge relative success or failure of the theoretical distributions to describe the experimental data. An extensive set of cloud droplet size distributions are measured under different aerosol injection rates. Five different aerosol injection rates are run both for size‐selected aerosol particles, and six aerosol injection rates are run for broad‐distribution, polydisperse aerosol particles. In relative comparison, the most favourable comparison to the measurements is the expression for stochastic condensation with size‐dependent droplet removal rate. However, even this optimal distribution breaks down for broad aerosol size distributions, primarily due to deviations from the measured large‐droplet tail. A possible explanation for the deviation is the Ostwald ripening effect coupled with deactivation/activation in polluted cloud conditions.more » « less
-
Kim, Y.; Moon, D.H. (Ed.)By using gravity/gauge correspondence, we construct a holographic model, constrained to mimic the lattice QCD equation of state at zero density, to investigate the temperature and baryon chemical potential dependence of the equation of state. We also obtained the energy loss of light and heavy partons within the hot and dense plasma represented by the heavy quark drag force, Langevin diffusion coefficients and jet quenching parameter at the critical point and across the first-order transition line predicted by the model.more » « less
-
Meier-Schellersheim, Martin (Ed.)We introduce a Stochastic Reaction-Diffusion-Dynamics Model (SRDDM) for simulations of cellular mechanochemical processes with high spatial and temporal resolution. The SRDDM is mapped into the CellDynaMo package, which couples the spatially inhomogeneous reaction-diffusion master equation to account for biochemical reactions and molecular transport within the Langevin Dynamics (LD) framework to describe dynamic mechanical processes. This computational infrastructure allows the simulation of hours of molecular machine dynamics in reasonable wall-clock time. We apply SRDDM to test performance of the Search-and-Capture of mitotic spindle assembly by simulating, in three spatial dimensions, dynamic instability of elastic microtubules anchored in two centrosomes, movement and deformations of geometrically realistic centromeres with flexible kinetochores and chromosome arms. Furthermore, the SRDDM describes the mechanics and kinetics of Ndc80 linkers mediating transient attachments of microtubules to the chromosomal kinetochores. The rates of these attachments and detachments depend upon phosphorylation states of the Ndc80 linkers, which are regulated in the model by explicitly accounting for the reactions of Aurora A and B kinase enzymes undergoing restricted diffusion. We find that there is an optimal rate of microtubule-kinetochore detachments which maximizes the accuracy of the chromosome connections, that adding chromosome arms to kinetochores improve the accuracy by slowing down chromosome movements, that Aurora A and kinetochore deformations have a small positive effect on the attachment accuracy, and that thermal fluctuations of the microtubules increase the rates of kinetochore capture and also improve the accuracy of spindle assembly.more » « less
-
Mucus, composed significantly of glycosylated mucins, is a soft and rheologically complex material that lines respiratory, reproductive, and gastrointestinal tracts in mammals. Mucus may present as a gel, as a highly viscous fluid, or as a viscoelastic fluid. Mucus acts as a barrier to the transport of harmful microbes and inhaled atmospheric pollutants to underlying cellular tissue. Studies on mucin gels have provided critical insights into the chemistry of the gels, their swelling kinetics, and the diffusion and permeability of molecular constituents such as water. The transport and dispersion of micron and sub-micron particles in mucin gels and solutions, however, differs from the motion of small molecules since the much larger tracers may interact with microstructure of the mucin network. Here, using brightfield and fluorescence microscopy, high-speed particle tracking, and passive microrheology, we study the thermally driven stochastic movement of 0.5–5.0 μm tracer particles in 10% mucin solutions at neutral pH, and in 10% mucin mixed with industrially relevant dust; specifically, unmodified limestone rock dust, modified limestone, and crystalline silica. Particle trajectories are used to calculate mean square displacements and the displacement probability distributions; these are then used to assess tracer diffusion and transport. Complex moduli are concomitantly extracted using established microrheology techniques. We find that under the conditions analyzed, the reconstituted mucin behaves as a weak viscoelastic fluid rather than as a viscoelastic gel. For small- to moderately sized tracers with a diameter of lessthan 2 μm, we find that effective diffusion coefficients follow the classical Stokes–Einstein relationship. Tracer diffusivity in dust-laden mucin is surprisingly larger than in bare mucin. Probability distributions of mean squared displacements suggest that heterogeneity, transient trapping, and electrostatic interactions impact dispersion and overall transport, especially for larger tracers. Our results motivate further exploration of physiochemical and rheological mechanisms mediating particle transport in mucin solutions and gels.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/https://doi.org/10.1016/j.jaerosci.2021.105746
