An efficient, localised approach for the simulation of elastic blood vessels using the lattice Boltzmann method

Many numerical studies of blood flow impose a rigid wall assumption due to the simplicity of its implementation compared to a full coupling with a solid mechanics model. In this paper, we present a localised method for incorporating the effects of elastic walls into blood flow simulations using the lattice Boltzmann method implemented by the open-source code HemeLB. We demonstrate that our approach is able to more accurately capture the flow behaviour expected in elastic walled vessels than ones with rigid walls. Furthermore, we show that this can be achieved with no loss of computational performance and remains strongly scalable on high performance computers. We finally illustrate that our approach captures the same trends in wall shear stress distribution as those observed in studies using a rigorous coupling between fluid dynamics and solid mechanics models to solve flow in personalised vascular geometries. These results demonstrate that our model can be used to efficiently and effectively represent flows in elastic blood vessels.