Tensor train versus Monte Carlo for the multicomponent Smoluchowski coagulation equation

作者: Sergey A. Matveev , Dmitry A. Zheltkov , Eugene E. Tyrtyshnikov , Alexander P. Smirnov

DOI: 10.1016/J.JCP.2016.04.025

关键词: Linear algebraApplied mathematicsMonte Carlo methodMathematicsGridTrapezoidal rule (differential equations)Smoluchowski coagulation equationCurse of dimensionalityKernel (statistics)ComputationMathematical analysisPhysics and Astronomy (miscellaneous)Computer Science Applications

摘要: In this paper we present a novel numerical algorithm for the space-homogeneous multicomponent (multidimensional) Smoluchowski coagulation equation, number of components is considered as dimensionality. The new methodology based on classical finite-difference predictor-corrector scheme. straightforward implementation scheme, however, one would have to compute and store prohibitively many values grid function at nodes multidimensional grid. We propose use special low-parametric representations functions well kernel. corresponding arrays are approximated by low-rank tensor-train decompositions reducing them combinations small low-dimensional arrays, eventually matrices which can fast algorithms linear algebra. Instead O ( N 2 d ) operations in method that requires only log ? operations, where per axis space components. work accelerate time-scheme trapezoidal rule computation integral operators. Thus, accuracy h + , step time step.

参考文章(28)
Eli Ben-Naim, Sidney Redner, Pavel L. Krapivsky, A Kinetic View of Statistical Physics ,(2010)
Dietrich Braess, Wolfgang Hackbusch, On the efficient computation of high-dimensional integrals and the approximation by exponential sums Springer, Berlin, Heidelberg. pp. 39- 74 ,(2009) , 10.1007/978-3-642-03413-8_3
Geethpriya Palaniswaamy, Sudarshan K. Loyalka, Direct simulation monte carlo aerosol dynamics : Coagulation and collisional sampling Nuclear Technology. ,vol. 156, pp. 29- 38 ,(2006) , 10.13182/NT06-A3771
A.A Lushnikov, Evolution of coagulating systems: III. Coagulating mixtures Journal of Colloid and Interface Science. ,vol. 54, pp. 94- 101 ,(1976) , 10.1016/0021-9797(76)90288-5
Kok Foong Lee, Robert I.A. Patterson, Wolfgang Wagner, Markus Kraft, Stochastic weighted particle methods for population balance equations with coagulation, fragmentation and spatial inhomogeneity Journal of Computational Physics. ,vol. 303, pp. 1- 18 ,(2015) , 10.1016/J.JCP.2015.09.031
Haibo Zhao, F. Einar Kruis, Chuguang Zheng, A differentially weighted Monte Carlo method for two-component coagulation Journal of Computational Physics. ,vol. 229, pp. 6931- 6945 ,(2010) , 10.1016/J.JCP.2010.05.031
Robin C. Ball, Colm Connaughton, Peter P. Jones, R. Rajesh, Oleg Zaboronski, Collective oscillations in irreversible coagulation driven by monomer inputs and large-cluster outputs. Physical Review Letters. ,vol. 109, pp. 168304- ,(2012) , 10.1103/PHYSREVLETT.109.168304
I. V. Oseledets, Tensor-Train Decomposition SIAM Journal on Scientific Computing. ,vol. 33, pp. 2295- 2317 ,(2011) , 10.1137/090752286
Shraddha Shekar, William J. Menz, Alastair J. Smith, Markus Kraft, Wolfgang Wagner, On a multivariate population balance model to describe the structure and composition of silica nanoparticles Computers & Chemical Engineering. ,vol. 43, pp. 130- 147 ,(2012) , 10.1016/J.COMPCHEMENG.2012.04.010
Vladimir A. Kazeev, Boris N. Khoromskij, Eugene E. Tyrtyshnikov, Multilevel Toeplitz Matrices Generated by Tensor-Structured Vectors and Convolution with Logarithmic Complexity SIAM Journal on Scientific Computing. ,vol. 35, ,(2013) , 10.1137/110844830