Numerical simulation of blood flow and pressure drop in the pulmonary arterial and venous circulation.

A novel multiscale mathematical and computational model of the pulmonary circulation is presented and used to analyse both arterial and venous pressure and flow. This work is a major advance over previous studies by Olufsen et al. (Ann Biomed Eng 28:1281-1299, 2012) which only considered the arterial circulation. For the first three generations of vessels within the pulmonary circulation, geometry is specified from patient-specific measurements obtained using magnetic resonance imaging (MRI). Blood flow and pressure in the larger arteries and veins are predicted using a nonlinear, cross-sectional-area-averaged system of equations for a Newtonian fluid in an elastic tube. Inflow into the main pulmonary artery is obtained from MRI measurements, while pressure entering the left atrium from the main pulmonary vein is kept constant at the normal mean value of 2 mmHg. Each terminal vessel in the network of 'large' arteries is connected to its corresponding terminal vein via a network of vessels representing the vascular bed of smaller arteries and veins. We develop and implement an algorithm to calculate the admittance of each vascular bed, using bifurcating structured trees and recursion. The structured-tree models take into account the geometry and material properties of the 'smaller' arteries and veins of radii ≥ 50 µm. We study the effects on flow and pressure associated with three classes of pulmonary hypertension expressed via stiffening of larger and smaller vessels, and vascular rarefaction. The results of simulating these pathological conditions are in agreement with clinical observations, showing that the model has potential for assisting with diagnosis and treatment for circulatory diseases within the lung.

Rarefaction and blood pressure in systemic and pulmonary arteries.

The effects of vascular rarefaction (the loss of small arteries) on the circulation of blood are studied using a multiscale mathematical model that can predict blood flow and pressure in the systemic and pulmonary arteries. We augmented a model originally developed for the systemic arteries (Olufsen et al. 1998, 1999, 2000, 2004) to (a) predict flow and pressure in the pulmonary arteries, and (b) predict pressure propagation along the small arteries in the vascular beds. The systemic and pulmonary arteries are modelled as separate, bifurcating trees of compliant and tapering vessels. Each tree is divided into two parts representing the `large' and `small' arteries. Blood flow and pressure in the large arteries are predicted using a nonlinear cross-sectional area-averaged model for a Newtonian fluid in an elastic tube with inflow obtained from magnetic resonance measurements. Each terminal vessel within the network of the large arteries is coupled to a vascular bed of small `resistance' arteries, which are modelled as asymmetric structured trees with specified area and asymmetry ratios between the parent and daughter arteries. For the systemic circulation, each structured tree represents a specific vascular bed corresponding to major organs and limbs. For the pulmonary circulation, there are four vascular beds supplied by the interlobar arteries. This manuscript presents the first theoretical calculations of the propagation of the pressure and flow waves along systemic and pulmonary large and small arteries. Results for all networks were in agreement with published observations. Two studies were done with this model. First, we showed how rarefaction can be modelled by pruning the tree of arteries in the microvascular system. This was done by modulating parameters used for designing the structured trees. Results showed that rarefaction leads to increased mean and decreased pulse pressure in the large arteries. Second, we investigated the impact of decreasing vessel compliance in both large and small arteries. Results showed, that the effects of decreased compliance in the large arteries far outweigh the effects observed when decreasing the compliance of the small arteries. We further showed that a decrease of compliance in the large arteries results in pressure increases consistent with observations of isolated systolic hypertension, as occurs in ageing.