A fluid-solid interaction problem in porous media

In this work, we derive asymptotic interface models for an elastic Muskat free boundary problem describing Darcy flow beneath an elastic membrane. In a weakly nonlinear regime of small interface steepness, we obtain nonlocal evolution equations that capture the free-boundary dynamics up to quadratic order. In the long-wave thin-film regime, we rewrite the kinematic condition in flux form, flatten the moving domain, and derive a lubrication-type equation. Moreover, we establish well-posedness for these models in suitable Wiener spaces.

💡 Research Summary

**
The paper investigates a coupled fluid–solid interaction problem in a porous medium, where a Darcy flow beneath an elastic membrane gives rise to a one‑phase Muskat free‑boundary system. The authors focus on two asymptotic regimes that simplify the dynamics: (i) a weakly nonlinear regime characterized by small interface slopes, and (ii) a long‑wave thin‑film regime where the depth of the porous layer is large compared to the horizontal wavelength.

In the first regime, the free surface is represented as a graph (h(x,t)) and the bulk Darcy velocity is expressed through a harmonic potential (\Phi). By expanding the Dirichlet‑to‑Neumann (DtN) operator around the flat equilibrium up to quadratic order in the slope parameter (\sigma), the authors obtain a nonlocal evolution equation involving the operator (\Lambda\tanh(\Lambda)) (with (\Lambda = \sqrt{-\Delta})). Two versions of the weakly nonlinear model are derived. The first retains the full quadratic nonlocal commutator, leading to an equation of the form (1.6) where the time derivative is acted upon by a nonlocal operator and the quadratic term contains another nonlocal commutator. The second model simplifies the quadratic term by replacing the time derivative of the height inside the nonlinearity with its leading‑order expression, resulting in a slightly more tractable equation (1.7)–(1.8). Both models incorporate the elastic restoring force (a Willmore‑type operator) and a tangential dissipative term, which together generate higher‑order spatial derivatives and nonlocal interactions.

The analytical core of the paper is the well‑posedness theory for these weakly nonlinear equations in Wiener spaces (A_s), defined via exponential weights on Fourier coefficients. By assuming zero‑mean initial data (h_0) with sufficiently small (A_1) norm, the authors prove existence and uniqueness of mild solutions. For the case (\lambda=0) (no fourth‑order regularization), local well‑posedness holds in (A_1); if additionally (h_0\in A_3), the solution extends globally and decays exponentially in the (A_0) norm (Theorem 5.1). When (\lambda>0), global well‑posedness is established in (A_3) with exponential decay in (A_0) (Theorem 5.2). The alternative model (1.7)–(1.8) enjoys the same global results (Theorem 5.4). The proofs rely on precise symbol estimates for the nonlocal operators, bilinear bounds in Wiener spaces, and a contraction‑mapping argument applied to the integral formulation of the equations.

In the second regime, the authors consider the long‑wave thin‑film limit (\delta\to\infty). They rewrite the kinematic condition in conservative (flux) form, flatten the moving domain onto a fixed strip, and perform a systematic asymptotic expansion. The resulting leading‑order lubrication equation (1.9) features a variable mobility (\mu(h)=\Theta B_{xx}B_t-\chi h-(\lambda/4)B_{xxxx}) that couples the time derivative of the height with the elastic term. This coupling produces a genuinely nonlinear elliptic operator acting on (B_t), a feature absent in classical thin‑film models.

The paper then establishes a global well‑posedness result for the thin‑film equation in the Wiener space (A_4). Assuming zero‑mean initial data with small (A_1) norm, a unique global mild solution (h\in C(