Numerical approximation with Nitsche's coupling of transient Stokes'/Darcy's flow problems applied to hemodynamics

We consider a time dependent coupled Stokes/Darcy flow problem and study an approximation method based on a unified finite element scheme complemented with implicit time stepping. Our finite element formulation relies on a weighing strategy in which the physical and discretization parameters are taken into account to robustly enforce interface and boundary conditions by means of the Nitsche method. We study absolute stability and convergence of the scheme, and discuss the algebraic properties of the associated discrete problem. Finally, we present numerical experiments confirming the predicted convergence behavior and algebraic properties, and report an application to the computational analysis of blood flow and plasma filtration in arteries after the implantation of a vascular graft.

[1]  Kent-André Mardal,et al.  Uniform preconditioners for the time dependent Stokes problem , 2006, Numerische Mathematik.

[2]  Erik Burman,et al.  A Domain Decomposition Method Based on Weighted Interior Penalties for Advection-Diffusion-Reaction Problems , 2006, SIAM J. Numer. Anal..

[3]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[4]  R. Rannacher,et al.  Finite element approximation of the nonstationary Navier-Stokes problem. I : Regularity of solutions and second-order error estimates for spatial discretization , 1982 .

[5]  Shigeru Tada,et al.  Internal elastic lamina affects the distribution of macromolecules in the arterial wall: a computational study. , 2004, American journal of physiology. Heart and circulatory physiology.

[6]  Kumbakonam R. Rajagopal,et al.  ON A HIERARCHY OF APPROXIMATE MODELS FOR FLOWS OF INCOMPRESSIBLE FLUIDS THROUGH POROUS SOLIDS , 2007 .

[7]  William J. Layton,et al.  Partitioned Time Stepping for a Parabolic Two Domain Problem , 2009, SIAM J. Numer. Anal..

[8]  P. Hansbo,et al.  A finite element method for domain decomposition with non-matching grids , 2003 .

[9]  Francesco Migliavacca,et al.  Numerical simulation of drug eluting coronary stents: Mechanics, fluid dynamics and drug release , 2009 .

[10]  P. Hansbo,et al.  A unified stabilized method for Stokes' and Darcy's equations , 2007 .

[11]  Jie Shen.  On error estimates of projection methods for Navier-Stokes equations: first-order schemes , 1992 .

[12]  V. Thomée.  Galerkin Finite Element Methods for Parabolic Problems (Springer Series in Computational Mathematics) , 2010 .

[13]  Miguel A. Fernández,et al.  Galerkin Finite Element Methods with Symmetric Pressure Stabilization for the Transient Stokes Equations: Stability and Convergence Analysis , 2008, SIAM J. Numer. Anal..

[14]  Panayot S. Vassilevski,et al.  Interior penalty preconditioners for mixed finite element approximations of elliptic problems , 1996, Math. Comput..

[15]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[16]  Paolo Zunino,et al.  A Finite Element Method Based on Weighted Interior Penalties for Heterogeneous Incompressible Flows , 2009, SIAM J. Numer. Anal..

[17]  Paolo Zunino,et al.  Mathematical models of mass transfer in the vascular walls , 2009 .

[18]  Carlo D'Angelo,et al.  Robust numerical approximation of coupled Stokes' and Darcy's flows applied to vascular hemodynamics and biochemical transport * , 2011 .

[19]  J. Cahouet,et al.  Some fast 3D finite element solvers for the generalized Stokes problem , 1988 .

[20]  Douglas N. Arnold,et al.  Preconditioning in Hdiv and applications , 1997, Math. Comput..