A semi-decoupled MAC scheme for the coupled fluid-poroelastic material interaction

https://doi.org/10.1016/j.camwa.2023.04.003Get rights and content

Abstract

In this paper, a semi-decoupled marker and cell (MAC) scheme is proposed for the Stokes-Biot system, in which the displacement of structure is split from the whole system using the time-lagging scheme, then the rest system composed by the velocity and pressure is treated in a similar way as the Stokes-Darcy coupling problem. It keeps the desirable qualities such as the mass conservation. The backward Euler scheme is used for the time discretization. We derive the stability and the error estimates for the fully discrete scheme, obtain the second order superconvergence for the discrete L2 norm for velocity, pressure and displacement and some terms of the H1 norm for displacement and velocity. Numerical experiments confirm well the theoretical analysis results and the robustness of parameters.

Introduction

The interaction between the incompressible viscous flow and the flow in saturated porous media with small deformable is important in many engineering and biological applications, such as reservoir engineering, groundwater flow, blood flow through vessels [1], [2], [3], [4]. We consider a fluid-poroelastic structure interaction (FPSI) problem, in which the fluid is described by the Stokes equations and the structure is provided by the fully dynamic Biot model, that is, coupled Stokes-Biot problem. The model features two different couplings, the Stokes-Darcy (SD) problem [5], [6], [7], [8], [9], [10] and the fluid-structure interaction (FSI) problem [11], [12], [13], [14]. It carries on all difficulties involving SD and FSI problems. Due to complexity of interface conditions, the analysis in SD cannot directly be used in FPSI. Thus, numerical methods for solving FPSI system paly a significant role.

The Stokes-Biot problem has been analyzed to be well-posed [15]. Concerning the computational perspective, different numerical approximates have attracted much attention of scholars. There are mainly two ways for solving the coupled Stokes-Biot flow model, including the monolithic method [16], [17], [18], [19], [20], [21], [22], [23] and partitioned approach [24], [25], [26], [27], [28], [29]. For example, Ambartsumyan et al. focused on the monolithic scheme for Stoke-Biot problem with the approximation of the mass conservation via a Lagrange multiplier [18]. Wen et al. constructed a strongly conservative discrete scheme [21], and proposed the interior penalty discontinuous finite element method [22]. Due to the monolithic approach requires solving a large linear system, it is computationally complex; obviously, to develop an effective and stable partitioned method to save computer storage and reduce executive time is natural and necessary. Buskač et al. developed a loosely coupled finite element scheme based on Lie operator splitting method to separate the Navier-Stokes equations from the Biot problem [24] and studied an alternative partitioned strategies for Stokes-Biot problems based on a Nitsche's method [25]. Moreover, Buskač [26] introduced a kinematically β-coupled scheme to split the fluid problem from the dynamic poroelasticity problem, and derived its stability. Recently, Kunwar [27] discussed the second-order time discretization optimization-based decoupling algorithms for Stokes-Biot system. Oyekole and Buskač [28] constructed two types of second-order loosely coupled schemes for Stokes/Biot problem, and proved the stability of them. Guo et al. [29] performed the stability and error estimates for a decoupled modified characteristic finite element scheme for Navier-Stokes/Biot problem. These partitioned methods are mainly based on the finite element methods where the interface conditions are subtly embedded in the weak form, which is not feasible in all applications. In this paper, we will concentrate on the partitioned approach based on the MAC finite difference scheme.

The MAC scheme has been a hot research topic for its fascinating properties such as satisfying the discrete incompressibility constraint, as well as locally conserving of the mass, momentum and energy [30], [31]. In recent years, MAC method has been used to solve Darcy-Forchheimer equation [32], Maxwell's equation [33], Stokes problem [34], linear elasticity problem [35], Stokes-Darcy problem [5], [7]. Rui et al. [5] constructed the second order MAC scheme for the SD system by using the techniques of integration by parts on some controlling volumes near the interface. We constructed a locking free finite difference method using the similar techniques for Stokes-Biot problem, established the stability and the error estimates [23]. However, it needs solve a large system which costs a lot of CPU time. The aim of this work is to develop and analyze an efficient semi-decoupled numerical solver for the Stokes-Biot system. More precisely, we split the displacement of structure from the whole system by designing a time-lagging scheme at each time step, then the residual system composed of velocity and pressure is similar to the coupled SD model. However, we need deal with the more complex interface conditions. In this paper, we rigorously and carefully derive the stability and error estimates by choosing the proper interface conditions. We obtain the second order convergence in the discrete L2 norm for the velocity, pressure and displacement and the second order superconvergence for some terms of the H1 norm for the velocity and displacement.

The paper is organized as follows. Section 2 is devoted to showing the Stokes-Biot equations and the semi-decoupled MAC scheme. The analysis of stability is derived in Section 3. An error estimate is carried out in Section 4. The numerical simulations are shown in Section 5. Throughout the paper we use C, whether it has subscrip or not, to represent a positive constant. It is independent of the spatial mesh size or the time step.

Section snippets

Model, notations and semi-decoupled MAC scheme

In this section, we consider the coupled Stokes-Biot model. Without loss of generality, we restrict the model to a two-dimensional (2D) geometrical model Ω=ΩfΩp. The fluid domain Ωf and poroelastic material region Ωp in our problem are given, respectively, byΩf:={(x,y)|0<x<a¯,0<y<b},Ωp:={(x,y)|a¯<x<a,0<y<b}, and the boundary Γf=ΩfΩ,Γp=ΩpΩ, the interface ΓI=ΩfΩp.

Let T be the final time. We model the free flow using the Stokes equations:σf(uf,pf)=ff in Ωf×(0,T],uf=0 in Ωf×(0,T],

The analysis of stability

In this section, we consider the discrete LBB condition and the stability analysis for semi-decoupled Stokes-Biot MAC scheme. Before proceeding, let us address the following property that will serve as auxiliary result for the stability and error analysis in Section 3.1. For simplicity, we present the preliminaries of discrete LBB condition in the Appendix. Here, we just show the result of discrete LBB condition. Then we derive the stability of scheme (2.23a), (2.23b), (2.23c), (2.23d), (2.23e)

The error estimates

In this section, we consider the error estimate of the semi-decoupled MAC scheme (2.23a)-(2.25d). First, by using Taylor's expansion, we can easily obtain the following results (Lemma 4.1, Lemma 4.3) for the pressure, velocity and displacement, which will be used in establishing the error analysis. In this section, the subscripts i,j will be dropped out, the grid point (xi,yj) still be denoted by γ. For some special points, we keep the subscripts. We use the standard Sobolev spaces that are

The numerical examples

In the numerical examples, we investigate accuracy of the proposed numerical scheme. We test Example 1, Example 2 to verify the convergence rates for the errors of velocity, pressure and displacement using the semi-decoupled MAC scheme. Here the results on non-uniform grids are similar to on uniform grids, so they are omitted.

Example 1

The exact solution is taken from Reference [21]. We set μ=ρ=λp=μp=α=s0=β=1, and κ=1. The domain Ω=[0,2]×[0,1], the Stokes domain Ωf=[0,1]×[0,1], poroelastic domain Ωp=[1,2]

Acknowledgements

This work is supported by the National Natural Science Foundation of China Grant No. 12131014.

References (41)

  • H. Rui et al.

    A locking-free finite difference method on staggered grids for linear elasticity problems

    Comput. Math. Appl.

    (2018)
  • J. Chen

    Time domain fundamental solution to Biot's complete equations of dynamic poroelasticity. Part II: three-dimensional solution

    Int. J. Solids Struct.

    (1994)
  • V. Calo et al.

    Multiphysics model for blood flow and drug transport with application to patient-specific coronary artery flow

    Comput. Mech.

    (2008)
  • N. Koshiba et al.

    Multiphysics simulation of blood flow and LDL transport in a porohyperelastic arterial wall model

    J. Biomech. Eng.

    (2007)
  • B. Tully et al.

    Coupling poroelasticity and CFD for cerebrospinal fluid hydrodynamics

    IEEE Trans. Biomed. Eng.

    (2009)
  • R. Zakerzadeh et al.

    Computational analysis of energy distribution of coupled blood flow and arterial defomation

    Int. J. Adv. Eng. Sci. Appl. Math.

    (2016)
  • H. Rui et al.

    A MAC scheme for coupled Stokes-Darcy equations on non-uniform grids

    J. Sci. Comput.

    (2020)
  • X. Li et al.

    Superconvergence of MAC scheme for a coupled free flow-porous media system with heat transport on non-uniform grids

    J. Sci. Comput.

    (2022)
  • M. Shiue et al.

    Convergence of the MAC scheme for the Stokes/Darcy coupling problem

    J. Sci. Comput.

    (2018)
  • Y. Cao et al.

    Parallel, non-iterative, multi-physics domain decomposition methods for time-dependent Stokes-Darcy systems

    Math. Comput.

    (2014)
  • Cited by (0)

    View full text