A fully averaged poroelastic Kirchhoff plate interacting with an incompressible, viscous fluid: analysis and numerical simulation
Pith reviewed 2026-05-20 12:50 UTC · model grok-4.3
The pith
A fully averaged poroelastic Kirchhoff plate model coupled to Stokes flow admits weak solutions and unique strong solutions for its regularized version.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the fully averaged poroelastic Kirchhoff plate equations posed on a codimension-one interface, when coupled to incompressible viscous fluid flow via continuity of normal velocities, Beavers-Joseph-Saffman slip, and force balance, possess weak solutions whose existence follows from energy estimates, while a regularized version of the coupled problem possesses a unique global strong solution whose existence follows from sectoriality of the spatial operator and maximal L^p-regularity of the resulting Cauchy problem; exponential decay holds for exponentially decaying data, and a surface finite-element method approximates the full Biot-Stokes system accurately in the薄结构限
What carries the argument
The fully averaged poroelastic Kirchhoff plate formulation placed on a codimension-one interface and coupled to the Stokes equations through kinematic and dynamic interface conditions.
If this is right
- Weak solutions exist for the linearly coupled fluid-structure system.
- A regularized version of the problem has a unique global-in-time strong solution.
- Solutions decay exponentially whenever the data decay exponentially.
- The surface finite-element scheme approximates the full Biot-Stokes system with high accuracy when the plate is thin.
Where Pith is reading between the lines
- The surface reduction may enable efficient simulation of flow through thin biological membranes or porous filters without meshing the plate interior.
- The sectoriality-plus-maximal-regularity technique could transfer to other thin-layer fluid-structure models that admit similar energy structures.
- Implementation simplicity on a single surface mesh may lower the barrier to exploring parameter studies in microfluidic or tissue-engineering contexts.
Load-bearing premise
The thin-structure regime in which averaging the poroelastic equations across the plate thickness produces an accurate reduction of the full three-dimensional Biot model.
What would settle it
A direct numerical comparison of the averaged surface model against the full three-dimensional Biot-Stokes system that reveals large discrepancies for plates of moderate thickness would falsify the practical accuracy claim.
Figures
read the original abstract
We study a new fully averaged poroelastic Kirchhoff plate model coupled with the flow of an incompressible, viscous fluid governed by the time-dependent Stokes equations. The fully averaged formulation offers several advantages over the classical Biot poroelastic plate model: both elastodynamic and pressure equations are posed on a codimension-one interface, the resulting numerical schemes are simpler to implement and computationally more efficient, and the fluid-structure coupling is more natural. We analyze a linearly coupled fluid-structure interaction problem with kinematic and dynamic interface conditions enforcing continuity of normal velocities, the Beavers-Joseph-Saffman slip in the tangential velocities, and balance of forces between the fluid and the poroelastic structure. We establish the existence of weak solutions using energy methods, and then prove global-in-time existence of a unique strong solution to a regularized version of the problem using sectoriality of the associated spatial operator and maximal $\mathrm{L}^p$-regularity of the resulting Cauchy problem. For data with exponential decay, we prove exponential decay of solutions. Finally, we develop a finite element method for the numerical approximation of the coupled system and show that it provides an excellent approximation of the full Biot-Stokes system in the thin-structure regime. The main advantage of this model lies in the remarkably simple implementation, as the poroelastic plate equations constitute a surface model bounding a bulk fluid domain. These results provide a rigorous analytical and computational framework for the study of coupled fluid-poroelastic structure interactions involving thin poroelastic interfaces modeled by the fully averaged Kirchhoff poroelastic plate equations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a fully averaged poroelastic Kirchhoff plate model coupled to the time-dependent Stokes equations for an incompressible viscous fluid. It establishes existence of weak solutions via energy methods, proves global-in-time existence and uniqueness of strong solutions for a regularized version of the problem via sectoriality of the spatial operator and maximal L^p-regularity, shows exponential decay for decaying data, and develops a finite element method that approximates the full Biot-Stokes system in the thin-structure regime.
Significance. If the results hold, the work supplies a simplified surface-based model for thin poroelastic fluid-structure interactions that retains analytical tractability while offering computational advantages over bulk Biot models. The proofs rely on standard energy estimates and maximal-regularity theory for coupled Stokes-structure systems, which are applied appropriately to the interface conditions (normal-velocity continuity, Beavers-Joseph-Saffman slip, force balance). The numerical demonstration of efficiency in the thin regime adds practical value, though it rests on the thin-structure approximation rather than a rigorous error analysis.
minor comments (3)
- The abstract states that the fully averaged formulation is 'computationally more efficient' and 'simpler to implement,' but the numerical section would benefit from explicit timing or degree-of-freedom comparisons with a standard Biot-Stokes discretization to quantify this advantage.
- In the section describing the regularized problem, the precise form of the regularization (e.g., added viscosity or smoothing of the interface) should be stated explicitly with an equation number so that readers can verify how it preserves the energy structure used for the weak-solution result.
- The claim of 'excellent approximation' of the full Biot-Stokes system in the thin-structure regime is supported only by numerical examples; adding a brief discussion of observed convergence rates or mesh-refinement studies would strengthen the computational section without altering the analytical core.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and the positive overall assessment. We appreciate the recommendation for minor revision and address the substantive points below.
read point-by-point responses
-
Referee: The numerical demonstration of efficiency in the thin regime adds practical value, though it rests on the thin-structure approximation rather than a rigorous error analysis.
Authors: We agree that the numerical section demonstrates the practical advantage of the fully averaged model by comparing it computationally to the full Biot-Stokes system in the thin limit, rather than supplying a rigorous a priori error estimate. Establishing such an estimate would require a separate, technically involved analysis that lies beyond the scope of the present work, which focuses on well-posedness of the new surface model and its numerical implementation. In the revised manuscript we have added a short paragraph in Section 6 clarifying this distinction and identifying a rigorous thin-limit error analysis as a natural direction for future research. revision: yes
Circularity Check
No significant circularity in derivation chain
full rationale
The derivation chain relies on standard energy methods for weak solutions and sectoriality plus maximal L^p-regularity for the regularized strong solutions, both of which are applied to the coupled Stokes-poroelastic system with explicitly stated interface conditions (normal velocity continuity, Beavers-Joseph-Saffman slip, force balance). These techniques are drawn from general PDE theory for fluid-structure interactions and do not reduce to self-definitions, fitted inputs renamed as predictions, or load-bearing self-citations. The exponential decay result follows directly from the energy structure, and the numerical finite-element claim is an approximation observation in the thin-structure regime rather than a derived prediction. The model is self-contained against external benchmarks of functional analysis.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The poroelastic plate is treated as a codimension-one interface with fully averaged elastodynamic and pressure equations.
- domain assumption Kinematic and dynamic interface conditions hold: continuity of normal velocities, Beavers-Joseph-Saffman slip, and force balance.
invented entities (1)
-
Fully averaged poroelastic Kirchhoff plate
no independent evidence
Reference graph
Works this paper leans on
-
[1]
A. Agresti and A. Hussein,MaximalL p-regularity andH ∞-calculus for block operator matrices and applications, J. Funct. Anal., 285 (2023), Paper No. 110146
work page 2023
-
[2]
Amann,Linear and Quasilinear Parabolic Problems
H. Amann,Linear and Quasilinear Parabolic Problems. Vol. I. Abstract Linear Theory, Monographs in Mathematics, vol. 89, Birkhäuser Boston, 1995
work page 1995
-
[3]
Amann,On the strong solvability of the Navier–Stokes equations, J
H. Amann,On the strong solvability of the Navier–Stokes equations, J. Math. Fluid Mech., 2 (2000), 16–98
work page 2000
-
[4]
I. Ambartsumyan, V. J. Ervin, T. Nguyen, and I. Yotov,A nonlinear Stokes–Biot model for the interaction of a non-Newtonian fluid with poroelastic media, ESAIM Math. Model. Numer. Anal., 53 (2019), 1915–1955
work page 2019
-
[5]
I. Ambartsumyan, E. Khattatov, I. Yotov, and P. Zunino,A Lagrange multiplier method for a Stokes–Biot fluid- poroelastic structure interaction problem, Numer. Math., 140 (2018), 513–553
work page 2018
-
[6]
J.-L. Auriault,Poroelastic Media. Homogenization and Porous Media, Interdiscip. Appl. Math., vol. 6, Springer, New York, 1997
work page 1997
- [7]
-
[8]
G. Avalos and J. T. Webster,Uniqueness of weak solutions for Biot-Stokes interactions, Pure Appl. Anal., 7 (2025), 1111–1139
work page 2025
- [9]
-
[10]
H. T. Banks, K. Bekele-Maxwell, L. Bociu, M. Noorman, and G. Guidoboni,Local sensitivity via the complex-step derivative approximation for 1D poro-elastic and poro-visco-elastic models, Math. Control Relat. Fields, 9 (2019), 623–642. FLUID-POROELASTIC STRUCTURE INTERACTION WITH A FULLY A VERAGED PLATE 53
work page 2019
- [11]
- [12]
- [13]
- [14]
-
[15]
T. Binz, M. Hieber, and A. Roy.fluid–structure interaction with porous media: the Beaver-Joseph condition in the strong sense, J. Differential Equations, 426 (2025), 660–689
work page 2025
-
[16]
M. A. Biot.Theory of elasticity and consolidation for a porous anisotropic solid, J. Appl. Phys., 26 (1955), 182–185
work page 1955
-
[17]
M. A. Biot.General theory of three-dimensional consolidation, J. Appl. Phys., 12 (1941), 155–164
work page 1941
- [18]
- [19]
- [20]
- [21]
- [22]
-
[23]
L. Bociu and M. Verri.Analysis of a coupled fluid–poroelastic structure interaction model, SIAM J. Appl. Math., 81 (2021), 2350–2375
work page 2021
-
[24]
L. Bociu and J. T. Webster.Nonlinear quasi-static poroelasticity, J. Differential Equations, 296 (2021), 242–278
work page 2021
- [25]
- [26]
- [27]
-
[28]
Bukač.Analysis of a fluid–structure interaction problem involving a poroelastic structure, Phys
M. Bukač.Analysis of a fluid–structure interaction problem involving a poroelastic structure, Phys. Rev. E, 94 (2016), Paper No. 033118
work page 2016
- [29]
- [30]
- [31]
- [32]
- [33]
- [34]
-
[35]
Cesmelioglu.Analysis of the coupled Navier–Stokes/Biot problem, J
A. Cesmelioglu.Analysis of the coupled Navier–Stokes/Biot problem, J. Math. Anal. Appl., 456 (2017), 970–991
work page 2017
-
[36]
Coussy.Poromechanics, John Wiley & Sons, New York, 2004
O. Coussy.Poromechanics, John Wiley & Sons, New York, 2004
work page 2004
- [37]
-
[38]
R. Denk, M. Hieber, and J. Prüss.R-boundedness, Fourier multipliers and problems of elliptic and parabolic type, Mem. Amer. Math. Soc., 166 (2003), no. 788
work page 2003
-
[39]
R. Denk, M. Hieber, and J. Prüss.OptimalLp-Lq-estimates for parabolic boundary value problems with inhomogeneous data, Math. Z., 257 (2007), 193–224
work page 2007
-
[40]
L. de Simon.Un’applicazione della teoria degli integrali singolari allo studio delle equazioni differenziali lineari astratte del prime ordine, Rend. Sem. Mat. Univ. Padova, 34 (1964), 205–223. [41]The FEniCS computing platform.https://fenicsproject.org/
work page 1964
-
[41]
J. L. Ferrín and A. Mikelić.Homogenizing the acoustic properties of a porous matrix containing an incompressible inviscid fluid, Math. Methods Appl. Sci., 26 (2003), 831–859. [43]The FreeFem computing platform.https://freefem.org/
work page 2003
-
[42]
Y. C. Fung.Biomechanics: Circulation, Springer, New York, 1997
work page 1997
-
[43]
C. Grandmont, M. Hillairet, and J. Lequeurre.Existence of local strong solutions to fluid-beam and fluid-rod interaction systems, Ann. Inst. H. Poincaré C Anal. Non Linéaire, 36 (2019), 1105–1149. 54 FELIX BRANDT, SUNČICA ČANIĆ, ANDREW SCHARF, AND JOSIP TAMBAČA
work page 2019
-
[44]
M. Hieber and T. Kashiwabara.Global strong well-posedness of the three-dimensional primitive equations inLp-spaces, Arch. Ration. Mech. Anal., 221 (2016), 1077–1115
work page 2016
-
[45]
M. Hieber and J. Saal.The Stokes equation in theLp-setting: well-posedness and regularity properties, In: Handbook of Mathematical Analysis in Mechanics of Viscous Fluids, Y. Giga and A. Novotný, eds., Springer, Cham, 2018, 117–206
work page 2018
-
[46]
J. M. Huyghe, T. Arts, D. H. van Campen, and R. S. Reneman.Porous medium finite element model of the beating left ventricle, Amer. J. Physiol., 262 (1992), 1256–1267
work page 1992
-
[47]
O. P. Iliev, A. E. Kolesov, and P. N. Vabishchevich.Numerical solution of plate poroelasticity problems, Transport in Porous Media, 115 (2016), 563–580
work page 2016
-
[48]
T. Kato.Perturbation Theory for Linear Operators, Reprint of the 1980 edition, Classics in Mathematics, Springer- Verlag, Berlin, 1995
work page 1980
-
[49]
R. Khot and R. Ruiz Baier.Virtual element methods for Biot-Kirchhoff poroelasticity, Math. Comp., 94 (2025), 1101– 1146
work page 2025
-
[50]
D. J. Korteweg.Über die Fortpflanzungsgeschwindigkeit des Schalles in elastischen Röhren, Annalen der Physik, 241 (1878), 525–542
-
[51]
J.Kuan, S.Čanić, andB.Muha.Fluid-poroviscoelastic structure interaction problem with nonlinear geometric coupling, J. Math. Pures Appl., 188 (2024), 345–445
work page 2024
-
[52]
J. Kuan, S. Čanić, and B. Muha.Existence of a weak solution to a regularized moving boundary fluid–structure interaction problem with poroelastic media. Comptes Rendus Mécanique, 351 (2023), 1–30
work page 2023
- [53]
-
[54]
P. Kunštek, M. Bukač, and B. Muha.Mass conservation in the validation of fluid-poroelastic structure interaction solvers, Appl. Math. Comput., 487 (2025), Paper No. 129081
work page 2025
-
[55]
P. C. Kunstmann and L. Weis.Perturbation theorem for maximalLp-regularity, Ann. Scuola Norm. Sup. Pisa Cl. Sci., 30 (2001), 415–435
work page 2001
-
[56]
A. Marciniak-Czochra and A. Mikelić.A rigorous derivation of the equations for the clamped Biot-Kirchhoff-Love poroelastic plate, Arch. Ration. Mech. Anal., 215 (2015), 1035–1062
work page 2015
-
[57]
A. Mikelić and J. Tambača.Derivation of a poroelastic flexural shell model.Multiscale Model. Simul., 14 (2016), 364–397
work page 2016
-
[58]
A. Mikelić and M. F. Wheeler.On the interface law between a deformable porous medium containing a viscous fluid and an elastic body, Math. Models Methods Appl. Sci., 22 (2012), Paper No. 1250031
work page 2012
-
[59]
P. Miotto and M. Bukač.A fluid–structure interaction problem with a semi-permeable poroelastic material, Com- put. Methods Appl. Mech. Engrg., 384 (2021), 113941
work page 2021
-
[60]
A. I. Moens.Die Pulskurve, E. J. Brill, Leiden, 1878
-
[61]
T. Nau.L p-Theory of Cylindrical Boundary Value Problems, PhD thesis, University of Konstanz, Springer Spektrum, 2012
work page 2012
-
[62]
Nau.TheL p-Helmholtz projection in finite cylinders, Czechoslovak Math
T. Nau.TheL p-Helmholtz projection in finite cylinders, Czechoslovak Math. J., 65 (2015), 119–134
work page 2015
-
[63]
Owczarek.A Galerkin method for Biot consolidation model, Math
S. Owczarek.A Galerkin method for Biot consolidation model, Math. Mech. Solids, 15 (2010), 42–56
work page 2010
-
[64]
O. Oyekola and M. Bukač.Second-order, loosely coupled methods for fluid-poroelastic material interaction, Nu- mer. Methods Partial Differential Equations, 36 (2020), 800–822
work page 2020
-
[65]
C. Parrow and M. Bukač.A Robin-Robin strongly coupled partitioned method for fluid-poroelastic structure interaction, J. Numer. Math., 33 (2025), 289–312
work page 2025
-
[66]
J. Prüss and G. Simonett.Moving Interfaces and Quasilinear Parabolic Evolution Equations, Monographs in Mathe- matics, vol. 105, Birkhäuser, 2016
work page 2016
-
[67]
E. Rohan and S. Naili.Homogenization of the fluid–structure interaction in acoustics of porous media perfused by viscous fluid, Z. Angew. Math. Phys., 71 (2020), Paper No. 137
work page 2020
- [68]
- [69]
-
[70]
R. E. Showalter.Diffusion in poro-elastic media, J. Math. Anal. Appl., 251 (2000), 310–340
work page 2000
-
[71]
R. E. Showalter.Poroelastic filtration coupled to Stokes flow, In: Lect. Notes Pure Appl. Math., Vol. 242, O. Imanuvilov et al., eds., Chapman & Hall/CRC, 2005, 229–241
work page 2005
-
[72]
R. E. Showalter and N. Su.Partially saturated flow in a poroelastic medium, Discrete Contin. Dyn. Syst. Ser. B, 1 (2001), 403–420
work page 2001
-
[73]
Terzaghi.Theoretical Soil Mechanics, John Wiley & Sons, 1943
K. Terzaghi.Theoretical Soil Mechanics, John Wiley & Sons, 1943
work page 1943
-
[74]
Tolksdorf.R-sectoriality of higher-order elliptic systems on general bounded domains, J
P. Tolksdorf.R-sectoriality of higher-order elliptic systems on general bounded domains, J. Evol. Equ., 18 (2018), 323–349
work page 2018
-
[75]
Tretter.Spectral Theory of Block Operator Matrices and Applications, Imperial College Press, 2008
C. Tretter.Spectral Theory of Block Operator Matrices and Applications, Imperial College Press, 2008
work page 2008
-
[76]
Triebel.Interpolation Theory, Function Spaces, Differential Operators, North-Holland, 1978
H. Triebel.Interpolation Theory, Function Spaces, Differential Operators, North-Holland, 1978
work page 1978
- [77]
-
[78]
Y. Wang, S. Čanić, M. Bukač, C. Blaha, and S. Roy.Mathematical and computational modeling of a poroelastic cell scaffold in a bioartificial pancreas, Fluids, 7 (2022), Paper No. 222. FLUID-POROELASTIC STRUCTURE INTERACTION WITH A FULLY A VERAGED PLATE 55
work page 2022
-
[79]
Weis.Operator-valued Fourier multiplier theorems and maximalLp-regularity, Math
L. Weis.Operator-valued Fourier multiplier theorems and maximalLp-regularity, Math. Ann., 319 (2001), 735–758
work page 2001
-
[80]
A. Ženíšek.The existence and uniqueness theorem in Biot’s consolidation theory, Aplikace Matematiky, 29 (1984), 194–211. Department of Mathematics, University of California at Berkeley, Berkeley, 94720, CA, USA. Email address:fbrandt@berkeley.edu Department of Mathematics, University of California at Berkeley, Berkeley, 94720, CA, USA. Email address:canic...
work page 1984
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.