FDA Benchmark Medical Device Flow Models for CFD Validation

Patient-specific prediction of arterial wall elasticity using medical image-informed in-silico simulations

Limitations in clinical cardiovascular research have driven the development of advanced simulations for patient-specific insights into arterial elasticity. However, uncertainties in model inputs, data resolution, and parameter estimation can compromise accuracy. Our research aimed to provide reliable estimates of the arterial wall elasticity non-invasively, where direct clinical measurement is difficult. By integrating patient-specific imaging with a simplified flow simulation model and uncertainty quantification, we sought to improve the reliability of these predictions as compared to the state-of-the-art. In a proof-of-concept study, we developed a simple area-averaged model of arterial hemodynamics, using Magnetic Resonance Angiogram (MRA)-derived geometries and input parameters based on the age, cuff blood pressure, and phase-contrast MRI data in five human subjects. This resulted in an in-silico model estimating the pressure and flow variations across the arterial-branches. Statistical uncertainties in the hemodynamic parameter predictions were quantified using non-intrusive Polynomial Chaos. Additionally, we developed a model to estimate the arterial elasticity by interlacing the results from fluid-structure interaction simulation for arterial hemodynamics with patient-specific clinical data. We found that the arterial elasticity values derived from our model, when used to predict the flowrates, closely matched the flow characteristics obtained from the patient-specific 4D flow MRI. The findings also showed zero or minimal positive/negative bias in our simulations, with no noticeable systematic error in predicting arterial elasticity values. Our results evidenced that accurate prediction of arterial wall elasticity is possible through use of an efficient simulation technique supplemented with clinically attainable imaging data. This has potential to predict cardiovascular-risk and guide individual patient management.

Hemodynamic evaluation of a novel double lumen cannula for left ventricle assist device system

BackgroundThe left ventricular assist device (LVAD) has been proven to be an effective therapy for providing temporary circulatory support. However, the use of this device can cause myocardial injury due to multiple insertions of various catheters.ObjectiveTherefore, this study aimed to evaluate the hemodynamic performance of a newly developed double-lumen catheter (DLC) for LVAD.MethodsTwo different LVAD DLC prototypes (a semi-circular and a concentric catheter) were designed based on the structure of venous DLC. Computational fluid dynamics (CFD) simulations were performed using the finite element method. The CFD results were confirmed through the testing of the 31 Fr prototype. The aorta is a large vessel with shear rates up to >300 s-1 and we used a reasonable approximation to model blood as a Newtonian fluid.ResultsAt a flow rate of 5 L/min, the semi-circular prototype achieved an infusion pressure of 74.68 mmHg, while the concentric prototype achieved an infusion pressure of 46.11 mmHg. The CFD results matched the experimental results with a mean percentage error of less than 7%. The peak wall shear stress in the semi-circular prototype (717.5 Pa) was higher than the hemolysis threshold (400 Pa), which could cause blood damage, and it also had a higher hemolysis index compared to concentric prototype. Moreover, both prototypes exhibited areas of blood stagnation and recirculation, suggesting a possible risk of thrombosis.ConclusionBoth prototypes of the LVAD DLC demonstrated similar blood flow rates. The semi-circular prototype showed superior infusion pressure compared to the concentric prototype, but had poorer hemolysis performance. However, the potential risk of thrombosis for both still exists. Therefore, further in vivo experiments are necessary to verify the safety and effectiveness of the LVAD DLC.

Structural and functional stenosis of the upper airway in Crouzon syndrome patients: A computational fluid dynamics analysis

OBJECTIVES

This study aimed to simulate the aerodynamics and to identify the spatial correlation between anatomical and functional stenoses in Crouzon syndrome patients.

METHODS

Six patients of Crouzon syndrome were included. Computational fluid dynamics (CFD) was utilized to simulate airflow dynamics, and characteristics, including the velocity, pressure intensity, wall shear stress, airflow resistance and streamline, were extracted for quantitative analysis both in overall and regionally. Structural stenosis was defined at the minimum cross-sectional area, while functional stenosis was identified at the point of maximum airflow velocity. The spatial distances between the Frankfurt plane and structural/functional stenosis were calculated and compared.

RESULTS

Structural stenosis occurred in the palatopharynx, while the highest inspiratory resistance and peak airflow velocity during expiration identified the glossopharynx as the functional stenosis site. A steep increase in negative pressure and a significant increase in wall shear stress could be observed surrounding the functional stenosis. The intensity and diffusion range of wall shear stress are positively correlated with age. Notably, the functional stenosis was consistently 5 mm below the structural stenosis (P < 0.05).

CONCLUSIONS

CFD effectively visualized both overall and regional aerodynamics of Crouzon syndrome, providing a novel method for functional airway evaluation. The spatial distributions of structural and functional stenoses did not strictly correspond; the structural stenosis was located on the palatopharynx, while the functional stenosis was on the glossopharynx. The wall shear stress worsens pathologically with age, aggravating functional stenosis to structural stenosis. Therefore, functional stenosis should also be addressed in airway management to ensure therapeutic effectiveness.

Computational fluid dynamics simulating of the FDA benchmark blood pump with different coefficient sets and scaler shear stress models used in the power-law hemolysis model

Hemolysis is the most important issue to consider in the design and optimization of blood-contacting devices. Although the use of Computational Fluid Dynamics (CFD) in hemolysis prediction studies provides convenience and has promising potential, it is an extremely challenging process. Hemolysis predictions with CFD depend on the mesh, implementation method, coefficient set, and scalar-shear-stress model. To this end, an attempt was made to find the combination that would provide the most accurate result in hemolysis prediction with the commonly cited power-law based hemolysis model. In the hemolysis predictions conducted using CFD on the Food and Drug Administration (FDA) benchmark blood pump, 3 different scalar-shear-stress models, and 5 different coefficient sets with the power-law based hemolysis model were used. Also, a mesh independence test based on hemolysis and pressure head was performed. The pressure head results of CFD simulations were compared with published pressure head of the FDA benchmark blood pump and a good agreement was observed. In addition, results of CFD-hemolysis predictions which are conducted with scalar-shear-stress model and coefficient set combinations were compared with experimental hemolysis data at three operating conditions such as 6-7 L/min flow rates at 3500 rpm rotational speeds and 6 L/min at 2500 rpm. One of the combinations of the scalar-shear-stress model and the coefficient set was found to be within the error limits of the experimental measurements, while all other combinations overestimated hemolysis.

Empirical and Computational Evaluation of Hemolysis in a Microfluidic Extracorporeal Membrane Oxygenator Prototype

Microfluidic devices promise to overcome the limitations of conventional hemodialysis and oxygenation technologies by incorporating novel membranes with ultra-high permeability into portable devices with low blood volume. However, the characteristically small dimensions of these devices contribute to both non-physiologic shear that could damage blood components and laminar flow that inhibits transport. While many studies have been performed to empirically and computationally study hemolysis in medical devices, such as valves and blood pumps, little is known about blood damage in microfluidic devices. In this study, four variants of a representative microfluidic membrane-based oxygenator and two controls (positive and negative) are introduced, and computational models are used to predict hemolysis. The simulations were performed in ANSYS Fluent for nine shear stress-based parameter sets for the power law hemolysis model. We found that three of the nine tested parameters overpredict (5 to 10×) hemolysis compared to empirical experiments. However, three parameter sets demonstrated higher predictive accuracy for hemolysis values in devices characterized by low shear conditions, while another three parameter sets exhibited better performance for devices operating under higher shear conditions. Empirical testing of the devices in a recirculating loop revealed levels of hemolysis significantly lower (<2 ppm) than the hemolysis ranges observed in conventional oxygenators (>10 ppm). Evaluating the model's ability to predict hemolysis across diverse shearing conditions, both through empirical experiments and computational validation, will provide valuable insights for future micro ECMO device development by directly relating geometric and shear stress with hemolysis levels. We propose that, with an informed selection of hemolysis parameters based on the shear ranges of the test device, computational modeling can complement empirical testing in the development of novel high-flow blood-contacting microfluidic devices, allowing for a more efficient iterative design process. Furthermore, the low device-induced hemolysis measured in our study at physiologically relevant flow rates is promising for the future development of microfluidic oxygenators and dialyzers.

Design and Analysis of a Polymeric Left Ventricular Simulator via Computational Modelling

Preclinical testing of medical devices is an essential step in the product life cycle, whereas testing of cardiovascular implants requires specialised testbeds or numerical simulations using computer software Ansys 2016. Existing test setups used to evaluate physiological scenarios and test cardiac implants such as mock circulatory systems or isolated beating heart platforms are driven by sophisticated hardware which comes at a high cost or raises ethical concerns. On the other hand, computational methods used to simulate blood flow in the cardiovascular system may be simplified or computationally expensive. Therefore, there is a need for low-cost, relatively simple and efficient test beds that can provide realistic conditions to simulate physiological scenarios and evaluate cardiovascular devices. In this study, the concept design of a novel left ventricular simulator made of latex rubber and actuated by pneumatic artificial muscles is presented. The designed left ventricular simulator is geometrically similar to a native left ventricle, whereas the basal diameter and long axis length are within an anatomical range. Finite element simulations evaluating left ventricular twisting and shortening predicted that the designed left ventricular simulator rotates approximately 17 degrees at the apex and the long axis shortens around 11 mm. Experimental results showed that the twist angle is 18 degrees and the left ventricular simulator shortens 5 mm. Twist angles and long axis shortening as in a native left ventricle show it is capable of functioning like a native left ventricle and simulating a variety of scenarios, and therefore has the potential to be used as a test platform.

Optimization of a centrifugal blood pump designed using an industrial method through experimental and numerical study

With improved treatment of coronary artery disease, more patients are surviving until heart failure occurs. This leads to an increase in patients needing devices for struggling with heart failure. Ventricular assist devices are known as the mainstay of these devices. This study aimed to design a centrifugal pump as a ventricular assist device. In order to design the pump, firstly, the geometrical parameters of the pump, including the gap distance, blade height, and position of the outlet relative to the blade, were investigated. Finally, the selected configuration, which had all the appropriate characteristics, both hydraulically and physiologically, was used for the rest of the study. The study of the blade, as the main component in energy transfer to the blood, in a centrifugal pump, has been considered in the present study. In this regard, the point-to-point design method, which is used in industrial applications, was implemented. The designer chooses the relationship between the blade angles at each radius in the point-to-point method. The present study selected logarithmic and second-order relations for designing the blade's profile. In total, 58 blades were examined in this study, which differed regarding blade inlet and outlet angles and the relationship between angle and radial position. ANSYS CFX 17.0 software was utilized to simulate blades' performances, and a benchmark pump provided by the US Food and Drug Administration (FDA) was used to validate the numerical simulations. Then, the selected impeller from the numerical investigation was manufactured, and its performance was compared experimentally with the FDA benchmark pump. A hydraulic test rig was also developed for experimental studies. The results showed that among the blades designed in this study, the blade with an input angle of 45° and an output angle of 55°, which is designed to implement a logarithmic relationship, has the best performance. The selected impeller configuration can increase the total head (at least by 20%) at different flow rates compared to the FDA pump.

Study on the optimal elastic modulus of flexible blades for right heart assist device supporting patients with single-ventricle physiologies

Background

Patients with single-ventricle physiologies continue to experience insufficient circulatory power after undergoing palliative surgeries. This paper proposed a right heart assist device equipped with flexible blades to provide circulatory assistance for these patients. The optimal elastic modulus of the flexible blades was investigated through numerical simulation.

Methods

A one-way fluid-structure interaction (FSI) simulation was employed to study the deformation of flexible blades during rotation and its impact on device performance. The process began with a computational fluid dynamics (CFD) simulation to calculate the blood pressure rise and the pressure on the blades' surface. Subsequently, these pressure data were exported for finite element analysis (FEA) to compute the deformation of the blades. The fluid domain was then recreated based on the deformed blades' shape. Iterative CFD and FEA simulations were performed until both the blood pressure rise and the blades' shape stabilized. The blood pressure rise, hemolysis risk, and thrombosis risk corresponding to blades with different elastic moduli were exhaustively evaluated to determine the optimal elastic modulus.

Results

Except for the case at 8,000 rpm with a blade elastic modulus of 40 MPa, the pressure rise associated with flexible blades within the studied range (rotational speeds of 4,000 rpm and 8,000 rpm, elastic modulus between 10 MPa and 200 MPa) was lower than that of rigid blades. It was observed that the pressure rise corresponding to flexible blades increased as the elastic modulus increased. Additionally, no significant difference was found in the hemolysis risk and thrombus risk between flexible blades of various elastic moduli and rigid blades.

Conclusion

Except for one specific case, deformation of the flexible blades within the studied range led to a decrease in the impeller's functionality. Notably, rotational speed had a more significant impact on hemolysis risk and thrombus risk compared to blade deformation. After a comprehensive analysis of blade compressibility, blood pressure rise, hemolysis risk, and thrombus risk, the optimal elastic modulus for the flexible blades was determined to be between 40 MPa and 50 MPa.

Incorporating Unresolved Stresses in Blood Damage Modeling: Energy Dissipation More Accurate Than Reynolds Stress Formulation

OBJECTIVE

Reynolds Averaged Navier Stokes (RANS) models are often used as the basis for modeling blood damage in turbulent flows. To predict blood damage by turbulence stresses that are not resolved in RANS, a stress formulation that represents the corresponding scales is required. Here, we compare two commonly employed stress formulations: a scalar stress representation that uses Reynolds stresses as a surrogate for unresolved fluid stresses, and an effective stress formulation based on energy dissipation.

METHODS

We conducted unsteady RANS simulations of the CentriMag blood pump with three different closure models and a Large Eddy Simulation (LES) for reference. We implemented both stress representations in all models and compared the resulting total stress distributions in Eulerian and Lagrangian frameworks.

RESULTS

The Reynolds-stress-based approach overestimated the contribution of unresolved stresses in RANS, with differences between closure models of up to several orders of magnitude. With the dissipation-based approach, the total stresses predicted with RANS deviated by about 50% from the LES reference, which was more accurate than only considering resolved stresses.

CONCLUSION

The Reynolds-stress-based formulation proved unreliable for estimating scalar stresses in our RANS simulations, while the dissipation-based approach provided an accuracy improvement over simply neglecting unresolved stresses.

SIGNIFICANCE

Our results suggest that dissipation-based inclusion of unresolved stresses should be the preferred choice for blood damage modeling in RANS.

Systematic analysis of non-intrusive polynomial chaos expansion to determine rotary blood pump performance over the entire operating range

This study applies non-intrusive polynomial chaos expansion (NIPCE) surrogate modeling to analyze the performance of a rotary blood pump (RBP) across its operating range. We systematically investigate key parameters, including polynomial order, training data points, and data smoothness, while comparing them to test data. Using a polynomial order of 4 and a minimum of 20 training points, we successfully train a NIPCE model that accurately predicts pressure head and axial force within the specified operating point range ([0-5000] rpm and [0-7] l/min). We also assess the NIPCE model's ability to predict two-dimensional velocity data across the given range and find good overall agreement (mean absolute error = 0.1 m/s) with a test simulation under the same operating condition. Our approach extends current NIPCE modeling of RBPs by considering the entire operating range and providing validation guidelines. While acknowledging computational benefits, we emphasize the challenge of modeling discontinuous data and its relevance to clinically realistic operating points. We offer open access to our raw data and Python code, promoting reproducibility and accessibility within the scientific community. In conclusion, this study advances comprehensive NIPCE modeling of RBP performance and underlines how critically NIPCE parameters and rigorous validation affect results.

Next Generation Development of Hybrid Continuous Flow Pediatric Total Artificial Heart Technology: Design-Build-Test

To address the unmet clinical need for pediatric circulatory support, we are developing an operationally versatile, hybrid, continuous-flow, total artificial heart ("Dragon Heart"). This device integrates a magnetically levitated axial and centrifugal blood pump. Here, we utilized a validated axial flow pump, and we focused on the development of the centrifugal pump. A motor was integrated to drive the centrifugal pump, achieving 50% size reduction. The motor design was simulated by finite element analysis, and pump design improvement was attained by computational fluid dynamics. A prototype centrifugal pump was constructed from biocompatible 3D printed parts for the housing and machined metal parts for the drive system. Centrifugal prototype testing was conducted using water and then bovine blood. The fully combined device ( i.e. , axial pump nested inside of the centrifugal pump) was tested to ensure proper operation. We demonstrated the hydraulic performance of the two pumps operating in tandem, and we found that the centrifugal blood pump performance was not adversely impacted by the simultaneous operation of the axial blood pump. The current iteration of this design achieved a range of operation overlapping our target range. Future design iterations will further reduce size and incorporate complete and active magnetic levitation.

Current progress toward isogeometric modeling of the heart biophysics

In this paper, we review a powerful methodology to solve complex numerical simulations, known as isogeometric analysis, with a focus on applications to the biophysical modeling of the heart. We focus on the hemodynamics, modeling of the valves, cardiac tissue mechanics, and on the simulation of medical devices and treatments. For every topic, we provide an overview of the methods employed to solve the specific numerical issue entailed by the simulation. We try to cover the complete process, starting from the creation of the geometrical model up to the analysis and post-processing, highlighting the advantages and disadvantages of the methodology.

Resistance valves in circulatory loops have a significant impact on in vitro evaluation of blood damage caused by blood pumps: a computational study

Background: Hemolysis and its complications are major concerns during the clinical application of blood pumps. In-vitro circulatory testing loops have been employed as the key procedure to evaluate the hemolytic and thrombogenic performance of blood pumps during the development phase and before preclinical in-vivo animal studies. Except for the blood damage induced by the pump under test, blood damage induced by loop components such as the resistance valve may affect the accuracy, reproducibility, and intercomparability of test results. Methods: This study quantitatively investigated the impact of the resistance valve on in vitro evaluation of blood damage caused by blood pumps under different operating points. A series of idealized tubing models under the resistance valve with different openings were created. Three pumps - the FDA benchmark pump, the HeartMate 3 LVAD, and the CH-VAD - were involved in hypothetical tests. Eight operating points were chosen to cover a relatively wide spectrum of testing scenarios. Computational fluid dynamics (CFD) simulations of the tubing and pump models were conducted at the same operating points. Results and Conclusion: Overall, hemolysis and platelet activation induced by a typical resistance valve are equivalent to 17%-45% and 14%-60%, respectively, of those induced by the pump itself. Both ratios varied greatly with flow rate, valve opening and pump models. Differences in blood damage levels between different blood pumps or working conditions can be attenuated by up to 45%. Thus, hemolysis and platelet activation induced by the resistance valve significantly affect the accuracy of in-vitro hemocompatibility evaluations of blood pumps. A more accurate and credible method for hemocompatibility evaluations of blood pumps will benefit from these findings.

Numerical study on the performance of mixed flow blood pump with superhydrophobic surface

To meet the clinical status of the wide application of percutaneous mechanical circulatory support, this paper selects the mixed flow blood pump applied with superhydrophobic surface as the research object. The Navier slip model was used to simulate the slip characteristics of superhydrophobic surface, and the effects of the blade wrap angle and the superhydrophobic surface on the performance of the mixed flow blood pump are studied by numerical simulation. The results show that (1) considering the head, hydraulic efficiency, and hemolysis index of the blood pump, the optimal value of the blade wrap angle of the mixed flow blood pump in this paper is 60°. (2) The hydraulic efficiency of the blood pump with superhydrophobic surface is improved, and the maximum growth rate is about 13.9%; superhydrophobic surface can reduce the hemolysis index of blood pump under various working conditions, and the maximum reduction rate of hemolysis index of blood pump is 22.9%. (3) The variation trends of blood pump head, hydraulic efficiency, and hemolysis index with the increased rotating speed before and after setting superhydrophobic slip boundary conditions are the same as their original variation trends.

In-Vitro Flow Validation of Third-Generation Ventricular Assist Devices: Feasibility and Challenges

Computational fluid dynamics (CFD) is a powerful tool for the in-silico evaluation of rotodynamic blood pumps (RBPs). Corresponding validation, however, is typically restricted to easily accessible, global flow quantities. This study showcased the HeartMate 3 (HM3) to identify feasibility and challenges of enhanced in-vitro validation in third-generation RBPs. To enable high-precision acquisition of impeller torques and grant access for optical flow measurements, the HM3 testbench geometry was geometrically modified. These modifications were reproduced in silico , and global flow computations validated along 15 operating conditions. The globally validated flow in the testbench geometry was compared with CFD-simulated flows in the original geometry to assess the impact of the necessary modifications on global and local hydraulic properties. Global hydraulic properties in the testbench geometry were successfully validated (pressure head: r = 0.999, root mean square error [RMSE] = 2.92 mmHg; torque: r = 0.996, RMSE = 0.134 mNm). In-silico comparison with the original geometry demonstrated good agreement ( r > 0.999, relative errors < 11.97%) of global hydraulic properties. Local hydraulic properties (errors up to 81.78%) and hemocopatibility predictions (deviations up to 21.03%), however, were substantially affected by the geometric modifications. Transferability of local flow measures derived on advanced in-vitro testbenches toward original pump designs is challenged by significant local effects associated with the necessary geometrical modifications.

Beyond CFD: Emerging methodologies for predictive simulation in cardiovascular health and disease

Physics-based computational models of the cardiovascular system are increasingly used to simulate hemodynamics, tissue mechanics, and physiology in evolving healthy and diseased states. While predictive models using computational fluid dynamics (CFD) originated primarily for use in surgical planning, their application now extends well beyond this purpose. In this review, we describe an increasingly wide range of modeling applications aimed at uncovering fundamental mechanisms of disease progression and development, performing model-guided design, and generating testable hypotheses to drive targeted experiments. Increasingly, models are incorporating multiple physical processes spanning a wide range of time and length scales in the heart and vasculature. With these expanded capabilities, clinical adoption of patient-specific modeling in congenital and acquired cardiovascular disease is also increasing, impacting clinical care and treatment decisions in complex congenital heart disease, coronary artery disease, vascular surgery, pulmonary artery disease, and medical device design. In support of these efforts, we discuss recent advances in modeling methodology, which are most impactful when driven by clinical needs. We describe pivotal recent developments in image processing, fluid-structure interaction, modeling under uncertainty, and reduced order modeling to enable simulations in clinically relevant timeframes. In all these areas, we argue that traditional CFD alone is insufficient to tackle increasingly complex clinical and biological problems across scales and systems. Rather, CFD should be coupled with appropriate multiscale biological, physical, and physiological models needed to produce comprehensive, impactful models of mechanobiological systems and complex clinical scenarios. With this perspective, we finally outline open problems and future challenges in the field.

Modeling flow in an in vitro anatomical cerebrovascular model with experimental validation

Acute ischemic stroke (AIS) is a leading cause of mortality that occurs when an embolus becomes lodged in the cerebral vasculature and obstructs blood flow in the brain. The severity of AIS is determined by the location and how extensively emboli become lodged, which are dictated in large part by the cerebral flow and the dynamics of embolus migration which are difficult to measure in vivo in AIS patients. Computational fluid dynamics (CFD) can be used to predict the patient-specific hemodynamics and embolus migration and lodging in the cerebral vasculature to better understand the underlying mechanics of AIS. To be relied upon, however, the computational simulations must be verified and validated. In this study, a realistic in vitro experimental model and a corresponding computational model of the cerebral vasculature are established that can be used to investigate flow and embolus migration and lodging in the brain. First, the in vitro anatomical model is described, including how the flow distribution in the model is tuned to match physiological measurements from the literature. Measurements of pressure and flow rate for both normal and stroke conditions were acquired and corresponding CFD simulations were performed and compared with the experiments to validate the flow predictions. Overall, the CFD simulations were in relatively close agreement with the experiments, to within ±7% of the mean experimental data with many of the CFD predictions within the uncertainty of the experimental measurement. This work provides an in vitro benchmark data set for flow in a realistic cerebrovascular model and is a first step towards validating a computational model of AIS.

Influence of Inlet Boundary Conditions on the Prediction of Flow Field and Hemolysis in Blood Pumps Using Large-Eddy Simulation

Inlet boundary conditions (BC) are one of the uncertainties which may influence the prediction of flow field and hemolysis in blood pumps. This study investigated the influence of inlet BC, including the length of inlet pipe, type of inlet BC (mass flow rate or experimental velocity profile) and turbulent intensity (no perturbation, 5%, 10%, 20%) on the prediction of flow field and hemolysis of a benchmark centrifugal blood pump (the FDA blood pump) and a commercial axial blood pump (Heartmate II), using large-eddy simulation. The results show that the influence of boundary conditions on integral pump performance metrics, including pressure head and hemolysis, is negligible. The influence on local flow structures, such as velocity distributions, mainly existed in the inlet. For the centrifugal FDA blood pump, the influence of type of inlet BC and inlet position on velocity distributions can also be observed at the diffuser. Overall, the effects of position of inlet and type of inlet BC need to be considered if local flow structures are the focus, while the influence of turbulent intensity is negligible and need not be accounted for during numerical simulations of blood pumps.

Modeling Flow in an In Vitro Anatomical Cerebrovascular Model with Experimental Validation

Acute ischemic stroke (AIS) is a leading cause of mortality that occurs when an embolus becomes lodged in the cerebral vasculature and obstructs blood flow in the brain. The severity of AIS is determined by the location and how extensively emboli become lodged, which are dictated in large part by the cerebral flow and the dynamics of embolus migration which are difficult to measure in vivo in AIS patients. Computational fluid dynamics (CFD) can be used to predict the patient-specific hemodynamics and embolus migration and lodging in the cerebral vasculature to better understand the underlying mechanics of AIS. To be relied upon, however, the computational simulations must be verified and validated. In this study, a realistic in vitro experimental model and a corresponding computational model of the cerebral vasculature are established that can be used to investigate flow and embolus migration and lodging in the brain. First, the in vitro anatomical model is described, including how the flow distribution in the model is tuned to match physiological measurements from the literature. Measurements of pressure and flow rate for both normal and stroke conditions were acquired and corresponding CFD simulations were performed and compared with the experiments to validate the flow predictions. Overall, the CFD simulations were in relatively close agreement with the experiments, to within ±7% of the mean experimental data with many of the CFD predictions within the uncertainty of the experimental measurement. This work provides an in vitro benchmark data set for flow in a realistic cerebrovascular model and is a first step towards validating a computational model of AIS.

Hemocompatibility and hemodynamic comparison of two centrifugal LVADs: HVAD and HeartMate3

Mechanical circulatory support using ventricular assist devices is a common technique for treating patients suffering from advanced heart failure. The latest generation of devices is characterized by centrifugal turbopumps which employ magnetic levitation bearings to ensure a gap clearance between moving and static parts. Despite the increasing use of these devices as a destination therapy, several long-term complications still exist regarding their hemocompatibility. The blood damage associated with different pump designs has been investigated profoundly in the literature, while the hemodynamic performance has been hardly considered. This work presents a novel comparison between the two main devices of the latest generation-HVAD and HM3-from both perspectives, hemodynamic performance and blood damage. Computational fluid dynamics simulations are performed to model the considered LVADs, and computational results are compared to experimental measurements of pressure head to validate the model. Enhanced performance and hemocompatibility are detected for HM3 owing to its design incorporating more conventional blades and larger gap clearances.

Mitigation effect of cell exclusion on blood damage in spiral groove bearings

Cell exclusion in spiral groove bearing (SGB) excludes red blood cells from high shear regions in the bearing gaps and potentially reduce haemolysis in rotary blood pumps. However, this mechanobiological phenomenon has been observed in ultra-low blood haematocrit only, whether it can mitigate blood damage in a clinically-relevant blood haematocrit remains unknown. This study examined whether cell exclusion in a SGB alters haemolysis and/or high-molecular-weight von Willebrand factor (HMW vWF) multimer degradation. Citrated human blood was adjusted to 35 % haematocrit and exposed to a SGB (n = 6) and grooveless disc (n = 3, as a non-cell exclusion control) incorporated into a custom-built Couette test rig operating at 2000RPM for an hour; shearing gaps were 20, 30, and 40 μm. Haemolysis was assessed via spectrophotometry and HMW vWF multimer degradation was detected with gel electrophoresis and immunoblotting. Haemolysis caused by the SGB at gaps of 20, 30 and 40 μm were 10.6 ± 3.3, 9.6 ± 2.7 and 10.5 ± 3.9 mg/dL.hr compared to 23.3 ± 2.6, 12.8 ± 3.2, 9.8 ± 1.8 mg/dL.hr by grooveless disc. At the same shearing gap of 20 µm, there was a significant reduced in haemolysis (P = 0.0001) and better preserved in HMW vWF multimers (p < 0.05) when compared SGB to grooveless disc. The reduction in blood damage in the SGB compared to grooveless disc is indicative of cell exclusion occurred at the gap of 20 µm. This is the first experimental study to demonstrate that cell exclusion in a SGB mitigates the shear-induced blood damage in a clinically-relevant blood haematocrit of 35 %, which can be potentially utilised in future blood pump design.

Results of the Interlaboratory Computational Fluid Dynamics Study of the FDA Benchmark Blood Pump

Computational fluid dynamics (CFD) is widely used to simulate blood-contacting medical devices. To be relied upon to inform high-risk decision making, however, model credibility should be demonstrated through validation. To provide robust data sets for validation, researchers at the FDA and collaborators developed two benchmark medical device flow models: a nozzle and a centrifugal blood pump. Experimental measurements of the flow fields and hemolysis were acquired using each model. Concurrently, separate open interlaboratory CFD studies were performed in which participants from around the world, who were blinded to the measurements, submitted CFD predictions of each benchmark model. In this study, we report the results of the interlaboratory CFD study of the FDA benchmark blood pump. We analyze the results of 24 CFD submissions using a wide range of different flow solvers, methods, and modeling parameters. To assess the accuracy of the CFD predictions, we compare the results with experimental measurements of three quantities of interest (pressure head, velocity field, and hemolysis) at different pump operating conditions. We also investigate the influence of different CFD methods and modeling choices used by the participants. Our analyses reveal that, while a number of CFD submissions accurately predicted the pump performance for individual cases, no single participant was able to accurately predict all quantities of interest across all conditions. Several participants accurately predicted the pressure head at all conditions and the velocity field in all but one or two cases. Only one of the eight participants who submitted hemolysis results accurately predicted absolute plasma free hemoglobin levels at a majority of the conditions, though most participants were successful at predicting relative hemolysis levels between conditions. Overall, this study highlights the need to validate CFD modeling of rotary blood pumps across the entire range of operating conditions and for all quantities of interest, as some operating conditions and regions (e.g., the pump diffuser) are more challenging to accurately predict than others. All quantities of interest should be validated because, as shown here, it is possible to accurately predict hemolysis despite having relatively inaccurate predictions of the flow field.

A CFD study on the interplay of torsion and vortex guidance by the mitral valve on the left ventricular wash-out making use of overset meshes (Chimera technique)

Cardiovascular disease often occurs with silent and gradual alterations of cardiac blood flow that can lead to the onset of chronic pathological conditions. Image-based patient-specific Computational Fluid Dynamics (CFD) models allow for an extensive quantification of the flow field beyond the direct capabilities of medical imaging techniques that could support the clinicians in the early diagnosis, follow-up, and treatment planning of patients. Nonetheless, the large and impulsive kinematics of the left ventricle (LV) and the mitral valve (MV) pose relevant modeling challenges. Arbitrary Lagrangian-Eulerian (ALE) based computational fluid dynamics (CFD) methods struggle with the complex 3D mesh handling of rapidly moving valve leaflets within the left ventricle (LV). We, therefore, developed a Chimera-based (overset meshing) method to build a patient-specific 3D CFD model of the beating LV which includes a patient-inspired kinematic model of the mitral valve (LVMV). Simulations were performed with and without torsion. In addition, to evaluate how the intracardiac LV flow is impacted by the MV leaflet kinematics, a third version of the model without the MV was generated (LV with torsion). For all model versions, six cardiac cycles were simulated. All simulations demonstrated cycle-to-cycle variations that persisted after six cycles but were albeit marginal in terms of the magnitude of standard deviation of velocity and vorticity which may be related to the dissipative nature of the numerical scheme used. The MV was found to have a crucial role in the development of the intraventricular flow by enhancing the direct flow, the apical washout, and the propagation of the inlet jet towards the apical region. Consequently, the MV is an essential feature in the patient-specific CFD modeling of the LV. The impact of torsion was marginal on velocity, vorticity, wall shear stress, and energy loss, whereas it resulted to be significant in the evaluation of particle residence times. Therefore, including torsion could be considered in patient-specific CFD models of the LV, particularly when aiming to study stasis and residence time. We conclude that, despite some technical limitations encountered, the Chimera technique is a promising alternative for ALE methods for 3D CFD models of the heart that include the motion of valve leaflets.

Smartphone-based particle image velocimetry for cardiovascular flows applications: A focus on coronary arteries

An experimental set-up is presented for the in vitro characterization of the fluid dynamics in personalized phantoms of healthy and stenosed coronary arteries. The proposed set-up was fine-tuned with the aim of obtaining a compact, flexible, low-cost test-bench for biomedical applications. Technically, velocity vector fields were measured adopting a so-called smart-PIV approach, consisting of a smartphone camera and a low-power continuous laser (30 mW). Experiments were conducted in realistic healthy and stenosed 3D-printed phantoms of left anterior descending coronary artery reconstructed from angiographic images. Time resolved image acquisition was made possible by the combination of the image acquisition frame rate of last generation commercial smartphones and the flow regimes characterizing coronary hemodynamics (velocities in the order of 10 cm/s). Different flow regimes (Reynolds numbers ranging from 20 to 200) were analyzed. The smart-PIV approach was able to provide both qualitative flow visualizations and quantitative results. A comparison between smart-PIV and conventional PIV (i.e., the gold-standard experimental technique for bioflows characterization) measurements showed a good agreement in the measured velocity vector fields for both the healthy and the stenosed coronary phantoms. Displacement errors and uncertainties, estimated by applying the particle disparity method, confirmed the soundness of the proposed smart-PIV approach, as their values fell within the same range for both smart and conventional PIV measured data (≈5% for the normalized estimated displacement error and below 1.2 pixels for displacement uncertainty). In conclusion, smart-PIV represents an easy-to-implement, low-cost methodology for obtaining an adequately robust experimental characterization of cardiovascular flows. The proposed approach, to be intended as a proof of concept, candidates to become an easy-to-handle test bench suitable for use also outside of research labs, e.g., for educational or industrial purposes, or as first-line investigation to direct and guide subsequent conventional PIV measurements.

Flow simulation-based particle swarm optimization for developing improved hemolysis models

The improvement and development of blood-contacting devices, such as mechanical circulatory support systems, is a life saving endeavor. These devices must be designed in such a way that they ensure the highest hemocompatibility. Therefore, in-silico trials (flow simulations) offer a quick and cost-effective way to analyze and optimize the hemocompatibility and performance of medical devices. In that regard, the prediction of blood trauma, such as hemolysis, is the key element to ensure the hemocompatibility of a device. But, despite decades of research related to numerical hemolysis models, their accuracy and reliability leaves much to be desired. This study proposes a novel optimization path, which is capable of improving existing models and aid in the development of future hemolysis models. First, flow simulations of three, turbulent blood flow test cases (capillary tube, FDA nozzle, FDA pump) were performed and hemolysis was numerically predicted by the widely-applied stress-based hemolysis models. Afterward, a multiple-objective particles swarm optimization (MOPSO) was performed to tie the physiological stresses of the simulated flow field to the measured hemolysis using an equivalent of over one million numerically determined hemolysis predictions. The results show that our optimization is capable of improving upon existing hemolysis models. However, it also unveils some deficiencies and limits of hemolysis prediction with stress-based models, which will need to be addressed in order to improve its reliability.

Examining the universality of the hemolysis power law model from simulations of the FDA nozzle using calibrated model coefficients

Computational fluid dynamics (CFD) is widely used to predict mechanical hemolysis in medical devices. The most popular hemolysis model is the stress-based power law model that is based on an empirical correlation between hemoglobin release from red blood cells (RBCs) and the magnitude of flow-induced stress and exposure time. Empirical coefficients are traditionally calibrated using data from experiments in simplified Couette-type blood-shearing devices with uniform-shear laminar flow and well-defined exposure times. Use of such idealized coefficients in simulations of real medical devices with complex hemodynamics is thought to be a primary reason for the historical inaccuracy of absolute hemolysis predictions using the power law model. Craven et al. (Biomech Model Mechanobiol 18:1005-1030, 2019) recently developed a CFD-based Kriging surrogate modeling approach for calibrating empirical coefficients in real devices that could potentially be used to more accurately predict absolute hemolysis. In this study, we use the FDA benchmark nozzle to investigate whether utilizing such calibrated coefficients improves the predictive accuracy of the standard Eulerian power law model. We first demonstrate the credibility of our CFD flow simulations by comparing with particle image velocimetry measurements. We then perform hemolysis simulations and compare the results with in vitro experiments. Importantly, the simulations use coefficients calibrated for the flow of a suspension of bovine RBCs through a small capillary tube, which is relatively comparable to the flow of bovine blood through the FDA nozzle. The results show that the CFD predictions of relative hemolysis in the FDA nozzle are reasonably accurate. The absolute predictions are, however, highly inaccurate with modified index of hemolysis values from CFD in error by roughly three orders of magnitude compared with the experiments, despite using calibrated model coefficients from a relatively similar geometry. We rigorously examine the reasons for the inaccuracy that include differences in the flow conditions in the hemolytic regions of each device and the lack of universality of the hemolysis power law model that is entirely empirical. Thus, while the capability to predict relative hemolysis is valuable for product development, further improvements are needed before the power law model can be relied upon to accurately predict the absolute hemolytic potential of a medical device.

Numerical study on the performance of centrifugal blood pump with superhydrophobic surface

AIM

In order to reduce the blood damage of an artificial heart pump and optimize its hydraulic performance, a centrifugal blood pump with superhydrophobic characteristics is proposed in this study.

METHODS

To study the influence of superhydrophobic surface characteristics on the performance of centrifugal blood pumps, the Navier slip model is used to simulate the slip characteristics of superhydrophobic surfaces, which is realized by the user defined function of ANSYS fluent. The user defined functions with different values of slip length are verified by two benchmark solutions of laminar flow and turbulence in the pipeline. The blood pump model adopts the designed centrifugal blood pump, and its head, hydraulic efficiency and hemolysis index are calculated. The Navier slip boundary condition (a constant slip-length of 50 μm) is applied to the walls of the blood pump impeller and a volute at different positions, and the influence of the superhydrophobic surface on the performance of the blood pump at the design point Q = 6 L/min was compared and analyzed.

RESULTS

The results show that the centrifugal blood pump model used in this paper has good blood compatibility and meets the design requirements; the superhydrophobic surface can significantly reduce the scalar shear stress in the blood pump. At the design point, when the slip length is 50 μm, the mass-average scalar shear stress in the impeller area and the volute area reduction rate is about 5.9%, the hydraulic efficiency growth rate is about 3.8%, the hemolysis index reduction rate is about 18.4%, and the pressure head changes little with a growth rate of 0.3%.

CONCLUSIONS

Centrifugal blood pumps with superhydrophobic surfaces can improve the efficiency of blood pumps and reduce hemolysis. Based on these encouraging results, vitro investigations for actual blood damage would be practicable.

Patient-specific changes in aortic hemodynamics is associated with thrombotic risk after fenestrated endovascular aneurysm repair with large diameter endografts

Background

The durability of fenestrated endovascular aneurysm repair (fEVAR) has been threatened by thrombotic complications. In the present study, we used patient-specific computational fluid dynamic (CFD) simulation to investigate the effect of the endograft diameter on hemodynamics after fEVAR and explore the hypothesis that diameter-dependent alterations in aortic hemodynamics can predict for thrombotic events.

Methods

A single-institutional retrospective study was performed of patients who had undergone fEVAR for juxtarenal aortic aneurysms. The patients were stratified into large diameter (34-36 mm) and small diameter (24-26 mm) endograft groups. Patient-specific CFD simulations were performed using three-dimensional paravisceral aortic models created from computed tomographic images with allometrically scaled boundary conditions. Aortic time-averaged wall shear stress (TAWSS) and residence time (RT) were computed and correlated with future thrombotic complications (eg, renal stent occlusion, development of significant intraluminal graft thrombus).

Results

A total of 36 patients (14 with a small endograft and 22 with a large endograft) were included in the present study. The patients treated with large endografts had experienced a higher incidence of thrombotic complications compared with small endografts (45.5% vs 7.1%; P = .016). Large endografts were associated with a lower postoperative aortic TAWSS (1.45 ± 0.76 dynes/cm² vs 3.16 ± 1.24 dynes/cm²; P < .001) and longer aortic RT (0.78 ± 0.30 second vs 0.34 ± 0.08 second; P < .001). In the large endograft group, a reduction >0.39 dynes/cm² in aortic TAWSS demonstrated discriminatory power for thrombotic complications (area under the receiver operating characteristic curve, 0.77). An increased aortic RT of ≥0.05 second had similar accuracy for predicting thrombotic complications (area under the receiver operating characteristic curve, 0.78). The odds of thrombotic complications were significantly higher if patients had met the hemodynamic threshold changes in aortic TAWSS (odds ratio, 7.0; 95% confidence interval, 1.1-45.9) and RT (odds ratio, 8.0; 95% confidence interval, 1.13-56.8).

Conclusions

Patient-specific CFD simulation of fEVAR in juxtarenal aortic aneurysms demonstrated significant endograft diameter-dependent differences in aortic hemodynamics. A postoperative reduction in TAWSS and an increased RT correlated with future thrombotic events after large-diameter endograft implantation. Patient-specific simulation of hemodynamics provides a novel method for thrombotic risk stratification after fEVAR.

The durability of fenestrated endovascular aneurysm repair (fEVAR) has been threatened by thrombotic complications. Using patient-specific computational flow simulation, the present retrospective study of 36 patients with juxtarenal aortic aneurysms treated with fEVAR identified several endograft diameter-dependent changes in aortic hemodynamics associated with thrombotic complications. A postoperative reduction in aortic wall shear stress and increased particle residence time correlated with the development of intraluminal graft thrombus and renal stent occlusion in patients treated with large diameter (>34 mm) endografts. These computationally estimated hemodynamic parameters could provide a novel method for patient-specific risk stratification for adverse events after fEVAR.

Simulation of Mechanical Heart Valve Dysfunction and the Non-Newtonian Blood Model Approach

The mechanical heart valve (MHV) is commonly used for the treatment of cardiovascular diseases. Nonphysiological hemodynamic in the MHV may cause hemolysis, platelet activation, and an increased risk of thromboembolism. Thromboembolism may cause severe complications and valve dysfunction. This paper thoroughly reviewed the simulation of physical quantities (velocity distribution, vortex formation, and shear stress) in healthy and dysfunctional MHV and reviewed the non-Newtonian blood flow characteristics in MHV. In the MHV numerical study, the dysfunction will affect the simulation results, increase the pressure gradient and shear stress, and change the blood flow patterns, increasing the risks of hemolysis and platelet activation. The blood flow passes downstream and has obvious recirculation and stagnation region with the increased dysfunction severity. Due to the complex structure of the MHV, the non-Newtonian shear-thinning viscosity blood characteristics become apparent in MHV simulations. The comparative study between Newtonian and non-Newtonian always shows the difference. The shear-thinning blood viscosity model is the basics to build the blood, also the blood exhibiting viscoelastic properties. More details are needed to establish a complete and more realistic simulation.

The influence of non-conformal grid interfaces on the results of large eddy simulation of centrifugal blood pumps

BACKGROUND

Computational fluid dynamics has been widely used to assist the design and evaluation of blood pumps. Discretization errors associated with computational grid may influence the credibility of numerical simulations. Non-conformal grid interfaces commonly exist in rotary machines, including rotary blood pumps. Should grid size across the interface differ greatly, large errors may occur.

METHODS

This study explored the effects of non-conformal grid interface on the prediction of the flow field and hemolysis in blood pumps using large eddy simulation (LES). Two benchmarks, a nozzle model and a centrifugal blood pump were chosen as test cases.

RESULTS

This study found that non-conformal grid interfaces with considerable change of grid sizes led to discontinuities of flow variables and brought errors to metrics such as pressure head (7%) and hemolysis (up to 14%).

CONCLUSIONS

The results on the full unstructured grid are more accurate with negligible changes of flow variables across the non-conformal grid interface. A full unstructured grid should be employed for centrifugal blood pumps to minimize the influence of non-conformal grid interfaces for LES simulations.

Hemolytic Footprint of Rotodynamic Blood Pumps

OBJECTIVE

In preclinical examinations, rotodynamic blood pumps (RBPs) are predominantly evaluated at design-point conditions. In clinical practice, however, they run at diversified modes of operation. This study aimed at extending current preclinical evaluation of hemolytic profiles in RBPs toward broader, clinically relevant ranges of operation.

METHODS

Two implantable RBPs the HeartMate 3 (HM3) and the HeartWare Ventricular Assist Device (HVAD) were analyzed at three pump speeds (HM3: 4300, 5600, 7000rpm; HVAD: 1800, 2760, 3600rpm) with three flow rates (1-9L/min) per speed setting. Hemolysis measurements were performed in heparinized bovine blood. The delta free hemoglobin (dfHb) and the normalized index of hemolysis (NIH) served as hemolytic measures. Statistical analysis was performed by multiple comparison of the 9 operating conditions. Moreover, computational fluid dynamics (CFD) was applied to provide mechanistic insights into the interrelation between hydraulics and hemolysis by correlating numerically computed hydraulic losses with in-vitro hemolytic measures.

RESULTS

In both devices, dfHb increased toward increasing speeds, particularly during low but also during high flow condition. By contrast, in both RBPs magnitudes of NIH were significantly elevated during low flow operation compared to high flow conditions (p<0.0036). Maps of hemolytic metrics revealed morphologically similar trends to in-silico hydraulic losses (r>0.793).

CONCLUSIONS

While off-design operation is associated with increased hemolytic profiles, the setting of different operating conditions render a preclinical prediction of clinical impact with current hemolysis metrics difficult.

SIGNIFICANCE

The identified increase in hemolytic measures during episodes of off-design operation is highlighting the need to consider worst-case operation during preclinical examinations.

An Accelerated Thrombosis Model for Computational Fluid Dynamics Simulations in Rotary Blood Pumps

Purpose

Thrombosis ranks among the major complications in blood-carrying medical devices and a better understanding to influence the design related contribution to thrombosis is desirable. Over the past years many computational models of thrombosis have been developed. However, numerically cheap models able to predict localized thrombus risk in complex geometries are still lacking. The aim of the study was to develop and test a computationally efficient model for thrombus risk prediction in rotary blood pumps.

Methods

We used a two-stage approach to calculate thrombus risk. The first stage involves the computation of velocity and pressure fields by computational fluid dynamic simulations. At the second stage, platelet activation by mechanical and chemical stimuli was determined through species transport with an Eulerian approach. The model was compared with existing clinical data on thrombus deposition within the HeartMate II. Furthermore, an operating point and model parameter sensitivity analysis was performed.

Results

Our model shows good correlation (R² > 0.93) with clinical data and identifies the bearing and outlet stator region of the HeartMate II as the location most prone to thrombus formation. The calculation of thrombus risk requires an additional 10–20 core hours of computation time.

Conclusion

The concentration of activated platelets can be used as a surrogate and computationally low-cost marker to determine potential risk regions of thrombus deposition in a blood pump. Relative comparisons of thrombus risk are possible even considering the intrinsic uncertainty in model parameters and operating conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00606-y.

Distribution and history of extensional stresses on vWF surrogate molecules in turbulent flow

The configuration of proteins is critical for their biochemical behavior. Mechanical stresses that act on them can affect their behavior leading to the development of decease. The von Willebrand factor (vWF) protein circulating with the blood loses its efficacy when it undergoes non-physiological hemodynamic stresses. While often overlooked, extensional stresses can affect the structure of vWF at much lower stress levels than shear stresses. The statistical distribution of extensional stress as it applies on models of the vWF molecule within turbulent flow was examined here. The stress on the molecules of the protein was calculated with computations that utilized a Lagrangian approach for the determination of the molecule trajectories in the flow filed. The history of the stresses on the proteins was also calculated. Two different flow fields were considered as models of typical flows in cardiovascular mechanical devises, one was a Poiseuille flow and the other was a Poiseuille–Couette flow field. The data showed that the distribution of stresses is important for the design of blood flow devices because the average stress can be below the critical value for protein damage, but tails of the distribution can be outside the critical stress regime.

Experimental Hemodynamics within the Penn State Fontan Circulatory Assist Device

For children born with a single functional ventricle, the Fontan operation bypasses the right ventricle by forming a four-way total cavopulmonary connection adapting the existing ventricle for the systemic circulation. However, upon adulthood, many Fontan patients exhibit low cardiac output and elevated venous pressure, eventually requiring a heart transplantation. Despite efforts to develop a Fontan pump or use an existing ventricular assist device for failing Fontan support, there is still no device designed or tested for subpulmonary support. Penn State University is developing a hydrodynamically levitated Fontan circulatory assist device (FCAD) for bridge-to-transplant or destination therapy. The FCAD hemodynamics, at both steady and pulsatile conditions for three pump operating conditions, were quantified using particle image velocimetry to determine the velocity magnitudes and Reynolds normal and shear stresses. Data were acquired at three planes (0 mm and ±25% of the radius) for the inferior and superior vena cavae inlets and the pulmonary artery outlet. The inlets had a blunt velocity profile that became skewed towards the collecting volute as fluid approached the rotor. At the outlet, regardless of the flow condition, a high-velocity jet exited the volute and moved downstream in a helical pattern. Turbulent stresses observed at the volute exit were influenced by the rotor's rotation. Regardless of inlet conditions, the pump demonstrated advantageous behavior for clinical use with a predictable flow field and a low risk of platelet adhesion and hemolysis based on calculated wall shear rates and turbulent stresses, respectively.

The effect of turbulence modelling on the assessment of platelet activation

Pathological platelet activation by abnormal shear stresses is regarded as a main clinical complication in recipients of cardiovascular mechanical devices. In order to improve their performance computational fluid dynamics (CFD) are used to evaluate flow fields and related shear stresses. CFD models are coupled with mathematical models that describe the relation between fluid dynamics variables, and in particular shear stresses, and the platelet activation state (PAS). These models typically use a Lagrangian approach to compute the shear stresses along possible platelet trajectories. However, in the case of turbulent flow, the choice of the proper turbulence closure is still debated for both concerning its effect on shear stress calculation and Lagrangian statistics. In this study different numerical simulations of the flow through a mechanical heart valve were performed and then compared in terms of Eulerian and Lagrangian quantities: a direct numerical simulation (DNS), a large eddy simulation (LES), two Reynolds-averaged Navier-Stokes (RANS) simulations (SST k-ω and RSM) and a "laminar" (no turbulence modelling) simulation. Results exhibit a large variability in the PAS assessment depending on the turbulence model adopted. "Laminar" and RSM estimates of platelet activation are about 60% below DNS, while LES is 16% less. Surprisingly, PAS estimated from the SST k- ω velocity field is only 8% less than from DNS data. This appears more artificial than physical as can be inferred after comparing frequency distributions of PAS and of the different Lagrangian variables of the mechano-biological model of platelet activation. Our study indicates how much turbulence closures may affect platelet activation estimates, in comparison to an accurate DNS, when assessing blood damage in blood contacting devices.

The effects of non-Newtonian blood modeling and pulsatility on hemodynamics in the food and drug administration's benchmark nozzle model

BACKGROUND

Computational fluid dynamics (CFD) is an important tool for predicting cardiovascular device performance. The FDA developed a benchmark nozzle model in which experimental and CFD data were compared, however, the studies were limited by steady flows and Newtonian models.

OBJECTIVE

Newtonian and non-Newtonian blood models will be compared under steady and pulsatile flows to evaluate their influence on hemodynamics in the FDA nozzle.

METHODS

CFD simulations were validated against the FDA data for steady flow with a Newtonian model. Further simulations were performed using Newtonian and non-Newtonian models under both steady and pulsatile flows.

RESULTS

CFD results were within the experimental standard deviations at nearly all locations and Reynolds numbers. The model differences were most evident at Re = 500, in the recirculation regions, and during diastole. The non-Newtonian model predicted blunter upstream velocity profiles, higher velocities in the throat, and differences in the recirculation flow patterns. The non-Newtonian model also predicted a greater pressure drop at Re = 500 with minimal differences observed at higher Reynolds numbers.

CONCLUSIONS

An improved modeling framework and validation procedure were used to further investigate hemodynamics in geometries relevant to cardiovascular devices and found that accounting for blood's non-Newtonian and pulsatile behavior can lead to large differences in predictions in hemodynamic parameters.

Validation of Numerically Predicted Shear Stress-dependent Dissipative Losses Within a Rotary Blood Pump

Computational fluid dynamics find widespread application in the development of rotary blood pumps (RBPs). Yet, corresponding simulations rely on shear stress computations that are afflicted with limited resolution while lacking validation. This study aimed at the experimental validation of integral hydraulic properties to analyze global shear stress resolution across the operational range of a novel RBP. Pressure head and impeller torque were numerically predicted based on Unsteady Reynolds-averaged Navier-Stokes (URANS) simulations and validated on a testbench with integrated sensor modalities (flow, pressure, and torque). Validation was performed by linear regression and Bland-Altman analysis across nine operating conditions. In power loss analysis (PLA), in silico hydraulic power losses were derived based on the validated hydraulic quantities and balanced with in silico shear-dependent dissipative power losses. Discrepancies among both terms provided a measure of in silico shear stress resolution. In silico and in vitro data correlated with low discordance in pressure (r = 0.992, RMSE = 1.02 mmHg), torque (r = 0.999, RMSE = 0.034 mNm), and hydraulic power losses (r = 0.990, RMSE = 0.015W). PLA revealed numerically predicted dissipative losses to be up to 34.4% smaller than validated computations of hydraulic losses. This study confirmed the suitability of URANS settings to predict integral hydraulic properties. However, numerical credibility was hampered by lacking resolution of shear-dependent dissipative losses.

Large eddy simulation as a fast and accurate engineering approach for the simulation of rotary blood pumps

An accurate representation of the flow field in blood pumps is important for the design and optimization of blood pumps. The primary turbulence modeling methods applied to blood pumps have been the Reynolds-averaged Navier-Stokes (RANS) or URANS (unsteady RANS) method. Large eddy simulation (LES) method has been introduced to simulate blood pumps. Nonetheless, LES has not been widely used to assist in the design and optimization of blood pumps to date due to its formidable computational cost. The purpose of this study is to explore the potential of the LES technique as a fast and accurate engineering approach for the simulation of rotary blood pumps. The performance of "Light LES" (using the same time and spatial resolutions as the URANS) and LES in two rotary blood pumps was evaluated by comparing the results with the URANS and extensive experimental results. This study showed that the results of both "Light LES" and LES are superior to URANS, in terms of both performance curves and key flow features. URANS could not predict the flow separation and recirculation in diffusers for both pumps. In contrast, LES is superior to URANS in capturing these flows, performing well for both design and off-design conditions. The differences between the "Light LES" and LES results were relatively small. This study shows that with less computational cost than URANS, "Light LES" can be considered as a cost-effective engineering approach to assist in the design and optimization of rotary blood pumps.

Computational evaluation of inferior vena cava filters through computational fluid dynamics methods

Numerical simulation is growing in its importance toward the design, testing and evaluation of medical devices. Computational fluid dynamics and finite element analysis allow improved calculation of stress, heat transfer, and flow to better understand the medical device environment. Current research focuses not only on improving medical devices, but also on improving the computational tools themselves. As methods and computer technology allow for faster simulation times, iterations and trials can be performed faster to collect more data. Given the adverse events associated with long-term inferior vena cava (IVC) filter placement, IVC filter design and device evaluation are of paramount importance. This work reviews computational methods used to develop, test, and improve IVC filters to ultimately serve the needs of the patient.

Human upper-airway respiratory airflow: In vivo comparison of computational fluid dynamics simulations and hyperpolarized 129Xe phase contrast MRI velocimetry

Computational fluid dynamics (CFD) simulations of respiratory airflow have the potential to change the clinical assessment of regional airway function in health and disease, in pulmonary medicine and otolaryngology. For example, in diseases where multiple sites of airway obstruction occur, such as obstructive sleep apnea (OSA), CFD simulations can identify which sites of obstruction contribute most to airway resistance and may therefore be candidate sites for airway surgery. The main barrier to clinical uptake of respiratory CFD to date has been the difficulty in validating CFD results against a clinical gold standard. Invasive instrumentation of the upper airway to measure respiratory airflow velocity or pressure can disrupt the airflow and alter the subject's natural breathing patterns. Therefore, in this study, we instead propose phase contrast (PC) velocimetry magnetic resonance imaging (MRI) of inhaled hyperpolarized 129Xe gas as a non-invasive reference to which airflow velocities calculated via CFD can be compared. To that end, we performed subject-specific CFD simulations in airway models derived from 1H MRI, and using respiratory flowrate measurements acquired synchronously with MRI. Airflow velocity vectors calculated by CFD simulations were then qualitatively and quantitatively compared to velocity maps derived from PC velocimetry MRI of inhaled hyperpolarized 129Xe gas. The results show both techniques produce similar spatial distributions of high velocity regions in the anterior-posterior and foot-head directions, indicating good qualitative agreement. Statistically significant correlations and low Bland-Altman bias between the local velocity values produced by the two techniques indicates quantitative agreement. This preliminary in vivo comparison of respiratory airway CFD and PC MRI of hyperpolarized 129Xe gas demonstrates the feasibility of PC MRI as a technique to validate respiratory CFD and forms the basis for further comprehensive validation studies. This study is therefore a first step in the pathway towards clinical adoption of respiratory CFD.

Stochastic simulation of the FDA centrifugal blood pump benchmark

In the present study, the effect of physical and operational uncertainties on the hydrodynamic and hemocompatibility characteristics of a centrifugal blood pump designed by the U.S. food and drug administration is investigated. Physical uncertainties include the randomness in the blood density and viscosity, while the operational uncertainties are composed of the pump rotational speed, mass flow rate, and turbulence intensity. The non-intrusive polynomial chaos expansion has been employed to conduct the uncertainty quantification analysis. Additionally, to assess each stochastic parameter's influence on the quantities of interest, the sensitivity analysis is utilized through the Sobol' indices. For numerical simulation of the pump's blood flow, the SST [Formula: see text] turbulence model and a power-law model of hemolysis were employed. The pump's velocity field is profoundly affected by the rotational speed in the bladed regions and the mass flow rate in other zones. Furthermore, the hemolysis index is dominantly sensitive to blood viscosity. According to the results, pump hydraulic characteristics (i.e., head and efficiency) show a more robust behavior than the hemocompatibility characteristics (i.e., hemolysis index) regarding the operational and physical uncertainties. Finally, it was found that the probability distribution function of the hemolysis index covers the experimental measurements.

Heuristic optimization of impeller sidewall gaps-based on the bees algorithm for a centrifugal blood pump by CFD

Optimization studies on blood pumps that require complex designs are gradually increasing in number. The essential design criteria of centrifugal blood pump are minimum shear stress with maximal efficiency. The geometry design of impeller sidewall gaps (blade tip clearance, axial gap, radial gap) is highly effective with regard to these two criteria. Therefore, unlike methods such as trial and error, the optimal dimensions of these gaps should be adjusted via a heuristic method, giving more effective results. In this study, the optimal gaps that can ensure these two design criteria with The Bees Algorithm (BA), which is a population-based heuristic method, are investigated. Firstly, a Computational Fluid Dynamics (CFD) analysis of sample pump models, which are selected according to the orthogonal array and pre-designed with different gaps, are performed. The dimensions of the gaps are optimized through this mathematical model. The simulation results for the improved pump model are nearly identical to those predicted by the BA. The improved pump model, as designed with the optimal gap dimensions so obtained, is able to meet the design criteria better than all existing sample pumps. Thanks to the optimal gap dimensions, it has been observed that compared to average values, it has provided a 42% reduction in aWSS and a 20% increase in efficiency. Moreover, original an approach to the design of impeller sidewall gaps was developed. The results show that computational costs have been significantly reduced by using the BA in blood pump geometry design.

Hemolysis estimation in turbulent flow for the FDA critical path initiative centrifugal blood pump

Hemolysis in medical devices and implants has been a primary concern over the past fifty years. Turbulent flow in particular can cause cell trauma and hemolysis in such devices. In this work, the effects of turbulence on red blood cell (RBC) damage are examined by simulating the flow field through a centrifugal blood pump that has been identified as a case study through the critical path initiative of the US Food and Drug Administration (FDA). In this study, a new model was employed to predict hemolysis in the turbulent flow environment in the pump selected for the FDA critical path initiative. The operating conditions for a centrifugal blood pump were specified by the FDA for rotational speeds of 2500 and 3500 rpm. The model is based on the analysis of the smaller eddies within the turbulent flow field, since it is assumed that turbulent flow eddies with sizes comparable to the dimensions of the RBCs lead to cell trauma. The Kolmogorov length scale of the velocity field is used to identify such small eddies. Using model parameters obtained in prior work through comparisons to capillary and jet flow, it is found that hemolysis for the 2500-rpm pump was predicted well, while hemolysis for the 3500-rpm pump was overpredicted. Results indicate refinement of the model and empirical constants with better experimental data is needed.

Dynamics of Blood Flows in Aortic Stenosis: Mild, Moderate, and Severe

Supraphysiologic high shear stresses created in calcific aortic stenosis (AS) are known to cause hemostatic abnormalities, however, the relationship between the complex blood flows over the severity of AS and hemostatic abnormalities still remains unclear. This study systematically characterized the blood flow in mild, moderate, and severe AS. A series of large eddy simulations (LES) validated by particle image velocimetry were performed on physiologically representative AS models with a peak physiologic flow condition of 18 liter per minute. Time-accurate velocity fields, transvalvular pressure gradient, and laminar viscous-and turbulent (or Reynolds) shear stresses (RSSmax) were evaluated for each degree of severity. The peak velocities of mild, moderate, and severe AS were on the order of 2.0, 4.0, and 8.0 m/s, respectively. Jet velocity in severe AS was highly skewed with extremely high velocity (as high as 8 m/s) and mainly traveled through the posterior aortic wall up to the aortic arch while still carrying a relatively high velocity, that is, >4 m/s. The mean laminar viscous wall shear stresses (WSS) for mild, moderate, and severe AS were on the order of 40, 100, and 180 Pa, respectively. The RSSmax were on the order of 260, 490, and 2,500 Pa for mild, moderate, and severe AS, respectively. This study may provide a link between altered flows in AS and hemostatic abnormalities such as acquired von Willebrand syndrome and hemolysis, thus, help diagnosing and timing of the treatment.

Validated Guidelines for Simulating Centrifugal Blood Pumps

PURPOSE

Rotary blood pumps (RBPs) employed as ventricular assist devices are developed to support the ventricles of patients suffering from heart failure. Computational Fluid Dynamics (CFD) is frequently used to predict the performance and haemocompatibility of these pumps during development, however different simulation techniques employed by various research groups result in inconsistent predictions. This inconsistency is further compounded by the lack of standardised model validation, thus it is difficult to determine which simulation techniques are accurate. To address these problems, the US Food and Drug Administration (FDA) proposed a simplified centrifugal RBP benchmark model. The aim of this paper was to determine simulation settings capable of producing accurate predictions using the published FDA results for validation.

METHODS

This paper considers several studies to investigate the impact of simulation options on the prediction of pressure and flow velocities. These included evaluation of the mesh density and interface position through steady simulations as well as time step size and turbulence models (k-ε realizable, k-ω SST, k-ω SST Intermittency, RSM ω-based, SAS and SBES) using a sliding mesh approach.

RESULTS

The most accurate steady simulation using the k-ω turbulence model predicted the pressure to within 5% of experimental results, however experienced issues with unphysical velocity fields. A more computationally expensive transient simulation that used the Stress-Blended Eddy Simulation (SBES) turbulence model provided a more accurate prediction of the velocity field and pressure rise to within experimental variation.

CONCLUSION

The findings of the study strongly suggest that SBES can be used to better predict RBP performance in the early development phase.

Patient-specific computational flow modelling for assessing hemodynamic changes following fenestrated endovascular aneurysm repair

Objective

This study aimed to develop an accessible patient-specific computational flow modelling pipeline for evaluating the hemodynamic performance of fenestrated endovascular aneurysm repair (fEVAR), with the hypothesis that computational flow modelling can detect aortic branch hemodynamic changes associated with fEVAR graft implantation.

Methods

Patients who underwent fEVAR for juxtarenal aortic aneurysms with the Cook ZFEN were retrospectively selected. Using open-source SimVascular software, preoperative and postoperative visceral aortic anatomy was manually segmented from computed tomography angiograms. Three-dimensional geometric models were then discretized into tetrahedral finite element meshes. Patient-specific pulsatile in-flow conditions were derived from known supraceliac aortic flow waveforms and adjusted for patient body surface area, average resting heart rate, and blood pressure. Outlet boundary conditions consisted of three-element Windkessel models approximated from physiologic flow splits. Rigid wall flow simulations were then performed on preoperative and postoperative models with the same inflow and outflow conditions. We used SimVascular's incompressible Navier-Stokes solver to perform blood flow simulations on a cluster using 72 cores.

Results

Preoperative and postoperative flow simulations were performed for 10 patients undergoing fEVAR with a total of 30 target vessels (20 renal stents, 10 mesenteric scallops). Postoperative models required a higher mean number of mesh elements to reach mesh convergence (3.2 ± 1.8 × 10⁶ vs 2.6 ± 1.1 × 10⁶; P = .005) with a longer mean computational time (10.3 ± 6.3 hours vs 7.8 ± 3.5 hours; P = .04) compared with preoperative models. fEVAR was associated with small but statistically significant increases in mean peak proximal aortic arterial pressure (140.3 ± 11.0 mm Hg vs 136.9 ± 8.7 mm Hg; P = .02) and peak renal artery pressure (131.6 ± 14.8 mm Hg vs 128.9 ± 11.8 mm Hg; P = .04) compared with preoperative simulations. No differences were observed in peak pressure in the celiac, superior mesenteric, or distal aortic arteries (P = .17-.96). When measuring blood flow, the only observed difference was an increase in peak renal flow rate after fEVAR (17.5 ± 3.8 mL/s vs 16.9 ± 3.5 mL/s; P = .04). fEVAR was not associated with changes in the mean pressure or the mean flow rate in the celiac, superior mesenteric, or renal arteries (P = .06-.98). Stenting of the renal arteries did not induce significant changes time-averaged wall shear stress in the proximal renal artery (23.4 ± 8.1 dynes/cm² vs 23.2 ± 8.4 dynes/cm²; P = .98) or distal renal artery (32.7 ± 13.9 dynes/cm² vs 29.6 ± 11.8 dynes/cm²; P = .23). In addition, computational visualization of cross-sectional velocity profiles revealed low flow disturbances associated with protrusion of renal graft fabric into the aortic lumen.

Conclusions

In a pilot study involving a selective cohort of patients who underwent uncomplicated fEVAR, patient-specific flow modelling was a feasible method for assessing the hemodynamic performance of various two-vessel fenestrated device configurations and revealed subtle differences in computationally derived peak branch pressure and blood flow rates. Structural changes in aortic flow geometry after fEVAR do not seem to affect computationally estimated renovisceral branch perfusion or wall shear stress adversely. Additional studies with invasive angiography or phase contrast magnetic resonance imaging are required to clinically validate these findings. (JVS–Vascular Science 2021;2:53-69.)

Clinical Relevance

Using a computational flow modelling for assessing the hemodynamic performance of fenestrated endovascular aneurysm repair (fEVAR), this real-world, patient-specific study included 10 participants and found that structural changes in aortic flow geometry after fEVAR did not seem to adversely impact estimated renal or visceral branch perfusion metrics (eg, peak and mean arterial pressure and flow rates) or wall shear stress. These findings overall support the ongoing clinical use of commercially available fEVAR devices for repair of juxtarenal aortic aneurysms, and provides a computational framework for future evaluation of fEVAR configurations in a preoperative or postoperative settings.

Application of physics-based flow models in cardiovascular medicine: Current practices and challenges

Personalized physics-based flow models are becoming increasingly important in cardiovascular medicine. They are a powerful complement to traditional methods of clinical decision-making and offer a wealth of physiological information beyond conventional anatomic viewing using medical imaging data. These models have been used to identify key hemodynamic biomarkers, such as pressure gradient and wall shear stress, which are associated with determining the functional severity of cardiovascular diseases. Importantly, simulation-driven diagnostics can help researchers understand the complex interplay between geometric and fluid dynamic parameters, which can ultimately improve patient outcomes and treatment planning. The possibility to compute and predict diagnostic variables and hemodynamics biomarkers can therefore play a pivotal role in reducing adverse treatment outcomes and accelerate development of novel strategies for cardiovascular disease management.