A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations

Simulation divergence due to backflow is a common, but not fully addressed, problem in three-dimensional simulations of blood flow in the large vessels. Because backflow is a naturally occurring physiologic phenomenon, careful treatment is necessary to realistically model backflow without artificially altering the local flow dynamics. In this study, we quantitatively compare three available methods for treatment of outlets to prevent backflow divergence in finite element Navier-Stokes solvers. The methods examined are (1) adding a stabilization term to the boundary nodes formulation, (2) constraining the velocity to be normal to the outlet, and (3) using Lagrange multipliers to constrain the velocity profile at all or some of the outlets. A modification to the stabilization method is also discussed. Three model problems, a short and long cylinder with an expansion, a right-angle bend, and a patient-specific aorta model, are used to evaluate and quantitatively compare these methods. Detailed comparisons are made to evaluate robustness, stability characteristics, impact on local and global flow physics, computational cost, implementation effort, and ease-of-use. The results show that the stabilization method offers a promising alternative to previous methods, with reduced effect on both local and global hemodynamics, improved stability, little-to-no increase in computational cost, and elimination of the need for tunable parameters.