King's College London

Research portal

A computationally efficient framework for the simulation of cardiac perfusion using a multi-compartment Darcy porous-media flow model

Research output: Contribution to journalArticle

Original languageEnglish
Pages (from-to)217-232
Number of pages16
JournalInternational Journal For Numerical Methods In Biomedical Engineering
Issue number2
Early online date18 Oct 2012
Publication statusPublished - Feb 2013

King's Authors


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. Copyright (c) 2012 John Wiley & Sons, Ltd.

View graph of relations

© 2018 King's College London | Strand | London WC2R 2LS | England | United Kingdom | Tel +44 (0)20 7836 5454