# A Posteriori Error Analysis for a Lagrange Multiplier Method for a Stokes/Biot Fluid-Poroelastic Structure Interaction Model

Authors: Houédanou Koffi Wilfrid
Journal: Abstract and Applied Analysis
 (2021)

---

## Abstract

In this work, we develop an a posteriori error analysis of a conforming mixed finite element method for solving the coupled problem arising in the interaction between a free fluid and a fluid in a poroelastic medium on isotropic meshes inℝd, d∈2,3. The approach utilizes a Lagrange multiplier method to impose weakly the interface conditions [Ilona Ambartsumyan et al., Numerische Mathematik, 140 (2): 513-553, 2018]. The a posteriori error estimate is based on a suitable evaluation on the residual of the finite element solution. It is proven that the a posteriori error estimate provided in this paper is both reliable and efficient. The proof of reliability makes use of suitable auxiliary problems, diverse continuous inf-sup conditions satisfied by the bilinear forms involved, Helmholtz decomposition, and local approximation properties of the Clément interpolant. On the other hand, inverse inequalities and the localization technique based on simplexe-bubble and face-bubble functions are the main tools for proving the efficiency of the estimator. Up to minor modifications, our analysis can be extended to other finite element subspaces yielding a stable Galerkin scheme.


---

## Body

## 1. Introduction

In this paper, we develop an a posteriori error analysis for solving the interaction of a free incompressible viscous Newtonian fluid with a fluid within a poroelastic medium. This is a challenging multiphysics problem with applications to predicting and controlling processes arising in groundwater flow in fractured aquifers, oil and gas extraction, arterial flows, and industrial filters. In these applications, it is important to model properly the interaction between the free fluid with the fluid within the porous medium and to take into account the effect of the deformation of the medium. For example, geomechanical effects play an important role in hydraulic fracturing, as well as in modeling phenomena such as subsidence and compaction.We adopt the Stokes equations to model the free fluid and the Biot system [1] for the fluid in the poroelastic media. In the latter, the volumetric deformation of the elastic porous matrix is complemented with the Darcy equation that describes the average velocity of the fluid in the pores. The model features two different kinds of coupling across the interface: Stokes-Darcy coupling [2–10] and fluid-structure interaction (FSI) [11–15].The well-posedness of the mathematical model based on the Stokes-Biot system for the coupling between a fluid and a poroelastic structure is studied in [16]. A numerical study of the problem, using a Navier-Stokes equations for the fluid, is presented in [11, 17], utilizing a variational multiscale approach to stabilize the finite element spaces. The problem is solved using both a monolithic and a partitioned approach, with the latter requiring subiterations between the two problems.Finite element analysis of an arbitrary Lagrangian-Eulerian method for Stokes/parabolic moving interface problem with jump coefficients has been studied in [18]. The authors in [19] study a numerical solution of the coupled system of the time-dependent Stokes and fully dynamic Biot equations. They establish stability of the scheme and derive error estimates for the fully discrete coupled scheme. Numerical errors and convergence rates for smooth problems as well as tests on realistic material parameters have been presented. In [20], Wen and He consider a strongly conservative discretization for the rearranged Stokes-Biot model based on interior penalty discontinuous Galerkin method and mixed finite element method. The existence and uniqueness of solution of the numerical scheme have been presented. Then, the analysis of stability and priori error estimates have been derived. The numerical examples under uniform meshes, which well validate the analysis of convergence and the strong mass conservation, are presented. A staggered finite element procedure for the coupled Stokes-Biot system with fluid entry resistance has been studied by Bergkamp et al. in [21] while Ambartsumyan et al. study in [22] flow and transport in fractured poroelastic media using Stokes flow in the fractures and the Biot model in the porous media. In [23], semidiscrete continuous-in-time approximation has been proposed for the weak coupled mixed formulation. For the discretization of the fluid velocity and pressure, the authors have used the finite elements which include the MINI elements, the Taylor-Hood elements, and the conforming Crouzeix-Raviart elements. For the discretization of the porous medium problem, they choose the spaces that include Raviart-Thomas and Brezzi-Douglas-Marini elements. An a priori error analysis is performed with some numerical tests confirming the convergence rates.A posteriori error estimators are computable quantities, expressed in terms of the discrete solution and of the data that measure the actual discrete errors without the knowledge of the exact solution. They are essential to design adaptive mesh refinement algorithms which aqui-distribute the computational effort and optimize the approximation efficiency. Since the pioneering work of Babuška and Rheinboldt [24–27], adaptive finite element methods based on a posteriori error estimates have been extensively investigated. To our best knowledge, there is no a posteriori error estimation for the Stokes/Biot fluid-poroelastic structure interaction model for finite element methods. Here, we develop such a posteriori error analysis for the semidiscrete conforming finite element methods. We have got a new family of a local indicator error ΘK (see Definition 5, Equation (48)) and global Θ (Equation (52)).The schedule of the paper is as follows. Section2 is devoted to notations and basic results that are used throughout the document. Our main results regarding a posteriori error analysis are stated in Section 3. We prove that our indicator errors are efficiency and reliability, and then, are optimal. The global inf-sup condition is the main tool yielding the reliability. In turn, the local efficiency result is derived using the technique of bubble function introduced by Verfürth [28] and used in similar context by Carstensen and Dolzmann [29]. Finally, this paper is summarized with further works in Section 4.


## 2. Preliminaries and Notations

### 2.1. Stokes-Biot Model Problem

We consider a multiphysics model problem for free fluid’s interaction with a flow in a deformable porous media, where the simulation domainΩ⊂ℝd, d=2,3, is a union of nonoverlapping regions Ωf and Ωp. Here, Ωf is a free fluid region with flow governed by the Stokes equations, and Ωp is a poroelastic material governed by the Biot system. For simplicity of notation, we assume that each region is connected. The extension to nonconnected regions is straightforward. Let Γfp=∂Ωf∩∂Ωp (see Figure 1).Figure 1
Global domainΩ consisting of the fluid region Ωf and the poroelastic media region Ωp separated by the interface Γfp.Letu∗,p∗ be the velocity-pressure pair in Ω∗, ∗=f,p, and let ηp be the displacement in Ωp. Let μ>0 be the fluid viscosity, let f∗ be the body force terms, and let q∗ be external source or sink terms. Let Duf and σfuf,pf denote, respectively, the deformation rate tensor and the stress tensor:
(1)Duf=12∇uf+∇ufT,σfuf,pf=−pfI+2μDuf.In the free fluid regionΩf, uf,pf satisfy the Stokes equations:
(2)−∇·σfuf,pf=ffinΩf×0,T,(3)∇·uf=qfinΩf×0,T,where T>0 is the final time. Let σeηp and σpηp,pp be the elastic and poroelastic stress tensors, respectively:
(4)σeηp=λp∇·ηpI+2μpDηp,σpηp,pp=σeηp−αppI,where 0<λmin≤λpx≤λmax and 0<μmin≤μpx≤μmax are the Lamé parameters and 0<α≤1 is the Biot-Willis constant. The poroelasticity region Ωp is governed by the quasistatic Biot system [23]:
(5)−∇·σpηp,pp=fp,μK−1up+∇pp=0inΩp×0,T,(6)∂∂ts0pp+α∇·ηp+∇·up=qpinΩp×0,T,where s0≥0 is a storage coefficient and K the symmetric and uniformly positive definite rock permeability tensor, satisfying, for some constants 0<kmin≤kmax,
(7)∀ξ∈ℝd,kminξTξ≤ξTKxξ≤kmaxξTξ,∀x∈Ωp.Following [1], the interface conditions on the fluid-poroelasticity interface Γfp are mass conservation, balance of stresses, and the Beavers-Joseph-Saffman (BJS) condition [30] modeling slip with friction:
(8)uf·nf+∂ηp∂t+up·np=0onΓfp×0,T,(9)−σfnf·nf=pp,σfnf+σpnp=0onΓfp×0,T,(10)−σfnf·τf,j=μαBJFKj−1uf−∂ηp∂t·τf,jonΓfp×0,T,where nf and np are the outward unit normal vectors to ∂Ωf and ∂Ωp, respectively, τf,j, 1≤j≤d−1, is an orthogonal system of unit tangent vectors on Γfp, Kj=Kτf,j·τf,j, and αBJS≥0 is an experimentally determined friction coefficient. We note that continuity of flux constraints the normal velocity of the solid skeleton, while the BJS condition accounts for its tangential velocity.The above system of equations needs to be complemented by a set of boundary and initial conditions. LetΓf=∂Ωf∩∂Ω and Γp=∂Ωp∩∂Ω. Let Γp=ΓpD∪ΓpN. We assume for simplicity homogeneous boundary conditions:
(11)uf=0onΓf×0,T,ηp=0onΓp×0,T,pp=0onΓpD×0,T,up·np=0onΓpN×0,T.To ovoid the issue with restricting the mean value of the pressure, we assume that∣ΓpD∣>0. We also assume that ΓpD is not adjacent to the interface Γfp, i.e., distΓpD,Γfp≥s>0. Nonhomogeneous displacement and velocity conditions can be handled in a standard way by adding suitable extensions of the boundary data. The pressure boundary condition is natural in the mixed Darcy formulation, so nonhomogeneous pressure data would lead to an additional boundary term. We further sat the initial conditions:
(12)ppx,0=pp,0x,ηpx,0=ηp,0xinΩp.Equations (2)–(10) consist of the model of the coupled Stokes and Biot flows problem that we will study below.


### 2.2. Weak Formulation

In this part, we first introduce some Sobolev spaces [31] and norms. If W is a bounded domain of ℝd and m is a nonnegative integer, the Sobolev space HmW=Wm,2W is defined in the usual way with the usual norm ∥·∥m,W and seminorm .m,W. In particular, H0W=L2W, and we write ∥·∥W for ∥·∥0,W. Similarly, we denote by ·,·W the L2WL2Wd or L2Wd×d inner product. For shortness if W is equal to Ω, we will drop the index Ω, while for any m≥0, ∥·∥m,å=∥·∥m,Ωå, .m,å=.m,Ωå and .,.å=·,·Ωå, for å=f,p. The space H0mΩ denotes the closure of C0∞Ω in HmΩ. Let HmΩd be the space of vector valued functions v=v1,⋯,vd with components vi in HmΩ. The norm and the seminorm on HmΩd are given by
(13)∥v∥m,Ω≔∑i=0d∥vi∥m,Ω21/2,vm,Ω≔∑i=0dvim,Ω21/2.For a connected open subset of the boundaryE⊂∂Ωf∪∂Ωp, we write .,.E for the L2E inner product (or duality pairing), that is, for scalar valued functions λ, σ one defines:
(14)λ,σE≔∫Eλσds.In the following, we derive a Lagrange multiplier type weak formulation of the system, which will be the basis for our finite element approximation. Let(15)Vf=vf∈H1Ωfd:vf=0onΓf,Wf=L2Ωf,Vp=vp∈H;Ωp:vp·np=0onΓpN,Wp=L2Ωp,Xp=ξp∈H1Ωpd:ξp=0onΓp,where H;Ωp is the space of L2Ωpd vectors with divergence in L2Ωp with a norm
(16)vH;Ωp2≔vΩp2+∇·vΩp2.We define the global velocity and pressure spaces as(17)V=v=vf,vp∈Vf×Vp,W=w=wf,wp∈Wf×Wp,with norms
(18)∥v∥V2≔∥vf∥Ωf2+∥vp∥H;Ωp2,∥w∥W2=∥wf∥Ωf2+∥wp∥Ωp2.The weak formulation is obtained by mutiplying the equations in each region by suitable test functions, integrating by parts the second-order terms in space, and utilizing the interface and boundary conditions.Let define(19)afuf,vf≔2μDuf,DvfΩf,apdup,vp≔μK−1up,vpΩp,apeηp,ξp≔2μpDηp,DξpΩp+λp∇·ηp,∇·ξpΩp,be the bilinear forms related to Stokes, Darcy, and the elasticity operator, respectively. Let
(20)b∗v,w≔−∇·v,wΩ∗.Integration by parts in (2) and the two equations in (5) lead to the interface term
(21)IΓfp=−σfnf,vfΓfp−σpnp,ξpΓfp+pp,vp·npΓfp.Using the first condition for balance of normal stress in (9), we set
(22)λ=−σfnf·nf=pponΓfp,which will be used as a Lagrange multiplier to impose the mass conservation interface condition (8). Utilizing the BJS condition (10) and the second condition for balance of stresses in (9), we obtain
(23)IΓfp=aBJSuf,∂tηp;vf,ξp+bΓvf,vp,ξp,λ,where
(24)aBJSuf,ηp;vf,ξp=∑j=1d−1μαBJSKj−1uf−ηp·τf,j,vf−ξp·τf,jΓfp,bΓvf,vp,ξp;μ=vf·nf+ξp+vp·np,μΓfp.For the well-posedness ofbΓ, we require that λ∈∧=Vp·np∣Γfp′. According to the normal trace theorem, since vp∈Vp⊂H;Ωp, then vp·np∈H−1/2∂Ωp. Furthermore, since vp·np=0 on ΓpN and disΓpD,Γfp≥s>0, then vp·np∈H−1/2Γfp (see, e.g., [32]). Therefore, we take ∧=H1/2Γfp.The Lagrange multiplier variational formulation is fort∈0,T, find uft∈Vf, pft∈Wf, upt∈Vp, ppt∈Wp, ηpt∈Xp, and λt∈∧ such that pp0=pp,0, ηp0=ηp,0, and for all vf∈Vf, wf∈Wf, vp∈Vp, wp∈Wp, ξp∈Xp, and μ∈∧,
(25)afuf,vf+apdup,vp+apeηp,ξp+aBJSuf,∂tηp;vf,ξp,(26)+bfvf,pf+bpvp,pp+αbpξp,pp+bΓvf,vp,ξp;λ=ff,vfΩf+fp,ξpΩps0∂tpp,wpΩp−αbp∂tηp,wp−bpup,wp−bfuf,wf=qf,wfΩf+qp,wpΩp,(27)bΓuf,up,∂tηp;μ=0.where we used the notation ∂t=∂/∂t.The assumptions on the fluid viscosityμ and the material coefficients K, λp, and μp imply that the bilinear forms af·,·, apd·,·, and ape·,· are coercive and continuous in the appropriate norms. In particular, there exist positive constants cf, cp, ce, Cf, Cp, Ce such that:
(28)cf∥vf∥H1Ωf2≤afvf,vf,∀vf∈Vf,(29)afvf,qf≤Cf∥vf∥H1Ωf∥qf∥H1Ωf,∀vf,qf∈Vf,(30)cp∥vp∥L2Ωp2≤apdvp,vp,∀vp∈Vp,(31)apdvp,qp≤Cp∥vp∥L2Ωp∥qp∥L2Ωp,∀vp,qp∈Vp,(32)ce∥ξp∥H1Ωp2≤apeξp,ξp∀ξp∈Xp,(33)apeξp,ζp≤Ce∥ξp∥H1Ωp∥ζp∥H1Ωp∀ξp,ζp∈Xp,where (28)–(29) and (32)–(33) hold true thanks to Poincaré inequality and (32)–(33) also relies on Korn’s inequality.In summary, from ([23], Corollary 3.1, Page 7), the following result holds:Theorem 1.
There exists a unique solutionuf,pf,up,pp,ηp,λ∈L∞0,T;Vf×L∞0,T;Wf×L∞0,T;Vp×W1,∞0,T;Wp×W1,∞0,T;Xp×L∞0,T;Λ to the problem (25)–(27).


### 2.3. Finite Element Discretization

LetThf and Thp be shape-regular and quasiuniform partition of Ωf and Ωp, respectively, both consisting of affine elements with maximal element diameter h. The two partitions may be nonmatching at the interface Γfp. For the discretization of the fluid velocity and pressure, we choose finite element spaces Vf,h⊂Vf and Wf,h⊂Wf, which are assumed to be inf-sup stable. Examples of such spaces include the MINI elements, the Taylor-Hood elements, and the conforming Crouzeix-Raviart elements. For the discretization of the porous medium problem, we choose Vp,h⊂Vp and Wp,h⊂Wp to be any of well-known inf-sup stable mixed finite element spaces, such as the Raviart-Thomas or the Brezzi-Douglas-Marini spaces. The global spaces are:
(34)Vh≔vh=vf,h,vp,h∈Vf,h×Vp,h,Wh=wh=wf,h,wp,h∈Wf,h×Wp,h.We employ a conforming Lagrangian finite element spaceXp,h⊂Xp to approximate the structure displacement. Note that the finite element spaces Vf,h, Vp,h, and Xp,h satisfy the prescribed homogeneous boundary conditions on the external boundaries. For the discrete Lagrange multiplier space, we take
(35)Λh=Vp,h·np∣Γfp.The semidiscrete continuous-in-time problem reads givenpp,h0 and ηp,h0, for t∈0,T, find uf,ht∈Vf,h, pf,ht∈Wf,h, up,ht∈Vp,h, pp,ht∈Wp,h, ηp,ht∈Xp,h, and λht∈Λh such that for all vf,h∈Vf,h, wf,h∈Wf,h, vp,h∈Vp,h, wp,h∈Wp,h, ξp,h∈Xp,h, and μh∈Λh,
(36)afuf,h,vf,h+apdup,h,vp,h+apeηp,h,ξp,h+aBJSuf,h,∂tηp,h;vf,h,ξp,h+bfvf,h,pf,h+bpvp,h,pp,h+αbpξp,h,pp,h+bΓvf,h,vp,h,ξp,h;λh=ff,h,vf,hΩf+fp,h,ξp,hΩp,(37)s0∂tpp,h,wp,hΩp−αbp∂tηp,h,wp,h−bpup,h,wp,h−bfuf,h,wf,h=qf,wf,hΩf+qp,wp,hΩp,(38)bΓuf,h,up,h,∂tηp,h;μh=0.We will takepp,h0 and ηp,h0 to be suitable projections of the initial data pp,0 and ηp,0.We introduce the errors for all variables:(39)ef≔uf−uf,h,ep≔up−up,h,es≔ηp−ηp,h,efp≔pf−pf,h,epp≔pp−pp,h,eλ≔λ−λh.The following results hold cf. [23]:Theorem 2 (A priori error estimation).
There exists a unique solutionuf,h,pf,h,up,h,pp,h,ηp,h,λh in L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh of the discrete problem (36)–(38), and if the solution uf,pf,up,pp,ηp,λ∈Vf×Wf×Vp×Wp×Xp×Λ of the continuous problem (25)–(27) is smooth enough, then we have:
(40)∥ef∥L20,T;H1Ωf+∥ep∥L20,T;L2Ωp+∥es∥L∞0,T;H1Ωp+∥ep∥L∞0,T;L2Ωp+∥ef∥L20,T;L2Ωf+∥eλ∥L20,T;Λh≤Cuf,pf,up,pp,ηp,λhr,r>0.ForWh=vf,h,wf,h,vp,h,wp,h,ξp,h,λh∈L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh, we can subtract (36)–(38) to (25)–(27) to obtain the Galerkin orthogonality relation for all t∈0,T:
(41)afeft,vf,ht+apdest,ξp,ht+aBJSeft,∂test;vf,ht,ξp,ht+bfvf,ht,eft+bpvp,ht,ept+αbpξpt,ept+bΓvft,vpt,ξpt,eλt+s0∂tept,wp,ht−αbp∂ttest,wp,ht−bpept,wp,ht−bfeft,wf,ht=0.


## 2.1. Stokes-Biot Model Problem

We consider a multiphysics model problem for free fluid’s interaction with a flow in a deformable porous media, where the simulation domainΩ⊂ℝd, d=2,3, is a union of nonoverlapping regions Ωf and Ωp. Here, Ωf is a free fluid region with flow governed by the Stokes equations, and Ωp is a poroelastic material governed by the Biot system. For simplicity of notation, we assume that each region is connected. The extension to nonconnected regions is straightforward. Let Γfp=∂Ωf∩∂Ωp (see Figure 1).Figure 1
Global domainΩ consisting of the fluid region Ωf and the poroelastic media region Ωp separated by the interface Γfp.Letu∗,p∗ be the velocity-pressure pair in Ω∗, ∗=f,p, and let ηp be the displacement in Ωp. Let μ>0 be the fluid viscosity, let f∗ be the body force terms, and let q∗ be external source or sink terms. Let Duf and σfuf,pf denote, respectively, the deformation rate tensor and the stress tensor:
(1)Duf=12∇uf+∇ufT,σfuf,pf=−pfI+2μDuf.In the free fluid regionΩf, uf,pf satisfy the Stokes equations:
(2)−∇·σfuf,pf=ffinΩf×0,T,(3)∇·uf=qfinΩf×0,T,where T>0 is the final time. Let σeηp and σpηp,pp be the elastic and poroelastic stress tensors, respectively:
(4)σeηp=λp∇·ηpI+2μpDηp,σpηp,pp=σeηp−αppI,where 0<λmin≤λpx≤λmax and 0<μmin≤μpx≤μmax are the Lamé parameters and 0<α≤1 is the Biot-Willis constant. The poroelasticity region Ωp is governed by the quasistatic Biot system [23]:
(5)−∇·σpηp,pp=fp,μK−1up+∇pp=0inΩp×0,T,(6)∂∂ts0pp+α∇·ηp+∇·up=qpinΩp×0,T,where s0≥0 is a storage coefficient and K the symmetric and uniformly positive definite rock permeability tensor, satisfying, for some constants 0<kmin≤kmax,
(7)∀ξ∈ℝd,kminξTξ≤ξTKxξ≤kmaxξTξ,∀x∈Ωp.Following [1], the interface conditions on the fluid-poroelasticity interface Γfp are mass conservation, balance of stresses, and the Beavers-Joseph-Saffman (BJS) condition [30] modeling slip with friction:
(8)uf·nf+∂ηp∂t+up·np=0onΓfp×0,T,(9)−σfnf·nf=pp,σfnf+σpnp=0onΓfp×0,T,(10)−σfnf·τf,j=μαBJFKj−1uf−∂ηp∂t·τf,jonΓfp×0,T,where nf and np are the outward unit normal vectors to ∂Ωf and ∂Ωp, respectively, τf,j, 1≤j≤d−1, is an orthogonal system of unit tangent vectors on Γfp, Kj=Kτf,j·τf,j, and αBJS≥0 is an experimentally determined friction coefficient. We note that continuity of flux constraints the normal velocity of the solid skeleton, while the BJS condition accounts for its tangential velocity.The above system of equations needs to be complemented by a set of boundary and initial conditions. LetΓf=∂Ωf∩∂Ω and Γp=∂Ωp∩∂Ω. Let Γp=ΓpD∪ΓpN. We assume for simplicity homogeneous boundary conditions:
(11)uf=0onΓf×0,T,ηp=0onΓp×0,T,pp=0onΓpD×0,T,up·np=0onΓpN×0,T.To ovoid the issue with restricting the mean value of the pressure, we assume that∣ΓpD∣>0. We also assume that ΓpD is not adjacent to the interface Γfp, i.e., distΓpD,Γfp≥s>0. Nonhomogeneous displacement and velocity conditions can be handled in a standard way by adding suitable extensions of the boundary data. The pressure boundary condition is natural in the mixed Darcy formulation, so nonhomogeneous pressure data would lead to an additional boundary term. We further sat the initial conditions:
(12)ppx,0=pp,0x,ηpx,0=ηp,0xinΩp.Equations (2)–(10) consist of the model of the coupled Stokes and Biot flows problem that we will study below.


## 2.2. Weak Formulation

In this part, we first introduce some Sobolev spaces [31] and norms. If W is a bounded domain of ℝd and m is a nonnegative integer, the Sobolev space HmW=Wm,2W is defined in the usual way with the usual norm ∥·∥m,W and seminorm .m,W. In particular, H0W=L2W, and we write ∥·∥W for ∥·∥0,W. Similarly, we denote by ·,·W the L2WL2Wd or L2Wd×d inner product. For shortness if W is equal to Ω, we will drop the index Ω, while for any m≥0, ∥·∥m,å=∥·∥m,Ωå, .m,å=.m,Ωå and .,.å=·,·Ωå, for å=f,p. The space H0mΩ denotes the closure of C0∞Ω in HmΩ. Let HmΩd be the space of vector valued functions v=v1,⋯,vd with components vi in HmΩ. The norm and the seminorm on HmΩd are given by
(13)∥v∥m,Ω≔∑i=0d∥vi∥m,Ω21/2,vm,Ω≔∑i=0dvim,Ω21/2.For a connected open subset of the boundaryE⊂∂Ωf∪∂Ωp, we write .,.E for the L2E inner product (or duality pairing), that is, for scalar valued functions λ, σ one defines:
(14)λ,σE≔∫Eλσds.In the following, we derive a Lagrange multiplier type weak formulation of the system, which will be the basis for our finite element approximation. Let(15)Vf=vf∈H1Ωfd:vf=0onΓf,Wf=L2Ωf,Vp=vp∈H;Ωp:vp·np=0onΓpN,Wp=L2Ωp,Xp=ξp∈H1Ωpd:ξp=0onΓp,where H;Ωp is the space of L2Ωpd vectors with divergence in L2Ωp with a norm
(16)vH;Ωp2≔vΩp2+∇·vΩp2.We define the global velocity and pressure spaces as(17)V=v=vf,vp∈Vf×Vp,W=w=wf,wp∈Wf×Wp,with norms
(18)∥v∥V2≔∥vf∥Ωf2+∥vp∥H;Ωp2,∥w∥W2=∥wf∥Ωf2+∥wp∥Ωp2.The weak formulation is obtained by mutiplying the equations in each region by suitable test functions, integrating by parts the second-order terms in space, and utilizing the interface and boundary conditions.Let define(19)afuf,vf≔2μDuf,DvfΩf,apdup,vp≔μK−1up,vpΩp,apeηp,ξp≔2μpDηp,DξpΩp+λp∇·ηp,∇·ξpΩp,be the bilinear forms related to Stokes, Darcy, and the elasticity operator, respectively. Let
(20)b∗v,w≔−∇·v,wΩ∗.Integration by parts in (2) and the two equations in (5) lead to the interface term
(21)IΓfp=−σfnf,vfΓfp−σpnp,ξpΓfp+pp,vp·npΓfp.Using the first condition for balance of normal stress in (9), we set
(22)λ=−σfnf·nf=pponΓfp,which will be used as a Lagrange multiplier to impose the mass conservation interface condition (8). Utilizing the BJS condition (10) and the second condition for balance of stresses in (9), we obtain
(23)IΓfp=aBJSuf,∂tηp;vf,ξp+bΓvf,vp,ξp,λ,where
(24)aBJSuf,ηp;vf,ξp=∑j=1d−1μαBJSKj−1uf−ηp·τf,j,vf−ξp·τf,jΓfp,bΓvf,vp,ξp;μ=vf·nf+ξp+vp·np,μΓfp.For the well-posedness ofbΓ, we require that λ∈∧=Vp·np∣Γfp′. According to the normal trace theorem, since vp∈Vp⊂H;Ωp, then vp·np∈H−1/2∂Ωp. Furthermore, since vp·np=0 on ΓpN and disΓpD,Γfp≥s>0, then vp·np∈H−1/2Γfp (see, e.g., [32]). Therefore, we take ∧=H1/2Γfp.The Lagrange multiplier variational formulation is fort∈0,T, find uft∈Vf, pft∈Wf, upt∈Vp, ppt∈Wp, ηpt∈Xp, and λt∈∧ such that pp0=pp,0, ηp0=ηp,0, and for all vf∈Vf, wf∈Wf, vp∈Vp, wp∈Wp, ξp∈Xp, and μ∈∧,
(25)afuf,vf+apdup,vp+apeηp,ξp+aBJSuf,∂tηp;vf,ξp,(26)+bfvf,pf+bpvp,pp+αbpξp,pp+bΓvf,vp,ξp;λ=ff,vfΩf+fp,ξpΩps0∂tpp,wpΩp−αbp∂tηp,wp−bpup,wp−bfuf,wf=qf,wfΩf+qp,wpΩp,(27)bΓuf,up,∂tηp;μ=0.where we used the notation ∂t=∂/∂t.The assumptions on the fluid viscosityμ and the material coefficients K, λp, and μp imply that the bilinear forms af·,·, apd·,·, and ape·,· are coercive and continuous in the appropriate norms. In particular, there exist positive constants cf, cp, ce, Cf, Cp, Ce such that:
(28)cf∥vf∥H1Ωf2≤afvf,vf,∀vf∈Vf,(29)afvf,qf≤Cf∥vf∥H1Ωf∥qf∥H1Ωf,∀vf,qf∈Vf,(30)cp∥vp∥L2Ωp2≤apdvp,vp,∀vp∈Vp,(31)apdvp,qp≤Cp∥vp∥L2Ωp∥qp∥L2Ωp,∀vp,qp∈Vp,(32)ce∥ξp∥H1Ωp2≤apeξp,ξp∀ξp∈Xp,(33)apeξp,ζp≤Ce∥ξp∥H1Ωp∥ζp∥H1Ωp∀ξp,ζp∈Xp,where (28)–(29) and (32)–(33) hold true thanks to Poincaré inequality and (32)–(33) also relies on Korn’s inequality.In summary, from ([23], Corollary 3.1, Page 7), the following result holds:Theorem 1.
There exists a unique solutionuf,pf,up,pp,ηp,λ∈L∞0,T;Vf×L∞0,T;Wf×L∞0,T;Vp×W1,∞0,T;Wp×W1,∞0,T;Xp×L∞0,T;Λ to the problem (25)–(27).


## 2.3. Finite Element Discretization

LetThf and Thp be shape-regular and quasiuniform partition of Ωf and Ωp, respectively, both consisting of affine elements with maximal element diameter h. The two partitions may be nonmatching at the interface Γfp. For the discretization of the fluid velocity and pressure, we choose finite element spaces Vf,h⊂Vf and Wf,h⊂Wf, which are assumed to be inf-sup stable. Examples of such spaces include the MINI elements, the Taylor-Hood elements, and the conforming Crouzeix-Raviart elements. For the discretization of the porous medium problem, we choose Vp,h⊂Vp and Wp,h⊂Wp to be any of well-known inf-sup stable mixed finite element spaces, such as the Raviart-Thomas or the Brezzi-Douglas-Marini spaces. The global spaces are:
(34)Vh≔vh=vf,h,vp,h∈Vf,h×Vp,h,Wh=wh=wf,h,wp,h∈Wf,h×Wp,h.We employ a conforming Lagrangian finite element spaceXp,h⊂Xp to approximate the structure displacement. Note that the finite element spaces Vf,h, Vp,h, and Xp,h satisfy the prescribed homogeneous boundary conditions on the external boundaries. For the discrete Lagrange multiplier space, we take
(35)Λh=Vp,h·np∣Γfp.The semidiscrete continuous-in-time problem reads givenpp,h0 and ηp,h0, for t∈0,T, find uf,ht∈Vf,h, pf,ht∈Wf,h, up,ht∈Vp,h, pp,ht∈Wp,h, ηp,ht∈Xp,h, and λht∈Λh such that for all vf,h∈Vf,h, wf,h∈Wf,h, vp,h∈Vp,h, wp,h∈Wp,h, ξp,h∈Xp,h, and μh∈Λh,
(36)afuf,h,vf,h+apdup,h,vp,h+apeηp,h,ξp,h+aBJSuf,h,∂tηp,h;vf,h,ξp,h+bfvf,h,pf,h+bpvp,h,pp,h+αbpξp,h,pp,h+bΓvf,h,vp,h,ξp,h;λh=ff,h,vf,hΩf+fp,h,ξp,hΩp,(37)s0∂tpp,h,wp,hΩp−αbp∂tηp,h,wp,h−bpup,h,wp,h−bfuf,h,wf,h=qf,wf,hΩf+qp,wp,hΩp,(38)bΓuf,h,up,h,∂tηp,h;μh=0.We will takepp,h0 and ηp,h0 to be suitable projections of the initial data pp,0 and ηp,0.We introduce the errors for all variables:(39)ef≔uf−uf,h,ep≔up−up,h,es≔ηp−ηp,h,efp≔pf−pf,h,epp≔pp−pp,h,eλ≔λ−λh.The following results hold cf. [23]:Theorem 2 (A priori error estimation).
There exists a unique solutionuf,h,pf,h,up,h,pp,h,ηp,h,λh in L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh of the discrete problem (36)–(38), and if the solution uf,pf,up,pp,ηp,λ∈Vf×Wf×Vp×Wp×Xp×Λ of the continuous problem (25)–(27) is smooth enough, then we have:
(40)∥ef∥L20,T;H1Ωf+∥ep∥L20,T;L2Ωp+∥es∥L∞0,T;H1Ωp+∥ep∥L∞0,T;L2Ωp+∥ef∥L20,T;L2Ωf+∥eλ∥L20,T;Λh≤Cuf,pf,up,pp,ηp,λhr,r>0.ForWh=vf,h,wf,h,vp,h,wp,h,ξp,h,λh∈L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh, we can subtract (36)–(38) to (25)–(27) to obtain the Galerkin orthogonality relation for all t∈0,T:
(41)afeft,vf,ht+apdest,ξp,ht+aBJSeft,∂test;vf,ht,ξp,ht+bfvf,ht,eft+bpvp,ht,ept+αbpξpt,ept+bΓvft,vpt,ξpt,eλt+s0∂tept,wp,ht−αbp∂ttest,wp,ht−bpept,wp,ht−bfeft,wf,ht=0.


## 3. A Posteriori Error Analysis

In order to solve the Stokes-Biot model problem by efficient adaptive finite element methods, reliable and efficient a posteriori error analysis is important to provide appropriated indicators. In this section, we first define the local and global indicators (Section3.1), and then, the lower and upper error bounds are derived (Sections 3.4 and 3.5).


### 3.1. Residual Error Estimators

The general philosophy of residual error estimators is to estimate an appropriate norm of the correct residual by terms that can be evaluated easier and that involve the data at hand. To this end, define the exact element residuals:Definition 3 (Exact Element Residuals).
Lett∈0,T and Wht=vf,ht,pf,ht,vp,ht,wp,ht,ξp,ht,λht∈Vf,h×Wf,h×Vp,h×Wp,h×Xp,h×Λh be an arbitrary finite element function. The exact element residuals over a triangle or tetrahedra K∈h and over E∈hΓfp are defined for all t∈0,T by:
(42)Rf,KWht=ff+∇·σfvf,ht,wf,ht,Rp,K,1Wht=fp+∇·σpξp,ht,wp,ht,Rp,K,2Wht=μK−1vp,ht+∇wp,ht,Rf,KWht=qf−∇·vf,ht,Rp,KWht=qp−∂ts0wp,ht+α∇·ξp,ht+∇·vp,ht,RE,pf,1Wht=vf,ht·nf,E+∂tξpt+vp,ht·np,E,RE,pf,2Wht=wp,h+σfvf,ht,wf,htnf,E·nf,E,RE,pf,3Wht=σfvf,ht,wf,htnf,E+σpξp,ht,wp,htnp,E,RE,pf,4Whtj=σfvf,h,wf,hτf,E,j+μαBJSKj−1vf,h−∂tξp,h·τf,E,j.As it is common, these exact residuals are replaced by some finite-dimensional approximation called approximate element residualr∗,K, r∗,K, ∗∈f,p, RE,pf,l, l∈1,2,3,4. This approximation is here achieved by projecting ff and qf on the space of piecewise constant functions in Ωf and piecewise ℙ1 functions in Ωp, more precisely for all K∈hf we take
(43)ff,K=pKff,qf,K=PKqf,withPK:L1K→ℝsuchthatpKϕ=1∣K∣∫Kϕxdx.While for allK∈hp, we take fp,K and qp as the unique element of ℙ1Kd, respectively ℙ1K such that:
(44)∫Kfp,Kx·qxdx=∫Kfpxqxdx,q∈ℙ1Kd,respectively,
(45)∫Kqp,Kxqxdx=∫Kqpxqxdx,∀q∈ℙ1K.Thereby, we define the approximate element residuals.Definition 4 [Approximate Element Residuals].
Lett∈0,T and Wht=vf,ht,pf,ht,vp,ht,wp,ht,ξp,ht,λht∈Vf,h×Wf,h×Vp,h×Wp,h×Xp,h×Λh be an arbitrary finite element function. Then, the approximate element residuals are defined for all t∈0,T by:
(46)rf,KWht=ff,K+∇·σfvf,ht,wf,ht,rp,K,1Wht=fp,K+∇·σpξp,ht,wp,ht,rp,K,2Wht=μK−1vp,ht+∇wp,ht,rf,KWht=qf,K−∇·vf,ht,rp,KWht=qp,K−∂ts0wp,ht+α∇·ξp,ht+∇·vp,ht.Next, introduce the gradient jump in normal direction by(47)JE,nE,fUh≔2μDuf,h−pf,hI·nEEifE∈hΩf,JE,nE,pUhJE,nE,pUh≔2μDηp,h−pp,hI·nEEifE∈hΩp.where I is the identity matrix of ℝd×d.Definition 5 [Residual Error Estimators].
LetUh=uf,h,pf,h,up,h,pp,h,ηp,h,λh be the finite element solution of the problem (36)–(38) in L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh. Then, the residual error estimator is locally defined by
(48)ΘKUh≔ΘK,f2Uh+ΘK,p2Uh+ΘK,pf2Uh1/2,where
(49)ΘK,f2Uh≔hK2∥rf,KUh∥L∞0,T;L2K2+∥rf,KUh∥L∞0,T;L2K2+∑E∈hΩfhE∥JE,nE,fUh∥L∞0,T;L2E2,(50)ΘK,p2Uh≔hK2∥rp,K,1Uh∥L∞0,T;L2K2+∥rp,K,2Uh∥L∞0,T;L2K2+hK2∥rp,K,2Uh∥L∞0,T;L2K2+∥rp,KUh∥L∞0,T;L2K2+∑E∈hΩphE∥JE,nE,pUh∥L∞0,T;L2E2,(51)ΘK,pf2Uh≔∑E∈hΓpfhE∥RE,pf,1Uh∥L∞0,T;L2E2+∑E∈hΓpfhE∥RE,pf,2Uh∥L∞0,T;L2E2+∑E∈hΓpfhE∥RE,pf,3Uh∥L∞0,T;L2E2+∑E∈hΓpfhE∑j=1d−1∥RE,pf,4j∥L∞0,T,L2E2.The global residual error estimator is given by(52)ΘUh≔∑K∈hΘKUh21/2.Furthermore, denote the local approximation terms by(53)ζK2≔ζK,f2+ζK,p2,where
(54)ζK,f2≔hK2∥Rf,KUh−rf,KUh∥L∞0,T;L2K2+hK2∥Rf,KUh−rf,KUh∥L∞0,T,L2K2ifK∈hf,ζK,p2≔hK2∥Rp,K,1Uh−rp,K,1Uh∥L∞0,T;L2K2+hK2∥Rp,KUh−rp,KUh∥L∞0,T,L2K2ifK∈hp.The global approximation term is defined by(55)ζ≔∑K∈hζK21/2.Remark 6.
The residual character of each term on the right-hand sides of (48)–(51) is quite clear since if Uh would be the exact solution of (2)–(10), then they would vanish.


### 3.2. Analytical Tools

#### 3.2.1. Inverse Inequalities

In order to derive the lower error bounds, we proceed similarly as in [29, 33] (see also [34]), by applying inverse inequalities, and the localization technique based on simplex-bubble and face-bubble functions. To this end, we recall some notation and introduce further preliminary results. Given K∈Th, and E∈K, we let bK and bE be the usual simplexe-bubble and face-bubble functions, respectively (see (1.5) and (1.6) in [28]). In particular, bK satisfies bK∈ℙ3K, suppbK⊆K, bK=0on∂K, and 0≤bK≤1onK. Similarly, bE∈ℙ2K, suppbE⊆ωE≔K′∈Th:E∈K′, bE=0on∂K\E and 0≤bE≤1inωE. We also recall from [35] that, given k∈ℕ, there exists an extension operator L:CE→CK that satisfies Lp∈ℙkK and Lp∣E=p,∀p∈ℙkE. A corresponding vectorial version of L, that is, the componentwise application of L, is denoted by L. Additional properties of bK, bE, and L are collected in the following lemma (see [35]):Lemma 7.
Givenk∈ℕ∗, there exist positive constants depending only on k and shape regularity of the triangulations (minimum angle condition), such that for each simplexe K and E∈K there hold
(56)∥ϕ∥K≲∥ϕbK1/2∥K≲∥ϕ∥K,∀ϕ∈ℙkK,ϕbK1,K≲hK−1∥ϕ∥K,∀ϕ∈ℙkK,∥ψ∥E≲∥bE1/2ψ∥E≲∥ψ∥E,∀ψ∈ℙkE,∥Lψ∥K+hELψ1,K≲hE1/2∥ψ∥E∀ψ∈ℙkE.Lemma 8 (Continuous trace inequality).
There exists a positive constantβ1>0 depending only on σ0 such that
(57)∥v∥∂K2≤β1∥v∥K∥v∥1,K,∀K∈Th,∀v∈H1Kd.


#### 3.2.2. Clément Interpolation Operator

In order to derive the upper error bounds, we introduce the Clément interpolation operatorICl0:H01Ω→PcbTh that approximates optimally nonsmooth functions by continuous piecewise linear functions:
(58)PcbTh≔v∈C0Ω¯:v∣K∈ℙ1K,∀K∈Thandv=0on∂Ω.In addition, we will make use of a vector valued version ofICl0, that is,ICl0:H01Ωd→PcbThd, which is defined componentwise by ICl0. The following lemma establishes the local approximation properties of ICl0 (and hence of ICl0), for a proof see ([36], Section 3).Lemma 9.
There exist constantsC1,C2>0, independent of h, such that for all v∈H01Ω there hold
(59)∥v−ICl0v∥K≤C1hK∥v∥1,ΔK,∀K∈Th,∥v−ICl0v∥E≤C2hE1/2∥v∥1,ΔE,∀E∈h,
whereΔK≔∪K′∈Th:K′∩K≠∅ and ΔE≔∪K′∈Th:K′∩E≠∅.


#### 3.2.3. Helmholtz Decomposition

Lemma 10 ([7]).
There existsCp>0 such that every vp∈H;Ωp can be decomposed as vp=w+β, where w∈H1Ωpd, β∈H1Ωp, ∫Ωpβxdx=0 and
(60)∥w∥1,Ωp+∥β∥1,Ωp≤Cp∥vp∥H;Ωp.


### 3.3. Error Equation

We setH=L∞0,T;Vf×L∞0,T;Wf×L∞0,T;Vp×W1,∞0,T;Wp×W1,∞0,T;Xp×L∞0,T;Λ and Hh=L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh.LetU=uf,pf,up,pp,ηp,λ∈H and W=vf,wf,vp,wp,ξp,μ∈H. For t∈0,T, we define the operator A by
(61)AUt,Wt≔afuft,vft+apdupt,vpt+bpvpt,ppt+αbpξpt,ppt+bΓvft,vpt,ξpt,λt+s0∂tppt,wptΩp−αbp∂tηpt,wpt−bpupt,wpt−bfuft,wft+bΓuft,upt,∂tηpt;μ,FWt≔ff,vftΩf+fp,ξpt+qf,wftΩf+qp,wptΩp.Then, the continuous problem (25)–(27) is equivalent to find U∈H such that for all t∈0,T, we have:
(62)AUt,Wt=FWt,∀W∈H.We define the discrete version by the same way: findUh∈Hh such that for all t∈0,T,
(63)AUht,Wht=FWht,∀Wh∈Hh.Since for allt∈0,T, and Wh∈Hh, AUt−Uht,Wht=0, then from (62), we obtain
(64)AUt−Uht,Wt=AUt−Uht,Wt−Wht=AUt,Wt−Wht−AUht,Wt−Wht=FWt−Wht−AUht,Wt−Wht=ff,vft−vf,htΩf+fp,ξpt−ξp,htΩp+qf,wft−wf,htΩf+qp,wpt−wp,htΩp−AUht,Wt−Wht=ff−ff,h,vft−vf,htΩf+fp−fp,h,ξpt−ξp,htΩp+qf−qf,h,wft−wf,htΩf+qp−qp,h,wpt−wp,htΩp+ff,h,vft−vf,htΩf+fp,h,ξpt−ξp,htΩp+qf,h,wft−wf,htΩf+qp,h,wpt−wp,htΩp−AUht,Wt−Wht.Now, we set,(65)AUt−Uht,Wt=∑K∈hAKUt−Uht,Wt.Then, from (64), we can write:
(66)AKUt−Uht,Wt≔ff−ff,h,vft−vf,htΩf∩K+fp−fp,h,ξpt−ξp,htΩp∩K+qf−qf,h,wft−wf,htΩf∩K+qp−qp,h,wpt−wp,htΩp∩K+ff,h,vft−vf,htΩf∩K+fp,h,ξpt−ξp,htΩp∩K+qf,h,wft−wf,htΩf∩K+qp,h,wpt−wp,htΩp∩K−AKUht,Wt−Wht=Rf,KUht−rf,KUht,vft−vf,htΩf∩K+Rf,KUht−rf,KUht,wft−wf,htΩf∩K+Rp,K,1Uht−rp,K,1Uht,ξpt−ξp,htΩp∩K+Rp,KUht−rp,KUht,wpt−wp,htΩp∩K+BKUht,Wt−Wht,where
(67)BKUht,Wt−Wht=ff,K,vft−vf,htΩf∩K+fp,K,ξpt−ξp,htΩp∩K+qf,K,wft−wf,htΩf∩K+qp,K,wpt−wp,htΩp∩K−AKUht,Wt−Wht.


### 3.4. Reliability of the A Posteriori Error Estimator

The first main result is given by the following theorem:Theorem 11 (Upper Error Bound).
LetU=uf,pf,up,pp,ηp,λ∈H be the exact solution and Uh=uf,h,pf,h,up,h,pp,h,ηp,h,λh∈Hh be the finite element solution. Then, there exist a positive constant Crel such that the error is bounded globally from above by:
(68)∥U−Uh∥Hh≤CrelΘUh+ζ.Proof.
Lett∈0;T. We consider the residue Equation (65), and we set Wt=vft,wft,vpt,wpt,ξpt,μt. We take Wht=vf,ht,0,vp,ht,0,ξp,ht,0 with vf,ht=ICl0vft and ξp,ht=ICl0ξpt. As vpt∈H;Ωp, then by Lemma 3.4, vpt admits the decomposition vpt=wpt+βt where wpt∈H1Ωpd and βt∈H1Ωp with ∫Ωpβtxdx=0 and ∥wp∥L∞0,T;H1Ωpd+∥β∥L∞0,T;H1Ωp≤Cp∥vp∥L∞0,T;Vp. We consider vp,ht=wp,ht+βp,ht with wp,ht=IClwpt and βp,ht=IClβpt. Thus, vpt−vp,ht=wpt−wp,ht+βpt−βp,ht. Now, we recall:
(69)BKUht,Wt−Wht=ff,K,vft−vf,htΩf∩K+fp,K,ξpt−ξp,htΩp∩K+qf,K,wft−wf,htΩf∩K+qp,K,wpt−wp,htΩp∩K−AKUht,Wt−Wht.
Therefore, integrate by parts element by element we may write:(70)BKUht,Wt−Wht=∑i=112ηi,K,where
(71)η1,K≔rf,KUht,vft−vf,htΩf∩K,η2,K≔rp,K,1Uht,ξpt−ξp,htΩp∩K,η3,K≔−rp,K,2Uht,wpt−wp,htΩp∩K,η4,K≔−rp,K,2Uht,βpt−βp,htΩp∩K,η5,K≔rf,KUht,wfΩf∩K,η6,K≔rp,KUht,wpΩp∩K,η7,K≔−∑E∈hΩf∩KJE,nE,fUht,vft−vf,htE,η8,K≔−∑E∈hΩp∩KJE,nE,pUht,ξpt−ξp,htE,η9,K≔∑E∈hΓfpRE,pf,1Uht,vft−vf,ht·nf,EE,η10,K≔∑E∈hΓfpRE,pf,2Uht,vft−vf,ht·nf,EE,η11,K≔∑E∈hΓfpRE,pf,3Uht,ξpt−ξp,htE,η12,K≔−∑E∈hΓfp∑j=1d−1RE,pf,4Uhtj,vft−vf,ht·τjE.
Thus, from Equations (66) and (70), we obtain
(72)AKUt−Uht,Wt≔Rf,KUht−rf,KUht,vft−vf,htΩf∩K+Rf,KUht−rf,KUht,wft−wf,htΩf∩K+Rp,K,1Uht−rp,K,1Uht,ξpt−ξp,htΩp∩K+Rp,KUht−rp,KUht,wpt−wp,htΩp∩K+∑i=112ηi,K.
Coercivity of operatorA leads to inf-sup condition:
(73)∥U−Uh∥Hh≤supW∈H∣AU−Uh,W∣∥W∥H.
By consequently, the identities (65) and (72), inf-sup condition of operator A (73), Cauchy-Schwarz inequality, estimation of Lemma 10, and the approximation properties of Lemma 9 imply the required estimate and finish the proof.


### 3.5. Efficiency of the A Posteriori Error Estimator

ForK∈h, we set HhK=L∞0,T;Vf,hK×L∞0,T;Wf,hK×L∞0,T;Vp,hK×W1,∞0,T;Wp,hK×W1,∞0,T;Xp,hK×L∞0,T;ΛhK,where
(74)Vf,hK≔vf,h∣K:vf,h∈Vf,h,Wf,hK≔wf,h∣K:wf,h∈Wf,h,Vp,hK≔vp,h∣K:vp,h∈Vp,h,Wp,hK≔wp,h∣K:wp,h∈Wp,h,Xp,hK≔ξp,h∣K:ξp,h∈Xp,h,ΛhK≔λh∣K:λh∈Λh.The error estimatorΘUh is considered efficient if it satisfies the following theorem.Theorem 12 (Lower Error Bound).
LetU=uf,pf,up,pp,ηp,λ∈H be the exact solution and Uh=uf,h,pf,h,up,h,pp,h,ηp,h,λh∈Hh be the finite element solution. Then, there exist a positive constant Ceff such that the error is bounded locally from below for all K∈h by:
(75)ΘKUh≤Ceff∥U−Uh∥HhK+∑K′⊂w~KζK′,where w~K is a finite union of neighborning elements of K and ∥·∥HhK is the product norm.Proof.
The lower bound are proved using standard elementwise integration by parts, namely, error equation of Section3.3 (i.e., identity (65) and Equation (72)) and some a inverse estimates of Lemma 7 (Cf. [37] for details).


## 3.1. Residual Error Estimators

The general philosophy of residual error estimators is to estimate an appropriate norm of the correct residual by terms that can be evaluated easier and that involve the data at hand. To this end, define the exact element residuals:Definition 3 (Exact Element Residuals).
Lett∈0,T and Wht=vf,ht,pf,ht,vp,ht,wp,ht,ξp,ht,λht∈Vf,h×Wf,h×Vp,h×Wp,h×Xp,h×Λh be an arbitrary finite element function. The exact element residuals over a triangle or tetrahedra K∈h and over E∈hΓfp are defined for all t∈0,T by:
(42)Rf,KWht=ff+∇·σfvf,ht,wf,ht,Rp,K,1Wht=fp+∇·σpξp,ht,wp,ht,Rp,K,2Wht=μK−1vp,ht+∇wp,ht,Rf,KWht=qf−∇·vf,ht,Rp,KWht=qp−∂ts0wp,ht+α∇·ξp,ht+∇·vp,ht,RE,pf,1Wht=vf,ht·nf,E+∂tξpt+vp,ht·np,E,RE,pf,2Wht=wp,h+σfvf,ht,wf,htnf,E·nf,E,RE,pf,3Wht=σfvf,ht,wf,htnf,E+σpξp,ht,wp,htnp,E,RE,pf,4Whtj=σfvf,h,wf,hτf,E,j+μαBJSKj−1vf,h−∂tξp,h·τf,E,j.As it is common, these exact residuals are replaced by some finite-dimensional approximation called approximate element residualr∗,K, r∗,K, ∗∈f,p, RE,pf,l, l∈1,2,3,4. This approximation is here achieved by projecting ff and qf on the space of piecewise constant functions in Ωf and piecewise ℙ1 functions in Ωp, more precisely for all K∈hf we take
(43)ff,K=pKff,qf,K=PKqf,withPK:L1K→ℝsuchthatpKϕ=1∣K∣∫Kϕxdx.While for allK∈hp, we take fp,K and qp as the unique element of ℙ1Kd, respectively ℙ1K such that:
(44)∫Kfp,Kx·qxdx=∫Kfpxqxdx,q∈ℙ1Kd,respectively,
(45)∫Kqp,Kxqxdx=∫Kqpxqxdx,∀q∈ℙ1K.Thereby, we define the approximate element residuals.Definition 4 [Approximate Element Residuals].
Lett∈0,T and Wht=vf,ht,pf,ht,vp,ht,wp,ht,ξp,ht,λht∈Vf,h×Wf,h×Vp,h×Wp,h×Xp,h×Λh be an arbitrary finite element function. Then, the approximate element residuals are defined for all t∈0,T by:
(46)rf,KWht=ff,K+∇·σfvf,ht,wf,ht,rp,K,1Wht=fp,K+∇·σpξp,ht,wp,ht,rp,K,2Wht=μK−1vp,ht+∇wp,ht,rf,KWht=qf,K−∇·vf,ht,rp,KWht=qp,K−∂ts0wp,ht+α∇·ξp,ht+∇·vp,ht.Next, introduce the gradient jump in normal direction by(47)JE,nE,fUh≔2μDuf,h−pf,hI·nEEifE∈hΩf,JE,nE,pUhJE,nE,pUh≔2μDηp,h−pp,hI·nEEifE∈hΩp.where I is the identity matrix of ℝd×d.Definition 5 [Residual Error Estimators].
LetUh=uf,h,pf,h,up,h,pp,h,ηp,h,λh be the finite element solution of the problem (36)–(38) in L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh. Then, the residual error estimator is locally defined by
(48)ΘKUh≔ΘK,f2Uh+ΘK,p2Uh+ΘK,pf2Uh1/2,where
(49)ΘK,f2Uh≔hK2∥rf,KUh∥L∞0,T;L2K2+∥rf,KUh∥L∞0,T;L2K2+∑E∈hΩfhE∥JE,nE,fUh∥L∞0,T;L2E2,(50)ΘK,p2Uh≔hK2∥rp,K,1Uh∥L∞0,T;L2K2+∥rp,K,2Uh∥L∞0,T;L2K2+hK2∥rp,K,2Uh∥L∞0,T;L2K2+∥rp,KUh∥L∞0,T;L2K2+∑E∈hΩphE∥JE,nE,pUh∥L∞0,T;L2E2,(51)ΘK,pf2Uh≔∑E∈hΓpfhE∥RE,pf,1Uh∥L∞0,T;L2E2+∑E∈hΓpfhE∥RE,pf,2Uh∥L∞0,T;L2E2+∑E∈hΓpfhE∥RE,pf,3Uh∥L∞0,T;L2E2+∑E∈hΓpfhE∑j=1d−1∥RE,pf,4j∥L∞0,T,L2E2.The global residual error estimator is given by(52)ΘUh≔∑K∈hΘKUh21/2.Furthermore, denote the local approximation terms by(53)ζK2≔ζK,f2+ζK,p2,where
(54)ζK,f2≔hK2∥Rf,KUh−rf,KUh∥L∞0,T;L2K2+hK2∥Rf,KUh−rf,KUh∥L∞0,T,L2K2ifK∈hf,ζK,p2≔hK2∥Rp,K,1Uh−rp,K,1Uh∥L∞0,T;L2K2+hK2∥Rp,KUh−rp,KUh∥L∞0,T,L2K2ifK∈hp.The global approximation term is defined by(55)ζ≔∑K∈hζK21/2.Remark 6.
The residual character of each term on the right-hand sides of (48)–(51) is quite clear since if Uh would be the exact solution of (2)–(10), then they would vanish.


## 3.2. Analytical Tools

### 3.2.1. Inverse Inequalities

In order to derive the lower error bounds, we proceed similarly as in [29, 33] (see also [34]), by applying inverse inequalities, and the localization technique based on simplex-bubble and face-bubble functions. To this end, we recall some notation and introduce further preliminary results. Given K∈Th, and E∈K, we let bK and bE be the usual simplexe-bubble and face-bubble functions, respectively (see (1.5) and (1.6) in [28]). In particular, bK satisfies bK∈ℙ3K, suppbK⊆K, bK=0on∂K, and 0≤bK≤1onK. Similarly, bE∈ℙ2K, suppbE⊆ωE≔K′∈Th:E∈K′, bE=0on∂K\E and 0≤bE≤1inωE. We also recall from [35] that, given k∈ℕ, there exists an extension operator L:CE→CK that satisfies Lp∈ℙkK and Lp∣E=p,∀p∈ℙkE. A corresponding vectorial version of L, that is, the componentwise application of L, is denoted by L. Additional properties of bK, bE, and L are collected in the following lemma (see [35]):Lemma 7.
Givenk∈ℕ∗, there exist positive constants depending only on k and shape regularity of the triangulations (minimum angle condition), such that for each simplexe K and E∈K there hold
(56)∥ϕ∥K≲∥ϕbK1/2∥K≲∥ϕ∥K,∀ϕ∈ℙkK,ϕbK1,K≲hK−1∥ϕ∥K,∀ϕ∈ℙkK,∥ψ∥E≲∥bE1/2ψ∥E≲∥ψ∥E,∀ψ∈ℙkE,∥Lψ∥K+hELψ1,K≲hE1/2∥ψ∥E∀ψ∈ℙkE.Lemma 8 (Continuous trace inequality).
There exists a positive constantβ1>0 depending only on σ0 such that
(57)∥v∥∂K2≤β1∥v∥K∥v∥1,K,∀K∈Th,∀v∈H1Kd.


### 3.2.2. Clément Interpolation Operator

In order to derive the upper error bounds, we introduce the Clément interpolation operatorICl0:H01Ω→PcbTh that approximates optimally nonsmooth functions by continuous piecewise linear functions:
(58)PcbTh≔v∈C0Ω¯:v∣K∈ℙ1K,∀K∈Thandv=0on∂Ω.In addition, we will make use of a vector valued version ofICl0, that is,ICl0:H01Ωd→PcbThd, which is defined componentwise by ICl0. The following lemma establishes the local approximation properties of ICl0 (and hence of ICl0), for a proof see ([36], Section 3).Lemma 9.
There exist constantsC1,C2>0, independent of h, such that for all v∈H01Ω there hold
(59)∥v−ICl0v∥K≤C1hK∥v∥1,ΔK,∀K∈Th,∥v−ICl0v∥E≤C2hE1/2∥v∥1,ΔE,∀E∈h,
whereΔK≔∪K′∈Th:K′∩K≠∅ and ΔE≔∪K′∈Th:K′∩E≠∅.


### 3.2.3. Helmholtz Decomposition

Lemma 10 ([7]).
There existsCp>0 such that every vp∈H;Ωp can be decomposed as vp=w+β, where w∈H1Ωpd, β∈H1Ωp, ∫Ωpβxdx=0 and
(60)∥w∥1,Ωp+∥β∥1,Ωp≤Cp∥vp∥H;Ωp.


## 3.2.1. Inverse Inequalities

In order to derive the lower error bounds, we proceed similarly as in [29, 33] (see also [34]), by applying inverse inequalities, and the localization technique based on simplex-bubble and face-bubble functions. To this end, we recall some notation and introduce further preliminary results. Given K∈Th, and E∈K, we let bK and bE be the usual simplexe-bubble and face-bubble functions, respectively (see (1.5) and (1.6) in [28]). In particular, bK satisfies bK∈ℙ3K, suppbK⊆K, bK=0on∂K, and 0≤bK≤1onK. Similarly, bE∈ℙ2K, suppbE⊆ωE≔K′∈Th:E∈K′, bE=0on∂K\E and 0≤bE≤1inωE. We also recall from [35] that, given k∈ℕ, there exists an extension operator L:CE→CK that satisfies Lp∈ℙkK and Lp∣E=p,∀p∈ℙkE. A corresponding vectorial version of L, that is, the componentwise application of L, is denoted by L. Additional properties of bK, bE, and L are collected in the following lemma (see [35]):Lemma 7.
Givenk∈ℕ∗, there exist positive constants depending only on k and shape regularity of the triangulations (minimum angle condition), such that for each simplexe K and E∈K there hold
(56)∥ϕ∥K≲∥ϕbK1/2∥K≲∥ϕ∥K,∀ϕ∈ℙkK,ϕbK1,K≲hK−1∥ϕ∥K,∀ϕ∈ℙkK,∥ψ∥E≲∥bE1/2ψ∥E≲∥ψ∥E,∀ψ∈ℙkE,∥Lψ∥K+hELψ1,K≲hE1/2∥ψ∥E∀ψ∈ℙkE.Lemma 8 (Continuous trace inequality).
There exists a positive constantβ1>0 depending only on σ0 such that
(57)∥v∥∂K2≤β1∥v∥K∥v∥1,K,∀K∈Th,∀v∈H1Kd.


## 3.2.2. Clément Interpolation Operator

In order to derive the upper error bounds, we introduce the Clément interpolation operatorICl0:H01Ω→PcbTh that approximates optimally nonsmooth functions by continuous piecewise linear functions:
(58)PcbTh≔v∈C0Ω¯:v∣K∈ℙ1K,∀K∈Thandv=0on∂Ω.In addition, we will make use of a vector valued version ofICl0, that is,ICl0:H01Ωd→PcbThd, which is defined componentwise by ICl0. The following lemma establishes the local approximation properties of ICl0 (and hence of ICl0), for a proof see ([36], Section 3).Lemma 9.
There exist constantsC1,C2>0, independent of h, such that for all v∈H01Ω there hold
(59)∥v−ICl0v∥K≤C1hK∥v∥1,ΔK,∀K∈Th,∥v−ICl0v∥E≤C2hE1/2∥v∥1,ΔE,∀E∈h,
whereΔK≔∪K′∈Th:K′∩K≠∅ and ΔE≔∪K′∈Th:K′∩E≠∅.


## 3.2.3. Helmholtz Decomposition

Lemma 10 ([7]).
There existsCp>0 such that every vp∈H;Ωp can be decomposed as vp=w+β, where w∈H1Ωpd, β∈H1Ωp, ∫Ωpβxdx=0 and
(60)∥w∥1,Ωp+∥β∥1,Ωp≤Cp∥vp∥H;Ωp.


## 3.3. Error Equation

We setH=L∞0,T;Vf×L∞0,T;Wf×L∞0,T;Vp×W1,∞0,T;Wp×W1,∞0,T;Xp×L∞0,T;Λ and Hh=L∞0,T;Vf,h×L∞0,T;Wf,h×L∞0,T;Vp,h×W1,∞0,T;Wp,h×W1,∞0,T;Xp,h×L∞0,T;Λh.LetU=uf,pf,up,pp,ηp,λ∈H and W=vf,wf,vp,wp,ξp,μ∈H. For t∈0,T, we define the operator A by
(61)AUt,Wt≔afuft,vft+apdupt,vpt+bpvpt,ppt+αbpξpt,ppt+bΓvft,vpt,ξpt,λt+s0∂tppt,wptΩp−αbp∂tηpt,wpt−bpupt,wpt−bfuft,wft+bΓuft,upt,∂tηpt;μ,FWt≔ff,vftΩf+fp,ξpt+qf,wftΩf+qp,wptΩp.Then, the continuous problem (25)–(27) is equivalent to find U∈H such that for all t∈0,T, we have:
(62)AUt,Wt=FWt,∀W∈H.We define the discrete version by the same way: findUh∈Hh such that for all t∈0,T,
(63)AUht,Wht=FWht,∀Wh∈Hh.Since for allt∈0,T, and Wh∈Hh, AUt−Uht,Wht=0, then from (62), we obtain
(64)AUt−Uht,Wt=AUt−Uht,Wt−Wht=AUt,Wt−Wht−AUht,Wt−Wht=FWt−Wht−AUht,Wt−Wht=ff,vft−vf,htΩf+fp,ξpt−ξp,htΩp+qf,wft−wf,htΩf+qp,wpt−wp,htΩp−AUht,Wt−Wht=ff−ff,h,vft−vf,htΩf+fp−fp,h,ξpt−ξp,htΩp+qf−qf,h,wft−wf,htΩf+qp−qp,h,wpt−wp,htΩp+ff,h,vft−vf,htΩf+fp,h,ξpt−ξp,htΩp+qf,h,wft−wf,htΩf+qp,h,wpt−wp,htΩp−AUht,Wt−Wht.Now, we set,(65)AUt−Uht,Wt=∑K∈hAKUt−Uht,Wt.Then, from (64), we can write:
(66)AKUt−Uht,Wt≔ff−ff,h,vft−vf,htΩf∩K+fp−fp,h,ξpt−ξp,htΩp∩K+qf−qf,h,wft−wf,htΩf∩K+qp−qp,h,wpt−wp,htΩp∩K+ff,h,vft−vf,htΩf∩K+fp,h,ξpt−ξp,htΩp∩K+qf,h,wft−wf,htΩf∩K+qp,h,wpt−wp,htΩp∩K−AKUht,Wt−Wht=Rf,KUht−rf,KUht,vft−vf,htΩf∩K+Rf,KUht−rf,KUht,wft−wf,htΩf∩K+Rp,K,1Uht−rp,K,1Uht,ξpt−ξp,htΩp∩K+Rp,KUht−rp,KUht,wpt−wp,htΩp∩K+BKUht,Wt−Wht,where
(67)BKUht,Wt−Wht=ff,K,vft−vf,htΩf∩K+fp,K,ξpt−ξp,htΩp∩K+qf,K,wft−wf,htΩf∩K+qp,K,wpt−wp,htΩp∩K−AKUht,Wt−Wht.


## 3.4. Reliability of the A Posteriori Error Estimator

The first main result is given by the following theorem:Theorem 11 (Upper Error Bound).
LetU=uf,pf,up,pp,ηp,λ∈H be the exact solution and Uh=uf,h,pf,h,up,h,pp,h,ηp,h,λh∈Hh be the finite element solution. Then, there exist a positive constant Crel such that the error is bounded globally from above by:
(68)∥U−Uh∥Hh≤CrelΘUh+ζ.Proof.
Lett∈0;T. We consider the residue Equation (65), and we set Wt=vft,wft,vpt,wpt,ξpt,μt. We take Wht=vf,ht,0,vp,ht,0,ξp,ht,0 with vf,ht=ICl0vft and ξp,ht=ICl0ξpt. As vpt∈H;Ωp, then by Lemma 3.4, vpt admits the decomposition vpt=wpt+βt where wpt∈H1Ωpd and βt∈H1Ωp with ∫Ωpβtxdx=0 and ∥wp∥L∞0,T;H1Ωpd+∥β∥L∞0,T;H1Ωp≤Cp∥vp∥L∞0,T;Vp. We consider vp,ht=wp,ht+βp,ht with wp,ht=IClwpt and βp,ht=IClβpt. Thus, vpt−vp,ht=wpt−wp,ht+βpt−βp,ht. Now, we recall:
(69)BKUht,Wt−Wht=ff,K,vft−vf,htΩf∩K+fp,K,ξpt−ξp,htΩp∩K+qf,K,wft−wf,htΩf∩K+qp,K,wpt−wp,htΩp∩K−AKUht,Wt−Wht.
Therefore, integrate by parts element by element we may write:(70)BKUht,Wt−Wht=∑i=112ηi,K,where
(71)η1,K≔rf,KUht,vft−vf,htΩf∩K,η2,K≔rp,K,1Uht,ξpt−ξp,htΩp∩K,η3,K≔−rp,K,2Uht,wpt−wp,htΩp∩K,η4,K≔−rp,K,2Uht,βpt−βp,htΩp∩K,η5,K≔rf,KUht,wfΩf∩K,η6,K≔rp,KUht,wpΩp∩K,η7,K≔−∑E∈hΩf∩KJE,nE,fUht,vft−vf,htE,η8,K≔−∑E∈hΩp∩KJE,nE,pUht,ξpt−ξp,htE,η9,K≔∑E∈hΓfpRE,pf,1Uht,vft−vf,ht·nf,EE,η10,K≔∑E∈hΓfpRE,pf,2Uht,vft−vf,ht·nf,EE,η11,K≔∑E∈hΓfpRE,pf,3Uht,ξpt−ξp,htE,η12,K≔−∑E∈hΓfp∑j=1d−1RE,pf,4Uhtj,vft−vf,ht·τjE.
Thus, from Equations (66) and (70), we obtain
(72)AKUt−Uht,Wt≔Rf,KUht−rf,KUht,vft−vf,htΩf∩K+Rf,KUht−rf,KUht,wft−wf,htΩf∩K+Rp,K,1Uht−rp,K,1Uht,ξpt−ξp,htΩp∩K+Rp,KUht−rp,KUht,wpt−wp,htΩp∩K+∑i=112ηi,K.
Coercivity of operatorA leads to inf-sup condition:
(73)∥U−Uh∥Hh≤supW∈H∣AU−Uh,W∣∥W∥H.
By consequently, the identities (65) and (72), inf-sup condition of operator A (73), Cauchy-Schwarz inequality, estimation of Lemma 10, and the approximation properties of Lemma 9 imply the required estimate and finish the proof.


## 3.5. Efficiency of the A Posteriori Error Estimator

ForK∈h, we set HhK=L∞0,T;Vf,hK×L∞0,T;Wf,hK×L∞0,T;Vp,hK×W1,∞0,T;Wp,hK×W1,∞0,T;Xp,hK×L∞0,T;ΛhK,where
(74)Vf,hK≔vf,h∣K:vf,h∈Vf,h,Wf,hK≔wf,h∣K:wf,h∈Wf,h,Vp,hK≔vp,h∣K:vp,h∈Vp,h,Wp,hK≔wp,h∣K:wp,h∈Wp,h,Xp,hK≔ξp,h∣K:ξp,h∈Xp,h,ΛhK≔λh∣K:λh∈Λh.The error estimatorΘUh is considered efficient if it satisfies the following theorem.Theorem 12 (Lower Error Bound).
LetU=uf,pf,up,pp,ηp,λ∈H be the exact solution and Uh=uf,h,pf,h,up,h,pp,h,ηp,h,λh∈Hh be the finite element solution. Then, there exist a positive constant Ceff such that the error is bounded locally from below for all K∈h by:
(75)ΘKUh≤Ceff∥U−Uh∥HhK+∑K′⊂w~KζK′,where w~K is a finite union of neighborning elements of K and ∥·∥HhK is the product norm.Proof.
The lower bound are proved using standard elementwise integration by parts, namely, error equation of Section3.3 (i.e., identity (65) and Equation (72)) and some a inverse estimates of Lemma 7 (Cf. [37] for details).


## 4. Discussion

In this paper, we have discussed a posteriori error estimates for a finite element approximation of the Stokes-Biot system. A residual type a posteriori error estimator is provided, that is both reliable and efficient.There are many issues to be addressed in this area such as the other types of a posteriori error estimates and related implementation of the adaptive finite element methods.Further, it is well known that an internal layer appears at the interfaceΓfp as the permeability tensor degenerates; in that case, anisotropic meshes have to be used in this layer (see for instance [10]). Hence, we intend to extend our results to such anisotropic meshes.


## 5. Nomenclature

(i)
Ω⊂ℝd,d∈2,3 bounded domain(ii)
Ωp: The poroelastic medium domain(iii)
Ωf=Ω\Ω¯d(iv)
Γfp=∂Ωf∩∂Ωp(v)
Γ∗=∂Ω∗\Γfp,∗=f,p(vi)
nf (resp. np): The unit outward normal vector along ∂Ωf (resp. ∂Ωp)(vii)
uf: The fluid velocity in Ωf(viii)
pf: The fluid pressure in Ωf(ix)
up,ηp: The fluid velocities in Ωp(x)
pp: The fluid pressure in Ωp(xi)
In2D, the curl of a scalar function w is given as usual by
(76)curlw≔∂w∂x2,−∂w∂x1Τ(xii)
In3D, the of a vector function w=w1,w2,w3 is given as usual by w≔∇×w, namely,
(77)curlw≔∂w3∂x2−∂w2∂x3,∂w1∂x3−∂w3∂x1,∂w2∂x1−∂w1∂x2(xiii)
ℙk: The space of polynomials of total degree not larger than k(xiv)
Th: Triangulation of Ω(xv)
Th∗: The corresponding induced triangulation of Ω∗, ∗∈f,p(xvi)
For anyK∈h, hK is the diameter of K, and ρK=2rK is the diameter of the largest ball inscribed into K(xvii)
h≔maxK∈τhhKandσh≔maxK∈τhhK/ρK(xviii)
εh: The set of all the edges or faces of the triangulation(xix)
εK: The set of all the edges (N=2) or faces (N=3) of an element K(xx)
εh≔⋃K∈hK(xxi)
NK: The set of all the vertices of an element K(xxii)
Nh≔⋃K∈hNK(xxiii)
ForA⊂Ω¯, EhA≔E∈Eh:E⊂A(xxiv)
ForE∈hΩ∗, we associate a unit vector nE,∗ such that nE,∗ is orthogonal to E and equals to the unit exterior normal vector to ∂Ω∗, ∗∈f,p(xxv)
ForE∈hΩ∗, ϕE is the jump across E in the direction of nE,∗(xxvi)
In order to avoid excessive use of constants, the abbreviationsx≲y and x~y stand for x⩽cy and c1x⩽y⩽c2x, respectively, with positive constants independent of x, y or h(xxvii)
∂t=∂/∂t(xxviii)
∥ϕ∥L20,T;X≔∫0T∥ϕs∥X2ds1/2(xxix)
∥ϕ∥L∞0,T;X≔sup0<t≤T∥ϕt∥X(xxx)
∥ϕ∥W1,∞0,T;X≔max∥ϕ∥L∞0,T;X,∥∂tϕ∥L∞0,T;X(xxxi)
∥·∥Hh: Product norm on Hh(xxxii)
∥·∥HhK: Local product norm on HhK.



---

*Source: 8877012-2021-03-04.xml*