Analysis of Transitional and Turbulent Flow Through the FDA Benchmark Nozzle Model Using Laser Doppler Velocimetry

The Accelerated Transcatheter Heart Valve Testing Environment: Loading, Motion, and Fluid Dynamics

Transcatheter aortic valve replacements (TAVRs) are an increasingly common treatment for aortic valve disease due to their minimally invasive delivery. As TAVR designs require thinner leaflets to facilitate catheter-based delivery, they experience greater leaflet operational stresses and potentially greater durability issues than conventional surgical valves. Yet, our understanding of TAVR durability remains largely unexplored. Currently, preclinical TAVR durability is evaluated within an ISO:5840 compliant accelerated wear tester (AWT) up to a required 200 × 106 cycles, corresponding to approximately five years in vivo. While AWTs use high cycle frequencies (10-20 Hz) to achieve realistic timeframes, the resulting valve loading behaviors and fluid dynamics are not representative of the in vivo environment and thus may not accurately predict failure mechanisms. Despite the importance of fatigue and failure predictions for replacement heart valves, surprisingly, little quantitative information exists on the dynamic AWT environment. To better understand this environment, we examined frequency and diastolic period effects for the first time using high-speed enface imaging and particle image velocimetry to quantify valve motion and flow, respectively, using a Durapulse™ AWT at frequencies of 10, 15, and 20 Hz. Regardless of operating condition, no waveform achieved a physiologically relevant transvalvular loading pressure, despite having an ISO compliant geometric orifice area (GOA) and waveform. General fluid mechanics were consistent with in vivo but the AWT geometry developed secondary flow structures, which could impact mechanical loading. Therefore, the nonphysiologic loading and variability induced by changes in operating condition must be carefully regulated to ensure physiologically relevant fatigue.

Computational Prediction of Thrombosis in Food and Drug Administration's Benchmark Nozzle

Thrombosis seriously threatens human cardiovascular health and the safe operation of medical devices. The Food and Drug Administration’s (FDA) benchmark nozzle model was designed to include the typical structure of medical devices. However, the thrombosis in the FDA nozzle has yet not been investigated. The objective of this study is to predict the thrombus formation process in the idealized medical device by coupling computational fluid dynamics and a macroscopic hemodynamic-based thrombus model. We developed the hemodynamic-based thrombus model by considering the effect of platelet consumption. The thrombus model was quantitatively validated by referring to the latest thrombosis experiment, which was performed in a backward-facing step with human blood flow. The same setup was applied in the FDA nozzle to simulate the thrombus formation process. The thrombus shaped like a ring was firstly observed in the FDA benchmark nozzle. Subsequently, the accuracy of the shear-stress transport turbulence model was confirmed in different turbulent flow conditions. Five scenarios with different Reynolds numbers were carried out. We found that turbulence could change the shape of centrosymmetric thrombus to axisymmetric and high Reynolds number blood flow would delay or even prevent thrombosis. Overall, the present study reports the thrombosis process in the FDA benchmark nozzle using the numerical simulation method, and the primary findings may shed light on the effect of turbulence on thrombosis.

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.

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.

Transcatheter Heart Valve Downstream Fluid Dynamics in an Accelerated Evaluation Environment

Transcatheter aortic valve replacements (TAVRs) provide minimally invasive delivery of bioprosthetic heart valves (BHVs) for the treatment of aortic valve disease. While surgical BHVs show efficacy for 8-10 years, long-term TAVR durability remains unknown. Pre-clinical testing evaluates BHV durability in an ISO:5840 compliant accelerated wear tester (AWT), yet, the design and development of AWTs and their accuracy in predicting in vivo performance, is unclear. As a result of limited knowledge on AWT environment and BHV loading, durability assessment of candidate valves remains fundamentally empirical. For the first time, high-speed particle image velocimetry quantified an ISO:5840 compliant downstream AWT velocity field, Reynolds stresses, and turbulence intensity. TAVR enface imaging quantified the orifice area and estimated the flow rate. When valve area and flow rate were at their maximum during peak systole (1.49 cm² and 16.05 L/min, respectively), central jet velocity, Reynolds normal and shear stress, and turbulence intensity grew to 0.50 m/s, 265.1 Pa, 124.6 Pa, and 37.3%, respectively. During diastole, unique AWT recirculation produced retrograde flow and the directional changes created eddies. These novel AWT findings demonstrated a substantially reduced valve fully loaded period and pressure not matching in vivo or in vitro studies, despite the comparable fluid environment and TAVR motion.

The effect of turbulence on transitional flow in the FDA's benchmark nozzle model using large-eddy simulation

The Food and Drug Administration's (FDA) benchmark nozzle model has been studied extensively both experimentally and computationally. Although considerable efforts have been made on validations of a variety of numerical models against available experimental data, the transitional flow cases are still not fully resolved, especially with regards to detailed comparison of predicted turbulence quantities with experimental measurements. This study aims to fill this gap by conducting large-eddy simulations (LES) of flow through the FDA's benchmark model, at a transitional Reynolds number of 2000. Numerical results are compared to previous interlaboratory experimental results, with an emphasis on turbulence characteristics. Our results show that the LES methodology can accurately capture laminar quantities throughout the model. In the pre-jet breakdown region, predicted turbulence quantities are generally larger than high resolution experimental data acquired with laser Doppler velocimetry. In the jet breakdown regions, where maximum Reynolds stresses occur, Reynolds shear stresses show excellent agreement. Differences of up to 4% and 20% are observed near the jet core in the axial and radial normal Reynolds stresses, respectively. Comparisons between viscous and Reynolds shear stresses show that peak viscous shear stresses occur in the nozzle throat reaching a value of 18 Pa in the boundary layer, whilst peak Reynolds shear stresses occur in the jet breakdown region reaching a maximum value of 87 Pa. Our results highlight the importance in considering both laminar and turbulent contributions towards shear stresses and that neglecting the turbulence effect can significantly underestimate the total shear force exerted on the fluid.

Evaluating medical device and material thrombosis under flow: current and emerging technologies

This review highlights the importance of flow in medical device thrombosis and explores current and emerging technologies to evaluate dynamic biomaterial Thrombosis in vitro.

The FDA nozzle benchmark: "In theory there is no difference between theory and practice, but in practice there is"

The utility of flow simulations relies on the robustness of computational fluid dynamics (CFD) solvers and reproducibility of results. The aim of this study was to validate the Oasis CFD solver against in vitro experimental measurements of jet breakdown location from the FDA nozzle benchmark at Reynolds number 3500, which is in the particularly challenging transitional regime. Simulations were performed on meshes consisting of 5, 10, 17, and 28 million (M) tetrahedra, with Δt  =  10^-5 seconds. The 5M and 10M simulation jets broke down in reasonable agreement with the experiments. However, the 17M and 28M simulation jets broke down further downstream. But which of our simulations are "correct"? From a theoretical point of view, they are all wrong because the jet should not break down in the absence of disturbances. The geometry is axisymmetric with no geometrical features that can generate angular velocities. A stable flow was supported by linear stability analysis. From a physical point of view, a finite amount of "noise" will always be present in experiments, which lowers transition point. To replicate noise numerically, we prescribed minor random angular velocities (approximately 0.31%), much smaller than the reported flow asymmetry (approximately 3%) and model accuracy (approximately 1%), at the inlet of the 17M simulation, which shifted the jet breakdown location closer to the measurements. Hence, the high-resolution simulations and "noise" experiment can potentially explain discrepancies in transition between sometimes "sterile" CFD and inherently noisy "ground truth" experiments. Thus, we have shown that numerical simulations can agree with experiments, but for the wrong reasons.

Steady Flow in a Patient-Averaged Inferior Vena Cava-Part I: Particle Image Velocimetry Measurements at Rest and Exercise Conditions

PURPOSE

Although many previous computational fluid dynamics (CFD) studies have investigated the hemodynamics in the inferior vena cava (IVC), few studies have compared computational predictions to experimental data, and only qualitative comparisons have been made. Herein, we provide particle image velocimetry (PIV) measurements of flow in a patient-averaged IVC geometry under idealized conditions typical of those used in the preclinical evaluation of IVC filters.

METHODS

Measurements are acquired under rest and exercise flow rate conditions in an optically transparent model fabricated using 3D printing. To ensure that boundary conditions are well-defined and to make follow-on CFD validation studies more convenient, fully-developed flow is provided at the inlets (i.e., the iliac veins) by extending them with straight rigid tubing longer than the estimated entrance lengths. Velocity measurements are then obtained at the downstream end of the tubing to confirm Poiseuille inflow boundary conditions.

RESULTS

Measurements in the infrarenal IVC reveal that flow profiles are blunter in the sagittal plane (minor axis) than in the coronal plane (major axis). Peak in-plane velocity magnitudes are 4.9 cm/s and 27 cm/s under the rest and exercise conditions, respectively. Flow profiles are less parabolic and exhibit more inflection points at the higher flow rate. Bimodal velocity peaks are also observed in the sagittal plane at the elevated flow condition.

CONCLUSIONS

The IVC geometry, boundary conditions, and infrarenal velocity measurements are provided for download on a free and publicly accessible repository at https://doi.org/10.6084/m9.figshare.7198703 . These data will facilitate future CFD validation studies of idealized, in vitro IVC hemodynamics and of similar laminar flows in vascular geometries.

Determination of Reynolds Shear Stress Level for Hemolysis

Reynolds shear stress (RSS) has served as a metric for the effect of turbulence on hemolysis. Forstrom (1969) and Sallam and Hwang (1984) determined the RSS threshold for hemolysis to be 50,000 and 4,000 dyne/cm, respectively, using a turbulent jet. Despite the order of magnitude discrepancy, the threshold by Sallam and Hwang has been frequently cited for hemolytic potential in blood pumps. We recreated a Sallam apparatus (SA) to resolve this discrepancy and provide additional data to be used in developing a more accurate hemolysis model. Hemolysis was measured over a large range of Reynolds numbers (Re) (Re = 1,000-80,000). Washed bovine red blood cells (RBCs) were injected into the free jet of phosphate buffered saline, and hemolysis was quantified using a percent hemolysis, Hp = h (100 - hematocrit [HCT])/Hb, where h (mg/dl) is free hemoglobin and Hb (mg/dl) is total hemoglobin. Reynolds shear stress was calculated using two-dimensional laser Doppler velocimetry. Reynolds shear stress of ≥30,000 dyne/cm corresponding to Re of ≥60,000 appeared to cause hemolysis (p < 0.05). This RSS is an order of magnitude greater than the RSS threshold that Sallam and Hwang suggested, and it is similar to Forstrom's RSS threshold. This study resolved a long-standing uncertainty regarding the critical values of RSS for hemolysis and may provide a foundation for a more accurate hemolysis model.

Computational Fluid Dynamics Assessment Associated with Transcatheter Heart Valve Prostheses: A Position Paper of the ISO Working Group

The governing international standard for the development of prosthetic heart valves is International Organization for Standardization (ISO) 5840. This standard requires the assessment of the thrombus potential of transcatheter heart valve substitutes using an integrated thrombus evaluation. Besides experimental flow field assessment and ex vivo flow testing, computational fluid dynamics is a critical component of this integrated approach. This position paper is intended to provide and discuss best practices for the setup of a computational model, numerical solving, post-processing, data evaluation and reporting, as it relates to transcatheter heart valve substitutes. This paper is not intended to be a review of current computational technology; instead, it represents the position of the ISO working group consisting of experts from academia and industry with regards to considerations for computational fluid dynamic assessment of transcatheter heart valve substitutes.

Experimental quantification of the fluid dynamics in blood-processing devices through 4D-flow imaging: A pilot study on a real oxygenator/heat-exchanger module

The performance of blood-processing devices largely depends on the associated fluid dynamics, which hence represents a key aspect in their design and optimization. To this aim, two approaches are currently adopted: computational fluid-dynamics, which yields highly resolved three-dimensional data but relies on simplifying assumptions, and in vitro experiments, which typically involve the direct video-acquisition of the flow field and provide 2D data only. We propose a novel method that exploits space- and time-resolved magnetic resonance imaging (4D-flow) to quantify the complex 3D flow field in blood-processing devices and to overcome these limitations. We tested our method on a real device that integrates an oxygenator and a heat exchanger. A dedicated mock loop was implemented, and novel 4D-flow sequences with sub-millimetric spatial resolution and region-dependent velocity encodings were defined. Automated in house software was developed to quantify the complex 3D flow field within the different regions of the device: region-dependent flow rates, pressure drops, paths of the working fluid and wall shear stresses were computed. Our analysis highlighted the effects of fine geometrical features of the device on the local fluid-dynamics, which would be unlikely observed by current in vitro approaches. Also, the effects of non-idealities on the flow field distribution were captured, thanks to the absence of the simplifying assumptions that typically characterize numerical models. To the best of our knowledge, our approach is the first of its kind and could be extended to the analysis of a broad range of clinically relevant devices.