УДК 624.29
Рассмотрены вопросы моделирования процесса седиментации частиц крови под действием неоднородных центробежных сил. Определены основные параметры оседания, а также показан характер движения границы плазма/эритроциты с учетом свойств крови животных
МОДЕЛИРОВАНИЕ ПРОЦЕССА ОСЕДАНИЯ ФЭК ПОД ДЕЙСТВИЕМ ЦЕНТРОБЕЖНЫХ СИЛ
Ю.Е. Мегель
Доктор технических наук, профессор* Контактный тел.: 717-42-63 E-mail Khstua@Lin.com
И.В. Возиян
Аспирант* *Кафедра кибернетики Харьковский национальный технический университет сельского
хозяйства им. П.Василенко ул. Артема 44, г.Харьков, 61002
А.И. Рыбалка
Кандидат технических наук, доцент Кафедра физики
Харьковский национальный университет радиоэлектроники
пр. Ленина, 14, Харьков 61166 Контактный тел.: 702-13-45 E-mail: fizika@kture.kharkov.ua
Актуальность темы
Оседание эритроцитов крови представляет собой сложный физико-химический процесс, зависящий от большого количества факторов: морфологической структуры самих эритроцитов, состава плазмы, в которой они взвешены, способности к агрегации, количества эритроцитов в исследуемом образце крови, а также геометрических размеров эритроцитарных образований.
Все эти факторы в совокупности определяют конечное значение параметра скорости оседания эритроцитов (СОЭ), характерное для того или иного животного, или птицы.
Исследование седиментации форменных элементов крови во времени довольно трудоемкий процесс, поэтому для ускорения седиментации эритроцитов и сокращения времени исследования может быть использовано поле неоднородных центробежных сил, при этом возникает необходимость осуществления предварительного расчета параметров предложенной методики и моделирование оседания под действием центробежных сил.
Анализ литературы
Не смотря на широкое применение в клинической диагностике показателя СОЭ, до сих пор не существует единой теории механизма оседания красной крови. Поэтому, разработан целый ряд моделей, позволяющих описать процесс седиментации эритроцитов.
Наиболее простая модель основана на законе Сток-са для оседания частиц в вязкой среде, согласно которому, частица c плотностью превышающей плотность среды, оседает под действием силы тяжести с постоянной скоростью. При этом скорость оседания пропорциональна квадрату радиуса частицы, разнице ее плотности и плотности среды, и обратно пропорциональна вязкости среды [1-2]. Следовательно, отличия в значениях СОЭ крови можно объяснить различиями вязкости плазмы, размеров оседающих частиц и их агрегатов, а также их плотностей. Постепенное замедление скорости оседания эритроцитов, наблюдаемое при длительном проведении опыта, объясняют увеличением вязкости, и плотности суспензии в нижней части капилляра [2]. Однако реальный ход оседания границы между красной кровью и чистой плазмой не
согласуется с моделями этого типа. Оседание всегда происходит неравномерно, скачками, периоды ускорения оседания сменяются замедлением. Какой-либо однообразной последовательности в чередовании ускорения и замедления нет.
Такое поведение системы не может быть объяснено в рамках моделей оседания эритроцитов, основанных на законе Стокса.
Существует двухфазная модель седиментации эритроцитов под действием гравитационных сил [35]. Модель основана на моделировании динамики процессов агрегации и уплотнения, т. к. общепризнано, что ускорение оседания в начальный момент определяется процессом агрегации, а последующее замедление вызвано уплотнением эритроцитов. При построении модели кровь рассматривают как двухфазную среду (эритроциты и плазма) и исходят из общих уравнений движения сплошной среды. Модель позволяет получить график оседания границы раздела слоя плазмы и слоя эритроцитов, определить скорость этого оседания, а также вычислить момент времени, начиная с которого сказывается влияние дна трубки и скорость оседания начинает замедляться.
На базе двухфазной модели разработаны различные модифицированные модели для оседания в неоднородном вдоль оси трубки поле сил, а также под действием магнитного поля [6-8].
При рассмотрении крови как двухфазной суспензии состоящей из частиц способных образовывать агрегаты и ньютоновской несжимаемой жидкости -плазмы крови, упускается тот факт, что образование агрегатов из эритроцитов в крови может сопровождаться захватом части несущей плазмы в состав агрегата. Поэтому, в качестве модернизации двухфазной модели разработана трехфазная модель седиментации эритроцитов в поле гравитационных сил [9].
Цель работы
Разработать математическую трехфазную модель оседания эритроцитов крови под действием неоднородных центробежных сил с учетом характерных особенностей крови вида животных и птицы.
Постановка задачи исследования
Поскольку при оседании эритроцитов агрегация имеет первостепенное значение, причем образование агрегатов сопровождается захватом части плазмы в состав агрегатов, то в основу построения модели положено рассмотрение процессов агрегации, уплотнения и запирания части плазмы в состав агрегатов. При этом кровь рассматривается как среда, состоящая из несущей несжимаемой жидкости (фаза 1) и взвешенных в ней агрегатов, состоящих из твердых несжимаемых частиц (фаза 2) и запертой жидкости (фаза 3). Кровь как суспензия характеризуется числовой концентрацией эритроцитов N объемной концентрацией Н и истинной концентрацией С.
Уравнения сохранения массы, импульсов и энергии фаз, а также уравнение для полного числа агрегатов имеют вид:
ЭР" + Эр"и"=ои & Эх ,
а duа Эр"
Р1Т=-дХ+ррс+*
НЕа дОа ра ^ = _£°_ + иа(р°Т + Яа) + №
Л Эх 1 с '
ЭN ЭШа
--1--— О
& Эх
(1) (2)
(3)
(4)
где ра - плотности фаз;
@а - скорости массообмена между фазами ( @2 =0, @3 -скорость запирания жидкости в агрегатах); иа - скорости фаз; ра - тензоры напряжений;
- объемные силы межфазового взаимодействия; Qа - поток энергии в фазе а ;
- межфазовые потоки энергии; N - числовая концентрация агрегатов; G - скорость образования агрегатов.
Осаждение эритроцитов крови является длительным процессом, который не заканчивается даже по истечении суток. Поэтому процесс оседания эритроцитов был ускорен за счет использования поля центробежных сил. Следовательно, время седиментации эритроцитов сокращено от нескольких часов до нескольких минут. Будем считать, что плотность жидкости (крови) рж, находящейся в капилляре с внутренним радиусом го и высотой заполнения капилляра кровью L оседает в осевом направлении со скоростью и, вращается вокруг продольной оси х-х с угловой скоростью ю, равной угловой скорости вращения диска центрифуги. В этом случае на жидкость действуют силы центробежного поля, создавая в ней центробежный напор рц.н, величина которого определяется по формуле:
Рц.н ="
Р,«2 (г2 - Г2)
2
(5)
где г - текущее расстояние частицы от оси вращения диска.
Из формулы (5) видно, что по мере увеличения г величина рц.н увеличивается, из-за чего на каждую частицу конечных размеров в центробежном поле действует выталкивающая сила Fв, направленная к оси вращения по радиусу. Величина выталкивающей силы для случая, когда частица имеет форму цилиндра, определяется следующим образом:
Рц.ц _ ^1рц.н1 ^2рц.н2 ,
(6)
где S1, S2 - площади поверхностей частиц, на которые действуют осредненные давления рц.н1 и рц.н2 соответственно.
Кроме выталкивающей силы Fв, на частицу имеющую массу тц действует центробежная сила, величину Fц можно определить по формуле:
1
F =—л^2Ьр ю2г
ц 4 е ге
Fс(х) - результирующая сила, действующая на частицу с объемом Vе в центробежном поле:
^ (х) = V. (р. -р, )ш2г;
(7)
ре, pf - плотность эритроцитов и плазмы крови соответственно;
Решая данную систему, процесс считают одномерным, выбрав ось х совпадающей с осью капилляра. Для оседания в конечной трубке высотой L, принимают за начало отсчета верхний конец столбца крови. В начальный момент времени, начальные и граничные условия имеют вид:
Н(0,х)=Н0 , с(0,х)=с0 , v(0,x)=v0 при 0 < х < L ,
где Н0, со, Vo - начальные объемные концентрации агрегатов, частиц и средний объем твердой фазы в агрегатах соответственно.
Считается, что в начальный момент времени (несколько секунд после установки капилляров в центрифугу) оседание отсутствует, следовательно, начальные скорости движения фаз принимаются равными и1=0, и2=0, и3=0.
Система уравнений (1-4) решается методом характеристик, что позволяет получить уравнения характеристик на начальной стадии процесса оседания, а также отображают динамику движения каждой из фаз во времени под действием поля центробежных сил и имеют следующий вид:
^ )
1
(1 - НГ ( §) ^ (р.-р£ )FC (х)
ац
(8)
(г2 + 2)Н _ 1 1 _ Н
(1 _ НГ (Н)^ (р. )FC (х)
ац
.(9)
Условия на характеристиках: dv Gv2
dt
С
-01+=0 Л ) 3С dt '
(10)
1 + (3г4 + 5)Н 3(1-Н)
аи _ нас _0 л с dt _ ",
где Н - средняя объемная концентрация агрегатов; С - средняя объемная концентрация частиц; и - средний объем твердой фазы в агрегатах; а, г1, г2 - постоянные положительные коэффициенты; ц - вязкость несущей жидкости;
На основе полученных выражений осуществляется расчет основных параметров оседания, а также динамику движения каждой из фаз во времени под действием поля центробежных сил. На рисунке 1 представлены характеристики семейств (5-6), которые берут начало в точках 0(0;0) и А(0;1), соответствующих вершине и дну капилляра.
X
1
ХВ1
Хв
ХВ2
ТВ1 ТВ ТВ2
а)
Хв1 Хв
ХВ2
1 А с
D
■ув
-
Тв Тв
б)
Рисунок 1. Характеристики системы: а — динамика СОЭ; б — изменение скорости движения границы плазма/эритроциты
Характеристики ОВD представляют собой графики движения границы раздела слоя плазмы и слоя эритроцитов, которые меняются в зависимости от начальных условий и состава исследуемого образца крови. На рисунке 1 приведены кривые ОВШ и ОВ2D для различных скоростей центрифугирования (ю1>ю<ю2). Под действием центробежной силы Fц частицы преодолевают сопротивление жидкости Fc и движутся в радиальном направлении (в направлении силовой линии поля), приобретая скорость и осаждения. Вместе с тем частица перемещается с жидкостью в капилляре, в направлении ко дну капилляра, со скоростью . Скорость оседания, используя уравнения характеристик, можно определить с учетом начальных условий:
(1 - Н0 )2 (р ^ ^ (х) 2
и = V-о> 1Не К) Л )(1+ кип3
(11)
Момент достижения максимума оседания, т. е. момент времени, начиная с которого сказывается наличие дна у капилляра и скорость начинает уменьшаться определяется согласно:
Т ш,
2кН0 (1 - Н0 )(ре-р, )Fc(x)
+1
-1
(12)
Т.о. в качестве диагностического показателя помимо стандартного значения СОЭ, может быть использован момент достижения максимальной скорости осе-
0
т
X
о
Т
Т
дания. Поскольку, экспериментальные исследования подтвердили, что при наличии патологических изменений в организме животного происходит изменение и в значении ТВ.
Выводы
В статье предложена трехфазная математическая модель позволяющая дополнительно учесть процесс запирания части плазмы крови в состав агрегатов при их образовании. Моделирование седиментации частиц крови осуществляется под действием поля неоднородных центробежных сил. Модель состоит из системы квазилинейных дифференциальных уравнений, решение которых методом характеристик позволяет получить характер динамики движения границы раздела слоя оседающих эритроцитов и слоя плазмы. Введение операции дополнительного дифференцирования позволяет получить кривую изменения скорости движения границы раздела фаз, которая имеет немонотонный характер и представляет собой многочисленные колебания скорости. На основе уравнений характеристик получены выражения для скорости седиментации эритроцитов, а также определен момент достижения максимума оседания.
Литература
1. Воейков В.Л. Физико-химические и физиологические аспекты реакции оседания эритроцитов. // Успехи физиологических наук. - 1998. - т. 29. №4. -с. 55-73.
2. Возиян И.В., Дацок О.М. Анализ моделей седиментации
эритроцитов. // Вестник НТУ ХПИ. - 2005. - Вып. 23.
- с. 162-169.
3. Лосев Е.С. К гидродинамической теории оседания эритро-
цитов. // Механика композитных материалов. - 1980.
- №1. - с. 136-141.
4. Регирер С.А. К вопросу о континуальных моделях суспен-
зий. - Прикл.математика и механика. - 1978. - т.42, №4.
- с.379-388.
5. Кизилова Н.Н., Черевко В.А. Моделирование агрегации и
оседания частиц суспензий в трубках в поле гравитационных и центробежных сил. //Материалы XV Между-нар. конф. по механике и современным прикладным программным системам. - М.: Вузовская книга. - 2007.
- С.264-265.
6. Дацок О.М., Жолонский Е.Н., Кизилова Н.Н. Двухфаз-
ная модель оседания эритроцитов в неоднородном поле сил.// Вестник ХГПУ. - 2002. - №135. - с. 36-45.
7. Кизилова Н.Н. Влияние некоторых физических полей на
механические процессы в биологических тканях: Авто-реф. Дисс..канд. физ.-мат.наук: 01.02.05/Харьк. Гос.ун-т.
- Харьков, 1993. - 18с.
8. Котельников А.А., Аристов А.А., Трифонова А.А. Моде-
лирование процесса оседания эритроцитов в капельной пробе // Современные техника и технологии: Труды XII Международной научно-практической конференции студентов и молодых ученых. - 27-31 марта 2006.
- Томск: Изд. ТПУ. - 2006 - Т. 1. - с. 324-326.
9. Лосев Е.С. Моделирование оседания агрегирующих ча-
стиц //Изв. АН СССР. МЖГ. - 1983. - №3. - с.71-78.