High fidelity blood flow in a patient-specific arteriovenous fistula
Output of a valveless Liebau pump with biologically relevant vessel properties and compression frequencies
Rubina Davtyan & Narine A. Sarvazyan
Adaptive constrained constructive optimisation for complex vascularisation processes
Gonzalo Daniel Maso Talou, Soroush Safaei, … Pablo Javier Blanco
Strongly Coupled Morphological Features of Aortic Aneurysms Drive Intraluminal Thrombus
D. Bhagavan, P. Di Achille & J. D. Humphrey
An inverse method for mechanical characterization of heterogeneous diseased arteries using intravascular imaging
Bharath Narayanan, Max L. Olender, … Farhad R. Nezami
The role of vascular complexity on optimal junction exponents
Jonathan Keelan & James P. Hague
5D Flow Tensor MRI to Efficiently Map Reynolds Stresses of Aortic Blood Flow In-Vivo
Jonas Walheim, Hannes Dillinger, … Sebastian Kozerke
Modeling early stage atherosclerosis in a primary human vascular microphysiological system
Xu Zhang, Muath Bishawi, … George A. Truskey
A neurovascular high-frequency optical coherence tomography system enables in situ cerebrovascular volumetric microscopy
Giovanni J. Ughi, Miklos G. Marosfoi, … Ajit S. Puri
Flow driven robotic navigation of microengineered endovascular probes
Lucio Pancaldi, Pietro Dirix, … Mahmut Selman Sakar
J. W. S. McCullough1 &
P. V. Coveney1,2
Scientific Reports volume 11, Article number: 22301 (2021) Cite this article
An arteriovenous fistula, created by artificially connecting segments of a patient's vasculature, is the preferred way to gain access to the bloodstream for kidney dialysis. The increasing power and availability of supercomputing infrastructure means that it is becoming more realistic to use simulations to help identify the best type and location of a fistula for a specific patient. We describe a 3D fistula model that uses the lattice Boltzmann method to simultaneously resolve blood flow in patient-specific arteries and veins. The simulations conducted here, comprising vasculatures of the whole forearm, demonstrate qualified validation against clinical data. Ongoing research to further encompass complex biophysics on realistic time scales will permit the use of human-scale physiological models for basic and clinical medicine.
The development of the virtual human to recreate an individual's specific physiological processes in a digital format is a major goal of computational biomedicine, heralding truly personalised medicine along with 'healthcasts'1,2,3,4,5,6,7. The increasing capability and power of computational software and hardware has allowed significant steps to be made towards this goal1,8. The development of such a 'digital doppelganger' is not a strict progression from low-level (e.g. genomic and DNA studies) to high-level (e.g. organ modelling or epidemiological phenomena) processes but both top-down and bottom-up9. The purpose of the present paper is to report a further advance to these efforts, specifically in the field of blood flow modelling. While 3D flow in large arterial structures has been simulated10,11, we demonstrate how scalable 3D modelling codes can be used to examine macroscopic blood flow in coupled arteries and veins at human scale. We show how the capillary beds linking the two networks in a virtual human can be captured in an efficient manner and applied to human-scale flow problems.
Blood flow is a complex and multiscale process that is fundamental to the operation of the body's physiological processes. From a simulation perspective however, the macroscopic nature of blood may be treated as a continuum fluid that behaves according to the Navier-Stokes equations for fluid flow. Many studies solve these equations for vascular networks through the use of highly simplified 1D representations. The specific properties of an individual's vascular network are used as input for nodal locations within such 1D models12,13,14,15. This approach has advantages relating to the relative ease in parameterisation of the model and the computational effort required to reach a solution. In comparison, 3D models permit a far higher fidelity study of the blood flow by means of an exact representation of the domain of interest and do not rely on major assumptions about flow behaviour (e.g. symmetry) that are made by 1D solvers10,11,16. However, 3D simulations are more demanding to execute. As supercomputers become more powerful, it is more tractable to study 3D flow of large, human-scale, vascular domains10,11,16. Studies of this nature are better equipped to reveal the subtle, patient-specific, phenomena that exist within complex vasculatures. To take advantage of increasing computational performance, a numerical method that demonstrates excellent parallel scaling capability whilst studying complex geometries is essential. In the development of the virtual human, such performance will allow an in silico model to generate necessary data and provide this to a patient's healthcare team in a timely fashion. The lattice Boltzmann method (LBM) used in this work is an ideal candidate for meeting these requirements due to its easily parallelised algorithm and often localised approach to solving the complex geometries or moving boundaries found at multiple scales in blood flow17,18,19.
Irrespective of whether a 1D or 3D model is used to represent the systemic blood vessels, there will remain sections of the vascular network, vital for physiological function, that must be described using sub-grid models. A typical instance of this is the flow through the capillary beds that feed the muscles and skin throughout the body. These vessels are too small to be resolved with typical scanning techniques used to obtain geometries of larger vessels. For instance, typical MRI scans used in hospitals can only resolve features larger than 1 mm20. Alternate scanning techniques such as computed tomography (0.5 mm) and ultrasound (0.15 mm) may achieve improved resolution though at the cost of lower contrast resolution21. The flow behaviour within these unresolved vessels becomes incorporated into the boundary conditions of macroscopic flow models. The development of a self-coupled version of the LBM code used in this study, HemeLB, has been used to show how this framework can simulate flow in human-scale arteries and veins16. Additionally, we provide further details on the modelling methods used in this current work in "Methods".
We provide a contextual setting for this work by examining 3D flow through a personalised representation of the major arteries and veins in the human forearm. For patients suffering from kidney failure who require haemodialysis to remove waste products from the bloodstream, an arteriovenous fistula is a common course of action to facilitate treatment22,23. This procedure often links the cephalic vein in the forearm to the radial artery, although other connections may also be selected. Bypassing the peripheral vasculature increases the blood flow through the vein, in turn increasing the size and strength of the vessel. Once matured, the vein is then used as an access point for blood to be removed for dialysis and returned to the body. The utilisation of in silico models prior to treatment permits surgeons to assess different treatment options prior to undertaking invasive procedures22,23. In this work we examine flow through the left forearm vessels illustrated in Fig. 1. The technique used to capture the 3D vessels was able to resolve vessels to a diameter of approximately 0.3 mm. Vascular structures smaller than this, such as the palmar arch, are represented within the sub-scale material model connecting the simulations of the arteries and veins. In our study, we firstly study flow in unmodified vessels to assess baseline parameters prior to treatment. Following this we demonstrate how the model can be used to assess a candidate fistula location and its impact on vessels throughout the arm.
Human scale vessels in the left forearm. Arteries are coloured red and the veins in blue. The technique used to capture the 3D vessels was able to resolve vessels to a diameter of approximately 0.3 mm. The vessels were truncated at the wrist as shown, the vessels in the hand were represented within the sub-scale material model connecting the simulations of the arteries and veins.
Within the literature, there exists a number of in silico models for studying arteriovenous fistula creation and its application to surgical planning and evaluation24,25,26,27,28,29. Of these, cases that capture large-scale vascular domains are typically 1D studies24,26,28 whilst 3D studies25,27,29 typically only focus on the region of the fistula and do not consider the wider arterial and venous networks throughout the arm. In our model we are able to combine these two features by efficiently simulating 3D resolved flow in detailed networks of large-scale vessels. This is an important step as factors dictating the success of a fistula depend on flow in both the created junction and throughout the remainder of the arm.
In what follows we discuss the formulation of our in silico model. We use this to analyse flow in normal arteries and veins within the left forearm and those that have been modified to incorporate an arteriovenous fistula. We compare our results to those reported in clinical studies and demonstrate excellent agreement.
We applied the sub-scale method discussed in "Sub-scale representation of capillary beds" to a large model of human blood vessels. This domain consists of the left forearm of the full human Yoon-sun model30 (Fig. 1), which is available as part of the Virtual Population (ViP) library31,32. The segmentation of the full human model domain meant that the finest vessels in this domain had a diameter of approximately 0.4 mm. We discretised these selected geometries to a resolution of \(\Delta x\) = 50 μm, resulting in a total of 6,128,855 lattice sites in the arterial geometry and 16,767,772 in the venous geometry. We conducted simulations with a time step of \(\Delta t\) = 3 μs, resulting in a LBM relaxation time of \(\tau\) = 0.5144. This was utilised in a two relaxation time collision kernel with a 'magic' parameter of \(\Lambda\) = 1/1233. These discretisation parameters were selected based on two key considerations. The first, governing \(\Delta x\), was that the finest vessels were sufficiently well resolved to capture the flow within them. The second was that \(\Delta t\) (and by corollary \(\tau\)) permitted the LBM method to stably solve flow within the expected physiological values. The values chosen meet our expectations for these criteria. All test cases were run for a total of 1.1 million time steps. This corresponds to 3.3 s of physical time being simulated and demonstrates the capacity of our model to efficiently represent multiple cardiac cycles in addition to a period of numerical 'warm-up' to prepare flow in the vessels. These were conducted on SuperMUC-NG (at the Leibniz Supercomputing Centre, Germany) and each took approximately 2 h of wall time to run on 50 nodes (2400 cores). In further work, we also examined a similar vascular set discretised to \(\Delta x\) = 25 μm; these contained 64,998,754 lattice sites in the arterial geometry and 166,652,137 in the venous geometry and took approximately 2.5 h to run on 1500 SuperMUC-NG nodes.
In the long term study of patients with arteriovenous fistula, a key observational parameter is how the vessels change shape in response to its insertion. Once this vascular remodelling has stabilised and blood flow is sufficient for dialysis to occur, the fistula is said to have matured. As our model is currently not able to capture this process dynamically, we look to examine the changes to the vascular flow patterns caused by the instantaneous creation of the fistula. These changes may serve as indicators for the success of maturation. A future study conducting such simulations on a population of patients would help to understand this further.
To demonstrate that our model captures the flow behaviour expected in both coupled arteries and veins as well as within an arteriovenous fistula we consider the following cases. The first is a 'healthy' case where we apply a typical flow waveform to the unmodified domains with a single capillary bed region (all venous inlets in the hand region are fed by the two arterial outlets). In the remaining cases we demonstrate how our model can capture an arteriovenous fistula. A common option for creating a fistula is to connect between the radial artery and cephalic vein34. This connection bypasses the capillaries in the hand and serves to increase blood flow in the cephalic vein. In the first of our examples, we modify the arterial and venous geometries to mimic a surgically created end-to-side fistula34 both before and after maturation. We explicitly couple these modifications together within the coupling map.
Inlet flow to the arterial geometry in all cases was provided according to the profile in Fig. 2a. This was generated based on data presented in Mynard and Smolich15 for the right brachial artery.
Healthy vessels
To represent healthy vessels, we leave the entire geometry unmodified. The capillary bed consists of a single bed constructed by both arterial outlets and the 12 venous inlets in the wrist region of the vascular structure.
In Fig. 2 we present the total flow rate entering and exiting the forearm domain being modelled. For the parameters demonstrated here, it can be observed that a steady-state flow response is obtained following the initial warm-up period.
Simulation inlet conditions used in all studies and the observed response in the healthy vessels. (a) details the flow velocity provided to the inlet of the arterial geometry for the forearm flow cases. The first 0.8s of the flow represents an initial warm-up period of flow within the system. (b) indicates the distribution of scaling weights applied to the flow velocity at the inlet of the arterial geometry. (c) compares the total inlet and outlet flow rate through the healthy vessels.
In the following sections we explicitly modify the vascular geometry to incorporate an arteriovenous fistula, firstly immediately after surgical connection and then following a period of maturation. We consider the end-to-side option of joining the cephalic vein to the radial artery. These are illustrated in Fig. 3 schematically alongside their realisation in the modified vascular domains.
An arteriovenous fistula is created by artificially connecting sections of the arterial and venous vasculatures to alter the distribution of blood flow through the vessels. In (a) we provide a schematic map of the geometric modifications made to represent an end-to-side fistula. In this configuration, the cephalic vein is terminated early and attached to the radial artery. In (b) the realisation of these features in pre- and post-maturation geometries is highlighted. In the arterial domain, the bifurcation to the ulnar artery is approximately 150 mm upstream of the anastomosis while the outlet boundary is approximately 20 mm downstream. In the venous domain, the first junction of another vessel (diameter approximately 0.9 mm) with the cephalic vein is approximately 45 mm downstream from the inlet coming from the anastomosis.
End-to-side arteriovenous fistula: pre-maturation
In this setting, we modify the venous geometry such that the cephalic vein is terminated ahead of a junction and mapped in a 1:1 fashion to an opening made in the radial artery. The cephalic vein is kept at the same diameter as the healthy vessels to capture flow pre-maturation. The diameter of the anastomosis connecting the two domains had a diameter of approximately 1.2 mm. This was determined by the diameter of the cephalic vein at the location where we created the fistula.
In Fig. 4, the total flow through the coupled domain again closely matches that provided at the input, indicating the conservation behaviour of the coupling scheme. The particular differences caused by the creation of the fistula can be seen more clearly in Fig. 4b which highlights the transient outflow of the major veins of the domain. Peak flow exiting via the basilic vein has roughly halved whilst the flow exiting the cephalic vein has increased by approximately two-thirds. Under this configuration, there is now almost equal flow exiting the two major vessels.
End-to-side arteriovenous fistula: post-maturation
In this setting, we modify the venous geometry such that the cephalic vein has dilated as a result of maturation. The mapping remains the same as before maturation with the transfer coefficients modified to reflect the change in size of the cephalic vein. In this dilated configuration, the diameter of the cephalic vein has increased by a factor of approximately 2.7. The arterial domain remained unchanged from the pre-maturation case. This dilated configuration was chosen to represent the change to the cephalic vein in a fully mature fistula. In practice, this state may take several weeks to reach. Whilst our model cannot simulate the transition between the pre- and post-maturation states it can compare the flow differences between them.
In Fig. 4, the change in similarities and differences in the flow patterns exiting the venous geometry can immediately be seen. During the periods of peak flow, the redistribution of the total flow compared to the healthy case is quite similar to that seen in the pre-maturation case. During the period of reduced driving flow, in the second half of the heartbeat profile, the flow exiting the cephalic vein reduces much more slowly then that seen in other cases. As the basilic vein profile remains similar to that seen in the pre-maturation case, this indicates that the larger cephalic vein is being supplemented by flow volume that would have exited the non-major veins.
Flow rates through the forearm blood vessels with an arteriovenous fistula in its pre- and post-maturation states. Figures (a) and (c) indicate the total inlet and outlet flow of the system. Figures (b) and (d) compare flow exiting the major veins to that in the 'healthy' vessels.
The use of 3D flow simulations permits far more detailed visualisation of flow than is possible in 1D studies. This can be particularly informative for flow in patient-specific geometries (as was highlighted in "Introduction"). We provide examples of flow through the arteries and veins of the healthy vessels and those with geometry-modified arteriovenous fistula in Fig. 5.
Snapshots of 3D flow after 700,000 time steps of simulation. In all cases the flow is represented by arrows oriented by velocity and coloured by velocity magnitude in m/s. In (a, b), we present flow in the healthy vessels where the preference for arterial flow in the radial artery is clearly evident due to the orientation of the vascular bifurcation. Figures (c, d) represent the flow in the pre- and post-maturation veins respectively. In the venous cases, the highest flow velocities do not occur in the cephalic or basilic veins but in particularly narrow peripheral vessels. The alteration of flow due to the sealed extension of the cephalic vein (noted in Fig. 3b) is obvious, particularly in the pre-maturation case. In a clinical setting, such simulation studies would allow such features to be noted and corrected ahead of any surgery being undertaken. In Figures (e, f) we provide zoomed images of flow features of particular interest. A key advantage of 3D simulation is to be able to generate and examine high fidelity representations of the simulated blood flow. For example, note the drop in flow velocity after the arterial exit in (e).
A particular challenge for the verification and validation of simulated 3D blood flow is the significant variation in the geometric and flow characteristics among humans. In the vessels examined in our study, the inlet diameter of the brachial artery providing flow to the full system was approximately 2.2 mm. In the healthy and pre-maturation configurations, the approximate diameter of the cephalic vein at the entry point of the fistula was 1.4 mm and at the exit plane 2.9 mm. In the matured vessels these measurements were 4.8 and 5.9 mm respectively. When compared to observational literature studies of flow in vessels of the forearm, the unmodified dimensions are particularly narrow but are not unusual35,36. Furthermore, the geometries we are studying were obtained from a female, where narrower vessels again are not unexpected36,37. A further point that may contribute to the narrow vessels in our geometries is that they were ultimately sourced from frozen sections of a cadaver that did not have any fixative or dyeing agents added to it30,38, possibly leading to a narrowing of vessels compared to a living person. It can also be noted from37 that the average flow rate in narrower vessels is also typically less than observed in larger vessels. The inlet velocity profile for our test cases (peak flow speed of 28.6 cm/s and average of 4.2 cm/s) is consistent with literature data39, though this is a measure that can vary widely. The maximum velocity in the brachial artery is approximately 40 cm/s, achieved at a narrow throat shortly after the inlet plane which is also consistent with literature data40. Given the narrow vessels present in our geometries, the brachial flow rate observed (9.8 ml/min) is naturally less then what may be reported in some cases within the literature but consistent with the lower bounds of these when stated standard deviations are taken into account37,41. In the healthy geometry, the flow exiting the cephalic vein averaged 3.0 ml/min whilst again very low, this is not inconsistent with literature data once the geometry and flow characteristics of the test case is taken into account42.
With the presence of a pre-maturation arteriovenous fistula, the flow at the exit plane of the cephalic vein increased by 50% to 4.5 ml/min. In the post-maturation configuration, this flow rate has further increased to 5.6 ml/min, a 90% increase compared to the healthy domain. Before attempting to compare this data with physical arteriovenous fistula, a number of points need to be taken into account. Firstly, these results were obtained with the arterial inlet velocity profile being identical to that used in the healthy case. This means that the increases in flow rates observed are solely due to the changes in vascular geometry caused by the creation of the fistula and the subsequent dilation of the cephalic vein. The 26% increase in outflow due to the dilation of the cephalic vein may be indicative of this enlarged vessel drawing flow from the other veins in the domain. In our study, the arterial domain was not dilated in response to the maturation of the fistula and the connection point on the radial vein remained at the same diameter (the dilation be accounted for in the connecting bed model). This meant that approximately 31% of the total arterial flow was passing through the fistula connection point, corresponding to 39% of the flow passing along the radial artery. With the cross-sectional area of the AVF connection point being only 31% of the radial outlet area, this would also suggest that the presence of the fistula enhances the flow exiting the arterial domain through this pathway.
However, compared to observed data, the flow through the fistula recorded in our setting is far lower than what would be considered successful in a clinical setting. Such examples of 'failed' fistula are not uncommon in practice, both in general and for patients with vessels of similar size to those tested here36. In Wedgwood et al.43, a study of 71 patients who had an arteriovenous fistula created between the radial artery and cephalic vein is presented. The authors present a correlation between the diameter of the vein used to create the fistula and the flow rate through it. As indicated in Fig. 6 the healthy and pre-maturation data is consistent with that the data presented in Wedgwood et al.43, whilst the post-maturation flow rate is significantly less than what would be expected for a vessel of this size. In this plot we have used the average of the inlet and outlet radii of the cephalic vein. It is noted by Van Tricht et al.34 that, for a radiocephalic arteriovenous fistula, a flow rate of at least 500 ml/min is required to be successful (though it has been noted that such a criteria varies globally22). This is around 100 times greater than the cephalic flow observed due to the physiological input to the cases presented here and indicates how significantly the creation of a fistula alters the vasculature of the arm. To achieve such changes, flow feeding the brachial artery would need to similarly increase through a combination of increased mean velocity (such as through a change in velocity profile) or vascular dilation. As an example, changing the inlet profile to a steady maximum input velocity of 11.4 cm/s (corresponds to a mean of 5.5 cm/s) yielded an average flow rate of 12.7 ml/min entering the brachial vein and 6.1 ml/min exiting the cephalic vein in the pre-maturation configuration. Interestingly, the post-maturation cephalic flow was 5.9 ml/min using this input. This indicates that the pulsatility of flow also has an impact on how flow is distributed through the venous domain. Flow through the fistula could be additionally achieved by increased relative area of the connection point of the cephalic vein to the radial artery. As can be seen in Ogawa et al.44, the brachial flow rate in patients with arteriovenous fistula can vary significantly between individual cases, making such behaviour difficult to predict in a modelling study without direct access to clinical data.
Comparison of cephalic vein flow rate calculated by our model and a correlation generated in Wedgwood et al.43. In (a), all data points of Wedgwood et al. have been indicated and the linear correlation extended to encompass the both the observed data and that obtained from our simulations. In (b) we zoom in on the low flow rate data in the vicinity of that observed in our simulation case.
These factors point to some of the limitations of our current modelling approach. The LBM can be sensitive to the choice of discretisation parameters when it comes to the stability and accuracy of simulations and the numerical factors driving these can act in opposition to each other33, particularly in the vicinity of physiologically realistic flows. Significant increase in the flow velocity using the current configuration would lead to inaccuracies and instabilities in the LBM. Whilst these could, in principle, be improved through a reduction in the grid spacing, this would need to be balanced against the challenges of geometry generation, data storage and longer simulation wall times. This should motivate further efforts in HemeLB to incorporate more sophisticated collision kernels that offer options to simulate challenging flow regimes33.
In our study, we have only considered a single set of patient-specific vessels. Increasing the number of domains studied (both pre- and post-maturation) would serve to further streamline our procedure for in silico fistula creation and allow us to develop further confidence in our approach for fistula simulation .
There are a number of further features of our current model that can be improved in future work. In particular, we plan to accurately integrate the behaviour of elastic walls into our model to better represent how they change in response to flow. We will also continue our work in extending the coupling approaches outlined here to enable circulatory flow to be captured within a full-human systemic geometry of arteries and veins that is ultimately coupled to a 3D model of the heart. The framework of the virtual human permits such increases in fidelity or performance to be readily introduced. The completion of these steps will constitute significant further progress along the pathway towards the virtual human. However, the culmination of this endeavour cannot be realised from a simulation perspective alone. The ongoing development of close collaborations between modellers and clinicians will be critical in ensuring that its significant benefits for predictive medicine, patient outcomes and general wellbeing become accessible to all members of society.
We employ the LBM-based solver HemeLB45 to simulate macroscopic blood flow in coupled arterial and venous geometries. Further technical information on the LBM can be found in the many textbooks that cover this topic in detail33,46,47,48,49. HemeLB is an open-source code that has been optimised for the study of flow in sparse geometries such as those typical of vascular networks. Since its inception50, HemeLB has been used to study problems including flow validation in cerebral vessels18, haemodynamic force distributions in retinal vessels51, wall shear stress analysis52,53,54 in cerebral arteries, and the transport of magnetic nanoparticles for drug delivery within the circle of Willis55. These studies demonstrate how HemeLB, and 3D modelling techniques more generally, are able to accurately capture complex physiological flow phenomena in ways that are not possible with 1D approaches and can be difficult or impossible to measure and assess directly within a patient. In this work, we implement velocity inlet boundary conditions and pressure outlet conditions utilising the methods described in previous publications using HemeLB56. We assume that solid walls are rigid with fluid-solid interfaces implemented using the technique of Bouzidi et al.57.
We build on the self-coupled HemeLB presented in our previous work16, which linked two instances of HemeLB in order to capture the respective flow in arterial and venous geometries. The scheme implemented here was designed to streamline the amount of data transferred between the two instances. In Fig. 7, we demonstrate the strong scaling of the coupled model described in McCullough et al.16 as applied to a pair of circle of Willis geometries used in this work for the scaling of a single component version of HemeLB. In this case, we run the simulation for 5050 steps and perform coupling every 500 steps. This demonstrates that the coupled model exhibits strong scaling behaviour up to 96,000 cores on the supercomputer SuperMUC-NG (at the Leibniz Supercomputing Centre, Germany).
Strong scaling performance of the self-coupled HemeLB implementation described in McCullough et al.16 to 96,000 cores of SuperMUC-NG. Ideal and 75% scaling lines are provided as a guide. The coupled simulation transferred data between simulations every 500 steps.
Sub-scale representation of capillary beds
In this work we outline one method designed to represent the capillary beds linking the arterial and venous geometries that cannot be fully resolved for macroscopic simulations. This begins by first mapping the arterial outlets to the venous inlets. The separate collections of these represent individual beds. For any given bed, the sub-scale model needs to be able to accommodate an arbitrary number of arteries and veins being connected. In most regions of the body, there are significantly more veins than arteries of simulatable scale. This means that most beds will have n arteries and m veins connected where typically \(n<m\). Local vascular features, however, may result in alternative configurations. In the current work we identify bed regions and use the sub-scale model to represent the capillaries connecting the arteries and veins located within these regions. Within this model, we use the average velocities from vessels on the opposite geometry to update those on the current set of vessels. The change in velocity here represents the action of the sub-scale capillary network on the flow. Coupling occurs at specified simulation intervals and we alternate the direction of coupling at each instance. As the LBM represents a quasi-incompressible fluid, we use the updated velocities on the arterial geometry to apply a force to the outlets based on the dynamic pressure at that outlet. This represents the volume of fluid within the sub-scale region. This forcing term is applied at the outlet using the approach of Guo et al.58.
The capillary bed approach used in this paper is based on resistor-capacitor models used in some existing models as capillary beds or windkessel boundary conditions. We have taken particular inspiration from the works of Mynard et al.15,59. In our approach we firstly determine the flow weighted velocity, \(v_D\) on the arterial and venous sides of a capillary bed
$$\begin{aligned} v_D = \sum _j\frac{A_jv_j^2}{A_jv_j}. \end{aligned}$$
The flow rate, Q, can then be determined by scaling \(v_D\), by the total area of the appropriate vessels. We then consider the capillary bed as the conceptual model presented in Fig. 8. Here the bed has capacitance, \({\bar{C}}\), and resistance, R. The change in pressure over the bed is computed as \(\Delta P = RQ_B\) where \(Q_B\) is the flow through the bed itself. Taking capacitance into account, the flow of the system can be described as
Conceptual layout of the resistor-capacitor based approach to representing the capillary beds. The arterial, \(Q_A\) and venous, \(Q_V\) flow through the bed is accumulated from and distributed to the connecting vessels by their relative portion of the total vessel area attached to the bed. The capacitance, \({\bar{C}}\), can be viewed as an average of arterial and venous vessel capacitances.
$$\begin{aligned} Q_A - {\bar{C}}\frac{d{\bar{P}}}{dt} = Q_V, \end{aligned}$$
where \({\bar{P}} = 0.5(P_A + P_V)\) is the average pressure of the system. We note that this can be recast as \({\bar{P}} = P_A - 0.5(P_A - P_V) = P_A - 0.5RQ_B\). We assume that the majority of the pressure change arises on the venous side of the system, leading to
$$\begin{aligned} \frac{d{\bar{P}}}{dt} \approx \frac{d}{dt} \left( -0.5RQ_B \right) = -0.5R\frac{dQ_B}{dt}. \end{aligned}$$
This then results in
$$\begin{aligned} Q_A + 0.5 {\bar{C}}R\frac{dQ_B}{dt} = Q_V, \end{aligned}$$
which can be used to update the flow rate on the receiving side of the capillary coupling. This is then distributed to the individual vessels based on the relative area of each. The implementation of this method in HemeLB can be obtained from https://github.com/UCL-CCS/HemePure_SelfCoupled/tree/ResistorCapacitorBeds.
Hunter, P. et al. A vision and strategy for the virtual physiological human: 2012 update. Interface Focus 3(2), 20130004. https://doi.org/10.1098/rsfs.2013.0004 (2013).
Muszkiewicz, A. et al. Variability in cardiac electrophysiology: Using experimentally-calibrated populations of models to move beyond the single virtual physiological human paradigm. Prog. Biophys. Mol. Biol. 120(1), 115–127. https://doi.org/10.1016/j.pbiomolbio.2015.12.002 (2016) (Recent Developments in Biophysics&amp; Molecular Biology of Heart Rhythm).
Biglino, G. et al. Computational modelling for congenital heart disease: How far are we from clinical translation?. Heart 103(2), 98–103. https://doi.org/10.1136/heartjnl-2016-310423 (2017).
Hoekstra, A. G., van Bavel, E., Siebes, M., Gijsen, F. & Geris, L. Virtual physiological human 2016: Translating the virtual physiological human to the clinic. Interface Focus 8(1), 20170067. https://doi.org/10.1098/rsfs.2017.0067 (2018).
Kim, Y. et al. Development of a Korean-specific virtual population for physiologically based pharmacokinetic modelling and simulation. Biopharm. Drug Disposit. 40(3–4), 135–150. https://doi.org/10.1002/bdd.2178 (2019).
Morton, S. E. et al. Predictive virtual patient modelling of mechanical ventilation: Impact of recruitment function. Ann. Biomed. Eng. 47(7), 1626–1641. https://doi.org/10.1007/s10439-019-02253-w (2019).
Hokken, T. W., Ribeiro, J. M., De Jaegere, P. P. & Van Mieghem, N. M. Precision medicine in interventional cardiology. Interv. Cardiac Rev. 15, 3. https://doi.org/10.15420/icr.2019.23 (2020).
Hunter, P. et al. A vision and strategy for the virtual physiological human in 2010 and beyond. Philos. Trans. R. Soc. A 368(1920), 2595–2614. https://doi.org/10.1098/rsta.2010.0048 (2010).
Noble, D. A theory of biological relativity: No privileged level of causation. Interface Focus 2(1), 55–64. https://doi.org/10.1098/rsfs.2011.0067 (2012).
Xiao, N., Humphrey, J. D. & Figueroa, C. A. Multi-scale computational model of three-dimensional hemodynamics within a deformable full-body arterial network. J. Comput. Phys. 244, 22–40. https://doi.org/10.1016/j.jcp.2012.09.016 (2013).
Article ADS MathSciNet PubMed MATH Google Scholar
Randles, A., Draeger, E. W. & Bailey, P. E. Massively parallel simulations of hemodynamics in the primary large arteries of the human vasculature. J. Comput. Sci. 9, 70–75. https://doi.org/10.1016/j.jocs.2015.04.003 (2015).
Sheng, C., Sarwal, S. N., Watts, K. C. & Marble, A. E. Computational simulation of blood flow in human systemic circulation incorporating an external force field. Med. Biol. Eng. Comput. 33(1), 8–17. https://doi.org/10.1007/BF02522938 (1995).
Qureshi, M. U. et al. Numerical simulation of blood flow and pressure drop in the pulmonary arterial and venous circulation. Biomech. Model. Mechanobiol. 13(5), 1137–1154. https://doi.org/10.1007/s10237-014-0563-y (2014).
Müller, L. O. & Toro, E. F. A global multiscale mathematical model for the human circulation with emphasis on the venous system. Int. J. Numer. Methods Biomed. Eng. 30(7), 681–725. https://doi.org/10.1002/cnm.2622 (2014).
Article MathSciNet Google Scholar
Mynard, J. P. & Smolich, J. J. One-dimensional haemodynamic modeling and wave dynamics in the entire adult circulation. Ann. Biomed. Eng. 43(6), 1443–1460 (2015).
McCullough, J. W. S. et al. Towards blood flow in the virtual human: Efficient self-coupling of HemeLB. J. R. Soc. Interface Focus 11, 20190119. https://doi.org/10.1098/rsfs.2019.0119 (2020).
Groen, D. et al. Analysing and modelling the performance of the HemeLB lattice-Boltzmann simulation environment. J. Comput. Sci. 4(5), 412–422. https://doi.org/10.1016/j.jocs.2013.03.002 (2013).
Groen, D. et al. Validation of patient-specific cerebral blood flow simulation using transcranial doppler measurements. Front. Physiol. 9, 721. https://doi.org/10.3389/fphys.2018.00721 (2018).
Latt, J. et al. Palabos: Parallel lattice Boltzmann solver. Comput. Math. Appl.https://doi.org/10.1016/j.camwa.2020.03.022 (2020).
Cleary, J. O. S. H. & Guimar aes, A. R. Magnetic resonance imaging. In Pathobiology of Human Disease (eds McManus, Linda M. & Mitchell, Richard N.) 3987–4004 (Academic Press, 2014). https://doi.org/10.1016/B978-0-12-386456-7.07609-7.
Lin, Eugene & Alessio, Adam. What are the basic concepts of temporal, contrast, and spatial resolution in cardiac CT?. J. Cardiovasc. Comput. Tomogr. 3(6), 403–408. https://doi.org/10.1016/j.jcct.2009.07.003 (2009).
Allon, M. & Robbin, M. L. Increasing arteriovenous fistulas in hemodialysis patients: Problems and solutions. Kidney Int. 62(4), 1109–1124. https://doi.org/10.1111/j.1523-1755.2002.kid551.x (2002).
Tullett, K. Renal Vascular Access: Having a Fistula for Haemodialysis. https://www.uhb.nhs.uk/Downloads/pdf/PiRenalVascularAccess.pdf. Accessed Dec 2019.
Bode, A. S. et al. Patient-specific computational modeling of upper extremity arteriovenous fistula creation: Its feasibility to support clinical decision-making. PLOS ONE 7(4), 1–8, 04. https://doi.org/10.1371/journal.pone.0034491 (2012).
Decorato, I. et al. Numerical simulation of the fluid structure interactions in a compliant patient-specific arteriovenous fistula. Int. J. Numer. Methods Biomed. Eng. 30(2), 143–159. https://doi.org/10.1002/cnm.2595 (2014).
Manini, S. et al. Computational model for simulation of vascular adaptation following vascular access surgery in haemodialysis patients. Comput. Methods Biomech. Biomed. Eng. 17(12), 1358–1367. https://doi.org/10.1080/10255842.2012.745857 (2014).
McGah, P. M., Leotta, D. F., Beach, K. W. & Aliseda, A. Effects of wall distensibility in hemodynamic simulations of an arteriovenous fistula. Biomech. Model. Mechanobiol. 13(3), 679–695 (2014).
Zonnebeld, N., Huberts, W., van Loon, M. M., Delhaas, T. & Tordoir, J. H. M. Preoperative computer simulation for planning of vascular access surgery in hemodialysis patients. J. Vasc. Access 18(1 suppl), S118–S124. https://doi.org/10.5301/jva.5000661 (2017).
Bai, Z. & Zhu, L. Simulation of blood flow past a distal arteriovenous-graft anastomosis at low Reynolds numbers. Phys. Fluids 31(9), 091902. https://doi.org/10.1063/1.5099635 (2019).
Human Model: Yoon-sun. https://itis.swiss/virtual-population/virtual-population/vip3/yoon-sun/. (2019).
Christ, A. et al. The Virtual Family: Development of surface-based anatomical models of two adults and two children for dosimetric simulations. Phys. Med. Biol. 55(2), N23. https://doi.org/10.1088/0031-9155/55/2/N01 (2010).
Article ADS PubMed Google Scholar
Gosselin, M. et al. Development of a new generation of high-resolution anatomical models for medical device evaluation: The Virtual Population 3.0. Phys. Med. Biol. 59(18), 5287. https://doi.org/10.1088/0031-9155/59/18/5287 (2014).
Krüger, T. et al. The Lattice Boltzmann Method: Principles and Practice (Springer, 2017). https://doi.org/10.1007/978-3-319-44649-3.
Van Tricht, I., De Wachter, D., Tordoir, J. & Verdonck, P. Hemodynamics and complications encountered with arteriovenous fistulas and grafts as vascular access for hemodialysis: A review. Ann. Biomed. Eng. 33(9), 1142–1157 (2005).
Lee, C. R. et al. Relation between digital peripheral arterial tonometry and brachial artery ultrasound measures of vascular function in patients with coronary artery disease and in healthy volunteers. Am. J. Cardiol. 109(5), 651–657. https://doi.org/10.1016/j.amjcard.2011.10.023 (2012).
Wilmink, T. & Corte-Real Houlihan, M. Diameter criteria have limited value for prediction of functional dialysis use of arteriovenous fistulas. Eur. J. Vasc. Endovasc. Surg. 56(4), 572–581. https://doi.org/10.1016/j.ejvs.2018.06.066 (2018).
Corretti, M. C., Plotnick, G. D. & Vogel, R. A. The effects of age and gender on brachial artery endothelium-dependent vasoactivity are stimulus-dependent. Clin. Cardiol. 18(8), 471–476. https://doi.org/10.1002/clc.4960180810 (1995).
Park, J. S. et al. Visible Korean human: Improved serially sectioned images of the entire body. IEEE Trans. Med. Imaging 24(3), 352–360. https://doi.org/10.1109/TMI.2004.842454 (2005).
Mitchell, G. F. et al. Local shear stress and brachial artery flow-mediated dilation. Hypertension 44(2), 134–139. https://doi.org/10.1161/01.HYP.0000137305.77635.68 (2004).
Zambanini, A. et al. Wave-energy patterns in carotid, brachial, and radial arteries: A noninvasive approach using wave-intensity analysis. Am. J. Physiol. Heart Circul. Physiol. 289(1), 270–276. https://doi.org/10.1152/ajpheart.00636.2003 (2005).
Chambers, J. C. et al. Association of Maternal endothelial dysfunction with preeclampsia. JAMA 285(12), 1607–1612. https://doi.org/10.1001/jama.285.12.1607 (2001).
Albayrak, R. et al. Hemodynamic changes in the cephalic vein of patients with hemodialysis arteriovenous fistula. J. Clin. Ultrasound 35(3), 133–137. https://doi.org/10.1002/jcu.20307 (2007).
Wedgwood, K. R., Wiggins, P. A. & Guillou, P. J. A prospective study of end-to-side vs. side-to-side arteriovenous fistulas for haemodialysis. Br. J. Surg. 71(8), 640–642. https://doi.org/10.1002/bjs.1800710831 (1984).
Ogawa, T., Matsumura, O., Matsuda, A., Hasegawa, H. & Mitarai, T. Brachial artery blood flow measurement: A simple and noninvasive method to evaluate the need for arteriovenous fistula repair. Dial. Transpl. 40(5), 206–210. https://doi.org/10.1002/dat.20565 (2011).
HemeLB. www.hemelb.org. (2019).
Succi, S. The Lattice Boltzmann Equation for Fluid Dynamics and Beyond (Oxford University Press, 2001).
MATH Google Scholar
Mohamad, A. A. Lattice Boltzmann Method: Fundamentals and Engineering Applications with Computer Codes (Springer, 2011).
Guo, Z. & Shu, C. Lattice Boltzmann Method and Its Applications in Engineering (World Scientific, 2013). https://doi.org/10.1142/8806.
Succi, S. The Lattice Boltzmann Equation: For Complex States of Flowing Matter (Oxford University Press, 2018).
Mazzeo, M. D. & Coveney, P. V. HemeLB: A high performance parallel lattice-Boltzmann code for large scale fluid flow in complex geometries. Comput. Phys. Commun. 178(12), 894–914. https://doi.org/10.1016/j.cpc.2008.02.013 (2008).
Article ADS MathSciNet CAS MATH Google Scholar
Bernabeu, M. O. et al. Computer simulations reveal complex distribution of haemodynamic forces in a mouse retina model of angiogenesis. J. R. Soc. Interface 11(99), 20140543. https://doi.org/10.1098/rsif.2014.0543 (2014).
Bernabeu, M. O. et al. Impact of blood rheology on wall shear stress in a model of the middle cerebral artery. Interface Focus 3(2), 20120094. https://doi.org/10.1098/rsfs.2012.0094 (2013).
Franco, C. A. et al. Dynamic endothelial cell rearrangements drive developmental vessel regression. PLOS Biol. 13(4), 1–19. https://doi.org/10.1371/journal.pbio.1002125 (2015).
Franco, C. A. et al. Non-canonical Wnt signalling modulates the endothelial shear stress flow sensor in vascular remodelling. eLife 5, e07727. https://doi.org/10.7554/eLife.07727 (2016).
Patronis, A. et al. Modeling patient-specific magnetic drug targeting within the intracranial vasculature. Front. Physiol. 9, 331. https://doi.org/10.3389/fphys.2018.00331 (2018).
Nash, R. W. et al. Choice of boundary condition for lattice-Boltzmann simulation of moderate-Reynolds-number flow in complex domains. Phys. Rev. E 89, 023303. https://doi.org/10.1103/PhysRevE.89.023303 (2014).
Bouzidi, M., Firdaouss, M. & Lallemand, P. Momentum transfer of a Boltzmann-lattice fluid with boundaries. Phys. Fluids 13(11), 3452–3459. https://doi.org/10.1063/1.1399290 (2001).
Article ADS CAS MATH Google Scholar
Guo, Z., Zheng, C. & Shi, B. Discrete lattice effects on the forcing term in the lattice Boltzmann method. Phys. Rev. E 65, 046308. https://doi.org/10.1103/PhysRevE.65.046308 (2002).
Mynard, J. P. & Valen-Sendstad, K. A unified method for estimating pressure losses at vascular junctions. Int. J. Numer. Methods Biomed. Eng. 31, 7. https://doi.org/10.1002/cnm.2717 (2015).
The authors acknowledge the Foundation for Research on Information Technologies in Society (IT'IS), particularly B. Lloyd and E. Neufeld, for providing the anatomical geometries used in this project as part of a collaboration within the CompBioMed Centre of Excellence.
We acknowledge funding support from European Commission CompBioMed Centre of Excellence (Grant Nos. 675451 and 823712). Support from the UK Engineering and Physical Sciences Research Council under the project 'UK Consortium on Mesoscale Engineering Sciences (UKCOMES)' (Grant No. EP/R029598/1) is gratefully acknowledged. We acknowledge funding support from MRC for a Medical Bioinformatics Grant (MR/L016311/1), and special funding from the UCL Provost. The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. (www.gauss-centre.eu) for funding this project by providing computing time on the GCS Supercomputer SuperMUC-NG at Leibniz Supercomputing Centre (www.lrz.de).
Centre for Computational Science, Department of Chemistry, University College London, London, UK
J. W. S. McCullough & P. V. Coveney
Informatics Institute, University of Amsterdam, Amsterdam, The Netherlands
P. V. Coveney
J. W. S. McCullough
J.M. developed the coupling model used here, devised, ran and analysed the simulations and drafted the initial version of the manuscript. P.C. supervised the project, provided scientific and technical advice and contributed to writing the manuscript.
Correspondence to P. V. Coveney.
McCullough, J.W.S., Coveney, P.V. High fidelity blood flow in a patient-specific arteriovenous fistula. Sci Rep 11, 22301 (2021). https://doi.org/10.1038/s41598-021-01435-8
About Scientific Reports
Guide to referees
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)