Transient, three-dimensional flow field simulation through a mechanical, trileaflet heart valve prosthesis.

Thromboembolic complications are one of the major challenges faced by designers and researchers in development of artificial organs with blood-contacting devices such as heart valve prostheses, especially mechanical valves. Besides increasing the thrombogenic potential, these valves change the hydrodynamic performance of the heart. In this study, the flow through a trileaflet, mechanical heart valve prosthesis was modeled with transient computational fluid dynamics to analyze flow patterns causing thrombus formations on valves. The valve was simulated under conditions of a test rig (THIA II), which was specially designed to analyze different valves with respect to thrombosis. The main goal of this study was to mimic the exact conditions of the test rig to be able to compare numerical and experimental results. The boundary conditions were obtained from experimental data as leaflet kinematics and pressure profiles. One complete cycle of the valve was simulated. Numerical flow and pressure results were analyzed and compared with experimental results. Shear stress and shear rates were determined with respect to thrombogenic potential, especially in the pivot regions, which seem to be the main influence for activation and deposition of thrombocytes. Approximately 0.7% of the blood volume moving through the fluid domain of the valve was exposed to shear rates high enough to cause platelet activation. However, shear rates of up to 20,000 s⁻¹ occurred in pivot regions. The pressure differences between the simulation and experimental data were approximately 2.5% during systole and increased up to 25% during diastole. The presented method, however, can be used to gain more information about the flow through different heart valve prostheses and, thus, improve the development process.

A validated CFD model to predict O2 and CO2 transfer within hollow fiber membrane oxygenators.

Hollow fiber oxygenators provide gas exchange to and from the blood during heart surgery or lung recovery. Minimal fiber surface area and optimal gas exchange rate may be achieved by optimization of hollow fiber shape and orientation (1). In this study, a modified CFD model is developed and validated with a specially developed micro membrane oxygenator (MicroMox). The MicroMox was designed in such a way that fiber arrangement and bundle geometry are highly reproducible and potential flow channeling is avoided, which is important for the validation. Its small size (V(Fluid)=0.04 mL) allows the simulation of the entire bundle of 120 fibers. A non-Newtonian blood model was used as simulation fluid. Physical solubility and chemical bond of O2 and CO2 in blood was represented by the numerical model. Constant oxygen partial pressure at the pores of the fibers and a steady state flow field was used to calculate the mass transport. In order to resolve the entire MicroMox fiber bundle, the mass transport was simulated for symmetric geometry sections in flow direction. In vitro validation was achieved by measurements of the gas transfer rates of the MicroMox. All measurements were performed according to DIN EN 12022 (2) using porcine blood. The numerical simulation of the mass transfer showed good agreement with the experimental data for different mass flows and constant inlet partial pressures. Good agreement could be achieved for two different fiber configurations. Thus, it was possible to establish a validated model for the prediction of gas exchange in hollow fiber oxygenators.

Flow distribution during cardiopulmonary bypass in dependency on the outflow cannula positioning.

Oxygen deficiency in the right brain is a common problem during cardiopulmonary bypass (CPB). This is linked to an insufficient perfusion of the carotid and vertebral artery. The flow to these vessels is strongly influenced by the outflow cannula position, which is traditionally located in the ascending aorta. Another approach however is to return blood via the right subclavian artery. A computational fluid dynamics (CFD) study was performed for both methods and validated by particle image velocimetry (PIV). A 3-dimensional computer aided design model of the cardiovascular (CV) system was generated from realtime computed tomography and magnetic resonance imaging data. Mesh generation (CFD) and rapid prototyping (PIV) were used for the further model creation. The simulations were performed assuming usual CPB conditions, and the same boundary conditions were applied for the PIV validation. The flow distribution was analyzed for 55 cannula positions inside the aorta and in relation to the distance between the cannula tip and the vertebral artery branch for subclavian cannulation. The study reveals that the Venturi effect due to the cannula jet appears to be the main reason for the loss in cerebral perfusion seen clinically. It provides a PIV-validated CFD method of analyzing the flow distribution in the CV system and can be transferred to other applications.

The impact of aortic/subclavian outflow cannulation for cardiopulmonary bypass and cardiac support: a computational fluid dynamics study.

Approximately 100 000 cases of oxygen deficiency in the brain occur during cardiopulmonary bypass (CPB) procedures each year. In particular, perfusion of the carotid and vertebral arteries is affected. The position of the outflow cannula influences the blood flow to the cardiovascular system and thus end organ perfusion. Traditionally, the cannula returns blood into the ascending aorta. But some surgeons prefer cannulation to the right subclavian artery. A computational fluid dynamics study was initially undertaken for both approaches. The vessel model was created from real computed tomography/magnetic resonance imaging data of young healthy patients. The simulations were run with usual CPB conditions. The flow distribution for different cannula positions in the aorta was studied, as well as the impact of the cannula tip distance to vertebral artery for the subclavian position. The study presents a fast method of analyzing the flow distribution in the cardiovascular system, and can be adapted for other applications such as ventricular assist device support. It revealed that two effects cause the loss of perfusion seen clinically: a vortex under the brachiocephalic trunk and low pressure regions near the cannula jet. The results suggest that cannulation to the subclavian artery is preferred if the cannula tip is sufficiently far away from the branch of the vertebral artery. For the aortic positions, however, the cannula should be injected from the left body side.

A validated computational fluid dynamics model to estimate hemolysis in a rotary blood pump.

A major part of developing rotary blood pumps requires the optimization of hemolytic properties of the entire pump. Application of a suited computational fluid dynamics (CFD)-based hemolysis model allows approximation of blood damage in an early phase of the design process. Thus, a drastic reduction of time- and cost- intensive hemolysis experiments can be achieved. For the MicroDiagonal Pump (MDP), still under development at Helmholtz-Institute in Aachen, Germany, different pump configurations have been analyzed, both numerically and experimentally. The CFD model of the pump has been successfully validated based on the comparison of the pressure head curves (H-Q curves), as discussed in a prior publication. In the present study, the authors focus on the development of a semiempiric blood damage model using the CFD and in vitro hemolysis data. On the one hand, mean key characteristics (shear stress and exposure time) and other characteristics affecting blood damage have been calculated based on numerical data. On the other hand, in vitro hemolysis tests have been accomplished in order to determine the hemolytic curves of two different pump configurations (with the same impeller but different tip clearances). Finally, a new function based on a general power law has been defined by means of the mean key characteristics. The unknown constants of the function have been determined by multidimensional regression analysis using the hemolytic curves. For the final validation of this new blood damage model, the calculated and the in vitro obtained hemolysis indices at the specific VAD operating point have been compared for all pump configurations. The comparison showed an excellent agreement, both qualitatively and quantitatively.

Comparison of hydraulic and hemolytic properties of different impeller designs of an implantable rotary blood pump by computational fluid dynamics.

A mixed-flow blood pump for long-term applications has been developed at the Helmholtz-Institute in Aachen, Germany. Central features of this implantable pump are a centrally integrated motor, a blood-immersed mechanical bearing, magnetic coupling of the impeller, and a shrouded impeller, which allows a relatively wide clearance. The aim of the study was a numerical analysis of hydraulic and hemolytic properties of different impeller design configurations. In vitro testing and numerical simulation techniques (computational fluid dynamics [CFD]) were applied to achieve a comprehensive overview. Pressure-flow charts were experimentally measured in a mock loop in order to validate the CFD data. In vitro hemolysis tests were performed at the main operating point of each impeller design. General flow patterns, pressure-flow charts, secondary flow rates, torque, and axial forces on the impeller were calculated by means of CFD. Furthermore, based on streak line techniques, shear stress (stress loading), exposure times, and volume percentage with critical stress loading have been determined. Comparison of CFD data with pressure head measurements showed excel-lent agreement. Also, impressive trend conformity was observed between in-vitro hemolysis results and numerical data. Comparison of design variations yielded clear trends and results. Design C revealed the best hydraulic and hemolytic properties and was chosen as the final design for the mixed-flow rotary blood pump.