Partitioning strategies for the interaction of a fluid with a poroelastic material based on a Nitsche's coupling approach

被引:68
作者
Bukac, M. [1 ]
Yotov, I. [2 ]
Zakerzadeh, R. [3 ]
Zunino, P. [3 ]
机构
[1] Univ Notre Dame, Dept Appl & Computat Math & Stat, Notre Dame, IN 46556 USA
[2] Univ Pittsburgh, Dept Math, Pittsburgh, PA 15260 USA
[3] Univ Pittsburgh, Dept Mech Engn & Mat Sci, Pittsburgh, PA 15261 USA
基金
美国国家科学基金会;
关键词
Fluid-structure interaction; Poroelasticity; Operator-splitting scheme; Nitsche's method; Preconditioning; FINITE-ELEMENT-METHOD; FLOW; STOKES; FRACTURES; STABILITY; SURFACE; MODELS;
D O I
10.1016/j.cma.2014.10.047
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
We develop a computational model to study the interaction of a fluid with a poroelastic material. The coupling of Stokes and Biot equations represents a prototype problem for these phenomena, which feature multiple facets. On one hand, it shares common traits with fluid-structure interaction. On the other hand it resembles the Stokes-Darcy coupling. For these reasons, the numerical simulation of the Stokes-Biot coupled system is a challenging task. The need of large memory storage and the difficulty to characterize appropriate solvers and related preconditioners for the equations at hand are typical shortcomings of classical discretization methods applied to this problem, such as the finite element method for spatial discretization and finite differences for time stepping. The application of loosely coupled time advancing schemes mitigates these issues, because it allows to solve each equation of the system independently with respect to the others, at each time step. In this work, we develop and thoroughly analyze a loosely coupled scheme for Stokes-Biot equations. The scheme is based on Nitsche's method for enforcing interface conditions. Once the interface operators corresponding to the interface conditions have been defined, time lagging allows us to build up a loosely coupled scheme with good stability properties. The stability of the scheme is guaranteed provided that appropriate stabilization operators are introduced into the variational formulation of each subproblem. The error of the resulting method is also analyzed, showing that splitting the equations pollutes the optimal approximation properties of the underlying discretization schemes. In order to restore good approximation properties, while maintaining the computational efficiency of the loosely coupled approach, we consider the application of the loosely coupled scheme as a preconditioner for the monolithic approach. Both theoretical insight and numerical results confirm that this is a promising way to develop efficient solvers for the problem at hand. (C) 2014 Elsevier B.V. All rights reserved.
引用
收藏
页码:138 / 170
页数:33
相关论文
共 42 条
  • [1] Coupling Biot and Navier-Stokes equations for modelling fluid-poroelastic media interaction
    Badia, Santiago
    Quaini, Annalisa
    Quarteroni, Alfio
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 228 (21) : 7986 - 8014
  • [2] Brezzi Franco, 1984, Efficient Solutions of Elliptic Systems, P11
  • [3] Bukac M., 2014, MODELING SIMULATION, V13
  • [4] Bukac M., 2014, NUMER METHODS PARTIA
  • [5] Bukac M., 2013, ANAL PARTITION UNPUB
  • [6] Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement
    Bukac, Martina
    Canic, Suncica
    Glowinski, Roland
    Tambaca, Josip
    Quaini, Annalisa
    [J]. JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 235 : 515 - 541
  • [7] Explicit strategies for incompressible fluid-structure interaction problems: Nitsche type mortaring versus Robin-Robin coupling
    Burman, Erik
    Fernandez, Miguel A.
    [J]. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2014, 97 (10) : 739 - 758
  • [8] Stabilization of explicit coupling in fluid-structure interaction involving fluid incompressibility
    Burman, Erik
    Fernandez, Miguel A.
    [J]. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2009, 198 (5-8) : 766 - 784
  • [9] FINITE ELEMENT APPROXIMATIONS FOR STOKES-DARCY FLOW WITH BEAVERS-JOSEPH INTERFACE CONDITIONS
    Cao, Yanzhao
    Gunzburger, Max
    Hu, Xiaolong
    Hua, Fei
    Wang, Xiaoming
    Zhao, Weidong
    [J]. SIAM JOURNAL ON NUMERICAL ANALYSIS, 2010, 47 (06) : 4239 - 4256
  • [10] Cao YZ, 2010, COMMUN MATH SCI, V8, P1