Pulse wave propagation in a model human arterial network: Assessment of 1-D numerical simulations against in vitro measurements

Manufacturing an artificial arterial tree using 3D printing

Models of the arterial network are useful in studying mechanical cardiac assist devices as well as complex pathological states that are difficult to investigate in-vivo otherwise. Earlier work of artificial arterial tree (AAT) have been constructed to include some of the major arteries and their branches for in-vitro experiments which focused on the aorta, using dipping or painting techniques, which resulted in inaccuracies and inconsistent wall thickness. Therefore, the aim of this work is to use 3D printing for manufacturing AAT based on physiologically correct dimensions of the largest 45 segments of the human arterial tree. A volume ratio mix of silicone rubber (98 %) and a catalyst (2 %) was used to create the walls of the AAT. To validate, the AAT was connected at its inlet to a piston pump that mimicked the heart and capillary tubes at the outlets that mimicked arterial resistances. The capillary tubes were connected to a reservoir that collected the water which was the fluid used in testing the closed-loop hydraulic system. Young's modulus of the AAT walls was determined using tensile testing of different segments of various wall thickness. The developed AAT produced pressure, diameter and flow rate waveforms that are similar to those observed in-vivo. The technique described here is low cost, may be used for producing arterial trees to facilitate testing mechanical cardiac assist devices and studying hemodynamic investigations.

Application and reduction of a nonlinear hyperelastic wall model capturing ex vivo relationships between fluid pressure, area, and wall thickness in normal and hypertensive murine left pulmonary arteries

Pulmonary hypertension is a cardiovascular disorder manifested by elevated mean arterial blood pressure (>20 mmHg) together with vessel wall stiffening and thickening due to alterations in collagen, elastin, and smooth muscle cells. Hypoxia-induced (type 3) pulmonary hypertension can be studied in animals exposed to a low oxygen environment for prolonged time periods leading to biomechanical alterations in vessel wall structure. This study introduces a novel approach to formulating a reduced order nonlinear elastic structural wall model for a large pulmonary artery. The model relating blood pressure and area is calibrated using ex vivo measurements of vessel diameter and wall thickness changes, under controlled pressure conditions, in left pulmonary arteries isolated from control and hypertensive mice. A two-layer, hyperelastic, and anisotropic model incorporating residual stresses is formulated using the Holzapfel-Gasser-Ogden model. Complex relations predicting vessel area and wall thickness with increasing blood pressure are derived and calibrated using the data. Sensitivity analysis, parameter estimation, subset selection, and physical plausibility arguments are used to systematically reduce the 16-parameter model to one in which a much smaller subset of identifiable parameters is estimated via solution of an inverse problem. Our final reduced one layer model includes a single set of three elastic moduli. Estimated ranges of these parameters demonstrate that nonlinear stiffening is dominated by elastin in the control animals and by collagen in the hypertensive animals. The pressure-area relation developed in this novel manner has potential impact on one-dimensional fluids network models of vessel wall remodeling in the presence of cardiovascular disease.

A variable heart rate multi-compartmental coupled model of the cardiovascular and respiratory systems

Current mathematical models of the cardiovascular system that are based on systems of ordinary differential equations are limited in their ability to mimic important features of measured patient data, such as variable heart rates (HR). Such limitations present a significant obstacle in the use of such models for clinical decision-making, as it is the variations in vital signs such as HR and systolic and diastolic blood pressure that are monitored and recorded in typical critical care bedside monitoring systems. In this paper, novel extensions to well-established multi-compartmental models of the cardiovascular and respiratory systems are proposed that permit the simulation of variable HR. Furthermore, a correction to current models is also proposed to stabilize the respiratory behaviour and enable realistic simulation of vital signs over the longer time scales required for clinical management. The results of the extended model developed here show better agreement with measured bio-signals, and these extensions provide an important first step towards estimating model parameters from patient data, using methods such as neural ordinary differential equations. The approach presented is generalizable to many other similar multi-compartmental models of the cardiovascular and respiratory systems.

Bio-inspired microfluidics: A review

Biomicrofluidics, a subdomain of microfluidics, has been inspired by several ideas from nature. However, while the basic inspiration for the same may be drawn from the living world, the translation of all relevant essential functionalities to an artificially engineered framework does not remain trivial. Here, we review the recent progress in bio-inspired microfluidic systems via harnessing the integration of experimental and simulation tools delving into the interface of engineering and biology. Development of "on-chip" technologies as well as their multifarious applications is subsequently discussed, accompanying the relevant advancements in materials and fabrication technology. Pointers toward new directions in research, including an amalgamated fusion of data-driven modeling (such as artificial intelligence and machine learning) and physics-based paradigm, to come up with a human physiological replica on a synthetic bio-chip with due accounting of personalized features, are suggested. These are likely to facilitate physiologically replicating disease modeling on an artificially engineered biochip as well as advance drug development and screening in an expedited route with the minimization of animal and human trials.

An anatomically detailed arterial-venous network model. Cerebral and coronary circulation

In recent years, several works have addressed the problem of modeling blood flow phenomena in veins, as a response to increasing interest in modeling pathological conditions occurring in the venous network and their connection with the rest of the circulatory system. In this context, one-dimensional models have proven to be extremely efficient in delivering predictions in agreement with in-vivo observations. Pursuing the increase of anatomical accuracy and its connection to physiological principles in haemodynamics simulations, the main aim of this work is to describe a novel closed-loop Anatomically-Detailed Arterial-Venous Network (ADAVN) model. An extremely refined description of the arterial network consisting of 2,185 arterial vessels is coupled to a novel venous network featuring high level of anatomical detail in cerebral and coronary vascular territories. The entire venous network comprises 189 venous vessels, 79 of which drain the brain and 14 are coronary veins. Fundamental physiological mechanisms accounting for the interaction of brain blood flow with the cerebro-spinal fluid and of the coronary circulation with the cardiac mechanics are considered. Several issues related to the coupling of arterial and venous vessels at the microcirculation level are discussed in detail. Numerical simulations are compared to patient records published in the literature to show the descriptive capabilities of the model. Furthermore, a local sensitivity analysis is performed, evidencing the high impact of the venous circulation on main cardiovascular variables.

A fast computational model for circulatory dynamics: effects of left ventricle-aorta coupling

The course of diseases such as hypertension, systolic heart failure and heart failure with a preserved ejection fraction is affected by interactions between the left ventricle (LV) and the vasculature. To study these interactions, a computationally efficient, biophysically based mathematical model for the circulatory system is presented. In a four-chamber model of the heart, the LV is represented by a previously described low-order, wall volume-preserving model that includes torsion and base-to-apex and circumferential wall shortening and lengthening, and the other chambers are represented using spherical geometries. Active and passive myocardial mechanics of all four chambers are included. The cardiac model is coupled with a wave propagation model for the aorta and a closed lumped-parameter circulation model. Parameters for the normal heart and aorta are determined by fitting to experimental data. Changes in the timing and magnitude of pulse wave reflections by the aorta are demonstrated with changes in compliance and taper of the aorta as seen in aging (decreased compliance, increased diameter and length), and resulting effects on LV pressure-volume loops and LV fiber stress and sarcomere shortening are predicted. Effects of aging of the aorta combined with reduced LV contractile force (failing heart) are examined. In the failing heart, changes in aortic properties with aging affect stroke volume and sarcomere shortening without appreciable augmentation of aortic pressure, and the reflected pressure wave contributes an increased proportion of aortic pressure.

The contribution of upper and lower body arterial vessels to the aortic root reflections: A one-dimensional computational study

BACKGROUND AND OBJECTIVES

Reflections measured at the aortic root are of physiological and clinical interest and thought to be composed of the superimposed reflections arriving from the upper and lower parts of the circulatory system. However, the specific contribution of each region to the overall reflection measurement has not been thoroughly examined. This study aims to elucidate the relative contribution of reflected waves arising from the upper and lower human body vasculature to those observed at the aortic root.

METHODS

We utilised a one-dimensional (1D) computational model of wave propagation to study reflections in an arterial model that included 37 largest arteries. A narrow Gaussian-shaped pulse was introduced to the arterial model from five distal locations: carotid, brachial, radial, renal, and anterior tibial. The propagation of each pulse towards the ascending aorta was computationally tracked. We calculated the reflected pressure and wave intensity at the ascending aorta in each case. The results are presented as a ratio of the initial pulse.

RESULTS

The findings of this study indicates that pressure pulses originated at the lower body can hardly be observed, while those originated from the upper body account for the largest portion of reflected waves seen at the ascending aorta.

CONCLUSIONS

Our study validates the findings of earlier studies, which demonstrated that human arterial bifurcations have a significantly lower reflection coefficient in the forward direction as compared to the backward direction. The results of this study underscore the need for further in-vivo investigations to provide a deeper understanding of the nature and characteristics of reflections observed in the ascending aorta, which can inform the development of effective strategies for the management of arterial diseases.

A network-based model of dynamic cerebral autoregulation

Cerebrovascular diseases continue to be one of the leading causes of morbidity and mortality in humans. Abnormalities in dynamic cerebral autoregulation (dCA) have been implicated in many of these disease conditions. Accurate models are therefore needed to better understand the complex pathophysiology behind impaired dCA. We thus present here a simple framework for modelling a vessel-driven network model of dCA in the microvasculature, as opposed to the conventional compartmental modelling approach. Network models incorporate the actual connectivity and anatomy of the vasculature, thereby allowing us to include and trace changes in the calibre and morphology of individual vessels, investigate the spatial specificity and heterogeneity of the various control mechanisms to help disentangle their contributions, and link the model parameters to the actual network physiology. The proposed control feedback mechanisms are incorporated at the level of the individual vessel, and the dynamic pressure and flow fields are solved for here within a simple vessel network. In response to an upstream pressure drop, the network is found to be able to recover cerebral blood flow (CBF) while exhibiting the characteristic autoregulatory behaviour in terms of changes in vessel calibre and the biphasic flow response. We assess the feasibility of our formulation in larger networks by comparing the simulation results to those obtained using a one-dimensional (1D) model of CBF applied to the same microvasculature network and find that our model results are in very good agreement with the 1D solution, while significantly reducing the computational cost, thus enabling more detailed models of network behaviour to be adopted in the future. Accurate and computationally feasible models of dCA that are more representative of the vasculature can help increase the translatability of haemodynamic models into the clinical environment, which would help develop more informed treatment guidelines for patients with cerebrovascular diseases.

Arterial wave dynamics preservation upon orthostatic stress: a modelling perspective

Pressure-flow travelling waves are a key topic for understanding arterial haemodynamics. However, wave transmission and reflection processes induced by body posture changes have not been thoroughly explored yet. Current in vivo research has shown that the amount of wave reflection detected at a central level (ascending aorta, aortic arch) decreases during tilting to the upright position, despite the widely proved stiffening of the cardiovascular system. It is known that the arterial system is optimized when in the supine position, i.e. propagation of direct waves is enabled and reflected waves are trapped, protecting the heart; however, it is not known whether this is preserved with postural changes. To shed light on these aspects, we propose a multi-scale modelling approach to inquire into posture-induced arterial wave dynamics elicited by simulated head-up tilting. In spite of remarkable adaptation of the human vasculature following posture changes, our analysis shows that, upon tilting from supine to upright: (i) vessel lumens at arterial bifurcations remain well matched in the forward direction, (ii) wave reflection at central level is reduced due to the backward propagation of weakened pressure waves produced by cerebral autoregulation, and (iii) backward wave trapping is preserved.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

Wave reflection quantification analysis and personalized flow wave estimation based on the central aortic pressure waveform

Pulse wave reflections reflect cardiac afterload and perfusion, which yield valid indicators for monitoring cardiovascular status. Accurate quantification of pressure wave reflections requires the measurement of aortic flow wave. However, direct flow measurement involves extra equipment and well-trained operator. In this study, the personalized aortic flow waveform was estimated from the individual central aortic pressure waveform (CAPW) based on pressure-flow relations. The separated forward and backward pressure waves were used to calculate wave reflection indices such as reflection index (RI) and reflection magnitude (RM), as well as the central aortic pulse transit time (PTT). The effectiveness and feasibility of the method were validated by a set of clinical data (13 participants) and the Nektar1D Pulse Wave Database (4,374 subjects). The performance of the proposed personalized flow waveform method was compared with the traditional triangular flow waveform method and the recently proposed lognormal flow waveform method by statistical analyses. Results show that the root mean square error calculated by the personalized flow waveform approach is smaller than that of the typical triangular and lognormal flow methods, and the correlation coefficient with the measured flow waveform is higher. The estimated personalized flow waveform based on the characteristics of the CAPW can estimate wave reflection indices more accurately than the other two methods. The proposed personalized flow waveform method can be potentially used as a convenient alternative for the measurement of aortic flow waveform.

Effect of lower extremity amputation on cardiovascular hemodynamic environment: An in vitro study

Lower extremity amputation (LEA) was associated with a greater risk of cardiovascular disease, but its hemodynamic mechanisms have not been fully studied. Therefore, to clarify the interrelationship between them, and figure out the potential pathogenesis, the exploration of the hemodynamic environment change of patients after LEA was premeditatedly executed. A near-physiological mock circulatory system (MCS) was employed in the present work to replicate the cardiovascular circulation after LEA in a short time and the unsteady-state numerical simulation was utilized as an auxiliary method to observe the changes of the hemodynamic environment inside the blood vessel. Higher severity of LEA leads to higher peripheral vascular impedance, higher blood pressure, and more obvious redistribution of blood perfusion volume. In addition, higher severity of LEA leads to lower wall shear stress (WSS), higher oscillatory shear index (OSI), and higher relative residence time (RRT) appeared in the infrarenal abdominal aorta and the iliac artery, while these changes are closely related to the higher probability of cardiovascular diseases. Results showed that different degrees of LEA (varying heights, unilateral/bilateral) have diverse effects on the patient's hemodynamic environment. This study explained the potential pathogenesis of cardiovascular diseases after LEA from a hemodynamic perspective and provided a certain reference value for the improvement of the cardiovascular hemodynamic environment and the prevention of cardiovascular diseases in lower extremity amputees.

The effect of cardiac properties on arterial pulse waves: An in-silico study

This study investigated the effects of cardiac properties variability on arterial pulse wave morphology using blood flow modelling and pulse wave analysis. A lumped-parameter model of the left part of the heart was coupled to a one-dimensional model of the arterial network and validated using reference pulse waveforms in turn verified by comparison with in vivo measurements. A sensitivity analysis was performed to assess the effects of variations in cardiac parameters on central and peripheral pulse waveforms. Results showed that left ventricle contractility, stroke volume, cardiac cycle duration, and heart valves impairment are determinants of central waveforms morphology, pulse pressure and its amplification. Contractility of the left atrium has negligible effects on arterial pulse waves. Results also suggested that it might be possible to infer left ventricular dysfunction by analysing the timing of the dicrotic notch and cardiac function by analysing PPG signals. This study has identified cardiac properties that may be extracted from in vivo central and peripheral pulse waves to assess cardiac function.

Automated generation of 0D and 1D reduced-order models of patient-specific blood flow

Three-dimensional (3D) cardiovascular fluid dynamics simulations typically require hours to days of computing time on a high-performance computing cluster. One-dimensional (1D) and lumped-parameter zero-dimensional (0D) models show great promise for accurately predicting blood bulk flow and pressure waveforms with only a fraction of the cost. They can also accelerate uncertainty quantification, optimization, and design parameterization studies. Despite several prior studies generating 1D and 0D models and comparing them to 3D solutions, these were typically limited to either 1D or 0D and a singular category of vascular anatomies. This work proposes a fully automated and openly available framework to generate and simulate 1D and 0D models from 3D patient-specific geometries, automatically detecting vessel junctions and stenosis segments. Our only input is the 3D geometry; we do not use any prior knowledge from 3D simulations. All computational tools presented in this work are implemented in the open-source software platform SimVascular. We demonstrate the reduced-order approximation quality against rigid-wall 3D solutions in a comprehensive comparison with N = 72 publicly available models from various anatomies, vessel types, and disease conditions. Relative average approximation errors of flows and pressures typically ranged from 1% to 10% for both 1D and 0D models, measured at the outlets of terminal vessel branches. In general, 0D model errors were only slightly higher than 1D model errors despite requiring only a third of the 1D runtime. Automatically generated ROMs can significantly speed up model development and shift the computational load from high-performance machines to personal computers.

Uncertainty quantification in cerebral circulation simulations focusing on the collateral flow: Surrogate model approach with machine learning

Collateral circulation in the circle of Willis (CoW), closely associated with disease mechanisms and treatment outcomes, can be effectively investigated using one-dimensional–zero-dimensional hemodynamic simulations. As the entire cardiovascular system is considered in the simulation, it captures the systemic effects of local arterial changes, thus reproducing collateral circulation that reflects biological phenomena. The simulation facilitates rapid assessment of clinically relevant hemodynamic quantities under patient-specific conditions by incorporating clinical data. During patient-specific simulations, the impact of clinical data uncertainty on the simulated quantities should be quantified to obtain reliable results. However, as uncertainty quantification (UQ) is time-consuming and computationally expensive, its implementation in time-sensitive clinical applications is considered impractical. Therefore, we constructed a surrogate model based on machine learning using simulation data. The model accurately predicts the flow rate and pressure in the CoW in a few milliseconds. This reduced computation time enables the UQ execution with 100 000 predictions in a few minutes on a single CPU core and in less than a minute on a GPU. We performed UQ to predict the risk of cerebral hyperperfusion (CH), a life-threatening condition that can occur after carotid artery stenosis surgery if collateral circulation fails to function appropriately. We predicted the statistics of the postoperative flow rate increase in the CoW, which is a measure of CH, considering the uncertainties of arterial diameters, stenosis parameters, and flow rates measured using the patients’ clinical data. A sensitivity analysis was performed to clarify the impact of each uncertain parameter on the flow rate increase. Results indicated that CH occurred when two conditions were satisfied simultaneously: severe stenosis and when arteries of small diameter serve as the collateral pathway to the cerebral artery on the stenosis side. These findings elucidate the biological aspects of cerebral circulation in terms of the relationship between collateral flow and CH.

Simulation of flow in an artery under pathological hemodynamic conditions: The use of a diagnostic disease descriptor

A numerical model for simulating and predicting blood flow dynamics in diseased arterial vessels has been developed. The time-dependent one-dimensional hyperbolic system of quasilinear partial differential equations which incorporates a diagnostic disease descriptor (kD) was used to simulate transient flow distribution for idealized healthy and diseased states. Blood flow simulations in the iliac arteries over about 125% of a cardiac cycle were generated and calibrated using the kD values from 0 to 3 representing hypothetical diseased states. Early results indicate that disease conditions induce abnormal flow in the artery, generating disorder and increased amplitude of blood pressure, flow and distensibility with increasing numerical values of the disease factor kD. More so, the prospective use of the kD-approach with documentation of in vivo adverse flow visualizations for diagnostic purposes was decisively discussed.

Wave Intensity Analysis Combined With Machine Learning can Detect Impaired Stroke Volume in Simulations of Heart Failure

Heart failure is treatable, but in the United Kingdom, the 1-, 5- and 10-year mortality rates are 24.1, 54.5 and 75.5%, respectively. The poor prognosis reflects, in part, the lack of specific, simple and affordable diagnostic techniques; the disease is often advanced by the time a diagnosis is made. Previous studies have demonstrated that certain metrics derived from pressure-velocity-based wave intensity analysis are significantly altered in the presence of impaired heart performance when averaged over groups, but to date, no study has examined the diagnostic potential of wave intensity on an individual basis, and, additionally, the pressure waveform can only be obtained accurately using invasive methods, which has inhibited clinical adoption. Here, we investigate whether a new form of wave intensity based on noninvasive measurements of arterial diameter and velocity can detect impaired heart performance in an individual. To do so, we have generated a virtual population of two-thousand elderly subjects, modelling half as healthy controls and half with an impaired stroke volume. All metrics derived from the diameter-velocity-based wave intensity waveforms in the carotid, brachial and radial arteries showed significant crossover between groups-no one metric in any artery could reliably indicate whether a subject's stroke volume was normal or impaired. However, after applying machine learning to the metrics, we found that a support vector classifier could simultaneously achieve up to 99% recall and 95% precision. We conclude that noninvasive wave intensity analysis has significant potential to improve heart failure screening and diagnosis.

A method of parameter estimation for cardiovascular hemodynamics based on deep learning and its application to personalize a reduced-order model

Precise model personalization is a key step towards the application of cardiovascular physical models. In this manuscript, we propose to use deep learning (DL) to solve the parameter estimation problem in cardiovascular hemodynamics. Based on the convolutional neural network (CNN) and fully connected neural network (FCNN), a multi-input deep neural network (DNN) model is developed to map the nonlinear relationship between measurements and the parameters to be estimated. In this model, two separate network structures are designed to extract the features of two types of measurement data, including pressure waveforms and a vector composed of heart rate (HR) and pulse transit time (PTT), and a shared structure is used to extract their combined dependencies on the parameters. Besides, we try to use the transfer learning (TL) technology to further strengthen the personalized characteristics of a trained-well network. For assessing the proposed method, we conducted the parameter estimation using synthetic data and in vitro data respectively, and in the test with synthetic data, we evaluated the performance of the TL algorithm through two individuals with different characteristics. A series of estimation results show that the estimated parameters are in good agreement with the true values. Furthermore, it is also found that the estimation accuracy can be significantly improved by a multicycle combination strategy. Therefore, we think that the proposed method has the potential to be used for parameter estimation in cardiovascular hemodynamics, which can provide an immediate, accurate, and sustainable personalization process, and deserves more attention in the future.

Pulsatility damping in the microcirculation: Basic pattern and modulating factors

Blood flow pulsatility is an important determinant of macro- and microvascular physiology. Pulsatility is damped largely in the microcirculation, but the characteristics of this damping and the factors that regulate it have not been fully elucidated yet. Applying computational approaches to real microvascular network geometry, we examined the pattern of pulsatility damping and the role of potential damping factors, including pulse frequency, vascular viscous resistance, vascular compliance, viscoelastic behavior of the vessel wall, and wave propagation and reflection. To this end, three full rat mesenteric vascular networks were reconstructed from intravital microscopic recordings, a one-dimensional (1D) model was used to reproduce pulsatile properties within the network, and potential damping factors were examined by sensitivity analysis. Results demonstrate that blood flow pulsatility is predominantly damped at the arteriolar side and remains at a low level at the venular side. Damping was sensitive to pulse frequency, vascular viscous resistance and vascular compliance, whereas viscoelasticity of the vessel wall or wave propagation and reflection contributed little to pulsatility damping. The present results contribute to our understanding of mechanical forces and their regulation in the microcirculation.

Impact of tapering of arterial vessels on blood pressure, pulse wave velocity, and wave intensity analysis using one-dimensional computational model

The angle of arterial tapering increases with ageing, and the geometrical changes of the aorta may cause an increase in central arterial pressure and stiffness. The impact of tapering has been primarily studied using frequency-domain transmission line theories. In this work, we revisit the problem of tapering and investigate its effect on blood pressure and pulse wave velocity (PWV) using a time-domain analysis with a 1D computational model. First, tapering is modelled as a stepwise reduction in diameter and compared with results from a continuously tapered segment. Next, we studied wave reflections in a combination of stepwise diameter reduction of straight vessels and bifurcations, then repeated the experiments with decreasing the length to physiological values. As the model's segments became shorter in length, wave reflections and re-reflections resulted in waves overlapping in time. We extended our work by examining the effect of increasing the tapering angle on blood pressure and wave intensity in physiological models: a model of the thoracic aorta and a model of upper thoracic and descending aorta connected to the iliac bifurcation. Vessels tapering inherently changed the ratio between the inlet and outlet cross-sectional areas, increasing the vessel resistance and reducing the compliance compared with non-tapered vessels. These variables influence peak and pulse pressure. In addition, it is well established that pulse wave velocity increases in an ageing arterial tree. This work provides confirmation that tapering induces reflections and offers an additional explanation to the observation of increased peak pressure and decreased diastolic pressure distally in the arterial tree.

A physiologically realistic virtual patient database for the study of arterial haemodynamics

This study creates a physiologically realistic virtual patient database (VPD), representing the human arterial system, for the primary purpose of studying the effects of arterial disease on haemodynamics. A low dimensional representation of an anatomically detailed arterial network is outlined, and a physiologically realistic posterior distribution for its parameters constructed through the use of a Bayesian approach. This approach combines both physiological/geometrical constraints and the available measurements reported in the literature. A key contribution of this work is to present a framework for including all such available information for the creation of virtual patients (VPs). The Markov Chain Monte Carlo (MCMC) method is used to sample random VPs from this posterior distribution, and the pressure and flow-rate profiles associated with each VP computed through a physics based model of pulse wave propagation. This combination of the arterial network parameters (representing a virtual patient) and the haemodynamics waveforms of pressure and flow-rates at various locations (representing functional response and potential measurements that can be acquired in the virtual patient) makes up the VPD. While 75,000 VPs are sampled from the posterior distribution, 10,000 are discarded as the initial burn-in period of the MCMC sampler. A further 12,857 VPs are subsequently removed due to the presence of negative average flow-rate, reducing the VPD to 52,143. Due to undesirable behaviour observed in some VPs-asymmetric under- and over-damped pressure and flow-rate profiles in left and right sides of the arterial system-a filter is proposed to remove VPs showing such behaviour. Post application of the filter, the VPD has 28,868 subjects. It is shown that the methodology is appropriate by comparing the VPD statistics to those reported in literature across real populations. Generally, a good agreement between the two is found while respecting physiological/geometrical constraints.

A proof of concept study for machine learning application to stenosis detection

This proof of concept (PoC) assesses the ability of machine learning (ML) classifiers to predict the presence of a stenosis in a three vessel arterial system consisting of the abdominal aorta bifurcating into the two common iliacs. A virtual patient database (VPD) is created using one-dimensional pulse wave propagation model of haemodynamics. Four different machine learning (ML) methods are used to train and test a series of classifiers—both binary and multiclass—to distinguish between healthy and unhealthy virtual patients (VPs) using different combinations of pressure and flow-rate measurements. It is found that the ML classifiers achieve specificities larger than 80% and sensitivities ranging from 50 to 75%. The most balanced classifier also achieves an area under the receiver operative characteristic curve of 0.75, outperforming approximately 20 methods used in clinical practice, and thus placing the method as moderately accurate. Other important observations from this study are that (i) few measurements can provide similar classification accuracies compared to the case when more/all the measurements are used; (ii) some measurements are more informative than others for classification; and (iii) a modification of standard methods can result in detection of not only the presence of stenosis, but also the stenosed vessel.

Graphical Abstract

Machine learning for detection of stenoses and aneurysms: application in a physiologically realistic virtual patient database

This study presents an application of machine learning (ML) methods for detecting the presence of stenoses and aneurysms in the human arterial system. Four major forms of arterial disease-carotid artery stenosis (CAS), subclavian artery stenosis (SAS), peripheral arterial disease (PAD), and abdominal aortic aneurysms (AAA)-are considered. The ML methods are trained and tested on a physiologically realistic virtual patient database (VPD) containing 28,868 healthy subjects, adapted from the authors previous work and augmented to include disease. It is found that the tree-based methods of Random Forest and Gradient Boosting outperform other approaches. The performance of ML methods is quantified through the [Formula: see text] score and computation of sensitivities and specificities. When using six haemodynamic measurements (pressure in the common carotid, brachial, and radial arteries; and flow-rate in the common carotid, brachial, and femoral arteries), it is found that maximum [Formula: see text] scores larger than 0.9 are achieved for CAS and PAD, larger than 0.85 for SAS, and larger than 0.98 for both low- and high-severity AAAs. Corresponding sensitivities and specificities are larger than 90% for CAS and PAD, larger than 85% for SAS, and larger than 98% for both low- and high-severity AAAs. When reducing the number of measurements, performance is degraded by less than 5% when three measurements are used, and less than 10% when only two measurements are used for classification. For AAA, it is shown that [Formula: see text] scores larger than 0.85 and corresponding sensitivities and specificities larger than 85% are achievable when using only a single measurement. The results are encouraging to pursue AAA monitoring and screening through wearable devices which can reliably measure pressure or flow-rates.

Multiscale Coupling of One-dimensional Vascular Models and Elastic Tissues

We present a computational multiscale model for the efficient simulation of vascularized tissues, composed of an elastic three-dimensional matrix and a vascular network. The effect of blood vessel pressure on the elastic tissue is surrogated via hyper-singular forcing terms in the elasticity equations, which depend on the fluid pressure. In turn, the blood flow in vessels is treated as a one-dimensional network. Intravascular pressure and velocity are simulated using a high-order finite volume scheme, while the elasticity equations for the tissue are solved using a finite element method. This work addresses the feasibility and the potential of the proposed coupled multiscale model. In particular, we assess whether the multiscale model is able to reproduce the tissue response at the effective scale (of the order of millimeters) while modeling the vasculature at the microscale. We validate the multiscale method against a full scale (three-dimensional) model, where the fluid/tissue interface is fully discretized and treated as a Neumann boundary for the elasticity equation. Next, we present simulation results obtained with the proposed approach in a realistic scenario, demonstrating that the method can robustly and efficiently handle the one-way coupling between complex fluid microstructures and the elastic matrix.

Influence of cross-sectional velocity profile on flow characteristics of arterial wall modeled as elastic and viscoelastic material

In the present work, blood flow behavior in a single artery and in arterial network is studied using time domain based one-dimensional wave propagation model retaining the nonlinear convective force. 1-D Navier-Stokes equation is used to model the flow behavior of the blood, using three unknown parameters: flow rate (q), cross-sectional area of artery (A) and pressure (p) based formulation. Three different approximate velocity profile functions across the cross-section namely modified flat, parabolic and the one proposed by Bessems are used to calculate the nonlinear convective force and the frictional force. Two different constitutive models, linear elastic model and standard linear solid model (Zener model) are used to model the arterial wall mechanical behavior. The system of partial differential equations is discretized using finite element and Crank Nicolson methods in space and time domains, respectively. Based on the formulation, an in house finite element code is developed to simulate flow behavior in both a single artery as well as in arterial network consisting of 20 small and large size arteries. Simulations are performed by enforcing a flow rate at the inlet and Windkessel model at the outlet. The results for elastic arterial wall model are found to be in good agreement with the results available in the literature. The flow rate/pressure predictions using different velocity profile functions are found to be nearly the same, however the Bessems velocity profile predicts more closer to 3D results compared to modified flat and parabolic profiles. Whereas, significant difference is found in the results predicted using elastic and viscoelastic artery wall models.

Arterial pulse wave propagation across stenoses and aneurysms: assessment of one-dimensional simulations against three-dimensional simulations and in vitro measurements

One-dimensional (1-D) arterial blood flow modelling was tested in a series of idealized vascular geometries representing the abdominal aorta, common carotid and iliac arteries with different sizes of stenoses and/or aneurysms. Three-dimensional (3-D) modelling and in vitro measurements were used as ground truth to assess the accuracy of 1-D model pressure and flow waves. The 1-D and 3-D formulations shared identical boundary conditions and had equivalent vascular geometries and material properties. The parameters of an experimental set-up of the abdominal aorta for different aneurysm sizes were matched in corresponding 1-D models. Results show the ability of 1-D modelling to capture the main features of pressure and flow waves, pressure drop across the stenoses and energy dissipation across aneurysms observed in the 3-D and experimental models. Under physiological Reynolds numbers (Re), root mean square errors were smaller than 5.4% for pressure and 7.3% for the flow, for stenosis and aneurysm sizes of up to 85% and 400%, respectively. Relative errors increased with the increasing stenosis and aneurysm size, aneurysm length and Re, and decreasing stenosis length. All data generated in this study are freely available and provide a valuable resource for future research.

A porous circulation model of the human brain for in silico clinical trials in ischaemic stroke

The advancement of ischaemic stroke treatment relies on resource-intensive experiments and clinical trials. In order to improve ischaemic stroke treatments, such as thrombolysis and thrombectomy, we target the development of computational tools for in silico trials which can partially replace these animal and human experiments with fast simulations. This study proposes a model that will serve as part of a predictive unit within an in silico clinical trial estimating patient outcome as a function of treatment. In particular, the present work aims at the development and evaluation of an organ-scale microcirculation model of the human brain for perfusion prediction. The model relies on a three-compartment porous continuum approach. Firstly, a fast and robust method is established to compute the anisotropic permeability tensors representing arterioles and venules. Secondly, vessel encoded arterial spin labelling magnetic resonance imaging and clustering are employed to create an anatomically accurate mapping between the microcirculation and large arteries by identifying superficial perfusion territories. Thirdly, the parameter space of the problem is reduced by analysing the governing equations and experimental data. Fourthly, a parameter optimization is conducted. Finally, simulations are performed with the tuned model to obtain perfusion maps corresponding to an open and an occluded (ischaemic stroke) scenario. The perfusion map in the occluded vessel scenario shows promising qualitative agreement with computed tomography images of a patient with ischaemic stroke caused by large vessel occlusion. The results highlight that in the case of vessel occlusion (i) identifying perfusion territories is essential to capture the location and extent of underperfused regions and (ii) anisotropic permeability tensors are required to give quantitatively realistic estimation of perfusion change. In the future, the model will be thoroughly validated against experiments.

On the poro-elastic models for microvascular blood flow resistance: An in vitro validation

Nowadays, adequate and accurate representation of the microvascular flow resistance constitutes one of the major challenges in computational haemodynamic studies. In this work, a theoretical, porous media framework, ultimately designed for representing downstream resistance, is presented and compared against an in vitro experimental results. The resistor consists of a poro-elastic tube, with either a constant or variable porosity profile in space. The underlying physics, characterizing the fluid flow through the porous media, is analysed by considering flow variables at different network locations. Backward reflections, originated in the reservoir of the in vitro model, are accounted for through a reflection coefficient imposed as an outflow network condition. The simulation results are in good agreement with the measurements for both the homogenous and heterogeneous porosity conditions. In addition, the comparison allows identification of the range of values representing experimental reservoir reflection coefficients. The pressure drops across the heterogeneous porous media increases with respect to the simpler configuration, whilst flow remains almost unchanged. The effect of some fluid network features, such as tube Young's modulus and fluid viscosity, on the theoretical results is also elucidated, providing a reference for the invitro and insilico simulation of different microvascular conditions.

Comparison of arterial wave intensity analysis by pressure-velocity and diameter-velocity methods in a virtual population of adult subjects

Pressure-velocity-based analysis of arterial wave intensity gives clinically relevant information about the performance of the heart and vessels, but its utility is limited because accurate pressure measurements can only be obtained invasively. Diameter-velocity-based wave intensity can be obtained noninvasively using ultrasound; however, due to the nonlinear relationship between blood pressure and arterial diameter, the two wave intensities might give disparate clinical indications. To test the magnitude of the disagreement, we have generated an age-stratified virtual population to investigate how the two dominant nonlinearities viscoelasticity and strain-stiffening cause the two formulations to differ. We found strong agreement between the pressure-velocity and diameter-velocity methods, particularly for the systolic wave energy, the ratio between systolic and diastolic wave heights, and older subjects. The results are promising regarding the introduction of noninvasive wave intensities in the clinic.

Conceptual Framework Development for a Double-Walled Aortic Stent-Graft to Manage Blood Pressure

A double-walled stent-graft (DWSG) design with a compressible gas layer was conceived with the goal of treating hypertension in patients receiving an aortic stent-graft. Early prototypes were developed to evaluate the design concept through static measurements from a finite element (FE) model and quasi-static inflation experiments, and through dynamic measurements from an in vitro flow loop and the three-element Windkessel model. The amount of gas in the gas layer and the properties of the flexible inner wall were the primary variables evaluated in this study. Properties of the inner wall had minimal effect on DWSG behavior, but increased gas charge led to increased fluid capacitance and larger reduction in peak and pulse pressures. In the flow loop, placement of the DWSG decreased pulse pressure by over 20% compared to a rigid stent-graft. Capacitance measurements were consistent across all methods, with the maximum capacitance estimated at 0.07 mL/mmHg for the largest gas charge in the 15 cm long prototype. Windkessel model predictions for in vivo performance of a DWSG placed in the aorta of a hypertensive patient showed pulse pressure reduction of 14% compared to a rigid stent-graft case, but pressures never returned to unstented values. These results indicate that the DWSG design has potential to be developed into a new treatment for hypertensive patients requiring an aortic intervention.

The Effects of Cerebral Vasospasm on Cerebral Blood Flow and the Effects of Induced Hypertension: A Mathematical Modelling Study

Background

Induced hypertension has been used to promote cerebral blood flow under vasospastic conditions although there is no randomised clinical trial to support its use. We sought to mathematically model the effects of vasospasm on the cerebral blood flow and the effects of induced hypertension.

Methods

The Anatomically Detailed Arterial Network (ADAN) model is employed as the anatomical substrate in which the cerebral blood flow is simulated as part of the simulation of the whole body arterial circulation. The pressure drop across the spastic vessel is modelled by inserting a specific constriction model within the corresponding vessel in the ADAN model. We altered the degree of vasospasm, the length of the vasospastic segment, the location of the vasospasm, the pressure (baseline mean arterial pressure [MAP] 90 mm Hg, hypertension MAP 120 mm Hg, hypotension), and the presence of collateral supply.

Results

Larger decreases in cerebral flow were seen for diffuse spasm and more severe vasospasm. The presence of collateral supply could maintain cerebral blood flow, but only if the vasospasm did not occur distal to the collateral. Induced hypertension caused an increase in blood flow in all scenarios, but did not normalise blood flow even in the presence of moderate vasospasm (30%). Hypertension in the presence of a complete circle of Willis had a marginally greater effect on the blood flow, but did not normalise flow.

Conclusion

Under vasospastic condition, cerebral blood flow varies considerably. Hypertension can raise the blood flow, but it is unable to restore cerebral blood flow to baseline.

Study protocol of validating a numerical model to assess the blood flow in the circle of Willis

INTRODUCTION

We developed a zero-dimensional (0D) model to assess the patient-specific haemodynamics in the circle of Willis (CoW). Similar numerical models for simulating the cerebral blood flow (CBF) had only been validated qualitatively in healthy volunteers by magnetic resonance (MR) angiography and transcranial Doppler (TCD). This study aims to validate whether a numerical model can simulate patient-specific blood flow in the CoW under pathological conditions.

METHODS AND ANALYSIS

This study is a diagnostic accuracy study. We aim to collect data from a previously performed prospective study that involved patients with aneurysmal subarachnoid haemorrhage (aSAH) receiving both TCD and brain Computerd Tomography angiography (CTA) at the same day. The cerebral flow velocities are calculated by the 0D model, based on the vessel diameters measured on the CTA of each patient. In this study, TCD is considered the gold standard for measuring flow velocity in the CoW. The agreement will be analysed using Pearson correlation coefficients.

ETHICS AND DISSEMINATION

This study protocol has been approved by the Medical Ethics Review Board of the University Medical Center Groningen: METc2019/103. The results will be submitted to an international scientific journal for peer-reviewed publication.

TRIAL REGISTRATION NUMBER

NL8114.

Dynamic Effects of Aortic Arch Stiffening on Pulsatile Energy Transmission to Cerebral Vasculature as A Determinant of Brain-Heart Coupling

Aortic stiffness increases with age and is a robust predictor of brain pathology including Alzheimer’s and other dementias. Aging causes disproportionate stiffening of the aorta compared with the carotid arteries, reducing protective impedance mismatches at their interface and affecting transmission of destructive pulsatile energy to the cerebral circulation. Recent clinical studies have measured regional stiffness within the aortic arch using pulse wave velocity (PWV) and have found a stronger association with cerebrovascular events than global stiffness measures. However, effects of aortic arch PWV on the transmission of harmful excessive pulsatile energy to the brain is not well-understood. In this study, we use an energy-based analysis of hemodynamic waves to quantify the effect of aortic arch stiffening on transmitted pulsatility to cerebral vasculature, employing a computational approach using a one-dimensional model of the human vascular network. Results show there exists an optimum wave condition—occurring near normal human heart rates—that minimizes pulsatile energy transmission to the brain. This indicates the important role of aortic arch biomechanics on heart-brain coupling. Our results also suggest that energy-based indices of pulsatility combining pressure and flow data are more sensitive to increased stiffness than using flow or pressure pulsatility indices in isolation.

Deep learning for prediction of cardiac indices from photoplethysmographic waveform: A virtual database approach

Deep learning methods combined with large datasets have recently shown significant progress in solving several medical tasks. However, collecting and annotating large datasets can be a very cumbersome and expensive task. We tackle these problems with a virtual database approach where training data is generated using computer simulations of related phenomena. Specifically, we concentrate on the following problem: can cardiovascular indices such as aortic elasticity, diastolic and systolic blood pressures, and blood flow from heart be predicted continuously using wearable photoplethysmographic sensors? We simulate the blood flow using a haemodynamic model consisting of the entire human circulation. Repeated evaluation of the simulator allows us to create a database of "virtual subjects" with size that is only limited by available computational resources. Using this database, we train neural networks to predict the cardiac indices from photoplethysmographic signal waveform. We consider two approaches: neural networks based on predefined input features and deep convolutional neural networks taking waveform directly as the input. The performance of the methods is demonstrated using numerical examples, thus carrying out a preliminary assessment of the approaches. The results show improvements in accuracy compared with the previous methods. The improvements are especially significant with indices related to aortic elasticity and maximum blood flow. The proposed approach would provide new means to measure cardiovascular health continuously, for example, with a simple wrist device.

Subject-specific simulation for non-invasive assessment of aortic coarctation: Towards a translational approach

We present a multi-scale CFD-based study conducted in a cohort of 11 patients with coarctation of the aorta (CoA). The study explores the potential for implementation of a workflow using non-invasive routinely collected medical imaging data and clinical measurements to provide a more detailed insight into local aortic haemodynamics in order to support clinical decision making. Our approach is multi-scale, using a reduced-order model (1D/0D) and an optimization process for the personalization of patient-specific boundary conditions and aortic vessel wall parameters from non-invasive measurements, to inform a more complex model (3D/0D) representing 3D aortic patient-specific anatomy. The reliability of the modelling approach is investigated by comparing 3D/0D model pressure drop estimation with measured peak gradients recorded during diagnostic cardiac catheterization and 2D PC-MRI flow rate measurements in the descending aorta. The current study demonstrated that the proposed approach requires low levels of user interaction, making it suitable for the clinical setting. The agreement between computed blood pressure drop and catheter measurements is 10  ±  8 mmHg at the coarctation site. The comparison between CFD derived and catheter measured pressure gradients indicated that the model has to be improved, suggesting the use of time varying pressure waveforms to further optimize the tuning process and modelling assumptions.

Quantification of blood flow patterns in the cerebral arterial circulation of individual (human) subjects

There is a growing research interest in quantifying blood flow distribution for the entire cerebral circulation to sharpen diagnosis and improve treatment options for cerebrovascular disease of individual patients. We present a methodology to reconstruct subject-specific cerebral blood flow patterns in accordance with physiological and fluid mechanical principles and optimally informed by in vivo neuroimage data of cerebrovascular anatomy and arterial blood flow rates. We propose an inverse problem to infer blood flow distribution across the visible portion of the arterial network that best matches subject-specific anatomy and a given set of volumetric flow measurements. The optimization technique also mitigates the effect of uncertainties by reconciling incomplete flow data and by dissipating unavoidable acquisition errors associated with medical imaging data.

Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data

Mathematical models are widely recognized as a valuable tool for cardiovascular diagnosis and the study of circulatory diseases, especially to obtain data that require otherwise invasive measurements. To correctly simulate body hemodynamics, the viscoelastic properties of vessels walls are a key aspect to be taken into account as they play an essential role in cardiovascular behavior. The present work aims to apply the augmented fluid-structure interaction system of blood flow to real case studies to assess the validity of the model as a valuable resource to improve cardiovascular diagnostics and the treatment of pathologies. Main contributions of the paper include the evaluation of viscoelastic tube laws, estimation of viscoelastic parameters and comparison of models with literature results and in-vivo experiments. The ability of the model to correctly simulate pulse waveforms in single arterial segments is verified using literature benchmark test cases, designed taking into account a simple elastic behavior of the wall in the upper thoracic aorta and in the common carotid artery. Furthermore, in-vivo pressure waveforms, extracted from tonometric measurements performed on four human common carotid arteries and two common femoral arteries, are compared to numerical solutions. It is highlighted that the viscoelastic damping effect of arterial walls is required to avoid an overestimation of pressure peaks. Finally, an effective procedure to estimate the viscoelastic parameters of the model is herein proposed, which returns hysteresis curves of the common carotid arteries dissipating energy fractions in line with values calculated from literature hysteresis loops in the same vessel.

On the accuracy of displacement-based wave intensity analysis: Effect of vessel wall viscoelasticity and nonlinearity

Recent studies showed that wave intensity analysis (WIA) provides clinically valuable information about local and global cardiovascular function. Wave intensity (WI) is computed as the product of the pressure change and the velocity change during short time intervals. The major limitation of WIA in clinical practice is the need for invasive pressure measurement. Since vessel wall displacement can be measured non-invasively, the usage of WI will be expanded if the vessel wall dilation is used instead of pressure in derivation of WI waveform. Our goal in this study is to investigate the agreement between wall displacement-based WI and the pressure-based WI for different vessel wall models including linear elastic, nonlinear and viscoelastic cases. The arbitrary Eulerian Lagrangian finite element method is employed to solve the coupled fluid-structure interaction (FSI). Our computational models also include two types of vascular disease-related cases with geometrical irregularities, aneurysm and stenosis. Our results show that for vessels with linear elastic wall, the displacement-based WI is almost identical to the pressure-based WI. The existence of vessel irregularities does not impact the accuracy of displacement-based WI. However, in a viscoelastic wall where there is a phase difference between pressure and vessel wall dilation, displacement-based WI deviated from pressure-based WI. The error associated with this phase difference increased nonlinearly with increasing viscosity. This results in a maximum error of 6.8% and 7.13% for a regular viscoelastic vessel wall and an irregular viscoelastic vessel wall, respectively. A separate analysis has also been performed on the agreement of backward and forward running waves extracted from a decomposition of the displacement-based and pressure-based WI. Our findings suggest that displacement-based WI is a reliable method of WIA for large central arteries that do not show viscoelastic behaviors. This can be clinically significant since the required information can be measured non-invasively.

Detection of arterial wall abnormalities via Bayesian model selection

Patient-specific modelling of haemodynamics in arterial networks has so far relied on parameter estimation for inexpensive or small-scale models. We describe here a Bayesian uncertainty quantification framework which makes two major advances: an efficient parallel implementation, allowing parameter estimation for more complex forward models, and a system for practical model selection, allowing evidence-based comparison between distinct physical models. We demonstrate the proposed methodology by generating simulated noisy flow velocity data from a branching arterial tree model in which a structural defect is introduced at an unknown location; our approach is shown to accurately locate the abnormality and estimate its physical properties even in the presence of significant observational and systemic error. As the method readily admits real data, it shows great potential in patient-specific parameter fitting for haemodynamical flow models.

Tapered vs. Uniform Tube-Load Modeling of Blood Pressure Wave Propagation in Human Aorta

In this paper, tapered vs. uniform tube-load models are comparatively investigated as mathematical representation for blood pressure (BP) wave propagation in human aorta. The relationship between the aortic inlet and outlet BP waves was formulated based on the exponentially tapered and uniform tube-load models. Then, the validity of the two tube-load models was comparatively investigated by fitting them to the experimental aortic and femoral BP waveform signals collected from 13 coronary artery bypass graft surgery patients. The two tube-load models showed comparable goodness of fit: (i) the root-mean-squared error (RMSE) was 3.3+/−1.1 mmHg in the tapered tube-load model and 3.4+/−1.1 mmHg in the uniform tube-load model; and (ii) the correlation was r = 0.98+/−0.02 in the tapered tube-load model and r = 0.98+/−0.01 mmHg in the uniform tube-load model. They also exhibited frequency responses comparable to the non-parametric frequency response derived from the aortic and femoral BP waveforms in most patients. Hence, the uniform tube-load model was superior to its tapered counterpart in terms of the Akaike Information Criterion (AIC). In general, the tapered tube-load model yielded the degree of tapering smaller than what is physiologically relevant: the aortic inlet-outlet radius ratio was estimated as 1.5 on the average, which was smaller than the anatomically plausible typical radius ratio of 3.5 between the ascending aorta and femoral artery. When the tapering ratio was restricted to the vicinity of the anatomically plausible typical value, the exponentially tapered tube-load model tended to underperform the uniform tube-load model (RMSE: 3.9+/−1.1 mmHg; r = 0.97+/−0.02). It was concluded that the uniform tube-load model may be more robust and thus preferred as the representation for BP wave propagation in human aorta; compared to the uniform tube-load model, the exponentially tapered tube-load model may not provide valid physiological insight on the aortic tapering, and its efficacy on the goodness of fit may be only marginal.

Pulse transit time estimation of aortic pulse wave velocity and blood pressure using machine learning and simulated training data

Recent developments in cardiovascular modelling allow us to simulate blood flow in an entire human body. Such model can also be used to create databases of virtual subjects, with sizes limited only by computational resources. In this work, we study if it is possible to estimate cardiovascular health indices using machine learning approaches. In particular, we carry out theoretical assessment of estimating aortic pulse wave velocity, diastolic and systolic blood pressure and stroke volume using pulse transit/arrival timings derived from photopletyshmography signals. For predictions, we train Gaussian process regression using a database of virtual subjects generated with a cardiovascular simulator. Simulated results provides theoretical assessment of accuracy for predictions of the health indices. For instance, aortic pulse wave velocity can be estimated with a high accuracy (r > 0.9) when photopletyshmography is measured from left carotid artery using a combination of foot-to-foot pulse transmit time and peak location derived for the predictions. Similar accuracy can be reached for diastolic blood pressure, but predictions of systolic blood pressure are less accurate (r > 0.75) and the stroke volume predictions are mostly contributed by heart rate.

Recently there has been a strong trend for self-monitoring of your cardiovascular health and new wearable sport trackers and mobile applications are coming to the market everyday. However, such solutions are mostly taking advantage of heart rate measurement. Other health indices such as blood pressure and pulse wave velocity reflecting to the condition of cardiovascular system would also be of great interest, but such solutions for continuous monitoring are barely existing or are at least unreliable. In this paper, we use computational modelling to assess theoretical capabilities of such measurements. We concentrate on predicting health indices using on pulse transmit time type of measurements. Such measurements could be carried out, for example, with photopletyshmography sensor or an optical sensor already found from several wearable sport trackers. We use cardiovascular modelling to create a database of “virtual subjects”, which is applied with machine learning to construct predictors for health indices. Our findings suggest that aortic pulse wave velocity and diastolic blood pressured could be predicted with a high accuracy, but predictions of systolic blood pressure are less accurate.

CardioFAN: open source platform for noninvasive assessment of pulse transit time and pulsatile flow in hyperelastic vascular networks

A profound analysis of pressure and flow wave propagation in cardiovascular systems is the key in noninvasive assessment of hemodynamic parameters. Pulse transit time (PTT), which closely relates to the physical properties of the cardiovascular system, can be linked to variations of blood pressure and stroke volume to provide information for patient-specific clinical diagnostics. In this work, we present mathematical and numerical tools, capable of accurately predicting the PTT, local pulse wave velocity, vessel compliance, and pressure/flow waveforms, in a viscous hyperelastic cardiovascular network. A new one-dimensional framework, entitled cardiovascular flow analysis (CardioFAN), is presented to describe the pulsatile fluid-structure interaction in the hyperelastic arteries, where pertaining hyperbolic equations are solved using a high-resolution total variation diminishing Lax-Wendroff method. The computational algorithm is validated against well-known numerical, in vitro and in vivo data for networks of main human arteries with 55, 37 and 26 segments, respectively. PTT prediction is improved by accounting for hyperelastic nonlinear waves between two arbitrary sections of the arterial tree. Consequently, arterial compliance assignments at each segment are improved in a personalized model of the human aorta and supra-aortic branches with 26 segments, where prior in vivo data were available for comparison. This resulted in a 1.5% improvement in overall predictions of the waveforms, or average relative errors of 5.5% in predicting flow, luminal area and pressure waveforms compared to prior in vivo measurements. The open source software, CardioFAN, can be calibrated for arbitrary patient-specific vascular networks to conduct noninvasive diagnostics.

A novel, FFT-based one-dimensional blood flow solution method for arterial network

In the present work, we propose an FFT-based method for solving blood flow equations in an arterial network with variable properties and geometrical changes. An essential advantage of this approach is in correctly accounting for the vessel skin friction through the use of Womersley solution. To incorporate nonlinear effects, a novel approximation method is proposed to enable calculation of nonlinear corrections. Unlike similar methods available in the literature, the set of algebraic equations required for every harmonic is constructed automatically. The result is a generalized, robust and fast method to accurately capture the increasing pulse wave velocity downstream as well as steepening of the pulse front. The proposed method is shown to be appropriate for incorporating correct convection and diffusion coefficients. We show that the proposed method is fast and accurate and it can be an effective tool for 1D modelling of blood flow in human arterial networks.

A novel porous media-based approach to outflow boundary resistances of 1D arterial blood flow models

In this paper we introduce a novel method for prescribing terminal boundary conditions in one-dimensional arterial flow networks. This is carried out by coupling the terminal arterial vessel with a poro-elastic tube, representing the flow resistance offered by microcirculation. The performance of the proposed porous media-based model has been investigated through several different numerical examples. First, we investigate model parameters that have a profound influence on the flow and pressure distributions of the system. The simulation results have been compared against the waveforms generated by three elements (RCR) Windkessel model. The proposed model is also integrated into a realistic arterial tree, and the results obtained have been compared against experimental data at different locations of the network. The accuracy and simplicity of the proposed model demonstrates that it can be an excellent alternative for the existing models.

An efficient full space-time discretization method for subject-specific hemodynamic simulations of cerebral arterial blood flow with distensible wall mechanics

A computationally inexpensive mathematical solution approach using orthogonal collocations for space discretization with temporal Fourier series is proposed to compute subject-specific blood flow in distensible vessels of large cerebral arterial networks. Several models of wall biomechanics were considered to assess their impact on hemodynamic predictions. Simulations were validated against in vivo blood flow measurements in six human subjects. The average root-mean-square relative differences were found to be less than 4.3% for all subjects with a linear elastic wall model. This discrepancy decreased further in a viscoelastic Kelvin-Voigt biomechanical wall. The results provide support for the use of collocation-Fourier series approach to predict clinically relevant blood flow distribution and collateral blood supply in large portions of the cerebral circulation at reasonable computational costs. It thus opens the possibility of performing computationally inexpensive subject-specific simulations that are robust and fast enough to predict clinical results in real time on the same day.

An In-Vitro Flow Study Using an Artificial Circle of Willis Model for Validation of an Existing One-Dimensional Numerical Model

A one-dimensional (1D) numerical model has been previously developed to investigate the hemodynamics of blood flow in the entire human vascular network. In the current work, an experimental study of water-glycerin mixture flow in a 3D-printed silicone model of an anatomically accurate, complete circle of Willis (CoW) was conducted to investigate the flow characteristics in comparison with the simulated results by the 1D numerical model. In the experiment, the transient flow and pressure waveforms were measured at 13 selected segments within the flow network for comparisons. In the 1D simulation, the initial parameters of the vessel network were obtained by a direct measurement of the tubes in the experimental setup. The results verified that the 1D numerical model is able to capture the main features of the experimental pressure and flow waveforms with good reliability. The mean flow rates measurement results agree with the predictions of the 1D model with an overall difference of less than 1%. Further experiment might be needed to validate the 1D model in capturing pressure waveforms.

A data-driven model to study utero-ovarian blood flow physiology during pregnancy

In this paper, we describe a mathematical model of the cardiovascular system in human pregnancy. An automated, closed-loop 1D-0D modelling framework was developed, and we demonstrate its efficacy in (1) reproducing measured multi-variate cardiovascular variables (pulse pressure, total peripheral resistance and cardiac output) and (2) providing automated estimates of variables that have not been measured (uterine arterial and venous blood flow, pulse wave velocity, pulsatility index). This is the first model capable of estimating volumetric blood flow to the uterus via the utero-ovarian communicating arteries. It is also the first model capable of capturing wave propagation phenomena in the utero-ovarian circulation, which are important for the accurate estimation of arterial stiffness in contemporary obstetric practice. The model will provide a basis for future studies aiming to elucidate the physiological mechanisms underlying the dynamic properties (changing shapes) of vascular flow waveforms that are observed with advancing gestation. This in turn will facilitate the development of methods for the earlier detection of pathologies that have an influence on vascular structure and behaviour.

Validation and Application of a Physically Nonlinear ID Computational Model for Bifurcated Arterial Networks

Reduced fluid-structure interaction models are the key component of hemodynamic simulation. In this work, a multi-purpose computational model applicable to specific physiological components such as arterial, venous and cerebrospinal fluid circulatory systems has been developed based on the Hamilton's variational principle. This model encompasses a viscous Newtonian fluid structure interaction (FSI) framework for the large compliant bifurcated arterial networks and its subsystems. This approach provides the groundworks for a correct formulation of reduced FSI models with an account for arbitrary non-linear viscoelastic properties of a compliant vascular tree. The hyperbolic properties of the derived mathematical model are analyzed and used to construct the Lax-Wendroff finite volume numerical scheme, with second order accuracy in time and space. The computational algorithm is validated against well-known numerical and in vitro experimental data reported in the literature for the case of human arterial trees, comprising 55 and 37 main arterial vessels. Utilizing the physics based nonlinear constitutive framework, this model can be adequately tested, calibrated and applied for patient-specific clinical diagnosis and prediction.