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

Abstract


The article describes a model of energy transfer in human upper respiratory tract, for which we estimate the hydrodynamic drag coefficient. The internal cavity of the nose represents a system of series- and parallel- connected channels. While modeling, we assumed the following points: (1) separation of central flow of inhaled air (hereinafter referred to as a stream filament) into four streams within nasal vestibule; (2) central distribution of stream filaments into communis, superior, medius, and inferior nasal meatus; (3) merging of the stream filaments into one in the choanae; (4) repetition of the pattern of redistribution of stream filaments under the reverse movement of air from choanae to the nasal vestibule. During the modeling of air flows inside the full-scale model of the nose the air flow velocity rate, which is necessary to calculate hydrodynamic drag coefficient, was determined for different parts of the nasal cavities, namely - nasal vestibule, choanae, superior, medius, inferior and communis nasal meatus. Full-scale model of the nose sinus was made of epoxy. For measuring velocity of air flow miniature bead heat-loss anemometers were placed inside the model: in the nasal vestibule, communis, superior, medius and inferior nasal meatus and choanae. The model’s energy transfer is based on the law of conservation for steady incompressible fluid flows (Bernoulli's equation) and Kirchhoff's rules. Solution of the Bernoulli's equation is founded on the following assumptions: (1) mass refers to the volume of a separate stream filament, and, since the density of the air as it moves inside the volume does not change, it can be considered identical in sections 1-1 (nasal vestibule) and 2-2 (choanae); (2) the share of the potential energy was not considered inside the nasal meatus due to the small size of the model; (3) the quantity of heat supplied to a unit mass in each section depends on the temperature of the heated walls, which is constant in the entire nasal cavity, therefore it is assumed that E heat 1= E heat 2 and is not considered in further calculations; (4) air flow velocity is maximum in the center of the flow, hence, due to the location of anemometer in this area, the maximum surface velocity and time velocity would coincide, causing the average surface integral of velocity change within the section can be replaced by time averaging; (5) because the density of the medium is constant, kinetic energy can be taken outside the integral sign; (6) air temperature inside the model can be considered approximately the same and equal to the temperature of the medium, then the values of the internal energy in sections 1-1 and 2-2 can be regarded as the same and not considered in the calculation of the hydrodynamic resistance coefficient; (7) the static pressure in the cross section depends on the air flow velocity, therefore, it as the velocity can be taken as an average value. The simulation results showed that: (1) during an intake of breath, hydrodynamic resistance coefficient was the lowest for the superior and inferior nasal meatus and was maximal in the communis and medius nasal meatus; (2) during exhalation the value of hydrodynamic resistance coefficient decreases because of the air movement in the communis and medius nasal meatus.

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

  1. Воронин А.А., Лукьянов Г.Н., Неронов Р.В. Моделирование воздушного потока в каналах нерегулярной формы // Научно-технический вестник информационных технологий, механики и оптики. 2013. № 3.
  2. Идельчик И.Е. Справочник по гидравлическим сопротивлениям. М., 1992.
  3. Нечипоренко А.С., Гарюк О.Г., Чмовж В.В. Критерий идентификации фаз дыхательного цикла // Вестник Национального технического университета Харьковский политехнический институт. Серия «Информатика и моделирование». 2013. № 19 (992).
  4. Пискунов Г.З. Причины роста распространенности заболеваний носа и околоносовых пазух // Российская ринология. 2009. № 1.
  5. Рассадина А.А., Лукьянов Г.Н., Усачев В.И. Натурная модель дыхательной системы человека // Вестник седьмой всероссийской научной конференции молодых ученых и аспирантов «Новые информационные технологии. Разработка и аспекты применения». Таганрог, 2004.
  6. Lukyanov G., Rassadina A., Usachev V. Comparison and the analysis of the processes of the movement of air through the human breathing system and its natural model // International Conference on Physics and Control, PhysCon 2005. St. Petersburg, 2005.

Statistics

Views

Abstract - 22

PDF (Russian) - 71

Cited-By


Article Metrics

Metrics Loading ...

Refbacks

  • There are currently no refbacks.




This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies