CC BY
49
УДК 612.13.001.575
MATHEMATICAL MODELING OF BLOOD FLOW IN BASILAR ARTERY BIFURCATION REGION*
S. V. Frolov1, S. V. Sindeev1, D. Liepsch2, A. Balasso3
Department "Biomedical Engineering", TSTU (1); ssindeev@yandex.ru;
Department of Mechanical Engineering, Munich University of Applied Sciences, Munich (Germany) (2); Interdisciplinary Research Laboratory of the Klinikum Rechts der Isar, Technical University of Munich, Munich (Germany) (3)
Key words and phrases: aneurysm; basilar artery; cardiovascular system; hemodynamics; mathematical model.
Abstract: Basilar artery aneurysm occurrence is associated with hemodynamics instability. For prediction of aneurysm various experimental and numerical methods are developed. However, the greatest interests are mathematical methods for computing the hemodynamic parameters at the bifurcation of the basilar artery. In this paper, we propose a partial derivatives mathematical model of hemodynamic to compute the three-dimensional blood flow in the basilar artery. The model is possible to compute the basic hemodynamic parameters of blood flow and find the value of wall shear stress, which has a significant influence on the formation and development of aneurysms of the basilar artery. The developed model can also be used in multiscale models of hemodynamics to combine the hemodynamic models of different levels of detail.
Introduction. Basilar artery, located in the pons, formed by the junction of the vertebral arteries. The main disease of basilar artery is aneurysm which is formed at the bifurcation region of basilar artery (Fig. 1). Currently, the causes of aneurysm development are not fully understood. It is assumed that the main factors contributing to
the formation of aneurysms are hemodynamic instability and geometrical parameters of the basilar artery [1]. As a result of hemodynamic instability in basilar artery there is a region with a low wall shear stress, which leads to the destruction of the inner layer of the vessel wall and the loss of elasticity [2]. As a result, in the vessel with a weakened wall aneurysm is formed. Especially dangerous is a rupture of the aneurysm, which leads to subarachnoid hemorrhage.
One of the topical tasks today is to develop methods for predicting the Fig. 1. Angiogram presents emergence and development of
the basilar artery aneurysm [3] aneurysms °f the basilar artery. The m°i5t
* The reported study was supported by the Supercomputing Center of Lomonosov Moscow State University [14].
promising method for such problems is the use of mathematical modeling of blood flow in the basilar artery.
Methods. Strictly speaking, the blood is a suspension, however, given that the diameter of the basilar artery is relatively large (> 1 mm), it can be accepted that blood is a Newtonian fluid. Vessel diameter size of 1 mm. is critical, because for vessels with diameters smaller than 1 mm observed Fahraeus-Lindqvist effect of or sigma effect [4, 5].
To describe the three-dimensional blood flow as incompressible Newtonian fluid we use the conservation momentum law [6]:
dux ~dt
dUy
duz
■ + uT
+ ux
- + u.
dux dx
duy
— + uy
dx y
■ + uT
duz
dt dx or in vector form:
- + u,
dux dy
du, _y
dy
duz
dy
- + uv
- + uz
dux
dz
duy
~dz
duz
dz
--v
- v
--v
( -.2
d ux
v dx2
C d 2
d u, dx 2
C d 2
d uz
dx2
d 2ux
d2u,
d \
dy2
d \ dz 2
d 2u ~~2
d 2u.
"iz2
+1 dP = f ■
p dx x'
1 dP . fy;
P dy
1 dP . + = fz p dz
— + uVu - vAu +1 VP = f, dt P
v =—, P
(1)
where u - blood velocity, m/s; v - kinematic blood viscosity, m /s; ^ - dynamic blood viscosity, Pa/s; p - blood density, kg/m3; P - blood pressure, Pa; f - external forces; t - time, s.
Under external forcesf in equation (1) is usually understood gravity. This value in the simulation of hemodynamics, tend to neglect (f = 0). Blood kinematic viscosity consider constant v = const. Thus, in equation (1) there are two unknowns: blood velocity u and pressure P. Consequently, it is needed to add one more equation. This equation is continuity equation
— + div pu = 0.
dt F
(2)
Since blood is modeled as incompressible Newtonian fluid, then p = const. Therefore, equation (2) can be written as
or in the vector form
du du du — + — +— = 0
dx dy dz
div u = 0.
(3)
Thus, equations (1), (3) formed a closed equation system, describing blood flow through basilar artery in 3D (3D hemodynamics model).
Consider the computational domain D (Fig. 2). Computational domain consists of basilar artery and two vessels in which basilar artery bifurcates. According to experimental data, length lBA of basilar artery vary from 24.8 to 38.5 mm. Diameter dBA of basilar artery vary from 2.7 to 4.28 mm. In bifurcation region basilar artery bifurcates on left and right posterior cerebral arteris. Bifurcation angle a vary from 30° to 180°. Length of left/right posterior cerebral artery (1^/lright) is (6.9±3.1) /
(6.8 ±2.7) mm respectively. Diameter of left /right posterior cerebral artery (dleft/ dright) is (2.2 ± 0.6)/(2.1 ± 0.7) mm respectively [7].
Left posterior cerebral artery
Basilar artery
D„
V
D
left, out
'BA
Dright, out
A
Right posterior cerebral artery
BA
V
Fig. 2. Geometrical 3D model of basilar artery
Basilar artery walls are elastic and can be stretched under the influence of moving blood on them. The increment in the radius of the vessel reaches approximately 15 % of the original [4]. However, for the basilar artery radius increment is negligible, therefore, the elasticity of the basilar artery can be neglected and considered them rigid.
For solving system of equations (1) - (3) it is necessary to apply proper initial and boundary conditions.
For the vessel wall Dw it is no-slip condition
ID,
= 0.
(4)
For inlet boundary Din, and also for outlet boundaries Dieft,out, Drightout
boundary conditions may be applied by researcher. However most promising method to apply boundary conditions is to obtain them from multiscale hemodynamics models [4, 8]. In this case simple hemodynamics models [9, 10] are used, for example global hemodynamics model of arterial tree model. Such models are less computationally costly, however only lumped parameters of blood flow can be found. As a rule, to get the boundary conditions for 3D model it is used results of artery tree model (1D model), converted according to special coupling algorithm (Fig. 3). Thus:
u ( X ¡„, y t )| D = A (uw ( x^ t) );
Din
u ( xleft, ou^ У, Zt) D = A (u1D (xleft, out, t)) ;
^left, out
u (xright, out,y,z,t) D = A (u1D (xright, out,t)),
right, out
where xin, xleft,out, xright out - coordinates along the vessel, where 3D and 1D models
coupling; u1D - value of blood velocity, computed by arterial tree model; A - coupling algorithm.
Fig. 3. Scheme of coupling of 3D basilar artery model and 1D arterial tree model
To apply proper initial conditions u0 = u(x, y, z, t) h P0 = P(x, y, z, t) when t = 0 it is needed to solve Stokes problem
[-v Auo +VPo = f;
\ (5)
[div u0 = 0.
In system (5) also as in the system (1), (3) influence of external forces is neglected (f = 0). The boundary conditions are the same as for system (1) - (3).
Results. To conduct numerical experiments computational domain D was constructed (Fig. 2). Parameters of computational domain are: lBA = 30 mm; dBA = 3 mm; a = 120°; l[eft = 7 mm; lright = 7 mm; d[eft = 2 mm; dright = 2 mm. Blood parameters are: v = 3.3-10-6 m2/s; ^ = 0.003 Pa/s; p = 1050 kg/m3.
For solving the mathematical model it was used Finite Element Method (FEM). Computational domain was divided on 18984 finite elements (tetrahedrons). In basilar artery Reynolds number Re vary from 200 to 600 [11], hence, it can be concluded that there is a laminar flow in basilar artery. For modeling of blood flow through basilar artery on inlet boundary Din velocity profile was applied according to Poiseuille's law [12] (Dirichlet condition):
u(^ У, ^ 0 D = Umax(t)
'\ - _Ky, z)2 ^
(ba/2)2
umax(t) = 2umn(t); (6)
r (y, z) = (y - y0)2 + (z - Z0)2,
where y0, z0 - coordinate of center of Din; umn - mean blood velocity; umax -maximum blood velocity.
umn(t) 0.564-
0.348
0.076 0.152 0.228 0.304 0.38 0.456 0.532 0.608 Fig. 4. Mean inlet blood flow velocity in basilar artery
0.684
t, s
According to [13] mean blood velocity in basilar artery is 0.397 m/s. For modeling it was used the following function of mean (Fig. 4).
For outlet boundaries Aeft,out> A-ight,out was applied free flowing conditions
(Neumann condition):
du dn
(7)
du
= 0;
Aeft,<
dn
= 0,
^right, out
where n - external normal to Dleft, out, Dright, out.
For vessel wall was applied the no-slip condition (4).
On С++ programming language was developed software for computing system of equations of mathematical model of blood flow through basilar artery (1), (3) - (7) with usage of MPI technology.
Solving equations (1), (3) is high computationally costly. Hence, for solving the Lomonosov supercomputer of MSU was used [14]. As of January 2014 Lomonosov supercomputer by performance takes:
- 1 place in Top50 list (supercomputers list of CIS) [15];
- 37 place in Top500 list (world supercomputers list) [16];
- 43 place in Graph500 (world supercomputers list, oriented on solving DIC (Data Intensive Computing) problems) [17].
Supercomputer consists of more than 5000 computational nodes with peak performance of 1.7 Pflops. Main processor for computational nodes is Intel Xeon 5570 Nehalem [18]. CPU frequency is 2.96 ГГц, number of cores - 4, number of threads - 8. Computational nodes are joined by QDR InfiniBand bus.
Software runs with following parameters: N = 256, n = 2048, where N - numbers of computational nodes; n - number of MPI-processes.
Time period of t e [0; 1] s was modeled with time discretization step dt = 0.0001s. Computation results show on Fig. 5.
d)
Fig. 5. Basilar artery bifurcation region:
a - pressure distribution at t = 0; b - velocity field at t = 0; c - basilar artery perpendicular slices at t = 0.133; d - basilar artery longitudinal slices at t = 0.133
Discussion. Thus, mathematical model of blood flow through basilar artery was developed. Using this model hemodynamics parameters and wall shear stress can be computed. This information can be used for forecasting of basilar artery aneurysm emergence and development.
References
1. Steiger H.J., Liepsch D.W., Poll A., Reulen H.J. Heart and Vessels, 1988, no. 4, pp. 162-169.
2. Boussel L., Rayz V., McCulloch C., Martin C., Acevedo-Bolton G., Lawton M., Higashida R., Smith W., Young W., Saloner D. Stroke, 2008, vol. 39, no. 11, pp. 2997-3002.
3. Aneurysms Radiology & Operative Images of Brain, available at: www.brain-aneurysm.com/ (accessed 13 January 2014).
4. Formaggia L., Quarteroni A., Veneziani A. Cardiovascular Mathematics. Modeling and simulation of cardiovascular system, Milan: Springer-Verlag, 2009, 522 p.
5. Baskurt, O.K., Hardeman, M.R., Rampling, M.W., Meiselman, H.J. (Eds.) Handbook of Hemorheology and Hemodynamics, Amsterdam: IOS Press, 2007, 468 p.
6. Loitsyanskii L.G. Mekhanika zhidkosti i gaza (Mechanics of fluid and gas), Moscow: Drofa, 2003, 840 p.
7. Kamath S. Journal of Anatomy, 1981, vol. 133, no. 3, pp. 419-423.
8. Frolov S.V., Sindeev S.V., Lischouk V.A., Gazizova D.Sh., Liepsch D., Balasso A. Вoprosy sovremennoi nauki i praktiki. Universitet imeni V.I. Vernadskogo, 2013, no. 4(48), pp. 46-53.
9. Frolov S.V., Sindeev S.V., Lischouk V.A., Gazizova D.Sh., Medvedeva S.A. Вoprosy sovremennoi nauki i praktiki. Universitet imeni V. I. Vernadskogo, 2012, no. 2(40), pp. 51-60.
10. Frolov S.V., Makoveev S.N., Gazizova D.Sh., Lischouk V.A. Transactions of the Tambov State Technical University, 2008, vol. 14, no. 4, pp. 892-902.
11. Ravensbergen J., Krijger J.K.B., Hillen B., Hoogstraten H.W. Journal of Fluid Mechanics, 1995, no. 304, pp. 119-141.
12. Westerhof N., Stergiopulos N., Noble M. Snapshots of Hemodynamics: An Aid for Clinical Research and Graduate Education, New York: Springer, 2010, 200 p.
13. Valencia A., Guzman A., Finol E., Amon C. Journal of Biomechanical Engineering, 2006, no. 4, pp. 516-526.
14. Voevodin Vl.V., Zhumatiy S.A., Sobolev S.I., Antonov A.S., Bryzgalov P.A., Nikitenko D.A., Stefanov K.S., Voevodin Vad.V. Open Systems, 2012, no. 7, pp. 36-39.
15. Top50 List, September 2013, available at: top50.supercomputers.ru (accessed 13 January 2014).
16. Top500 List, November 2013, available at: www.top500.org (accessed 13 January 2014).
17. Graph500 List, November 2013, available at: www.graph500.org (accessed 13 January 2014).
18. X5570 Intel Xeon Processor, available at: http://ark.intel.com/ (accessed 13 January 2014).
Математическое моделирование движения крови в области бифуркации базилярной артерии
С. В. Фролов1, С. В. Синдеев1, Д. Липш2, А. Балассо3
Кафедра «Биомедицинская техника», ФГБОУ ВПО «ТГТУ»; ssindeev@yandex.ru;
кафедра механики, Мюнхенский университет прикладных наук, г. Мюнхен (Германия) (2); Междисциплинарная исследовательская лаборатория клиники Рехтс дер Изар, Технический университет Мюнхена, г. Мюнхен (Германия) (3)
Ключевые слова и фразы: аневризма; базилярная артерия; гемодинамика; математическая модель; сердечно-сосудистая система.
Аннотация: Возникновение аневризмы базилярной артерии связано с нарушением гемодинамики. Для прогнозирования возникновения аневризмы применяются различные методы натурных и численных экспериментов. Однако наибольший интерес представляют собой математические методы расчета гемо-динамических параметров в области бифуркации базилярной артерии. В работе предложена модель гемодинамики в частных производных для расчета трехмерного течения крови по базилярной артерии. С помощью модели возможен расчет основных гемодинамических параметров течения крови и вычисление пристеночного напряжения сдвига, которое оказывает существенное влияние на образование и развитие аневризмы базилярной артерии. Также возможно использование разработанной модели в многомасштабных моделях гемодинамики, позволяющих объединить модели гемодинамики разного уровня детализации.
Список литературы
1. Hemodynamic Stress in Terminal Saccular Aneurysms: A laser-Doppler Study / H. J. Steiger [at al.] // Heart and Vessels. - 1988. - № 4. - P. 162 - 169.
2. Aneurysm Growth Occurs at Region of Low Wall Shear Stress: Patient-Specific Correlation of Hemodynamics and Growth in a Longitudinal Study / L. Boussel [at al.] // Stroke. - 2008. - Vol. 39, № 11. - P. 2997-3002.
3. Aneurysms Radiology & Operative Images of Brain. - URL: www.brain-aneurysm.com.
4. Formaggia, L. Cardiovascular Mathematics. Modeling and Simulation of Cardiovascular System / L. Formaggia, A. Quarteroni, A. Veneziani. - Milan : Springer-Verlag, 2009. - 522 p.
5. Handbook of Hemorheology and Hemodynamics / O.K. Baskurt [at al.]. -Amsterdam : IOS Press, 2007. - 468 p.
6. Лойцянский, Л. Г. Механика жидкости и газа / Л. Г. Лойцянский. - М. : Дрофа, 2003. - 840 с.
7. Kamath, S. Observations on the Length and Diameter of Vessels Forming the Cicle of Willis / S. Kamath // Journal of Anatomy. - 1981. - Vol. 133, No. 3. - P. 419 - 423.
8. Development of Multiscale Hemodynamics Model for Research of Basilar Artery Circulation / S. V. Frolov [at al.] // Вопр. соврем. науки и практики. Ун-т им. В. И. Вернадского. - 2013. - № 4(48). - С. 46 - 53.
9. Четырехкамерная модель сердечно-сосудистой системы человека / С. В. Фролов [и др.] // Вопр. соврем. науки и практики. Ун-т им. В. И. Вернадского. - 2012. - № 2(40). - С. 51 - 60.
10. Модель сердечно-сосудистой системы, ориентированная на интенсивную терапию / С. В. Фролов [и др.] // Вестн. Тамб. гос. техн. ун-та. - 2008. - Т. 14, № 4. -С. 892 - 902.
11. Merging flows in an arterial confluence: the vertebrobasilar junction / J. Ravensbergen [at al.] // Journal of Fluid Mechanics. - 1995. - No. 304. - P. 119 - 141.
12. Westerhof, N. Snapshots of Hemodynamics: An Aid for Clinical Research and Graduate Education / N. Westerhof, N. Stergiopulos, M. Noble. - New York : Springer, 2010. - 200 p.
13. Blood Flow Dynamics in Saccular Aneurysm Models of the Basilar Artery / A. Valencia [at al.] // Journal of Biomechanical Engineering. - 2006. - No. 4. -P. 516 - 526.
14. Практика суперкомпьютера «Ломоносов» / Вл. В. Воеводин [и др.] // Открытые системы. СУБД. - 2012. - № 7. - C. 36 - 39.
15. %p50 : 19-я ред. рейтинга. - Режим доступа : top50.supercomputers.ru (дата обращения: 13.01.2014).
16. TOP500 : November 2013. - URL : www.top500.org.
17. Graph 500 List : November 2013. - URL : www.graph500.org.
18. X5570 Intel Xeon Processor. - URL : http://ark.intel.com/.
Mathematische Modellierung der Bewegung des Blutes auf dem Gebiet der Bifuraktion der Basilararterie
Zusammenfassung: Das Entstehen die Aneurysmen der Basilararterie ist mit der Verletzung der Hämodynamik verbunden. Für die Prognostizierung des Entstehens die Aneurysmen werden verschiedene Methoden der Außen- und numerischen Experimente verwandt. Jedoch stellen das meiste Interesse die mathematischen Methoden der Berechnung der hemodynamischen Parameter auf dem Gebiet der Bifuraktion der Basilatarterie dar. In der Arbeit ist das Modell der Hämodynamik in den privaten Ableitungen für die Berechnung der dreimeßigen Strömungen des Blutes nach der Basilararterie angeboten. Mit Hilfe des Modells ist die Berechnung der hämodynamischen Hauptparameter der Strömung des Blutes und die Berechnung der Wandanstrengungen der Verschiebung möglich, die den wesentlichen Einfluss auf die Bildung und die Entwicklung die Aneurysmen der Basilararterie leistet. Auch ist die Nutzung des entwickelten Modells in den mehrgroßzügigen Modellen der Hämodynamik, zulassend möglich, die Modelle der Hämodynamik verschiedenen Niveaus der Detaillierung zu vereinigen.
Modélage mathématique du courant de sang dans le domaine de la bifurcation de l'artère basilaire
Résumé: L'occurence de l'anévrysme de l'artère basilaire est liée à l'instabilité de l'hémodynamique. Pour la prévision de l'occurence de l'anévrysme sont employées de différentes méthodes des expériments naturels et numériques. Le plus grand intérêt présentent les méthodes mathématiques du calcul des paramètres hémodynamiques dans le domaine de la bifurcation de l'artère basilaire. Est proposé le modèle de l'hémodynamique dans les dérivées particulières pour le calcul du courant de sang par une artère basilaire. A l'aide de ce modèle est possible de calculer des essentiels paramètres hémodynamiques du courant de sang ainsi que la tension du décalage qui influence sur la formation et le développement de l'anévrysme de l'artère basilaire. Il est possible d'utiliser le modèle élaboré dans les modèles de grande échelle de l'hémodynamique permettant d'unifier les modèles de l'hémodynamique de différents niveaux des détails.
Авторы: Фролов Сергей Владимирович - доктор технических наук, профессор, заведующий кафедрой «Биомедицинская техника»; Синдеев Сергей Вячеславович - аспирант кафедры «Биомедицинская техника», ФГБОУ ВПО «ТГТУ»; Липш Дитер - PhD, профессор кафедры механики, Мюнхенский университет прикладных наук, г. Мюнхен (Германия); Балассо Андреа - PhD, профессор, инженер междисциплинарной исследовательской лаборатории неврологического центра клиники Рехтс дер Изар, Технический университет Мюнхена, г. Мюнхен (Германия).
Рецензент: Литовка Юрий Владимирович - доктор технических наук, профессор кафедры «Системы автоматизированной поддержки принятия решений», ФГБОУ ВПО «ТГТУ».
