A benchmark study of numerical schemes for one-dimensional arterial blood flow modelling

Modelling blood flow in the circle of Willis in continuous flow left ventricular assist devices: possible relevance to strokes

Purpose

Despite significant improvements in the design and performance of continuous flow left ventricular assist devices (CFLVADs), one of the most important reasons hampering further penetration of this technology is the occurrence of adverse events, especially strokes. One of the well-known risk factors for strokes is hypertension which is particularly common in patients undergoing a CFLVAD implant. While the device is implanted in the heart, strokes happen due to pathology in the brain and we hypothesised that modelling the blood flow in the circle of Willis might shed light on the causation of strokes in this situation.The aim of the study was two-fold:1. What is the reason for hypertension in CFLVADs? Are there physical factors at play, besides neurohumoral mechanisms?2. Do anatomical factors in the circle of Willis play a role in the causation of strokes in these patients?

Methods

The circle of Willis is often incomplete and has a number of anatomical variations, the commonest being the absence of the posterior communicating artery. Hypertension is common after CFLVAD implantation and is also a well-known risk factor for strokes. We examined the blood pressure in the cerebral circulation with pulsatile and non-pulsatile flow for identical conditions and the effect of the absence of the posterior communicating artery on regional cerebral blood flow and pressure. One-dimensional blood flow model was used, taking into account wave propagation and reflections and physiological data obtained from anatomically detailed arterial network (ADAN86) which has data from 86 arteries including detailed cerebral network.

Results

The mean arterial pressure was significantly higher in the non-pulsatile blood flow of CFLVADs compared to pulsatile flow, for identical conditions, across all arteries. With increasing imparted pulsatility to CFLVAD flow, the mean arterial pressure progressively decreased. Isolated absence of the posterior communicating artery had no effect on the flow as well as pressure in the middle cerebral artery. However, when combined with the absence of flow in the ipsilateral carotid artery, the flow as well as the pressure decreased very significantly in both continuous and pulsatile flow situations.

Conclusions

Physiologically significant pulsatility in CFLVADs can have important clinical advantages in lowering of blood pressure which can lead to lower incidence of strokes, pump thrombosis, gastrointestinal (GI) bleeds, and aortic incompetence. Patient-specific anatomical variations in the circle of Willis, especially the absence of the posterior communicating artery, can have important consequences in regional cerebral perfusion under some circumstances.

Graphical abstract

Uncertainty Quantification and Sensitivity Analysis for Non-invasive Model-Based Instantaneous Wave-Free Ratio Prediction

The main objectives of this work are to validate a 1D-0D unsteady solver with a distributed stenosis model for the patient-specific estimation of resting haemodynamic indices and to assess the sensitivity of instantaneous wave-free ratio (iFR) predictions to uncertainties in input parameters. We considered 52 patients with stable coronary artery disease, for which 81 invasive iFR measurements were available. We validated the performance of our solver compared to 3D steady-state and transient results and invasive measurements. Next, we used a polynomial chaos approach to characterise the uncertainty in iFR predictions based on the inputs associated with boundary conditions (coronary flow, compliance and aortic/left ventricular pressures) and vascular geometry (radius). Agreement between iFR and the ratio between cardiac cycle averaged distal and aortic pressure waveforms (restingP d / P a $$ {P}_d/{P}_a $$ ) obtained through 1D-0D and 3D models was satisfactory, with a bias of 0.0-0.005 (±0.016-0.026). The sensitivity analysis showed that iFR estimation is mostly affected by uncertainties in vascular geometry and coronary flow (steady-state parameters). In particular, our 1D-0D method overestimates invasive iFR measurements, with a bias of -0.036 (±0.101), indicating that better flow estimates could significantly improve our modelling pipeline. Conversely, we showed that standard pressure waveforms could be used for simulations, since the impact of uncertainties related to inlet-pressure waveforms on iFR prediction is negligible. Furthermore, while compliance is the most relevant transient parameter, its effect on iFR estimates is negligible compared to that of vascular geometry and flow. Finally, we observed a strong correlation between iFR and restingP d / P a $$ {P}_d/{P}_a $$ , suggesting that steady-state simulations could replace unsteady simulations for iFR prediction.

Noninvasive estimation of central blood pressure through fluid-structure interaction modeling

Central blood pressure (cBP) is considered a superior indicator of cardiovascular fitness than brachial blood pressure (bBP). Even though bBP is easy to measure noninvasively, it is usually higher than cBP due to pulse wave amplification, characterized by the gradual increase in peak systolic pressure during pulse wave propagation. In this study, we aim to develop an individualized transfer function that can accurately estimate cBP from bBP. We first construct a three-dimensional, patient-specific model of the upper limb arterial system using fluid-structure interaction simulations, incorporating variable material properties and complex boundary conditions. Then, we develop an analytical brachial-aortic transfer function based on novel solutions for compliant vessels. The accuracy of this transfer function is successfully validated against numerical simulation results, which effectively reproduce pulse wave propagation and amplification, with key hemodynamic parameters falling within the range of clinical measurements. Further analysis of the transfer function reveals that cBP is a linear combination of bBP and aortic flow rate in the frequency domain, with the coefficients determined by vessel geometry, material properties, and boundary conditions. Additionally, bBP primarily contributes to the steady component of cBP, while the aortic flow rate is responsible for the pulsatile component. Furthermore, local sensitivity analysis indicates that the lumen radius is the most influential parameter in accurately estimating cBP. Although not directly applicable clinically, the proposed transfer function enhances understanding of the underlying physics-highlighting the importance of aortic flow and lumen radius-and can guide the development of more practical transfer functions.

openBF: an open-source finite volume 1D blood flow solver

Computational simulations are widely adopted in cardiovascular biomechanics because of their capability of producing physiological data otherwise impossible to measure with non-invasive modalities.Objective.This study presents openBF, a computational library for simulating the blood dynamics in the cardiovascular system.Approach.openBF adopts a one-dimensional viscoelastic representation of the arterial system, and is coupled with zero-dimensional windkessel models at the outlets. Equations are solved by means of the finite-volume method and the code is written in Julia. We assess its predictions by performing a multiscale validation study on several domains available from the literature.Main results.At all scales, which range from individual arteries to a population of virtual subjects, openBF's solution show excellent agreement with the solutions from existing software. For reported simulations, openBF requires low computational times.Significance.openBF is easy to install, use, and deploy on multiple platforms and architectures, and gives accurate prediction of blood dynamics in short time-frames. It is actively maintained and available open-source on GitHub, which favours contributions from the biomechanical community.

An In Silico Modelling Approach to Predict Hemodynamic Outcomes in Diabetic and Hypertensive Kidney Disease

Early diagnosis of kidney disease remains an unmet clinical challenge, preventing timely and effective intervention. Diabetes and hypertension are two main causes of kidney disease, can often appear together, and can only be distinguished by invasive biopsy. In this study, we developed a modelling approach to simulate blood velocity, volumetric flow rate, and pressure wave propagation in arterial networks of ageing, diabetic, and hypertensive virtual populations. The model was validated by comparing our predictions for pressure, volumetric flow rate and waveform-derived indexes with in vivo data on ageing populations from the literature. The model simulated the effects of kidney disease, and was calibrated to align quantitatively with in vivo data on diabetic and hypertensive nephropathy from the literature. Our study identified some potential biomarkers extracted from renal blood flow rate and flow pulsatility. For typical patient age groups, resistive index values were 0.69 (SD 0.05) and 0.74 (SD 0.02) in the early and severe stages of diabetic nephropathy, respectively. Similar trends were observed in the same stages of hypertensive nephropathy, with a range from 0.65 (SD 0.07) to 0.73 (SD 0.05), respectively. Mean renal blood flow rate through a single diseased kidney ranged from 329 (SD 40, early) to 317 (SD 38, severe) ml/min in diabetic nephropathy and 443 (SD 54, early) to 388 (SD 47, severe) ml/min in hypertensive nephropathy, showing potential as a biomarker for early diagnosis of kidney disease. This modelling approach demonstrated its potential application in informing biomarker identification and facilitating the setup of clinical trials.

Computational modelling of cardiovascular pathophysiology to risk stratify commercial spaceflight

For more than 60 years, humans have travelled into space. Until now, the majority of astronauts have been professional, government agency astronauts selected, in part, for their superlative physical fitness and the absence of disease. Commercial spaceflight is now becoming accessible to members of the public, many of whom would previously have been excluded owing to unsatisfactory fitness or the presence of cardiorespiratory diseases. While data exist on the effects of gravitational and acceleration (G) forces on human physiology, data on the effects of the aerospace environment in unselected members of the public, and particularly in those with clinically significant pathology, are limited. Although short in duration, these high acceleration forces can potentially either impair the experience or, more seriously, pose a risk to health in some individuals. Rather than expose individuals with existing pathology to G forces to collect data, computational modelling might be useful to predict the nature and severity of cardiovascular diseases that are of sufficient risk to restrict access, require modification, or suggest further investigation or training before flight. In this Review, we explore state-of-the-art, zero-dimensional, compartmentalized models of human cardiovascular pathophysiology that can be used to simulate the effects of acceleration forces, homeostatic regulation and ventilation-perfusion matching, using data generated by long-arm centrifuge facilities of the US National Aeronautics and Space Administration and the European Space Agency to risk stratify individuals and help to improve safety in commercial suborbital spaceflight.

Computational haemodynamics for pulmonary valve replacement by means of a reduced fluid-structure interaction model

Pulmonary valve replacement (PVR) consists of substituting a patient's original valve with a prosthetic one, primarily addressing pulmonary valve insufficiency, which is crucially relevant in Tetralogy of Fallot repairment. While extensive clinical and computational literature on aortic and mitral valve replacements is available, PVR's post-procedural haemodynamics in the pulmonary artery and the impact of prosthetic valve dynamics remain significantly understudied. Addressing this gap, we introduce a reduced Fluid-Structure Interaction (rFSI) model, applied for the first time to the pulmonary valve. This model couples a three-dimensional computational representation of pulmonary artery haemodynamics with a one-degree-of-freedom model to account for valve structural mechanics. Through this approach, we analyse patient-specific haemodynamics pre and post PVR. Patient-specific geometries, reconstructed from CT scans, are virtually equipped with a template valve geometry. Boundary conditions for the model are established using a lumped-parameter model, fine-tuned based on clinical patient data. Our model accurately reproduces patient-specific haemodynamic changes across different scenarios: pre-PVR, six months post-PVR, and a follow-up condition after a decade. It effectively demonstrates the impact of valve implantation on sustaining the diastolic pressure gradient across the valve. The numerical results indicate that our valve model is able to reproduce overall physiological and/or pathological conditions, as preliminary assessed on two different patients. This promising approach provides insights into post-PVR haemodynamics and prosthetic valve effects, shedding light on potential implications for patient-specific outcomes.

Global sensitivity analysis with multifidelity Monte Carlo and polynomial chaos expansion for vascular haemodynamics

Computational models of the cardiovascular system are increasingly used for the diagnosis, treatment, and prevention of cardiovascular disease. Before being used for translational applications, the predictive abilities of these models need to be thoroughly demonstrated through verification, validation, and uncertainty quantification. When results depend on multiple uncertain inputs, sensitivity analysis is typically the first step required to separate relevant from unimportant inputs, and is key to determine an initial reduction on the problem dimensionality that will significantly affect the cost of all downstream analysis tasks. For computationally expensive models with numerous uncertain inputs, sample-based sensitivity analysis may become impractical due to the substantial number of model evaluations it typically necessitates. To overcome this limitation, we consider recently proposed Multifidelity Monte Carlo estimators for Sobol' sensitivity indices, and demonstrate their applicability to an idealized model of the common carotid artery. Variance reduction is achieved combining a small number of three-dimensional fluid-structure interaction simulations with affordable one- and zero-dimensional reduced-order models. These multifidelity Monte Carlo estimators are compared with traditional Monte Carlo and polynomial chaos expansion estimates. Specifically, we show consistent sensitivity ranks for both bi- (1D/0D) and tri-fidelity (3D/1D/0D) estimators, and superior variance reduction compared to traditional single-fidelity Monte Carlo estimators for the same computational budget. As the computational burden of Monte Carlo estimators for Sobol' indices is significantly affected by the problem dimensionality, polynomial chaos expansion is found to have lower computational cost for idealized models with smooth stochastic response.

Age and sex-dependent sensitivity analysis of a common carotid artery model

The common carotid artery (CCA) is an accessible and informative site for assessing cardiovascular function which makes it a prime candidate for clinically relevant computational modelling. The interpretation of supplemental information possible through modelling is encumbered by measurement uncertainty and population variability in model parameters. The distribution of model parameters likely depends on the specific sub-population of interest and delineation based on sex, age or health status may correspond to distinct ranges of typical parameter values. To assess this impact in a 1D-CCA-model, we delineated specific sub-populations based on age, sex and health status and carried out uncertainty quantification and sensitivity analysis for each sub-population. We performed a structured literature review to characterize sub-population-specific variabilities for eight model parameters without consideration of health status; variations for a healthy sub-populations were based on previously established references values. The variabilities of diameter and distensibility found in the literature review differed from those previously established in a healthy population. Model diameter change and pulse pressure were most sensitive to variations in distensibility, while pressure was most sensitive to resistance in the Windkessel model for all groups. Uncertainties were lower when variabilities were based on a healthy sub-population; however, the qualitative distribution of sensitivity indices was largely similar between the healthy and general population. Average sensitivity of the pressure waveform showed a moderate dependence on age with decreasing sensitivity to distal resistance and increasing sensitivity to distensibility and diameter. The female population was less sensitive to variations in diameter but more sensitive to distensibility coefficient than the male population. Overall, as hypothesized input variabilities differed between sub-populations and resulted in distinct uncertainties and sensitivities of the 1D-CCA-model outputs, particularly over age for the pressure waveform and between males and females for pulse pressure.

Tube law parametrization using in vitro data for one-dimensional blood flow in arteries and veins: TUBE LAW PARAMETRIZATION IN ARTERIES AND VEINS

The deformability of blood vessels in one-dimensional blood flow models is typically described through a pressure-area relation, known as the tube law. The most used tube laws take into account the elastic and viscous components of the tension of the vessel wall. Accurately parametrizing the tube laws is vital for replicating pressure and flow wave propagation phenomena. Here, we present a novel mathematical-property-preserving approach for the estimation of the parameters of the elastic and viscoelastic tube laws. Our goal was to estimate the parameters by using ovine and human in vitro data, while constraining them to meet prescribed mathematical properties. Results show that both elastic and viscoelastic tube laws accurately describe experimental pressure-area data concerning both quantitative and qualitative aspects. Additionally, the viscoelastic tube law can provide a qualitative explanation for the observed hysteresis cycles. The two models were evaluated using two approaches: (i) allowing all parameters to freely vary within their respective ranges and (ii) fixing some of the parameters. The former approach was found to be the most suitable for reproducing pressure-area curves.

A Modular Framework for Implicit 3D-0D Coupling in Cardiac Mechanics

In numerical simulations of cardiac mechanics, coupling the heart to a model of the circulatory system is essential for capturing physiological cardiac behavior. A popular and efficient technique is to use an electrical circuit analogy, known as a lumped parameter network or zero-dimensional (0D) fluid model, to represent blood flow throughout the cardiovascular system. Due to the strong physical interaction between the heart and the blood circulation, developing accurate and efficient numerical coupling methods remains an active area of research. In this work, we present a modular framework for implicitly coupling three-dimensional (3D) finite element simulations of cardiac mechanics to 0D models of blood circulation. The framework is modular in that the circulation model can be modified independently of the 3D finite element solver, and vice versa. The numerical scheme builds upon a previous work that combines 3D blood flow models with 0D circulation models (3D fluid - 0D fluid). Here, we extend it to couple 3D cardiac tissue mechanics models with 0D circulation models (3D structure - 0D fluid), showing that both mathematical problems can be solved within a unified coupling scheme. The effectiveness, temporal convergence, and computational cost of the algorithm are assessed through multiple examples relevant to the cardiovascular modeling community. Importantly, in an idealized left ventricle example, we show that the coupled model yields physiological pressure-volume loops and naturally recapitulates the isovolumic contraction and relaxation phases of the cardiac cycle without any additional numerical techniques. Furthermore, we provide a new derivation of the scheme inspired by the Approximate Newton Method of Chan (1985), explaining how the proposed numerical scheme combines the stability of monolithic approaches with the modularity and flexibility of partitioned approaches.

Efficient sensitivity analysis for biomechanical models with correlated inputs

In most variance-based sensitivity analysis (SA) approaches applied to biomechanical models, statistical independence of the model input is assumed. However, often the model inputs are correlated. This might alter the interpretation of the SA results, which may severely impact the guidance provided during model development and personalization. Potential reasons for the infrequent usage of SA techniques that account for input correlation are the associated high computational costs, especially for models with many parameters, and the fact that the input correlation structure is often unknown. The aim of this study was to propose an efficient correlated global sensitivity analysis method by applying a surrogate model-based approach. Furthermore, this article demonstrates how correlated SA should be interpreted and how the applied method can guide the modeler during model development and personalization, even when the correlation structure is not entirely known beforehand. The proposed methodology was applied to a typical example of a pulse wave propagation model and resulted in accurate SA results that could be obtained at a theoretically 27,000× lower computational cost compared to the correlated SA approach without employing a surrogate model. Furthermore, our results demonstrate that input correlations can significantly affect SA results, which emphasizes the need to thoroughly investigate the effect of input correlations during model development. We conclude that our proposed surrogate-based SA approach allows modelers to efficiently perform correlated SA to complex biomechanical models and allows modelers to focus on input prioritization, input fixing and model reduction, or assessing the dependency structure between parameters.

Learning reduced-order models for cardiovascular simulations with graph neural networks

Reduced-order models based on physics are a popular choice in cardiovascular modeling due to their efficiency, but they may experience loss in accuracy when working with anatomies that contain numerous junctions or pathological conditions. We develop one-dimensional reduced-order models that simulate blood flow dynamics using a graph neural network trained on three-dimensional hemodynamic simulation data. Given the initial condition of the system, the network iteratively predicts the pressure and flow rate at the vessel centerline nodes. Our numerical results demonstrate the accuracy and generalizability of our method in physiological geometries comprising a variety of anatomies and boundary conditions. Our findings demonstrate that our approach can achieve errors below 3% for pressure and flow rate, provided there is adequate training data. As a result, our method exhibits superior performance compared to physics-based one-dimensional models while maintaining high efficiency at inference time.

Advanced waveform analysis of the photoplethysmogram signal using complementary signal processing techniques for the extraction of biomarkers of cardiovascular function

Introduction

Photoplethysmogram signals from wearable devices typically measure heart rate and blood oxygen saturation, but contain a wealth of additional information about the cardiovascular system. In this study, we compared two signal-processing techniques: fiducial point analysis and Symmetric Projection Attractor Reconstruction, on their ability to extract new cardiovascular information from a photoplethysmogram signal. The aim was to identify fiducial point analysis and Symmetric Projection Attractor Reconstruction indices that could classify photoplethysmogram signals, according to age, sex and physical activity.

Methods

Three datasets were used: an in-silico dataset of simulated photoplethysmogram waves for healthy male participants (25-75 years old); an in-vivo dataset containing 10-min photoplethysmogram recordings from 57 healthy subjects at rest (18-39 or > 70 years old; 53% female); and an in-vivo dataset containing photoplethysmogram recordings collected for 4 weeks from a single subject, in daily life. The best-performing indices from the in-silico study (5/48 fiducial point analysis and 6/49 Symmetric Projection Attractor Reconstruction) were applied to the in-vivo datasets.

Results

Key fiducial point analysis and Symmetric Projection Attractor Reconstruction indices, which showed the greatest differences between groups, were found to be consistent across datasets. These indices were related to systolic augmentation, diastolic peak positioning and prominence, and waveform variability. Both fiducial point analysis and Symmetric Projection Attractor Reconstruction techniques provided indices that supported the classification of age and physical activity, but not sex.

Conclusions

Both fiducial point analysis and Symmetric Projection Attractor Reconstruction techniques demonstrated utility in identifying cardiovascular differences between individuals and within an individual over time. Future research should investigate the potential utility of these techniques for extracting information on fitness and disease, to support healthcare-decision making.

Mechanistic insights on age-related changes in heart-aorta-brain hemodynamic coupling using a pulse wave model of the entire circulatory system

Age-related changes in aortic biomechanics can impact the brain by reducing blood flow and increasing pulsatile energy transmission. Clinical studies have shown that impaired cardiac function in patients with heart failure is associated with cognitive impairment. Although previous studies have attempted to elucidate the complex relationship between age-associated aortic stiffening and pulsatility transmission to the cerebral network, they have not adequately addressed the effect of interactions between aortic stiffness and left ventricle (LV) contractility (neither on energy transmission nor on brain perfusion). In this study, we use a well-established and validated one-dimensional blood flow and pulse wave computational model of the circulatory system to address how age-related changes in cardiac function and vasculature affect the underlying mechanisms involved in the LV-aorta-brain hemodynamic coupling. Our results reveal how LV contractility affects pulsatile energy transmission to the brain, even with preserved cardiac output. Our model demonstrates the existence of an optimal heart rate (near the normal human heart rate) that minimizes pulsatile energy transmission to the brain at different contractility levels. Our findings further suggest that the reduction in cerebral blood flow at low levels of LV contractility is more prominent in the setting of age-related aortic stiffening. Maintaining optimal blood flow to the brain requires either an increase in contractility or an increase in heart rate. The former consistently leads to higher pulsatile power transmission, and the latter can either increase or decrease subsequent pulsatile power transmission to the brain.NEW & NOTEWORTHY We investigated the impact of major aging mechanisms of the arterial system and cardiac function on brain hemodynamics. Our findings suggest that aging has a significant impact on heart-aorta-brain coupling through changes in both arterial stiffening and left ventricle (LV) contractility. Understanding the underlying physical mechanisms involved here can potentially be a key step for developing more effective therapeutic strategies that can mitigate the contributions of abnormal LV-arterial coupling toward neurodegenerative diseases and dementia.

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 novel physics-based model for fast computation of blood flow in coronary arteries

Blood flow and pressure calculated using the currently available methods have shown the potential to predict the progression of pathology, guide treatment strategies and help with postoperative recovery. However, the conspicuous disadvantage of these methods might be the time-consuming nature due to the simulation of virtual interventional treatment. The purpose of this study is to propose a fast novel physics-based model, called FAST, for the prediction of blood flow and pressure. More specifically, blood flow in a vessel is discretized into a number of micro-flow elements along the centerline of the artery, so that when using the equation of viscous fluid motion, the complex blood flow in the artery is simplified into a one-dimensional (1D) steady-state flow. We demonstrate that this method can compute the fractional flow reserve (FFR) derived from coronary computed tomography angiography (CCTA). 345 patients with 402 lesions are used to evaluate the feasibility of the FAST simulation through a comparison with three-dimensional (3D) computational fluid dynamics (CFD) simulation. Invasive FFR is also introduced to validate the diagnostic performance of the FAST method as a reference standard. The performance of the FAST method is comparable with the 3D CFD method. Compared with invasive FFR, the accuracy, sensitivity and specificity of FAST is 88.6%, 83.2% and 91.3%, respectively. The AUC of FFR_FAST is 0.906. This demonstrates that the FAST algorithm and 3D CFD method show high consistency in predicting steady-state blood flow and pressure. Meanwhile, the FAST method also shows the potential in detecting lesion-specific ischemia.

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.

A noninvasive method of estimating patient-specific left ventricular pressure waveform

BACKGROUND AND OBJECTIVE

The left ventricular pressure waveform is indispensable for the construction of the pressure strain loop when investigating coronary artery disease (CAD) patients. In previous studies by others, exclusion of CAD patients has not allowed a reliable estimation of the left ventricular pressure waveform from the pressure strain loop of these patients. To remedy this, we propose a patient-specific noninvasive method for the estimation of left ventricular pressure.

METHODS

A simplified systemic circulation model consisting primarily of a single fiber model and a 1D simulation of the arterial tree was used. Sensitivity analysis based on the Morris method was performed to select a subset of the important parameters. Following this, the important parameter subset and the set of all the parameters were identified in the model using the pressure waveform of a peripheral artery as input, in a two-step process. In addition, the left ventricular pressure waveform was estimated using the set of all parameters.

RESULTS

Reducing the size of the parameter subset significantly decreases the computational cost of parameter optimization in the first step and greatly simplifies the identification of the full parameter set in the second step. Comparison with the reference left ventricular pressure waveform from CAD patients, showed that the proposed method provides a good estimate of the reference left ventricular pressure waveform. The correlation coefficients between the estimated and reference were r = 0.907, r = 0.904 and r = 0.780 for systolic blood pressure, pulse pressure and mean blood pressure, respectively.

CONCLUSIONS

This work may provide a convenient surrogate for the estimation of the left ventricular pressure waveform.

Modelling collateral flow and thrombus permeability during acute ischaemic stroke

The presence of collaterals and high thrombus permeability are associated with good functional outcomes after an acute ischaemic stroke. We aim to understand the combined effect of the collaterals and thrombus permeability on cerebral blood flow during an acute ischaemic stroke. A cerebral blood flow model including the leptomeningeal collateral circulation is used to simulate cerebral blood flow during an acute ischaemic stroke. The collateral circulation is varied to capture the collateral scores: absent, poor, moderate and good. Measurements of the transit time, void fraction and thrombus length in acute ischaemic stroke patients are used to estimate thrombus permeability. Estimated thrombus permeability ranges between 10^-7 and 10^-4 mm². Measured flow rates through the thrombus are small and the effect of a permeable thrombus on brain perfusion during stroke is small compared with the effect of collaterals. Our simulations suggest that the collaterals are a dominant factor in the resulting infarct volume after a stroke.

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.

A coupling strategy for a first 3D-1D model of the cardiovascular system to study the effects of pulse wave propagation on cardiac function

A key factor governing the mechanical performance of the heart is the bidirectional coupling with the vascular system, where alterations in vascular properties modulate the pulsatile load imposed on the heart. Current models of cardiac electromechanics (EM) use simplified 0D representations of the vascular system when coupling to anatomically accurate 3D EM models is considered. However, these ignore important effects related to pulse wave transmission. Accounting for these effects requires 1D models, but a 3D-1D coupling remains challenging. In this work, we propose a novel, stable strategy to couple a 3D cardiac EM model to a 1D model of blood flow in the largest systemic arteries. For the first time, a personalised coupled 3D-1D model of left ventricle and arterial system is built and used in numerical benchmarks to demonstrate robustness and accuracy of our scheme over a range of time steps. Validation of the coupled model is performed by investigating the coupled system's physiological response to variations in the arterial system affecting pulse wave propagation, comprising aortic stiffening, aortic stenosis or bifurcations causing wave reflections. Our first 3D-1D coupled model is shown to be efficient and robust, with negligible additional computational costs compared to 3D-0D models. We further demonstrate that the calibrated 3D-1D model produces simulated data that match with clinical data under baseline conditions, and that known physiological responses to alterations in vascular resistance and stiffness are correctly replicated. Thus, using our coupled 3D-1D model will be beneficial in modelling studies investigating wave propagation phenomena.

Euler-Lagrange numerical simulation of improved magnetic drug delivery in a three-dimensional CT-based carotid artery bifurcation

BACKGROUND AND OBJECTIVE

Magnetic drug targeting (MDT) is a promising method to improve the therapy efficiency for cardiovascular diseases (CVDs) and cancers. In MDT, therapeutic agents are bonded to superparamagnetic iron oxide nanoparticle (SPION) cores and then are guided toward the damaged tissue through a magnetic field. Fundamentally, it's vital to steer the SPIONs to the desired location to increase the capture efficiency at the target lesion. Hence, the present study aims to enhance the drug delivery to the desired branch in a carotid bifurcation. Besides, it is tried to decrement the particles' entry to the unwanted outlet by using four different magnet configurations (with a maximum magnetic flux density of 0.4 T) implanted adjacent to the artery wall. Also, the effect of particles' diameter -ranging from 200 nm to 2 µm- on the drug delivery performance is studied in the four cases.

METHODS

The Eulerian-Lagrangian approach with one-way coupling is employed for numerical simulation of the problem using the finite element method (FEM). The dominant forces acting on particles are drag and magnetophoretic. A computed tomography (CT) model of the carotid bifurcation is adopted to have a 3D realistic geometry. The blood flow is considered to be laminar, incompressible, pulsatile, and non-Newtonian. Boundary conditions are applied using the three-element Windkessel equation.

RESULTS

Results are presented in terms of velocity, pressure, magnetic field flux density, wall shear stress, and streamlines. Also, the number of particles in each direction is presented for the four studied cases. The results show that using proper magnets configurations makes it possible to guide more particles to the desired branch (up to 4%) while preventing particles from entering the unwanted branch (up to 13%). By defining connectivity between oscillatory shear index (OSI) value and magnetic drug delivery efficacy, it becomes clear that places with lower OSI values are more proper to place the magnets than areas with higher OSI values.

CONCLUSIONS

It was observed that increasing the diameter of particles does not necessarily result in a higher drug delivery efficiency. The configuration of the magnets and the size of particles are the main affecting parameters that should be chosen precisely to meet the most efficient drug delivery performance. Magnetic drug targeting (MDT) is a promising method to improve the therapy efficiency for cardiovascular diseases (CVDs) and cancers. Fundamentally, it's vital to steer the superparamagnetic iron oxide nanoparticles (SPIONs) to the target lesion location to increase the capture efficiency. Hence, the present study aims to enhance the drug delivery to the desired branch in a 3D carotid bifurcation. Besides, it is tried to decrement the particles' entry to the unwanted outlet by using four different magnet configurations implanted adjacent to the artery wall. The Eulerian-Lagrangian approach with one-way coupling is employed for numerical simulation of the problem using the finite element method (FEM). The dominant forces acting on particles are drag and magnetophoretic. The blood flow is laminar, incompressible, pulsatile, and non-Newtonian. The results show that it is possible to guide more particles to the desired branch (up to 4%) while preventing particles from entering the unwanted branch (up to 13%). By defining connectivity between oscillatory shear index (OSI) value and magnetic drug delivery efficacy, it becomes clear that places with lower OSI values are more proper to place the magnets than areas with higher OSI values.

Improving the accuracy and robustness of carotid-femoral pulse wave velocity measurement using a simplified tube-load model

Arterial stiffness, as measured by pulse wave velocity, for the early non-invasive screening of cardiovascular disease is becoming ever more widely used and is an independent prognostic indicator for a variety of pathologies including arteriosclerosis. Carotid-femoral pulse wave velocity (cfPWV) is regarded as the gold standard for aortic stiffness. Existing algorithms for cfPWV estimation have been shown to have good repeatability and accuracy, however, further assessment is needed, especially when signal quality is compromised. We propose a method for calculating cfPWV based on a simplified tube-load model, which allows for the propagation and reflection of the pulse wave. In-vivo cfPWV measurements from 57 subjects and numerical cfPWV data based on a one-dimensional model were used to assess the method and its performance was compared to three other existing approaches (waveform matching, intersecting tangent, and cross-correlation). The cfPWV calculated using the simplified tube-load model had better repeatability than the other methods (Intra-group Correlation Coefficient, ICC = 0.985). The model was also more accurate than other methods (deviation, 0.13 ms⁻¹) and was more robust when dealing with noisy signals. We conclude that the determination of cfPWV based on the proposed model can accurately and robustly evaluate arterial stiffness.

Steady-state solutions of one-dimensional equations of non-Newtonian hemodynamics

This paper is devoted to obtaining and analysis of steady-state solutions of one-dimensional equations for the simulation of blood flow when the non-Newtonian nature of blood is taken into account. The models, based on the rheological relations, widely used for the blood, are considered. The expressions for the nonlinear frictional term are presented. For the Power Law, Simplified Cross, and Quemada models, the exact integrals of the nonlinear ordinary differential equation, obtained from the averaged momentum equation, are obtained. It is demonstrated that several solutions exist for every rheological model, but the physically relevant solutions can be selected by the appropriate value of Mach number. The effects of the velocity profile and the value of hematocrit on the steady-state solutions are analyzed. It is demonstrated that the flattening of the velocity profile, which is typical for the blood, leads to the diminishing of the length of the interval, where the solution exists. The same effect is observed when the hematocrit value is increased.

A mathematical model of maternal vascular growth and remodeling and changes in maternal hemodynamics in uncomplicated pregnancy

The maternal vasculature undergoes tremendous growth and remodeling (G&R) that enables a > 15-fold increase in blood flow through the uterine vasculature from conception to term. Hemodynamic metrics (e.g., uterine artery pulsatility index, UA-PI) are useful for the prognosis of pregnancy complications; however, improved characterization of the maternal hemodynamics is necessary to improve prognosis. The goal of this paper is to develop a mathematical framework to characterize maternal vascular G&R and hemodynamics in uncomplicated human pregnancies. A validated 1D model of the human vascular tree from the literature was adapted and inlet blood flow waveforms at the ascending aorta at 4 week increments from 0 to 40 weeks of gestation were prescribed. Peripheral resistances of each terminal vessel were adjusted to achieve target flow rates and mean arterial pressure at each gestational age. Vessel growth was governed by wall shear stress (and axial lengthening in uterine vessels), and changes in vessel distensibility were related to vessel growth. Uterine artery velocity waveforms generated from this model closely resembled ultrasound results from the literature. The literature UA-PI values changed significantly across gestation, increasing in the first month of gestation, then dramatically decreasing from 4 to 20 weeks. Our results captured well the time-course of vessel geometry, material properties, and UA-PI. This 1D fluid-G&R model captured the salient hemodynamic features across a broad range of clinical reports and across gestation for uncomplicated human pregnancy. While results capture available data well, this study highlights significant gaps in available data required to better understand vascular remodeling in pregnancy.

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.

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.

Parameter estimation to study the immediate impact of aortic cross-clamping using reduced order models

Aortic cross-clamping is a common strategy during vascular surgery, however, its instantaneous impact on hemodynamics is unknown. We, therefore, developed two numerical models to estimate the immediate impact of aortic clamping on the vascular properties. To assess the validity of the models, we recorded continuous invasive pressure signals during abdominal aneurysm repair surgery, immediately before and after clamping. The first model is a zero-dimensional (0D) three-element Windkessel model, which we coupled to a gradient-based parameter estimation algorithm to identify patient-specific parameters such as vascular resistance and compliance. We found a 10% increase in the total resistance and a 20% decrease in the total compliance after clamping. The second model is a nine-artery network corresponding to an average human body in which we solved the one-dimensional (1D) blood flow equations. With a similar parameter estimation method and using the results from the 0D model, we identified the resistance boundary conditions of the 1D network. Determining the patient-specific total resistance and the distribution of peripheral resistances through the parameter estimation process was sufficient for the 1D model to accurately reproduce the impact of clamping on the pressure waveform. Both models gave an accurate description of the pressure wave and had a high correlation (R² > .95) with experimental blood pressure data.

Five years of cardio-ankle vascular index (CAVI) and CAVI0: how close are we to a pressure-independent index of arterial stiffness?

Pulse wave velocity, a common metric of arterial stiffness, is an established predictor for cardiovascular events and mortality. However, its intrinsic pressure-dependency complicates the discrimination of acute and chronic impacts of increased blood pressure on arterial stiffness. Cardio-ankle vascular index (CAVI) represented a significant step towards the development of a pressure-independent arterial stiffness metric. However, some potential limitations of CAVI might render this arterial stiffness metric less pressure-independent than originally thought. For this reason, we later introduced CAVI0. Nevertheless, advantages of one approach over the other are left debated. This review aims to shed light on the pressure (in)dependency of both CAVI and CAVI0. By critically reviewing results from studies reporting both CAVI and CAVI0 and using simple analytical methods, we show that CAVI0 may enhance the pressure-independent assessment of arterial stiffness, especially in the presence of large inter-individual differences in blood pressure.

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.

Comparison of Non-Newtonian Models of One-Dimensional Hemodynamics

The paper is devoted to the comparison of different one-dimensional models of blood flow. In such models, the non-Newtonian property of blood is considered. It is demonstrated that for the large arteries, the small parameter is observed in the models, and the perturbation method can be used for the analytical solution. In the paper, the simplified nonlinear problem for the semi-infinite vessel with constant properties is solved analytically, and the solutions for different models are compared. The effects of the flattening of the velocity profile and hematocrit value on the deviation from the Newtonian model are investigated.

A 1D-3D Hybrid Model of Patient-Specific Coronary Hemodynamics

PURPOSE

Coronary flow is affected by evolving events such as atherosclerotic plaque formation, rupture, and thrombosis, resulting in myocardial ischemia and infarction. Highly resolved 3D hemodynamic data at the stenosis is essential to model shear-sensitive thrombotic events in coronary artery disease.

METHODS

We developed a hybrid 1D-3D simulation framework to compute patient-specific coronary hemodynamics efficiently. A 1D model of the coronary flow is coupled to an image-based 3D model of the region of interest. This framework affords the advantages of reduced-order modeling, decreasing the global computational cost, without sacrificing the accuracy of the quantities of interest.

RESULTS

We validated our 1D-3D model against full 3D coronary simulations in healthy and diseased conditions. Our results showed good agreement between the 3D and the 1D-3D models while reducing the computational cost by 40-fold compared to the 3D simulation. The 1D-3D model predicted left/right coronary flow distribution within 3% and provided an accurate estimation of fractional flow reserve and wall shear stress distribution at the stenosis comparable to the 3D simulation.

CONCLUSION

Savings in computational cost may be significant in situations with changing geometry, such as growing thrombosis. Also, this approach would allow quantifying the time-dependent effect of thrombotic growth and occlusion on the global coronary circulation.

In silico trials for treatment of acute ischemic stroke: Design and implementation

An in silico trial simulates a disease and its corresponding therapies on a cohort of virtual patients to support the development and evaluation of medical devices, drugs, and treatment. In silico trials have the potential to refine, reduce cost, and partially replace current in vivo studies, namely clinical trials and animal testing. We present the design and implementation of an in silico trial for treatment of acute ischemic stroke. We propose an event-based modelling approach for the simulation of a disease and injury, where changes to the state of the system (the events) are assumed to be instantaneous. Using this approach we are able to combine a diverse set of models, spanning multiple time scales, to model acute ischemic stroke, treatment, and resulting brain tissue injury. The in silico trial is designed to be modular to aid development and reproducibility. It provides a comprehensive framework for application to any potential in silico trial. A statistical population model is used to generate cohorts of virtual patients. Patient functional outcomes are also predicted with a statistical model, using treatment and injury results and the patient's clinical parameters. We demonstrate the functionality of the event-based modelling approach and trial framework by running proof of concept in silico trials. The proof of concept trials simulate the same cohort of patients twice: once with successful treatment (successful recanalisation) and once with unsuccessful treatment (unsuccessful treatment). Ways to overcome some of the challenges and difficulties in setting up such an in silico trial are discussed, such as validation and computational limitations.

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.

Global sensitivity analysis for assessing the parameters importance and setting a stopping criterion in a biomedical inverse problem

This paper shows how to obtain in addition to the standard deviations available after a data assimilation procedure based on the ensemble Kalman filter, an apportioning of the total uncertainty in the outputs of a patient-specific blood flow model into small portions of uncertainty due to input parameters. Statistical indicators generally used for identifying the importance of numerical parameters, namely the Sobol' first order and total indices, are introduced and discussed. These allow the identification of the importance rank of the different input parameters for the patient-specific blood flow model, as well as the influence of the interactions between these parameters on the model output variance. The results show that knowing the importance rank of the model input parameters during the assimilation procedure is useful to avoid unnecessary over-solving and to find a suitable stopping criterion in clinical situations where faster diagnosis is always requested. Indeed, the work permits to reduce typically by a factor of six the time to solution and most importantly with very limited extra calculation using already available information.

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.

Uncertainty quantification implementations in human hemodynamic flows

BACKGROUND AND OBJECTIVE

Human hemodynamic modeling is usually influenced by uncertainties occurring from a considerable unavailability of information linked to the boundary conditions and the physical properties used in the numerical models. Calculating the effect of these uncertainties on the numerical findings along the cardiovascular system is a demanding process due to the complexity of the morphology of the body and the area dynamics. To cope with all these difficulties, Uncertainty Quantification (UQ) methods seem to be an ideal tool.

RESULTS

This study focuses on analyzing and summarizing some of the recent research efforts and directions of implementing UQ in human hemodynamic flows by analyzing 139 research papers. Initially, the suitability of applying this approach is analyzed and demonstrated. Then, an overview of the most significant research work in various fields of biomedical hemodynamic engineering is presented. Finally, it is attempted to identify any possible forthcoming directions for research and methodological progress of UQ in biomedical sciences.

CONCLUSION

This review concludes that by finding the best statistical methods and parameters to represent the propagated uncertainties, while achieving a good interpretation of the interaction between input-output, is crucial for implementing UQ in biomedical sciences.

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.

Modelling the leptomeningeal collateral circulation during acute ischaemic stroke

A novel model of the leptomeningeal collateral circulation is created by combining data from multiple sources with statistical scaling laws. The extent of the collateral circulation is varied by defining a collateral vessel probability. Blood flow and pressure are simulated using a one-dimensional steady state blood flow model. The leptomeningeal collateral vessels provide significant flow during a stroke. The pressure drop over an occlusion predicted by the model ranges between 60 and 85 mmHg depending on the extent of the collateral circulation. The linear transport of contrast material was simulated in the circulatory network. The time delay of peak contrast over an occlusion is 3.3 s in the model, and 2.1 s (IQR 0.8-4.0 s) when measured in dynamic CTA data of acute ischaemic stroke patients. Modelling the leptomeningeal collateral circulation could lead to better estimates of infarct volume and patient outcome.

Complementing sparse vascular imaging data by physiological adaptation rules

Mathematical modeling of pressure and flow waveforms in blood vessels using pulse wave propagation (PWP) models has tremendous potential to support clinical decision making. For a personalized model outcome, measurements of all modeled vessel radii and wall thicknesses are required. In clinical practice, however, data sets are often incomplete. To overcome this problem, we hypothesized that the adaptive capacity of vessels in response to mechanical load could be utilized to fill in the gaps of incomplete patient-specific data sets. We implemented homeostatic feedback loops in a validated PWP model to allow adaptation of vessel geometry to maintain physiological values of wall stress and wall shear stress. To evaluate our approach, we gathered vascular MRI and ultrasound data sets of wall thicknesses and radii of central and arm arterial segments of 10 healthy subjects. Reference models (i.e., termed RefModel, n = 10) were simulated using complete data, whereas adapted models (AdaptModel, n = 10) used data of one carotid artery segment only, and the remaining geometries in this model were estimated using adaptation. We evaluated agreement between RefModel and AdaptModel geometries, as well as that between pressure and flow waveforms of both models. Limits of agreement (bias ± 2 SD of difference) between AdaptModel and RefModel radii and wall thicknesses were 0.2 ± 2.6 mm and -140 ± 557 µm, respectively. Pressure and flow waveform characteristics of the AdaptModel better resembled those of the RefModels as compared with the model in which the vessels were not adapted. Our adaptation-based PWP model enables personalization of vascular geometries even when not all required data are available.NEW & NOTEWORTHY To benefit personalized pulse wave propagation (PWP) modeling, we propose a novel method that, instead of relying on extensive data sets on vascular geometries, incorporates physiological adaptation rules. The developed vascular adaptation model adequately predicted arterial radius and wall thickness compared with ultrasound and MRI estimates, obtained in humans. Our approach could be used as a tool to facilitate personalized modeling, notably in case of missing data, as routinely found in clinical settings.

Coupling one-dimensional arterial blood flow to three-dimensional tissue perfusion models for in silico trials of acute ischaemic stroke

An acute ischaemic stroke is due to the sudden blockage of an intracranial blood vessel by an embolized thrombus. In the context of setting up in silico trials for the treatment of acute ischaemic stroke, the effect of a stroke on perfusion and metabolism of brain tissue should be modelled to predict final infarcted brain tissue. This requires coupling of blood flow and tissue perfusion models. A one-dimensional intracranial blood flow model and a method to couple this to a brain tissue perfusion model for patient-specific simulations is presented. Image-based patient-specific data on the anatomy of the circle of Willis are combined with literature data and models for vessel anatomy not visible in the images, to create an extended model for each patient from the larger vessels down to the pial surface. The coupling between arterial blood flow and tissue perfusion occurs at the pial surface through the estimation of perfusion territories. The coupling method is able to accurately estimate perfusion territories. Finally, we argue that blood flow can be approximated as steady-state flow at the interface between arterial blood flow and tissue perfusion to reduce the cost of organ-scale simulations.

In vitro Biomodels in Stenotic Arteries to Perform Blood Analogues Flow Visualizations and Measurements: A Review

Cardiovascular diseases are one of the leading causes of death globally and the most common pathological process is atherosclerosis. Over the years, these cardiovascular complications have been extensively studied by applying in vivo, in vitro and numerical methods (in silico). In vivo studies represent more accurately the physiological conditions and provide the most realistic data. Nevertheless, these approaches are expensive, and it is complex to control several physiological variables. Hence, the continuous effort to find reliable alternative methods has been growing. In the last decades, numerical simulations have been widely used to assess the blood flow behavior in stenotic arteries and, consequently, providing insights into the cardiovascular disease condition, its progression and therapeutic optimization. However, it is necessary to ensure its accuracy and reliability by comparing the numerical simulations with clinical and experimental data. For this reason, with the progress of the in vitro flow measurement techniques and rapid prototyping, experimental investigation of hemodynamics has gained widespread attention. The present work reviews state-of-the-art in vitro macro-scale arterial stenotic biomodels for flow measurements, summarizing the different fabrication methods, blood analogues and highlighting advantages and limitations of the most used techniques.

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.

Towards enabling a cardiovascular digital twin for human systemic circulation using inverse analysis

An exponential rise in patient data provides an excellent opportunity to improve the existing health care infrastructure. In the present work, a method to enable cardiovascular digital twin is proposed using inverse analysis. Conventionally, accurate analytical solutions for inverse analysis in linear problems have been proposed and used. However, these methods fail or are not efficient for nonlinear systems, such as blood flow in the cardiovascular system (systemic circulation) that involves high degree of nonlinearity. To address this, a methodology for inverse analysis using recurrent neural network for the cardiovascular system is proposed in this work, using a virtual patient database. Blood pressure waveforms in various vessels of the body are inversely calculated with the help of long short-term memory (LSTM) cells by inputting pressure waveforms from three non-invasively accessible blood vessels (carotid, femoral and brachial arteries). The inverse analysis system built this way is applied to the detection of abdominal aortic aneurysm (AAA) and its severity using neural networks.

Experimental validation of computational fluid dynamics for solving isothermal and incompressible viscous fluid flow

In order to validate a computational method for solving viscous fluid flows, experiments are carried out in an eccentric cylindrical cavity showing various flow formations over a range of Reynolds numbers. Especially, in numerical solution approaches for isothermal and incompressible flows, we search for simple experimental data for evaluating accuracy as well as performance of the computational method. Verification of different computational methods is arduous, and analytic solutions are only obtained for simple geometries like a channel flow. Clearly, a method is expected to predict different flow patterns within a cavity. Thus, we propose a configuration generating different flow formations depending on the Reynolds number and make the experimental results freely available in order to be used as an assessment criterion to demonstrate the reliability of a new computational approach.

Mechanics of the dicrotic notch: An acceleration hypothesis

The dicrotic notch is a prominent and distinctive feature of the pressure waveform in the central arteries. It is universally used to demarcate the end of systole and the beginning of diastole in these arteries. Despite its importance clinically, no physical mechanism for the formation of the dicrotic notch has been demonstrated convincingly. We first explore a mechanism based on the reflection of a backward wavefront from the aortic valve at the time of closure. This hypothesis is rejected on the basis of experimental evidence from measurements made in dogs. A new hypothesis is presented involving the acceleration of the aortic valve apparatus at the time of valve closure. This hypothesis is supported by new calculations of the acceleration of the aortic valve apparatus during the cardiac cycle based on computed tomography scans in man.

Cardiovascular models for personalised medicine: Where now and where next?

The aim of this position paper is to provide a brief overview of the current status of cardiovascular modelling and of the processes required and some of the challenges to be addressed to see wider exploitation in both personal health management and clinical practice. In most branches of engineering the concept of the digital twin, informed by extensive and continuous monitoring and coupled with robust data assimilation and simulation techniques, is gaining traction: the Gartner Group listed it as one of the top ten digital trends in 2018. The cardiovascular modelling community is starting to develop a much more systematic approach to the combination of physics, mathematics, control theory, artificial intelligence, machine learning, computer science and advanced engineering methodology, as well as working more closely with the clinical community to better understand and exploit physiological measurements, and indeed to develop jointly better measurement protocols informed by model-based understanding. Developments in physiological modelling, model personalisation, model outcome uncertainty, and the role of models in clinical decision support are addressed and 'where-next' steps and challenges discussed.