Existence, asymptotic behaviour and convergence of a generalised 3D Muskat problem in stable regime

1. Problem Statement

The paper investigates a generalised three-dimensional $\alpha$-Muskat problem, which models the evolution of the interface between two incompressible fluids with different densities in a porous medium.

• Physical Context: The system describes the interface $x_3 = f(x_1, x_2, t)$ separating two domains $\Omega_1$ (denser fluid, $\rho_2$) and $\Omega_2$ (lighter fluid, $\rho_1$). The problem is considered in the stable regime ($\rho_2 > \rho_1$), meaning the denser fluid is below the lighter one.
• Mathematical Formulation: The model modifies the classical Darcy's law by introducing a fractional derivative operator. The velocity field $u$ is related to the density $\rho$ via a pseudo-differential operator of order $\alpha \in [0, 1)$:
  $$u + \nabla P = -(0, 0, (-\Delta)^{\alpha/2} \rho)$$
  This leads to a contour equation for the interface $f$:
  $$\partial_t f(x, t) = \frac{\rho_2 - \rho_1}{4\pi} \text{PV} \int_{\mathbb{R}^2} \frac{(\nabla f(x) - \nabla f(x-y)) \cdot y}{[|y|^2 + (f(x) - f(x-y))^2]^{\frac{3+\alpha}{2}}} dy$$
  For simplicity, the constant factor is normalized to 1.
• Goal: To establish well-posedness, global existence, decay rates, and convergence properties of the solution $f$ as the parameter $\alpha$ varies, specifically focusing on the 3D case which had not been fully addressed for this generalised model.

2. Methodology

The authors employ a combination of energy estimates, Fourier analysis, maximum principles, and asymptotic analysis.

• Energy Estimates: For local well-posedness, the authors use standard $L^2$ and $H^k$ energy methods. They multiply the contour equation by $f$ (or its derivatives) and integrate, carefully handling the singular integral terms via integration by parts and splitting the domain into $|y|<1$ and $|y|>1$.
• Fourier Analysis & Norms: To prove global existence, they introduce a specific weighted $L^1$-type norm in Fourier space: $\|f\|_s = \int_{\mathbb{R}^2} |\xi|^s |\hat{f}(\xi)| d\xi$. This allows them to control the nonlinearity using Taylor expansions of the kernel function.
• Maximum Principles: The authors derive maximum principles for the $L^\infty$ norms of the solution $f$ and its gradient $\nabla f$. This involves analyzing the evolution of the maximum point of the function using the Rademacher theorem and properties of hypergeometric functions.
• Convergence Analysis: To study the limit $\alpha \to 0^+$, they define the difference $g_\alpha = f_\alpha - f$ and derive an evolution equation for the $L^2$ norm of the difference, utilizing uniform bounds on the solutions to apply Grönwall's inequality.

3. Key Contributions and Results

A. Local-in-Time Well-Posedness

• Result: For any $\alpha \in [0, 1)$ and initial data $f_0 \in H^k(\mathbb{R}^2)$ with $k \ge 4$, there exists a time $T > 0$ such that the system possesses a unique solution $f \in C^1([0, T]; H^k(\mathbb{R}^2))$.
• Significance: This establishes the mathematical foundation for the generalised 3D model, extending previous 2D results to three dimensions with fractional dissipation.

B. Global-in-Time Existence (Strong Solutions)

• Result: For $\alpha \in [0, 1/2)$, if the initial data satisfies a smallness condition $\|f_0\|_1 < k_0(\alpha)$ (where $\|f\|_1$ is the Fourier-weighted norm defined in Eq. 3.1), then a unique global solution exists for all $t > 0$.
• Constraint: The global existence is currently restricted to $\alpha < 1/2$. The authors note that the condition relies on ensuring the set $(\alpha, 1-\alpha)$ is non-empty, and the case $\alpha \ge 1/2$ remains an open problem.
• Mechanism: The proof relies on showing that the norm $\|f\|_1$ is non-increasing (a maximum principle in Fourier space) under the smallness condition, preventing blow-up.

C. Maximum Principles and Asymptotic Decay

• $L^\infty$ Decay of $f$: The authors prove that $\|f(t)\|_{L^\infty} \le \|f_0\|_{L^\infty}$. Furthermore, if the initial data does not change sign (e.g., $f_0 \ge 0$), the solution decays algebraically in time:
  $$\|f(t)\|_{L^\infty} \le \frac{\|f_0\|_{L^\infty}}{[1 + C t]^{\frac{2}{1+\alpha}}}$$
• $L^\infty$ Decay of $\nabla f$: Under the assumption that $\|\nabla f_0\|_{L^\infty}$ is sufficiently small (specifically $< \sqrt{\frac{1+\alpha-\epsilon}{5+\alpha+\epsilon}}$) and $\|\nabla f\|_{L^1}$ is bounded, the gradient also decays:
  $$\|\nabla f(t)\|_{L^\infty} \le \frac{\|\nabla f_0\|_{L^\infty}}{[1 + C t]^{\frac{2}{1+\alpha}}}$$
• Significance: These results characterize the long-time behavior of the interface, showing that the interface flattens out over time, with the decay rate dependent on the fractional parameter $\alpha$.

D. Convergence as $\alpha \to 0^+$

• Result: The solutions $f_\alpha$ of the generalised problem converge to the solution $f$ of the classical Muskat problem ($\alpha = 0$) in the $L^2$ and $H^k$ norms as $\alpha \to 0^+$.
• Extension: Under the condition that $f_0$ does not change sign, the convergence also holds in the $L^1$ norm.
• Significance: This rigorously validates the generalised model as a perturbation of the classical Muskat problem, ensuring consistency in the limit of vanishing fractional dissipation.

4. Significance of the Work

1. First 3D Generalised Analysis: This is the first work to address the generalised Muskat problem with fractional dissipation ($\alpha > 0$) in three dimensions. Previous literature largely focused on 2D or the classical 3D case ($\alpha=0$).
2. Handling Non-Local Operators: The paper successfully adapts energy methods to handle the non-local nature of the fractional Laplacian in the velocity law, which introduces significant technical challenges compared to the local Darcy law.
3. Stability and Decay: By establishing global existence for small data and proving explicit decay rates, the paper provides a comprehensive picture of the stability of the interface in the stable regime.
4. Mathematical Rigor: The derivation of the maximum principles for both the function and its gradient, along with the convergence proofs, adds depth to the understanding of active scalar equations with drift velocities defined by pseudo-differential operators.

In summary, this paper bridges the gap between classical porous media flow and fractional diffusion models in 3D, providing rigorous existence, uniqueness, and long-time behavior results for the generalised Muskat interface.