Outflow boundary conditions for blood flow in arterial trees

Differential sensitivities to blood pressure variations in internal carotid and intracranial arteries: a numerical approach to stroke prediction

Stroke remains a global health concern, necessitating early prediction for effective management. Atherosclerosis-induced internal carotid and intra cranial stenosis contributes significantly to stroke risk. This study explores the relationship between blood pressure and stroke prediction, focusing on internal carotid artery (ICA) branches: middle cerebral artery (MCA), anterior cerebral artery (ACA), and their role in hemodynamics. Computational fluid dynamics (CFD) informed by the Windkessel model were employed to simulate patient-specific ICA models with introduced stenosis. Central to our investigation is the impact of stenosis on blood pressure, flow velocity, and flow rate across these branches, incorporating Fractional Flow Reserve (FFR) analysis. Results highlight differential sensitivities to blood pressure variations, with M1 branch showing high sensitivity, ACA moderate, and M2 minimal. Comparing blood pressure fluctuations between ICA and MCA revealed heightened sensitivity to potential reverse flow compared to ICA and ACA comparisons, emphasizing MCA's role. Blood flow adjustments due to stenosis demonstrated intricate compensatory mechanisms. FFR emerged as a robust predictor of stenosis severity, particularly in the M2 branch. In conclusion, this study provides comprehensive insights into hemodynamic complexities within major intracranial arteries, elucidating the significance of blood pressure variations, flow attributes, and FFR in stenosis contexts. Subject-specific data integration enhances model reliability, aiding stroke risk assessment and advancing cerebrovascular disease understanding.

An analytical method informed by clinical imaging data for estimating outlet boundary conditions in computational fluid dynamics analysis of carotid artery blood flow

Stroke occur mainly due to arterial thrombosis and rupture of cerebral blood vessels. Previous studies showed that blood flow-induced wall shear stress is an essential bio marker for estimating atherogenesis. It is a common practice to use computational fluid dynamics (CFD) simulations to calculate wall shear stress and to quantify blood flow. Reliability of predicted CFD results greatly depends on the accuracy of applied boundary conditions. Previously, the boundary conditions were estimated by varying values so that they matched the clinical data. It is applicable upon the availability of clinical data. Meanwhile, in most cases all that can be accessed are arterial geometry and inflow rate. Consequently, there is a need to devise a tool to estimate boundary values such as resistance and compliance of arteries. This study proposes an analytical framework to estimate the boundary conditions for a carotid artery based on the geometries of the downstream arteries available from clinical images.

Central arterial pressure and patient-specific model parameter estimation based on radial pressure measurements

One-dimensional arterial flow simulations are suitable to estimate the aortic pressure from peripheral measurements in a patient-specific arterial network. This study introduces a reduction of the system parameters, and a novel calculation method to estimate the patient-specific set and the aortic curve based on radial applanation tonometry. Peripheral and aortic pressure curves were measured in patients, optimization were carried out. The aortic pressure curves were reproduced well, with an overestimation of the measured Systolic and Mean Blood Pressures on average by 0.6 and 0.5 mmHg respectively, and the Root Mean Square Difference of the curves was 3 mmHg on average.

First blood: An efficient, hybrid one- and zero-dimensional, modular hemodynamic solver

Low-dimensional (1D or 0D) models can describe the whole human blood circulation, for example, 1D distributed parameter model for the arterial network and 0D concentrated models for the heart or other organs. This paper presents a combined 1D-0D solver, called first_blood, that solves the governing equations of fluid dynamics to model low-dimensional hemodynamic effects. An extended method of characteristics is applied here to solve the momentum, and mass conservation equations and the viscoelastic wall model equation, mimicking the material properties of arterial walls. The heart and the peripheral lumped models are solved with a general zero-dimensional (0D) nonlinear solver. The model topology can be modular, that is, first_blood can solve any 1D-0D hemodynamic model. To demonstrate the applicability of first_blood, the human arterial system, the heart and the peripherals are modelled using the solver. The simulation time of a heartbeat takes around 2 s, that is, first_blood requires only twice the real-time for the simulation using an average PC, which highlights the computational efficiency. The source code is available on GitHub, that is, it is open source. The model parameters are based on the literature suggestions and on the validation of output data to obtain physiologically relevant results.

Wall shear stress indicators influence the regular hemodynamic conditions in coronary main arterial diseases: cardiovascular abnormalities

Computational hemodynamic (CH) characteristics play a central role in the onset and expansion of atherosclerotic plaques in the coronary main arteries. This study has explored the effects of hemodynamic properties especially coronary arterial wall tangential stresses on various healthy and diseased patient-based coronary artery models based on coronary computed tomography angiography (CCTA) imaging. The key components of the work are the CCTA image acquisition, accurate three-dimensional (3 D) model segmentation, reconstruction, appropriate grid generation, CH simulations, and analysis of the results by using open-source techniques. The CH simulation results have produced hemodynamic variables, including velocity magnitude (VM), mean arterial pressure difference, wall shear stress (WSS), time-averaged WSS (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and finally, computational fractional flow reserve (cFFR), that allow the pathophysiological conditions in patient-based coronary models. The VM, mean pressure difference, and WSS indices have yielded consistent simulation results for predicting the severity conditions of coronary diseases. We have compared our cFFR results with the published results and observed that the WSS indices were a good alternative approach for measuring the severity of coronary lesions. The CH results allow a medical expert to estimate the severity of a lumen area and stenosis physiological blood flow conditions in a non-invasive way.

Parametric analysis of an efficient boundary condition to control outlet flow rates in large arterial networks

Substantial effort is being invested in the creation of a virtual human-a model which will improve our understanding of human physiology and diseases and assist clinicians in the design of personalised medical treatments. A central challenge of achieving blood flow simulations at full-human scale is the development of an efficient and accurate approach to imposing boundary conditions on many outlets. A previous study proposed an efficient method for implementing the two-element Windkessel model to control the flow rate ratios at outlets. Here we clarify the general role of the resistance and capacitance in this approach and conduct a parametric sweep to examine how to choose their values for complex geometries. We show that the error of the flow rate ratios decreases exponentially as the resistance increases. The errors fall below 4% in a simple five-outlets model and 7% in a human artery model comprising ten outlets. Moreover, the flow rate ratios converge faster and suffer from weaker fluctuations as the capacitance decreases. Our findings also establish constraints on the parameters controlling the numerical stability of the simulations. The findings from this work are directly applicable to larger and more complex vascular domains encountered at full-human scale.

In Vitro Flow Chamber Design for the Study of Endothelial Cell (Patho)Physiology

In the native vasculature, flowing blood produces a frictional force on vessel walls that affects endothelial cell function and phenotype. In the arterial system, the vasculature's local geometry directly influences variations in flow profiles and shear stress magnitudes. Straight arterial sections with pulsatile shear stress have been shown to promote an athero-protective endothelial phenotype. Conversely, areas with more complex geometry, such as arterial bifurcations and branch points with disturbed flow patterns and lower, oscillatory shear stress, typically lead to endothelial dysfunction and the pathogenesis of cardiovascular diseases. Many studies have investigated the regulation of endothelial responses to various shear stress environments. Importantly, the accurate in vitro simulation of in vivo hemodynamics is critical to the deeper understanding of mechanotransduction through the proper design and use of flow chamber devices. In this review, we describe several flow chamber apparatuses and their fluid mechanics design parameters, including parallel-plate flow chambers, cone-and-plate devices, and microfluidic devices. In addition, chamber-specific design criteria and relevant equations are defined in detail for the accurate simulation of shear stress environments to study endothelial cell responses.

Evaluation and verification of patient-specific modelling of type B aortic dissection

Quantitative assessment of the complex hemodynamic environment in type B aortic dissection (TBAD) through computational fluid dynamics (CFD) simulations can provide detailed insights into the disease and its progression. As imaging and computational technologies have advanced, methodologies have been developed to increase the accuracy and physiological relevance of CFD simulations. This study presents a patient-specific workflow to simulate blood flow in TBAD, utilising the maximum amount of in vivo data available in the form of CT images, 4D-flow MRI and invasive Doppler-wire pressure measurements, to implement the recommended current best practice methodologies in terms of patient-specific geometry and boundary conditions. The study aimed to evaluate and verify this workflow through detailed qualitative and quantitative comparisons of the CFD and in vivo data. Based on data acquired from five TBAD patients, a range of essential model inputs was obtained, including inlet flow waveforms and 3-element Windkessel model parameters, which can be utilised in further studies where in vivo flow data is not available. Local and global analysis showed good consistency between CFD results and 4D-MRI data, with the maximum velocity in the primary entry tear differing by up to 0.3 m/s, and 80% of the analysed regions achieving moderate or strong correlations between the predicted and in vivo velocities. CFD predicted pressures were generally well matched to the Doppler-wire measurements, with some deviation in peak systolic values. Overall, this study presents a validated comprehensive workflow with extensive data for CFD simulation of TBAD.

Self-Similar Functional Circuit Models of Arteries and Deterministic Fractal Operators: Theoretical Revelation for Biomimetic Materials

In this paper, the self-similar functional circuit models of arteries are proposed for bioinspired hemodynamic materials design. Based on the mechanical-electrical analogous method, the circuit model can be utilized to mimic the blood flow of arteries. The theoretical mechanism to quantitatively simulate realistic blood flow is developed by establishing a fractal circuit network with an infinite number of electrical components. We have found that the fractal admittance operator obtained from the minimum repeating unit of the fractal circuit can simply and directly determine the blood-flow regulation mechanism. Furthermore, according to the operator algebra, the fractal admittance operator on the aorta can be represented by Gaussian-type convolution kernel function. Similarly, the arteriolar operator can be described by Bessel-type function. Moreover, by the self-similar assembly pattern of the proposed model, biomimetic materials which contain self-similar circuits can be designed to mimic physiological or pathological states of blood flow. Studies show that the self-similar functional circuit model can efficiently describe the blood flow and provide an available and convenient structural theoretical revelation for the preparation of in vitro hemodynamic bionic materials.

On the relevance of boundary conditions and viscosity models in blood flow simulations in patient-specific aorto-coronary bypass models

Physiologically realistic results are the aim of every blood flow simulation. This is not different in aorto-coronary bypasses where the properties of the coronary circulation may significantly affect the relevance of the performed simulations. By considering three patient-specific bypass geometries, the present article focuses on two aspects of the coronary blood flow - its phasic flow pattern and its behaviour affected by blood rheology. For the phasic flow property, a multiscale modelling approach is chosen as a means to assess the ability of five different types of coronary boundary conditions (mean arterial pressure, Windkessel model and three lumped parameter models) to attain realistic coronary haemodynamics. From the analysed variants of boundary conditions, the best option in terms of physiological characteristics and its potential for use in patient-based simulations, is utilised to account for the effect of shear-dependent viscosity on the resulting haemodynamics and wall shear stress stimulation. Aside from the Newtonian model, the blood rheology is approximated by two non-Newtonian models in order to determine whether the choice of a viscosity model is important in simulations involving coronary circulation. A comprehensive analysis of obtained results demonstrated notable superiority of all lumped parameter models, especially in comparison to the constant outlet pressure, which regardless of bypass type gave overestimated and physiologically misleading results. In terms of rheology, it was noted that blood in undamaged coronary arteries behaves as a Newtonian fluid, whereas in vessels with atypical lumen geometry, such as that of anastomosis or stenosis, its shear-thinning behaviour should not be ignored.

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.

Application of 0D model of blood flow to study vessel abnormalities in the human systemic circulation: An in-silico study

In this paper, a multi-compartment 0D model of the blood flow is considered to study the vessel abnormalities (stenoses and aneurysms) in the human systemic circulation (SC). In the complete SC, different levels of stenosis and aneurysms are artificially created by decreasing and increasing the vessel diameters respectively and their effects on pressure and flow are studied using sensitivity analysis (SA). Normalized local sensitivity analysis (LSA) is used to study the impact of stenosis and aneurysms on pressure and flow wave pattern. Furthermore, global sensitivity analysis (GSA), Sobol’s method is used to quantify the overall influence of stenoses and aneurysms in the complete SC.

The results of global sensitivity analysis revealed that the impact of both stenoses and aneurysms is strong within the individual structures (arm, legs, carotid bifurcation, aorta), while, aortic stenoses and aneurysms have effect on almost all downstream nodes. Moreover, the study could be useful for medical doctors, teachers and students to observe the hemodynamical changes in the SC with respect to vessel abnormalities, which could further help in making any clinical decision for patients having different levels of vessel abnormalities in any part or structure of the SC.

Computational Fluid Dynamics and Aortic Dissections: Panacea or Panic?

This paper reviews the methodology, benefits and limitations associated with computational flow dynamics (CFD) in the field of vascular surgery. Combined with traditional imaging of the vasculature, CFD simulation enables accurate characterisation of real-time physiological and haemodynamic parameters such as wall shear stress. This enables vascular surgeons to understand haemodynamic changes in true and false lumens, and exit and re-entry tears. This crucial information may facilitate triaging decisions. Furthermore, CFD can be used to assess the impact of stent graft treatment, as it provides a haemodynamic account of what may cause procedure-related complications. Efforts to integrate conventional imaging, individual patient data and CFD are paramount to its success, given its potential to replace traditional registry-based, population-averaged data. Nonetheless, methodological limitations must be addressed before clinical implementation. This must be accompanied by further research with large sample sizes, to establish the association between haemodynamic patterns as observed by CFD and progression of aortic dissection.

Data assimilation for identification of cardiovascular network characteristics

A method to estimate the hemodynamics parameters of a network of vessels using an Ensemble Kalman filter is presented. The elastic moduli (Young's modulus) of blood vessels and the terminal boundary parameters are estimated as the solution of an inverse problem. Two synthetic test cases and a configuration where experimental data are available are presented. The sensitivity analysis confirms that the proposed method is quite robust even with a few numbers of observations. The simulations with the estimated parameters recovers target pressure or flow rate waveforms at given specific locations, improving the state-of-the-art predictions available in the literature. This shows the effectiveness and efficiency of both the parameter estimation algorithm and the blood flow model.

Patient-Specific, Multi-Scale Modeling of Neointimal Hyperplasia in Vein Grafts

Neointimal hyperplasia is amongst the major causes of failure of bypass grafts. The disease progression varies from patient to patient due to a range of different factors. In this paper, a mathematical model will be used to understand neointimal hyperplasia in individual patients, combining information from biological experiments and patient-specific data to analyze some aspects of the disease, particularly with regard to mechanical stimuli due to shear stresses on the vessel wall. By combining a biochemical model of cell growth and a patient-specific computational fluid dynamics analysis of blood flow in the lumen, remodeling of the blood vessel is studied by means of a novel computational framework. The framework was used to analyze two vein graft bypasses from one patient: a femoro-popliteal and a femoro-distal bypass. The remodeling of the vessel wall and analysis of the flow for each case was then compared to clinical data and discussed as a potential tool for a better understanding of the disease. Simulation results from this first computational approach showed an overall agreement on the locations of hyperplasia in these patients and demonstrated the potential of using new integrative modeling tools to understand disease progression.

Impact of Using Conventional Inlet/Outlet Boundary Conditions on Haemodynamic Metrics in a Subject-Specific Rabbit Aorta

Computational fluid dynamics is a tool capable of accurately measuring metrics currently used to predict the behaviour of cardiovascular diseases. This study quantifies the impact various commonly used inlet and outlet boundary conditions have on various shear rate–based haemodynamic metrics currently used for predicting the localisation of cardiovascular diseases. Simulations are conducted on an accurately represented rabbit aorta configuration and comparison has been made against available in vivo data. The boundary conditions studied include two different inlet profiles, three outlet boundary conditions, and steady-state versus pulsatile flow cases. Large variations were found in the results, particularly when using different outlet boundary conditions, and the discrepancies were evaluated both quantitatively and qualitatively. The results clearly highlight the importance of the type of boundary condition used when simulating complex cardiovascular models. By restricting the attention to the flow within the aorta and the intercostal branches, the results suggest that prescribing transient simulation and fully developed flow at the inlet are not required. Furthermore, assuming the widely accepted low wall shear stress theory of Caro, it was found that Murray’s law–based outlet boundary condition returns the most physiologically accurate results when compared to in vivo data.

Uncertainty Quantification in a Patient-Specific One-Dimensional Arterial Network Model: EnKF-Based Inflow Estimator

Successful clinical use of patient-specific models for cardiovascular dynamics depends on the reliability of the model output in the presence of input uncertainties. For 1D fluid dynamics models of arterial networks, input uncertainties associated with the model output are related to the specification of vessel and network geometry, parameters within the fluid and wall equations, and parameters used to specify inlet and outlet boundary conditions. This study investigates how uncertainty in the flow profile applied at the inlet boundary of a 1D model affects area and pressure predictions at the center of a single vessel. More specifically, this study develops an iterative scheme based on the ensemble Kalman filter (EnKF) to estimate the temporal inflow profile from a prior distribution of curves. The EnKF-based inflow estimator provides a measure of uncertainty in the size and shape of the estimated inflow, which is propagated through the model to determine the corresponding uncertainty in model predictions of area and pressure. Model predictions are compared to ex vivo area and blood pressure measurements in the ascending aorta, the carotid artery, and the femoral artery of a healthy male Merino sheep. Results discuss dynamics obtained using a linear and a nonlinear viscoelastic wall model.

Patient-specific structural effects on hemodynamics in the ischemic lower limb artery

Lower limb peripheral artery disease is a prevalent chronic non-communicable disease without obvious symptoms. However, the effect of ischemic lower limb peripheral arteries on hemodynamics remains unclear. In this study, we investigated the variation of the hemodynamics caused by patient-specific structural artery characteristics. Computational fluid dynamic simulations were performed on seven lower limb (including superficial femoral, deep femoral and popliteal) artery models that were reconstructed from magnetic resonance imaging. We found that increased wall shear stress (WSS) was mainly caused by the increasing severity of stenosis, bending, and branching. Our results showed that the increase in the WSS value at a stenosis at the bifurcation was 2.7 Pa. In contrast, the isolated stenosis and branch caused a WSS increase of 0.7 Pa and 0.5 Pa, respectively. The WSS in the narrow popliteal artery was more sensitive to a reduction in radius. Our results also demonstrate that the distribution of the velocity and pressure gradient are highly structurally related. At last, Ultrasound Doppler velocimeter measured result was presented as a validation. In conclusion, the distribution of hemodynamics may serve as a supplement for clinical decision-making to prevent the occurrence of a morbid or mortal ischemic event.

Non-dimensional physics of pulsatile cardiovascular networks and energy efficiency

In Nature, there exist a variety of cardiovascular circulation networks in which the energetic ventricular load has both steady and pulsatile components. Steady load is related to the mean cardiac output (CO) and the haemodynamic resistance of the peripheral vascular system. On the other hand, the pulsatile load is determined by the simultaneous pressure and flow waveforms at the ventricular outlet, which in turn are governed through arterial wave dynamics (transmission) and pulse decay characteristics (windkessel effect). Both the steady and pulsatile contributions of the haemodynamic power load are critical for characterizing/comparing disease states and for predicting the performance of cardiovascular devices. However, haemodynamic performance parameters vary significantly from subject to subject because of body size, heart rate and subject-specific CO. Therefore, a 'normalized' energy dissipation index, as a function of the 'non-dimensional' physical parameters that govern the circulation networks, is needed for comparative/integrative biological studies and clinical decision-making. In this paper, a complete network-independent non-dimensional formulation that incorporates pulsatile flow regimes is developed. Mechanical design variables of cardiovascular flow systems are identified and the Buckingham Pi theorem is formally applied to obtain the corresponding non-dimensional scaling parameter sets. Two scaling approaches are considered to address both the lumped parameter networks and the distributed circulation components. The validity of these non-dimensional number sets is tested extensively through the existing empirical allometric scaling laws of circulation systems. Additional validation studies are performed using a parametric numerical arterial model that represents the transmission and windkessel characteristics, which are adjusted to represent different body sizes and non-dimensional haemodynamic states. Simulations demonstrate that the proposed non-dimensional indices are independent of body size for healthy conditions, but are sensitive to deviations caused by off-design disease states that alter the energetic load. Sensitivity simulations are used to identify the relationship between pulsatile power loss and non-dimensional characteristics, and optimal operational states are computed.