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.

Multi-scale model of lumen formation via inverse membrane blebbing mechanism during sprouting angiogenesis process.

Cancer is one of the leading causes of mortality and morbidity among people worldwide. Cancer appears as solid tumors in many cases. Angiogenesis is the growth of blood vessels from the existing vasculature and is one of the imperative processes in tumor growth. Another vital phenomenon for formation and functionality of this vasculature network is lumen formation. The results of recent studies indicate the importance of blood pressure in this mechanism. Computational modeling can study these processes in different scales. Hence, wide varieties of these models have been proposed during recent years. In this research, a multi-scale model is developed for the angiogenesis process. In the extracellular scale, the growth factor concentration is calculated via the reaction diffusion equation. At the cellular scale, growth, migration, and the adhesion of endothelial cells are modeled by the Potts cellular model. At the intra-cellular scale by considering biochemical signals, a Boolean network model describes migration, division, or apoptosis of endothelial cells. A stochastic model developed for lumen formation via inverse membrane blebbing mechanism. A CFD simulation was also used to investigate the role of pulsated blood pressure in the inverse membrane blebbing mechanism. The lumen formation model shows stochastic behavior in blebs expansion and lumen expansion. Comparing the stochastic model's results with the CFD simulation also shows the vital role of pressure pulse and the topology of the blebs in bleb retraction.

Design of an aortic polymeric valve with asymmetric leaflets and evaluation of its performance by finite element method.

BACKGROUND: The geometry of leaflets plays a significant role in prosthetic valves' (PVs) performance. Typically, natural aortic valves have three unequal leaflets, which differ in size. The present study aims to design an asymmetric tri-leaflet polymeric valve with one large and two small leaflets based on commissure lengths and leaflet eccentricities. METHODS: Eccentricity was related to commissure lengths based on the deformation of the free margins for the fully-opened state of leaflets. The polystyrene-block-polyethylene-polypropylene-block-polystyrene polymer characterized the material properties of the leaflets. The Finite Element Method (FEM) was used to evaluate performance parameters, including maximum geometric orifice area (GOA), average GOA, maximum von Mises stress, and leaflet's coaptation surface area (CSA). RESULTS: Asymmetric valves with no eccentricity provided a low level of GOA because the asymmetric form of small leaflets caused them to close faster than the large leaflet, leading to a sudden drop in the GOA during systole. As the radial curve tends towards a straight line, an undesirable coaptation occurs, and peak stress increases despite higher GOAs. A new radial curve consisting of two straight lines connected by an arc that provided 25.64 mm2 coaptation surface area (CAS) and 117.54 mm2 average GOA, was proposed to improve coaptation and GOA. CONCLUSION: The radial curve of leaflets affects the valve's performance more than other geometric parameters. The combination of straight lines and arcs for radial curves was selected as the reference model for asymmetric valves with one large and two small leaflets.

Investigation the effect of geometry and position of polymeric heart valves on hemodynamic with fluid-structure interaction numerical method.

Computational modeling and numerical simulation of heart valve dynamics incorporating both fluid dynamics and valve structural replications has been challenging. In this study, we developed a double-coupled fluid-structure interaction (FSI) model using arbitrary lagrangian eulerian(ALE) and steered adaptive mesh(SAM). So we were looking to simulate transcatheter aortic valve (TAV) hemodynamic performance throughout entire cardiac cycles [1]. To reach this object, semi-real geometry of aorta and aortic polymeric valves has been created. At model inlet, left ventricular pressure and at the model outlet the elastic porous tube have been considered. Nonlinear finite element way and Sparse solver was utilized to couple fluid and solid equation. Consequently, extensive and comparative simulation were performed to investigate the impact of valve elasticity and valve positions on hemodynamics and solid parameters. Effective Orifice Area(EOA) also has been calculated [1]. The simulation results indicated that the lower of the elastic modulus cause to increase the EOA. Furthermore, the result of valve position showed, whenever the valve is closer to sinuses, a greater EOA and lower stresses impose on the leaflet are achievable.

Numerical study of hemodynamics in a complete coronary bypass with venous and arterial grafts and different degrees of stenosis.

Cardiovascular diseases are among the leading causes of death in the world. The coronary blockage is one of most common types of these diseases that in the majority of cases has been treated by bypass surgery. In the bypass surgery, a graft is implemented to alter the blocked coronary and allow the blood supply process. The hemodynamic characteristics of the bypass strongly depend on the geometry and mechanical properties of the graft. In the present study, the fluid-structure interaction (FSI) analysis is conducted to investigate the bypass performance for a thoracic artery as well as a saphenous vein graft. Blood flow introduces a pressure on the walls of the graft which behaves as a hyperelastic material. A complete coronary bypass with stenosis degrees of 70% and 100% is modeled. To consider the nonlinear stress-strain behavior of the grafts, a five parameter Mooney-Rivlin hyperplastic model is implemented for the structural analysis and blood is assumed to behave as a Newtonian fluid. The simulations are performed for a structured grid to solve the governing equations using finite element method (FEM). The results show that wall shear stress (WSS) for saphenous vein is larger than that of thoracic artery while the total deformation of the thoracic artery is larger compared to the saphenous vein. Also, for the venous grafts or lower stenosis degree, the oscillatory shear index (OSI) is higher at both left and right anastomoses meaning that venous grafts as well as lower degree of stenosis are more critical in terms of restenosis.

Local hemodynamic analysis of the C-Pulse Device by 3D fluid-structure interaction simulation.

Background: C-Pulse is a new, nonblood contacting device based on the concept of counter-pulsation that is designed for long-term implantation. However, there is a lack of comprehensive investigation of the pressure and velocity fields under the action of C-Pulse. Aim: In this paper, we aim to conduct a numerical simulation of the underlying mechanism of the device in order to analyze its performance and related undesirable issues. Materials & methods: A 3D finite element model is utilized to simulate the mechanism of the blood pumping. Results & conclusion: The simulation well reproduced the essential characteristics of the C-Pulse. Preliminary results were in a reasonable range while a couple of irregular flow patterns were identified.

Simulating glioblastoma growth consisting both visible and invisible parts of the tumor using a diffusion-reaction model followed by resection and radiotherapy.

Glioblastoma is known to be among one of the deadliest brain tumors in the world today. There have been major improvements in the detection of cancerous cells in the twenty-first century. However, the threshold of detection of these cancerous cells varies in different scanning techniques such as magnetic resonance imaging (MRI) and computed tomography (CT). The growth of these tumors and different treatments have been modeled to assist medical experts in better predictions of the related tumor growth and in the selection of more accurate treatments. In clinical terms the tumor consisted of two parts known as the visible part, which is the part of the tumor that is above the threshold of the detecting device and the invisible part, which is below the detecting threshold. In this study, the common reaction-diffusion model of tumor growth is used to simulate the growth of the glioblastoma tumor. Also resection and radiotherapy have been modeled as methods to prevent the growth of the tumor. The results demonstrate that although the selected treatments were effective in reducing the number of cancerous cells to under the threshold of detection, they did not eliminate all cancerous cells and if no further treatments were applied, the cancerous cells would spread and become malignant again. Although previous studies have suggested that the ratio of proliferation to diffusion could describe the malignancy of the tumor, this study in addition shows the importance of each of the coefficients regarding the malignancy of the tumor.

Personalized image-based tumor growth prediction in a convection-diffusion-reaction model.

Inter-individual heterogeneity of tumors leads to non-effectiveness of unique therapy plans. This issue has caused a growing interest in the field of personalized medicine and its application in tumor growth evaluation. Accordingly, in this paper, a framework of personalized medicine is presented for growth prediction of brain glioma tumors. A convection-diffusion-reaction model is used as the patient-specific tumor growth model which is associated with multimodal magnetic resonance images (MRIs). Two parameters of intracellular area fraction (ICAF) and metabolic rate have been used to incorporate the physiological data obtained from medical images into the model. The framework is tested on the data of two cases of glioma tumors to document the approach; parameter estimation is made using particle swarm optimization (PSO) and genetic algorithm (GA) and the model is evaluated by comparing the predicted tumors with the observed tumors in terms of root mean square error of the ICAF maps (IRMSE), relative area difference (RAD) and Dice's coefficient (DC). Results show the differences of IRMSE, RAD and DC in 4.1 ∓ 1.15%, 0.099 ∓ 0.041 and 85.5 ∓ 7.5%, respectively. Survival times are estimated by assuming the tumor radius of 35 mm as the fatal burden. Results confirm that less-diffusive tumors lead to higher survival times. The represented framework makes it possible to personally predict the growth behavior of glioma tumors only based on patients' routine MRIs and provides a basis for modeling the personalized therapy and walking in the path of personalized medicine.

Computational fluid dynamics-based study of possibility of generating pulsatile blood flow via a continuous-flow VAD.

Until recent years, it was almost beyond remedy to save the life of end-stage heart failure patients without considering a heart transplant. This is while the need for healthy organs has always far exceeded donations. However, the evolution of VAD technology has certainly changed the management of these patients. Today, blood pumps are designed either pulsatile flow or continuous flow, each of which has its own concerns and limitations. For instance, pulsatile pumps are mostly voluminous and hardly can be used for children. On the other hand, the flow generated by continuous-flow pumps is in contrast with pulsatile flow of the natural heart. In this project, having used computational fluid dynamics, we studied the possibility of generating pulsatile blood flow via a continuous-flow blood pump by adjusting the rotational speed of the pump with two distinct patterns (sinusoidal and trapezoidal), both of which have been proposed and set based on physiological needs and blood flow waveform of the natural heart. An important feature of this study is setting the outlet pressure of the pump similar to the physiological conditions of a patient with heart failure, and since these axial pumps are sensitive to outlet pressures, more secure and reliable results of their performance are achieved. Our results show a slight superiority of a sinusoidal pattern compared to a trapezoidal one with the potential to achieve an adequate pulsatile flow by precisely controlling the rotational speed.

Numerical simulation of osteocyte cell in response to directional mechanical loadings and mechanotransduction analysis: Considering lacunar-canalicular interstitial fluid flow.

The osteocyte cell is a bone cell that also functions as a bone mechanosensor. In this work, a three-dimensional (3D) fluid-structure interaction (FSI) model of an osteocyte cell under different mechanical loading conditions was used to obtain a better understanding of osteocyte cell behavior under different physiological conditions. In the current study, both fluid and solid parts of osteocyte cell were considered in order to allow for more accurate results. Five different loading conditions have been applied to the osteocyte cell, and consequently the different interstitial fluid flow velocities and shear stresses have been investigated. Furthermore, using a mathematical model, the change in the stimulus function value with shear stress and NO enzyme was revealed. This work suggests that changes in osteocyte morphology and direction of loadings affect cell stimulation. It was found that cell is mostly stimulated and expanded in the direction experiencing the most shear stress. Finally, the amount of cell stimulation was shown quantitatively and there was strong dependency between stimulus function, shear stress, calcium, and NO concentration.

A comparative study of thermal effects of 3 types of laser in eye: 3D simulation with bioheat equation.

The application of laser in ophthalmology and eye surgery is so widespread that hardly can anyone deny its importance. On the other hand, since the human eye is an organ susceptible to external factors such as heat waves, laser radiation rapidly increases the temperature of the eye and therefore the study of temperature distribution inside the eye under laser irradiation is crucial; but the use of experimental and invasive methods for measuring the temperature inside the eye is typically high-risk and hazardous. In this paper, using the three-dimensional finite element method, the distribution of heat transfer inside the eye under transient condition was studied through three different lasers named Nd:Yag, Nd:Yap and ArF. Considering the metabolic heat and blood perfusion rate in various regions of the eye, numerical solution of space-time dependant Pennes bioheat transfer equation has been applied in this study. Lambert-Beer's law has been used to model the absorption of laser energy inside the eye tissues. It should also be mentioned that the effect of the ambient temperature, tear evaporation rate, laser power and the pupil diameter on the temperature distribution have been studied. Also, temperature distribution inside the eye after applying each laser and temperature variations of six optional regions as functions of time have been investigated. The results show that these radiations cause temperature rise in various regions, which will in turn causes serious damages to the eye tissues. Investigating the temperature distribution inside the eye under the laser irradiation can be a useful tool to study and predict the thermal effects of laser radiation on the human eye and evaluate the risk involved in performing laser surgery.

A new design and computational fluid dynamics study of an implantable axial blood pump.

Considering small thoracic space, using implantable ventricular assist device requires reduction in a pump size. Among many available blood pumps, axial blood pumps have attracted greatly because of their small size. In this article, a new miniature axial blood pump has been designed and studied which can be easily implanted in the human body. In this design, the pump overall length decreased by a little increasing in the pump diameter, and new blade geometry is used to produce a streamlined, idealized, and nonobstructing blood flow path in the pump. By means of computational fluid dynamic, the flow pattern through the pump has been predicted and overall pump performance and efficiency has been computed. Also, to ensure a reliable VAD design, two methods for checking wall shear stress were used to confirm that this pump wouldn't cause serious blood damage.

Comparison of blood flow velocity through the internal carotid artery based on Doppler ultrasound and numerical simulation.

Doppler ultrasound is a usual non-invasive method to estimate the stenosis percentage in large arteries such as carotid by measuring maximum velocity of blood flow. Based on clinical investigations, because of vessel wall motions, Doppler positioning and angle correction, some errors can arise in Doppler results which lead to incorrect diagnosis. The aim of this study was to compare the results of Doppler test and the numerical simulation of blood flow in the same case. For this evaluation, two patients including an 87-year-old man and a 72-year-old woman suffering from stenosis in the internal carotid artery were selected. First, clinical information of each patient such as CT-Angio scan images and Doppler ultrasound results on different locations of the stenosed artery were obtained. Then, the geometries were reconstructed and numerical simulations were carried out using ANSYS software. Results showed that the velocity profile of Doppler test and numerical simulation were in good agreement at the regions of pre-and post-stenosis. However, the value of maximum velocity at the stenotic region had significant differences.