Уравнение состояния крови при течении в мелких сосудах Текст научной статьи по специальности «Ветеринарные науки»

УДК 532.5.031+612.13

А.Е. Медведев

Уравнение состояния крови при течении в мелких сосудах

A.E. Medvedev

The Equation of the State of a Blood Flow in Small Vessels

Предложена двухфазная модель крови для описания течения как в крупных, так и в мелких кровеносных сосудах. На основе данной модели дано объяснение давно известным особенностям (эффектам) течения крови в сосудах: зависимости показателя гематокрита от диаметра сосуда; существованию бесклеточного слоя плазмы вблизи стенки сосуда; тупому (по сравнению с профилем течения Пуазейля) профилю скорости крови; зависимости вязкости крови от диаметра сосуда. Получены зависимости скорости и вязкости крови от диаметра кровеносного сосуда.

Ключевые слова: показатель гематокрита, течение Пуазейля, математическая модель, двухфазное течение, относительная вязкость.

При течении крови в мелких сосудах (менее 200 мкм) наблюдается зависимость реологических свойств крови от размера сосуда - кажущаяся вязкость крови падает с уменьшением диаметра кровеносного сосуда. Для описания течения крови в мелких сосудах необходимо создание уравнения состояния, зависящее (кроме остальных параметров) от диаметра сосуда.

Особенности течения крови. Кровь представляет собой суспензию, состоящую из плазмы (вязкая несжимаемая жидкость) и эритроцитов (двояковогнутые деформируемые диски размером 8^2,5 мкм). Одной из основных характеристик крови является показатель гематокрита Н - объемное содержание эритроцитов. Течение крови в сосудах отличается особенностями (эффектами): I) зависимость показателя гематокрита от диаметра сосуда (эффект Фареуса); II) существование пристеночного слоя плазмы без эритроцитов; III) тупой (по сравнению с профилем течения Пуазейля) профиль скорости крови; IV) вязкость крови падает с уменьшением размера сосуда (эффект Фареуса-Линдквиста.

Модель течения крови. Рассмотрим кровь как суспензию, состоящую из двух несжимаемых фаз. Первая фаза - плазма крови, вторая - эритроциты. Относительная вязкость суспензии зависит от концентрации и, согласно формуле Эйнштейна, имеет вид

^ = =1 + п (ш2 )ти2, (1)

The paper presents the two-phase model to describe blood flow in large and in small blood vessels. Based on this model the study give an explanation for long-known features (effects) of blood flow in the vessels: dependence of the hematocrit (packed cell volume) on the diameter of the vessel, existence of cell-free plasma layer near the vessel wall, obtuse (as compared with the profile of Poiseuille flow) velocity profile of a blood; dependence of the blood viscosity on the diameter of the vessel. The researcher determines dependences of the blood rate and viscosity on the diameter of a blood vessel.

Key words: hematocrit, Poiseuille flow, mathematical model, two-phase flow, relative viscosity.

где m2 = H - объемная доля эритроцитов (локальный показатель гематокрита); д, д - динамическая вязкость крови и плазмы соответственно. Для крови зависимость n от m2 зададим в виде n(m2) = = 3 +5m2.

Известно [2], что эритроциты неравномерно распределены по сечению сосуда - объемная доля эритроцитов m2 монотонно убывает от максимума m20 на оси сосуда до нуля на стенке. Решение уравнений, аналогичных уравнениям Пуазейля, но с переменной вязкостью, дает формулы для скорости крови

5

w(5) = wm„ [м (i)-м (5)], м (5) = J25/л (5)d5, (2)

0

где 5 = r/r0 - безразмерный радиус, wmax - максимальная скорость течения Пуазейля.

Скорость крови w (5) имеет более тупой профиль, по сравнению с параболическим решением Пуазейля wP (5) = wmx (i - 52 ). Это связано с тем, что

концентрация эритроцитов m2 и относительная вязкость п имеют максимум на оси сосуда и минимум -на стенке. В силу этого w (5) < wP (5) по всему сечению сосуда, лишь на стенке - w(1) = wp(\) = 0.

Для простоты примем, что распределение объемной доля эритроцитов m2 по сечению сосуда задается ступенчатой функцией:

Уравнение состояния крови при течении в мелких сосудах

(3)

\т20 при 0 1 -И ,

т2 (£) = <

[0 при 1 - к <| < 1,

где к - относительная толщина пристеночного слоя плазмы.

Эффект образования пристеночного слоя связан с поперечной миграцией эритроцитов при движении по сосуду. Для нахождения уравнения состояния крови были взяты экспериментальные данные по зависимости показателя гематокрита от диаметра сосуда (рис. 1а).

Задача нахождения уравнения состояния крови сводится к решению алгебраического уравнения на толщину пристеночного слоя к = к(ё0) и объемной доли эритроцитов т20 = т20(^0):

1 -[1 - 1/^(т20 )] X2 =

(4)

= 2Ф( Нв ) -[1 - 1/2П(т20)) х},

где Х = (1 - к )2 = НВ НВ )/т20 ; Ф( НВ ) -

функция, аппроксимирующая экспериментальные данные на рисунке 1а.

а Ь

Рис. 1. Зависимость отношения показателей гематокрита НТ/ Нв (а) и относительной наблюдаемой вязкости Г|ге1 (Ь) от диаметра сосуда ^. Значки - экспериментальные данные [1]. Линии (а) - аппроксимация экспериментальных точек [1],

линии (Ь) - результаты расчета вязкости пге1 по формуле (5)

После определения к(й?0) и т20(й?0) находим зависимость относительной вязкости П и скорости V крови от диаметра кровеносного сосуда й?0.

Относительная наблюдаемая вязкость (вязкость, наблюдаемая в капиллярных вискозиметрах) определяется как отношение пуазейлевского расхода воды Qp к расходу крови ^:

Пге1 = Qp / Q■ (5)

0.8

0.7

0.6-

0.5-

i

0.4

0.3-

0.2-

0.1

0.0

20

Проведено сравнение с известными экспериментальными данными [1-3] по относительной наблюдаемой вязкости Пге1 (рис. 1Ь), по толщине пристеночного слоя к (рис. 2а) и профилю продольной скорости крови w (рис. 2Ь). Как видно из рисунков 1 и 2, несмотря на грубое приближение профиля локального гематокрита ступенчатой функцией (3), результаты расчета по модели находятся в пределах погрешности экспериментальных измерений.

d-, (мкм)

а

r (мкм) b

Точки■

Рис. 2. a - зависимость относительной толщины пристеночного слоя плазмы h от диаметра сосуда d0.

- эксперименты из [2]. Сплошные линии - расчет по предложенной модели; б - сравнение экспериментального (точки из [3]) и расчетного по формуле (2) (сплошная красная кривая) распределения скорости крови в стеклянной трубке диаметром 54,2 мкм (HD = 0,335, градиент давления - dp/dz = 38090 Па/м3). Пунктирная кривая - скорость течения Пуазейля

Выводы. Получено уравнение состояние крови - зависимость вязкости крови, показателя гематокрита и толщины пристеночного слоя от диаметра

сосуда. Данные зависимости имеют единый вид для сосудов всех размеров и переходят в формулы течения Пуазейля при больших диаметрах сосудов.

Библиографический список

1. Pries A.R., Secomb T.W. Blood flow in microvascular networks // Handbook of Physiology: Microcirculation / ed. R.F. Tuma, W.N. Dura, K. Ley. - Academ Press, 2008.

2. Sharan M., Popel A.S. A two-phase model for flow of blood in narrow tubes with increased effective viscosity near the wall // Biorheology. - 2001. - Vol. 38.

3. Long D.S., Smith M.L., Pries A.R. et al. Micro-viscometry reveals reduced blood viscosity and altered shear rate and shear stress profiles in microvessels after hemo-dilution // Proc. Natl. - Acad. Sci. USA. 2004. - Vol. 101, №27.