Modelling Circulation in Blood Vessel Aneurysms

Introduction. A two-dimensional hydrodynamic problem in the «stream function–vorticity» variables is numerically solved in an open rectangular cavity simulating blood flow in a blood vessel aneurysm. Two solution algorithms are proposed for Reynolds numbers Re < 1 and for Re ≥ 1. Materials and Methods. To accelerate the numerical solution with an explicit finite-difference scheme for the vorticity dynamics equation, the initial condition damping method, the n-fold splitting method of the explicit finite-differences cheme ( n = 100, 200), and the symmetry plane of the rectangular cavity–aneurysm were employed. In the splitting method, the maximum time step proportional to the square of the spatial step was used without violating the spectral stability of the explicit scheme in the vorticity equation. On half of the rectangular aneurysm, symmetric solutions were considered with a uniform 100 × 50 grid and equal steps h 1 = h 2 = 0.01. The inverse matrix for solving the Poisson equation in the «stream function–vorticity» variables with a finite number of elementary operations was computed using the MSIMSL library. Results. The numerical solution showed that the number and location of circulation regions in the aneurysm at small Reynolds numbers depend on the ratio of the vessel diameter to the aneurysm diameter. At small values of this parameter, the aneurysm contains a single large vortex that narrows the vessel lumen in the case of thrombus formation inside the aneurysm. The narrowing of the blood flow tube inside the aneurysm reaches 34%. It was found that the formation of the hydrodynamic structure in the aneurysm occurs in a time negligible (0.002%) compared to the period between pulsation waves (1 s). For the first time, a boundary condition with fourth-order accuracy was proposed to relate velocity, vorticity, and stream function. Discussion. The approximation of the equations in systems (4) and (22) has sixth-order accuracy at interior nodes and fourth-order accuracy at boundary nodes. The problem was also solved for blood motion in arteries at high Reynolds numbers (Re = 1500). The solution shows that in the aneurysm symmetry plane a chain of connected vortices is formed with alternating signs of vorticity, carried by the blood flow along the vessel. Conclusion. The initial–boundary value problems (4), (22) formulated in this work make it possible to qualitatively model blood flow in aneurysms of capillaries, arterioles, and arteries at low and high velocities, as well as blood motion in elements of medical equipment.