Full Text

Введение В современном обществе наблюдается устойчивый рост сердечно-сосудистых заболеваний, развивающихся часто у людей трудоспособного возраста и порой приводящих к летальному исходу. В основе этих заболеваний чаще всего лежит нарушение анатомической и функциональной целостности кровеносных сосудов и, как следствие, нарушение тока крови по пораженным сосудам[1, 4, 8]. Повсеместное внедрение компьютерных технологий послужило толчком к стремительному развитию медицины. Появляются уникальные возможности для исследования, моделирования, выявления и прогнозирования причин многих заболеваний, в том числе сердечно-сосудистой системы и системы кровообращения [13, 26]. Система кровообращения человека является достаточно сложной, и множество факторов влияет на ее функционирование [6, 10]. Сбои в работе кровеносной системы имеют значимые негативные последствия [12, 15]. Особо уязвимыми являются места сужения, расширения и разветвления сосудов, где возникновение аневризмы или стеноза приводит к нарушению циркуляции крови [32]. Помимо прочего, при определенных ситуациях возможно увеличение аневризмы вплоть до ее разрыва, что может повлечь за собой самые неблагоприятные последствия[9]. Кровь является гетерогенной жидкостью, и в силу сложного состава (эритроциты, тромбоциты, лейкоциты, плазма, макромолекулы различных белков) и наличия, вследствие этого, особых реологических свойств (аномалия вязкости, псевдопластичность, тиксотропия и другие), проявляющихся при течении крови, ее можно отнести к неньютоновским жидкостям. Красные кровяные тельца, называемые эритроцитами, отвечают за транспортировку кислорода в ткани; тромбоциты - за регуляцию активности свертывающей системы. Всем форменным элементам крови свойственно деформироваться, собираться в скопления, что вносит значительное влияние на поведение потока крови [14, 17]. В упрощенном представлении кровь можно рассматривать как суспензию форменных элементов в физиологическом растворе. Макромолекулы белков способны собираться в клубки, распутываться и изменять ориентацию в потоке, что соответствует изменению конформации этих макромолекул в процессе переноса с потоком крови [3, 34]. С помощью численного моделирования стало возможным получение представлений о течении крови в сосудах разных диаметров с использованием различных моделей жидкостей. Изучение кровообращения как в организме в целом, так и на отдельном участке невозможно без анализа изменений реологических свойств крови. Для каждого отдельного случая эти свойства зависят от многих факторов: некоторые виды заболеваний, эмоциональное состояние и т.д. [17]. Проявление таких свойств, как аномалия вязкости, наличие продольной вязкости, зависящей от скорости продольного деформирования, и упругих свойств позволяет любую жидкую среду считать упруговязкой жидкостью. Класс упруговязких жидкостей довольно широк, и отнести к нему можно в равной степени и биологические жидкости, и полимерные растворы. Существует огромное количество моделей, адекватно описывающих поведение полимерных растворов. Выбор подходящей модели для описания течения крови зависит не только от заявленных свойств потока, но и от межмолекулярного взаимодействия. Кровь можно представить как суспензию форменных элементов в плазме. Неньютоновские свойства, проявляемые при течении крови по сосудам, могут быть отнесены к форменным элементам. Именно им свойственно образовывать скопления, которые ориентируются в потоке. В то же время образовавшиеся скопления не склеиваются между собой. При выборе модели это необходимо учитывать. Таким образом, взаимодействие эритроцитов и макромолекул белка с плазмой может быть соотнесено с взаимодействием макромолекул полимера и растворителя. Сложность моделирования течения крови, прежде всего, вызвана наличием особых ее свойств. При наличии заболеваний концентрация форменных элементов изменяется, соответственно, меняются и реологические свойства крови. Таким образом, реологическое поведение крови может быть описано в рамках соотношения FENE-P [2], полученного для моделирования течений разбавленных полимерных растворов. Сходство в поведении макромолекул полимера и белков позволяет применять это соотношение и для моделирования течения крови. Целью настоящей работы является исследование эффекта потери симметрии течения крови в сливающемся симметричном элементе кровеносной системы с аневризмой с помощью реологического конститутивного соотношения FENE-P. Постановка задачи Аневризмы представляют собой выпячивание артериальной стенки. Как правило, они локализуются в области разветвления/слияния сосуда или в устьях крупных ветвей, отходящих от артерий. Наиболее часто патология встречается в переднем отделе артериального круга большого мозга (виллизиев круг). Дегенеративные и атеросклеротические изменения стенки сосуда влияют на возникновение аневризм, однако стоит отметить, что они могут быть и врожденными [46]. Форма аневризмы варьируется от еле заметной припухлости до раздутого пузыря. Наиболее часто встречаются мешотчатые формы, а наиболее редко - веретенообразные. По величине различают аневризмы от 3 мм и более 25 мм. Наличие такой патологии может стать причиной внутричерепного кровоизлияния в любой период жизни человека вне зависимости от возраста[5-7, 22]. В качестве модели кровеносного сосуда с аневризмой рассматривается сливающийся канал с двумя видами каверн: прямоугольной и круглой, где высота каверны в два раза меньше ширины канала. В соответствии с тем, что стенки сосудов головного мозга деформируются незначительно [22], принято следующее допущение: рассматривается плоское течение в канале с твердыми стенками. Для представления сливающегося сосуда выбран плоский Т-образный канал с двумя видами каверн (рис. 1). Сетка расчетной области - неравномерная, со сгущением 1:300 на подходе к центральной области, с количеством узлов 40×20 в рукавах канала (рис. 2). Рис.1. Схема канала Рис. 2. Расчетные сетки с количеством узлов 60×20 и 40×20 соответственно В тестовом режиме были проведены расчеты и на других сетках с разным количеством узлов и разными сгущениями. Как показали тестовые расчеты, для формирования профиля скорости вблизи области разветвления длину канала следует выбирать не менее 10 ширин. Изотермические течения неньютоновских упруговязких жидкостей описываются уравнениями движения и неразрывности: (1) (2) где - девиатор напряжений, - плотность жидкости. Для ньютоновской жидкости девиатор напряжений может быть записан в следующем виде: (3) Ньютоновская жидкость характеризуется постоянной сдвиговой и продольной вязкостью, а также отсутствием упругих свойств. В данном случае рассматривается течение жидкости, которое описывается моделью FENE-P. Были приняты следующие начальные и граничные условия. На входе в канал (Г1) задается постоянство скоростей: (4) Для формирования установившегося профиля скорости, т.е. для предотвращения влияния заданного профиля скорости во входном сечении на течение в центральной части канала, длина входной части канала выбрана как 10 ширин. При таком выборе, как показали тестовые расчеты, течение жидкости в центральной части канала не зависит от задания профиля скорости во входном сечении. На выходе из канала (Г2) задаются условия , , (5) Данные граничные условия задают определенное значение давления на выходе из канала и определяют плоское течение. Для установления течения длина входного и выходного каналов, в соответствии с тестовыми расчетами, также была выбрана как 10 ширин. На твердых стенках (Г3) задается условие прилипания (6) Метод решения Для моделирования течения крови в местах слияния кровеносной системы использован метод контрольных объемов, который обладает свойством консервативности, характеризующимся выполнением интегральных балансовых соотношений. В настоящей работе результаты получены с помощью программной среды Open Source Field Operation and Manipulation (Open FOAM) [49]. При помощи метода Гаусса интегралы по контрольному объему сводятся к поверхностным, а значение искомой функции на поверхности ячейки получается путем интерполяции ее значений в центроидах соседних ячеек. Для конвективных слагаемых дифференциальных уравнений применяется процедура QUICK. Для диффузионных слагаемых при дискретизации оператора Лапласа необходимо вычислить нормальные градиенты скорости на поверхности ячеек с учетом скоростей в центрах соседних ячеек по симметричной схеме второго порядка. Для дискретизации производной по времени использовалась полунеявная схема Эйлера первого порядка точности. Шаг по координате выбирается исходя из выполнения критерия устойчивости схемы Co < 1[49]. Рис. 3. Распределение горизонтальной компоненты скорости по центральной линии входной части канала Сходимость решения обеспечивается путем подтверждения полученных результатов на сгущающихся сетках. В тестовом режиме были проведены расчеты на равномерных и неравномерных сетках с различным количеством узлов. Было установлено, что полученные решения на сетках 60×20 и 40×20 отличаются меньше чем на 3 % (рис. 3). Поэтому сетка 40×20 была выбрана как наиболее подходящая. Определяющие уравнения Как уже отмечено выше, для описания поведения крови в сливающемся элементе кровеносной системы использована реологическая модель FENE-P, предсказывающая свойства, характерные реальному току крови [45]. Обзор моделей Максвелла, Олдроида-Б и FENE-P Неньютоновское поведение крови [21, 37, 41] обусловлено спонтанным склеиванием эритроцитов в «монетные столбики». Данные агрегации образуют трехмерные структуры даже при отсутствии потока. При низких скоростях сдвига высокая вязкость крови зависит от размера и взаимодействия «монетных столбиков». При высоких скоростях сдвига наблюдается разрушение агрегаций, приводящих к снижению вязкости. Было предложено множество моделей для описания реологического поведения крови [18, 35, 36, 38, 39, 44, 50, 51]. Модели, применяющиеся для описания биомеханического вязкоупругого поведения крови, обычно являются обобщенными моделями Максвелла с временем релаксации и вязкостью сдвига, зависящей от агрегации/разрушения «монетных столбиков» эритроцитов. В случае стационарного движения жидкости при скорости сдвига напряжение будет пропорционально данной скорости с коэффициентом , который называется вязкостью. В отличие от чисто упругих веществ, у вязкоупругого вещества есть и эластичный, и вязкий компоненты. Вязкость вязкоупругой материи позволяет веществу растягиваться в зависимости от времени [27]. Вязкоупругость - это молекулярная перестановка. Во время напряжения вязкоупругого материала, такого как полимер, части длинной полимерной цепочки меняют позиции. Это движение (или перестановка) называется ползучестью. Уравнение Максвелла, описывающее вязкоупругое поведение материала, может быть записано в виде (7) где - вязкость; - скорость сдвига; - напряжение; - модуль упругости; t - время. При прекращении сдвига напряжения экспоненциально снижаются за характерное время λ, называемое временем релаксации. В дальнейшем модель Максвелла была обобщена для произвольных деформаций с помощью модели Олдроида (Oldroyd-Bmodel) и модели FENE-P (аббревиатура от finitely extensible nonlinear elastic model, индекс P означает фамилию автора: Петерлин). Их сходство заключается в том, что обе модели полимерную молекулу рассматривают как гибкую гантель, представляющую собой бусинки, соединенные между собой пружинкой (рис. 4). Модель Oldroyd-B предсказывает упругие свойства жидкости (конечное время релаксации), а модель FENE-P - помимо упругих свойств предсказывает наличие Рис. 4. Представление «монетного столбика» эритроцитов в виде модели из бусинок, соединённых пружинкой ( - длина вектора, положение центра масс определяется как аномальной вязкости (зависимость вязкости от сдвиговой скорости) и продольной вязкости (в продольных течениях зависимость продольной вязкости от скорости продольного деформирования). Применение модели FENE-P Несмотря на то что течение крови в больших артериях сердечно-сосудистой системы человека может быть адекватно описано с использованием уравнений Навье-Стокса, реалистичное моделирование при течении в более мелких сосудах требует, чтобы надлежащим образом учитывалось изменение вязкости в радиальном направлении. Макроскопическое моделирование крови в узких трубках часто включало представление течения крови как концентрированного ядра красных кровяных телец, окруженного более узким кольцом с меньшей вязкостью вблизи стенки сосуда [42, 43]. Трехмерное моделирование с использованием метода immersed finite element method было выполнено в работе [31] для исследования зависимости вязкости крови от скорости сдвига, влияния деформируемости клеток на вязкость. Однако сравнения между экспериментами и прогнозируемой кажущейся вязкостью при уменьшении диаметра трубки не проводились. Подход в работе [20] заключался в том, чтобы получить и использовать модели сплошной среды на макроуровне для течения крови с учeтом отклика на мезоуровне [36, 40]. Owensetal. [36] показали, что предсказания петель гистерезиса сдвигового напряжения в вискозиметре Куэтта находились в тесном согласии с экспериментальными данными [23]. Фанг и Оуэнс [28] использовали микроструктурную модель Оуэнса [36] для исследования устойчивого, колебательного и пульсирующего потока крови в трубке радиусом 0,43 мм. Упругие эффекты были наиболее выраженными при низких скоростях и амплитудах течения, что соответствовало данным из работы [47]. В работе [25] было произведено сравнение моделей Олдроида-Б, FENE-P и модели Оуэнса для моделирования течения крови. Были показаны пределы применимости моделей в зависимости от числа Вайсенберга. Также было показано, что модели учитывают нелинейные вязкоупругие эффекты, присущие крови. Методика нахождения параметров модели FENE-P с использованием реологических экспериментальных исследований представлена в работе [48]. Приложения для исследования течений жидкости FENE-Pпри ламинарном и турбулентном режимах представлены в работах [33, 40]. Исходя из вышесказанного, поскольку мы рассматриваем течение в артериях головного мозга, размер которых меньше по сравнению с крупными артериями, и существует необходимость учeта вязкоупругих эффектов, которые, как было показано в работах [19, 24, 30], оказывают значительное влияние на характер течения (см. раздел «Результаты»), в данной работе применяется модель FENE-P. Определяющие соотношения FENE-P Для модели FENE-P, в соответствии с принципом расщепления напряжений, общее напряжение можно записать в виде суммы (8) где - неньютоновская составляющая напряжения; - ньютоновская составляющая напряжения. Конститутивное реологическое соотношение запишется в следующем виде: (9) (10) (11) где - вектор скорости; - характерное время релаксации напряжений; - вязкость крови при нулевой сдвиговой скорости (); - динамическая вязкость неньютоновской составляющей жидкости при нулевой сдвиговой скорости, - динамическая вязкость плазмы; - вектор конфигурации (макромолекулы белка представляются в виде механической модели гибкой гантели: две бусинки, соединенные между собой пружинкой заданной жесткости (), где и - радиусы-векторы, определяющие положение бусинок гантели); - тензор конфигурации (), где - диадное произведение векторов конфигурации, - осреднение по ансамблю, где - вероятность того, что случайно выбранная макромолекула имеет заданный размер, находящийся в интервале от до ;- тензор скоростей деформации (), где - процедура транспонирования. С помощью стандартной процедуры приведения к характерным масштабам уравнения движения записываются в безразмерном виде и содержат следующие безразмерные величины: где - характерная скорость; - характерный линейный масштаб; - длина вектора конфигурации макромолекулы, находящейся в равновесном состоянии; - максимально возможная длина вектора конфигурации. Число Вайсенберга We характеризует отношение упругих свойств жидкости к вязким, число Рейнольдса Re - отношение инерционных свойств к вязким, коэффициент ретардации β - отношение вязкости растворителя к нулевой сдвиговой вязкости, параметр L2 характеризует степень распутывания макромолекул белка в потоке. Известно, что значения средней скорости крови в различных частях кровеносной системы могут варьироваться от 0,0007 до 0,7 м/св зависимости от расположения элемента кровеносной систем. Соответственно, и число Рейнольдса может изменяться от 0,001 до 4500 [7]. Для исследования модели течения крови в сливающемся элементе кровеносной системы использовалась феноменологическая аналогия, связанная с близкими реологическими свойствами крови и разбавленного полимерного раствора. Для решения задачи в данном модельном случае были выбраны следующие значения как физических, так и безразмерных параметров. Характерный поперечный размер элемента кровеносной системы может изменяться от долей миллиметра до нескольких сантиметров [11]. В настоящей работе характерный размер l принят как 0,0002 м. Параметры жидкости, используемые при вычислениях Параметр Величина vx 0,0005 м/с 1000 кг/м3 22,9∙10-3 кг/м·с 0,1 Re 0,04 L2 10; 100 We 0,01; 0,6 Результаты Биологические жидкости, такие как плазма крови или кровь с пониженной концентрацией макромолекул, могут быть условно отнесены к ньютоновской жидкости [34]. В предельном случае при значении числа Вайсенберга, стремящемся к нулю (), модель упруговязкой жидкости FENE-P хорошо описывает течение ньютоновской жидкости [3, 7, 46]. Параметр We = 0,01 означает, что упругие свойства малы по сравнению с вязкими. Для ньютоновской жидкости характерно то, что в рукавах канала при течении жидкости присутствуют только касательные напряжения, а возникающие в средней части нормальные напряжения являются следствием поворота потока жидкости. На рис. 5 представлены линии тока и первая разность нормальных напряжений в случае течения ньютоновской жидкости. Из этого рисунка можно видеть, что течение остается симметричным, а внутри каверны возникают две застойные области. Причиной такого поведения крови может стать и недостаточное содержание в крови эритроцитов. Для обеих форм каверн можно наблюдать аналогичные картины распределения нормальных напряжений (см. рис. 5). Первая разность нормальных напряжений, определяемая как, является одной из характеристик потока жидкости. Для неньютоновских жидкостей в случае сдвигового течения эта величина, в отличие от ньютоновских, имеет ненулевое значение. Очертания линий графиков распределения N1 по сечениям одинаковы для обеих форм канала. Исходя из этого можно предположить, что в данном случае форма каверны практически не влияет на распределение исследуемой характеристики потока ньютоновской жидкости. На рис. 7-8 представлены результаты моделирования течения модельной жидкости с достаточным содержанием макромолекул (We = 0,6) в канале с двумя видами каверн. Из рис. 7, а и рис 8, а можно видеть, что значение степени распутывания макромолекулы, равное L2 = 10 при We = 0,6, предсказывает сохранение симметричной формы течения модельной жидкости и незначительное сужение сходящегося потока в месте разветвления канала. Точка стагнации потока смещается внутрь каверны. Сравнивая рис. 5, рис. 7, а и рис. 8, а, можно предположить, что в случае L2 = 10 степень распутывания макромолекулы оказывается недостаточной для проявления потери симметрии потока. И, как следствие, полученная картина течения в этом случае близка к ньютоновскому случаю. Картина течения при L2 = 100 и We = 0,6 отличается от предыдущего случая. Можно видеть (см. рис. 7, б и рис. 8, б), что поток крови из одного рукава полностью занимает каверну. Кроме того, в канале с прямоугольной каверной (см. рис. 7, б) наблюдается асимметричное течение крови (асимметрия потока вправо), тогда как в канале со скругленной каверной - влево. а б Рис. 6. Распределение N1 по линии симметрии в горизонтальном (а) и вертикальном (б) сечениях а б Рис. 7. Линии тока при течении крови (We = 0,6; L2 = 10 (а) и We = 0,6; L2 = 100 (б)) в канале с прямоугольной формой каверны Возникающая асимметрия в одну из сторон не является характеристикой потока и равновероятно может быть направлена как в одну, так и в другую сторону. Моделирование, проведенное на разных сетках, на разных вычислительных машинах, может показывать отклонение как в одну, так и в другую сторону. В таких случаях картины с отклонением в разные стороны должны быть зеркальным отображением друг друга (рис. 9, а, б и рис. 10). На рис. 9 представлено распределение вертикальной компоненты скорости vy на расстоянии одной ширины канала от нижних угловых точек для канала со скругленной (см. рис. 9, а) и прямоугольной (см. рис. 9, б) формой каверны. Из подробного анализа рис. 9, а и рис. 9, б следует, что профиль скорости для обоих случаев несимметричный. В канале присутствуют асимметрия потока и асимметричное распределение скорости в сечении выходного канала вблизи места ветвления. Несмотря на это, на достаточном удалении от области слияния потоков (четыре ширины канала) линия раздела потоков жидкости проходит по линии симметрии выходного канала (см. рис. 9, в). Это свидетельствует о том, что на этом расстоянии от места слияния потоков макромолекулы успевают отрелаксировать, а их влияние на поток уменьшается. Феноменологически кровь становится ньютоновской жидкостью. И как следствие этого, графики распределения vy для канала со скругленной и прямоугольной формами совпадают. Таким образом, можно констатировать, что неньютоновские свойства крови проявляются только в местах разветвления кровеносной системы. Для составления более полного представления о поведении модели крови в сливающихся сосудах с аневризмой на рис. 10 представлено распределение разности главных напряжений в сечении, проходящем по центральной линии входной части канала. Величина () определяет степень оптической неоднородности потока, вызванной изменением конфигурации макромолекул [29]. Величины (), полученные как для канала с прямоугольной каверной, так и для канала со скругленной каверной, показывают схожие распределения, за исключением того, что пик величины () смещается вправо или влево в зависимости от формы каверны. Диапазон изменения величины разности главных напряжений одинаков для этих случаев. Зоной особого интереса является область разветвления канала. Из анализа рис. 10 можно сделать следующий вывод. В случае потери симметрии (We = 0,6; L 2= 100) пики () практически на порядок превышают максимальные значения (), полученные для We = 0,01 и We = 0,6; L2 = 10. При значениях параметров We = 0,6; L2 = 10 изменение максимальных значений разности главных напряжений превышает значения для ньютоновской жидкости (We = 0,01) примерно в полтора раза (см. рис. 10). Таким образом, в тех случаях, когда концентрация макромолекул в крови увеличена, а в кровеносной системе присутствуют аневризмы, возможно нарушение симметрии потока. При этом в центральной части канала, как следует из рис. 10, возникает значительное увеличение нормальных напряжений, что может привести как к сильному деформированию форменных элементов крови, так и к негативному воздействию на стенки кровеносной системы. Заключение В данной статье рассмотрено математическое моделирование течения крови в сливающемся кровеносном сосуде с аневризмой. Моделирование течения жидкости проводилось с помощью реологического конститутивного соотношения FENE-P, предсказывающего свойства, соответствующие реальной биологической жидкости, а именно, аномалию вязкости, переменную продольную вязкость и конечное время релаксации напряжений. Показано, что при определенном наборе значений параметров Weи L2, при постоянных значениях числа Рейнольдса и коэффициента ретардации поток крови в сливающемся симметричном элементе кровеносной системы теряет симметрию. Если принять значение числа We = 0,01, то для любого значения L2 картина течения будет симметричной, при этом увеличение степени распутывания макромолекулы L2 ведет только к увеличению застойных зон и большему сужению основного потока. При значениях этих параметров (We = 0,6 и L2 = 100) течение крови, помимо заявленных ранее эффектов, показывает нарушение симметричной формы потока в сливающемся элементе кровеносной системы. В зависимости от типа каверны асимметрия потока может возникать в разные стороны. Однако распределения характеристик потока в канале будут зеркальным отображением друг друга, а диапазоны изменения значения будут почти одинаковыми. Потеря симметрии тока крови для симметричного элемента кровеносной системы с аневризмой равновероятно возможна как в одну, так и в другую сторону. Физическое обоснование эффекта потери симметрии заключается во взаимодействии макромолекулы достаточной длины с основным током крови. В местах искривления потока возникают нормальные напряжения, обусловленные сливающейся формой канала. Под действием нормальных напряжений конформация макромолекул изменяется, а сама цепочка ориентируется в потоке. С другой стороны, ориентированные цепочки макромолекул влияют на изменение потока.

About the authors

F. Kh Tazyukov

E. R Kutuzova

B. A Snigerev

F. A Garifullin

References

Statistics

Views

Abstract - 93

PDF (Russian) - 209

Refbacks

• There are currently no refbacks.