Computational haemodynamics: geometry and non-newtonian modelling using spectral/hp element methods

The increased availability of large scale computational facilities is enabling the use of computational modelling for the treatment of complex problems of industrial interest. Computational haemodynamics is an area where improved computational power and algorithmic development have permitted the application in biomechanics of techniques originally devised for direct numerical simulations in others areas of scientific interest. In this paper we present a highly accurate algorithm for the numerical solution of the incompressible Navier–Stokes equations which employs a spectral/hp element spatial discretisation and a high-order time splitting. Two important developments discussed in this paper are the coupling of these techniques with unstructured mesh generation, used in many everyday engineering applications, and the introduction of a non-Newtonian viscosity model to account for the shear thinning properties of blood. The proposed method will be applied to the analysis of shear stress distributions in flows of interest in haemodynamics. Received: 5 May 1999 / Accepted: 21 september 1999