A tangential regularization method for backflow stabilization in hemodynamics

In computational simulations of fluid flows, instabilities at the Neumann boundaries may appear during backflow regime. It is widely accepted that this is due to the incoming energy at the boundary from the convection term, which cannot be controlled when the velocity field is unknown. Hence, we propose a stabilized formulation based on a local regularization of the fluid velocity along the tangential directions on the Neumann boundaries. The stabilization term is proportional to the amount of backflow, and does not require any further assumption on the velocity profile. The performance of the method is assessed on a two- and three-dimensional Womersley flows, as well as considering a hemodynamic physiological regime in a patient-specific aortic geometry.

[1]  Volker Gravemeier,et al.  A novel formulation for Neumann inflow boundary conditions in biomechanics , 2012, International journal for numerical methods in biomedical engineering.

[2]  Charles A. Taylor,et al.  Augmented Lagrangian method for constraining the shape of velocity profiles at outlet boundaries for three-dimensional finite element simulations of blood flow , 2009 .

[3]  Thomas J. R. Hughes,et al.  Conservation properties for the Galerkin and stabilised forms of the advection–diffusion and incompressible Navier–Stokes equations , 2005 .

[4]  Alexandre Joel Chorin,et al.  On the Convergence of Discrete Approximations to the Navier-Stokes Equations* , 1989 .

[5]  R. Temam.  Une méthode d'approximation de la solution des équations de Navier-Stokes , 1968 .

[6]  A. Marsden,et al.  A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations , 2011 .

[7]  L. Formaggia,et al.  On the stability of the coupling of 3D and 1D fluid-structure interaction models for blood flow simulations , 2007 .

[8]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[9]  Thomas J. R. Hughes,et al.  Patient-specific isogeometric fluid–structure interaction analysis of thoracic aortic blood flow due to implantation of the Jarvik 2000 left ventricular assist device , 2009 .

[10]  Alejandro F. Frangi,et al.  Efficient computational fluid dynamics mesh generation by image registration , 2007, Medical Image Anal..

[11]  C. Bertoglio.  Forward and Inverse Problems in Fluid-Structure Interaction. Application to Hemodynamics , 2012 .

[12]  H. Si,et al.  Adaptive tetrahedral mesh generation by constrained Delaunay refinement , 2008 .

[13]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: VI. Convergence analysis of the generalized SUPG formulation for linear time-dependent multi-dimensional advective-diffusive systems , 1987 .

[14]  Rolf Rannacher,et al.  ARTIFICIAL BOUNDARIES AND FLUX AND PRESSURE CONDITIONS FOR THE INCOMPRESSIBLE NAVIER–STOKES EQUATIONS , 1996 .

[15]  Alfio Quarteroni,et al.  Numerical Treatment of Defective Boundary Conditions for the Navier-Stokes Equations , 2002, SIAM J. Numer. Anal..

[16]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[17]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .