Two-Phase Free Interface Problems

• Two-phase free interface problems are defined by regions of differing phases separated by an unknown free boundary with discontinuities in physical and geometric quantities.
• Advanced variational methods and PDE techniques reveal sharp regularity results, including C¹,α bounds for 0 < γ < 1 and Log-Lipschitz continuity for γ = 0.
• The analysis has practical implications in cavitation, fluid vortex dynamics, and reaction-diffusion, employing innovative comparison methods and stability estimates.

A two-phase free interface problem involves regions of different physical or mathematical phases separated by an a priori unknown interface—referred to as the "free boundary"—across which physical or geometric quantities may have discontinuities or singular behaviors. The interface's evolution and regularity are governed simultaneously by the bulk equations in each phase and specific coupling (jump or transmission) conditions prescribed on the interface. These problems arise in nonlinear elliptic and parabolic PDEs, with singular or degenerate structures arising from non-differentiable potentials, discontinuous coefficients, or degenerate absorption terms. Prominent applications include cavitation, two-phase fluid mechanics, phase transitions, and combustion. The analysis of two-phase free interface problems requires advanced PDE theory to address singularities at the interface and to obtain sharp regularity estimates for both solutions and the free boundary itself.

1. Variational Formulation and Singular Energies

A general two-phase free interface problem is formulated as the minimization of variational energies with singular, non-differentiable terms. A canonical example is the minimization of: $$J_\gamma(v) = \int_\Omega \left( |\nabla v(x)|^p + F_\gamma(v(x)) + f(x) v(x) \right) dx$$ where $v$ is sought in an appropriate Sobolev space with boundary data, $p > 1$, and $F_\gamma$ is a non-smooth potential that encodes the two-phase nature: $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$ with $v^+ = \max\{v, 0\}$, $v^- = \min\{v, 0\}$, and in the limiting case $\gamma = 0$, $$F_0(v) = \lambda_+ \chi_{\{v>0\}} + \lambda_- \chi_{\{v<0\}}$$ where $\chi$ denotes the characteristic function. The parameters $v$0 and $v$1 modulate the strength and discontinuity of absorption/jump effects between phases.

The associated Euler–Lagrange equations become singular along the free interface $v$2. For $v$3 minimizer: $v$4 where $v$5 is the $v$6-Laplacian. For $v$7, the right-hand side becomes distributional, and the interface condition emerges as a flux jump: $v$8 with $v$9 the unit normal to the interface.

This framework captures a wide range of problems, including degenerate absorption, obstacle-type systems, and models with non-smooth double-well potentials.

2. Degenerate Operators and Interface Singularities

Degenerate operators, most notably the $p > 1$0-Laplacian with $p > 1$1, induce geometric and analytic challenges distinct from linear cases. Regularity theory for $p > 1$2-harmonic functions provides only $p > 1$3 continuity, with the exponent $p > 1$4 dependent on $p > 1$5 and the spatial dimension. The presence of a degenerate absorption term $p > 1$6, especially as $p > 1$7, causes the right-hand side of the Euler–Lagrange equation to blow up as the interface is approached, and standard tools, such as monotonicity formulas, are unavailable. The interface is not prescribed and must be found as part of the solution—often an implicit, highly non-smooth set along which the differential operator degenerates or becomes singular.

The case $p > 1$8 corresponds to the sharpest singular regime, typical of cavitation, jet, or vortex-core problems, where the model enforces the absorption/discontinuity directly via the indicator functions. Such constructions appear in heterogeneous Prandtl–Batchelor-type models in fluid mechanics (vortex regions with prescribed vorticity surrounded by irrotational flow), in chemical diffusion, and in other free boundary phenomena.

3. Sharp Regularity Theory for Minimizers

The regularity of minimizers in these problems is highly sensitive to the degeneracy and the absorption singularity. The principal findings are:

• For $p > 1$9, every local minimizer is $F_\gamma$0 in the interior, with the sharp Hölder exponent $F_\gamma$1 dependent on dimension $F_\gamma$2, $F_\gamma$3, and $F_\gamma$4:

$F_\gamma$5

where $F_\gamma$6 is the optimal exponent for $F_\gamma$7-harmonic functions and $F_\gamma$8 describes the integrability of $F_\gamma$9.

• For $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$0, minimizers are universally Log-Lipschitz continuous:

$$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$1

at the borderline integrability $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$2.

These regularity results are sharp—that is, optimal with respect to the structure of the coefficients, operator degeneracy, and singularity strength. The approach leverages comparison with $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$3-harmonic functions, Campanato-type arguments, and the fine analysis of singular absorption, extending beyond purely variational minimality to encompass viscosity/weak formulations for the interface conditions.

The distinction between the $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$4 and Log-Lipschitz regimes reflects a phase transition in the interface regularity as the strength of the singular potential changes ($$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$5).

4. Physical Models, Applications, and Interface Geometry

Such two-phase free interface problems, as formulated above, are central to multiple physical phenomena, including:

• Heterogeneous cavitation and fluid vortex dynamics: The limiting case ($$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$6) yields Prandtl–Batchelor-type models where constant vorticity regions (jets or cavities) are embedded within irrotational flow, and the interface is determined as a sharp vortex boundary.
• Reaction–diffusion in porous media: The absorption term models concentrated sources or sinks where chemical density undergoes sudden disappearance or ignition at the interface.
• Obstacle-type systems with degenerate nonlinearities: The potential encodes non-smooth phase constraints, typical of contact or phase separation problems.

The improved regularity of both solutions and interfaces is pivotal for analyzing interface stability, shape optimization, and geometric evolution. For instance, in cavitation and vortex core detection, the sharp $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$7 or Log-Lipschitz continuity ensures geometric control of the interface, excluding fractal or pathological behaviors in relevant parameter regimes.

5. Novel Analytical Techniques and Stability

The derivation of these sharp regularity results introduces several innovations:

• Quantitative $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$8 regularity via effective comparison functions: Using suitable $$F_\gamma(v) = \lambda_+ (v^+)^{1+\gamma} + \lambda_- (v^-)^{1+\gamma}, \qquad (0 < \gamma \leq 1)$$9-harmonic comparison and Campanato's method, the regularity up to the free boundary is achieved, showing that singular absorption does not destroy interior gradient bounds for positive $v^+ = \max\{v, 0\}$0.
• Asymptotically optimal Log-Lipschitz bounds at the degenerate endpoint ($v^+ = \max\{v, 0\}$1), not attainable via classical energy estimates.
• Stability and convergence of minimizers as $v^+ = \max\{v, 0\}$2: Minimizers of $v^+ = \max\{v, 0\}$3 converge to minimizers of the limiting cavitation energy $v^+ = \max\{v, 0\}$4, demonstrating the robustness of the approach across the parameter space.
• Lower gradient (non-degeneracy) bounds close to the free boundary, ensuring sharp transitions between phases.
• Novel comparison arguments and BMO-based tools that compensate for the lack of monotonicity formulas in the degenerate regime.

These methods are essential for extending regularity theory from variational to fully nonlinear, degenerate, or singular free boundary PDEs.

6. Broader Implications, Limitations, and Future Directions

The analytical framework and regularity theory established generalize and strengthen classical results, providing explicit dependence of regularity exponents on dimension, operator degeneracy, and absorption singularity. Notably, the results are novel even for linear ($v^+ = \max\{v, 0\}$5) cases—a context where the interplay between non-differentiable potentials and interface geometry had not been fully resolved.

Limitations stem from the intrinsic singularity: for example, in the strongest singular regime ($v^+ = \max\{v, 0\}$6), one can at best expect Log-Lipschitz continuity; full Lipschitz or higher regularity is out of reach without further structural assumptions.

This body of work underpins subsequent advances in the study of more complex or higher-dimensional two-phase degenerate free boundary models, including those with variable exponent or nonlocal (fractional) operators, furnishing a unifying analytic foundation for the geometric and functional analysis of such interfaces.

---

Key Regularity Table

Parameter regime	Regularity of minimizer	Dependence of exponent
$v^+ = \max\{v, 0\}$7	$v^+ = \max\{v, 0\}$8	$v^+ = \max\{v, 0\}$9 explicit in $v^- = \min\{v, 0\}$0, $v^- = \min\{v, 0\}$1, $v^- = \min\{v, 0\}$2
$v^- = \min\{v, 0\}$3	Log-Lipschitz	Universal

This classification highlights the singular limit's effect on regularity and encapsulates the main analytical insight provided by (Leitão et al., 2012).