Abstract:We study a class of discrete models in which a collection of particles evolves in time following the gradient flow of an energy depending on the cell areas of an associated Laguerre (i.e. a weighted Voronoi) tessellation. We consider the high number of cell limit of such systems and, using a modulated energy argument, we prove convergence towards smooth solutions of nonlinear diffusion PDEs of porous medium type.
What problem does this paper attempt to address?
### The problems the paper attempts to solve
The paper aims to study a class of discrete models in which the evolution of a set of particles over time follows an energy - gradient flow related to the cell areas of the associated Laguerre (weighted Voronoi) tessellation. Specifically, the paper focuses on the limiting behavior of these systems as the number of particles tends to infinity and uses the modulated energy argument method to prove the convergence of such discrete models to the smooth solutions of nonlinear diffusion partial differential equations (PDEs).
#### Main problem description
1. **Particle system and tessellation**:
- Given a compact domain \(\Omega\subset\mathbb{R}^d\) and a finite number \(N\) of particles, the position of each particle is parameterized by the vector \(X=(x_1,\dots,x_N)\in(\mathbb{R}^d)^N\).
- The interaction between particles is described by the Laguerre tessellation, which is a weighted Voronoi tessellation.
2. **Energy function**:
- The energy function \(F_{\varepsilon}(X)\) is defined as:
\[
F_{\varepsilon}(X):=\inf_{L\in L_N(\Omega)}\sum_i\int_{L_i}\frac{\vert x - x_i\vert^2}{2\varepsilon}\,dx+\sum_iC_i(\vert L_i\vert)
\]
where \(C_i(a)\) is a function depending on the cell volume, usually adopting a nonlinear spring model:
\[
C_i(a)=K(a)\frac{\vert a-\bar{a}\vert^2}{2}
\]
Here \(K(a)\) is the bulk modulus and \(\bar{a}\) is the target volume.
3. **Particle dynamics**:
- The dynamics of the particle positions is determined by the gradient flow of the energy \(F_{\varepsilon}\), specifically:
\[
\dot{X}(t)=-\nabla_{m_0}F_{\varepsilon}(X(t))
\]
where \(\nabla_{m_0}\) represents the gradient with respect to the weighted \(l_2\) metric.
4. **Continuous limit**:
- The main result of the paper is to prove that the discrete solution \(\bar{\mu}_N\) converges to the smooth solution \(\rho\) of the nonlinear diffusion equation (such as the porous medium equation). Specifically, for a given initial condition \(X_0\) and appropriate parameters \(\delta_N\) and \(\varepsilon\), the error estimate is:
\[
\max_{t\in[0,T]}\int_{\Omega}U(\bar{\mu}_N(t,\cdot)\vert\rho(t,\cdot))\leq C\left(\frac{\delta_N^2}{\varepsilon}+\varepsilon^{p - 1}\right)
\]
where \(p=\min\{2,\alpha\}\), \(U(r\vert s)=U(r)-U(s)-U'(s)(r - s)\).
### Summary
The core problem of the paper is to study the behavior of the discrete particle system in the high - particle - number limit and prove its convergence to the solutions of the continuous nonlinear diffusion equation. By introducing an appropriate energy function and the modulated energy method, the author successfully builds a bridge between the discrete model and the continuous model and provides a rigorous mathematical proof.