Assessing hydrodynamic resistence coefficient within energetic model of the human upper respiratory tract

Full Text

Исследование движения воздушных потоков внутри полостей носа, представляющих сложную нерегулярную структуру каналов и пазух, открывает путь к пониманию сложных процессов в дыхательной системе человека. Внутренняя архитектоника носа определяет проходимость носовых каналов. Нарушения ее сечений являются проявлением ряда серьезных заболеваний, таких как риниты, в том числе аллергические, синуситы, полипы и др. [4]. Важную роль в определении проходимости верхних дыхательных путей ринология отводит оценке коэффициента носового сопротивления, объективно отражающего нарушения аэродинамики дыхательного цикла. Коэффициент носового сопротивления является аналогом коэффициента гидродинамического сопротивления крови, предложенного французским физиологом Пуазейлем для ламинарного течения крови в сосуде постоянного сечения, и определяется по величине перепада давления и расхода внутри трубки, подводимой к органам дыхания. Проведенные исследования, например [3; 6], показали, что движение воздуха внутри полостей носа и околоносовых пазух не является однозначно ламинарным и сопровождается наличием на некоторых участках турбулентных вихрей. В этой связи возникает необходимость привлечения новых, более чувствительных к турбулентному потоку моделей. Кроме того, важность представляет оценка коэффициента на различных участках внутреннего сечения носа, что, в связи с особенностями строения носа, возможно только при применении моделирования. Предлагаемая вниманию читателя вычислительная модель основана на применении закона сохранения энергии при распределении воздушных потоков внутри полостей носа человека. Необходимые в расчете коэффициента скорости воздушного потока на различных участках полостей носа, а именно - в области преддверия носа, хоанах, верхнем, среднем, нижнем и общем носовых ходах, были определены в процессе моделирования движения воздушных потоков внутри натурной модели носа [5]. Коэффициент носового сопротивления в вычислительной модели заменен на коэффициент гидродинамического сопротивления, являющийся безразмерной величиной и обладающий аналогичным физическим смыслом [2]. Внутренняя полость носа рассмотрена как система, состоящая из последовательного и параллельного соединения каналов (рис. 1). В энергетической модели принято: · в области преддверия носа, в некоторой точке, удаленной от входа в нос, центральная струйка тока делится на четыре части; · после некоторого неравновесного положения образовавшиеся струйки центрально распределяются внутри общего, верхнего, среднего и нижнего носовых ходов (в этих местах были расположены анемометры внутри натурной модели носа для определения скорости воздушного потока); · при дальнейшем движении воздуха внутри полости носа происходит объединение струек в одну, также центральную для области хоан; · та же картина перераспределения повторяется при обратном движении воздуха от области хоан к области преддверия носа. Рис. 1. Нос как система последовательно и параллельно соединенных каналов Индексы при расходе G соответствуют: 1 - общему носовому ходу, 2 - среднему носовому ходу, 3 - верхнему носовому ходу, 4 - нижнему носовому ходу. Индексы при площади поперечного сечения канала F соответствуют: 11 - в области преддверия носа, 21 - в общем носовом ходу, 22 - в среднем носовом ходу, 23 - в верхнем носовом ходу, 24 - в нижнем носовом ходу, 33 - в хоанах На рис. 1 величина G соответствует расходу вдыхаемого и выдыхаемого воздуха, величина F соответствует площадям сечения внутри каналов. Общее гидродинамическое сопротивление системы можно представить формулой: , (1) где ζ11 - коэффициент гидродинамического сопротивления области преддверия носа, ζ2 - суммарный коэффициент гидродинамического сопротивления областей среднего, нижнего, верхнего и общего носовых ходов, ζ33 - коэффициент гидродинамического сопротивления области хоан. Для двух сечений центрального положения струек тока выполняется закон сохранения энергии: (2) В этом уравнении масса относится к единице объема отдельно взятой струйки, и, так как плотность воздуха при его движении внутри объема не меняется, ее можно считать одинаковой в сечениях 1-1 и 2-2. Отношение N·t отражает работу легких и является произведением силы дыхательных движений легких («мощности легких») на время, затрачиваемое на переход из сечения 1-1 в сечение 2-2. ΔH выражает энергию, преобразующуюся в дополнительное тепло при диссипации. Произведения: , , , характеризуют кинетическую энергию, статическое давление, долю потенциальной энергии и внутреннюю тепловую энергию в i-том носовом ходе соответственно. Поскольку размеры модели невелики, доля потенциальной энергии внутри носовых ходов не меняется, поэтому в дальнейших вычислениях она не рассматривалась. Количество теплоты, подводимое к единице массы в каждом сечении, зависит от температуры нагретых стенок, которая постоянна для всей полости носа . Поэтому можно считать, что Eтепл1=Eтепл2 и в дальнейших вычислениях не рассматривать. Тогда уравнение Бернулли с учетом изменения скорости, температуры и давления внутри сечений: , (3) где ΔNобщ - общая мощность, теряемая на участках 1-1, 2-2 и характеризующая величину механической энергии, превращаемой в теплоту, Вт. В общем случае, коэффициент гидродинамического сопротивления может быть определен из соотношения: (4) Гидродинамическое сопротивление внутри модели можно рассчитать по формулам (3) и (4), приняв следующие допущения: 1) скорость воздушного потока максимальна в центре потока, следовательно, при условии расположения анемометра в этой области при моделировании дыхания максимальное значение поверхностной и временной скоростей совпадут, тогда средний поверхностный интеграл изменения скорости внутри рассматриваемого сечения можно заменить осреднением по времени: i. ; 2) долю кинетической энергий можно вынести за знак интеграла, т.к. плотность среды постоянна; 3) температуру воздуха внутри модели можно считать примерно одинаковой и равной температуре среды, тогда значения внутренней энергии в сечениях 1-1 и 2-2 можно считать одинаковыми и при расчете коэффициента гидравлического сопротивления не рассматривать; 4) статическое давление в сечении зависит от скорости воздушного потока. Его, как и скорость, можно принять как некоторую усредненную величину. (В экспериментальных исследованиях не рассматривалась). Введя среднюю скорость , где Т - промежуток времени, называемый периодом осреднения, уравнение Бернулли можно представить как следующее соотношение: , (5) где ρ - массовая плотность воздуха, кг·с2/м4; , - средние скорости воздушного потока в выбранных сечениях рассматриваемого и последующего за ним каналов, м/с. Коэффициент гидравлического сопротивления ζi: . (6) Коэффициент гидродинамического сопротивления в области преддверия носа, нижнем и верхнем ходах и в области хоан был оценен через отношение: , (7) где z - гидравлическое сопротивление, определяемое через расход . (8) Для нижнего и верхнего носовых ходов расход определялся по формуле: , (9) где G - расход в области преддверия носа, а также в хонах, G1 - в общем носовом ходу, G2 - в среднем носовом ходу, G3 - в верхнем носовом ходу, G4 - в нижнем. Общее гидравлическое сопротивление при рассмотрении носа как системы последовательно параллельно соединенных каналов (рис. 2): (10) Рис. 2. Блок-схема гидравлических сопротивлений внутри полостей модели носа Обобщенное сопротивление для нижнего и верхнего ходов ζнижн,верхн были найдены через площадь, равную сумме площадей этих каналов. Обобщенное сопротивление «преддверие - общий ход - ζпох», можно найти как сумму сопротивлений «преддверие - общий ход - ζпреддв-общий», «общий ход - хоаны - ζобщий-хоаны»: (11) Также можно найти обобщенное сопротивление «преддверие - средний ход», ζпсх: , (12) где ζпреддверие-средний - коэффициент гидравлического сопротивления области «преддверие - средний ход», ζсредний-хоаны - коэффициент гидравлического сопротивления области средний ход - хоаны. Результаты оценки гидродинамического сопротивления в энергетической модели верхних дыхательных путей представлены в таблице ниже. Таблица 1 Оценка коэффициента гидродинамического сопротивления внутри модели носа через отношение энергий ζ, вдох ζ, выдох Область преддверия носа 0,27 0,28 Область преддверия - общий ход 0,75 0,29 Область преддверия - средний ход 0,83 0,28 Общий ход - хоаны 0,07 0,94 Средний ход - хоаны 0,15 0,80 Верхний и нижний ходы 0,04 0,04 Хоаны 0,20 0,20 Общее гидродинамическое сопротивление 0,49 0,50 Сравнение коэффициентов гидродинамического сопротивления показало: 1. При вдохе воздушный поток испытывает минимальное сопротивление в области верхнего и нижнего носовых ходов, поэтому основной поток вдыхаемого воздуха должен устремиться именно в эти каналы. 2. Величины гидродинамического сопротивления при переходе из общего и среднего носовых ходов в хоаны оказались меньше, чем при переходе воздушного потока из области преддверия носа в эти носовые ходы, поэтому при выдохе движение воздушного потока по носовым ходам не будет испытывать такого сопротивления, как при вдохе. Возможно, что такое движение воздушного потока способствует обогащению пазух носа, которые в основном выходят в средний носовой ход, теплым выдыхаемым воздухом. 3. Величина общего гидродинамического сопротивления при вдохе оказалась равной величине общего гидродинамического сопротивления при выдохе. Дополнительно к проведенному исследованию добавим, что нами в работе [1] показано наличие ярко выраженной вихревой природы для потока вдыхаемого воздуха. Отмечено, что наиболее сложные траектории движения получены для нижнего носового хода. Минимальное значение гидродинамического сопротивления в области нижнего носового хода при наличии вихрей [1] способствует увеличению контакта вдыхаемого потока воздуха со стенками носа, что обеспечивает согревание воздуха в носу.

About the authors

Gennadij Nikolaevich Lukyanov

St. Petersburg National Research University of Information Technologies, Mechanics and Optics

Head of Electronics Department

Anna Aleksandrovna Rassadina

St. Petersburg National Research University of Information Technologies, Mechanics and Optics

Assistant Professor at the Electronics Department

References

Supplementary files