Table of Contents
Fetching ...

Well-Posedness of the Free Boundary Incompressible Porous Media Equation

Mickaël Latocca, Huy Q. Nguyen

TL;DR

This work tackles local well-posedness for the free boundary incompressible porous media equation, a quasilinear IPM system with a moving interface between fluid and dry regions under gravity. The authors identify a stability functional $\mathcal{T}(f)=\gamma(f)-G[f]\Gamma(f)$ and prove local existence-uniqueness for small density perturbations around stratified states whenever $\inf_x \mathcal{T}(f_0)(x)>0$, using a robust combination of Dirichlet-Neumann operator analysis, para-differential calculus, and a domain-flattening technique. Central contributions include tame fractional Sobolev bounds for operators linking boundary data to Poisson solutions in Sobolev domains and contraction estimates that enable a fixed-point construction in Sobolev spaces, covering both infinite and finite depth geometries. The results advance the mathematical understanding of Muskat-like free boundary problems with variable density and provide tools potentially applicable to related stability analyses and numerical schemes for porous-media flows.

Abstract

We consider the free boundary incompressible porous media equation which describes the dynamics of a density transported by a Darcy flow in the field of gravity, with a free boundary between the fluid region and the dry region above it. For any stratified density state, we identify a stability condition for the initial free boundary. Under this condition, we prove that small localized perturbations of the stratified density lead to unique local-in-time solutions in Sobolev spaces. Our proof involves analytic ingredients that are of independent interest, including tame fractional Sobolev estimates for operators that map the Dirichlet boundary function and the forcing function of Poisson's equation to its solution in domains of Sobolev regularity.

Well-Posedness of the Free Boundary Incompressible Porous Media Equation

TL;DR

This work tackles local well-posedness for the free boundary incompressible porous media equation, a quasilinear IPM system with a moving interface between fluid and dry regions under gravity. The authors identify a stability functional and prove local existence-uniqueness for small density perturbations around stratified states whenever , using a robust combination of Dirichlet-Neumann operator analysis, para-differential calculus, and a domain-flattening technique. Central contributions include tame fractional Sobolev bounds for operators linking boundary data to Poisson solutions in Sobolev domains and contraction estimates that enable a fixed-point construction in Sobolev spaces, covering both infinite and finite depth geometries. The results advance the mathematical understanding of Muskat-like free boundary problems with variable density and provide tools potentially applicable to related stability analyses and numerical schemes for porous-media flows.

Abstract

We consider the free boundary incompressible porous media equation which describes the dynamics of a density transported by a Darcy flow in the field of gravity, with a free boundary between the fluid region and the dry region above it. For any stratified density state, we identify a stability condition for the initial free boundary. Under this condition, we prove that small localized perturbations of the stratified density lead to unique local-in-time solutions in Sobolev spaces. Our proof involves analytic ingredients that are of independent interest, including tame fractional Sobolev estimates for operators that map the Dirichlet boundary function and the forcing function of Poisson's equation to its solution in domains of Sobolev regularity.

Paper Structure

This paper contains 43 sections, 62 theorems, 830 equations, 2 figures.

Key Result

Theorem 1.4

Let $d\geqslant 1$ and $\gamma \in \mathfrak{G}$, defined in defi:gamma. Let $s> \frac{3}{2}+\frac{d}{2}$ and $(f_0,g_0) \in H^s(\mathbb{R}^d)\times H^s(\Omega_{f_0})$, where $\Omega_{f_0}$ is defined by either eq.domain (infinite depth) or eq.domain.finite (finite depth). Assume that In the finite-depth case we also assume that $b(x)=-H+b_0(x)$ with $H>0$ and $b_0 \in H^{s+\frac{1}{2}}(\mathbb{R

Figures (2)

  • Figure 1: The domain $\Omega_{f(t)}$ in the finite-depth case.
  • Figure 2: The domains $\mathcal{U}_i$

Theorems & Definitions (118)

  • Definition 1.2
  • Remark 1.3
  • Theorem 1.4
  • Remark 1.5
  • Remark 1.6
  • Theorem 2.1: ABZNP
  • Theorem 2.2
  • Lemma 2.3
  • proof
  • Remark 2.4
  • ...and 108 more