A Nitsche method for Navier--Stokes/generalized poroelasticity interface problems

Aparna Bansal; Dwijendra Narain Pandey; Nicolas A. Barnafi; Ricardo Ruiz-Baier

arxiv: 2508.09388 · v2 · pith:VG4DTPUOnew · submitted 2025-08-12 · 🧮 math.NA · cs.NA

A Nitsche method for Navier--Stokes/generalized poroelasticity interface problems

Aparna Bansal , Nicolas A. Barnafi , Dwijendra Narain Pandey , Ricardo Ruiz-Baier This is my paper

Pith reviewed 2026-05-21 23:34 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords Nitsche methodNavier-Stokes equationsporoelasticityfinite element methodinterface problemsmonolithic discretizationerror estimates

0 comments

The pith

A Nitsche formulation yields a stable monolithic finite element scheme for Navier-Stokes and generalized poroelasticity interface problems.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a unified monolithic finite element discretization for the time-dependent coupled Navier-Stokes and generalized poroelastic flow problem. It handles the fluid-structure interface through a Nitsche-type weak enforcement that avoids the need for conforming meshes. Well-posedness of the discrete problem follows from differential-algebraic equations theory combined with the Banach fixed-point theorem, while stability and a priori error estimates hold for a penalty parameter chosen independently of mesh size. This setup produces convergent approximations that capture the coupled dynamics without mesh-dependent tuning.

Core claim

The authors construct a monolithic implicit time-stepping finite element scheme that incorporates a Nitsche-type formulation to enforce the interface conditions between the Navier-Stokes fluid and the generalized poroelastic medium. The resulting discrete system is shown to be well-posed, stable, and convergent, with the stability and error bounds controlled by a properly chosen penalty parameter that remains independent of the mesh size.

What carries the argument

Nitsche-type formulation, a technique for weakly enforcing interface conditions between fluid and poroelastic domains within a single monolithic discretization.

If this is right

The fully discrete scheme remains well-posed for arbitrary mesh sizes.
Stability is obtained with a single penalty parameter that does not vary with refinement.
A priori error estimates hold uniformly for the velocity, pressure, and displacement fields.
The method reproduces the coupled dynamics accurately in the presented numerical tests.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The mesh-independent penalty choice may reduce the need for parameter retuning when the same code is applied to different material regimes.
The monolithic structure could simplify the design of adaptive time-stepping strategies for long-time simulations of the interface problem.
Similar Nitsche treatments might be tested on other multiphysics couplings that mix incompressible flow with deformable porous media.

Load-bearing premise

The penalty parameter can be chosen once to guarantee stability and convergence for the entire coupled system without depending on mesh size or requiring hidden restrictions on material parameters.

What would settle it

Numerical experiments in which the observed convergence rate falls below the predicted order when the mesh is successively refined while the penalty parameter is held fixed would falsify the stability and error claims.

Figures

Figures reproduced from arXiv: 2508.09388 by Aparna Bansal, Dwijendra Narain Pandey, Nicolas A. Barnafi, Ricardo Ruiz-Baier.

**Figure 6.** Figure 6 [PITH_FULL_IMAGE:figures/full_fig_p021_6.png] view at source ↗

read the original abstract

We consider a time-dependent coupled Navier--Stokes/generalized poroelastic flow problem and propose a unified and monolithic finite element discretization based on implicit time stepping. To handle the fluid-structure interface we employ a Nitsche-type formulation. The resulting discrete problem is shown to be well-posed using the theory of differential-algebraic equations (DAEs) and the Banach fixed-point theorem. We prove stability and derive a priori error estimates for the fully discrete scheme. The stability and convergence of the method are ensured by a properly chosen penalty parameter independent of the mesh size. Numerical tests are presented to confirm the theoretical convergence rates and to illustrate the ability of the method to capture the coupled dynamics accurately.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This extends Nitsche coupling to a monolithic Navier-Stokes/generalized poroelasticity system with DAE well-posedness and error estimates; the mesh-independent penalty claim is plausible but needs explicit checks against material coefficients.

read the letter

The main contribution is a monolithic finite element scheme that couples Navier-Stokes to generalized poroelasticity via Nitsche terms at the interface, using implicit time stepping throughout. Well-posedness of the semi-discrete problem follows from DAE theory, the convective nonlinearity is handled by a Banach fixed-point argument, and they derive stability plus a priori error estimates for the fully discrete version. Numerical tests back the convergence rates and show the coupled dynamics working in practice. This is a direct, honest extension of existing interface techniques rather than a new paradigm, and the choice to treat the whole thing monolithically avoids splitting errors that often appear in these problems. The citation pattern draws on standard Nitsche and poroelasticity references without circularity or self-referential fitting. The analysis looks formally grounded in the tools they invoke. One soft spot is the penalty parameter. The abstract asserts a choice independent of mesh size that guarantees stability and convergence, but Nitsche lower bounds routinely absorb viscosity, permeability, Lamé parameters, and sometimes the time step through the contraction constant. If those dependencies are not ruled out explicitly or if the bound requires hidden restrictions on the coefficients, the claim would need qualification. The paper does not appear to introduce invented entities or load-bearing assumptions that contradict its own equations. This is aimed at computational mechanicians who need a workable tool for fluid-poroelastic interface problems in applications such as tissue mechanics or reservoir simulation. Readers who already work with monolithic schemes or DAE-based analysis will find the error estimates and tests directly usable. It has enough formal structure and reproducible numerics to deserve peer review rather than a desk reject.

Referee Report

2 major / 2 minor

Summary. The paper proposes a monolithic Nitsche-type finite element discretization with implicit time stepping for the time-dependent Navier-Stokes equations coupled to a generalized poroelasticity model across an interface. Well-posedness of the semi-discrete system is established via differential-algebraic equation (DAE) theory combined with the Banach fixed-point theorem to treat the convective nonlinearity; stability and a priori error estimates are derived for the fully discrete scheme, with the Nitsche penalty parameter asserted to be independent of mesh size. Numerical experiments are included to verify convergence rates and illustrate coupled dynamics.

Significance. If the well-posedness, stability, and error analysis hold with a truly mesh-independent penalty, the work supplies a rigorous, unfitted interface treatment for a challenging fluid-poroelastic coupling that appears in biomechanics and reservoir modeling. The combination of DAE theory for the monolithic system and explicit a priori estimates would strengthen the literature on monolithic FSI schemes.

major comments (2)

[§3 and abstract] §3 (well-posedness analysis) and the statement in the abstract: the lower bound on the Nitsche penalty parameter γ is claimed to be independent of mesh size, yet the standard inverse-trace estimates in the Nitsche consistency terms necessarily absorb the viscosity μ, permeability κ, Biot-Willis coefficient, and Lamé parameters. The DAE index-2 structure and the contraction constant of the Banach fixed-point map may further introduce explicit dependence on the time-step size Δt; the manuscript must exhibit the precise lower bound (or prove its independence from Δt and material ratios) to substantiate the central stability claim.
[§4] Theorem on a priori error estimates (likely §4): the error bound is stated to hold for a properly chosen penalty independent of h, but the proof sketch does not clarify whether the hidden constants remain uniform when the fluid-to-solid parameter contrast becomes large or when Δt is not sufficiently small relative to the material coefficients.

minor comments (2)

[Numerical experiments] The numerical section should tabulate the specific values of γ used for each mesh size and material parameter set to allow direct verification of mesh independence.
[§2] Notation for the generalized poroelasticity model (storage coefficient, Biot-Willis tensor) should be introduced once in §2 and used consistently thereafter.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major comment below and indicate the changes planned for the revised version.

read point-by-point responses

Referee: [§3 and abstract] §3 (well-posedness analysis) and the statement in the abstract: the lower bound on the Nitsche penalty parameter γ is claimed to be independent of mesh size, yet the standard inverse-trace estimates in the Nitsche consistency terms necessarily absorb the viscosity μ, permeability κ, Biot-Willis coefficient, and Lamé parameters. The DAE index-2 structure and the contraction constant of the Banach fixed-point map may further introduce explicit dependence on the time-step size Δt; the manuscript must exhibit the precise lower bound (or prove its independence from Δt and material ratios) to substantiate the central stability claim.

Authors: We thank the referee for highlighting this point. In the well-posedness analysis of Section 3 the lower bound on γ is obtained after absorbing all consistency terms via inverse-trace inequalities; the resulting threshold is independent of the mesh size h by construction. The bound does, however, depend on the physical coefficients (μ, κ, Lamé parameters, Biot-Willis coefficient) and, through the DAE index-2 theory and the contraction mapping constant, may also depend on Δt. We will revise Theorem 3.1 to display the explicit lower bound in terms of these quantities, add a remark clarifying the independence from h, and update the abstract and the statement in Section 3 accordingly. This makes the dependence transparent while preserving the central claim that the penalty need not grow with mesh refinement. revision: yes
Referee: [§4] Theorem on a priori error estimates (likely §4): the error bound is stated to hold for a properly chosen penalty independent of h, but the proof sketch does not clarify whether the hidden constants remain uniform when the fluid-to-solid parameter contrast becomes large or when Δt is not sufficiently small relative to the material coefficients.

Authors: We agree that the dependence of the hidden constants should be stated more explicitly. The a priori estimates in Section 4 are derived under the assumption that γ satisfies the lower bound from the well-posedness result; the constants are independent of h but may grow with the material-parameter contrast and with 1/Δt. We will insert a remark after the main error theorem noting this dependence and emphasizing that, for any fixed set of physical coefficients and time-step size, the estimates remain uniform with respect to h. No alteration of the theorem statement itself is required, but the proof sketch will be expanded to indicate where the parameter dependence enters. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation relies on external DAE theory and fixed-point arguments

full rationale

The paper's central well-posedness argument invokes standard theory of differential-algebraic equations for the semi-discrete system and the Banach fixed-point theorem to treat the convective term. Stability and a priori estimates are derived from this framework with a penalty parameter selected to satisfy a mesh-independent lower bound via standard inverse and trace inequalities. No load-bearing step reduces the main result to a self-definition, a fitted input renamed as prediction, or a self-citation chain; the cited mathematical tools are independent of the present construction and the penalty choice is not shown to be tautological within the paper's own equations.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard mathematical tools for DAEs and fixed-point arguments plus a mesh-independent penalty parameter whose selection is not derived from first principles within the abstract.

free parameters (1)

Nitsche penalty parameter
Chosen to ensure stability and convergence independent of mesh size; its specific value or selection rule is not derived in the abstract.

axioms (2)

standard math Theory of differential-algebraic equations applies to the semi-discrete system
Invoked to establish well-posedness of the discrete problem.
standard math Banach fixed-point theorem guarantees unique solution
Used for existence and uniqueness after the DAE reduction.

pith-pipeline@v0.9.0 · 5662 in / 1385 out tokens · 39246 ms · 2026-05-21T23:34:09.909870+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

[1]

Ambartsumyan, V

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), pp. 1915–1955

work page 2019
[2]

Ambartsumyan, E

I. Ambartsumyan, E. Khattatov, T. Nguyen, and I. Yotov, Flow and transport in fractured poroelastic media, GEM Int. J. Geomath., 10 (2019), pp. Paper No. 11, 34

work page 2019
[3]

Ambartsumyan, E

I. Ambartsumyan, E. Khattatov, I. Yotov, and P. Zunino, A Lagrange multiplier method for a Stokes-Biot fluid-poroelastic structure interaction model , Numer. Math., 140 (2018), pp. 513–553

work page 2018
[4]

Badia, A

S. Badia, A. Quaini, and A. Quarteroni , Coupling Biot and Navier-Stokes equations for modelling fluid-poroelastic media interaction , J. Comput. Phys., 228 (2009), pp. 7986– 8014

work page 2009
[5]

Badia and F

S. Badia and F. Verdugo , Gridap: an extensible finite element toolbox in julia , Journal of Open Source Software, 5 (2020), p. 2520

work page 2020
[6]

Bansal, N

A. Bansal, N. A. Barnafi, D. N. Pandey, and R. Ruiz-Baier , A lagrange multiplier- based method for stokes-linearized poro-hyperelastic interface problems , arXiv preprint arXiv:2407.13684, (2024)

work page arXiv 2024
[7]

Barnafi, P

N. Barnafi, P. Zunino, L. Ded`e, and A. Quarteroni, Mathematical analysis and numerical approximation of a general linearized poro-hyperelastic model , Comput. Math. Appl., 91 (2021), pp. 202–228

work page 2021
[8]

Bazilevs, M.-C

Y. Bazilevs, M.-C. Hsu, D. J. Benson, S. Sankaran, and A. L. Marsden , Computational fluid-structure interaction: methods and application to a total cavopulmonary connection , Comput. Mech., 45 (2009), pp. 77–89

work page 2009
[9]

Boffi, F

D. Boffi, F. Brezzi, and M. Fortin , Mixed finite element methods and applications , vol. 44 of Springer Series in Computational Mathematics, Springer, Heidelberg, 2013

work page 2013
[10]

K. E. Brenan, S. L. Campbell, and L. R. Petzold, Numerical solution of initial value prob- lems in differential-algebraic equations , North-Holland Publishing Co., New York, 1989

work page 1989
[11]

S. C. Brenner and L. R. Scott , The mathematical theory of finite element methods , vol. 15 of Texts in Applied Mathematics, Springer, New York, third ed., 2008

work page 2008
[12]

Bukaˇc, I

M. Bukaˇc, I. Yotov, R. Zakerzadeh, and P. Zunino, Partitioning strategies for the interac- tion of a fluid with a poroelastic material based on a Nitsche’s coupling approach , Comput. Methods Appl. Mech. Engrg., 292 (2015), pp. 138–170

work page 2015
[13]

Bukaˇc, I

M. Bukaˇc, I. Yotov, and P. Zunino , An operator splitting approach for the interaction be- tween a fluid and a multilayered poroelastic structure, Numer. Methods Partial Differential Equations, 31 (2015), pp. 1054–1100

work page 2015
[14]

Burman and M

E. Burman and M. A. Fern´andez, Stabilization of explicit coupling in fluid-structure interac- tion involving fluid incompressibility , Comput. Methods Appl. Mech. Engrg., 198 (2009), pp. 766–784

work page 2009
[15]

, Explicit strategies for incompressible fluid-structure interaction problems: Nitsche type mortaring versus Robin-Robin coupling , Internat. J. Numer. Methods Engrg., 97 (2014), pp. 739–758

work page 2014
[16]

Burtschell, P

B. Burtschell, P. Moireau, and D. Chapelle, Numerical analysis for an energy-stable total discretization of a poromechanics model with inf-sup stability, Acta Math. Appl. Sin. Engl. Ser., 35 (2019), pp. 28–53

work page 2019
[17]

Y. Cao, M. Gunzburger, F. Hua, and X. Wang, Coupled Stokes-Darcy model with Beavers- Joseph interface boundary condition , Commun. Math. Sci., 8 (2010), pp. 1–25

work page 2010
[18]

Caucao, T

S. Caucao, T. Li, and I. Yotov, A multipoint stress-flux mixed finite element method for the Stokes-Biot model, Numer. Math., 152 (2022), pp. 411–473

work page 2022
[19]

Cesmelioglu and P

A. Cesmelioglu and P. Chidyagwai, Numerical analysis of the coupling of free fluid with a poroelastic material, Numer. Methods Partial Differential Equations, 36 (2020), pp. 463– NITSCHE METHODS FOR NS/POROELASTICITY INTERFACE PROBLEMS 23 494

work page 2020
[20]

Cesmelioglu, J

A. Cesmelioglu, J. J. Lee, and S. Rhebergen, Hybridizable discontinuous Galerkin methods for the coupled Stokes-Biot problem , Comput. Math. Appl., 144 (2023), pp. 12–33

work page 2023
[21]

D’Angelo and P

C. D’Angelo and P. Zunino , Numerical approximation with Nitsche’s coupling of transient Stokes’/Darcy’s flow problems applied to hemodynamics , Appl. Numer. Math., 62 (2012), pp. 378–395

work page 2012
[22]

Discacciati, E

M. Discacciati, E. Miglio, and A. Quarteroni , Mathematical and numerical models for coupling surface and groundwater flows , vol. 43, 2002, pp. 57–74

work page 2002
[23]

M. A. Fern ´andez, Incremental displacement-correction schemes for the explicit coupling of a thin structure with an incompressible fluid , C. R. Math. Acad. Sci. Paris, 349 (2011), pp. 473–477

work page 2011
[24]

G. P. Galdi and R. Rannacher, eds., Fundamental trends in fluid-structure interaction, vol. 1 of Contemporary Challenges in Mathematical Fluid Dynamics and Its Applications, World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2010

work page 2010
[25]

Galvis and M

J. Galvis and M. Sarkis, Non-matching mortar discretization analysis for the coupling Stokes- Darcy equations, Electron. Trans. Numer. Anal., 26 (2007), pp. 350–384

work page 2007
[26]

G. N. Gatica, R. Oyarz´ua, and F.-J. Sayas, Analysis of fully-mixed finite element methods for the Stokes-Darcy coupled problem , Math. Comp., 80 (2011), pp. 1911–1948

work page 2011
[27]

Z. Ge, J. Pang, and J. Cao , Multiphysics mixed finite element method with Nitsche’s tech- nique for Stokes-poroelasticity problem, Numer. Methods Partial Differential Equations, 39 (2023), pp. 544–576

work page 2023
[28]

S. Guo, Y. Sun, Y. Wang, X. Yue, and H. Zheng , A Fully Parallelizable Loosely Coupled Scheme for Fluid-Poroelastic Structure Interaction Problems , SIAM J. Sci. Comput., 47 (2025), pp. B951–B975

work page 2025
[29]

Kashiwabara, I

T. Kashiwabara, I. Oikawa, and G. Zhou , Penalty method with Crouzeix-Raviart approx- imation for the Stokes equations under slip boundary condition , ESAIM Math. Model. Numer. Anal., 53 (2019), pp. 869–891

work page 2019
[30]

J. Kuan, S. ˇCani´c, and B. Muha , Fluid-poroviscoelastic structure interaction problem with nonlinear geometric coupling, J. Math. Pures Appl. (9), 188 (2024), pp. 345–445

work page 2024
[31]

Nitsche, ¨Uber ein Variationsprinzip zur L¨ osung von Dirichlet-Problemen bei Verwendung von Teilr¨ aumen, die keinen Randbedingungen unterworfen sind, Abh

J. Nitsche, ¨Uber ein Variationsprinzip zur L¨ osung von Dirichlet-Problemen bei Verwendung von Teilr¨ aumen, die keinen Randbedingungen unterworfen sind, Abh. Math. Sem. Univ. Hamburg, 36 (1971), pp. 9–15

work page 1971
[32]

Quarteroni and A

A. Quarteroni and A. Valli , Numerical approximation of partial differential equations , vol. 23 of Springer Series in Computational Mathematics, Springer-Verlag, Berlin, 1994

work page 1994
[33]

Ruiz-Baier, M

R. Ruiz-Baier, M. Taffetani, H. D. Westermeyer, and I. Yotov, The Biot-Stokes coupling using total pressure: formulation, analysis and application to interfacial flow in the eye , Comput. Methods Appl. Mech. Engrg., 389 (2022), pp. Paper No. 114384, 30

work page 2022
[34]

L. R. Scott and S. Zhang , Finite element interpolation of nonsmooth functions satisfying boundary conditions, Math. Comp., 54 (1990), pp. 483–493

work page 1990
[35]

Seboldt, O

A. Seboldt, O. Oyekole, J. Tamba ˇca, and M. Buka ˇc, Numerical modeling of the fluid- porohyperelastic structure interaction, SIAM J. Sci. Comput., 43 (2021), pp. A2923–A2948

work page 2021
[36]

R. E. Showalter, Poroelastic filtration coupled to Stokes flow , in Control theory of partial differential equations, vol. 242 of Lect. Notes Pure Appl. Math., Chapman & Hall/CRC, Boca Raton, FL, 2005, pp. 229–241

work page 2005
[37]

J. Yu, Y. Sun, F. Shi, and H. Zheng , Nitsche’s type stabilized finite element method for the fully mixed Stokes-Darcy problem with Beavers-Joseph conditions , Appl. Math. Lett., 110 (2020), pp. 106588, 8

work page 2020
[38]

Yu and Y

J. Yu and Y. Zhang , Nitsche’s type stabilization for the fully mixed Navier-Stokes/Darcy problem, J. Appl. Anal. Comput., 11 (2021), pp. 1481–1493. Appendix A. Proof of Theorem 5.3. Proof. We introduce the errors for all variables and split them into approximation and discretization errors: en f := uS,n f − uS,n f,h = (uS,n f − If,huS,n f ) + (If,huS,n f...

work page 2021

[1] [1]

Ambartsumyan, V

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), pp. 1915–1955

work page 2019

[2] [2]

Ambartsumyan, E

I. Ambartsumyan, E. Khattatov, T. Nguyen, and I. Yotov, Flow and transport in fractured poroelastic media, GEM Int. J. Geomath., 10 (2019), pp. Paper No. 11, 34

work page 2019

[3] [3]

Ambartsumyan, E

I. Ambartsumyan, E. Khattatov, I. Yotov, and P. Zunino, A Lagrange multiplier method for a Stokes-Biot fluid-poroelastic structure interaction model , Numer. Math., 140 (2018), pp. 513–553

work page 2018

[4] [4]

Badia, A

S. Badia, A. Quaini, and A. Quarteroni , Coupling Biot and Navier-Stokes equations for modelling fluid-poroelastic media interaction , J. Comput. Phys., 228 (2009), pp. 7986– 8014

work page 2009

[5] [5]

Badia and F

S. Badia and F. Verdugo , Gridap: an extensible finite element toolbox in julia , Journal of Open Source Software, 5 (2020), p. 2520

work page 2020

[6] [6]

Bansal, N

A. Bansal, N. A. Barnafi, D. N. Pandey, and R. Ruiz-Baier , A lagrange multiplier- based method for stokes-linearized poro-hyperelastic interface problems , arXiv preprint arXiv:2407.13684, (2024)

work page arXiv 2024

[7] [7]

Barnafi, P

N. Barnafi, P. Zunino, L. Ded`e, and A. Quarteroni, Mathematical analysis and numerical approximation of a general linearized poro-hyperelastic model , Comput. Math. Appl., 91 (2021), pp. 202–228

work page 2021

[8] [8]

Bazilevs, M.-C

Y. Bazilevs, M.-C. Hsu, D. J. Benson, S. Sankaran, and A. L. Marsden , Computational fluid-structure interaction: methods and application to a total cavopulmonary connection , Comput. Mech., 45 (2009), pp. 77–89

work page 2009

[9] [9]

Boffi, F

D. Boffi, F. Brezzi, and M. Fortin , Mixed finite element methods and applications , vol. 44 of Springer Series in Computational Mathematics, Springer, Heidelberg, 2013

work page 2013

[10] [10]

K. E. Brenan, S. L. Campbell, and L. R. Petzold, Numerical solution of initial value prob- lems in differential-algebraic equations , North-Holland Publishing Co., New York, 1989

work page 1989

[11] [11]

S. C. Brenner and L. R. Scott , The mathematical theory of finite element methods , vol. 15 of Texts in Applied Mathematics, Springer, New York, third ed., 2008

work page 2008

[12] [12]

Bukaˇc, I

M. Bukaˇc, I. Yotov, R. Zakerzadeh, and P. Zunino, Partitioning strategies for the interac- tion of a fluid with a poroelastic material based on a Nitsche’s coupling approach , Comput. Methods Appl. Mech. Engrg., 292 (2015), pp. 138–170

work page 2015

[13] [13]

Bukaˇc, I

M. Bukaˇc, I. Yotov, and P. Zunino , An operator splitting approach for the interaction be- tween a fluid and a multilayered poroelastic structure, Numer. Methods Partial Differential Equations, 31 (2015), pp. 1054–1100

work page 2015

[14] [14]

Burman and M

E. Burman and M. A. Fern´andez, Stabilization of explicit coupling in fluid-structure interac- tion involving fluid incompressibility , Comput. Methods Appl. Mech. Engrg., 198 (2009), pp. 766–784

work page 2009

[15] [15]

, Explicit strategies for incompressible fluid-structure interaction problems: Nitsche type mortaring versus Robin-Robin coupling , Internat. J. Numer. Methods Engrg., 97 (2014), pp. 739–758

work page 2014

[16] [16]

Burtschell, P

B. Burtschell, P. Moireau, and D. Chapelle, Numerical analysis for an energy-stable total discretization of a poromechanics model with inf-sup stability, Acta Math. Appl. Sin. Engl. Ser., 35 (2019), pp. 28–53

work page 2019

[17] [17]

Y. Cao, M. Gunzburger, F. Hua, and X. Wang, Coupled Stokes-Darcy model with Beavers- Joseph interface boundary condition , Commun. Math. Sci., 8 (2010), pp. 1–25

work page 2010

[18] [18]

Caucao, T

S. Caucao, T. Li, and I. Yotov, A multipoint stress-flux mixed finite element method for the Stokes-Biot model, Numer. Math., 152 (2022), pp. 411–473

work page 2022

[19] [19]

Cesmelioglu and P

A. Cesmelioglu and P. Chidyagwai, Numerical analysis of the coupling of free fluid with a poroelastic material, Numer. Methods Partial Differential Equations, 36 (2020), pp. 463– NITSCHE METHODS FOR NS/POROELASTICITY INTERFACE PROBLEMS 23 494

work page 2020

[20] [20]

Cesmelioglu, J

A. Cesmelioglu, J. J. Lee, and S. Rhebergen, Hybridizable discontinuous Galerkin methods for the coupled Stokes-Biot problem , Comput. Math. Appl., 144 (2023), pp. 12–33

work page 2023

[21] [21]

D’Angelo and P

C. D’Angelo and P. Zunino , Numerical approximation with Nitsche’s coupling of transient Stokes’/Darcy’s flow problems applied to hemodynamics , Appl. Numer. Math., 62 (2012), pp. 378–395

work page 2012

[22] [22]

Discacciati, E

M. Discacciati, E. Miglio, and A. Quarteroni , Mathematical and numerical models for coupling surface and groundwater flows , vol. 43, 2002, pp. 57–74

work page 2002

[23] [23]

M. A. Fern ´andez, Incremental displacement-correction schemes for the explicit coupling of a thin structure with an incompressible fluid , C. R. Math. Acad. Sci. Paris, 349 (2011), pp. 473–477

work page 2011

[24] [24]

G. P. Galdi and R. Rannacher, eds., Fundamental trends in fluid-structure interaction, vol. 1 of Contemporary Challenges in Mathematical Fluid Dynamics and Its Applications, World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2010

work page 2010

[25] [25]

Galvis and M

J. Galvis and M. Sarkis, Non-matching mortar discretization analysis for the coupling Stokes- Darcy equations, Electron. Trans. Numer. Anal., 26 (2007), pp. 350–384

work page 2007

[26] [26]

G. N. Gatica, R. Oyarz´ua, and F.-J. Sayas, Analysis of fully-mixed finite element methods for the Stokes-Darcy coupled problem , Math. Comp., 80 (2011), pp. 1911–1948

work page 2011

[27] [27]

Z. Ge, J. Pang, and J. Cao , Multiphysics mixed finite element method with Nitsche’s tech- nique for Stokes-poroelasticity problem, Numer. Methods Partial Differential Equations, 39 (2023), pp. 544–576

work page 2023

[28] [28]

S. Guo, Y. Sun, Y. Wang, X. Yue, and H. Zheng , A Fully Parallelizable Loosely Coupled Scheme for Fluid-Poroelastic Structure Interaction Problems , SIAM J. Sci. Comput., 47 (2025), pp. B951–B975

work page 2025

[29] [29]

Kashiwabara, I

T. Kashiwabara, I. Oikawa, and G. Zhou , Penalty method with Crouzeix-Raviart approx- imation for the Stokes equations under slip boundary condition , ESAIM Math. Model. Numer. Anal., 53 (2019), pp. 869–891

work page 2019

[30] [30]

J. Kuan, S. ˇCani´c, and B. Muha , Fluid-poroviscoelastic structure interaction problem with nonlinear geometric coupling, J. Math. Pures Appl. (9), 188 (2024), pp. 345–445

work page 2024

[31] [31]

Nitsche, ¨Uber ein Variationsprinzip zur L¨ osung von Dirichlet-Problemen bei Verwendung von Teilr¨ aumen, die keinen Randbedingungen unterworfen sind, Abh

J. Nitsche, ¨Uber ein Variationsprinzip zur L¨ osung von Dirichlet-Problemen bei Verwendung von Teilr¨ aumen, die keinen Randbedingungen unterworfen sind, Abh. Math. Sem. Univ. Hamburg, 36 (1971), pp. 9–15

work page 1971

[32] [32]

Quarteroni and A

A. Quarteroni and A. Valli , Numerical approximation of partial differential equations , vol. 23 of Springer Series in Computational Mathematics, Springer-Verlag, Berlin, 1994

work page 1994

[33] [33]

Ruiz-Baier, M

R. Ruiz-Baier, M. Taffetani, H. D. Westermeyer, and I. Yotov, The Biot-Stokes coupling using total pressure: formulation, analysis and application to interfacial flow in the eye , Comput. Methods Appl. Mech. Engrg., 389 (2022), pp. Paper No. 114384, 30

work page 2022

[34] [34]

L. R. Scott and S. Zhang , Finite element interpolation of nonsmooth functions satisfying boundary conditions, Math. Comp., 54 (1990), pp. 483–493

work page 1990

[35] [35]

Seboldt, O

A. Seboldt, O. Oyekole, J. Tamba ˇca, and M. Buka ˇc, Numerical modeling of the fluid- porohyperelastic structure interaction, SIAM J. Sci. Comput., 43 (2021), pp. A2923–A2948

work page 2021

[36] [36]

R. E. Showalter, Poroelastic filtration coupled to Stokes flow , in Control theory of partial differential equations, vol. 242 of Lect. Notes Pure Appl. Math., Chapman & Hall/CRC, Boca Raton, FL, 2005, pp. 229–241

work page 2005

[37] [37]

J. Yu, Y. Sun, F. Shi, and H. Zheng , Nitsche’s type stabilized finite element method for the fully mixed Stokes-Darcy problem with Beavers-Joseph conditions , Appl. Math. Lett., 110 (2020), pp. 106588, 8

work page 2020

[38] [38]

Yu and Y

J. Yu and Y. Zhang , Nitsche’s type stabilization for the fully mixed Navier-Stokes/Darcy problem, J. Appl. Anal. Comput., 11 (2021), pp. 1481–1493. Appendix A. Proof of Theorem 5.3. Proof. We introduce the errors for all variables and split them into approximation and discretization errors: en f := uS,n f − uS,n f,h = (uS,n f − If,huS,n f ) + (If,huS,n f...

work page 2021