CC BY
20
DOI: 10.15593/RZhBiomeh/2017.3.02 УДК 532:612.3
к Российский
Журнал / Биомеханики
www.biomech.ru
ТЕЧЕНИЕ МНОГОКОМПОНЕНТНОЙ СМЕСИ В ЖЕЛУДКЕ И ДВЕНАДЦАТИПЕРСТНОЙ КИШКЕ С УЧЕТОМ ФУНКЦИОНАЛЬНЫХ НАРУШЕНИЙ: РЕЗУЛЬТАТЫ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ ДЛЯ ОПРЕДЕЛЕНИЯ КИСЛОТНОСТИ
1 Федеральный научный центр медико-профилактических технологий управления рисками здоровью населения, Россия, 614045, Пермь, ул. Монастырская, 82, email: [email protected]
2 Кафедра математического моделирования систем и процессов, Пермский национальный исследовательский политехнический университет, Россия, 614990, Пермь, Комсомольский проспект, 29
Аннотация. Работа посвящена некоторым аспектам развития математической модели для описания течения смеси в антродуоденуме. Детально представлены материалы для численного исследования распределения уровня рН в антродуоденальном отделе желудочно-кишечного тракта. В состав однофазной многокомпонентной смеси включены растворенные на молекулярном уровне компоненты: соляная кислота, гидрокарбонат натрия, вода и продукт реакции кислоты и щелочи. Учитываются функциональные нарушения четырех участков антродуоденума по выполнению секреторной и моторной функций. Математическая постановка задачи включает в себя уравнения сохранения массы и импульса для компонент смеси с учетом массовых источников за счет секреции компонент и реакции нейтрализации. На стенках антродуоденума с учетом перистальтического движения в предположении прилипания частиц среды задаются кинематические граничные условия. На остальных границах исследуемой области (сечения входа/выхода) задаются значения массовых долей компонент первой фазы, давление и условия равенства нулю касательных компонент вектора напряжений. Для определения уровня рН смеси через концентрации компонент получены два алгебраических уравнения третьего и четвертого порядков. Идентификация параметров модели выполнена на основе литературных данных и проведенных ультразвуковых исследований. В расчетных сценариях исследуется влияние изначально повышенной кислотности (например, за счет принятой кислой жидкости) на динамическое распределение кислотности в тракте в зависимости от объема и pH жидкости, наличия функциональных нарушений секреции щелочи, моторики антрума. Показано, что напитки с pH = 2,3 обладают повреждающим потенциалом для слизистого покрова стенки тракта, особенно при наличии нарушений секреции щелочи. Предложенные в статье подходы в дальнейшем могут быть использованы для предоставления индивидуальных рекомендаций по потреблению напитков и пищи с повышенной кислотностью.
Ключевые слова: численное моделирование, многокомпонентная смесь, повышенная кислотность, уровень рН, антродуоденум, функциональные нарушения.
© Камалтдинов М.Р., Трусов П.В., Зайцева Н.В., 2017
Камалтдинов Марат Решидович, к.ф.-м.н., старший научный сотрудник отдела математического моделирования систем и процессов, Пермь
Трусов Петр Валентинович, д.ф.-м.н., профессор, заведующий кафедрой математического моделирования систем и процессов, главный научный сотрудник отдела математического моделирования систем и процессов, Пермь
Зайцева Нина Владимировна, академик РАН, д.м.н., профессор, директор, Пермь
М.Р. Камалтдинов1, П.В. Трусов1'2, Н.В. Зайцева1
Введение
Повышенная кислотность в желудочно-кишечном тракте может приводить к эрозивным эффектам на эмали и дентине зубов [20, 66] и к образованию язв в желудке и двенадцатиперстной кишке [44, 57]. Чаще всего язвообразование наблюдается в антродуоденуме [1, 10], т.е. в антруме (нижнем отделе желудка, расположенном ближе к кишечнику) и в дуоденуме (двенадцатиперстной кишке).
Традиционными подходами к определению кислотности в желудочно-кишечном тракте являются экспериментальные методы. В 50-70-х гг. прошлого столетия во всем мире активно велись разработки беспроводной радиокапсулы, позволяющей определять различные характеристики среды в пищеварительном тракте, в том числе давление, температуру [49] и уровень рН [2]. Тем не менее данный метод имел ряд недостатков, связанных с невозможностью управления капсулой, сложностью определения места ее нахождения в тракте и невозможностью сопутствующего забора гистологического материала. В этой связи использование радиокапсулы ограничилось только научными целями без широкого внедрения в медицинскую практику, исключением являются рН-капсулы для диагностики гастроэзофагеальной рефлюксной болезни (заболеваний пищевода) [50]. Исследования в данном направлении возобновились в нашем веке [16, 21, 35, 40, 55, 67], в частности, разрабатываются и внедряются в практику видеокапсулы для эндоскопии тонкого и толстого кишечника [3], помогающие диагностировать болезнь Крона и некоторые другие заболевания пищеварительного тракта. Следует отметить перспективность применения радиокапсульных технологий для идентификации и верификации математических моделей течения среды в желудочно-кишечном тракте в случае преодоления упомянутых выше недостатков, так как несомненным преимуществом капсулы являются ее небольшой размер и, как следствие, малое влияние на изучаемые процессы.
Наиболее достоверным методом измерения кислотности в пищеводе, желудке и дуоденуме является рН-метрия, в основе которой лежит введение через нос или рот пациента рН-зонда, представляющего собой трубку с проводом [60]. Зонд передает электрические сигналы на регистрирующее устройство, записывающее результаты измерений уровня кислотности в нескольких точках желудочно-кишечного тракта. Другим вариантом измерений является откачка порций содержимого тракта для последующего анализа [15]. Помимо измерения концентрации кислоты, кратковременная рН-метрия позволяет определить функциональные возможности желудка и дуоденума по секреции и нейтрализации кислоты, для этого используются стимуляция гормонами и щелочные нагрузочные тесты [27, 53]. Суточная рН-метрия применяется для более длительных исследований [26, 56]: пациент в течение 24 часов носит зонд с регистрирующим прибором, который позволяет определить кислотность среды в желудочно-кишечном тракте на различных стадиях пищеварения в привычных повседневных условиях. К недостаткам метода относится ограниченное количество точек замера, сложность определения точных мест забора анализа, а также физические и психологические неудобства обследуемого. Кроме того, само наличие инородного предмета в полости тракта способно исказить результаты измерений.
Помимо непосредственных методов измерения кислотности в желудочно-кишечном тракте, существуют косвенные неинвазивные методы, например, дыхательные тесты [17, 58]. Однако неинвазивные методы не дают точных количественных оценок, что ограничивает их практическое использование [27].
К общим недостаткам экспериментальных методов измерения кислотности в тракте относятся значительная погрешность измерения, существенные временные и трудовые затраты, сложность замера в реальных условиях при различных режимах
потребления пищи. В этой связи целесообразно применять математические подходы для оценки и прогнозирования уровня pH в желудочно-кишечном тракте, которые обладают рядом преимуществ: существенная экономия материальных ресурсов, возможность легкого управления и быстрого прогнозирования, включения/выключения отдельных факторов. Кроме того, использование моделей на основе дифференциальных уравнений в частных производных позволяет оценивать и прогнозировать кислотность в любой точке области моделирования.
С развитием электронно-вычислительных мощностей и методов диагностики, позволяющих получать снимки внутренних органов с высоким разрешением, развиваются модели CFD (Computational Fluid Dynamics) для описания процессов течения жидкости со сложной реологией в разветвленных биоканалах различной формы и кривизны [41]. Тем не менее применительно к пищеварительной системе потенциал использования подходов математики и механики раскрыт не в полной мере. Одни из первых моделей для описания течения в пищеводе [45], желудке [52] и кишечнике [34] разработаны около 15-20 лет назад. В работах последних лет, посвященных трехмерному моделированию течения в замкнутом желудке, численно исследовано влияние вязкости жидкости [23], плотности и объемной доли частиц второй фазы [68] на характеристики процесса течения, намечены пути дальнейшего развития моделей - учет секреции желудочного сока и его взаимодействия с пищей [24]. В работах по двумерному моделированию течения в желудке рассмотрено распространение фермента пепсина в антральном отделе [39], двухфазное течение с учетом эвакуации в кишечник [28]. Модели пищеварения в кишечнике ограничены 1D- и 2D-постановками, но более глубоки по содержанию, например, позволяют учитывать биохимические реакции [63]. Тем не менее исследователями преимущественно рассматриваются только моторика тракта, течение однофазной неньютоновской жидкости [18, 25, 46, 48, 59]. Таким образом, количество трехмерных моделей, описывающих индивидуальное строение желудочно-кишечного тракта, довольно ограничено. В указанных работах не рассматриваются кислотность содержимого тракта, биохимические реакции, секреция пищеварительных желез и растворение пищи. Кроме того, не используются особенности пищеварительных процессов при наличии патологий.
Между тем ненормативный уровень кислотности в тракте, обусловленный, в том числе, недостаточным уровнем секреции щелочи для нейтрализации соляной кислоты, может вызывать повреждения слизистой оболочки тракта и существенно влиять на пищеварительные процессы. Например, в экспериментах показано, что внутривенное введение гидрокарбонатов способствует уменьшению вероятности образования язвы у крыс при остром воздействии HCl, а введение веществ, блокирующих выработку щелочи, приводит к увеличению вероятности язвообразования [61]. В экспериментальных исследованиях обнаружено, что при перфузии через дуоденум свиньи в течение 15 минут раствора соляной кислоты концентрацией 0,03 и 0,1 моль/л (pH = 1,5 и pH = 1) поражается 16 и 57 % площади слизистой оболочки соответственно [14]. Кроме того, в этой же работе показано, что при внутривенном введении вазоактивного интестинального пептида, стимулирующего секрецию гидрокарбонатов слизистой оболочкой, поражается только 7 % площади слизистой при воздействии соляной концентрацией 0,03 моль/л.
В силу того что организм человека представляет собой набор сложных биологических систем с наличием множества взаимосвязей, при разработке подходов для описания отдельных органов и систем необходимо стремиться к учету взаимодействий с другими системами. Указанное положение соответствует принципам фундаментальной концепции международного проекта Physiome, посвященного
проблемам описания физиологии всего организма человека [7, 31]. Для описания процессов накопления повреждений в органах и системах под воздействием факторов различной природы коллективом исследователей, в том числе авторами данной статьи, разрабатывается многоуровневая математическая модель [12]. На верхнем уровне модели описывается осредненное взаимодействие систем организма между собой. На мезоуровне (уровне органа) более детально рассматриваются повреждения участков антродуоденума под воздействием химических веществ, в том числе соляной кислоты из полости тракта. Предложена математическая модель для описания течения многофазной многокомпонентной смеси в трехмерном антродуоденуме с учетом нарушений основных функций, наличия массовых источников за счет межфазного обмена, секреции и всасывания компонент смеси, в численных экспериментах исследовано влияние функциональных нарушений и свойств фаз, описывающих пищу, на характеристики процесса течения [64].
Разработанная модель для описания течения смеси в антродуоденуме позволяет определять массовые концентрации компонент смеси, в том числе соляной и угольной кислот, гидрокарбоната натрия. Для идентификации и верификации получаемых результатов, оценок степени повреждения слизистой оболочки тракта встает задача определения уровня pH смеси, так как именно этот показатель используется в экспериментах. Взаимодействие соляной кислоты и гидрокарбонатов в желудочно-кишечном тракте осуществляется на молекулярном уровне, при этом в избытке остается либо смесь соляной и угольной кислот, либо угольной кислоты и гидрокарбонатов. Приближенные формулы для расчета уровня pH в зависимости от концентрации компонент дают значительные погрешности при малых концентрациях. В этой связи целесообразно разработать подходы, позволяющие получать более точные решения. Таким образом, статья посвящена проблемам развития математической модели антродуоденума для определения кислотности смеси в полости тракта. В данной работе детально изложены материалы для численного исследования распределения уровня pH в антродуоденуме в зависимости от наличия функциональных нарушений.
Концептуальная постановка
Рассматривается стадия пищеварения, когда содержимое желудка начинает переходить в кишечник, при этом температура смеси в желудке достигла температуры тела человека. Исследуемые процессы в первом приближении полагаются изотермическими. Так как в данной работе основное внимание уделено определению уровня pH в тракте, чтобы дополнительно не перегружать текст статьи, сделано существенное упрощение разработанной модели. Рассматривается течение однофазной многокомпонентной жидкости в антродуоденуме с растворенными на молекулярном
уровне компонентами (индекс i = 0, I ): соляной кислотой ( i = 0), продуктом реакции соляной кислоты и гидрокарбоната натрия (углекислым газом и хлоридом натрия) ( i = 1), гидрокарбонатом натрия ( i = 2), водой ( i = 3).
В модели учитывается реакция нейтрализации между соляной кислотой (HCl) и гидрокарбонатом натрия ( NaHCO3 ), в результате реакции образуется хлорид натрия NaCl и угольная кислота H2CO3, которая может разлагаться на воду H2O и диоксид углерода CO. Продукт реакции предполагается состоящим из двух компонент - воды
и смеси хлорида натрия и диоксида углерода, которые растворены на молекулярном уровне. Для описания химической реакции в уравнениях сохранения массы вводятся
массовые источники для компонент смеси. В силу малых концентраций взаимодействующих веществ тепловой эффект реакции не учитывается.
В антродуоденуме осуществляется секреция соляной кислоты и гидрокарбоната натрия. Секреторные процессы описываются с помощью объемных массовых источников вблизи границы области моделирования (стенки тракта). Так как в данной постановке компоненты пищи в смеси отсутствуют, интенсивность массовых источников кислоты постоянна и не зависит ни от каких факторов. При описании секреции гидрокарбоната натрия учитывается функциональность по передаче управляющего сигнала через нейроэндокринные механизмы. Предполагается, что рецепторы определяют уровень кислотности в теле желудке и передают команды на изменение интенсивности массового источника гидрокарбоната натрия в антродуоденуме.
На рассматриваемой стадии пищеварения для антрума характерно распространение перистальтических волн по направлению к кишечнику со значительной степенью сжатия мышц стенки (амплитудой). В модели предполагается, что моторика антродуоденума является периодической, учитываются перистальтические волны в антруме, дуоденуме, а также сокращение/сжатие пилорического сфинктера. Параметры волнового процесса полагаются не зависящими от свойств содержимого тракта.
На основе консультаций со специалистами граница области моделирования разбита на несколько подобластей по функциональному и геометрическому признаку [64] (рис. 1). В области тела желудка ЗЦ2) вводятся массовые источники соляной
кислоты. В нижней части желудка (в антруме) ЗЦ3), примыкающей к области пилорического отверстия, секретируется гидрокарбонат натрия. В области пилорического отверстия и двенадцатиперстной кишки также введены массовые источники гидрокарбоната натрия. В центральной части двенадцатиперстной кишки выделена зона ЗЦ5), которая имитирует место выхода протоков от поджелудочной
железы и печени, секреты которых также содержат гидрокарбонат натрия. , представляют собой входную и выходную границы исследуемой области,
Рис. 1. Выделение подобластей по функциональному и геометрическому признаку: ЗЦ2) - зона секреции кислоты; ЗЦ3), , - зоны секреции
гидрокарбоната натрия; ЗЦ^, ЗЦ6) - сечения входа/выхода расчетной
области [64]
в первом приближении принимаемые плоскими сечениями, свободными от касательных нагрузок.
Каждый 1-й участок стенки тракта обладает собственной поврежденностью ^(0(т) е [0; 1] по выполнению т-й функции. Функциональность участка стенки тракта
связана с поврежденностью формулой )= 1 - Е(1)(т), функциональность, отличная
от 1, является количественной мерой снижения функций тракта. В рамках данной статьи рассматриваются функциональные нарушения по секреторной (т = 1) и моторной функциям (т = 2) тракта.
Предполагается, что при нарушении секреторной функции снижается
интенсивность массового источника соляной кислоты в области . С медицинской точки зрения гипосекреция может соответствовать заболеванию атрофическим гастритом различной степени тяжести [38]. Например, при ^ 0,5 интенсивность
массового источника снижена в два раза, что соответствует средней тяжести заболевания. В случае тяжелых форм заболевания скорость секреции кислоты может снижаться до нескольких процентов от нормального физиологического уровня [29]. В рамках работы не рассматривается гиперсекреция кислоты, которая, скорее, является следствием нарушений в регуляции, в то время как пониженная секреция может быть вызвана повреждением секреторных клеток. Аналогично в модели рассматривается только снижение интенсивности массового источника гидрокарбоната натрия, т.е. нарушения (I = 3, 5) антрума, дуоденума, поджелудочной железы и печени.
Снижение уровня секреции щелочи в дуоденуме может быть связано с рядом заболеваний. Например, у пациентов с язвой дуоденума скорость секреции гидрокарбоната натрия в двенадцатиперстной кишке почти в два-три раза ниже, чем у здоровых людей [32]. Возможно, именно снижение функции секреции приводит к истончению слизистого слоя тракта и образованию язвы. Секреция гидрокарбоната натрия с соком поджелудочной железы также может уменьшаться вследствие заболеваний [19], возраста [65] или под действием негативных факторов, например, курения [36].
Нарушения моторной функции выражаются в снижении амплитуды перистальтических волн в антруме и дуоденуме, в снижении степени сжатия пилорического сфинктера, что может быть вызвано слабостью мышц циркулярного слоя стенки тракта. Другие нарушения моторики, например, изменение периодичности сокращений, не рассматриваются.
Математическая постановка
Для описания среды в антродуоденуме используем диффузионное приближение для гомогенной смеси. Предположим, что компоненты растворены и взаимодействуют на молекулярном уровне, при этом скорости относительного движения учитываются только для определения концентраций компонент, а динамическими и инерционными эффектами за счет относительного движения можно пренебречь. Относительные движения компонент определяются диффузионными механизмами. Кроме того, в уравнения сохранения массы необходимо включить массовые источники за счет процессов секреции, реакции между компонентами. Таким образом, уравнения сохранения массы для компонент смеси могут быть записаны в виде
д - — -(р7(0) + У- (руГ(г)) = -У-1(0 + + Г ер, I е[0; Т), I = 0, I, (1)
где г - радиус-вектор пространственных точек, м; р - плотность среды, кг/м ; О = ОидО - замкнутая область (внутренность всей области О и ее граница дО); - 1-я граница области, I = 1, Ь; О) - внутренность области, прилегающая к 1-й
границе, I = 1, Ь; V - скорость частиц среды, м/с; 7(г) - массовая доля г-й компоненты, X ^ = 1, г = 0, I; 1(г) - вектор интенсивности потока массы г-й компоненты за счет
I
диффузионных процессов, кг/(м2-с), г = 0, I; - интенсивность источника массы г-й компоненты за счет реакций между компонентами, кг/(м3-с), г = 0, I; £(г) -интенсивность массовых источников г-й компоненты первой фазы в области за счет секреции, кг/(м3х), г = 0, I.
Содержимое антродуоденума является гомогенной смесью - каждая компонента занимает весь ее объем, поэтому при описании компонент оперируем с массовой долей компоненты 7(г):
X 7(0 = 1, / = 0,7. (2)
г
Уравнение сохранения импульса для вязкой многокомпонентной жидкости имеет вид
д
— + У- (рvv) = -Ур + V•T + рg, г £ О, t £[0; Т), (3)
дt
где р - давление, Па; Т(у) - девиаторная часть тензора напряжений Коши, Па, которую для вязкой несжимаемой жидкости можно записать в виде
Т = Л(^ + ^)Т), г £ О, (4)
где ^ - сдвиговая вязкость, Па-с.
Для ламинарного течения вектор интенсивности потока массы г-й компоненты за
счет диффузионных процессов может быть представлен в виде
1(0 =-рК(1 )VY,), (5)
где К) - коэффициент диффузии г-й компоненты в полости желудочно-кишечного тракта, м2/с, г = 0, 7 . В первом приближении К(г) предполагается одинаковым для всех компонент.
Слагаемое источника массы за счет реакции нейтрализации пропорционально произведению концентраций соляной кислоты и гидрокарбоната натрия:
Щ) = кШ2№(г)С(0)С(2) , г = 0 3 , (6)
где к(0)(2) - константа скорости реакции, м /(с-кмоль); М(0 - молярная масса г-й
компоненты первой фазы, кг/кмоль, г = 0, 7; С(;) - молярная концентрация г-й
компоненты первой фазы в полости антродуоденума, кмоль/м3, г = 0, 7.
Соляная кислота (г = 0) и гидрокарбонат натрия (г = 2) расходуются (массовые стоки) в результате реакции (6), при этом образуются (массовые источники) вода
(/ = 3), смесь углекислого газа и хлорида натрия (/ = 1), поэтому слагаемое для
реагентов принимается отрицательным, а для продуктов реакции - положительным.
Молярную концентрацию 1-й компоненты можно вычислить, используя соотношение
С,, = 7,)р /М,). (7)
Массовые источники компонент за счет секреции зависят от участка тракта, поэтому при дальнейшем описании используем дополнительный индекс области I для источника £(г). В областях, примыкающих к границам, массовый источник за счет секреции компонент задается в виде
)(1) = Е(!)(1)ЯЦхг),Г е Ц1),1 = 2, 5, (8)
где - интенсивность массовых источников 1-й компоненты первой фазы в 1-й
подобласти , кг/(м3х); ^^^ - интенсивность массовых источников 1-й компоненты
первой фазы в 1-й подобласти при отсутствии функциональных
нарушений, кг/(м •с).
Таким образом, секреторный поток пропорционален функциональности, при )= 0 секреция полностью отсутствует, при }= 1 интенсивность источника
соответствует секреции без функциональных нарушений.
Так как в рассматриваемой постановке питательные вещества в антродуоденуме отсутствуют, то интенсивность массовых источников кислоты определяется соотношением
Xм = ?0 г = 3 4 5 (9)
Х(0)(2) ^(0)(2) , 1 3, 4, 5 , (9)
где £(°оХ2) - базовый уровень секреции кислоты в области (при отсутствии
растворенных компонент пищи), кг/(м •с).
Скорость секреции гидрокарбоната натрия определяется соотношением
8(2)(1) [Р(1)а(1)^(0) ](2 ^ (2)(1) + [Р(1)а(1)^(0) ](2)
_ 0 F S (2)(lP(1)a(1)^(0)](2) . г
S(2)(l) _S(2)(l) ^ F(l)(ne) , г т^ i ' 1 J' 4 5 ' (10)
где - базовый уровень секреции гидрокарбоната натрия (при отсутствии кислоты) в области , кг/(м •с); ^(г)(пе) - функциональность нервной системы
по передаче сигнала в область , ^(/)(пе) е [0; 1]; 5(2)(/) - константа скорости секреции гидрокарбоната натрия в области , кг/(м •с); £ '(2)(г) - вторая константа скорости секреции гидрокарбоната натрия в области , кг/м ; [р(1)а(1^(о) Ъ) - средняя массовая концентрация соляной кислоты вблизи стенки тракта Ц2), кг/м .
Соотношение вида (10) является аналогом уравнения Михаэлиса-Ментен и применяется для описания воздействий с насыщением. Соотношение (10) отражает известные экспериментальные факты [13], при этом при усилении секреции кислоты в области Ц2) происходит повышение интенсивности секреции гидрокарбоната
в областях , , для нейтрализации кислоты. В (10) присутствует
сомножитель ^г)(пе), описывающий влияние функциональности нервной системы
на передачу управляющего сигнала, так как в этом случае секреция щелочи
и измерение рецепторами кислотности среды осуществляются в разных областях антродуоденума.
Скорость секреции воды в верхней части желудка пропорциональна
скорости секреции кислоты £(3)(2) = ^зхг^охг). Скорость секреции воды в областях Ц3),
, пропорциональна скоростям секреции гидрокарбонатов £(2)(г}. Таким образом,
в первом приближении предполагается постоянство концентраций кислоты и щелочи в секретах желез.
На стенках антродуоденума с учетом перистальтического движения в предположении прилипания частиц среды задаются кинематические граничные условия и условия нулевого потока компонент через границу:
дУ ___
у(г, г) = = 0, I = 0, I, г е[0; Т), г(г) , I = 2, 5, (11)
дп )
где г(г) е - радиус-вектор материальной точки стенки антродуоденальной области желудочно-кишечного тракта. Амплитуда сокращений пропорциональна функциональности области }. На остальных границах исследуемой области
(сечения входа/выхода) задаются значения массовых долей компонент первой фазы, давление и условия равенства нулю касательных компонент вектора напряжений:
У(0 = У™, п-а-п = р, 1 -(п-а-п)п = 0, 1 = п-а, г е[0; Т), г едЦ0, I = 1, 6. (12) Система уравнений дополняется начальными условиями в виде
Y,.(/, r) = К?,, v(t, r) = v0, r gQ , t = 0.
(13)
Определение уровня ph смеси
В медицинских исследованиях для оценки функции желудка применяется щелочной тест - пациенту через рот вводят значительную концентрацию гидрокарбоната натрия и наблюдают за изменением уровня pH [5]. В этом случае уровень pH рассчитывают по приближенным формулам через значение концентрации вещества, находящегося в избытке, при этом концентрация образовавшейся углекислоты не учитывается, так как она (при быстром введении значительных концентраций гидрокарбоната) распадается на воду и углекислый газ (диссоциацией которых пренебрегаем, вероятно, при таких концентрациях углекислый газ выводится из раствора). Формулы для приближенного расчета pH при избытке соляной кислоты и гидрокарбоната натрия соответственно записываются в виде
pH = - lg
C(0) +
2 + 4K
(0) ^ 4K(w)
(14)
pH = - lg
K(d) • (K(d )C + K )
(K(2) C(2) ^ K(w))
41)
K(d) + C
K(1) ^ C(2)
где
K (d) K(1)
К,^? - вторая константа диссоциации угольной кислоты, моль/л;
K,
( w)
2/ 2
произведение воды, моль /л .
(15)
- первая константа диссоциации угольной кислоты, моль/л;
- ионное
2
В реальных условиях компоненты могут иметь малые концентрации, поэтому при расчете уровня рН для задачи течения в антродуодемуме целесообразно использовать более точные формулы с учетом концентрации угольной кислоты.
При выводе расчетных формул используются следующие основные положения.
1. В дальнейших выкладках предполагается равенство активности ионов концентрациям, поэтому уровень рН вычисляется через логарифм концентрации ионов оксония:
рН = - 1в[а(НзО+)] = - 1в[ЫзО+], (16)
где а(Н3О+) - активность ионов оксония; [Н3О+ ] - концентрация ионов оксония.
2. Предполагается независимость констант диссоциации угольной кислоты и воды от параметров среды, в основном в силу гипотезы об изотермичности рассматриваемых процессов. Таким образом, уровень рН будет определяться только концентрацией компонент в антродуоденуме.
3. Вместо первой константы диссоциации угольной кислоты используется значение кажущейся константы, учитывающей равновесие угольной кислоты с гидратом диоксида углерода.
4. Вторая ступень ионизации угольной кислоты в водном растворе не рассматривается, так как ее учет не влияет на точность оценки рН при любых концентрациях.
5. Поскольку время реакции между кислотой и гидрокарбонатом натрия намного меньше рассматриваемых временных масштабов, полагается, что в любой точке среды присутствует только одна из этих компонент, находящихся в избытке. В этой связи вывод формулы рН проведен для двух случаев, в первом случае для водного раствора гидрокарбоната натрия и угольной кислоты, во втором - для водного раствора соляной и угольной кислот.
Для вывода расчетной формулы рН смеси гидрокарбоната натрия с концентрацией С(2) и угольной кислоты с концентрацией С(1) рассматриваются
уравнения реакций, описывающие ионизацию аниона НСО- (который образовался вследствие диссоциации гидрокарбоната натрия), угольной кислоты по первой ступени и автопротолиз воды:
НСО- + НО <—> НСО + ОН-,
HCO3 + ho <—> и3°+ + CO3 , H2CO3 + HO <-> H3O+ + HCO3,
HO+HO <—> H3O+ + OH 3.
(17)
Уравнения материального баланса для системы (17) имеют вид С(2) =[НСОз ] + [ИСОз] + [СО32- ],
(18)
С(1)=[НзСОз]з + [НСО- ]з,
где [НСО3 ], [НСО] - концентрации ионов в первых двух реакциях (17); [Н2СО3]2 , [НСО- ]2 - концентрации ионов в третьей реакции (17).
В силу того что концентрация частиц, отдавших протоны, равна концентрации частиц, присоединивших протоны, уравнение электронейтральности записывается в виде
[И2С03 ] + [И30+ ] = [ОН" ] + [С032^ ] + [ИС03 ]2. (19)
Уравнения дополняются соотношениями для констант диссоциации кислоты и ионного произведения воды:
) = к«) (нсо)=[Н0+ ] [ИСО" ] / [НСО ], (2°)
к™ = к ^) (ИСО3) = [И30+ ] [С023 ] / [НСО" ], (21)
К(К) = [И30+ ][0И-]. (22) Используя второе уравнение (18), получаем:
[И2га3]=[итд + [И2С03]2 = [ИС0Д + С(2) -[ИС0" ]2. (23)
Выразим разность [ИС0] "[ИС0- ]2 из (23) и подставим в (19):
[И2С03 ] - С(1) + [И30+ ] = [0И- ] + [С03- ]. (24)
Из выражений для констант диссоциации угольной кислоты (20)-(21)
и уравнений материального баланса (18) можно получить соотношения для [И2С03], [С03- ]:
[H2CO'] = [HO* f * OS)' ] * KK [H3°' ' (25)
C * C
[CO2- ] =_^_^_K(d) K(d) (26)
[C°3 ] [H°*]2 * ] * KdKd K(1) K(2) • (26)
Подставляя (25), (26), а также [0И ] из (22) в формулу (24), после преобразований можно получить уравнение четвертой степени относительно [И30+ ]:
[Из0+ ]4 + (С(2) + К^ДО ]3 + (К*К2) -стК* -К(К))[Из0+ ]2 + +(-С(2)К™К((2)) - 2Са)К<?К((3)) - К(к)К<? )[Из0+ ] - К(к)К<?К((3)) = 0.
В растворе кислот ионизация угольной кислоты протекает по первой ступени, осуществляется автопротолиз воды, а соляная кислота с концентрацией С(0) практически полностью подвергается ионизации:
ИС1 + И20 <-> И30+ + С1- . (28)
Уравнение электронейтральности раствора кислот с учетом С(0) = с (ИС1) = [С1- ]:
[Из0+ ] = [0И- ] + [ИС0- ] + С(0). (29)
Используя соотношения для первой константы диссоциации и ионного произведения воды, после преобразований можно получить кубическое уравнение относительно [И30+ ]:
[Из0+ ]з + (К* -С(0))[Из0+ ]2 -(С(0)К^ + СтК™ + К(К))[Из0+ ] -К^К(К) = 0. (3°)
Для решения уравнений (27) и (30) в численной реализации модели использован метод касательных. Вычисление десятичного логарифма с обратным знаком от концентрации ионов оксония в соответствии с формулой (16) приводит к итоговому значению рН раствора. Следует заметить, что присутствие в растворе соли №С1 не оказывает влияния на уровень рН.
Идентификация параметров
Основные трудности идентификации вызывает определение пространственно распределенных реологических характеристик содержимого желудка на различных стадиях пищеварения. На данном этапе параметры секреции компонентов определены на основе данных литературы по максимальной и базовой скоростям секреции [13] (табл. 1). Реакция ионного обмена протекает практически мгновенно: £(ох2) > 104 м3/(кмоль-с), в расчетах использовано минимальное значение [5].
Коэффициенты диффузии К^ для компонент первой фазы приняты равными
Ы0 м/с [5]. Константы диссоциации принимаются равными К^) = 4,9-10-7 моль/л,
= 5,62-Ш-11 моль/л, значение ионного произведения воды равно
Км = 5,62-Ш-11 моль2/л2 (при температуре 35 °С) [9].
Базовые параметры моторики (при отсутствии функциональных нарушений) заданы на основе известных литературных данных [52] и результатов ультразвукового исследования [11]: периодичность - 18 с, ширина волны - 0,02 м, амплитуда - 0,009 м. Волна инициируется в теле желудка и распространяется к пилорическому отверстию со скоростью 0,0022 м/с в течение 38 с. Параметры волн в двенадцатиперстной кишке: периодичность - 9 с, ширина волны - 0,04 м, амплитуда - 0,0035 м, волна инициируется вблизи пилорического отверстия и распространяется со скоростью 0,005 м/с в течение 36 с. Открытие и закрытие пилорического отверстия осуществляется в течение 2 с каждые 18 с. Таким образом, при заданных параметрах моторика антрального отдела, пилорического отверстия и двенадцатиперстной кишки синхронизирована. Геометрия расчетной области меняется периодически.
Результаты и обсуждение
Алгоритм реконструкции трехмерной формы антродуоденума по результатам ультразвуковых исследований, а также алгоритм перепостроения расчетной сетки для описания перистальтических волн подробно изложены в предыдущих публикациях авторов [11, 4]. Представленные в данной статье результаты получены на примере среднестатистической по геометрическим характеристикам (размерам, форме) антродуоденальной области желудочно-кишечного тракта.
Исследуется влияние изначально повышенной кислотности (например, за счет принятой кислой жидкости) на динамическое распределение кислотности в тракте в зависимости от объема и рН жидкости, наличия функциональных нарушений секреции щелочи, моторики антрума (табл. 2).
В первом сценарии рассматривается случай употребления напитка с нейтральной кислотностью (рН = 7) объемом 170 мл, например, воды. Во втором и третьем сценарии рассматривается случай употребления напитка с повышенной кислотностью (рН = 2,3) объемом 100 и 170 мл соответственно. Такой уровень кислотности наблюдается у лимонного сока, «Кока-колы» и подобных газированных напитков [54]. В четвертом сценарии анализируется распределение кислотности
Таблица 1
Таблица параметров секреции, оценка по [13]
Обозначение параметра Значение параметра
Л(0)(2) ~0,404-10-7/К(2) кг/(м3-с)
Л(2)(3) ~1,75-10-7К(3) кг/(м3-с)
5(3)(3) ~7-10-7К(3) кг/(м3-с)
5 '(2)(3) ~0,5 кг/м3
5(2)(4) ~0,278-10-7К(4) кг/(м3-с)
5(2)(4) ~1,112-10-7К(4) кг/(м3-с)
5 (2)(4) ~0,5 кг/м3
5(2)(5) ~1,972-10-7/К(4) кг/(м3-с)
5(2)(5) ~7,888-10-7/К(4) кг/(м3-с)
5 (2)(5) ~0,5 кг/м3
5(3)(2) ~171,3
5(3)(3) ~395,5
5(3)(4) ~395,5
5(3)(5) ~78,3
Таблица 2
Параметры для сценариев модели
№ п/п Фунциональные нарушения секреции щелочи, ), 1 = 3,5 Фунциональные нарушения моторики антрума, ^^ Объем напитка, мл Кислотность напитка, pH
1 1 1 170 7
2 1 1 100 2,3
3 1 1 170 2,3
4 1 1 170 3,5
5 0,4 1 170 2,3
6 1 0,75 170 2,3
в желудочно-кишечном тракте при употреблении напитка с уровнем рН = 3,5 объемом 170 мл, что соответствует, например, яблочному, апельсиновому соку или другим фруктовым сокам и некоторым газированным напиткам. В сценариях 1-4 какие-либо функциональные нарушения в антродуоденуме отсутствуют, в частности,
интенсивность секреции соляной кислоты и гидрокарбоната натрия соответствует физиологическим уровням у здорового человека. В сценарии 5 рассматривается случай употребления напитка (рН = 2,3, объем 170 мл) при наличии функциональных
нарушений секреции гидрокарбоната натрия (= 0,4, I = з, 5). В сценарии 6
анализируется влияние нарушений моторики антрума ( ^(3)(3) = 0,75) на характеристики
течения жидкости в антродуоденуме при употреблении кислого напитка (рН = 2,3, объем 170 мл).
На рис. 2 отмечены четыре точки мониторинга в пристеночном слое тела желудка (точка 1), пилорического отдела (точка 2), луковицы двенадцатиперстной кишки (точка 3), в начальном сегменте нисходящей части дуоденума (точка 4).
Анализируя результаты, мы можем видеть, что во всех сценариях в полости тела желудка кислотность смеси составляет от 2 до 3,5 рН (рис. 2). Полученные результаты соответствуют физиологическим нормативам. В сценариях 1, 4 в пристеночных слоях полости антрума и двенадцатиперстной кишки образуется область меньшей кислотности (от 5 до 7 рН) для защиты от негативного воздействия (рис. 2, а, г).
В сценариях 2, 3, 5, 6 можно наблюдать области повышенной кислотности (менее 4 рН) в полости пилорического отдела желудка и дуоденума (рис. 2, б, в, д, е), что может оказывать негативное воздействие на состояние слизистой оболочки тракта. Следует отметить, что зоны с низким рН имеют различную локализацию: в сценариях 2, 3 аномально низкие рН наблюдаются в области луковицы двенадцатиперстной кишки, в сценарии 5 высокая кислотность образуется во всей верхней части дуоденума и в начальном сегменте нисходящей части двенадцатиперстной кишки, а в сценарии 6 можно видеть только небольшой участок повышенной кислотности в области луковицы двенадцатиперстной кишки.
В первом сценарии можно наблюдать, что кислотность в приграничных слоях желудка находится в пределах физиологической нормы (рис. 3, а). В пилорическом отделе желудка и луковице двенадцатиперстной кишки наблюдается небольшое повышение кислотности вследствие поступления кислотного содержимого желудка в дуоденум.
Во втором и третьем сценарии можно наблюдать значительно повышенную кислотность в пилорическом отделе желудка, луковице двенадцатиперстной кишки и в нисходящей части двенадцатиперстной кишки вследствие наличия кислого напитка (рис. 3, б, в, г). Постепенно (после 3 минут от приема напитка) уровень рН в пристеночном слое пилорического отдела желудка начинает медленно возрастать вследствие нейтрализации кислоты и перехода кислого напитка в кишечник (см. рис. 3, б). В силу поступления напитка в кишечник в области луковицы дуоденума кислотность остается очень высокой (2,3 рН) (см. рис. 3, в). В области нисходящей части двенадцатиперстной кишки уровень рН начинает медленно возрастать (спустя 3 минуты после приема напитка) вследствие нейтрализации кислоты гидрокарбонатом натрия за счет массового источника в области I = 5 (выход протоков от поджелудочной железы и печени) (см. рис. 3, г). Напиток большего объема (сценарий 2) быстрее поступает в кишечник, формируя низкие уровни рН = 2,94 уже в первые 90 с после начала моделируемой стадии пищеварения (см. рис. 3, в).
В четвертом сценарии уровень рН в пристеночном слое тела желудка и пилорического отдела близок к значениям, полученным в первом сценарии при употреблении воды (см. рис. 3, а, б). В этом случае механизмы защиты слизистой антрума и пилорического отдела справляются с нейтрализацией агрессивного воздействия кислоты. Вследствие поступления кислого содержимого в кишечник рН в дуоденуме понижается до 4 в течение 5,5 минут после приема напитка, оставаясь
Точка 4
6.73 5.15
3.58 О X 2=00 Точка 2
Точка 1
2.00
Точка 2 Точка 1
д
Рис. 2. Уровень рН, t = 330 с: а - сценарий 1; б - сценарий 2; в - сценарий 3; г - сценарий 4; д - сценарий 5; е - сценарий 6
значительно выше значений рН, полученных в сценарии 3 (см. рис. 3, в). В области нисходящей части двенадцатиперстной кишки рН восстанавливается до 5,7 вследствие нейтрализации кислоты (см. рис. 3, г), что очень близко к значениям, полученным в сценарии 1 (рН = 5,9).
При нарушении функции секреции гидрокарбоната натрия (сценарий 5) пилорический отдел желудка не может нейтрализовать аномально избыточную кислотность (рН = 2,5) на протяжении всего времени моделирования (см. рис. 3, б). В этом случае следует ожидать большую вероятность локализации повреждений слизистого покрова в желудке по сравнению с другими сценариями. Кроме того, если в сценариях 2 и 3 зона низких значений рН в верхней части двенадцатиперстной кишки формируется преимущественно вблизи вогнутой стенки дуоденума, то в сценарии 5 высокая кислотность смеси обнаруживается вблизи вогнутой и выпуклой стенок двенадцатиперстной кишки.
0 100 200 300 t,c ^"Сценарий 1 ^"Сценарий 2 ^"Сценарий 3 ^"Сценарий 4 ^"Сценарий 5 ^"Сценарий 6
а
0 100 200 300 tc ^—Сценарий 1 ^—Сценарий 2' ^"Сценарий 3 ^"Сценарий 4 ^"Сценарий 5 ^"Сценарий 6
б
0 100 200 300 t,c ^"Сценарий 1 ^—Сценарий 2 ^—Сценарий 3 ^—Сценарий4 ^"Сценарий 5 ^"Сценарий 6
0 100 200 300 t,c ^-Сценарий 1 ^—Сценарий 2 ^"Сценарий 3 ^"Сценарий 4 ^"Сценарий 5 ^"Сценарий 6
в г
Рис. 3. Динамика уровня рН вблизи стенки тракта в различных частях желудочно-кишечного тракта: а - тело желудка (точка 1); б - пилорический отдел (точка 2); в - луковица двенадцатиперстной кишки (точка 3); г -нисходящая часть двенадцатиперстной кишки (точка 4)
В сценарии 6 вследствие наличия нарушений моторной функции антрума переход кислого содержимого желудка в кишечник осуществляется медленней, чем в сценариях 1-5. Например, при открытии пилорического отверстия скорость течения в наименьшем сечении сфинктера составляет 0,01 м/с в сценарии 5 и 0,019 м/с в сценарии 1. В этой связи в сценарии 6 кишечник успешно нейтрализует кислоту, поэтому аномально низкие рН наблюдаются только в точке мониторинга № 3 (в небольшом участке луковицы двенадцатиперстной кишки).
В целом можно отметить, что при отсутствии функциональных нарушений гидрокарбонат натрия значительно эффективней нейтрализует избыточную кислотность в антродуоденуме при употреблении заданного объема напитка с рН > 3,5. Напитки с рН = 2,3 обладают существенно большим повреждающим потенциалом для слизистой оболочки стенки желудочно-кишечного тракта, так как в сценарии 3 в течение всего расчетного времени уровень кислотности в антродуоденуме остается критически высоким. При наличии нарушений секреции щелочи ситуация дополнительно усугубляется расширением площади потенциального поражения стенки тракта. Негативное воздействие кислоты может выражаться в подавлении способностей к восстановлению клеток слизистого слоя [37].
Полученные результаты по распределению кислотности в тракте согласуются с известными данными экспериментов: у здоровых людей средний уровень pH в теле желудка равен 2,7 [30], а при употреблении «Кока-колы» достоверно увеличивается время повышенной кислотности в дуоденуме [51].
Заключение
Разработанная математическая модель для описания процесса течения в антродуоденуме позволяет оценивать и прогнозировать уровень кислотности в полости тракта. Уровень pH в антродуоденуме может изменяться в зависимости от индивидуальных параметров формы и моторики тракта, наличия секреторных функциональных нарушений, объема и состава принятой пищи. В связи с этим предложенные в статье подходы в дальнейшем могут быть использованы для предоставления индивидуальных рекомендаций по потреблению напитков и пищи с повышенной кислотностью.
Используя разработанную модель, можно выделить области аномально повышенной кислотности в антральном отделе желудке и в дуоденуме. Однако для предоставления прогнозных количественных оценок развития нарушений требуется более детальное рассмотрение механизмов повреждения и восстановления стенки тракта, в том числе при воздействиях соляной кислоты, химических веществ, содержащихся в пище, лекарственных препаратов. Вполне возможно, что однократное химическое воздействие не приведет к развитию патологий. Например, в экспериментальных исследованиях показана возможность быстрой регенерации клеток слизистой стенки желудка [33, 62] и дуоденума [22, 47]. В этой связи требуется учитывать периодичность принятия пищи для рассмотрения циклов повреждения-восстановления.
Для более адекватного представления процесса течения пищи необходимо развивать модели многофазного течения с учетом упругих свойств стенок канала, регуляторных взаимодействий с другими органами и системами, в частности, с двенадцатиперстной кишкой и билиарной системой [6, 42, 43]. Актуальной проблемой остается исследование влияния реологических свойств биожидкости на характеристики процесса течения.
Благодарность
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 16-01-00126_a.
Список литературы
1. Авраменко А.А., Гоженко А.И., Гойдык В.С. К вопросу о локализации и численности язвенных дефектов, которые образуются у больных хроническим хеликобактериозом и при экспериментальном моделировании на крысах // Актуальные проблемы транспортной медицины. -2008. - T. 12, № 2. - С. 124-127.
2. Бабский Е.Б., Сорин А.М., Давыдов С.Н. Приборы эндорадиозондирования. - М.: Наука, 1975. -176 с.
3. Домарев Л.В., Русских А.Е., Шитиков Е.А. Капсульная колоноскопия. Первый опыт использования // Колопроктология. - 2015. - T. 51, № 1. - С. 27-31.
4. Камалтдинов М.Р. Трехмерное моделирование моторики антродуоденальной области пищеварительного тракта для задач оценки риска здоровью при пероральной экспозиции химических веществ // Анализ риска здоровью. - 2014. - № 2. - С. 68-77.
5. Краткий справочник физико-химических величин / под ред. А.А. Равделя, А.М. Пономаревой. - Изд. десятое, испр. и доп. - СПб.: Иван Федоров, 2003. - 117 с.
6. Кучумов А.Г. Математическое моделирование перистальтического течения литогенной желчи через проток при рубцовом стенозе, рассматриваемый в виде трубки с сужающимися стенками конечной длины // Российский журнал биомеханики. - 2016. - Т. 20, № 2. - C. 96-115.
7. Кучумов А.Г., Няшин Ю.И., Самарцев В.А., Гаврилов В.А., Менар М. Биомеханический подход к моделированию билиарной системы как шаг в направлении к построению виртуальной модели физиологии человека // Российский журнал биомеханики. - 2011. - Т. 15, № 2. - C. 32-48.
8. Невский Д.И., Ракитин А.Б., Ракитин Б.В. Анализ внутрижелудочных рН-грамм при щелочном тесте // Медицинская техника. - 2008. - № 1. - С. 23-27.
9. Новый справочник химика и технолога. Химическое равновесие. Свойства растворов. - СПб.: Профессионал, 2004. - 998 с.
10. Осадчий В.А., Буканова Т.Ю. Клинико-морфологические и патогенетические особенности эрозивно-язвенных поражений гастродуоденальной зоны у больных с различной тяжестью хронической сердечной недостаточности, ассоциированной с ишемической болезнью сердца // Сердечная недостаточность. - 2014. - Т. 15. - T. 87, № 6. - С. 374-381.
11. Трусов П.В., Зайцева Н.В., Камалтдинов М.Р. Течение в антродуоденальной области пищеварительного тракта: математическая модель и некоторые результаты // Математическая биология и биоинформатика. - 2015. - № 1. - С. 34-53.
12. Трусов П.В., Зайцева Н.В., Кирьянов Д.А., Камалтдинов М.Р., Цинкер М.Ю., Чигвинцев В.М., Ланин Д.В. Математическая модель эволюции функциональных нарушений в организме человека с учетом внешнесредовых факторов // Математическая биология и биоинформатика. - 2012. - № 2. -С. 589-610.
13. Физиология человека / под ред. В.М. Покровского, Г.Ф. Коротько. - М.: Медицина, 1997. -Т. 1. - 448 с.; Т. 2. - 368 с.
14. Ainsworth M.A., Fenger C., Svendsen P., Schaffalitzky De Muckadell O.B. Effect of stimulation of mucosal HCO-secretion on acid-induced injury to porcine duodenal mucosa // Scandinavian Journal of Gastroenterology. - 1993. - Vol. 28. - P. 1091-1097.
15. Boeykens K., Steeman E., Duysburgh I. Reliability of pH measurement and the auscultatory method to confirm the position of a nasogastric tube // International Journal of Nursing Studies. - 2014. - Vol. 51. -P. 1427-1433.
16. Chang J.H., Choi M.G., Yim D.S., Cho Y.K., Park J.M., Lee I.S., Kim S.W., Chung I.S. A novel placement method of the bravo wireless ph monitoring capsule for measuring intragastric pH // Dig. Dis. Sci. - 2009. -Vol. 54. - P. 578-585.
17. Clough M.R., Axon A.T. The сalcium сarbonate breath test, a noninvasive test of stimulated gastric acid secretion: preliminary communication // Eur. J. Gastroenterol. Hepatol. - 2009. - Vol. 21, № 3. -P. 266-272.
18. De Loubens S., Lentle R.G., Love R.J., Hulls C., Janssen P.W.M. Fluid mechanical consequences of pendular activity, segmentation and pyloric outflow in the proximal duodenum of the rat and the guinea pig // Journal of the Royal Society Interface. - 2013. - Vol. 10. - Article ID 20130027.
19. DiMagno E.P., Malagelada J.R., Go V.L. The relationships between pancreatic ductal obstruction and pancreatic secretion in man // Mayo Clinic Proceedings. - 1979. - Vol. 54, № 3. - P. 157-162.
20. Ehlen L.A., Marshall T.A., Qian F., Wefel J.S., Warren J.J. Acidic beverages increase the risk of in vitro tooth erosion // Nutr. Res. - 2008. - Vol. 28, № 5. - P. 299-303.
21. Farmer A.D., Scott S.M., Hobson A.R. Gastrointestinal motility revisited: The wireless motility capsule // United European Gastroenterol. J. - 2013. - Vol. 1, № 6. - P. 413-421.
22. Feil W., Wenzl E., Vattay P., Starlinger M., Sogukoglu T., Schiessel R. Repair of rabbit duodenal mucosa after acid injury in vivo and in vitro // Gastroenterology. - 1987. - Vol. 92. - P. 1973-1986.
23. Ferrua M.J., Xue Z., Singh R.P. On the kinematics and efficiency of advective mixing during gastric digestion - A numerical analysis // Journal of Biomechanics. - 2014. - Vol. 47. - P. 3664-3673.
24. Ferrua M.J., Singh R.P. Computational modelling of gastric digestion: current challenges and future directions // Current Opinion in Food Science. - 2015. - Vol. 4. - P. 116-123.
25. Fullard L.A., Lammers W.J., Ferrua M.J. Advective mixing due to longitudinal and segmental contractions in the ileum of the rabbit // Journal of Food Engineering. - 2015. - Vol. 160. - P. 1-10.
26. Georgia D., Stamatina T., Maria N., Konstantinos A., Konstantinos F., Emmanouil L., Georgios Z., Dimitrios T. 24-h multichannel intraluminal impedance pH-metry 1 year after laparocopic sleeve gastrectomy: an objective assessment of gastroesophageal reflux disease // Obesity Surgery. - 2017. -Vol. 27. - P. 749-753.
27. Ghosh T., Lewis D.I., Axon A.T.R., Everett S.M. Review article: methods of measuring gastric acid secretion // Aliment. Pharmacol. Ther. - 2011. - Vol. 33. - P. 768-781.
28. Hao S., Wang B., Wang Y. Density-dependent gastroretentive microparticles motion in human gastric emptying studied using computer simulation // European Journal of Pharmaceutical Sciences. - 2015. -Vol. 70. - P. 72-81.
29. Haruma K., Kamada T., Kawaguchi H., Okamoto S., Yoshihara M., Sumii K., Inoue M., Kishimoto S., Kajiyama G., Miyoshi A. Effect of age and Helicobacter pylori infection on gastric acid secretion // Journal of Gastroenterology and Hepatology. - 2000. - Vol. 15. - P. 277-283.
30. Haruma K., Mihara M., Okamoto E., Kusunoki H., Hananoki M., Tanaka S., Yoshihara M., Sumii K., Kajiyama G. Eradication of Helicobacter pylori increases gastric acidity in patients with atrophic gastritis of the corpus-evaluation of 24-h pH monitoring // Alimentary Pharmacology & Therapeutics. - 1999. -Vol. 13. - P. 155-162.
31. Hunter P., Robbins P., Noble D. The IUPS human Physiome Project // Pflugers Archiv - European Journal of Physiology. - 2002. - Vol. 445. - P. 1-9.
32. Isenberg J.I., Selling J.A., Hogan D.L., Koss M.A. Impaired proximal duodenal mucosal bicarbonate secretion in patients with duodenal ulcer // The New England Journal of Medicine. - 1987. - Vol. 316. -P. 374-379.
33. Ito S., Lacy E.R., Rutten M.J., Critchlow J., Silen W. Rapid repair of injured gastric mucosa // Scand. J. Gastroenterol. Suppl. - 1984. - Vol. 101. - P. 87-95.
34. Jeffrey B., Udaykumar H.S., Schulze K.S. Flow fields generated by peristaltic reflex in isolated guinea pig ileum: impact of contraction depth and shoulders // American Journal of Physiology - Gastrointestinal and Liver Physiology. - 2003. - Vol. 285. - P. G907-G918.
35. Johannessen E.A., Wang L., Reid S.W.J., Cumminga D.R.S, Cooper J.M. Implementation of radiotelemetry in a lab-in-a-pill format // Lab. Chip. - 2006. - Vol. 6. - P. 39-45.
36. Kadiyala V., Lee L.S., Banks P.A., Suleiman S., Paulo J.A., Wang W., Jessica Rosenblum J., Sainani N.I., Mortele K., Conwell D.L. Cigarette smoking impairs pancreatic duct cell bicarbonate secretion // Journal of the Pancreas. - 2013. - Vol. 14, № 1. - P. 31-38.
37. Kapicioglu S., Baki A., Tekelioglu Y., Arslan M., Sari M., Ovali E. The inhibiting effect of cola on gastric mucosal cell cycle proliferation in humans // Scand. J. Gastroenterol. - 1998. - Vol. 33, № 7. -P. 701-703.
38. Katelaris P.H., Seow F., Lin B.P., Napoli J., Ngu M.C., Jones D.B. Effect of age, Helicobacter pylori infection, and gastritis with atrophy on serum gastrin and gastric acid secretion in healthy men // Gut. -1993. - Vol. 34. - P. 1032-1037.
39. Kozu H., Kobayashi I., Nakajima M., Uemura K., Sato S., Ichikawa S. Analysis of flow phenomena in gastric contents induced by human gastric peristalsis using CFD // Food Biophysics. - 2010. - Vol. 5. -P. 330-336.
40. Koziolek M., Grimm M., Becker D., Iordanov V., Zou H., Shimizu J., Wanke C., Garbacz G., Weitschies W. Investigation of pH and temperature profiles in the GI tract of fasted human subjects using the intellicap® system // J. Pharm. Sci. - 2015. - Vol. 104. - P. 2855-2863.
41. Kuchumov A.G., Gilev V., Popov V., Samartsev V., Gavrilov V. Non-newtonian flow of pathological bile in the biliary system: Experimental investigation and CFD simulations // Korea Aust. Rheol. J. - 2014. -Vol. 26. - P. 81-90.
42. Kuchumov A.G., Nyashin Y.I., Samarcev V.A., Gavrilov V.A. Modelling of the pathological bile flow in the duct with a calculus // Acta Bioeng. Biomech. - 2013. - Vol. 15. - P. 9-17.
43. Kuchumov A, Tuktamyshev V., Kamaltdinov M. Peristaltic flow of lithogenic bile in the vateri's papilla as non-newtonian fluid in the finite-length tube: analytical and numerical results for reflux study and optimization // Lekar a Technika - Clinician and Technology. - 2017. - Vol. 47, № 2. - P. 35-42.
44. Lam S.K. Pathogenesis and pathophysiology of duodenal ulcer // Clinics in Gastroenterology. - 1984. -Vol. 13, № 2. - P. 447-472.
45. Li M., Brasseur J.G., Dodds W.J. Analysis of normal and abnormal esophageal transport using computer simulation // American Journal of Physiology - Gastrointestinal and Liver Physiology. - 1994. - Vol. 266. -P. 525-543.
46. Lim Y.F., de Loubens C., Love R.J., Lentle R.G., Janssen P.W.M. Flow and mixing by small intestine villi // Food and Function. - 2015. - Vol. 6. - P. 1787-1795.
47. Livingston E.H., Passaro E.P. Jr, Miller J., Guth P.H. Spectrum of injury produced in the duodenum by perfusion with luminal acid in the rat // Gastroenterology. - 1992. - Vol. 103. - P. 481-489.
48. Love R.J., Lentle R.G., Asvarujanon P., Hemar Y., Stafford K.J. An expanded finite element model of the intestinal mixing of digesta // Food Digestion. - 2013. - Vol. 4. - P. 26-35.
49. Mackay S., Jacobson B. Endoradiosonde // Nature. - 1957. - Vol. 179. - P. 1239-1240.
50. Maerten P., Ortner M., Michetti P., Dorta G. Wireless capsule pH monitoring: does it fulfil all expectations? // Digestion. - 2007. - Vol. 76. - P. 235-240.
51. McCloy R.F., Greenberg G.R., Baron J.H. Duodenal pH in health and duodenal ulcer disease: effect of a meal, Coca-Cola, smoking, and cimetidine // Gut. - 1984. - Vol. 25, № 4. - P. 386-392.
52. Pal A., Indireshkumar K., Schwizer W., Abrahamsson B., Fried M., Brasseur J.G. Gastric flow and mixing studied using computer simulation // Proceedings of the Royal Society B. - 2004. - Vol. 271. -P. 2587-2594.
53. Phan J., Benhammou J.N., Pisegna J.R. Gastric hypersecretory states: investigation and management // Curr. Treat. Options Gastro. - 2015. - Vol. 13. - P. 386-397.
54. Reddy A., Norris D.F., Momeni S.S., Waldo B., Ruby J.D. The pH of beverages in the United States // J. Am. Dent. Assoc. - 2016. - Vol. 147, № 4. - P. 255-263.
55. Saad, R.J. The wireless motility capsule: a one-stop shop for the evaluation of GI motility disorders // Curr. Gastroenterol. Rep. - 2016. - Vol. 18. - P. 14.
56. Sammon A.M., Ndebia E.J., Umapathy E., Iputo J.E. 24-Hour measurement of gastric pH in rural South Africa // Gastroenterology Research and Practice. - 2015. - Vol. 2015. - Article ID 658106. - 5 pp.
57. Shay H. Etiology and pathology of gastric and duodenal ulcer/ H. Shay, D.C.H. Sun // Bockus H.L. Gastroenterology. - Philadelphia-London, 1968. - Vol. 1. - P. 420-465.
58. Shinkai H., Iijima K., Koike T., Nakagawa K., Maejima R., Endo H., Ara N., Asano N., Imatani A., Ohara S., Shimosegawa T. Calcium carbonate breath test for non-invasive estimation of gastric acid secretion // Tohoku J. Exp. Med. - 2014. - Vol. 232, № 4. - P. 255-261.
59. Sinnott M.D., Cleary P.W., Harrison S.M. Peristaltic transport of a particulate suspension in the small intestine // Applied Mathematical Modelling. - 2017. - Vol. 44. - P. 143-159.
60. Solana M.J., López-Herce J., Sánchez C., Cordero M., Urbano J., Botrán M., Del Castillo J., Bellón J.M., Carrillo A. Comparison of continuous pH-meter and intermittent pH paper monitoring of gastric pH in critically ill children // Eur. J. Gastroenterol. Hepatol. - 2012. - Vol. 24, № 1. - P. 33-36.
61. Starlinger M., Jakesz R., Matthews J.B., Yoon C., Schiessel R. The relative importance of HCO3 - and blood flow in the protection of rat gastric mucosa during shock // Gastroenterology. - 1981. - Vol. 81, № 4. -P. 732-735.
62. Svanes K., Ito S., Takeuchi K., Silen W. Restitution of the surface epithelium of the in vitro frog gastric mucosa after damage with hyperosmolar sodium chloride // Gastroenterology.- 1982. - Vol. 82. -P. 1409-1426.
63. Taghipoor M., Barles G., Georgelin C., Licois J.R., Lescoat P. Digestion modeling in the small intestine: Impact of dietary fiber // Mathematical Biosciences. - 2014. - Vol. 258. - P. 101-112.
64. Trusov P.V., Zaitseva N. V., Kamaltdinov M.R. A multiphase flow in the antroduodenal portion of the gastrointestinal tract: a mathematical model // Computational and Mathematical Methods in Medicine. -2016. - Vol. 2016. - Article ID 5164029. - 18 pp.
65. Vellas B., Balas D., Moreau J., Bouisson M., Senegas-Balas F., Guidet M., Ribet A. Exocrine pancreatic secretion in the elderly // International Journal of Pancreatology. - 1988. - Vol. 3. - P. 497-502.
66. Wang Y.L., Chang C.C., Chi C.W., Chang H.H., Chiang Y.C., Chuang Y.C., Chang H.H., Huang G.F., Liao Y.S., Lin C.P. Erosive potential of soft drinks on human enamel: an in vitro study // J. Formos. Med. Assoc. - 2014. - Vol. 113, No. 11. - P. 850-856.
67. Weinstein D.H., deRijke S., Chow C.C., Foruraghi L., Zhao X., Wright E.C., Whatley M., Maass-Moreno R., Chen C.C., Wank S.A. A new method for determining gastric acid output using a wireless pH-sensing capsule // Alimentary Pharmacology and Therapeutics. - 2013. - Vol. 37. -P. 1198-1209.
68. Xue Z., Ferrua M.J., Singh R.P. Computational fluid dynamics modeling of granular flow in human stomach // Alimentos hoy. - 2012. - Vol. 21. - P. 3-14.
MULTI-COMPONENT MIXTURE FLOW IN THE STOMACH AND DUODENUM ALLOWING FOR FUNCTIONAL DISORDERS: RESULTS OF NUMERIC MODELLING FOR DETERMINING ACIDITY
M.R. Kamaltdinov, P.V. Trusov, N.V. Zaytseva (Perm, Russia)
This paper covers certain aspects of mathematic model development for describing mixture flow in the antroduodenum. It details data for numeric research on pH level distribution in the antroduodenal section of the gastrointestinal tract. One-phase multi-component mixture contains components which are dissolved at the molecular level:
hydrochloric acid, sodium hydrocarbonate, water, and a product of the reaction between acid and alkaline. The model allows for functional disorders in four antroduodenum sections which result in poorer secretory and motor functions performance. Mathematical tasking includes mass and momentum conservation equations for mixture components taking account mass sources due to components secretion and neutralization reaction. Kinematic boundary conditions are set on antroduodenum walls allowing for peristaltic motion (no slip condition). As for other boundaries of the examined area (inlet/outlet sections), here we set values of the 1st phase mass fractions, pressure, and conditions under which tangent components of stress vector are equal to zero. To determine pH level of the mixture via components concentrations, we obtained two algebraic equations, cubic and quatric ones. Model parameters were identified on the basis of data taken from literature and performed ultrasound research. We discuss various calculated scenarios where we examined impacts exerted by initially increased acidity (for example, due to acidic fluid intake) on dynamic acidity spread in the tract depending on fluid volume and pH level, occurence of functional disorders in alkaline secretion and in antrum motility. We showed that drinks with pH = 2.3 can potentially damage mucous coat of the tract walls, especially if alkaline secretion disorders occur. The approaches we suggest in the paper can in future be used for working out individual recommendations on consuming food and drinks with increased acidity.
Key words: numeric modelling, multi-component mixture, increased acidity, pH level, antroduodenum, functional disorders.
Получено 26 июля 2017
