A numerical study on large-time asymptotics of the Lifshitz-Slyozov system

We numerically investigate the long time behavior of solutions of the Lifshitz-Slyozov system. We propose a numerical scheme in which the numerical dissipation is controlled in such a way that the results for large time are meaningful. In this respect, we find the long time behavior to crucially depend on the distribution of largest aggregates present in the solution. This fact proved, in some particular cases in [23], was difficult to obtain with previous numerical schemes in the engineering literature leading to wrong statements. We propose a numerical scheme in which we can observe and quantify the equilibration rates towards the right asymptotic profile. Moreover, this system appears to be a very interesting test problem for any anti dissipative scheme for conservation laws.