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

Abstract

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.

In silico coronary wave intensity analysis: application of an integrated one-dimensional and poromechanical model of cardiac perfusion

“…2012; Hyde et al. 2014), which can distinguish the pressure and flow across the multiscale hierarchy of distal vessels.…”

Spatio-Temporal Simulation of First Pass Drug Perfusion in the Liver

“…Our combined model was inspired by [24], where a lobule-scale perfusion model in more physical detail and also allowing for deformation of the porous medium is introduced. A model for cardiac perfusion using very similar modeling techniques is presented in [25]. It considers multiple geometric scales, but no draining vascular system and no metabolization.…”

An automatic service for the personalization of ventricular cardiac meshes

“…The methodology proposed here is aimed at mechanical simulations, with application also including continuum approaches for perfusion modelling [38]. The definition of minimum information required to personalize a mesh will depend on the anatomical accuracy needed for the posterior simulation or analysis.…”