Sergey A. Matveev,Ilya Tretyak
Abstract:We propose an efficient implementation of the numerical tensor-train (TT) based algorithm solving the multicomponent coagulation equation preserving the nonnegativeness of solution. Unnatural negative elements in the constructed approximation arise due to the errors of the low-rank decomposition and discretization scheme. In this work, we propose to apply the rank-one corrections in the TT-format proportional to the minimal negative element. Such an element can be found via application of the global optimization methods that can be fully implemented within efficient operations in the tensor train format. We incorporate this trick into the time-integration scheme for the multicomponent coagulation equation and also use it for post-processing of the stationary solution for the problem with the source of particles.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the non - negativity of the numerical solution of the multi - component Smoluchowski coagulation equation. Specifically, the author proposes an efficient algorithm based on the Tensor Train (TT) format to solve the multi - component coagulation equation and ensure the non - negativity of the solution. Due to the errors of low - rank decomposition and discretization schemes, unnatural negative elements may appear in the constructed approximate solution. To solve this problem, the author proposes to apply a rank - one correction in the TT format to eliminate these negative elements.
### Main problems
1. **Numerical solution of multi - component coagulation equation**:
- Coagulation/aggregation refers to the aggregation process of interacting particles and has wide applications in nature and technology, such as raindrop formation, rubber production, and nanoparticle synthesis.
- The basic Smoluchowski equation describes the evolution of the particle number density function \(n(v, t)\), where \(v\) represents the volume or mass of the particle and \(t\) is time. For multi - component systems, this equation can be generalized to a multi - dimensional form.
2. **Maintaining the non - negativity of the solution**:
- During the numerical solution process, due to the errors of low - rank decomposition and discretization schemes, negative values may appear in the solution, which is unreasonable because the particle concentration should be non - negative.
- The author proposes to solve this problem by applying a rank - one correction in the TT format, that is, adding a rank - one correction term proportional to the smallest negative element, thereby ensuring the non - negativity of the solution.
### Solutions
- **Rank - one correction**: Find the smallest negative element through a global optimization method and apply a rank - one correction in the TT format to ensure the non - negativity of the solution.
- **Time integration scheme**: Incorporate this correction method into the time integration scheme and use it to process the steady - state solution post - processing of the multi - component coagulation problem with a source term.
### Mathematical expressions
The basic form of the multi - component Smoluchowski equation is as follows:
\[
\frac{\partial n(v, t)}{\partial t}=\frac{1}{2} \int_{0}^{v_{1}} \cdots \int_{0}^{v_{d}} K(v - u; u) n(v - u, t) n(u, t) \, du_{1} \cdots du_{d}-n(v, t) \int_{0}^{\infty} \cdots \int_{0}^{\infty} K(u; v) n(u, t) \, du_{1} \cdots du_{d}+q(v, t)
\]
where \(n(v, t)\) is the particle concentration function, \(K(u, v)\) is the symmetric reaction rate kernel function, and \(q(v, t)\) is the particle source term.
### Conclusions
The author verifies the effectiveness of this method through numerical experiments and shows its applicability in different dimensions and different coagulation kernels. The results show that this method can effectively solve the multi - component coagulation equation while ensuring the non - negativity of the solution and has a relatively low computational complexity.
### References
- [1] P. L. Krapivsky, S. Redner, and E. Ben - Naim, *A Kinetic View of Statistical Physics*. Cambridge University Press, 2010.
- [2] F. Leyvraz, "Scaling theory and exactly solved models in the kinetics of irreversible aggregation," *Phys. Reports*, vol. 383, pp. 95–212, 2003.
- [3] G. Falkovich, A. Fouxon, and M. G. Stepanov, "Acceleration of rain initiation by cloud turbulence," *Nature*, vol. 41