Large friction limit of the almost pressureless Euler-Poisson system

1. Problem Statement and Motivation

The System:
The paper investigates the almost pressureless Euler-Poisson (EP) system with repulsive force and large friction. The system models the dynamics of a plasma (charged particles) in an electric field generated by a constant background charge $M$. The governing equations (in a rescaled time variable $\tau = t/\varepsilon$) are:
$$\begin{aligned} \partial_\tau \rho + \frac{1}{\varepsilon} \text{div}(\rho u) &= 0, \\ \partial_\tau (\rho u) + \frac{1}{\varepsilon} \text{div}(\rho u \otimes u) + \frac{1}{\varepsilon^{1-\alpha}} \nabla p + \frac{1}{\varepsilon^2} \rho u + \frac{1}{\varepsilon} \rho \nabla \phi &= 0, \\ -\Delta \phi &= \rho - M, \end{aligned}$$
where $p = \rho^\gamma$ ($\gamma > 1$), $\varepsilon \ll 1$ is a small parameter representing the inverse of the friction coefficient, and $\alpha > 0$ is a fixed parameter. The force term $\rho \nabla \phi$ is repulsive (unlike the attractive force in self-gravitating systems).

The Limit:
As $\varepsilon \to 0$ (the large friction limit), the system formally converges to the hyperbolic-elliptic Keller-Segel (KS) system of consumption type:
$$\partial_\tau \sigma - \text{div}(\sigma \nabla \psi) = 0, \quad -\Delta \psi = \sigma - M,$$
where $\sigma$ represents the density and $\psi$ the potential. This limit system is known to exhibit finite-time blow-up (singularity formation) for general initial data, particularly in higher dimensions, due to the aggregation mechanism.

The Challenge:
The primary mathematical challenge is to establish the global existence of smooth solutions for the EP system in the large friction regime and prove that these solutions converge to the global smooth solutions of the limit KS system. This is non-trivial because:

1. The EP system is hyperbolic and generally prone to shock formation.
2. The limit KS system (without vacuum) is also prone to blow-up for general data.
3. The limit is singular: the time derivatives are not uniformly bounded with respect to $\varepsilon$ due to an initial layer where the friction rapidly damps the velocity.

2. Methodology

The author employs a direct singular limit method rather than the relative entropy method used in previous works. The key methodological steps include:

• Reformulation via Keller-Segel Map:
  The velocity $u$ is decomposed using an ansatz: $u = \varepsilon v_\rho + \varepsilon^\alpha w$, where $v_\rho = -\nabla(-\Delta)^{-1}(\rho - M)$ is the velocity field associated with the KS limit. This transforms the original EP system into a new system for $(\rho, w)$ where the singular structure is separated, and the KS structure is explicit.
  $$\begin{aligned} \partial_\tau \rho + \frac{1}{\varepsilon^{1-\alpha}} \text{div}(\rho w) + \text{div}(\rho v_\rho) &= 0, \\ \varepsilon^\alpha \rho \partial_\tau w + \dots &= \text{Source terms involving } \varepsilon. \end{aligned}$$

• Uniform-in-$\varepsilon$ Estimates:
  The core of the proof involves deriving energy estimates that are uniform with respect to $\varepsilon$.

  • Energy Functionals: The author constructs a total energy functional $E_{total}$ involving $H^3$ norms of $\rho - M$ and weighted norms of $w$.
  • Dissipation: A dispersion functional $D_{total}$ is defined to capture the damping effects.
  • Handling Nonlinearity: A crucial step is controlling the nonlinear terms, specifically the term $\|\nabla \rho\|_{L^4}$. In 3D, the smallness of the initial gradient $\|\nabla \rho_0\|_{L^4}$ is required to close the global estimates. In 1D, this condition is not needed due to the specific structure of the 1D potential.

• Spectral Analysis:
  The paper includes a spectral analysis of the linearized operator around the equilibrium $(M, 0)$. This confirms that the equilibrium is linearly stable and justifies the ansatz for the velocity scaling, showing that the velocity perturbation is of order $O(\varepsilon) + O(\varepsilon^{\alpha+1})$, which is consistent with the ansatz used (though the initial data allows for slightly larger perturbations).

• Compactness and Convergence:
  Using the uniform bounds and the Aubin-Lions compactness lemma, the author proves that as $\varepsilon \to 0$, the solutions $(\rho, w)$ converge to the solution $(\sigma, v)$ of the limit KS system.

3. Key Contributions and Results

1. Global Smooth Solutions for the EP System:
The paper establishes the existence of unique global-in-time classical solutions for the almost pressureless EP system with large friction.

• Initial Data Requirements:
  • 3D: Requires smallness of the initial velocity in $H^3$ and the initial density gradient in $L^4$ ($\|\nabla \rho_0\|_{L^4}$). Crucially, the initial density perturbation $\|\rho_0 - M\|_{H^3}$ does not need to be small, provided the non-vacuum condition holds ($\rho_0 > 0$).
  • 1D: Requires only smallness of the initial velocity in $H^2$; no smallness condition on $\|\nabla \rho_0\|_{L^4}$ is needed.
• Mechanism: The global existence is driven by the large damping/friction effect, which suppresses the formation of singularities (shocks) that would otherwise occur in the pressureless limit.

2. Convergence to the Limit System:
The paper proves that the global solutions of the EP system converge to the global smooth solutions of the hyperbolic-elliptic KS system of consumption type as $\varepsilon \to 0$.

• Regularity: The convergence holds in strong topologies ($C(0, T; H^2_{loc})$ for density, weak convergence for time derivatives).
• No Singularity in the Limit: Unlike the general KS system which may blow up, the limit solutions obtained here are global and smooth because they are limits of the globally stable EP solutions.

3. Analysis of Vacuum and 1D Asymptotics:
The paper investigates the behavior when vacuum (zero density) is present, specifically in 1D.

• Vacuum Shrinkage: If the initial data contains a vacuum interval, the KS flow causes this interval to shrink exponentially to a point as $t \to \infty$.
• Loss of Regularity: While the solution approaches the equilibrium $M$ almost everywhere, the Sobolev norms (specifically $H^2$) of the solution grow exponentially in time due to the vacuum. This explains why the non-vacuum condition is essential for the global uniform estimates and the validity of the large friction limit in the current framework.

4. Significance

• Resolution of Singularity Formation: The work demonstrates that the large friction limit acts as a regularizing mechanism, preventing the finite-time blow-up typically associated with the pressureless Euler-Poisson and Keller-Segel systems.
• Relaxed Initial Data: Unlike previous results on the EP system which required small perturbations for both density and velocity, this paper allows for large density perturbations (as long as they are bounded away from vacuum) and only requires small velocity perturbations. This is a significant advancement in the theory of global smooth solutions for plasma models.
• Methodological Innovation: The use of the "Keller-Segel map" ansatz to reformulate the system allows for a direct singular limit analysis where the regularity of the EP and KS solutions matches exactly. This avoids the need for the limit solution to be more regular than the primitive system, a common constraint in relative entropy methods.
• Physical Insight: The results provide a rigorous mathematical justification for using the Keller-Segel model to describe plasma dynamics in high-friction regimes, provided the density remains non-vacuous.

In summary, this paper provides a rigorous proof of the global stability of the almost pressureless Euler-Poisson system under large friction, establishes its convergence to a global smooth Keller-Segel system, and clarifies the critical role of the non-vacuum condition in maintaining regularity.