We present a method to efficiently simulate coronary perfusion in subject-specific models of the heart within clinically relevant time frames. Perfusion is modelled as a Darcy porous-media flow, where the permeability tensor is derived from homogenization of an explicit anatomical representation of the vasculature. To account for the disparity in length scales present in the vascular network, in this study, this approach is further refined through the implementation of a multi-compartment medium where each compartment encapsulates the spatial scales in a certain range by using an effective permeability tensor. Neighbouring compartments then communicate through distributed sources and sinks, acting as volume fluxes. Although elegant from a modelling perspective, the full multi-compartment Darcy system is computationally expensive to solve. We therefore enhance computational efficiency of this model by reducing the N-compartment system of Darcy equations to N pressure equations, and N subsequent projection problems to recover the Darcy velocity. The resulting 'reduced' Darcy formulation leads to a dramatic reduction in algebraic-system size and is therefore computationally cheaper to solve than the full multi-compartment Darcy system. A comparison of the reduced and the full formulation in terms of solution time and memory usage clearly highlights the superior performance of the reduced formulation. Moreover, the implementation of flux and, specifically, impermeable boundary conditions on arbitrarily curved boundaries such as epicardium and endocardium is straightforward in contrast to the full Darcy formulation. Finally, to demonstrate the applicability of our methodology to a personalized model and its solvability in clinically relevant time frames, we simulate perfusion in a subject-specific model of the left ventricle.

A method to extract myocardial coronary permeabilities appropriate to parameterise a continuum porous perfusion model using the underlying anatomical vascular network is developed. Canine and porcine whole-heart discrete arterial models were extracted from high-resolution cryomicrotome vessel image stacks. Five parameterisation methods were considered that are primarily distinguished by the level of anatomical data used in the definition of the permeability and pressure-coupling fields. Continuum multi-compartment porous perfusion model pressure results derived using these parameterisation methods were compared quantitatively via a root-mean-square metric to the Poiseuille pressure solved on the discrete arterial vasculature. The use of anatomical detail to parameterise the porous medium significantly improved the continuum pressure results. The majority of this improvement was attributed to the use of anatomically-derived pressure-coupling fields. It was found that the best results were most reliably obtained by using porosity-scaled isotropic permeabilities and anatomically-derived pressure-coupling fields. This paper presents the first continuum perfusion model where all parameters were derived from the underlying anatomical vascular network.