Optimizing cerebral perfusion and hemodynamics during cardiopulmonary bypass through cannula design combining in silico, in vitro and in vivo input.

Cardiopulmonary bypass (CPB) is a standard technique for cardiac surgery, but comes with the risk of severe neurological complications (e.g. stroke) caused by embolisms and/or reduced cerebral perfusion. We report on an aortic cannula prototype design (optiCAN) with helical outflow and jet-splitting dispersion tip that could reduce the risk of embolic events and restores cerebral perfusion to 97.5% of physiological flow during CPB in vivo, whereas a commercial curved-tip cannula yields 74.6%. In further in vitro comparison, pressure loss and hemolysis parameters of optiCAN remain unaffected. Results are reproducibly confirmed in silico for an exemplary human aortic anatomy via computational fluid dynamics (CFD) simulations. Based on CFD simulations, we firstly show that optiCAN design improves aortic root washout, which reduces the risk of thromboembolism. Secondly, we identify regions of the aortic intima with increased risk of plaque release by correlating areas of enhanced plaque growth and high wall shear stresses (WSS). From this we propose another easy-to-manufacture cannula design (opti2CAN) that decreases areas burdened by high WSS, while preserving physiological cerebral flow and favorable hemodynamics. With this novel cannula design, we propose a cannulation option to reduce neurological complications and the prevalence of stroke in high-risk patients after CPB.

Tipless transseptal cannula concept combines improved hemodynamic properties and risk-reduced placement: An in silico proof-of-concept.

As a leading cause of death worldwide, heart failure is a serious medical condition in which many critically ill patients require temporary mechanical circulatory support (MCS) as a bridge-to-recovery or bridge-to-decision. In many cases, the TandemHeart system is used to unload the left heart by draining blood from the left atrium (LA) to the femoral artery via a transseptal multistage cannula. However, even though the correct positioning of the cannula is crucial for a safe treatment, the long cannula tip currently used in transseptal cannulas complicates positioning, making the cannula vulnerable to displacement during MCS. To overcome these limitations, we propose the development of a new tipless transseptal cannula with improved hemodynamic properties. We discuss the tipless cannula concept by comparing it to the common multistage cannula concept using computational fluid dynamics simulations and assess the flow field in the LA, the wall shear stresses (WSS), and the pressure loss. Across the two distinct time points of end-systole and end-diastole and two drainage flow rates of 3.5 and 5.0 L/min, we find a more homogeneous inlet flow pattern for the tipless cannula concept, accompanied by a remarkably reduced area of platelet-activating WSS (up to 10-times smaller area compared to the multistage cannula). Moreover, pressure loss is up to 14.5% lower in the tipless cannula concept, confirming overall improved hemodynamic properties of the tipless cannula concept. Finally, a diameter-dependent study reveals that lower WSS and pressure losses can be further reduced by large-lumen designs for any simulation setting. Overall, our results suggest that a tipless cannula concept remedies the crucial disadvantages of a long-tip multistage cannula by reducing the risk of misplacement, and it furthermore promotes optimized hemodynamics. With this successful proof-of-concept, we underscore the potential for and encourage the realization of further experimental investigations regarding the development of a tipless transseptal cannula for MCS.

Crucial Aspects for Using Computational Fluid Dynamics as a Predictive Evaluation Tool for Blood Pumps.

The suitability of computational fluid dynamics (CFD) as a regulatory tool for safety assessment of medical devices is still limited: A lack of standardized validation and evaluation methods impairs the quantitative comparability and reliability of simulation studies, particularly regarding the assessment of hemocompatibility. This study investigated important aspects of validation and verification for three common turbulence modeling approaches (laminar, k-ω shear stress transport [SST] and stress-blended eddy simulation [SBES]) and three different mesh refinements. Simulation results for pressure head, characteristic velocity, and shear stress for the benchmark blood pump model of the Food and Drug Administration critical path initiative were compared with its published experimental results. For the highest mesh resolution, all three models predicted the hydraulic pump characteristics with a relative deviation averaged over six operating conditions below 6.1%. In addition, the SBES model showed an accurate agreement of the characteristic velocity field in the pump's diffusor region (relative error <2.9%), while the laminar and SST model calculated significantly elevated and deviating velocity amplitudes (>43.6%). The ability to quantify shear stress is fundamental for the prediction of blood damage. In this respect, this study demonstrated that: 1) a close agreement and validation of both pressure head and characteristic velocity was feasible and 2) the shear stress quantification demanded higher near-wall mesh resolutions, although such high resolutions were not required for the validation of only pressure heads or velocity. Hence, a mesh verification analysis for shear stresses may prove significant for the development of credible CFD blood damage predictions in the future.

Mesh Sensitivity Analysis for Quantitative Shear Stress Assessment in Blood Pumps Using Computational Fluid Dynamics.

The reduction of excessive, nonphysiologic shear stresses leading to blood trauma can be the key to overcome many of the associated complications in blood recirculating devices. In that regard, computational fluid dynamics (CFD) are gaining in importance for the hydraulic and hemocompatibility assessment. Still, direct hemolysis assessments with CFD remain inaccurate and limited to qualitative comparisons rather than quantitative predictions. An underestimated quantity for improved blood damage prediction accuracy is the influence of near-wall mesh resolution on shear stress quantification in regions of complex flows. This study investigated the necessary mesh refinement to quantify shear stress for two selected, meshing sensitive hotspots within a rotary centrifugal blood pump (the blade leading edge and tip clearance gap). The shear stress in these regions is elevated due to presence of stagnation points and the flow around a sharp edge. The nondimensional mesh characteristic number y+, which is known in the context of turbulence modeling, underestimated the maximum wall shear stress by 60% on average with the recommended value of 1, but was found to be exact below 0.1. To evaluate the meshing related error on the numerical hemolysis prediction, three-dimensional simulations of a generic centrifugal pump were performed with mesh sizes from 3 × 106 to 30 × 106 elements. The respective hemolysis was calculated using an Eulerian scalar transport model. Mesh insensitivity was found below a maximum y+ of 0.2 necessitating 18 × 106 mesh elements. A meshing related error of up to 25% was found for the coarser meshes. Further investigations need to address: (1) the transferability to other geometries and (2) potential adaptions on blood damage estimation models to allow better quantitative predictions.

Fluid-structure interaction of a pulsatile flow with an aortic valve model: A combined experimental and numerical study.

The complex fluid-structure interaction problem associated with the flow of blood through a heart valve with flexible leaflets is investigated both experimentally and numerically. In the experimental test rig, a pulse duplicator generates a pulsatile flow through a biomimetic rigid aortic root where a model of aortic valve with polymer flexible leaflets is implanted. High-speed recordings of the leaflets motion and particle image velocimetry measurements were performed together to investigate the valve kinematics and the dynamics of the flow. Large eddy simulations of the same configuration, based on a variant of the immersed boundary method, are also presented. A massively parallel unstructured finite-volume flow solver is coupled with a finite-element solid mechanics solver to predict the fluid-structure interaction between the unsteady flow and the valve. Detailed analysis of the dynamics of opening and closure of the valve are conducted, showing a good quantitative agreement between the experiment and the simulation regarding the global behavior, in spite of some differences regarding the individual dynamics of the valve leaflets. A multicycle analysis (over more than 20 cycles) enables to characterize the generation of turbulence downstream of the valve, showing similar flow features between the experiment and the simulation. The flow transitions to turbulence after peak systole, when the flow starts to decelerate. Fluctuations are observed in the wake of the valve, with maximum amplitude observed at the commissure side of the aorta. Overall, a very promising experiment-vs-simulation comparison is shown, demonstrating the potential of the numerical method.

Hemodynamic analysis of outflow grafting positions of a ventricular assist device using closed-loop multiscale CFD simulations: Preliminary results.

Subclavian arteries are a possible alternate location for left ventricular assist device (LVAD) outflow grafts due to easier surgical access and application in high risk patients. As vascular blood flow mechanics strongly influence the clinical outcome, insights into the hemodynamics during LVAD support can be used to evaluate different grafting locations. In this study, the feasibility of left and right subclavian artery (SA) grafting was investigated for the HeartWare HVAD with a numerical multiscale model. A 3-D CFD model of the aortic arch was coupled to a lumped parameter model of the cardiovascular system under LVAD support. Grafts in the left and right SA were placed at three different anastomoses angles (90°, 60° and 30°). Additionally, standard grafting of the ascending and descending aorta was modelled. Full support LVAD (5l/min) and partial support LVAD (3l/min) in co-pulsation and counter-pulsation mode were analysed. The grafting positions were investigated regarding coronary and cerebral perfusion. Furthermore, the influence of the anastomosis angle on wall shear stress (WSS) was evaluated. Grafting of left or right subclavian arteries has similar hemodynamic performance in comparison to standard cannula positions. Angularity change of the graft anastomosis from 90° to 30° slightly increases the coronary and cerebral blood flow by 6-9% while significantly reduces the WSS by 35%. Cannulation of the SA is a feasible anastomosis location for the HVAD in the investigated vessel geometry.

Investigation of hemodynamics during cardiopulmonary bypass: A multiscale multiphysics fluid-structure-interaction study.

Neurological complications often occur during cardiopulmonary bypass (CPB). Hypoperfusion of brain tissue due to diminished cerebral autoregulation (CA) and thromboembolism from atherosclerotic plaque reduce the cerebral oxygen supply and increase the risk of perioperative stroke. To improve the outcome of cardiac surgeries, patient-specific computational fluid dynamic (CFD) models can be used to investigate the blood flow during CPB. In this study, we establish a computational model of CPB which includes cerebral autoregulation and movement of aortic walls on the basis of in vivo measurements. First, the Baroreflex mechanism, which plays a leading role in CA, is represented with a 0-D control circuit and coupled to the 3-D domain with differential equations as boundary conditions. Additionally a two-way coupled fluid-structure interaction (FSI) model with CA is set up. The wall shear stress (WSS) distribution is computed for the whole FSI domain and a comparison to rigid wall CFD is made. Constant flow and pulsatile flow CPB is considered. Rigid wall CFD delivers higher wall shear stress values than FSI simulations, especially during pulsatile perfusion. The flow rates through the supraaortic vessels are almost not affected, if considered as percentages of total cannula output. The developed multiphysic multiscale framework allows deeper insights into the underlying mechanisms during CPB on a patient-specific basis.

Numerical prediction of thrombus risk in an anatomically dilated left ventricle: the effect of inflow cannula designs.

BACKGROUND: Implantation of a rotary blood pump (RBP) can cause non-physiological flow fields in the left ventricle (LV) which may trigger thrombosis. Different inflow cannula geometry can affect LV flow fields. The aim of this study was to determine the effect of inflow cannula geometry on intraventricular flow under full LV support in a patient specific model. METHODS: Computed tomography angiography imaging of the LV was performed on a RBP candidate to develop a patient-specific model. Five inflow cannulae were evaluated, which were modelled on those used clinically or under development. The inflow cannulae are described as a crown like tip, thin walled tubular tip, large filleted tip, trumpet like tip and an inferiorly flared cannula. Placement of the inflow cannula was at the LV apex with the central axis intersecting the centre of the mitral valve. Full support was simulated by prescribing 5 l/min across the mitral valve. Thrombus risk was evaluated by identifying regions of stagnation. Rate of LV washout was assessed using a volume of fluid model. Relative haemolysis index and blood residence time was calculated using an Eulerian approach. RESULTS: The inferiorly flared inflow cannula had the lowest thrombus risk due to low stagnation volumes. All cannulae had similar rates of LV washout and blood residence time. The crown like tip and thin walled tubular tip resulted in relatively higher blood damage indices within the LV. CONCLUSION: Changes in intraventricular flow due to variances in cannula geometry resulted in different stagnation volumes. Cannula geometry does not appreciably affect LV washout rates and blood residence time. The patient specific, full support computational fluid dynamic model provided a repeatable platform to investigate the effects of inflow cannula geometry on intraventricular flow.

A numerical framework to investigate hemodynamics during endovascular mechanical recanalization in acute stroke.

Ischemic stroke, caused by embolism of cerebral vessels, inflicts high morbidity and mortality. Endovascular aspiration of the blood clot is an interventional technique for the recanalization of the occluded arteries. However, the hemodynamics in the Circle of Willis (CoW) are not completely understood, which results in medical misjudgment and complications during surgeries. In this study we establish a multiscale description of cerebral hemodynamics during aspiration thrombectomy. First, the CoW is modeled as a 1D pipe network on the basis of computed tomography angiography (CTA) scans. Afterwards, a vascular occlusion is placed in the middle cerebral artery and the relevant section of the CoW is transferred to a 3D computational fluid dynamic (CFD) domain. A suction catheter in different positions is included in the CFD simulations. The boundary conditions of the 3D domain are taken from the 1D domain to ensure system coupling. A Eulerian-Eulerian multiphase simulation describes the process of thrombus aspiration. The physiological blood flow in the 1D and 3D domains is validated with literature data. Further on, it is proved that domain reduction and pressure coupling at the boundaries are an appropriate method to reduce computational costs. Future work will apply the developed framework to various clinical questions.

In vitro flow investigations in the aortic arch during cardiopulmonary bypass with stereo-PIV.

The cardiopulmonary bypass is related to complications like stroke or hypoxia. The cannula jet is suspected to be one reason for these complications, due to the sandblast effect on the vessel wall. Several in silico and in vitro studies investigated the underlying mechanisms, but the applied experimental flow measurement techniques were not able to address the highly three-dimensional flow character with a satisfying resolution. In this work in vitro flow measurements in a cannulated and a non-cannulated aortic silicone model are presented. Stereo particle image velocimetry measurements in multiple planes were carried out. By assembling the data of the different measurement planes, quasi 3D velocity fields with a resolution of~1.5×1.5×2.5 mm(3) were obtained. The resulting velocity fields have been compared regarding magnitude, streamlines and vorticity. The presented method shows to be a suitable in vitro technique to measure and address the three-dimensional aortic CPB cannula flow with a high temporal and spatial resolution.

Development of a numerical pump testing framework.

It has been shown that left ventricular assist devices (LVADs) increase the survival rate in end-stage heart failure patients. However, there is an ongoing demand for an increased quality of life, fewer adverse events, and more physiological devices. These challenges necessitate new approaches during the design process. In this study, computational fluid dynamics (CFD), lumped parameter (LP) modeling, mock circulatory loops (MCLs), and particle image velocimetry (PIV) are combined to develop a numerical Pump Testing Framework (nPTF) capable of analyzing local flow patterns and the systemic response of LVADs. The nPTF was created by connecting a CFD model of the aortic arch, including an LVAD outflow graft to an LP model of the circulatory system. Based on the same geometry, a three-dimensional silicone model was crafted using rapid prototyping and connected to an MCL. PIV studies of this setup were performed to validate the local flow fields (PIV) and the systemic response (MCL) of the nPTF. After validation, different outflow graft positions were compared using the nPTF. Both the numerical and the experimental setup were able to generate physiological responses by adjusting resistances and systemic compliance, with mean aortic pressures of 72.2-132.6 mm Hg for rotational speeds of 2200-3050 rpm. During LVAD support, an average flow to the distal branches (cerebral and subclavian) of 24% was found in the experiments and the nPTF. The flow fields from PIV and CFD were in good agreement. Numerical and experimental tools were combined to develop and validate the nPTF, which can be used to analyze local flow fields and the systemic response of LVADs during the design process. This allows analysis of physiological control parameters at early development stages and may, therefore, help to improve patient outcomes.

Analysis of emboli and blood flow in the ophthalmic artery to understand retinal artery occlusion.

Retinal artery occlusion (RAO) is a common ocular vascular occlusive disorder that may lead to partial or complete retinal ischemia with sudden visual deterioration and visual field defects. Although RAO has been investigated since 1859, the main mechanism is still not fully understood. While hypoperfusion of the ophthalmic artery (OA) due to severe stenosis of the internal carotid artery might lead to RAO, emboli are assumed to be the main reason. Intra-arterial thrombolysis is not a sufficient treatment for RAO, and current research is mainly focused on risk factors. In this study, a computational fluid dynamic model is presented to analyse flow conditions and clot behaviour at the junction of the internal carotid artery and OA based on a realistic geometry from a RAO patient. Clot diameters varied between 5 and 200 µm, and the probability of clots reaching the OA or being washed into the brain was analysed. Results show sufficient blood flow and perfusion pressure at the end of OA. The probability that clots from the main blood flow will to be washed into the brain is 7.32 ± 1.08%. A wall shear stress hotspot is observed at the curvature proximal to the internal carotid artery/OA junction. Clots released from this hotspot have a higher probability of causing RAO. The occurrence of such patient-specific pathophysiologies will have to be considered in the future.

Design modifications and computational fluid dynamic analysis of an outflow cannula for cardiopulmonary bypass.

Cardiopulmonary bypass is a well-established technique during open heart surgeries. However, neurological complications due to insufficient cerebral oxygen supply occur and the severe consequences must not be neglected. Recent computational fluid dynamics (CFD) studies showed that during axillary cannulation the cerebral perfusion is strongly affected by the distance between the cannula tip and the vertebral artery branch. In this study we use two modifications of the cannula design to analyze the flow characteristics by means of CFD and experimental validation with particle image velocimetry (PIV). One approach applies a spin to the blood stream with a helical surface in the cannula cross section. Another approach uses radial bores in a constricted cannula tip to split the outflow jet. The additional helicity improves the perfusion of the cerebral vessels and suppresses the blood suction in the right vertebral artery observed with a standard cannula. The cannula with a helix throughout the entire length changes the blood flow from 2124 to 32 mL/min in comparison with an unmodified design and has the lowest prediction of blood damage. Separating the blood stream does not deliver satisfying results. The PIV measurements validate the simulations and correspond with the velocity distribution as well as vortex locations.

Numerical washout study of a pulsatile total artificial heart.

PURPOSE: For blood pumps with long term indication, blood stagnation can result in excessive thromboembolic risks for patients. This study numerically investigates the washout performance of the left pump chamber of a pulsatile total artificial heart (TAH) as well as the sensitivity of the rotational orientation of the inlet bileaflet mechanical heart valve (MHV) on blood stagnation. METHODS: To quantitatively evaluate the washout efficiency, a fluid-structure interaction (FSI) simulation of the artificial heart pumping process was combined with a blood washout model. Four geometries with different orientations (0°, 45°, 90° and 135°) of the inlet valve were compared with respect to washout performance. RESULTS: The calculated flow field showed a high level of agreement with particle image velocimetry (PIV) measurements. Almost complete washout was achievable after three ejection phases. Remains of old blood in relation to the chamber volume was below 0.6% for all configurations and were mainly detected opposite to the inlet and outlet port at the square edge where the membrane and the pump chamber are connected. Only a small variation in the washout efficiency and the general flow field was observed. An orientation of 0° showed minor advantages with respect to blood stagnation and recirculation. CONCLUSIONS: Bileaflet MHVs were demonstrated to be only slightly sensitive to rotation regarding the washout performance of the TAH. The proposed numerical washout model proved to be an adequate tool to quantitatively compare different configurations and designs of the artificial organ regarding the potential for blood stagnation where experimental measurements are limited.

A multiscale 0-D/3-D approach to patient-specific adaptation of a cerebral autoregulation model for computational fluid dynamics studies of cardiopulmonary bypass.

Neurological complication often occurs during cardiopulmonary bypass (CPB). One of the main causes is hypoperfusion of the cerebral tissue affected by the position of the cannula tip and diminished cerebral autoregulation (CA). Recently, a lumped parameter approach could describe the baroreflex, one of the main mechanisms of cerebral autoregulation, in a computational fluid dynamics (CFD) study of CPB. However, the cerebral blood flow (CBF) was overestimated and the physiological meaning of the variables and their impact on the model was unknown. In this study, we use a 0-D control circuit representation of the Baroreflex mechanism, to assess the parameters with respect to their physiological meaning and their influence on CBF. Afterwards the parameters are transferred to 3D-CFD and the static and dynamic behavior of cerebral autoregulation is investigated. The parameters of the baroreflex mechanism can reproduce normotensive, hypertensive and impaired autoregulation behavior. Further on, the proposed model can mimic the effects of anesthetic agents and other factors controlling dynamic CA. The CFD simulations deliver similar results of static and dynamic CBF as the 0-D control circuit. This study shows the feasibility of a multiscale 0-D/3-D approach to include patient-specific cerebral autoregulation into CFD studies.

Development of a hemodynamically optimized outflow cannula for cardiopulmonary bypass.

The jet of the outflow cannula is a potential risk for patients undergoing cardiopulmonary bypass (CPB), because increased jet velocities lead to altered flow conditions and might furthermore mobilize atherosclerotic plaques from calcified aortas. The cannula jet is therefore among the main reasons for cerebral hypoxia and stroke in CPB patients. In the past, we developed a validated computational fluid dynamics (CFD) model to analyze flow conditions during CPB as dependent on cannulation and support modalities. This model is now applied to develop a novel CPB outflow cannula to reduce the jet effect and increase cerebral blood flow. The Multi-Module Cannula (MMC) is based on a generic elbow cannula that was iteratively improved. It features an inner wall to smoothly guide the blood as well as an elliptically shaped outlet diffuser. During standard CPB conditions of 5 L/min, the pressure drop over the MMC is 61 mm Hg, compared with 68 mm Hg with a standard cannula. The maximum velocities are decreased from 3.7 m/s to 3.3 m/s. In the cannula jet of the MMC, the velocities are reduced further, down to 1.6 m/s. The cerebral blood flow is typically reduced during CPB. Using the MMC, however, it reaches almost physiological values at 715 mL/min. These results suggest that the MMC outperforms standard CPB cannulas. Further design improvements and improved insertion techniques are under consideration.

Implementation of intrinsic lumped parameter modeling into computational fluid dynamics studies of cardiopulmonary bypass.

Stroke and cerebral hypoxia are among the main complications during cardiopulmonary bypass (CPB). The two main reasons for these complications are the cannula jet, due to altered flow conditions and the sandblast effect, and impaired cerebral autoregulation which often occurs in the elderly. The effect of autoregulation has so far mainly been modeled using lumped parameter modeling, while Computational Fluid Dynamics (CFD) has been applied to analyze flow conditions during CPB. In this study, we combine both modeling techniques to analyze the effect of lumped parameter modeling on blood flow during CPB. Additionally, cerebral autoregulation is implemented using the Baroreflex, which adapts the cerebrovascular resistance and compliance based on the cerebral perfusion pressure. The results show that while a combination of CFD and lumped parameter modeling without autoregulation delivers feasible results for physiological flow conditions, it overestimates the loss of cerebral blood flow during CPB. This is counteracted by the Baroreflex, which restores the cerebral blood flow to native levels. However, the cerebral blood flow during CPB is typically reduced by 10-20% in the clinic. This indicates that either the Baroreflex is not fully functional during CPB, or that the target value for the Baroreflex is not a full native cerebral blood flow, but the plateau phase of cerebral autoregulation, which starts at approximately 80% of native flow.

Implementation of cerebral autoregulation into computational fluid dynamics studies of cardiopulmonary bypass.

Peri- or postoperative neurological complications are among the main risks for patients undergoing extracorporeal circulatory support (ECC). Two of the main reasons are an increased risk for strokes and altered flow conditions leading to cerebral hypoperfusion. This is strongly affected by cerebral autoregulation, which is the body's intrinsic ability to provide sufficient cerebral blood flow (CBF) despite changes in cerebral perfusion pressure (CPP). This complex mechanism has been mainly neglected in numerical studies, which have often been applied for analysis of ECC. In this study, a mathematical model is presented to implement cerebral autoregulation into computational fluid dynamics (CFD) studies. CFD simulations of cardiopulmonary bypass (CPB) were performed in a 3D model of the cardiovascular system, with flow variations between 4.5-6 L/min. Cerebral outlets were modeled using an equation to calculate CBF based on CPP. Assuming full regulation, CBF was kept constant for CPP between 80 and 120 mm Hg. A deviation in CBF of 20% occurred for CPP between 55-80 mm Hg and 120-145 mm Hg, respectively. The level of regulation was varied to take possible impairment of cerebral autoregulation into account. Furthermore, chronic hypertension was modeled by increasing the baseline CPP. Results indicate that even for full autoregulation, CBF is decreased during CPB. It is even lower for impaired autoregulation and hypertensive patients, demonstrating the strong impact of autoregulation on CBF. It is therefore imperative to include this mechanism into CFD studies. The presented model can help to improve CPB support conditions based on patient-specific autoregulation parameters.

Mimicking of cerebral autoregulation by flow-dependent cerebrovascular resistance: a feasibility study.

Understanding circulatory autoregulation is essential for improving physiological control of rotary blood pumps and support conditions during cardiopulmonary bypass (CPB). Cerebral autoregulation (CAR), arguably the most critical, is the body's intrinsic ability to maintain sufficient cerebral blood flow (CBF) despite changes in aortic perfusion pressure. It is therefore imperative to include this mechanism into computational fluid dynamics (CFD), particle image velocimetry (PIV), or mock circulation loop (MCL) studies. Without such inclusions, potential losses of CBF are overestimated. In this study, a mathematical model to mimic CAR is implemented in a MCL- and PIV-validated CFD model. A three-dimensional model of the human vascular system was created from magnetic resonance imaging records. Numerical flow simulations were performed for physiological conditions and CPB. The inlet flow was varied between 4.5 and 6 L/min. Arterial outlets were modeled using vessel-specific, flow-dependent cerebrovascular resistances (CVRs), resulting in a variation of the pressure drop between 0 and 80mmHg. CBF is highly dependent on the level of CAR during CPB. By varying the CVR parameters up to the beginning of plateau phase, it can be regulated between 0 and 80% of physiological CBF. So while implementing autoregulation, CBF remains unchanged during a simulated native cardiac output of 5L/min or CPB support of 6L/min. Neglecting CAR, constant backflow from the brain occurs for some cannula positions. Using flow-dependent CVR, CBF returns to its baseline at a rate of recovery of 0.25s. Results demonstrate that modeling of CAR by flow-dependent CVR delivers feasible results. The presented method can be used to optimize physiological control of assist devices dependent upon different levels of CAR representing different patients.