Simon J. Sonntag

Fluid-Structure Interaction Model of a Percutaneous Aortic Valve: Comparison with an In Vitro Test and Feasibility Study in a Patient-Specific Case.

Transcatheter aortic valve replacement (TAVR) represents an established recent technology in a high risk patient base. To better understand TAVR performance, a fluid–structure interaction (FSI) model of a self-expandable transcatheter aortic valve was proposed. After an in vitro durability experiment was done to test the valve, the FSI model was built to reproduce the experimental test. Lastly, the FSI model was used to simulate the virtual implant and performance in a patient-specific case. Results showed that the leaflet opening area during the cycle was similar to that of the in vitro test and the difference of the maximum leaflet opening between the two methodologies was of 0.42%. Furthermore, the FSI simulation quantified the pressure and velocity fields. The computed strain amplitudes in the stent frame showed that this distribution in the patient-specific case is highly affected by the aortic root anatomy, suggesting that the in vitro tests that follow standards might not be representative of the real behavior of the percutaneous valve. The patient-specific case also compared in vivo literature data on fast opening and closing characteristics of the aortic valve during systolic ejection. FSI simulations represent useful tools in determining design errors or optimization potentials before the fabrication of aortic valve prototypes and the performance of tests.

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.

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.

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.

Effects of Interaction Between Ventricular Assist Device Assistance and Autoregulated Mock Circulation Including Frank-Starling Mechanism and Baroreflex.

A mock heart circulation loop (MHCL) is a hydraulic model simulating the human circulatory system. It allows in vitro investigations of the interaction between cardiac assist devices and the human circulatory system. In this study, a preload sensitive MHCL, the MHCLAUTO , was developed to investigate the interaction between the left ventricle and left ventricular assist devices (LVADs). The Frank-Starling mechanism was modeled by regulating the stroke volume (SV) based on the measured mean diastolic left atrial pressure (MLAPdiast ). The baroreflex autoregulation mechanism was implemented to maintain a constant mean aortic pressure (MAP) by varying ventricular contractility (Emax ), heart rate (HR), afterload/systemic vascular resistance (SVR) and unstressed venous volume (UVV). The DP3 blood pump (Medos Medizintechnik GmbH) was used to simulate the LVAD. Characteristic parameters were measured in pathological conditions both with and without LVAD to assess the hemodynamic effect of LVAD on the MHCLAUTO . The results obtained from the MHCLAUTO show a high correlation to literature data. The study demonstrates the possibility of using the MHCLAUTO as a research tool to better understand the physiological interactions between cardiac implants and human circulation.

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.

In Vitro Evaluation of Intra-Aneurysmal, Flow-Diverter-Induced Thrombus Formation: A Feasibility Study

BACKGROUND AND PURPOSE: Intracranial aneurysm treatment by flow diverters aims at triggering intra-aneurysmal thrombosis. By combining in vitro blood experiments with particle imaging velocimetry measurements, we investigated the time-resolved thrombus formation triggered by flow diverters. MATERIALS AND METHODS: Two test setups were built, 1 for particle imaging velocimetry and 1 for blood experiments, both generating the same pulsatile flow and including a silicone aneurysm model. Tests without flow diverters and with 2 different flow-diverter sizes (diameter: 4.5 and 4.0 mm) were performed. In the blood experiments, the intra-aneurysmal flow was monitored by using Doppler sonography. The experiments were stopped at 3 different changes of the spatial extent of the signal. RESULTS: No thrombus was detected in the aneurysm model without the flow diverter. Otherwise, thrombi were observed in all aneurysm models with flow diverters. The thrombi grew from the proximal side of the aneurysm neck with fibrin threads connected to the flow diverter and extending across the aneurysm. The thrombus resulting from the 4.0-mm flow diverter grew along the aneurysm wall as a solid and organized thrombus, which correlates with the slower velocities near the wall detected by particle imaging velocimetry. The thrombus that evolved by using the 4.5-mm flow diverter showed no identifiable growing direction. The entire thrombus presumably resulted from stagnation of blood and correlates with the central vortex detected by particle imaging velocimetry. CONCLUSIONS: We showed the feasibility of in vitro investigation of time-resolved thrombus formation in the presence of flow diverters.

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.

Fluid Dynamics in Rotary Piston Blood Pumps

Mechanical circulatory support can maintain a sufficient blood circulation if the native heart is failing. The first implantable devices were displacement pumps with membranes. They were able to provide a sufficient blood flow, yet, were limited because of size and low durability. Rotary pumps have resolved these technical drawbacks, enabled a growing number of mechanical circulatory support therapy and a safer application. However, clinical complications like gastrointestinal bleeding, aortic insufficiency, thromboembolic complications, and impaired renal function are observed with their application. This is traced back to their working principle with attenuated or non-pulsatile flow and high shear stress. Rotary piston pumps potentially merge the benefits of available pump types and seem to avoid their complications. However, a profound assessment and their development requires the knowledge of the flow characteristics. This study aimed at their investigation. A functional model was manufactured and investigated with particle image velocimetry. Furthermore, a fluid–structure interaction computational simulation was established to extend the laboratory capabilities. The numerical results precisely converged with the laboratory measurements. Thus, the in silico model enabled the investigation of relevant areas like gap flows that were hardly feasible with laboratory means. Moreover, an economic method for the investigation of design variations was established.

Combined Computational and Experimental Approach to Improve the Assessment of Mitral Regurgitation by Echocardiography

Mitral regurgitation (MR) is one of the most frequent valvular heart diseases. To assess MR severity, color Doppler imaging (CDI) is the clinical standard. However, inadequate reliability, poor reproducibility and heavy user-dependence are known limitations. A novel approach combining computational and experimental methods is currently under development aiming to improve the quantification. A flow chamber for a circulatory flow loop was developed. Three different orifices were used to mimic variations of MR. The flow field was recorded simultaneously by a 2D Doppler ultrasound transducer and Particle Image Velocimetry (PIV). Computational Fluid Dynamics (CFD) simulations were conducted using the same geometry and boundary conditions. The resulting computed velocity field was used to simulate synthetic Doppler signals. Comparison between PIV and CFD shows a high level of agreement. The simulated CDI exhibits the same characteristics as the recorded color Doppler images. The feasibility of the proposed combination of experimental and computational methods for the investigation of MR is shown and the numerical methods are successfully validated against the experiments. Furthermore, it is discussed how the approach can be used in the long run as a platform to improve the assessment of MR quantification.