Biot–Kirchhoff Equations in Poroelastic Plates

• Biot–Kirchhoff equations are a class of reduced poroelastic plate models that couple a Darcy-type pressure equation with a Kirchhoff–Love bending law for thin porous structures.
• They derive from 3D Biot poroelasticity via dimension reduction and homogenization, incorporating mixed boundary conditions to capture fluid–structure interactions.
• Recent studies employ rigorous mathematical analysis and advanced discretization methods to ensure strong convergence and reliable error estimates in numerical approximations.

Searching arXiv for recent and foundational papers on Biot–Kirchhoff equations and related poroelastic plate models. First search: Biot–Kirchhoff / poroelastic plate / Kirchhoff-Love thin-layer literature. Biot–Kirchhoff equations denote, in the most precise sense used in current arXiv literature, a class of reduced poroelastic plate models in which a Biot-type pressure or Darcy law is coupled to a Kirchhoff or Kirchhoff–Love plate equation for a thin porous structure. The term is not uniform across the literature: some works explicitly analyze “Biot–Kirchhoff poroelasticity,” some derive a “Biot-plate-system” whose structural part is of Kirchhoff–Love type, some study a “clamped Biot–Kirchhoff–Love poroelastic plate,” and some recent fluid–structure interaction papers formulate bulk–surface or fully averaged poroelastic Kirchhoff plate models instead (Khot et al., 2023, Gahn, 2024, Marciniak-Czochra et al., 2012, Dassi et al., 4 Aug 2025, Brandt et al., 17 May 2026).

1. Terminological scope and model classes

In the surveyed papers, the phrase refers most directly to thin poroelastic plates rather than to standard bulk Biot poroelasticity. The common feature is a fourth-order plate equation for transverse deflection, or a Kirchhoff–Love plate reduction, coupled with a pressure equation representing pore-fluid storage, diffusion, Darcy flow, or a pressure moment across thickness. This suggests that “Biot–Kirchhoff equations” is best treated as an umbrella label for a family of thin-structure poroelastic reductions rather than as a single canonical PDE system (Khot et al., 2023, Marciniak-Czochra et al., 2012, Gahn, 2024, Dassi et al., 4 Aug 2025, Brandt et al., 17 May 2026).

Model family	Primary unknowns	Characteristic coupling
Biot–Kirchhoff poroelasticity	$u$, $p$	fourth-order plate equation coupled with second-order pressure equation
Clamped Biot–Kirchhoff–Love plate	$w^0$, $w_3^0$, $\pi^0$	2D Navier stretching, bending, and pressure diffusion in thickness
Biot-plate-system	$p_0$, $\widetilde{u}_1$, $u_0^n$	generalized Darcy law coupled to a Kirchhoff-Love-type plate equation
Fully averaged poroelastic Kirchhoff plate	$w$, $[q]$, $p$0	all plate unknowns posed on the same mid-surface
Stokes/Biot–Kirchhoff bulk–surface model	$p$1, $p$2, $p$3, $p$4	3D Stokes flow coupled to a surface poroelastic plate

A central distinction in this literature is between plate reductions that retain only mid-surface variables and reductions that keep an explicit thickness variable for pressure. The former include the surface Biot–Kirchhoff and fully averaged models; the latter include the clamped Biot–Kirchhoff–Love derivation, where the pressure equation remains parabolic only in the vertical direction (Marciniak-Czochra et al., 2012, Brandt et al., 17 May 2026, Dassi et al., 4 Aug 2025).

2. Three-dimensional origins and thin-plate reduction

A foundational route to Biot–Kirchhoff equations starts from three-dimensional quasi-static Biot poroelasticity in a thin domain and then passes to the limit as thickness tends to zero. In the clamped plate derivation, the starting body is a thin poroelastic plate, the lateral boundary is clamped, and the authors choose Terzaghi’s time corresponding to the plate thickness. The limit displacement acquires the Kirchhoff–Love structure, the in-plane stretching is described by 2D Navier’s linear elasticity equations, the transverse deflection satisfies a bending equation, and the pore pressure remains as an additional unknown coupled to both stretching and bending (Marciniak-Czochra et al., 2012).

The same thin-structure logic appears in a different microscopic setting in the derivation of a thin periodically heterogeneous poroelastic layer. There the physical domain is

$p$5

both the layer thickness and the microstructure period are of order $p$6, the fluid is governed by quasistatic Stokes equations, and the solid by linear elasticity. The coupling is through continuity of velocities and continuity of normal stress at the fluid–solid interface. Passing to the limit combines homogenization of the fluid-solid microstructure with dimension reduction from a three-dimensional layer to a plate, and the result is a coupled Biot-plate-system whose structural part is of Kirchhoff-Love type (Gahn, 2024).

In both derivations, the thin-plate limit changes the form of poroelastic coupling. The three-dimensional volumetric interaction of standard Biot theory is replaced by reduced couplings involving membrane strain, bending curvature, averaged pressure, or pressure moments across thickness. In the clamped Biot–Kirchhoff–Love plate, the pressure equation is parabolic only in the vertical direction; in the thin periodic layer, the pressure equation becomes a generalized Darcy law on the mid-surface $p$7 (Marciniak-Czochra et al., 2012, Gahn, 2024).

3. Canonical governing equations

A representative Biot–Kirchhoff plate model is the scaled system

$p$8

$p$9

Here $w^0$0 is the averaged-through-thickness deflection and $w^0$1 is the first moment of fluid pressure relative to the solid. The coupling is symmetric in energy: $w^0$2 enters the plate equation through $w^0$3, while $w^0$4 enters the pressure equation through $w^0$5. This is the form explicitly analyzed under the title “Biot-Kirchhoff poroelasticity” (Khot et al., 2023).

The same paper uses mixed boundary conditions. On the clamped with zero-flux part $w^0$6,

$w^0$7

while on the simply supported / pressure-Dirichlet part $w^0$8,

$w^0$9

The corresponding weak formulation is posed in

$w_3^0$0

and the mixed bilinear form is coercive in a parameter-weighted norm because the off-diagonal coupling terms cancel exactly (Khot et al., 2023).

The clamped Biot–Kirchhoff–Love derivation produces a more structured reduced system. The in-plane part is a 2D Navier-type system for $w_3^0$1, the pressure fluctuation $w_3^0$2 satisfies

$w_3^0$3

and the transverse displacement satisfies a bending equation containing

$w_3^0$4

The new term is the first moment of pressure through the thickness, and it acts as a pressure-induced bending moment (Marciniak-Czochra et al., 2012).

A third canonical form arises in the homogenized thin-layer limit. The effective pressure equation is

$w_3^0$5

while the plate equations are an in-plane membrane equilibrium and a fourth-order bending equation for $w_3^0$6. Here the proper identification is a generalized Darcy law coupled to a Kirchhoff-Love plate equation including the Darcy pressure (Gahn, 2024).

4. Bulk–surface and fully averaged formulations

Recent work extends Biot–Kirchhoff equations from isolated plates to coupled fluid–structure systems. In the Stokes/Biot–Kirchhoff bulk–surface model, the plate lies on a flat surface $w_3^0$7, the bulk fluid is governed by Stokes flow in $w_3^0$8, and the surface unknowns are the plate normal displacement $w_3^0$9 and the first moment of the pore pressure head $\pi^0$0. After time semi-discretisation, the surface equations are

$\pi^0$1

$\pi^0$2

with

$\pi^0$3

The coupling conditions are

$\pi^0$4

together with a tangential Beavers–Joseph–Saffman–Jones condition (Dassi et al., 4 Aug 2025).

A closely related but distinct reduction is the fully averaged poroelastic Kirchhoff plate interacting with an incompressible, viscous fluid. The classical Biot plate model recalled there has a mixed-dimensional character: $\pi^0$5 lives on the mid-surface $\pi^0$6, while $\pi^0$7 and $\pi^0$8 live in the three-dimensional plate $\pi^0$9. The fully averaged model replaces the three-dimensional pressure field by the pressure jump

$p_0$0

and the transverse average

$p_0$1

so that both elastodynamics and pressure variables are defined on the same codimension-one manifold $p_0$2. The plate subsystem contains a Kirchhoff bending term $p_0$3, a pressure-jump contribution $p_0$4, and surface equations for $p_0$5 and $p_0$6 (Brandt et al., 17 May 2026).

These models show two distinct reductions of Biot plate theory. The Stokes/Biot–Kirchhoff bulk–surface model retains a pressure-head moment $p_0$7 and tangential Darcy flow on $p_0$8. The fully averaged model replaces the three-dimensional pore pressure by $p_0$9 and $\widetilde{u}_1$0, and the authors stress that this preserves the energy structure of the original Biot plate model while making coupling to a bulk fluid more natural (Dassi et al., 4 Aug 2025, Brandt et al., 17 May 2026).

5. Mathematical analysis and discretisation

Rigorous derivation is a defining feature of this literature. The clamped Biot–Kirchhoff–Love plate paper proves strong convergence of the three-dimensional solid displacement, fluid pressure and total poroelastic stress to the solution of the reduced plate equations. The proof proceeds through uniform estimates, compactness, identification of weak limits, and correction terms that upgrade weak convergence to strong convergence (Marciniak-Czochra et al., 2012).

The thin periodically heterogeneous layer paper proves microscopic well-posedness, compactness and convergence of the microscopic sequence to limit fields, and identifies the unique weak solution of the effective Biot-plate system. The analysis uses two-scale convergence in thin periodic structures, compactness in perforated thin domains, extension operators for thin perforated domains, Korn and Poincaré inequalities with explicit $\widetilde{u}_1$1-dependence, a Bogovskii operator in thin perforated domains, and rescaled test functions adapted to plate kinematics (Gahn, 2024).

For the continuous Biot–Kirchhoff plate problem on polygonal meshes, conforming and nonconforming virtual element methods of arbitrary polynomial degree have been developed. The governing equations consist of one fourth-order equation for the transverse displacement of the middle surface coupled with a second-order equation for the pressure head relative to the solid with mixed boundary conditions. The analysis proves a priori error estimates in the best-approximation form and derives residual-based reliable and efficient a posteriori error estimates in appropriate norms, with robustness with respect to the main model parameters. A notable ingredient is the construction of enrichment operators, referred to as a conforming companion operator in the context of finite element methods, connecting nonconforming virtual element spaces to continuous Sobolev spaces (Khot et al., 2023).

The bulk–surface Stokes/Biot–Kirchhoff model is analyzed as a double perturbed saddle-point problem. Unique solvability of the continuous formulation is proved using Fredholm's theory for compact operators and the Babuška–Brezzi approach for saddle-point problems with penalty. The proposed stable virtual element method establishes a discrete inf-sup condition under a small mesh assumption through a Fortin interpolant that requires only $\widetilde{u}_1$2-regularity for the Stokes problem, and the monolithic discrete formulation is shown to be well posed with optimal convergence in the energy norm (Dassi et al., 4 Aug 2025).

The fully averaged poroelastic Kirchhoff plate coupled to Stokes flow is studied by energy methods at the weak-solution level and by sectoriality of the associated spatial operator together with maximal $\widetilde{u}_1$3-regularity for a regularized strong-solution problem. The paper proves existence of weak solutions, global-in-time existence of a unique strong solution to a regularized version of the problem, exponential decay of solutions for data with exponential decay, and develops a finite element method that provides an excellent approximation of the full Biot-Stokes system in the thin-structure regime (Brandt et al., 17 May 2026).

6. Distinctions from neighboring theories

A recurrent source of confusion is that not every Biot paper and not every Kirchhoff paper is about Biot–Kirchhoff equations. The high-frequency Biot-JKD model is a Biot-type poroelastic wave system with dynamic tortuosity and an augmented first-order PDE–ODE system that removes explicit memory convolutions; it is explicitly described as not using the phrase “Biot-Kirchhoff equations” and not being a Kirchhoff plate or shell theory (Ou et al., 2018).

Another distinct tradition is the classical Kirchhoff equations for rigid body motion in an ideal incompressible fluid. That literature concerns the Lie algebra $\widetilde{u}_1$4, the impulsive moment $\widetilde{u}_1$5, the impulsive force $\widetilde{u}_1$6, and Hamiltonian rigid-body dynamics, not Biot poroelasticity. It is therefore relevant only to the Kirchhoff half of the phrase and not to the poroelastic plate models surveyed here (Dragovic et al., 2011).

Standard quasi-static Biot poroelasticity, even when coupled to Stokes flow or equipped with advanced solvers and discretisations, is also not in itself a Biot–Kirchhoff theory. This includes stabilized finite element and iterative solver treatments of the classical Biot equations, total-pressure formulations for Biot–Stokes coupling, and stable cell-centered finite volume discretisations for the Biot equations. Those works address the Biot side of poroelasticity but do not introduce a Kirchhoff plate reduction (Cai et al., 2017, Ruiz-Baier et al., 2021, Nordbotten, 2015).

The most accurate use of the term in the available arXiv literature is therefore restricted to thin porous plates or thin poroelastic interfaces where one finds, in exact or derived form, a Biot-type pressure equation together with a Kirchhoff or Kirchhoff–Love bending law. Within that scope, the phrase covers the pressure-moment plate model, the clamped Biot–Kirchhoff–Love reduction, the homogenized Biot-plate-system of Kirchhoff-Love type, and recent bulk–surface or fully averaged interface formulations (Khot et al., 2023, Marciniak-Czochra et al., 2012, Gahn, 2024, Dassi et al., 4 Aug 2025, Brandt et al., 17 May 2026).