An unconditionally stable second-order Robin partitioned method for fluid-poroelastic structure interaction
Abstract
We study an unsteady fluid-poroelastic structure interaction (FPSI) problem, in which the free fluid region is governed by the incompressible Stokes equations and the poroelastic region by the fully dynamic Biot system. The original physical coupling conditions are equivalently reformulated as Robin-type interface conditions and an auxiliary interface variable is introduced as common Robin data shared by the two subsystems. This reformulation decomposes the coupled FPSI problem into fluid and poroelastic subproblems posed on their respective subdomains. The problem is discretized in space by the finite element method (FEM), while the backward Euler scheme and the BDF2 scheme combined with the second-order interface extrapolation are employed for the temporal discretization. By exploiting the resulting discrete dynamic relations, we derive an update formula for the auxiliary interface variable and obtain fully discrete first- and second-order Robin partitioned schemes in which the two subproblems are solved independently at each time step. The unconditional stability of both schemes is rigorously established. Numerical experiments verify the temporal convergence orders of the proposed schemes and investigate the influence of the Robin parameter. A classical pressure-wave benchmark verifies the accuracy of the proposed methods, while supplementary nonlinear moving-domain simulations illustrate their applicability to more complex FPSI configurations.
Related
- An arbitrary-order time-sub-cycling framework for conservative partitioned fluid–structure interaction
- High-Order Unconditionally Stable IMEX BDF--FEM Schemes for the Hodgkin--Huxley Model on a Branched Axon
- Unconditionally L2-stable explicit--implicit--null filtered Runge--Kutta methods for general nonlinear degenerate parabolic equations
- A linear, second-order, unconditionally energy-stable scalar-auxiliary-variable scheme for the CHQ viscoelastic phase-separation model
- Weak Approximation Properties Induced by Unconditionally p-Compact Sets
BibTeX
@misc{he2026unconditionally,
title = {{An unconditionally stable second-order Robin partitioned method for fluid-poroelastic structure interaction}},
author = {Wenlong He and Thomas Wick and Xiaohe Yue and Jiwei Zhang and Haibiao Zheng},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.25350},
eprint = {2609.25350},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.25350},
}
RIS
TY - UNPB
TI - An unconditionally stable second-order Robin partitioned method for fluid-poroelastic structure interaction
AU - Wenlong He
AU - Thomas Wick
AU - Xiaohe Yue
AU - Jiwei Zhang
AU - Haibiao Zheng
PY - 2026
DA - 2026/09/21
PB - arXiv
DO - 10.48550/arxiv.2609.25350
UR - https://arxiv.org/abs/2609.25350
AB - We study an unsteady fluid-poroelastic structure interaction (FPSI) problem, in which the free fluid region is governed by the incompressible Stokes equations and the poroelastic region by the fully dynamic Biot system. The original physical coupling conditions are equivalently reformulated as Robin-type interface conditions and an auxiliary interface variable is introduced as common Robin data shared by the two subsystems. This reformulation decomposes the coupled FPSI problem into fluid and poroelastic subproblems posed on their respective subdomains. The problem is discretized in space by the finite element method (FEM), while the backward Euler scheme and the BDF2 scheme combined with the second-order interface extrapolation are employed for the temporal discretization. By exploiting the resulting discrete dynamic relations, we derive an update formula for the auxiliary interface variable and obtain fully discrete first- and second-order Robin partitioned schemes in which the two subproblems are solved independently at each time step. The unconditional stability of both schemes is rigorously established. Numerical experiments verify the temporal convergence orders of the proposed schemes and investigate the influence of the Robin parameter. A classical pressure-wave benchmark verifies the accuracy of the proposed methods, while supplementary nonlinear moving-domain simulations illustrate their applicability to more complex FPSI configurations.
ER -