Stationary mass distribution and nonlocality in models of coalescence and shattering

Pattern formation by turbulent cascades

Fully developed turbulence is a universal and scale-invariant chaotic state characterized by an energy cascade from large to small scales at which the cascade is eventually arrested by dissipation1-6. Here we show how to harness these seemingly structureless turbulent cascades to generate patterns. Pattern formation entails a process of wavelength selection, which can usually be traced to the linear instability of a homogeneous state7. By contrast, the mechanism we propose here is fully nonlinear. It is triggered by the non-dissipative arrest of turbulent cascades: energy piles up at an intermediate scale, which is neither the system size nor the smallest scales at which energy is usually dissipated. Using a combination of theory and large-scale simulations, we show that the tunable wavelength of these cascade-induced patterns can be set by a non-dissipative transport coefficient called odd viscosity, ubiquitous in chiral fluids ranging from bioactive to quantum systems8-12. Odd viscosity, which acts as a scale-dependent Coriolis-like force, leads to a two-dimensionalization of the flow at small scales, in contrast with rotating fluids in which a two-dimensionalization occurs at large scales4. Apart from odd viscosity fluids, we discuss how cascade-induced patterns can arise in natural systems, including atmospheric flows13-19, stellar plasma such as the solar wind20-22, or the pulverization and coagulation of objects or droplets in which mass rather than energy cascades23-25.

Lattice models for ballistic aggregation: Cluster-shape-dependent exponents

We study ballistic aggregation on a two-dimensional square lattice, where particles move ballistically in between momentum and mass conserving coalescing collisions. Three models are studied based on the shapes of the aggregates: In the first the aggregates remain point particles, in the second they retain the fractal shape at the time of collision, and in the third they assume a spherical shape. The exponents describing the power-law temporal decay of number of particles and energy as well as dependence of velocity correlations on mass are determined using large-scale Monte Carlo simulations. It is shown that the exponents are universal only for the point-particle model. In the other two cases, the exponents are dependent on the initial number density and correlations vanish at high number densities. The fractal dimension for the second model is close to 1.49.

Anomalous aggregation regimes of temperature-dependent Smoluchowski equations

Temperature-dependent Smoluchowski equations describe the ballistic agglomeration. In contrast to the standard Smoluchowski equations for the evolution of cluster densities, with constant rate coefficients, the temperature-dependent equations describe both-the evolution of the densities as well as cluster temperatures, which determine the agglomeration rates. To solve these equations, we develop a Monte Carlo technique based on the low-rank approximation for the aggregation kernel. Using this highly effective approach, we perform a comprehensive study of the kinetic phase diagram of the system and reveal a few surprising regimes, including permanent temperature growth and "density separation" regime, with a large gap in the size distribution for middle-size clusters. We perform scaling analysis and classify the aggregation kernels for the temperature-dependent equations. Furthermore, we conjecture the lack of gelation in such systems. The results of the scaling theory agree well with the simulation data.

Stochastic Theory of Discrete Binary Fragmentation—Kinetics and Thermodynamics

We formulate binary fragmentation as a discrete stochastic process in which an integer mass k splits into two integer fragments j, k−j, with rate proportional to the fragmentation kernel Fj,k−j. We construct the ensemble of all distributions that can form in fixed number of steps from initial mass M and obtain their probabilities in terms of the fragmentation kernel. We obtain its partition function, the mean distribution and its evolution in time, and determine its stability using standard thermodynamic tools. We show that shattering is a phase transition that takes place when the stability conditions of the partition function are violated. We further discuss the close analogy between shattering and gelation, and between fragmentation and aggregation in general.

Nonextensive Supercluster States in Aggregation with Fragmentation

Systems evolving through aggregation and fragmentation may possess an intriguing supercluster state (SCS). Clusters constituting this state are mostly very large, so the SCS resembles a gelling state, but the formation of the SCS is controlled by fluctuations and in this aspect, it is similar to a critical state. The SCS is nonextensive, that is, the number of clusters varies sublinearly with the system size. In the parameter space, the SCS separates equilibrium and jamming (extensive) states. The conventional methods, such as, e.g., the van Kampen expansion, fail to describe the SCS. To characterize the SCS we propose a scaling approach with a set of critical exponents. Our theoretical findings are in good agreement with numerical results.

Hopf bifurcation in addition-shattering kinetics

In aggregation-fragmentation processes, a steady state is usually reached. This indicates the existence of an attractive fixed point in the underlying infinite system of coupled ordinary differential equations. The next simplest possibility is an asymptotically periodic motion. Never-ending oscillations have not been rigorously established so far, although oscillations have been recently numerically detected in a few systems. For a class of addition-shattering processes, we provide convincing numerical evidence for never-ending oscillations in a certain region U of the parameter space. The processes which we investigate admit a fixed point that becomes unstable when parameters belong to U and never-ending oscillations effectively emerge through a Hopf bifurcation.

Statistical Mechanics of Discrete Multicomponent Fragmentation

We formulate the statistics of the discrete multicomponent fragmentation event using a methodology borrowed from statistical mechanics. We generate the ensemble of all feasible distributions that can be formed when a single integer multicomponent mass is broken into fixed number of fragments and calculate the combinatorial multiplicity of all distributions in the set. We define random fragmentation by the condition that the probability of distribution be proportional to its multiplicity, and obtain the partition function and the mean distribution in closed form. We then introduce a functional that biases the probability of distribution to produce in a systematic manner fragment distributions that deviate to any arbitrary degree from the random case. We corroborate the results of the theory by Monte Carlo simulation, and demonstrate examples in which components in sieve cuts of the fragment distribution undergo preferential mixing or segregation relative to the parent particle.

Role of energy in ballistic agglomeration

We study a ballistic agglomeration process in the reaction-controlled limit. Cluster densities obey an infinite set of Smoluchowski rate equations, with rates dependent on the average particle energy. The latter is the same for all cluster species in the reaction-controlled limit and obeys an equation depending on densities. We express the average energy through the total cluster density that allows us to reduce the governing equations to the standard Smoluchowski equations. We derive basic asymptotic behaviors and verify them numerically. We also apply our formalism to the agglomeration of dark matter.

Steady oscillations in aggregation-fragmentation processes

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 simulation results when θ≡ν-μ<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 2 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.