Численные исследования течения воздуха в дыхательном тракте человека Текст научной статьи по специальности «Физика»

Мезо-, нано-, биомеханика и механика природных процессов Вестник Нижегородского университета им. Н.И. Лобачевского, 2011, № 4 (2), с. 410-412

УДК 533.6.011:51+612.215.4

ЧИСЛЕННЫЕ ИССЛЕДОВАНИЯ ТЕЧЕНИЯ ВОЗДУХА В ДЫХАТЕЛЬНОМ ТРАКТЕ ЧЕЛОВЕКА

© 2011 г. В.Л. Ганимедов, М.И. Мучная, А. С. Садовский, В.Н. Шепеленко

Институт теоретической и прикладной механики им. С. А. Христиановича СО РАН, Новосибирск

ganim@itam.nsc.ru

Поступила в редакцию 16.05.2011

Построение геометрической модели дыхательного тракта и расчеты в рамках уравнений Навье - Сто -кса проводились с помощью программного пакета «Fluent» и его геометрического препроцессора «Gambit». Основные результаты расчетов следующие. Получено, что возвратное течение, возникающее в нижней части широкой носовой полости, не препятствуя дыхательной функции, играет позитивную роль, заключающуюся в прогреве вдыхаемого воздуха. Показано, что течение в трахее при вдохе представляет собой закрученный поток с образованием нескольких вихрей. При выдохе на структуру и параметры течения воздуха в трахее основное влияние оказывает геометрия бронхиального дерева и, в меньшей степени, строение носовой полости. Численная имитация клинического метода риноманометрии показала, что его результаты не отражают реальной зависимости «перепад давления - объемный расход» при дыхании в естественных условиях.

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

Геометрия носовой полости человека моделировалась по серии томограмм, выполненных в параллельных коронарных сечениях, полученных на рентгеновском томографе с шагом 2-^5 мм. Построение геометрии модели и расчеты в рамках уравнений Навье - Стокса проводились с помощью программного пакета «Fluent» и его геометрического препроцессора

«Gambit» (Лицензия ID:KMMG-SK-LO1). На рис. 1 приведено изображение поверхности одной из рассматриваемых носовых полостей.

К сожалению, в представленных медиками томограммах не Рис. 1 имелось информации о

носоглотке и трахее пациентов. Поэтому после носовой полости поверхность дыхательного тракта аппроксимировалась элементами стандартных фигур.

Перечислим основные результаты проведенных исследований.

1) При широком основном носовом канале слабое возвратное течение, возникающее в нижней части носовой полости, не препятствуя дыхательной функции, играет позитивную роль,

заключающуюся в прогреве вдыхаемого воздуха. Особенно это должно быть важно при низких наружных температурах. На рис. 2 приведены поля температуры (а) и продольной ком -поненты скорости (б) в коронарных сечениях.

б) Рис. 2

2) Создана математическая модель используемой в клинических исследованиях процедуры передней активной риноманометрии (ПАР) (с помощью аппарата фирмы НОМОТН), регистрирующего значение объемного расхода воздуха через нос в зависимости от респираторного усилия. Проведена серия параметрических расчетов. Получено, что течение в носовых полостях при выполнении процедуры ПАР имеет струйный характер и зависит от того, под каким углом установлен оператором носовой адаптер. На рис. 3 приведена картина визуализации течения по продольной компоненте скорости У2 в правой носовой полости при угле установки носового адаптера, равном 22°, относительно продольной оси.

Рис. 3

Анализ результатов приводит к выводу, что показания ПАР дают только косвенную информацию о возможностях транспортной функции носа, и риноманограммы не отражают реальной зависимости «перепад давления - объемный расход» при дыхании в естественных условиях.

3) Численные исследования показали, что течение в трахее при вдохе в моделях, построенных для пяти различных геометрий носовой полости конкретных людей, представляет собой закрученный поток с образованием нескольких вихрей. На рис. 4 приведен пример построенной геометрической модели дыхательных путей, на рис. 5 - распределение векторов скорости при вдохе в двух характерных сечениях этой модели: Б-Б (рис. 5а) и В-В (рис. 5б).

а) б)

Рис. 5

Количество и интенсивность вихрей полностью определяется профилем скорости в сечении, где соединяются в одну правая и левая носовые полости, этот профиль зависит от реальной геометрии обоих носовых полостей. При выдохе на структуру и параметры течения воздуха в трахее основное влияние оказывает геометрия бронхиального дерева и, в меньшей степени, строение носовой полости.

Список литературы

1. Фомин В.М. и др. Исследование течения воздуха в носовой полости человека // ПМТФ. 2010. Т. 51, №2 С. 107-115.

NUMERICAL STUDY OF AIR FLOW IN THE HUMAN RESPIRATORY TRACT V.L. Ganimedov, M.I. Muchnaya, A.S. Sadovsky, V.N. Shepelenko

Construction of a geometric model and the calculations in the framework of the Navier - Stokes equations were performed using the software package «Fluent» and its geometrical preprocessor «Gambit». The main results of the calculations are as follows. It was found that the reverse flow that occurs in the lower part of the broad nasal cavity does not prevent respiratory function and plays a positive role, consisting in heating of inhaled air. It is shown that the swirling flow is generated in the trachea during the inspiration with the formation of several vortices. During the expiration the air flow structure and parameters in the trachea are influenced mainly by the geometry of the bronchial tree and, to a lesser extent, by the nasal cavity configuration. The numerical simulation of a clinic rhinomanometry method has shown that it does not reflect the real "the pressure drop - volume flow rate" relationship for breathing in natural conditions.

Keywords: numerical simulation, nasal cavity, coronal sections, trachea, swirling flow, volume flow rate.