-
Steady oscillations in aggregation-fragmentation processes
Authors:
N. V. Brilliantov,
W. Otieno,
S. A. Matveev,
A. P. Smirnov,
E. E. Tyrtyshnikov,
P. L. Krapivsky
Abstract:
We report surprising steady oscillations in aggregation-fragmentation processes. Oscillating solutions are observed for the class of aggregation kernels $K_{i,j} = i^νj^μ + j^νi^μ$ homogeneous in masses $i$ and $j$ of merging clusters and fragmentation kernels, $F_{ij}=λK_{ij}$, with parameter $λ$ quantifying the intensity of the disruptive impacts. We assume a complete decomposition (shattering)…
▽ More
We report surprising steady oscillations in aggregation-fragmentation processes. Oscillating solutions are observed for the class of aggregation kernels $K_{i,j} = i^νj^μ + j^νi^μ$ homogeneous in masses $i$ and $j$ of merging clusters and fragmentation kernels, $F_{ij}=λK_{ij}$, with parameter $λ$ quantifying the intensity of the disruptive impacts. We assume a complete decomposition (shattering) of colliding partners into monomers. We show that an assumption of a steady-state distribution of cluster sizes, compatible with governing equations, yields a power law with an exponential cutoff. This prediction agrees with simulations results when $ θ\equiv ν-μ<1$. For $ θ=ν-μ>1$, however, the densities exhibit an oscillatory behavior. While these oscillations decay for not very small $λ$, they become steady if $θ$ is close to two and $λ$ is very small. Simulation results lead to a conjecture that for $ θ<1$ the system has a stable fixed point, corresponding to the steady-state density distribution, while for any $θ>1 $ there exists a critical value $λ_c(θ)$, such that for $λ< λ_c(θ)$, the system has an attracting limit cycle. This is rather striking for a closed system of Smoluchowski-like equations, lacking any sinks and sources of mass.
△ Less
Submitted 17 July, 2023;
originally announced July 2023.
-
Aggregation kinetics at sedimentation: the impact of particles diffusion
Authors:
Rishat R. Zagidullin,
Alexander P. Smirnov,
Sergey A. Matveev,
Nikolai V. Brilliantov
Abstract:
We investigate the aggregation kinetics of sedimenting particles theoretically and numerically, using the advection-diffusion equation. Agglomeration, caused by both transport mechanisms (diffusion and advection), is important for small particles, like primary ash or soot particles in atmosphere, and large particles of equal or close size, where the advection mechanism is weak. For small Peclet nu…
▽ More
We investigate the aggregation kinetics of sedimenting particles theoretically and numerically, using the advection-diffusion equation. Agglomeration, caused by both transport mechanisms (diffusion and advection), is important for small particles, like primary ash or soot particles in atmosphere, and large particles of equal or close size, where the advection mechanism is weak. For small Peclet numbers, which quantify the relative importance of diffusion and advection, we obtain the aggregation rates, as an expansion in Peclet numbers. For large Peclet numbers we use purely ballistic aggregation rates. Combining these results we obtain the rational approximant for the whole range of Peclet numbers. We also compute the aggregation rates by numerically solving the advection-diffusion equation. The results of the numerical simulations are in excellent agreement with the analytical theory for the studied Peclet numbers, varying by four orders of magnitude.
△ Less
Submitted 3 November, 2022; v1 submitted 29 December, 2021;
originally announced December 2021.
-
Aggregation in non-uniform systems with advection and localized source
Authors:
Rishat R. Zagidullin,
Alexander P. Smirnov,
Sergey A. Matveev,
Nikolai V. Brilliantov,
P. L. Krapivsky
Abstract:
We explore analytically and numerically agglomeration driven by advection and localized source. The system is inhomogeneous in one dimension, viz. along the direction of advection. We analyze a simplified model with mass-independent advection velocity, diffusion coefficient, and reaction rates. We also examine a model with mass-dependent coefficients describing aggregation with sedimentation. For…
▽ More
We explore analytically and numerically agglomeration driven by advection and localized source. The system is inhomogeneous in one dimension, viz. along the direction of advection. We analyze a simplified model with mass-independent advection velocity, diffusion coefficient, and reaction rates. We also examine a model with mass-dependent coefficients describing aggregation with sedimentation. For the simplified model, we obtain an exact solution for the stationary spatially dependent agglomerate densities. In the model describing aggregation with sedimentation, we report a new conservation law and develop a scaling theory for the densities. For numerical efficiency we exploit the low-rank approximation technique; this dramatically increases the computational speed and allows simulations of large systems. The numerical results are in excellent agreement with the predictions of our theory.
△ Less
Submitted 28 December, 2021; v1 submitted 25 December, 2021;
originally announced December 2021.
-
Oscillations in aggregation-shattering processes
Authors:
S. A. Matveev,
P. L. Krapivsky,
A. P. Smirnov,
E. E. Tyrtyshnikov,
N. V. Brilliantov
Abstract:
We observe never-ending oscillations in systems undergoing aggregation and collision-controlled shattering. Specifically, we investigate aggregation-shattering processes with aggregation kernels K_{i,j} = (i/j)^a+(j/i)^a and shattering kernels F_{i,j}=λK_{i,j}, where i and j are cluster sizes and parameter λquantifies the strength of shattering. When 0<a<1/2, there are no oscillations and the syst…
▽ More
We observe never-ending oscillations in systems undergoing aggregation and collision-controlled shattering. Specifically, we investigate aggregation-shattering processes with aggregation kernels K_{i,j} = (i/j)^a+(j/i)^a and shattering kernels F_{i,j}=λK_{i,j}, where i and j are cluster sizes and parameter λquantifies the strength of shattering. When 0<a<1/2, there are no oscillations and the system monotonically approaches to a steady state for all values of λ; in this region we obtain an analytical solution for the stationary cluster size distribution. Numerical solutions of the rate equations show that oscillations emerge in the 1/2<a<1 range. When the shattering rate is sufficiently large oscillations decay and eventually disappear, while for λ<λ_c(a) oscillations apparently persist forever. Thus never-ending oscillations can arise in closed aggregation-shattering processes without sinks and sources of particles.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
Low-rank approximation in the numerical modeling of the Farley-Buneman instability in ionospheric plasma
Authors:
S. V. Dolgov,
A. P. Smirnov,
E. E. Tyrtyshnikov
Abstract:
We consider the numerical modeling of the Farley-Buneman instability development in the earth's ionosphere plasma. The ion behavior is governed by the kinetic Landau equation in the four-dimensional phase space, and since the finite difference discretization on a tensor product grid is used, this equation becomes the most computationally challenging part of the scheme. To relax the complexity and…
▽ More
We consider the numerical modeling of the Farley-Buneman instability development in the earth's ionosphere plasma. The ion behavior is governed by the kinetic Landau equation in the four-dimensional phase space, and since the finite difference discretization on a tensor product grid is used, this equation becomes the most computationally challenging part of the scheme. To relax the complexity and memory consumption, an adaptive model reduction using the low-rank separation of variables, namely the Tensor Train format, is employed.
The approach was verified via the prototype MATLAB implementation. Numerical experiments demonstrate the possibility of efficient separation of space and velocity variables, resulting in the solution storage reduction by a factor of order tens.
△ Less
Submitted 27 August, 2013;
originally announced August 2013.