Personalized Numerical Cardiovascular Model with Weight Growth for Evaluating Pediatric Left Ventricular Assist Devices: Derivation from an Experimental Mock Circulatory Loop.

Pediatric patients with heart failure have limited treatment options because of a shortage of donor hearts and compatible left ventricular assist devices (LVADs). To address this issue, our group is developing an implantable pediatric LVAD for patients weighing 5-20 kg, capable of accommodating different physiological hemodynamic conditions as patients grow. To evaluate LVAD prototypes across a wide range of conditions, we developed a numerical cardiovascular model, using data from a mock circulatory loop (MCL) and patient-specific elastance functions. The numerical MCL was validated against experimental MCL results, showing good agreement, with differences ranging from 0 to 11%. The numerical model was also tested under left heart failure conditions and showed a worst-case difference of 16%. In an MCL study with a pediatric LVAD, a pediatric dataset was obtained from the experimental MCL and used to tune the numerical MCL. Then, the numerical model simulated LVAD flow by using an HQ curve obtained from the LVAD's impeller. When the numerical MCL was validated against the experimental MCL, hemodynamic differences ranged between 0 and 9%. These findings suggest that the numerical model can replicate various physiological conditions and impeller designs, indicating its potential as a tool for developing and optimizing pediatric LVADs.

Blood-contacting magnetic levitation bearing design using computational fluid dynamics for haemocompatibility.

The Hemocompatibility Assessment Platform (HAP) is a testing rig that will allow for the evaluation of blood trauma caused by individual components of rotary blood pumps including the NeoVAD - a proposed paediatric Left Ventricular Assist Device (LVAD). It is important that the HAP itself is only minimally haemolytic such that the plasma free haemoglobin measured can be assumed to come from the test component. In this study, Computational Fluid Dynamics simulations have been carried out to inform the design of a magnetically levitated motor bearing gap. Simulations show that issues with the original design, namely stagnation regions and large recirculation zones can be mitigated with the introduction of a pipe that introduces blood-flow to the centre of the bearing and disrupts the secondary flow patterns that cause these issues.Clinical relevance- The consequent reduction in shear exposure time will reduce heamolsyis from the HAP. The redesign of the bearing will result in reduced baseline blood trauma from the HAP, thus allowing quantification of test component haemolysis and will therefore aid the design of future paediatric LVADs.

Dissipated energy and efficiency as objective functions for the design of the NeoVAD rotary blood pump.

Minimising haemolytic blood damage is an important objective when designing rotary blood pumps, however, calculating haemolysis can be computationally expensive and inaccurate. Efficiency and dissipated energy are much more easily calculable hydraulic parameters in the design and analysis of rotary blood pumps and although there is work to suggest that efficiency is not a good indicator of haemocompatibility, i.e. more efficient pumps do not necessarily cause less damage, there is recent speculation that dissipated energy can act as an easily calculable haemolysis analogue.This study shows that for design purposes, optimising for maximum efficiency and minimum dissipated energy are functionally the same as they are inherently and closely linked. Moreover a demonstration of rotary blood pump design has been completed using the NeoVAD paediatric left ventricular assist device optimising for both objective functions. The resulting designs appear similar in rotor blade shape and are similar in hydraulic performance.Clinical relevance- This reinforces the direct link between efficiency and dissipated energy when analysing rotary blood pumps at a given design operating point. This raises questions either of the claim that efficiency cannot be used as an easily calculable analogue for haemolysis or the validity of dissipated energy to act in this same manner.

Machine learning based on computational fluid dynamics enables geometric design optimisation of the NeoVAD blades.

The NeoVAD is a proposed paediatric axial-flow Left Ventricular Assist Device (LVAD), small enough to be implanted in infants. The design of the impeller and diffuser blades is important for hydrodynamic performance and haemocompatibility of the pump. This study aimed to optimise the blades for pump efficiency using Computational Fluid Dynamics (CFD), machine learning and global optimisation. Meshing of each design typically included 6 million hexahedral elements and a Shear Stress Transport turbulence model was used to close the Reynolds Averaged Navier-Stokes equations. CFD models of 32 base geometries, operating at 8 flow rates between 0.5 and 4 L/min, were created to match experimental studies. These were validated by comparison of the pressure-flow and efficiency-flow curves with those experimentally measured for all base prototype pumps. A surrogate model was required to allow the optimisation routine to conduct an efficient search; a multi-linear regression, Gaussian Process Regression and a Bayesian Regularised Artificial Neural Network predicted the optimisation objective at design points not explicitly simulated. A Genetic Algorithm was used to search for an optimal design. The optimised design offered a 5.51% increase in efficiency at design point (a 20.9% performance increase) as compared to the best performing pump from the 32 base designs. An optimisation method for the blade design of LVADs has been shown to work for a single objective function and future work will consider multi-objective optimisation.

Fluid-structure interaction modelling of a positive-displacement Total Artificial Heart.

For those suffering from end-stage biventricular heart failure, and where a heart transplantation is not a viable option, a Total Artificial Heart (TAH) can be used as a bridge to transplant device. The Realheart TAH is a four-chamber artificial heart that uses a positive-displacement pumping technique mimicking the native heart to produce pulsatile flow governed by a pair of bileaflet mechanical heart valves. The aim of this work was to create a method for simulating haemodynamics in positive-displacement blood pumps, using computational fluid dynamics with fluid-structure interaction to eliminate the need for pre-existing in vitro valve motion data, and then use it to investigate the performance of the Realheart TAH across a range of operating conditions. The device was simulated in Ansys Fluent for five cycles at pumping rates of 60, 80, 100 and 120 bpm and at stroke lengths of 19, 21, 23 and 25 mm. The moving components of the device were discretised using an overset meshing approach, a novel blended weak-strong coupling algorithm was used between fluid and structural solvers, and a custom variable time stepping scheme was used to maximise computational efficiency and accuracy. A two-element Windkessel model approximated a physiological pressure response at the outlet. The transient outflow volume flow rate and pressure results were compared against in vitro experiments using a hybrid cardiovascular simulator and showed good agreement, with maximum root mean square errors of 15% and 5% for the flow rates and pressures respectively. Ventricular washout was simulated and showed an increase as cardiac output increased, with a maximum value of 89% after four cycles at 120 bpm 25 mm. Shear stress distribution over time was also measured, showing that no more than [Formula: see text]% of the total volume exceeded 150 Pa at a cardiac output of 7 L/min. This study showed this model to be both accurate and robust across a wide range of operating points, and will enable fast and effective future studies to be undertaken on current and future generations of the Realheart TAH.

MRI-based strain measurements reflect morphological changes following myocardial infarction: A study on the UK Biobank cohort.

In a porcine experimental model of myocardial infarction, a localised, layer-specific, circumferential left ventricular strain metric has been shown to indicate chronic changes in ventricular function post-infarction more strongly than ejection fraction. This novel strain metric might therefore provide useful prognostic information clinically. In this study, existing clinical volume indices, global strains, and the novel, layer-specific strain were calculated for a large human cohort to assess variations in ventricular function and morphology with age, sex, and health status. Imaging and health data from the UK Biobank were obtained, including healthy volunteers and those with a history of cardiovascular illness. In total, 710 individuals were analysed and stratified by age, sex and health. Significant differences in all strain metrics were found between healthy and unhealthy populations, as well as between males and females. Significant differences in basal circumferential strain and global circumferential strain were found between healthy males and females, with males having smaller absolute values for both (all p ≤ 0.001). There were significant differences in the functional variables left ventricular ejection fraction, end-systolic volume, end-systolic volume index and mid-ventricular circumferential strain between healthy and unhealthy male cohorts aged 65-74 (all p ≤ 0.001). These results suggest that whilst regional circumferential strains may be useful clinically for assessing cardiovascular health, care must be taken to ensure critical values are indexed correctly to age and sex, due to the differences in these values observed here.

Video-based valve motion combined with computational fluid dynamics gives stable and accurate simulations of blood flow in the Realheart total artificial heart.

BACKGROUND: Patients with end-stage, biventricular heart failure, and for whom heart transplantation is not an option, may be given a Total Artificial Heart (TAH). The Realheart® is a novel TAH which pumps blood by mimicking the native heart with translation of an atrioventricular plane. The aim of this work was to create a strategy for using Computational Fluid Dynamics (CFD) to simulate haemodynamics in the Realheart®, including motion of the atrioventricular plane and valves. METHODS: The accuracies of four different computational methods for simulating fluid-structure interaction of the prosthetic valves were assessed by comparison of chamber pressures and flow rates with experimental measurements. The four strategies were: prescribed motion of valves opening and closing at the atrioventricular plane extrema; simulation of fluid-structure interaction of both valves; prescribed motion of the mitral valve with simulation of fluid-structure interaction of the aortic valve; motion of both valves prescribed from video analysis of experiments. RESULTS: The most accurate strategy (error in ventricular pressure of 6%, error in flow rate of 5%) used video-prescribed motion. With the Realheart operating at 80 bpm, the power consumption was 1.03 W, maximum shear stress was 15 Pa, and washout of the ventricle chamber after 4 cycles was 87%. CONCLUSIONS: This study, the first CFD analysis of this novel TAH, demonstrates that good agreement between computational and experimental data can be achieved. This method will therefore enable future optimisation of the geometry and motion of the Realheart®.

Discrete responses of erythrocytes, platelets, and von Willebrand factor to shear.

Despite decades of technological advancements in blood-contacting medical devices, complications related to shear flow-induced blood trauma are still frequently observed in clinic. Blood trauma includes haemolysis, platelet activation, and degradation of High Molecular Weight von Willebrand Factor (HMW vWF) multimers, all of which are dependent on the exposure time and magnitude of shear stress. Specifically, accumulating evidence supports that when blood is exposed to shear stresses above a certain threshold, blood trauma ensues; however, it remains unclear how various constituents of blood are affected by discrete shears experimentally. The aim of this study was to expose blood to discrete shear stresses and evaluate blood trauma indices that reflect red cell, platelet, and vWF structure. Citrated human whole blood (n = 6) was collected and its haematocrit was adjusted to 30 ± 2% by adding either phosphate buffered saline (PBS) or polyvinylpyrrolidone (PVP). Viscosity of whole blood was adjusted to 3.0, 12.5, 22.5 and 37.5 mPa·s to yield stresses of 3, 6, 9, 12, 50, 90 and 150 Pa in a custom-developed shearing system. Blood samples were exposed to shear for 0, 300, 600 and 900 s. Haemolysis was measured using spectrophotometry, platelet activation using flow cytometry, and HMW vWF multimer degradation was quantified with gel electrophoresis and immunoblotting. For tolerance to 300, 600 and 900 s of exposure time, the critical threshold of haemolysis was reached after blood was exposed to 90 Pa for 600 s (P < 0.05), platelet activation and HMW vWF multimer degradation were 50 Pa for 600 s and 12 Pa for 300 s respectively (P < 0.05). Our experimental results provide simultaneous comparison of blood trauma indices and thus also the relation between shear duration and magnitude required to induce damage to red cells, platelets, and vWF. Our results also demonstrate that near-physiological shear stress (<12 Pa) is needed in order to completely avoid any form of blood trauma. Therefore, there is an urgent need to design low shear-flow medical devices in order to avoid blood trauma in this blood-contacting medical device field.

Hemocompatibility Assessment Platform Drive System Design: Trade-off between Motor Performance and Hemolysis.

Left ventricular assist devices (LVADs) have long been used to treat adults with heart failure, but LVAD options for pediatric patients with heart failure are lacking. Despite the urgent need for long-term, implantable pediatric LVADs, design challenges such as hemolysis, pump thrombosis, and bleeding persist. We have developed a Hemocompatibility Assessment Platform (HAP) to identify blood trauma from individual LVAD components. A HAP would aid in refining pump components before in vivo testing, thereby preventing unnecessary animal sacrifice and reducing development time and cost. So that the HAP does not confound hemolysis data, the HAP drive system consists of an enlarged air-gap motor coupled to a magnetic levitation system. Although it is known that an enlarged air gap motor will have diminished performance, while the larger gap in the motor will cause less blood damage, the trade-offs are not fully characterized. Therefore, in this study we evaluated these trade-offs to determine an optimal rotor diameter for the HAP drive motor. The motor performance was characterized with an experimental method by determining the torque constant for the HAP drive motor with varied rotor diameters. The torque threshold was set as 10 mNm to achieve a nominal current of 3.5A. Hemolysis in the HAP drive motor gap was estimated by calculating scalar shear stress generated in the HAP motor gap analytically and numerically. A design criterion of 30 Pa was selected for scalar shear stress to achieve minimal hemolysis and platelet activation in the HAP drive system.Clinical Relevance- We evaluated a Hemocompatibility Assessment Platform for developing LVAD prototypes that can best balance motor performance and hemocompatibility. This design method can assist with optimizing the drive system during the research stage and illustrates how motor geometry can be tuned to reduce blood trauma.

Comparison of the within-reader and inter-vendor agreement of left ventricular circumferential strains and volume indices derived from cardiovascular magnetic resonance imaging.

PURPOSE: Volume indices and left ventricular ejection fraction (LVEF) are routinely used to assess cardiac function. Ventricular strain values may provide additional diagnostic information, but their reproducibility is unclear. This study therefore compares the repeatability and reproducibility of volumes, volume fraction, and regional ventricular strains, derived from cardiovascular magnetic resonance (CMR) imaging, across three software packages and between readers. METHODS: Seven readers analysed 16 short-axis CMR stacks of a porcine heart. Endocardial contours were manually drawn using OsiriX and Simpleware ScanIP and repeated in both softwares. The images were also contoured automatically in Circle CVI42. Endocardial global, apical, mid-ventricular, and basal circumferential strains, as well as end-diastolic and end-systolic volume and LVEF were compared. RESULTS: Bland-Altman analysis found systematic biases in contour length between software packages. Compared to OsiriX, contour lengths were shorter in both ScanIP (-1.9 cm) and CVI42 (-0.6 cm), causing statistically significant differences in end-diastolic and end-systolic volumes, and apical circumferential strain (all p<0.006). No differences were found for mid-ventricular, basal or global strains, or left ventricular ejection fraction (all p<0.007). All CVI42 results lay within the ranges of the OsiriX results. Intra-software differences were found to be lower than inter-software differences. CONCLUSION: OsiriX and CVI42 gave consistent results for all strain and volume metrics, with no statistical differences found between OsiriX and ScanIP for mid-ventricular, global or basal strains, or left ventricular ejection fraction. However, volumes were influenced by the choice of contouring software, suggesting care should be taken when comparing volumes across different software.

Normal fluid stresses are prevalent in rotary ventricular assist devices: A computational fluid dynamics analysis.

Despite the evolution of ventricular assist devices, ventricular assist device patients still suffer from complications due to the damage to blood by fluid dynamic stress. Since rotary ventricular assist devices are assumed to exert mainly shear stress, studies of blood damage are based on shear flow experiments. However, measurements and simulations of cell and protein deformation show normal and shear stresses deform, and potentially damage, cells and proteins differently. The aim was to use computational fluid dynamics to assess the prevalence of normal stress, in comparison with shear stress, in rotary ventricular assist devices. Our calculations showed normal stresses do occur in rotary ventricular assist devices: the fluid volumes experiencing normal stress above 10 Pa were 0.011 mL (0.092%) and 0.027 mL (0.39%) for the HeartWare HVAD and HeartMate II (HMII), and normal stresses over 100 Pa were present. However, the shear stress volumes were up to two orders of magnitude larger than the normal stress volumes. Considering thresholds for red blood cell and von Willebrand factor deformation by normal and shear stresses, the fluid volumes causing deformation by normal stress were between 2.5 and 5 times the size of those causing deformation by shear stress. The exposure times to the individual normal stress deformation regions were around 1 ms. The results clearly show, for the first time, that while blood within rotary ventricular assist devices experiences more shear stress at much higher magnitudes as compared with normal stress, there is sufficient normal stress exposure present to cause deformation of, and potentially damage to, the blood components. This study is the first to quantify the fluid stress components in real blood contacting devices.

Experimental measurement and numerical modelling of dye washout for investigation of blood residence time in ventricular assist devices.

Ventricular assist devices have become the standard therapy for end-stage heart failure. However, their use is still associated with severe adverse events related to the damage done to the blood by fluid dynamic stresses. This damage relates to both the stress magnitude and the length of time the blood is exposed to that stress. We created a dye washout technique which combines experimental and numerical approaches to measure the washout times of ventricular assist devices. The technique was used to investigate washout characteristics of three commercially available and clinically used ventricular assist devices: the CentriMag, HVAD and HeartMate II. The time taken to reach 5% dye concentration at the outlet (T05) was used as an indicator of the total residence time. At a typical level of cardiac support, 5 L/min and 100 mmHg, T05 was 0.93, 0.28 and 0.16 s for CentriMag, HVAD and HeartMate II, respectively, and increased to 5.06, 1.64 and 0.96 s for reduced cardiac support of 1 L/min. Regional variations in washout characteristics are described in this article. While the volume of the flow domain plays a large role in the differences in T05 between the ventricular assist devices, after standardising for ventricular assist device volume, the secondary flow path was found to increase T05 by 35%. The results explain quantitatively, for the first time, why the CentriMag, which exerts low shear stress magnitude, has still been found to cause acquired von Willebrand Syndrome in patients.

Ultrasound imaging velocimetry with interleaved images for improved pulsatile arterial flow measurements: a new correction method, experimental and in vivo validation.

Blood velocity measurements are important in physiological science and clinical diagnosis. Doppler ultrasound is the most commonly used method but can only measure one velocity component. Ultrasound imaging velocimetry (UIV) is a promising technique capable of measuring two velocity components; however, there is a limit on the maximum velocity that can be measured with conventional hardware which results from the way images are acquired by sweeping the ultrasound beam across the field of view. Interleaved UIV is an extension of UIV in which two image frames are acquired concurrently, allowing the effective interframe separation time to be reduced and therefore increasing the maximum velocity that can be measured. The sweeping of the ultrasound beam across the image results in a systematic error which must be corrected: in this work, we derived and implemented a new velocity correction method which accounts for acceleration of the scatterers. We then, for the first time, assessed the performance of interleaved UIV for measuring pulsatile arterial velocities by measuring flows in phantoms and in vivo and comparing the results with spectral Doppler ultrasound and transit-time flow probe data. The velocity and flow rate in the phantom agreed within 5-10% of peak velocity, and 2-9% of peak flow, respectively, and in vivo the velocity difference was 9% of peak velocity. The maximum velocity measured was 1.8 m s-1, the highest velocity reported with UIV. This will allow flows in diseased arteries to be investigated and so has the potential to increase diagnostic accuracy and enable new vascular research.

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 Darcy's 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.

Ultrasound imaging velocimetry: effect of beam sweeping on velocity estimation.

As an emerging flow-mapping tool that can penetrate deep into optically opaque media such as human tissue, ultrasound imaging velocimetry has promise in various clinical applications. Previous studies have shown that errors occur in velocity estimation, but the causes have not been well characterised. In this study, the error in velocity estimation resulting from ultrasound beam sweeping in image acquisition is quantitatively investigated. The effects on velocity estimation of the speed and direction of beam sweeping relative to those of the flow are studied through simulation and experiment. The results indicate that a relative error in velocity estimation of up to 20% can be expected. Correction methods to reduce the errors under steady flow conditions are proposed and evaluated. Errors in flow estimation under unsteady flow are discussed.

Comparison and experimental validation of fluid dynamic numerical models for a clinical ventricular assist device.

With the recent advances in computer technology, computational fluid dynamics (CFDs) has become an important tool to design and improve blood-contacting artificial organs, and to study the device-induced blood damage. Commercial CFD software packages are readily available, and multiple CFD models are provided by CFD software developers. However, the best approach of using CFD effectively to characterize fluid flow and to predict blood damage in these medical devices remains debatable. This study aimed to compare these CFD models and provide useful information on the accuracy of each model in modeling blood flow in circulatory assist devices. The laminar and five turbulence models (Spalart-Allmaras, k-ε (k-epsilon), k-ω (k-omega), SST [Menter's Shear Stress Transport], and Reynolds Stress) were implemented to predict blood flow in a clinically used circulatory assist device, the CentriMag centrifugal blood pump. In parallel, a transparent replica of the CentriMag pump was constructed and selected views of the flow fields were measured with digital particle image velocimetry (DPIV). CFD results were compared with the DPIV experimental results. Compared with the experiment, all the selected CFD models predicted the flow pattern fairly well except the area of the outlet. However, quantitatively, the laminar model results were the most deviated from the experimental data. On the other hand, k-ε renormalization group theory models and Reynolds Stress model are the most accurate. In conclusion, for the circulatory assist devices, turbulence models provide more accurate results than the laminar model. Among the selected turbulence models, k-ε and Reynolds Stress Method models are recommended.

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.

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(α)τ(ß)) 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 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.

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.