VaMpy: A Python Package to Solve 1D Blood Flow Problems

Alexandra K. Diem; Neil W. Bressloff

doi:10.5334/jors.159

Abstract

Finite-differences methods such as the Lax-Wendroff method (LW) are commonly used to solve 1D models of blood flow. These models solve for blood flow and lumen area and are useful in disease research, such as hypertension and atherosclerosis, where flow and pressure are good indicators for the presence of disease. Despite the popularity of the LW method to solve the blood flow equations, no implementation of a LW solver for these equations has been published and made publicly available. This leads to the reimplementation of the same methods within different research groups and makes verification of results more difficult. The Vascular Modelling in Python (VaMpy) toolkit is a Python package that aims to fill this gap. It implements Richtmyer’s two-step Lax-Wendroff scheme to solve 1D model equations of blood flow in arterial trees and aims at facilitating the solution of blood flow problems for various medical applications.

Funding statement: The development of this software was supported by an EPSRC Doctoral Training Centre grant (EP/G03690X/1).

References

Olufsen M S “Numerical Simulation and Experimental Validation of Blood Flow in Arteries with Structured-Tree Outflow Condition” Annals of Biomedical Engineering 2000 28 11 1281 1299 10.1114/1.1326031 URL: http://link.springer.com/article/10.1114/1.1326031
Open DOI Search in Google Scholar Back to article
Sherwin S J “Computational modelling of 1D blood flow with variable mechanical properties and its application to the simulation of wave propagation in the human arterial system” International Journal for Numerical Methods In Fluids 2003 43 6–7 673 700 10.1002/fld.543 URL: http://onlinelibrary.wiley.com/doi/10.1002/fld.543/abstract
Open DOI Search in Google Scholar Back to article
Alastruey J “Modelling the circle of Willis to assess the effects of anatomical variations and occlusions on cerebral flows” Journal of Biomechanics 2007 40 8 1794 1805 10.1016/j.jbiomech.2006.008
Open DOI Search in Google Scholar Back to article
Kolachalama V “Predictive Haemodynamics in a One-Dimensional Carotid Artery Bifurcation. Part I Application to Stent Design” IEEE Transactions on Biomedical Engineering 2007 54 5 802 812 10.1109/TBME.2006.889188 URL: http://ieeexplore.ieee.org/document/4155000/
Open DOI Search in Google Scholar Back to article
Cousins W Gremaud P A “Impedance boundary conditions for general transient hemodynamics” International Journal for Numerical Methods in Biomedical Engineering 2014 30 11 1249 1313 10.1002/cnm.2658 URL: http://onlinelibrary.wiley.com/doi/10.1002/cnm.2658/abstract
Open DOI Search in Google Scholar Back to article
Diem A K “Prediction of Perivascular Drainage of Ab from the Brain Using Computational Modelling: Implications for Alzheimer’s Disease”. PhD thesis 2016 University of Southampton
Search in Google Scholar Back to article
LeVeque R J Numerical Methods for Conservation Laws 1992 2nd Basel, Switzerland Birkhäuser Verlag 122 135 10.1007/978-3-0348-8629-1_12
Open DOI Search in Google Scholar Back to article
Richtmyer R D “A Survey of Difference Methods for Non-Steady Fluid Dynamics” NCAR Technical Notes 1963 63 2 10.5065/D67P8WCQ URL: http://opensky.ucar.edu/islandora/object/technotes:49
Open DOI Search in Google Scholar Back to article
Smith N P Pullan A J Hunter P J “An Anatomically Based Model of Transient Coronary Blood Flow in the Heart” SIAM Journal on Applied Mathematics 2002 62 3 990 1018 0036-1399 10.1137/S0036139999355199 URL: http://epubs.siam.org/doi/abs/10.1137/S0036139999355199
Open DOI Search in Google Scholar Back to article
Azer K Peskin C S “A One-dimensional Model of Blood Flow in Arteries with Friction and Convection Based on the Womersley Velocity Profile” Cardiovascular Engineering 2007 7 2 51 73 https://link.springer.com/article/10.1007/s10558-007-9031-y 10.1007/s10558-007-9031-y
Open DOI Search in Google Scholar Back to article
Itu L M Suciu C “Analysis of outflow boundary condition implementations for 1D blood flow models” Proceedings of the 3rd International Conference on E-Health and Bioengineering 2011 4 24 27 URL: http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=6150403
Search in Google Scholar Back to article
Mynard J P Nithiarasu P “A 1D arterial blood flow model incorporating ventricular pressure, aortic valve and regional coronary flow using the locally conservative Galerkin (LCG) method” Communications in Numerical Methods in Engineering 2008 24 367 417 20407939 10.1002/cnm.1117 URL: http://onlinelibrary.wiley.com/doi/10.1002/cnm.1117/full
Open DOI Search in Google Scholar Back to article
Manini S “pyNS: An Open-Source Framework for 0D Haemodynamic Modelling” Annals of Biomedical Engineering 2015 43 6 1461 1473 15739686 10.1007/s10439-014-1234-y URL: https://link.springer.com/article/10.1007/s10439-014-1234-y
Open DOI Search in Google Scholar Back to article
Devault K “Blood Flow in the Circle of Willis Modeling and Calibration” Multiscale Modeling & Simulation 2008 7 2 888 909 10.1137/07070231X URL: http://epubs.siam.org/doi/abs/10.1137/07070231X
Open DOI Search in Google Scholar Back to article
Olufsen M S “Modeling of the Arterial System with Reference to an Anesthesia Simulator”. PhD Thesis 1998 Denmark University of Roskilde URL: http://rudar.ruc.dk/handle/1800/744
Search in Google Scholar Back to article
Chnafa C “Improved reduced-order modelling of cerebrovascular flow distribution by accounting for arterial bifurcation pressure drops” Journal of Biomechanicsiomechanics 2017 51 83 88 10.1016/j.jbiomech.2016.12.004
Open DOI Search in Google Scholar Back to article
Bakker E N T P “Lymphatic Clearance of the Brain Perivascular, Paravascular and Significance for Neurodegenerative Diseases” Cellular and Molecular Neurobiology 2016 36 2 181 194 10.1007/s10571-015-0273-8 URL: https://link.springer.com/article/10.1007/s10571-015-0273-8
Open DOI Search in Google Scholar Back to article

VaMpy: A Python Package to Solve 1D Blood Flow Problems

Abstract

Paradigm

My account