M. Ertan Taskin

The use of computational fluid dynamics in the development of ventricular assist devices

Progress in the field of prosthetic cardiovascular devices has significantly contributed to the rapid advancements in cardiac therapy during the last four decades. The concept of mechanical circulatory assistance was established with the first successful clinical use of heart-lung machines for cardiopulmonary bypass. Since then a variety of devices have been developed to replace or assist diseased components of the cardiovascular system. Ventricular assist devices (VADs) are basically mechanical pumps designed to augment or replace the function of one or more chambers of the failing heart. Computational Fluid Dynamics (CFD) is an attractive tool in the development process of VADs, allowing numerous different designs to be characterized for their functional performance virtually, for a wide range of operating conditions, without the physical device being fabricated. However, VADs operate in a flow regime which is traditionally difficult to simulate; the transitional region at the boundary of laminar and turbulent flow. Hence different methods have been used and the best approach is debatable. In addition to these fundamental fluid dynamic issues, blood consists of biological cells. Device-induced biological complications are a serious consequence of VAD use. The complications include blood damage (haemolysis, blood cell activation), thrombosis and emboli. Patients are required to take anticoagulation medication constantly which may cause bleeding. Despite many efforts blood damage models have still not been implemented satisfactorily into numerical analysis of VADs, which severely undermines the full potential of CFD. This paper reviews the current state of the art CFD for analysis of blood pumps, including a practical critical review of the studies to date, which should help device designers choose the most appropriate methods; a summary of blood damage models and the difficulties in implementing them into CFD; and current gaps in knowledge and areas for future work.

A Quantitative Comparison of Mechanical Blood Damage Parameters in Rotary Ventricular Assist Devices: Shear Stress, Exposure Time and Hemolysis Index

Ventricular assist devices (VADs) have already helped many patients with heart failure but have the potential to assist more patients if current problems with blood damage (hemolysis, platelet activation, thrombosis and emboli, and destruction of the von Willebrand factor (vWf)) can be eliminated. A step towards this goal is better understanding of the relationships between shear stress, exposure time, and blood damage and, from there, the development of numerical models for the different types of blood damage to enable the design of improved VADs. In this study, computational fluid dynamics (CFD) was used to calculate the hemodynamics in three clinical VADs and two investigational VADs and the shear stress, residence time, and hemolysis were investigated. A new scalar transport model for hemolysis was developed. The results were compared with in vitro measurements of the pressure head in each VAD and the hemolysis index in two VADs. A comparative analysis of the blood damage related fluid dynamic parameters and hemolysis index was performed among the VADs. Compared to the centrifugal VADs, the axial VADs had: higher mean scalar shear stress (sss); a wider range of sss, with larger maxima and larger percentage volumes at both low and high sss; and longer residence times at very high sss. The hemolysis predictions were in agreement with the experiments and showed that the axial VADs had a higher hemolysis index. The increased hemolysis in axial VADs compared to centrifugal VADs is a direct result of their higher shear stresses and longer residence times. Since platelet activation and destruction of the vWf also require high shear stresses, the flow conditions inside axial VADs are likely to result in more of these types of blood damage compared with centrifugal VADs.

Study of Flow‐Induced Hemolysis Using Novel Couette‐Type Blood‐Shearing Devices

To assist the development and application of blood-contacting medical devices, two novel flow-through Couette-type blood-shearing devices have been developed to study the quantitative relationship between blood damage indexes and flow-dependent parameters. One device is an axial flow-through Couette-type device supported by a pair of pin bearings adapted from the adult Jarvik 2000 blood pump. The other is a centrifugal flow-through Couette-type device supported with magnetic bearings adapted from the CentriMag blood pump. In both devices, a rotor spindle was used to replace the original impeller blades so that a small gap was created between the housing and the rotating spindle surface. Computational fluid dynamics simulations have shown that a uniform, high shear stress region can be generated inside the small gap while the shear stresses elsewhere are relatively low. The possibility of secondary blood damage caused by mechanical seals was eliminated due to the use of a magnetic rotor system. Blood flow through the gap was driven by an externally pressurized reservoir. By adjusting the rotational speed and blood flow rate, shear-induced hemolysis was quantified at a matrix of exposure time (0.039 to 1.48 s) and shear stress (50 to 320 Pa). All of the experiments were conducted at room temperature using heparinized ovine blood with a hematocrit value of 30%. The measured hemolysis levels were much lower than those published in the literature, and the overestimation of those earlier studies may be attributable to device-related secondary blood-damaging effects. A new set of coefficients for the power law model was derived from the regression of the experimental data.

Evaluation of Eulerian and Lagrangian models for hemolysis estimation

Hemolysis caused by flow-induced mechanical damage to red blood cells is still a problem in medical devices such as ventricular assist devices (VADs), artificial lungs, and mechanical heart valves. A number of different models have been proposed by different research groups for calculating the hemolysis, and of these, the power law–based models (HI(%)=Ct&agr;&tgr;&bgr;) have proved the most popular because of their ease of use and applicability to a wide range of devices. However, within this power law category of models there are a number of different implementations. The aim of this work was to evaluate different power law–based models by calculating hemolysis in a specifically designed shearing device and a clinical VAD, and comparing the estimated results with experimental measurements of the hemolysis in these two devices. Both the Eulerian scalar transport and all the Lagrangian models had fairly large percentage of errors compared with the experiments (minimum Eulerian 91% and minimum Lagrangian 57%) showing they could not accurately predict the magnitude of the hemolysis. However, the Eulerian approach had large correlation coefficients (>0.99) showing that this method can predict relative hemolysis, which would be useful in comparative analysis, for example, for ranking different devices or for design optimization studies.

Computational characterization of flow and hemolytic performance of the UltraMag blood pump for circulatory support

The Levitronix UltraMag blood pump is a next generation, magnetically suspended centrifugal pump and is designed to provide circulatory support for pediatric and adult patients. The aim of this study is to investigate the hemodynamic and hemolytic characteristics of this pump using the computational fluid dynamics (CFD) approach. The computational domain for CFD analysis was constructed from the three-dimensional geometry (3D) of the UltraMag blood pump and meshed into 3D tetrahedral/hybrid elements. The governing equations of fluid flow were computationally solved to obtain a blood flow through the blood pump. Further, hemolytic blood damage was calculated by solving a scalar transport equation where the scalar variable and the source term were obtained utilizing an empirical power-law correlation between the fluid dynamic variables and hemolysis. To obtain mesh independent flow solution, a comparative examination of vector fields, hydrodynamic performance, and hemolysis predictions were carried out. Different sizes of tetrahedral and tetrahedral/hexahedral mixed hybrid models were considered. The mesh independent solutions were obtained by a hybrid model. Laminar and SST κ-ω turbulence flow models were used for different operating conditions. In order to pinpoint the most significant hemolytic region, the flow field analysis was coupled to the hemolysis predictions. In summary, computational characterization of the device was satisfactorily carried out within the targeted operating conditions of the device, and it was observed that the UltraMag blood pump can be safely operated for its intended use to create a circulatory support for both pediatric and adult-sized patients.

Computational Fluid Dynamics Analysis of Thrombosis Potential in Left Ventricular Assist Device Drainage Cannulae

Cannulation is necessary when blood is removed from the body, for example in hemodialysis, cardiopulmonary bypass, blood oxygenators, and ventricular assist devices. Artificial blood contacting surfaces are prone to thrombosis, especially in the presence of stagnant or recirculating flow. In this work, computational fluid dynamics was used to investigate the blood flow fields in three clinically available cannulae (Medtronic DLP 12, 16, and 24 F), used as drainage for pediatric circulatory support and to calculate parameters that may be indicative of thrombosis potential. The results show that using the 24 F cannula below flow rates of about 0.75 L/min produces hemodynamic conditions, which may increase the risk of blood clotting within the cannula. No reasons are indicated for not using the 12 or 16 F cannulae with flow rates between 0.25 and 3.0 L/min.

Computational Design and In Vitro Characterization of an Integrated Maglev Pump-Oxygenator

For the need for respiratory support for patients with acute or chronic lung diseases to be addressed, a novel integrated maglev pump-oxygenator (IMPO) is being developed as a respiratory assist device. IMPO was conceptualized to combine a magnetically levitated pump/rotor with uniquely configured hollow fiber membranes to create an assembly-free, ultracompact system. IMPO is a self-contained blood pump and oxygenator assembly to enable rapid deployment for patients requiring respiratory support or circulatory support. In this study, computational fluid dynamics (CFD) and computer-aided design were conducted to design and optimize the hemodynamics, gas transfer, and hemocompatibility performances of this novel device. In parallel, in vitro experiments including hydrodynamic, gas transfer, and hemolysis measurements were conducted to evaluate the performance of IMPO. Computational results from CFD analysis were compared with experimental data collected from in vitro evaluation of the IMPO. The CFD simulation demonstrated a well-behaved and streamlined flow field in the main components of this device. The results of hydrodynamic performance, oxygen transfer, and hemolysis predicted by computational simulation, along with the in vitro experimental data, indicate that this pump-lung device can provide the total respiratory need of an adult with lung failure, with a low hemolysis rate at the targeted operating condition. These detailed CFD designs and analyses can provide valuable guidance for further optimization of this IMPO for long-term use.

Computational model‐based design of a wearable artificial pump‐lung for cardiopulmonary/respiratory support

Mechanical ventilation and extracorporeal membrane oxygenation are the only immediate options available for patients with respiratory failure. However, these options present significant shortcomings. To address this unmet need for respiratory support, innovative respiratory assist devices are being developed. In this study, we present the computational model-based design, and analysis of functional characteristics and hemocompatibility performance, of an innovative wearable artificial pump-lung (APL) for ambulatory respiratory support. Computer-aided design and computational fluid dynamics (CFD)-based modeling were utilized to generate the geometrical model and to acquire the fluid flow field, gas transfer, and blood damage potential. With the knowledge of flow field, gas transfer, and blood damage potential through the whole device, design parameters were adjusted to achieve the desired specifications based on the concept of virtual prototyping using the computational modeling in conjunction with consideration of the constraints on fabrication processes and materials. Based on the results of the CFD design and analysis, the physical model of the wearable APL was fabricated. Computationally predicted hydrodynamic pumping function, gas transfer, and blood damage potential were compared with experimental data from in vitro evaluation of the wearable APL. The hydrodynamic performance, oxygen transfer, and blood damage potential predicted with computational modeling, along with the in vitro experimental data, indicated that this APL meets the design specifications for respiratory support with excellent biocompatibility at the targeted operating condition.

Computational Study of the Blood Flow in Three Types of 3D Hollow Fiber Membrane Bundles

The goal of this study is to develop a computational fluid dynamics (CFD) modeling approach to better estimate the blood flow dynamics in the bundles of the hollow fiber membrane based medical devices (i.e., blood oxygenators, artificial lungs, and hemodialyzers). Three representative types of arrays, square, diagonal, and random with the porosity value of 0.55, were studied. In addition, a 3D array with the same porosity was studied. The flow fields between the individual fibers in these arrays at selected Reynolds numbers (Re) were simulated with CFD modeling. Hemolysis is not significant in the fiber bundles but the platelet activation may be essential. For each type of array, the average wall shear stress is linearly proportional to the Re. For the same Re but different arrays, the average wall shear stress also exhibits a linear dependency on the pressure difference across arrays, while Darcys law prescribes a power-law relationship, therefore, underestimating the shear stress level. For the same Re, the average wall shear stress of the diagonal array is approximately 3.1, 1.8, and 2.0 times larger than that of the square, random, and 3D arrays, respectively. A coefficient C is suggested to correlate the CFD predicted data with the analytical solution, and C is 1.16, 1.51, and 2.05 for the square, random, and diagonal arrays in this paper, respectively. It is worth noting that C is strongly dependent on the array geometrical properties, whereas it is weakly dependent on the flow field. Additionally, the 3D fiber bundle simulation results show that the three-dimensional effect is not negligible. Specifically, velocity and shear stress distribution can vary significantly along the fiber axial direction.

THE EFFECT OF IMPELLER POSITION ON CFD CALCULATIONS OF BLOOD FLOW IN MAGNETICALLY LEVITATED CENTRIFUGAL BLOOD PUMPS

Cardiovascular disease is the leading cause of mortality globally. Among various forms of cardiovascular disease, heart failure (HF) affects 5.7 million patients in the United States with about 670,000 new patients diagnosed for the first time annually (1). The fatality rate for HF is high, with one in five people dying within 1 year (1). The number of deaths has increased (1) despite advances in surgical treatment and new pharmaceutical therapies. Many therapies are available to treat patients with HF, including lifestyle changes, medications, transcatheter interventions and surgery. However, despite optimal medical and surgical therapies, some patients still do not improve and the available therapies fail to control their symptoms; for them, cardiac transplantation may be the only treatment option. However, only approximately 2300 donor hearts become available each year resulting in around 2200 transplants (1), or only about 6% of the estimated 35,000 US patients who would benefit from a heart actually receiving a transplant. To address the need to support the circulation in patients with end-stage HF a wide variety of mechanical circulatory support devices (MCSDs) have been developed over the past four decades. These MCSDs have been developed as a bridge to transplant, a bridge to recovery, and as an end stage treatment. They can be implanted as a ventricular assist device (VAD) to support the left ventricle (LVAD) or the right ventricle (RVAD) or two devices are used to support both left and right ventricles (Bi-VAD).Copyright