[Paper Review] Dynamical scaling in Smoluchowski's coagulation equations: uniform convergence

Govind Menon, Robert L. Pego|ArXiv.org|Jun 24, 2003

Coagulation and Flocculation Studies14 references4 citations

TL;DR

This paper establishes uniform convergence of solutions to self-similar forms in Smoluchowski's coagulation equations for solvable kernels $K=2$, $x+y$, and $xy$, using Fourier inversion and Laplace transform methods with complex characteristics. It proves uniform convergence of densities to exponential-tailed self-similar solutions under optimal moment and regularity conditions, extending classical probabilistic limit theorems to coagulation dynamics.

ABSTRACT

We consider the approach to self-similarity (or dynamical scaling) in Smoluchowski's coagulation equations for the solvable kernels K(x,y)=2, x+y and xy. We prove the uniform convergence of densities to the self-similar solution with exponential tails under the regularity hypothesis that a suitable moment have an integrable Fourier transform. For the discrete equations we prove uniform convergence under optimal moment hypotheses. Our results are completely analogous to classical local convergence theorems for the normal law in probability theory. The proofs rely on the Fourier inversion formula and the solution by the method of characteristics for the Laplace transform.

Motivation & Objective

• To rigorously establish uniform convergence of solution densities to self-similar forms in Smoluchowski’s coagulation equations for solvable kernels.
• To extend classical probabilistic limit theorems—such as the central limit theorem—to the dynamics of coagulation via analytical methods.
• To weaken regularity and decay assumptions in prior results, achieving optimal moment hypotheses for convergence.
• To unify the treatment of continuous and discrete equations under a common framework using similarity variables and rescaling.
• To provide a rigorous foundation for the dynamical scaling hypothesis widely used in physics and chemistry, especially near gelation times.

Proposed method

• Utilizes the Fourier inversion formula to analyze convergence in the frequency domain, enabling uniform bounds on solution densities.
• Applies the method of characteristics in the right half-plane to derive strong decay estimates for the Laplace transform of the solution.
• Employs exact solution formulas for the Laplace transform of the density, derived from the method of characteristics.
• Uses a change of variables to map the multiplicative kernel $K=xy$ to the additive kernel $K=x+y$, reducing the problem to a known case.
• Defines similarity variables that rescale time and mass to isolate self-similar behavior near gelation or long-time limits.
• Applies uniform convergence criteria via integrability of the Fourier transform of $x^3 n_0(x)$, ensuring convergence in sup-norm.

Experimental results

Research questions

• RQ1 Under what conditions does the solution density of Smoluchowski’s coagulation equation converge uniformly to a self-similar profile?
• RQ2 How can classical probabilistic limit theorems be adapted to prove uniform convergence in coagulation dynamics?
• RQ3 What is the minimal regularity and moment condition required for uniform convergence to self-similar solutions with exponential tails?
• RQ4 How does the convergence behavior differ for the kernels $K=2$, $x+y$, and $xy$, especially near gelation?
• RQ5 Can the convergence to self-similar form for the $K=xy$ kernel be rigorously established as $t \to T_{\rm gel}$ using rescaling and transform methods?

Key findings

• For $K=2$ and $x+y$, the paper proves uniform convergence of solution densities to self-similar solutions with exponential tails under the assumption that the Fourier transform of $x^3 n_0(x)$ is integrable.
• The result for $K=2$ improves upon Kreer and Penrose by weakening the decay and regularity assumptions on initial data to almost optimal moment and continuity conditions.
• For $K=xy$, uniform convergence of rescaled densities to the self-similar profile $\hat{n}_{*,2}(\hat{x}) = \frac{1}{\sqrt{2\pi \hat{x}^5}} e^{-\hat{x}/2}$ is established as $t \to 1^-$, the gelation time.
• The rescaling $\hat{x} = (1-t)^2 x$, $\hat{n}(t,\hat{x}) = (1-t)^{-5} n(t,x)$ preserves the second moment and captures the divergent mass flux near gelation.
• In the discrete case, uniform convergence is proven via the relation $h l n_l^{\rm mul}(t) = (1-t)^{-1} n_l^{\rm add}(\log(1-t)^{-1})$, with similar sup-norm convergence in similarity variables.
• The self-similar profile for $K=xy$ is shown to emerge universally from initial data with finite second and third moments, confirming the scaling hypothesis in the limit.