Poroplate Dynamics: Fluid-Structure Coupling

• Poroplate dynamics is the study of thin porous plate deformation coupled with transverse pressure diffusion, derived rigorously from the 3D Biot system.
• The model uses asymptotic analysis to reduce complex interactions into a 2.5D framework that captures key features of elastic and fluid coupling.
• Applications span layered filtration, biological tissue mechanics, and engineered composites, with semigroup theory ensuring dynamic well-posedness.

Poroplate dynamics concerns the mathematical modeling, analysis, and physical interpretation of systems where the elastic deformation of thin porous plates is coupled with the perfusive dynamics of saturating fluids. The “poroplate” paradigm emerges as a rigorous limit of the three-dimensional (3D) Biot poroelasticity system under thin-plate scaling, yielding a so-called “2.5D” model that retains out-of-plane pressure diffusion and its coupling to transverse plate deformation. Poroplate dynamics is central to the understanding of layered filtration, fluid-structure interaction in biological tissues, and the mechanics of engineered sandwich materials, especially where plates act as interfaces between bulk flow domains or poroelastic substrates. The resulting partial differential equation (PDE) systems feature highly non-trivial mathematical structure: degenerate parabolicity, non-autonomous coefficients, strong interface coupling, and non-classical boundary and transmission laws.

1. Mathematical Formulation and Plate-Biot Coupling

The derivation of poroplate dynamics begins with the 3D Biot system on a domain $$\Omega_p^\varepsilon = \omega \times (-\varepsilon h, +\varepsilon h)$$, where $\omega \subset \mathbb{R}^2$ represents the in-plane mid-surface and $\varepsilon \ll 1$ the aspect ratio (thinness). The equations governing the solid displacement $u^\varepsilon$ and pore pressure $p^\varepsilon$ are: $$\begin{cases} \rho_s \partial_{tt} u^\varepsilon - \nabla \cdot (C \varepsilon(u^\varepsilon)) + \alpha \nabla p^\varepsilon = F^\varepsilon\ \partial_t(c_0 p^\varepsilon + \alpha \nabla \cdot u^\varepsilon) - \nabla \cdot (k \nabla p^\varepsilon) = G^\varepsilon \end{cases}$$ After nondimensionalization and formal asymptotic analysis as $\varepsilon \to 0$, the out-of-plane displacement $w(x_1, x_2, t)$ (Euler-Bernoulli plate kinematics) and fluid pressure $p(x_1, x_2, z, t)$ become the principal unknowns. The reduced “2.5D” system for the plate-poroelastic regime is: \begin{equation} \begin{cases} \rho_p \partial_{tt} w + D \Delta^2 w + \alpha \Delta \displaystyle\int_{-h}^h z\, p(\cdot, z, t)\,dz = f(x, t)\ \partial_t[c_p p - \alpha z \Delta w] - \partial_z[k(x, z, t) \partial_z p] = g(x, z, t) \end{cases} \end{equation} where $D$ is flexural rigidity, $\rho_p$ is plate density, $c_p$ is fluid storage, and $\alpha$ is the Biot-Willis constant. Notably, the pressure diffusion acts only in the transverse ($z$) direction, a distinguishing feature of poroplate models (Gurvich et al., 2021).

2. Interface Problems and Multilayer Coupling

Poroplate dynamics is further generalized in multilayered systems, where a poroelastic plate serves as an interface between bulk free-flow regions (Stokes), bulk Biot poroelastic domains, or both. In these setups, the governing equations for each region are supplemented by physically justified interface conditions:

• Conservation of mass: continuity of normal Darcy flux across the interface,
• Stress balance: continuity of normal stresses (including contributions from plate bending and fluid pressure),
• Tangential slip: Beavers–Joseph–Saffman-type slip conditions for fluid velocity on the plate interface.

The Mikelić plate model at the interface introduces additional operators: $$\mathcal{K}(p_p) = \int_{-h/2}^{h/2} p_p(x', s)\, ds,\quad \widetilde{\mathcal{K}}(q)(x', s) = s\,q(x')$$ and the governing PDEs: $$\begin{cases} \rho_p w_{tt} + D \Delta^2 w + \gamma w + \alpha_p \Delta[\mathcal{K}(p_p)] = F_p\ [c_p p_p - \alpha_p \widetilde{\mathcal{K}}(\Delta w)]_t - \partial_s(k_p \partial_s p_p) = 0 \end{cases}$$ with precise interface coupling to bulk Biot and Stokes equations via mass and stress continuity, as established in (Avalos et al., 1 Feb 2026).

3. Functional Analysis, Operator Framework, and Weak Solution Theory

Solutions to poroplate PDEs are formulated in anisotropic Sobolev spaces adapted to the plate (e.g., $W = H^2(\omega) \cap H^1_0(\omega)$) and pressure fields ($V$ for pressure with normal-derivative vanishing boundaries). For the quasi-static case ($\rho_p = 0$), the PDE system is equivalently written as an implicit evolution equation: $\partial_t(c_p I + B)p + A(t)p = g^*$ where $A(t): V \to V'$ is the time-dependent (possibly degenerate) diffusion operator in $z$, and $B$ is a bounded memory operator arising from the bending-pressure coupling: $$B = \beta \mathcal{K}^* \Delta^{-1} \mathcal{K},\quad \beta = \alpha^2/D$$ Weak solutions are defined via duality pairings and integral constraints involving test functions from the relevant Sobolev classes. Existence is established using the Lions–Showalter theorem for implicit evolution equations under minimal regularity and coercivity assumptions on $k(x, z, t)$. Uniqueness requires additional time-regularity of the permeability function (Gurvich et al., 2021).

4. Semigroup Theory and Dynamic Well-posedness

When inertia is retained ($\rho_p > 0$), the system is recast as a first-order evolution in a Hilbert phase space, encompassing the plate displacement, velocity, and pressure: $$X = (H^2(\omega) \cap H^1_0(\omega)) \times L^2(\omega) \times L^2(\Omega)$$ The infinitesimal generator $\mathbf{A}$ encapsulates dynamic plate, pressure diffusion, and coupling: $$\mathbf{A} = \begin{pmatrix} 0 & I & 0 \ -D \Delta^2 & 0 & -\alpha \Delta \int_{-h}^h z(\cdot, z)dz \ -\alpha z \Delta & 0 & -\partial_z(k \partial_z) \end{pmatrix}$$ With suitable domain conditions encoding both physical and mathematical boundary/interface laws, dissipativity and $m$-dissipativity are established. The Lumer–Phillips theorem guarantees $\mathbf{A}$ generates a contraction $C_0$-semigroup, yielding unique mild and strong solutions (Gurvich et al., 2021, Avalos et al., 1 Feb 2026). A parallel theory extends to multilayered systems, where the state space incorporates all bulk and interface variables, and the full system generator is shown to admit a semigroup representation under strong interface transmission constraints.

5. Nonlinear Dynamics and Perturbation Theory

Nonlinear structural effects, such as von Kármán-type large deflection, can be incorporated into the plate equation as locally Lipschitz perturbations. For instance, the addition

$f(w) = [w, v(w)] + F_0$

with Airy stress $v(w)$ subject to biharmonic constraints, leads to a nonlinear operator $\mathcal{F}$ acting on the phase space. By Pazy’s theorem, $A + \mathcal{F}$ (with $A$ the linear generator) continues to generate a local-in-time mild solution. Global-in-time estimates are obtained via energy bounds incorporating plate potential energy $\Pi(w)$, with Grönwall’s lemma ensuring dissipativity (Avalos et al., 1 Feb 2026).

6. Physical Interpretation, Regularity, and Stability

Poroplate dynamics provides mechanistic insight into several physical effects:

• Bending operator $D\Delta^2$: resists curvature, governing plate flexibility.
• Pressure-moment coupling $\alpha \Delta \int z p\,dz$: mediates feedback between fluid perfusion and elastic response; fluid content changes can induce or resist plate motion.
• Degenerate parabolicity: pressure diffuses only through plate thickness ($z$), isolated from in-plane diffusion, modeling filtration-like behaviors.
• Dissipation channels: arise both from bulk diffusion (in the Biot and Stokes regions), and from plate pressure diffusion and slip drag at interfaces. In systems with a plate interface, additional regularizing and stabilizing effects are observed due to the enhanced dissipation and coupling structure, as spectral analysis of the generator demonstrates altered decay rates and potential stabilization across frequency regimes (Avalos et al., 1 Feb 2026).
• Comparison to Biot–Stokes systems: The introduction of a poroplate interface in multilayer filtration increases both the mathematical complexity and the physical richness of the problem, with extra buffering and enhanced damping due to plate bending and perfusion.

7. Outlook and Open Problems

The rigorous treatment of poroplate dynamics has catalyzed advances in interface modeling for multi-physics filtration and fluid-structure systems. The establishment of semigroup solvability for linear and (via perturbation theory) certain nonlinear problems unlocks future stability and regularity analyses. Key open issues include:

• Systematic comparison of spectral gaps and decay in different interface configurations,
• Long-time regularity with less restrictive assumptions on coefficients (e.g., permeability variability),
• Analysis of strong nonlinearities or non-monotone couplings in biological or synthetic composite settings.

References: (Gurvich et al., 2021, Avalos et al., 1 Feb 2026)