О влиянии реологических и миогенных факторов на кровоток в резистивном сосуде

Автор: Шадрина Н.Х.

Журнал: Российский журнал биомеханики @journal-biomech

Статья в выпуске: 4 (62) т.17, 2013 года.

Бесплатный доступ

Рассмотрено квазиодномерное течение крови в одиночном сосуде с учетом эффектов Фареуса–Линдквиста, Фареуса и реакций сосудистой стенки на механические стимулы. В стенке сосуда рассматриваются три слоя: эндотелиальный, гладкомышечный и адвентиция. Предполагается, что диаметр сосуда регулируется двумя параметрами: концентрацией свободных ионов кальция и концентрацией оксида азота (NO) в гладкомышечном слое. Первый из них зависит от среднего окружного напряжения в стенке сосуда, NO продуцируется эндотелием под действием касательных напряжений на стенке. Для описания реологических факторов используется решение задачи о двухслойном течении в цилиндрической трубке. Исследован вклад реологических факторов в изменения кровотока при варьировании входного давления в пассивном сосуде, сосуде с реакцией Бейлисса и сосуде, обладающем чувствительностью как к реакции Бейлисса, так и касательному напряжению на стенке. При расчетах использовались данные публикаций, описывающих результаты экспериментов на мелких церебральных артериях крыс. Для сосудов этого типа показано, что при одинаковых условиях на концах сосуда радиус и кровоток в активном сосуде на десятки процентов меньше, чем в пассивном сосуде во всем исследованном диапазоне давлений, причем это различие монотонно увеличивается с ростом входного давления. Расход при давлении выше 60 мм рт. ст. снижается почти на 90%. Реологические факторы на эти различия практически не влияют. Рост кровотока за счет NO-зависимой реакции в активном сосуде может превысить 12%, эта величина зависит от давления и определяется уровнем сосудистого тонуса. Вклад реологических эффектов также зависит от давления, их влияние на кровоток в пассивном и активном сосудах относительно мало и находится в пределах нескольких процентов. Реологические факторы уменьшают снижение кровотока в «бейлиссовском» сосуде по сравнению с пассивным и повышают его в сосуде с двумя реакциями по сравнению с «бейлиссовским».

Еще

Резистивный сосуд, реакция бейлисса, no-зависимая реакция, пристеночный слой плазмы, эффект фареуса–линдквиста, эффект фареуса

Короткий адрес: https://sciup.org/146216118

IDR: 146216118

Текст научной статьи О влиянии реологических и миогенных факторов на кровоток в резистивном сосуде

Сопротивление резистивных сосудов, к которым относят мелкие артерии и артериолы, определяется сосудистым тонусом и кажущейся вязкостью крови, зависящей от показателя гематокрита, распределения и свойств форменных элементов крови, вязкости плазмы, напряжения сдвига на стенке сосуда [5]. На тонус сосуда, помимо нервных и химических факторов, влияют механические стимулы – колебания трансмурального давления и напряжения сдвига на стенке сосуда [12]. Известно,

Шадрина Нажия Хабибуллаевна, к.ф.-м.н., научный сотрудник Института физиологии имени И.П. Павлова, Санкт-Петербург что ширина свободного от форменных элементов пристеночного слоя плазмы меняется с изменением напряжения сдвига на стенке сосуда [19, 28]. В свою очередь, плазменный слой влияет на транспорт оксида азота и диффузию его в стенке, а следовательно, и на регуляцию тонуса [10]. Немаловажное значение в формировании плазменного слоя имеют деформируемость и агрегируемость эритроцитов [29], влияющие и на вязкость крови. Сложные взаимосвязи сосудистых реакций и реологических факторов обсуждались не раз [10, 16], имеются единичные примеры качественного рассмотрения условий, при которых проявляется компенсирующее или усиливающее влияние пристеночного слоя на сосудистое сопротивление [7]. Тем не менее количественные оценки не производились ни экспериментально, ни теоретически.

Ранее были получены [2, 8, 9] оценки изменений радиуса и кровотока в ответ на колебания трансмурального давления (разности между давлениями на внутренней и наружной поверхности стенки сосуда) и изменения касательного напряжения сдвига на внутренней стенке сосуда. Первый из этих эффектов, известный как эффект Бейлисса, приводит к расширению сосуда при снижении трансмурального давления и сужению при увеличении давления. Эффект обусловлен зависимостью сократительной активности гладких мышц, регулируемой концентрацией свободных ионов кальция в гладкомышечных клетках, от напряжений в сосудистой стенке. Реакция на изменения напряжения сдвига, или NO-зависимая реакция, проявляется в виде увеличения/уменьшения сосудистого просвета при росте/снижении напряжения сдвига на стенке сосуда. Эта реакция вызвана выделением оксида азота (NO) эндотелиальными клетками, выстилающими изнутри сосудистую стенку, под действием сдвигового напряжения. Считается, что NO, продуцируемый эндотелием, диффундирует в гладкомышечный слой и снижает сократительную активность гладкомышечных клеток. Поскольку целью работ [2, 8, 9] было исследование роли каждой из двух сосудистых реакций в регуляции кровотока, кровь считалась ньютоновской жидкостью, зависимость вязкости крови от радиуса сосуда не учитывалась. В предлагаемой работе вводятся в рассмотрение реологические факторы, а именно зависимость кажущейся вязкости крови от радиуса (эффект Фареуса– Линдквиста) и отличие показателя гематокрита в сосуде (динамического гематокрита) от гематокрита в сливном резервуаре (эффект Фареуса). Сравнивается величина кровотока при учете и без рассмотрения двух этих эффектов.

Постановка задачи

Рассматривается течение крови в сосуде, внутренняя поверхность которого мало отличается от кругового цилиндра и слабо меняется вдоль оси. Полагается, что радиус сосуда мал по сравнению с длиной. Известно, что течение в резистивных сосудах характеризуется малыми значениями числа Рейнольдса (Re << 1) и параметра Уомерсли W = R^ ro/v <  0,5 ( R - внутренний радиус сосуда; ш - угловая частота пульсаций градиента давления; ν – кинематическая вязкость крови). Эти условия позволяют применить квазистационарные уравнения гидродинамической теории смазки для всей области потока в сосуде (см., например, [4]).

д Q дп R 2  -          п R 4 д p

— +----= 0, Q =----.

д x     д t                  8 n app д x

Здесь Q – расход; x – продольная координата; t – время; η app – кажущаяся вязкость крови; p – трансмуральное давление. Далее в тексте по умолчанию под словом «давление» будет подразумеваться именно трансмуральное давление. Давление на внешней поверхности предполагается постоянным и равным, как правило, атмосферному давлению.

Предполагается, что регуляция просвета сосуда осуществляется концентрацией свободных ионов кальция в гладкомышечных клетках C m и средней концентрацией оксида азота C в гладкомышечном слое. При этом постулируется, что C m однозначным образом зависит от среднего окружного напряжения в стенке: с = pR / h ( h - толщина стенки сосуда), а скорость производства NO эндотелием пропорциональна модулю напряжения сдвига τ w на стенке.

Уравнения для радиуса и концентрации кальция имеют следующий вид [2, 9]:

Л , ^ + R = Л , dp + ( Fpas ( p) - Fact (Cm ) )(Ф( C) + 1) ,(2)

а-С^ = — Cm + *) + P^p.(3)

9 t9 t

Здесь ЛR, Лp, а, в   - постоянные; функция F (p) = Rmax - a —max0

p + a описывает зависимость статического радиуса пассивного сосуда от давления, Rmax ,R0 ,a – постоянные величины, определяемые аппроксимацией экспериментальных кривых, приведенных в работе [17]. Функция Fact (Cm ) отражает изменения радиуса активного сосуда, вызванные сокращением гладких мышц, она получена сплайн-аппроксимацией данных работы [17]. Множитель (Ф(C) +1) в уравнении (2) описывает регуляцию радиуса оксидом азота. Зависимость Ф( C) получена на основе литературных данных (см. [2]) и имеет вид Ф(C) = 0[1 + tanh(n lg( 4 C))]. Функция у(с) в (3) представляет собой сплайн-аппроксимацию результатов, представленных в работе [17] (см. [2]).

Концентрация оксида азота C находится из уравнений, описывающих радиальную диффузию NO из эндотелия. Полагается, что стенка сосуда состоит из трех концентрических слоев: эндотелия, гладкомышечного и наружного (адвентиции), ширина которых меняется при сокращении/релаксации гладкомышечных клеток. Содержимое сосуда разбивается на две области: внутреннее «ядро» с эритроцитами и наружный пристеночный слой плазмы. Уравнения диффузии записываются для каждой области внутри сосуда и его стенки:

9 C i = D i 5

91    r 9 r

5 C ri

5 r

- k i C/- , i = 1, 2, ..., 5.

Здесь i = 1 (0 <  r R 1 ) соответствует эритроцитарному ядру, i = 2 ( R 1 r R 2) -кольцевому плазменному слою, i = 3 ( R 2 r R 3 ) - эндотелию, i = 4 ( R 3 r R 4 ) -гладкомышечному слою, i = 5 ( R 4 r R 5 ) - адвентиции и прилегающему слою окружающей ткани шириной 6; R 2 = R , R 5 = R + h + 5 . На границе между областями выполняются условия непрерывности концентраций и потоков

C.(t,R,x) = Ci+1(t,Ri,x),    i = 1, ..., 4,(5)

9 C ( t , R , x )         9 C i 41 ( t , R i , x )

------1-------= Di+1------5--------, i =1; 4.

9 r9 r

На границах эндотелия с плазмой крови и гладкомышечным слоем условия непрерывности потока суть

5 c- ( t , R , x ) - D . 5 C + 1 ( t , R , x ) =dt| , i = 2;3; t= R d p .

9r           i + 1        9r 1 w|               w 2 9x

Значения параметров

Параметр

Значение

Единицы измерения

Источник

Λ R

40

c

[9]

Λ p

1,47∙10–7

м3∙с∙Н–1

[9]

a

3360,49

Н∙м–2

[9]

R*

98∙10–6

м

[9]

R max

116,4∙10–6

м

[9]

R 0

48,92∙10–6

м

[9]

θ

0,06

аппроксимация данных [15]

ξ

18,75∙104

м3∙моль–1

оценка по данным [15]

α

50

с

[9]

β

2,25∙10–7

моль∙с∙м–1∙Н–1

[9]

ζ

0,25∙10–9

моль∙с–1∙Н–1

см. [8]

Ζ v

8,6579∙109

Н∙с∙м–5

[8]

S w

1,112∙10–12 ( R* )2

м2

оценка по данным [17, 21]

S 3

1,21∙10–13 S w

м2

[11]

S 4

6,94∙10–13 S w

м2

[11]

D i , i = 1 ^ 3

3,3∙10–9

м2 с–1

[2]

D , , i = 4;5

8,5∙10–10

м2∙с–1

[8]

k 2

0,26

моль–1∙м3∙с–1

[18]

k 3

0,26

моль–1∙м3∙с–1

[18]

k 4

50

моль–1∙м3∙с–1

[27]

k 5

1

с–1

[26]

X 1 , X 5

1

[26]

X 2 , Х з

2

[18]

X 4

2

[27]

В уравнениях (4)–(7) C i , D i , k i , χ i концентрация, коэффициент диффузии, скорость и порядок реакции NO в i- й области, D i , k i , χ i , ζ – постоянные величины, их значения приведены в таблице.

В центре сосуда выполняется условие симметрии

8 C i ( t ,0, x ) = 0

5 r

На внешней границе r = R + h + 5 поток считается нулевым:

8 С 5 ( t , R 5, x ) = 0 d r

Управляющий параметр C находится осреднением концентрации NO в гладкомышечном слое:

C ( t , x ) =

n ( R 2 - R 2)

2 n j 0

R 4 ( t , x )

J C 4( t , r , x ) rdrd ф .

R 3 ( t , x )

Условие несжимаемости стенки и ее слоев при отсутствии продольных растяжений сводится к соотношениям

S w =n[ ( R + h )2 - R 2 ] = const, S i =n ( R i - R i - ) = const, i = 3; 4 .          (11)

Если считать кровь ньютоновской жидкостью, то n app = const, и после задания давлений на входном p in и выходном p out концах сосуда система уравнений (1)-(11) становится замкнутой. Такая постановка задачи была изложена ранее в работе [2]. Далее давление на концах сосуда будет задаваться в следующем виде:

Pn = Pa ,  Pout = Pv — ZvQv , где pa, pv – константы, Qv и Zv – расход и сопротивление на выходе из сосуда.

При учете эффекта Фареуса–Линдквиста необходимо ввести в рассмотрение зависимость кажущейся вязкости крови от внутреннего радиуса сосуда. На основе анализа большого числа публикаций и собственных исследований течения суспензий эритроцитов в стеклянных трубках авторы работы [22] вывели эмпирические соотношения, связывающие кажущуюся вязкость крови с внутренним диаметром трубки D , показателем гематокрита в сливном резервуаре H D и вязкостью плазмы η pl :

П rel

app n pl

1 + ( n rel ,0,45

- 1)

(1 H d ) M 1 (1 - 0,45) M - 1

Здесь

Пrel,0,45 = 220 exp(-1,3D) + 3,2 - 2,44 exp(-0,06D) ,(14)

M = (0,8 + exp(-0,0753D )(-1 +-----——) +-----——.

1 + 10 - 11 D 12    1 + 10 - 11 D 12

Той же группой авторов на основе анализа литературных данных была получена эмпирическая формула для вычисления динамического гематокрита H T [23]:

-^ = HD + (1 - HD )(1 +1,7 exp(-0,35 D) - 0,6 exp(-0,01D).(16)

HD app

П о

П о

R 1

П с

R

+1 -

R 1

R

- 1

H C

HD '\.:: к f R_ ) + 2 f R_ ) 1 _(R ) Н с   П о   п С I R ) I R ) I R )

Решение получено при условии, что скорости и касательные напряжения на границе между слоями непрерывны. Показатели гематокрита H D и H T определяются соотношениями

RR

HD = Q J rht (r) u (r) dr , HT = R J rht (r) dr , где ht (r), u(r) – локальные значения показателя гематокрита и скорости.

Выражение для вязкости п С получено в работе [24] из эмпирической формулы (13) в предположении, что п С = П app при D ^ 0:

п ,.. (1 - Н с Г0,8 - 1

1 + 2, 2                 л о • пpl           (1 - 0,45)-0,8 -1

Воспользуемся также предложенным авторами работы [24] способом учета неоднократно описанных в экспериментальных [19, 20, 28] и теоретических работах [7, 13] колебаний ширины плазменного слоя. Поскольку флуктуационные радиальные смещения эритроцитов вызывают дополнительную диссипацию энергии вблизи стенки, полагается, что п 0 больше вязкости плазмы и является функцией толщины пристеночного слоя и показателя гематокрита.

Таким образом, течение крови в резистивном сосуде с учетом эффектов Фареуса и Фареуса–Линдквиста описывается уравнениями (1)–(4), (10), (11), (13)–(20) и граничными условиями (5)-(9), (12). Начальные условия при t = 0 берутся из экспериментальных данных работы [17].

Методы решения и значения параметров

Система уравнений, описывающих течение ньютоновской жидкости, решалась численно описанным ранее методом [2, 9]. После перехода к безразмерным переменным и уравнениям в конечных разностях производилась линеаризация последних. Затем для каждого момента времени t = t j из конечно-разностных аналогов уравнений (1–3) и условия (12) находились значения p ( t j , x ), R ( t j , x ), C m ( t j , x ) при C = C ( x , t j_ i ) из предыдущего временного слоя. После исключения из уравнений переменных R ( t j , x ), C m ( t j , x ) вычислялись итерации p ( x , t j ) методом прогонки. Итерации проводились до тех пор, пока модуль разности между s- й и ( s + 1)-й итерациями p ( t j , x ) не оказывался меньше 10 - 7 , затем из линейных соотношений находились R ( t j , x ), C m ( t j , x ). Для нахождения C ( t j , x ) решалась задача о радиальной диффузии (4–9), (11).

После перехода к новым переменным т = t, p = r I R1, i = 1; p = (r - Ri-1)I( Ri - Ri-1), i = 2; 3; 4; p = (r - R4)I 5, i = 5, где 0 < p < 1, методом параметрической прогонки [3] в каждой из пяти областей определялась концентрация оксида азота Ci(т, р, x), т = tj, и по формуле (10) вычислялся параметр C(tj, x).

При рассмотрении неньютоновских свойств крови сначала для каждого момента времени t = t j описанным выше методом вычислялись p ( t j , x ), R ( t j , x ), C m ( t j , x ) при C = C ( x , t j- 1 ), n app ( x , t j- 1 ). Зная R ( t j , x ), из эмпирических соотношений (13-16) определяли n app ( x , t j ), H t ( x , t j ) и решали уравнения (17-20) относительно R 1 ( t j , x ), Н С ( t j , x ), n 0 и П c • Уравнения (17-20) можно свести к одному неявному уравнению относительно переменной H С ( t j , x ):

Hc = Y (Нс ),(21)

Y ( Н с ) =

_____________________ app HT I? n app - Н n C ( Н с ) +   n app - Hd n C ( H c )   - 4 n app ( Hd I H t - 2 ) -n с ( H c )

Ht           yv Ht7

Здесь

Пс(Hc ) = npi [1 - A + A (1 - Hc)-“]. A =      2'2.

( 1 - 0,45 ) - 1

Уравнение (21) решается методом итераций. Оценки показали, что отображение Hc 2 = Y ( Нс 1 ) является сжимающим, т.е. решение удовлетворяет критерию единственности [1]. После нахождения Н с определяются R 1 ( t j , x ), n 0 и n С и рассматриваются уравнения диффузии (4–9), (11) для момента времени t j с использованием вычисленных для t = tj значений R 1 ( t j , x ), НС ( t j , x ). Последний параметр влияет на коэффициент взаимодействия NO с гемоглобином в области 1 (0 <  r R 1 ).

Счет осуществлялся в среде Turbo Pascal. Полагалось, что npi = 0,0158 П, Hd = 0,45, длина сосуда l = 20R*, где R* - радиус пассивного сосуда при p = 60 мм рт. ст. Параметр к 1, характеризующий скорость реакции NO в центральной части сосуда, пропорционален HС, скорости взаимодействия NO c эритроцитарным гемоглобином и концентрации гемоглобина в одном эритроците. Если принять, что при Нс = 0,45 к1 = 1230 с-1 [14], то при изменении динамического гематокрита к 1 (с-1) пересчитывается по формуле к ,= 1230-Нс-.

1        0,45

Значения остальных параметров приведены в таблице. Для описания свойств сосудистой стенки использовались параметры и функциональные зависимости, полученные в экспериментах на мелких мозговых артериях крыс [17]. При вычислениях задавалось p a , полагалось, что p a - p v = 1 мм рт. ст., Z v = const.

Результаты вычислений представлены в виде относительных изменений сравниваемых величин, выраженных в процентах. Расчеты проводились для точки, равноудаленной от концов сосуда.

Результаты и их обсуждение

На рис. 1 представлена зависимость относительных изменений радиуса, расхода и кажущейся вязкости крови, вычисленных для «бейлиссовского» сосуда с учетом реологических эффектов, по сравнению с их значениями, полученными без рассмотрения реологических факторов, от входного давления. Относительные изменения каждого параметра Y вычисляются как A Y = ( Y +- Y - )/ Y - - 100%, где Y + и Y - - соответственно значения параметра Y , вычисленные с учетом эффектов Фареуса и Фареуса–Линдквиста и без их рассмотрения. В случае, когда реологические эффекты не рассматривались, кажущаяся вязкость крови бралась равной среднему в рассмотренном диапазоне p in значению вязкости, вычисленной с учетом этих эффектов.

Реологические эффекты практически не влияют на радиус сосуда: A R = 0 во всем рассмотренном диапазоне давлений на входе в сосуд (см. рис. 1, линия 1 ). Поскольку для приближения к физиологическим условиям разность давлений на концах сосуда бралась много меньше входного давления, изменения радиуса при переходе от статического режима к динамическому практически отсутствуют. Внутренний радиус сосуда главным образом определяется входным давлением, и его зависимость от p a аналогична статической кривой на рис. 2. Вклад эффекта Фареуса– Линдквиста в изменения кажущейся вязкости и расхода зависит от входного давления (см. рис. 1, линии 2 , 3 ). Если сравнить график функции Ap app ( p a ) (см. рис. 1, линия 3 ) и зависимость статического радиуса активного сосуда от давления (см. рис. 2), то обнаруживается положительная корреляция относительных изменений вязкости и статического радиуса. В результате зависимость относительных изменений расхода от давления имеет противоположную направленность (см. рис. 1, линия 2 ). Тот факт, что эффект Фареуса–Линдквиста не приводит к видимым изменениям внутреннего радиуса, но вызывает изменения кажущейся вязкости, можно объяснить высокой чувствительностью кажущейся вязкости крови даже к очень малым колебаниям диаметра.

Аналогичная картина относительных изменений Q и р app в зависимости от p a следует и из расчетов для сосуда с обеими механогенными реакциями (рис. 3). Как и для сосуда с реакцией Бейлисса, сначала производились вычисления с учетом реологических факторов. Кажущаяся вязкость при отсутствии реологических эффектов определялась осреднением значений кажущейся вязкости, вычисленных

Рис. 1. Зависимость относительных изменений R (линия 1 ), Q ( 2 ) и η a pp ( 3 ) за счет эффекта Фареуса–Линдквиста от входного давления в сосуде с реакцией Бейлисса

Рис. 2. Статическая зависимость внутреннего радиуса сосуда от трансмурального давления для мелкой мозговой артерии крысы (по данным [17])

Рис. 3. Зависимость относительных изменений R (линия 1), Q (2), ηapp (3) и C (4), вызванных реологическими факторами, от входного давления в сосуде с обеими реакциями в рассматриваемом диапазоне входных давлений с учетом реологических факторов и обеих реакций на механические стимулы. Если для зависимостей ΔQ(pa) и ∆ηapp (pa) на рис. 1 и 3 наблюдается лишь некоторое количественное отличие, то зависимость ∆R(pa) отличается качественно. Наиболее заметно различие проявляется при низких значениях давления: за счет реологических факторов вслед за уменьшением концентрации NO происходит снижение радиуса. В диапазоне 40–90 мм рт. ст. наличие реологических эффектов приводит к меньшим изменениям ΔC и менее заметным изменениям радиуса. Следует отметить, что изменения радиуса достаточно малы: в рассмотренных сосудах максимальное изменение, обусловленное реологическими факторами, составляет 0,1%. Если в сосуде с реакцией Бейлисса относительные изменения расхода за счет реологических эффектов связаны с изменениями кажущейся вязкости, то при наличии NO-зависимой реакции относительные изменения расхода обусловлены изменениями как радиуса, так и ηapp . Количественные оценки показали, что относительные изменения расхода по абсолютной величине не превышают 2% в сосуде с реакцией Бейлисса и 2,5% в сосуде с двумя реакциями.

В пассивном сосуде при отсутствии активации гладкомышечных клеток зависимость относительных изменений кажущейся вязкости, вызванных эффектом Фареуса–Линдквиста, от входного давления имеет монотонный характер: с ростом давления кажущаяся вязкость, как и статический радиус, увеличивается (рис. 4), кровоток, соответственно, снижается. Диапазон изменений ΔQ в пассивном и активном сосудах практически совпадает. Как и в сосуде с реакцией Бейлисса, реологические эффекты не влияют на радиус пассивного сосуда.

Сопоставление внутреннего радиуса и кровотока в сосуде с реакцией Бейлисса и пассивном сосуде, а также в сосуде с обеими реакциями и «бейлиссовском» сосуде при тех же условиях на концах сосуда проводилось с помощью следующих параметров:

A R B _ pas = ( R b - R p„ ) / R pas ' 100% - A Q b _ , = ( Q b Q . ) / Q p. -100% ,

AR NO _ B =( Rno + B - Rb ) / Rb '100%, AQno _ B =( Qno + B - Qb ) / Qb-100%, где нижние индексы «pas», «B», «NO+B» соответствуют пассивному, «бейлиссовскому» и сосуду с обеими реакциями.

На рис. 5 показаны изменения A RB pas , A Qb pas в зависимости от давления на входе в сосуд, полученные при вычислениях с учетом (линии 1 , 2 ) и без учета (линии 3 , 4 ) реологических эффектов. Снижение радиуса сосуда с активными миоцитами относительно радиуса пассивного сосуда (рис. 5, линии 1 , 3 ) составляет от 24 до 44%,

Рис. 4. Зависимость относительных изменений R (линия 1 ), Q ( 2 ) и η app ( 3 ), вызванных эффектом Фареуса–Линдкиста, от входного давления в пассивном сосуде

Рис. 5. Изменения R (линии 1 , 3 ) и Q ( 2 , 4 ) в «бейлиссовском» сосуде относительно тех же параметров в пассивном сосуде, вычисленные с учетом ( 1 , 2 ) и без учета ( 3 , 4 ) реологических факторов

Рис. 6. Изменения R (1, 3) и Q (2, 4) в сосуде с двумя реакциями относительно тех же параметров в «бейлиссовским» сосуде, вычисленные с учетом (1, 2) и без учета (3, 4) реологических факторов а снижение кровотока может достигать 90% (см. рис. 3, линии 2, 4). Реологические факторы не влияют на ∆RB_pas , но уменьшают различие расходов в «бейлиссовском» и пассивном сосудах, это уменьшение не превышает 2%.

В сосуде с обеими реакциями при тех же условиях на концах сосуда радиус и расход увеличиваются по сравнению с «бейлиссовским» сосудом (рис. 6). Если максимальное увеличение радиуса составляет менее 4%, то рост расхода может доходить до 14%. Отметим явно выраженную немонотонность зависимостей R NO B , QNO _ B от входного давления по сравнению с параметрами RB _ pas , QB _ pas . Снижение R NO B и QNO B происходит в диапазоне давлений, соответствующем падающему участку статической зависимости R ( p ). Как и на рис. 5, наличие эффектов Фареуса и Фареуса–Линдквиста заметно не сказывается на величине R NO B (линии 1 и 3 на рис. 6 совпадают) и приводит к небольшому снижению QNO B , причем величина этого снижения зависит от давления и не превышает 2%.

Расчеты для мелкой мозговой артерии крысы показали, что при тех же условиях на концах сосуда радиус и кровоток в активном сосуде на десятки процентов меньше, чем в пассивном, на всем исследованном диапазоне давлений, причем это различие увеличивается с ростом входного давления (см. рис. 5). При давлении выше 60 мм рт. ст. расход снижается почти на 90%. Реологические факторы на эти различия практически не влияют. Рост кровотока за счет NO-зависимой реакции в активном сосуде может превысить 12%, эта величина зависит от давления и, в конечном счете, при неизменной разности давлений на концах сосуда определяется уровнем сосудистого тонуса, что соответствует выводам экспериментальных работ [6, 25]. Вклад реологических факторов также зависит от давления, но он относительно мал. Эти факторы уменьшают снижение кровотока в «бейлиссовском» сосуде по сравнению с пассивным и повышают его в сосуде с двумя реакциями по сравнению с «бейлиссовским».

Степень влияния сосудистых реакций существенным образом зависит от индивидуальных свойств сосуда: от его положения в сосудистой иерархии, свойств его стенки, вида животного. Поэтому количественные оценки могут быть различны для разных сосудов. Все вычисления выполнены применительно к мелким мозговым артериям. Рассмотрение именно этих сосудов было обусловлено тем, что только для этих сосудов удалось найти максимальный набор необходимых экспериментальных данных.

Заключение

Рассмотрено изменение кровотока в одиночном сосуде под действием механических (трансмурального давления и напряжения сдвига на сосудистой стенке) и реологических факторов (эффектов Фареуса и Фареуса–Линдквиста). Оценки показали существенное влияние сосудистых реакций на сопротивление: снижение кровотока за счет реакции Бейлисса может достигать 90%, NO-зависимая реакция приводит к повышению кровотока, которое может доходить до 14% относительно кровотока в «бейлиссовском» сосуде. Изменение кровотока за счет рассмотренных реологических факторов не превышает 2,5%. Оценки сделаны применительно к мелким мозговым артериям. Степень влияния сосудистых реакций существенным образом зависит от индивидуальных свойств сосуда, поэтому количественные оценки для разных сосудов могут отличаться.

Список литературы О влиянии реологических и миогенных факторов на кровоток в резистивном сосуде

  • Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. -М.: БИНОМ. Лаборатория знаний, 2008. -413 с.
  • Бучин В.А., Шадрина Н.Х. О регуляции просвета резистивного кровеносного сосуда механическими стимулами//Известия РАН. Механика жидкости и газа. -2010. -№ 2. -С. 51-63.
  • Воеводин А.Ф., Шугрин С.М. Численные методы расчета одномерных систем. -Новосибирск: Наука, 1981. -208 с.
  • Егоров В.А., Регирер С.А., Шадрина Н.Х. Течение крови в микрососудистой сети мышцы при регуляторных реакциях. Квазистационарные задачи//Известия РАН. Механика жидкости и газа. -1993. -№ 1. -С. 137-145.
  • Левтов В.А., Регирер С.А., Шадрина Н.Х. Реология крови. -М.: Медицина, 1982. -272 с.
  • Мелькумянц А.М., Балашов С.А. Механочувствительность артериального эндотелия. -Тверь: Триада, 2005. -205 с.
  • Регирер С.А., Шадрина Н.Х. О совместном влиянии тонуса стенки и пристенного плазматического слоя на сопротивление микрососуда//Тез. докл. V Всерос. конф. по биомеханике, 29 мая -2 июня. -Н. Новгород, 2000. -С. 21.
  • Шадрина Н.Х. О пульсирующем течении крови в резистивном сосуде//Известия РАН. Механика жидкости и газа. -2012. -№ 6. -С. 75-86.
  • Шадрина Н.Х., Бучин В.А. О моделировании реакции резистивного сосуда на давление//Биофизика. -2009. -Т. 54, вып. 2. -С. 267-273.
  • Bascurt O.K., Yalcin O., Meiselman H.J. Hemorheology and vascular control mechanisms//Clin. Hemorheol. Microcirc. -2004. -Vol. 30. -P. 169-178.
  • Baumbach G.L, Walmsley J.G., Hart M.N. Composition and mechanics of cerebral arterioles in hypertensive rats//Am. J. Pathology. -1988. -Vol. 133, No. 3. -P. 464-471.
  • Davis M.J. Perspective: Physiological role(s) of the vascular myogenic response//Microcirculation. -2012. -Vol. 19. -P. 99-114.
  • Fedosov D.A., Caswel B., Popel A.S., Karniadakis G.E. Blood flow and cell-free layer in microvessels//Microcirculation. -2010. -Vol. 17, No. 8. -P. 615-628.
  • Kavdia M., Popel A.S. Contribution of nNOS and eNOS-derived NO to microvascular smooth muscle NO exposure//J. Appl. Physiol. -2004. -Vol. 97, No. 1. -P. 293-301.
  • Kelm M., Feelisch M., Spahr R., Piper H., Noack E., Schrader J. Quantitative and kinetic characterization of nitric oxide and EDRF released from cultured endothelial cells//Biochem. Biophys. Res. Commun. -1988. -Vol. 154, No. 1. -P. 236-244.
  • Kim S., Ong P.K, Yalcin O., Intaglietta M., Johnson P.C. The cell-free layer in microvascular blood flow//Biorheology. -2009. -Vol. 46, No. 3. -P. 181-189.
  • Knot H.J., Nelson M.T. Regulation of arterial diameter and wall [Ca2+] in cerebral arteries of rat by membrane potential and intravascular pressure//J. Physiology. -1998. -Vol. 508, No. 1. -P. 199-209.
  • Lewis R.S., Deen W.M. Kinetics of the reaction of nitric oxide with oxygen in aqueous solutions//Chem. Res. Toxicol. -1994. -Vol. 7, No. 4. -P. 568-574.
  • Namgung B., Ong P.K., Johnson P.C., Kim S. Effect of cell-free layer variation on arteriolar wall shear stress//Ann. Biomed. Eng. -2011. -Vol. 39, No. 1. -P. 359-366.
  • Ong P.K., Namgung B., Johnson P.C., Kim S. Effect of erythrocyte aggregation and flow rate on cell-free layer formation in arterioles//Am. J. Physiol. Heart Circ. Physiol. -2010. -Vol. 298. -P. H1870-H1878.
  • Osol G., Halpern W. Myogenic properties of cerebral blood vessels from normotensive and hypertensive rats//Am. J. Physiol. Heart Circ. Physiol. -1985. -Vol. 249, No. 5. -P. H914-H921.
  • Pries A.R., Neuhaus D., Gaehtgens P. Blood viscosity in tube flow: dependence on diameter and hematocrit//Am. J. Physiol. -1992. -Vol. 263, No. 6. -P. H1770-H1778.
  • Pries A.R., Secomb T.W., Gaehtgens P., Gross J.F. Blood flow in microvascular networks: experiments and simulation//Circ. Res. -1990. -Vol. 67. -P. 826-834.
  • 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. -P. 415-428.
  • Thorin-Trescases N., Bevan J.A. High levels of myogenic tone antagonize the dilator response of small rabbit cerebral arteries//Stroke. -1998. -Vol. 29, No. 6. -P. 1194-1200.
  • Tsoukias N.M., Kavdia M., Popel A.S. A theoretical model of nitric oxide transport in arterioles: frequency-vs. amplitude-dependent control of cGMP formation//Am. J. Physiol. Heart Circ. Physiol. -2004. -Vol. 286. -P. H1043-H1056.
  • Vaughn M.W., Kuo L., Liao J.C. Estimation of nitric oxide production and reaction rates in tissue by use of a mathematical model//Am. J. Physiol. Heart Circ. Physiol. -1998. -Vol. 274, No. 6. -P. H2163-H2176.
  • Yin X., Zhang J. Cell-free layer and wall shear stress variation in microvessels//Biorheology. -2012. -Vol. 49. -P. 261-270.
  • Zhang J., Johnson P.C., Popel A.S. Effects of erythrocyte deformability and aggregation on the cell free layer and apparent viscosity of microscopic blood flows//Microvasc. Res. -2009. -Vol. 77. -P. 265-272.
Еще
Статья научная