Piotr Reorowicz

Numerical simulations of the blood flow in the patient-specific arterial cerebral circle region

The Cerebral Circle Region, also known as the Circle of Willis (CoW), is a loop of arteries that form arterial connections between supply arteries to distribute blood throughout the cerebral mass. Among the population, only 25% to 50% have a complete system of arteries forming the CoW. 3D time-varying simulations for three different patient-specific artery anatomies of CoW were performed in order to gain a better insight into the phenomena existing in the cerebral blood flow. The models reconstructed on the basis of computer tomography images start from the aorta and include the largest arteries that supply the CoW and the arteries of CoW. Velocity values measured during the ultrasound examination have been compared with the results of simulations. It is shown that the flow in the right anterior artery in some cases may be supplied from the left internal carotid artery via the anterior communicating artery. The investigations conducted show that the computational fluid dynamic tool, which provides high resolution in both time and space domains, can be used to support physicians in diagnosing patients of different ages and various anatomical arterial structures.

Particle image velocimetry tests on pediatric 45-cc and 30-cc ventricle assist devices: effects of heart rate on VAD operation

Background This study investigated flow analysis inside pediatric ventricle assist devices (VADs) designed and manufactured at the Foundation for Cardiac Surgery Development (FRK), Zabrze, Poland. The main goal of the experiment was to define the minimal heart rate admissible in clinical practice. Methods The flow was directed by mechanical, single-disc valves developed at the Lodz University of Technology, Institute of Turbomachinery in Lodz, Poland. VAD operation conditions under different heart rates were analyzed. Measurements were performed on Religa PED pediatric VADs (45 cm3 and 30 cm3) with a particle image velocimetry (PIV) system. Results Due to the PIV method used, the measurements were made without interference of the measuring system onto the flow structure in the investigated channel, as the measurement procedure is noninvasive. During the investigations conducted in different measurement planes, the majority of the flow volume in the chamber was observable. Conclusions The measurements at different heart rates demonstrated a significant influence of this parameter on the flow nature in the heart ventricle. Additionally, it was found that the heart rate affected the operation of heart valves in the VAD.

Numerical investigations of the unsteady blood flow in the end-to-side arteriovenous fistula for hemodialysis

PURPOSE The aim of this study was to investigate the blood flow in the end-to-side arteriovenous (a-v) fistula, taking into account its pulsating nature and the patient-specific geometry of blood vessels. Computational Fluid Dynamics (CFD) methods were used for this analysis. METHODS DICOM images of the fistula, obtained from the angio-computed tomography, were a source of the data applied to develop a 3D geometrical model of the fistula. The model was meshed, then the ANSYS CFX v. 15.0 code was used to perform simulations of the flow in the vessels under analysis. Mesh independence tests were conducted. The non-Newtonian rheological model of blood and the Shear Stress Transport model of turbulence were employed. Blood vessel walls were assumed to be rigid. RESULTS Flow patterns, velocity fields, the volume flow rate, the wall shear stress (WSS) propagation on particular blood vessel walls were shown versus time. The maximal value of the blood velocity was identified in the anastomosis - the place where the artery is connected to the vein. The flow rate was calculated for all veins receiving blood. CONCLUSIONS The blood flow in the geometrically complicated a-v fistula was simulated. The values and oscillations of the WSS are the largest in the anastomosis, much lower in the artery and the lowest in the cephalic vein. A strong influence of the mesh on the results concerning the maximal and area-averaged WSS was shown. The relation between simulations of the pulsating and stationary flow under time-averaged flow conditions was presented.

Methods for determination of stagnation in pneumatic ventricular assist devices

Background: A pneumatic paediatric ventricular assist device developed at the Foundation of Cardiac Surgery Development, Zabrze, equipped with valves based on J. Moll’s design, with later modifications introduced at the Institute of Turbomachinery, Lodz University of Technology, was tested numerically and experimentally. The main aim of those investigations was to detect stagnation zones within the ventricular assist device and indicate advantages and limitations of both approaches. Methods: In the numerical transient test, a motion of the diaphragm and discs was simulated. Two different methods were used to illustrate stagnation zones in the ventricular assist device. The flow pattern inside the chamber was represented by velocity contours and vectors to validate the results using images obtained in the laser particle image velocimetry experiment. Results: The experimental light-based method implied problems with proper illumination of regions in the wall vicinity. High-resolution flow data and other important parameters as stagnation regions or flow patterns in regions not accessible for light in the particle image velocimetry method are derived in the numerical solution. However, computations of a single case are much more time-consuming if compared to a single experiment conducted on a well-calibrated stand. Conclusion: The resulting main vortexes in the central part of the pump chamber and the velocity magnitudes are correlated in both methods, which are complementary and when used together offer better insight into the flow structure inside the ventricular assist device and enable a deeper analysis of the results.