Schauder-type Estimates and Well-posedness for Nonlocal Quasilinear Evolution Equations in Fluid Dynamics

1. Problem Statement

The paper addresses the fundamental challenge of establishing rigorous well-posedness (existence, uniqueness, and stability) for nonlocal quasilinear parabolic evolution equations arising in fluid dynamics. These equations typically feature:

• Nonlocality: The highest-order dissipative terms are pseudo-differential operators (e.g., fractional Laplacians $\Lambda^s$) rather than local differential operators.
• Quasilinearity: The coefficients of the principal operator depend on the solution itself (specifically its derivatives), making the equation fully nonlinear in the highest order.
• Critical Regularity: The goal is to prove well-posedness in scaling-critical spaces (specifically Besov spaces $\dot{B}^1_{\infty,\infty}$ and Lipschitz spaces $\dot{W}^{1,\infty}$), which are the lowest regularity spaces where the scaling invariance of the equation is preserved.

The authors focus on three specific physical models:

1. The 2D Muskat Equation with Surface Tension: Describing the interface between two immiscible fluids in a porous medium.
2. The 2D Peskin Problem: Modeling an elastic string immersed in a 2D Stokes fluid.
3. The 3D Peskin Problem: Modeling an elastic membrane immersed in a 3D Stokes fluid.

2. Methodology

The core of the paper is the development of a unified analytical framework based on Schauder-type estimates for linear parabolic systems with variable coefficients.

A. Kernel-Adapted Freezing-Coefficient Method

Traditional Schauder theory for variable-coefficient equations often treats coefficient variations as lower-order perturbations, which usually restricts results to short time intervals or requires smallness of the coefficients. The authors introduce a novel kernel-adapted freezing-coefficient method:

• Freezing: At a reference point $x_0$, the variable coefficients are frozen to create a constant-coefficient operator.
• Explicit Representation: Instead of treating the difference as a perturbation immediately, they derive an explicit integral representation formula for the solution using the fundamental kernel (Green's function) of the frozen operator.
• Evaluation: The solution is evaluated at the physical point $x$. The error terms (residuals) arising from the coefficient variation are explicitly controlled within the leading-order dynamics.
• Advantage: This approach avoids artificial spatial truncations and provides robust control of residual terms, allowing for estimates that are valid globally in time (under smallness conditions) and in critical spaces.

B. Time-Weighted Critical Spaces

The analysis is conducted in time-weighted Hölder/Besov spaces. The norms are defined as:
$$\|u\|_T = \sup_{t \in [0,T]} \left( \|u(t)\|_{\dot{C}^b} + t^{\frac{m+\kappa}{s}} \|u(t)\|_{\dot{C}^{m+\kappa+b}} \right)$$
This weighting allows the solution to exhibit instantaneous smoothing (gaining regularity for $t>0$) even if the initial data is rough (e.g., in $\dot{B}^1_{\infty,\infty}$). This is crucial for handling critical data where standard fixed-point arguments in Lipschitz spaces fail.

C. Strategy for Quasilinear Equations

To handle the quasilinear nature:

1. Linearization: The equation is linearized around a smooth reference profile $\Phi$.
2. A Priori Estimates: The new Schauder estimates are applied to the linearized system to derive bounds on the solution in terms of the initial data and the reference profile.
3. Approximation & Compactness: Since the critical space $\dot{B}^1_{\infty,\infty}$ is not a Banach algebra and does not support standard contraction mapping directly, the authors first prove a "relaxed" version in Lipschitz spaces ($\dot{W}^{1,\infty}$). They then construct a sequence of smooth approximations, apply the uniform a priori estimates, and pass to the limit using compactness arguments to obtain solutions in the critical Besov space.

3. Key Contributions

1. General Schauder Estimates (Theorem 1.1)

The authors establish sharp Schauder-type estimates for linear parabolic systems driven by variable-coefficient nonlocal pseudo-differential operators of order $s > 0$.

• Scope: Applies to operators with coefficients depending on time, space, and the solution's derivatives.
• Result: Provides bounds in critical time-weighted Hölder/Besov spaces. Crucially, if the coefficient norm is sufficiently small, the estimates hold globally in time.

2. Well-posedness for the Muskat Equation with Surface Tension

• Result: Proves local well-posedness for arbitrary initial data in $W^{1,\infty}$ and global well-posedness for small initial data in the critical space $\dot{B}^1_{\infty,\infty}$.
• Significance: This is the first result establishing well-posedness for the Muskat equation with surface tension in the largest critical space. It handles the case where the heavier fluid is below (stable) and the case where the heavier fluid is above (unstable, yielding only local existence).
• Mechanism: Surface tension provides a third-order dissipative term ($\Lambda^3$), which dominates the dynamics. The authors show that the nonlocal quasilinear structure fits their framework, overcoming the difficulty that the principal coefficient depends on the solution's slope.

3. Well-posedness for the 2D and 3D Peskin Problems

• 2D Case: Proves local well-posedness in the critical Besov space $\dot{B}^1_{\infty,\infty}$. They also extend this to global well-posedness and asymptotic stability for solutions close to circular equilibrium states, utilizing the specific spectral properties of the linearized operator around a circle.
• 3D Case: Proves local well-posedness in critical spaces.
• Innovation: The 3D case is particularly challenging due to the lack of a Hilbert transform structure (present in 2D) and complex geometry. The authors use stereographic projection to map the problem to $\mathbb{R}^2$, reformulating the nonlocal operator into a form compatible with their Schauder estimates. They rigorously control the error terms introduced by the coordinate transformation and the curvature of the sphere.

4. Main Results Summary

Model	Space	Result	Conditions
2D Muskat (with $\sigma > 0$)	$\dot{B}^1_{\infty,\infty}$	Global Well-posedness	Small initial data; $\rho_0 = 0$ (no gravity jump) or stable regime.
2D Muskat (with $\sigma > 0$)	$W^{1,\infty}$	Local Well-posedness	Arbitrary initial data; general $\rho_0$.
2D Peskin	$\dot{B}^1_{\infty,\infty}$	Local Well-posedness	Arc-chord condition satisfied.
2D Peskin	$\dot{B}^1_{\infty,\infty}$	Global Well-posedness	Small perturbation of a circular equilibrium.
3D Peskin	$W^{1,\infty}$ / Critical	Local Well-posedness	Arc-chord condition satisfied; uses stereographic projection.

5. Significance and Impact

• Unified Framework: The paper provides a single, robust analytical framework capable of handling diverse nonlocal fluid dynamics problems (Muskat, Peskin) that were previously treated with disparate, model-specific techniques.
• Critical Regularity: By working in $\dot{B}^1_{\infty,\infty}$, the results are optimal with respect to scaling. This is a significant improvement over previous results that required higher regularity (e.g., $H^s$ with $s > d/2 + 1$) or subcritical spaces.
• Handling Nonlocality: The "kernel-adapted freezing-coefficient" method offers a new tool for analyzing variable-coefficient nonlocal operators, potentially applicable to other areas like kinetic theory (Boltzmann/Landau equations) and fractional diffusion.
• Geometric Complexity: The successful extension to the 3D Peskin problem demonstrates that the framework can handle high-dimensional geometric complexities and non-Euclidean manifolds without relying on specific symmetries (like the Hilbert transform in 2D).

In conclusion, this work represents a major advancement in the mathematical theory of free boundary problems and fluid-structure interactions, establishing the well-posedness of these complex systems at the lowest possible regularity threshold.