Критическое давление потери устойчивости в складных трубках, актуальное для биомедицинских потоков
ДомДом > Новости > Критическое давление потери устойчивости в складных трубках, актуальное для биомедицинских потоков

Критическое давление потери устойчивости в складных трубках, актуальное для биомедицинских потоков

Jun 12, 2023

Научные отчеты, том 13, Номер статьи: 9298 (2023) Цитировать эту статью

77 Доступов

Подробности о метриках

Поведение спавшихся или стенозированных сосудов в организме человека можно изучать с помощью упрощенных геометрических фигур, таких как складная трубка. Целью данной работы является определение величины критического давления выпучивания разборной трубы с использованием теории фазового перехода Ландау. Методика основана на реализации экспериментально подтвержденной численной 3D модели разборной трубы. Критическое давление потери устойчивости оценивается для различных значений геометрических параметров системы путем рассмотрения связи между внутристеночным давлением и площадью центрального сечения как функции параметра порядка системы. Результаты показывают зависимость критических давлений потери устойчивости от геометрических параметров разборной трубы. Выведены общие безразмерные уравнения для критических давлений потери устойчивости. Преимущество этого метода в том, что он не требует каких-либо геометрических предположений, а основан исключительно на наблюдении, что выпучивание складной трубы можно рассматривать как фазовый переход второго рода. Исследованные геометрические и упругие параметры пригодны для биомедицинского применения, при этом особый интерес представляет изучение бронхиального дерева при патофизиологических состояниях, таких как астма.

Возможность изучения транспорта масс в организме человека, будь то воздух или кровь, с помощью математических и числовых моделей представляет собой один из наиболее плодотворных примеров соединения медицины и техники. Применение моделей вычислительной гидродинамики (CFD), взаимодействия жидкости со структурой (FSI) и аэроакустики значительно улучшило понимание патофизиологических состояний системы кровообращения1, дыхательной системы2, 3, процесса голосового воспроизведения4 и цереброваскулярной системы5, среди другие. Достоверность результатов, полученных с помощью таких численных моделей, должна быть подтверждена экспериментальными кампаниями в каждом конкретном случае. Разнообразие и геометрическая сложность сосудов человека могут сделать этот важный шаг чрезвычайно трудным. В связи с этим упрощенные модели, такие как складные трубки6,7,8, все еще широко используются как в числовых моделях, так и в клинической практике. Несмотря на упрощенную геометрию, феноменология складной трубки достаточно богата, чтобы охватить наиболее важные физические механизмы коллапса сосудов9. Динамика складной трубки существенно зависит от так называемого внутримурального давления, которое определяется как разница давлений между внутренней частью (просветом) и внешней частью трубки. При наличии потока жидкости необходимо учитывать дополнительный вклад, обусловленный ускорением потока вблизи сужения, что приводит к образованию области отрицательного статического давления. По мере увеличения внешнего давления (т. е. внутристеночного давления становится отрицательным) трубка начинает разрушаться. При критическом значении внутристеночного давления трубка испытывает явление коробления, в результате чего поперечное сечение становится двухлепестковым (см. рис. 1). Такая величина называется критическим давлением деформации и играет важную роль в оценке и диагностике многих патологий, связанных со стенозом и сужением10,11,12. В этой конфигурации небольшие изменения внутримурального давления могут привести к большим изменениям площади просвета. Если внешнее давление продолжает увеличиваться (или внутреннее давление продолжает уменьшаться из-за ускорения потока), внутренние стенки трубки соприкоснутся друг с другом (см. рис. 1), что в конечном итоге приведет к полному закрытию просвета. Точная и индивидуальная оценка критического давления коробления позволяет принимать более обоснованные клинические решения. Примером может служить критическое давление выпучивания глотки у пациентов, страдающих обструктивным апноэ во сне (СОАС). СОАС является наиболее распространенной патологией в спектре нарушений дыхания во сне13. Пациенты, страдающие СОАС, испытывают повторяющийся коллапс глотки во время сна, вызывая апноэ, которое серьезно влияет на качество жизни пациентов. Оценка тяжести патологии и выбор лечения сильно зависят от значений критического давления выпучивания глотки14. Однако его оценка требует, чтобы пациенты ночевали в больнице и находились под постоянным наблюдением, что приводит к весьма навязчивому опыту для пациента и к высокому экономическому воздействию на систему здравоохранения15.

0\) is the maximum value of the external pressure and \(t\in [0,\tau ]\). A sensitivity study in terms of the ratio \(p_{max}/\tau\) is discussed in "Sensitivity analysis" . At each time step, the value of \(p_{intr} = -p_{ext}\) is recorded. The central cross-section of the tube will be the most collapsed since it is the furthest radial cross-section from the constrained faces. To determine the value of the area, the radial coordinates \((r_i^j, \vartheta _i^j)\) of the corresponding deformed perimeter are recorded at each time-step, where the index \(i=1,\dots ,N\) labels the mesh elements and N is the total number of mesh elements on the perimeter. The index \(j=1,\dots ,M\) corresponds to the j-th time step and M is the total number of time steps needed for the external pressure to reach the value \(p_{max}\) in Eq. (3). The area \(A^j\) of the central cross-section can be then computed as:/p>0\) and \(c_2 >0\) are constants, \(\xi\) is a thermodynamic variable, and \(\xi _{crit}\) is the critical value for \(\xi\) at which the transition occurs. The minimization of the Landau potential yields an equation for the order parameter:/p>0\) and \(c_2>0\) are two free parameters, \(\tilde{p}=-p_{intr}\), and \(\beta >0\) is an additional free parameter, usually called critical exponent. The results of this analysis have been obtained according to the following procedure: /p>