Full Text

Реология в воздушно-капельной передаче инфекций Воздушно-капельный путь передачи инфекции является основным не только для актуального в настоящий момент коронавируса, но и для многих других инфекционных заболеваний, таких как ангина, дифтерия, корь, краснуха, ветряная оспа и др. При этих заболеваниях возбудители попадают в воздух с капельками слюны или лёгочной слизи при кашле, чихании, разговоре. Ротовая жидкость (или смешанная слюна) состоит из секретов слюнных желез (основной компонент) и множества микродобавок, присутствующих в ротовой полости и дыхательной системе. Воздушно-капельный путь передачи инфекции состоит из двух стадий. Сначала при дыхании, кашле, чихании, человека в окружающее воздушное пространство поступает воздух изо рта и легких человека. При движении воздуха через ротовую полость происходит захват и вовлечение в движение находящейся здесь жидкости. В результате выходящий воздух содержит множество мельчайших капель ротовой жидкости, которые в случае болезни человека могут содержать патогенные агенты. На второй стадии другой человек, вдыхая воздух с инфицированными каплями, может принести инфекцию в свой организм и подвергнуться заражению. Для препятствия поступления капель от больного человека в окружающую атмосферу, а также для защиты от проникновения капель в организм здорового человека используются различные фильтры и маски, задача которых задерживать движение капель. При чихании выделяется большое количество (до 50 000) капель слизи и слюны. Скорость капель может быть очень высокой, до 40 м/с (разные авторы дают значения 4,5-100 м/с), капли разлетаются на расстояние до 4-6 м. Судьбу выброшенных капель можно примерно описать следующим образом: капли диаметром более 5 мкм быстро оседают под действием силы тяжести, а капли менее 5 мкм могут длительное время оставаться в воздухе и переноситься воздушными потоками. Когда относительно крупная капля жидкости сталкивается с твёрдой поверхностью или попадает на материал маски или фильтра, она может разрушиться на несколько более мелких капель, создавая таким образом новый источник опасной инфекции в виде нового массива капель. Как ранее установлено, формирование всплеска и разрушение капли жидкости при столкновении с препятствием зависит от её реологических (объёмных и поверхностных) свойств [18]. В случае столкновения капель воды с препятствием при высоких значениях ударных чисел Рейнольдса и Вебера (Rei = rvidi/m, Wei = rvi2di/g, r - плотность жидкости; vi - скорость капли перед соударением; di - диаметр капли; m - вязкость жидкости; g - поверхностное натяжение жидкости) процесс контролируется исключительно инерцией и поверхностным натяжением [18, 35]. Эксперименты показали, что капля при соударении сохраняет сплошность при низких числах Вебера и разрушается при высоких числах Вебера. Переход от соударения с сохранением сплошности к соударению с разрушением капли (разбрызгиванию) происходит в диапазоне чисел Вебера Wei*Î(137, 206) (рис. 1, а, б) [32]. Исследования [10, 22, 26, 29, 30, 38] выявили упругость ротовой жидкости, притом, что её сдвиговая вязкость всего в 1-4 раза превосходит вязкость воды [9, 14]. Об упругости ротовой жидкости свидетельствует способность при растяжении формировать устойчивые утончающиеся нити. Время жизни таких нитей измеряется секундами, в то время как нить, сформированная водой, распадается за доли миллисекунды [3]. Опыты с ударным разрушением упругих растворов высокополимеров показали, что упругость радикально меняет устойчивость капель при соударении. Капли становятся значительно более прочными [33]. При столкновении с краёв всплеска срываются вторичные капли, которые, однако, при достаточной упругости жидкости не отрываются от всплеска, а оказываются связанными с ним тонкими утончающимися нитями. Под действием упругости нитей вторичные капли возвращаются назад, сливаясь с жидкостью исходной капли, при этом никакие свободные вторичные капли не формируются (рис. 1, в). Капля сохраняет сплошность. Более того, при достаточной упругости край всплеска может вообще не разрушаться, а виде жидкого тора сначала расширяться, а затем схлопываться, сливаясь в единую каплю (рис. 1, г). Исследования [33] показали, что определяющими параметрами, описывающими переход к разрушению упруговязкой капли при столкновении, наряду с уже упомянутым числом Вебера Wei = rvi2di/g, являются число Деборы Dei = θvi/di и упруго-капиллярное число Gdi/γ, где θ - время релаксации; G - модуль упругости. а б в г Рис. 1. Наблюдения сверху за разрушением капель жидкостей при падении на дискообразное препятствие диаметром 4 мм (выглядит как черный круг в центре кадра) [18, 32]. Прошло 6 мс после начала соударения: а, б - всплеск капли воды диаметром 2,8 мм при падении с высоты 18 и 27 см (Wei = 137 и 206) соответственно; в, г - аналогичные всплески капель растворов высокомолекулярных полиэтиленоксида и полиакриламида концентрацией 100 млн-1. Высота падения капли 65 см Реологические параметры θ и G - константы реологических уравнений состояния (РУС) Олдройда (вариант Oldroyd-B) или Максвелла (вариант Upper-Convected Maxwell UCM) [1, 24, 28]. Оба уравнения оказываются эквивалентными для описания быстрых деформаций растяжений при многократной вытяжке (см. раздел Теоретическая модель), как это имеет место при разрушении капель и утончении капиллярных нитей. Другим аспектом влияния упругости жидкости на передачу инфекций является способность упругости влиять на относительно медленное движение жидкости через каналы сложной формы или пористую среду [11, 24], что имеет место при защите от инфицированных капель при помощи различного рода фильтров. Те же параметры РУС Олдройда (вариант Oldroyd-B) и Максвелла (вариант Upper-Convected Maxwell UCM), а именно θ и G, управляют упругими эффектами в данном классе течений. Настоящая работа была предпринята, чтобы ответить на вопрос: какие конкретные значения реологических параметров ротовой жидкости θ и G можно использовать для моделирования столкновения капли ротовой жидкости с тем или иным препятствием или при движении через сложные среды. До сих пор реологические тесты ротовой жидкости проводились только для определения лишь одного параметра, а именно, времени релаксации θ [4, 27, 36, 38]. Вместе с тем для того, чтобы иметь полное реологическое описание свойств ротовой жидкости, необходимо установить значения модуля упругости G. Пара θ и G позволяет моделировать не только распад самоутончающихся нитей, но и другие явления, в которых участвует ротовая жидкость, такие как разрушения капель при столкновении с препятствиями и движение в фильтрах и каналах сложной формы. Экспериментальные наблюдения Данный раздел посвящён методу реологического испытания ротовой жидкости при помощи утончающейся капиллярной нити данной жидкости [23]. Проведены эксперименты по наблюдению за распадом капель ротовой жидкости in situ. Ротовая жидкость - водный раствор различных органических и минеральных компонентов, суммарная концентрация которых ca. 0,5% [16]. Среди основных компонентов - смесь секретов слюнных желез. Среди желез выделяют большие: околоушную, подчелюстную, подъязычную и ряд малых слюнных желез (размером 1-5 мм) - щечных, молярных, губных, язычных твердого и мягкого неба. Сами названия свидетельствуют о разнообразии источников слюнных секретов и, следовательно, о возможной неоднородности формируемой смеси. Плотность ротовой жидкости практически не отличается от плотности воды r = 1001…1017 кг/м3. Ротовая жидкость отбиралась при помощи пинцета непосредственно с губ испытуемого. На современном этапе автор проводил опыты на себе (мужчина 65 лет). Использовалась простейшая экспериментальная процедура. Сразу после отбора проба ротовой жидкости растягивалась между ножками пинцета. Далее сформированная жидкая нить утончалась сама по себе без какого-либо внешнего вмешательства. Процесс записывался при помощи видеокамеры смартфона iPhone 6 в режиме 30 кадров в секунду. Возможный режим 240 кадров в секунду не выявил преимущества и в экспериментах использовался редко. Использовалось разрешение кадров Full HD 1920´1080. Для получения размерных зависимостей в поле записи помещалась прозрачная линейка. Фотография использованного оборудования представлена на рис. 2. Следует отметить, что использование смартфона для визуализации утончения капиллярных нитей полимерных жидкостей, по-видимому, впервые продемонстрировано в [30], что показало возможность существенного удешевления данной экспериментальной методики. Рис. 2. Экспериментальное оборудование для исследования реологии упругих жидкостей: смартфон, пинцет и линейка Типичная видеозапись представлена на рис. 3. Видеоданные анализировались с помощью приложения QuickTime Player, хотя возможно использование и других приложений, таких как, например, ImageJ или Phantom. 3/0 4/0,0333 5/0,0667 6/0,100 8/0,167 10/0,233 15/0,400 30/0,900 Рис. 3. Видеозапись формирования и утончения нити ротовой жидкости (опыт 7114, также см. рис. 4, 5). Время возрастает слева направо и сверху вниз. На врезках указаны номер кадра/время утончения нити (секунды) Нить утончается во времени. На некоторой стадии возможен вторичный распад, когда появляется серия вторичных капелек. Измерения диаметра нити на видеокадрах позволили построить зависимость утончения нити ротовой жидкости во времени d = d(t). Появление вторичных капелек на поздних стадиях распада при измерениях игнорировалось. Начало отсчёта времени было выбрано в момент раздвижения ножек пинцета (кадр 3). Масштаб изображений определялся по прозрачной линейке. На рис. 4 приведена типичная зависимость. Данные свидетельствуют, что сразу после начала растяжения капля жидкости теряет устойчивость и быстро утончается под действием поверхностного натяжения от начального диаметра db до некоторого промежуточного диаметра ds в некоторый момент времени ts. Эта фаза деформирования соответствует растяжению капли обычной ньютоновской жидкости [5]. На этой стадии растяжения силы поверхностного натяжения уравновешиваются силами инерции жидкости и внутренними напряжениями. Внутренние напряжения сначала имеют вязкую природу, но по мере утончения проявляются упругие напряжения, вызванные быстрым деформированием жидкости. Так ведут себя растворы высокополимеров, в которых растворитель играет роль вязкой жидкости, а упругость формируется гигантским растяжением упругих гибкоцепных макромолекул [5]. В определённый момент (качественно обозначаемый ts) упругость начинает доминировать над всеми прочими силами и процесс распада нити резко замедляется. Дальнейшее утончение нити приобретает экспоненциальный характер, как об этом свидетельствуют данные на рис. 4. t, c 1,0 0,5 0 ts 1 0,1 d, мм ds db Рис. 4. Утончение нити ротовой жидкости - зависимость диаметра нити от времени. Сплошная кривая - аппроксимация экспериментальных точек (опыт 7114) аналитической зависимостью (6) на участке t > 0 Теоретическая модель Данный раздел представляет методику извлечения реологических характеристик ротовой жидкости из наблюдений за распадом нити этой жидкости. Исследования реологических свойств ротовой жидкости при растяжении с помощью метода утончающейся жидкости [4] показали, что реологическое поведение ротовой жидкости при интенсивном растяжении подобно поведению растворов полимеров [5]. Для описания реологии растворов полимеров при интенсивном растяжении апробированы РУС Олдройда (вариант Oldroyd-B) [24, 28] и Максвелла (вариант Upper-Convected Maxwell UCM) [1, 24, 28]. РУС Олдройда (вариант Oldroyd-B) описывается уравнениями s = -pI + 2msE + τ, τ = G(A - 1); dA/dt = A×Ñv + ÑvT×A - (A - I)/q, где s - тензор напряжений; p - давление; ms - вязкость растворителя; G - модуль упругости; E - тензор скоростей деформации; A - тензор упругих деформаций; I - единичный тензор; t - время; Ñv - градиент скорости; q - время релаксации. Рассматривая течение растяжения в направлении 1 с накоплением больших упругих деформаций (A1 >> 1), получим зависимость напряжений от кинематики растяжения A1 = l2exp(-t/q), t1 = GA1, (1) где l - степень удлинения образца в направлении 1, l >> 1. РУС Максвелла (вариант Upper-Convected Maxwell UCM) представлена уравнениями s = -pI + τ; dτ/dt = τ·E + τ·E - τ/θ + 2μE/θ, где s - тензор напряжений; τ - тензор избыточных напряжений; p - давление; I - единичный тензор; μ - элонгационная вязкость; E - тензор скоростей деформаций; q - время релаксации; t - время. В случае больших упругих деформаций (l2 >> 1, t1 >> μ/θ) t1 = (μ/θ) l2exp(-t/q), (2) где для упруговязких полимерных жидкостей отношение μ/θ представляет собой модуль упругости G [21]. Таким образом, в случае доминирования упругих напряжений осевая компонента тензора упругих напряжений (1) и осевая компонента тензора избыточных напряжений (2) совпадают, что свидетельствует об эквивалентности РУС Олдройда и Максвелла для описания течений упругих жидкостей с большими деформациями. Механическим аналогом данного РУС является последовательное соединение вязкого и упругого элементов [12] - рис. 5, а. θ G 2γ/d d λ τ1 а б Рис. 5. Модели упруговязкой реологии (а) и утончающейся нити (б) Как показали прямые экспериментальные измерения величины упругого натяжения нити раствора высокополимера на стадии экспоненциального утончения (t > ts) [2, 6], упругая компонента t1 приблизительно совпадает с величиной капиллярного давления в нити t1 » 2g/d (рис. 5, б), где g - поверхностное натяжение. Полагая для степени удлинения l = (db/d)2, из (1) (или (2)) получим зависимость утончения нити на экспоненциальной стадии 2g/d = G(db/d)4exp(-t/q), которую можно представить в виде d = ds exp(-(t-ts)/3q), (3) где величины уравнения (3) связаны между собой условием (2g/G)(ds3/db4)exp(ts/q) = 1. (4) Начальная быстрая стадия утончения нити t < ts моделировалась линейной зависимостью c коэффициентом пропорциональности α d = db - αt, (5) которая качественно соответствует утончению капиллярной нити вязкой жидкости [3-5, 23, 37]. В настоящей работе формула (5) использовалась лишь для кинематических оценок. Реологические параметры жидкости определялись путём аппроксимации экспериментальной зависимости d = d(t) (см., например, рис. 4) аналитической зависимостью, построенной на основе суперпозиции теоретических кривых (3) и (5). d = (db - αt)H(ts - t) + ds exp(-(t - ts)/3q)H(t - ts), (6) где H - функция Хевисайда; ds = db - αts в силу непрерывности теоретической зависимости d = d(t). Данная кривая непрерывна и состоит из двух ветвей, первая из которых описывает быструю начальную стадию утончения нити (5) при t < ts, а вторая меленную экспоненциальную (3) при t > ts. Сшивка решений в точке ts, ds определяет границу перехода одного режима течения в другой. Аппроксимация осуществлялась с помощью приложения Origin 6.1 (Раздел Analysis, NonLinear Curve Fitting). Определению подлежали параметры α, ts, q. В вычислениях аппроксимирующая функция задавалась как (для опыта 7114) (2.0 - ALF·x)/(1 + exp(-2·100·(-x + ts))) + + (2.0 - ALF·ts)·exp(-(x - ts)/3/TETA)/(1 + exp(-2·100·(x - ts))), в которой 2.0 = db - измеренный в данном опыте параметр db (см. рис. 3 и 4); ALF, ts, TETA - искомые параметры α, ts, q соответственно; x - независимая переменная, отвечающая времени t; 1/(1 + exp(-2·100·(-x + ts))) и 1/(1 + exp(-2·100·(x - ts))) гладкое моделирование функций Хевисайда H(ts - t) и H(t - ts). В результате численных вычислений находились искомые неизвестные параметры и строилась аппроксимирующая кривая (см. пример на рис. 4). Обработаны видеонаблюдения 8 опытов, результаты представлены на графике рис. 6 и в таблице. Время релаксации q определялось непосредственно в результате аппроксимации экспериментальных точек, а модуль упругости G рассчитывался по формуле (4). В вычислениях поверхностное натяжение ротовой жидкости полагалось g = 0,050 Н/м [9, 15]. Измеренные реологические данные (см. таблицу) характеризуются чрезвычайно высоким разбросом, хотя беглый взгляд на рис. 5 не дает оснований априори предполагать такую особенность. Две причины способствуют разбросу выходных данных, особенно заметному для модуля упругости. Во-первых, как ранее отмечалось, ротовая жидкость формируется путём смешения секретов нескольких слюнных желез, каждая из которых может формировать секрет с уникальными свойствами. Другая причина - сильная зависимость от входных данных. Так, модуль упругости определяется путём возведения в третью степень промежуточного диаметра ds и в четвёртую степень начального диаметра db. Поэтому небольшая ошибка в значениях ds, db многократно усиливается в процессе возведения в степени. Такая ситуация близка к случаю некорректно поставленной задачи [8]. С учётом приведенных обстоятельств целесообразно в качестве усредненных данных использовать медианы измеренных значений [13]. «Оценки медианы более робастны, её оценивание может быть более предпочтительным для распределений с так называемыми тяжёлыми хвостами. Главное преимущество медианы в устойчивости к выбросам - аномально большим или малым значениям». Использование приложения Origin 6.1 к анализу данных таблицу приводит к следующим медианным значениям реологических параметров ротовой жидкости <G> = 0,10785 Па, <q> = 0,32883 с. (7) 0 1 2 3 4 5 6 7 8 9 10 11 12 0,1 1 d, мм t + C, c Рис. 6. Зависимости изменения диаметра нити от времени по данным 8 опытов (символы и номера опытов показаны на врезке). Для наглядности кривые сдвинуты по времени с шагом 1 с. Горизонтальные участки на некоторых кривых показывают размер капли до начала растяжения. Сплошные кусочно-линейные зависимости - аппроксимации экспериментальных точек аналитической зависимостью (6) Результаты обработки данных Номер опыта db, мм ds, мм ts, с G, Па q, c 2222 2,777 0,4977 0,04795 0,2328 0,41364 2223 1,6627 0,5658 0,14098 4,1957 0,24684 7032 2,8485 0,4308 0,08505 0,1440 0,49996 7035 1,319 0,4320 0,24703 4,8745 0,40877 7066 2,336 0,3022 0,07169 0,1241 0,24563 7069 3,9 0,4182 0,1533 0,0856 0,41975 7102 3,0 0,3448 0,0199 0,0538 0,32747 7111 2,0 0,4257 0,04266 0,5730 0,24713 7114 2,0 0,2353 0,03873 0,0916 0,33019 7117 3,0 0,3252 0,03535 0,0476 0,3073 7120 2,477 0,3047 0,03585 0,0836 0,33761 Сравнивая параметры (7) с результатами измерения аналогичных параметров в водных растворах полиакриламида (ПАА) и полиэтиленоксида (ПЭО) [33], можно заключить, что близкие к ротовой жидкости параметры упругости <G> = 0,10785 Па демонстрируют растворы ПАА и ПЭО концентраций 252 и 12449 млн-1, а релаксации <q> = 0,32883 с - растворы ПАА и ПЭО концентрациями 2102 и 19320 млн-1 соответственно. Это означает, что можно подобрать полимеры и соответствующую концентрацию раствора, чтобы приготовить растворы полимеров, моделирующие реальную ротовую жидкость. В частности, раствор ПЭО концентрацией 1-2 % достаточно близок по своей реологии к реальной ротовой жидкости. Заключение Важные следствия вытекают из выполненной работы. Во-первых, впервые определен/оценен полный набор численных значений реологических характеристик ротовой жидкости, а, именно, модуля упругости G и времени релаксации q. (Предыдущие измерения касались только времени релаксации.) Установление физических характеристик ротовой жидкости позволяет выбрать среди возможных растворов полимеров такие, реология которых совпадает с реологией реальной ротовой жидкости. Использование выбранных модельных полимерных жидкостей позволяет экспериментально воспроизводить в лаборатории гидродинамические процессы, протекающие с реальной ротовой жидкостью. Во-вторых, набор G и q «замыкает» реологические уравнения состояний Олдройда-Б или Максвелла с верхней конвективной производной, которые описывают быстрые растяжения упруговязких жидкостей. Именно такие процессы протекают в ротовой жидкости в ходе выброса капель при кашле, чихании, разговоре. Дальнейшая судьба капель, а именно их взаимодействие (в том числе коалесценция и разбрызгивание) с материалом защитных масок, фильтров и других препятствий, также может быть описана с использованием данных моделей. В частности, найденные физические данные помогут предсказать, что будет с каплями в окружающем пространстве, т.е. упадут ли капли на пол, консолидируются ли в защитных средствах, либо разрушатся до уровня аэрозоля и станут удобными объектами для переноса патогенных агентов воздушными конвективными потоками? Описание взаимодействия упруговязких капель с препятствиями может оказаться полезным при разработке новых средств защиты от инфекций. Наконец, на примере ротовой жидкости показана возможность проведения реологических тестов биологических упруговязких жидкостей при помощи простейших (фактически бытовых) инструментов - линейки, смартфона и пинцета. Прогресс в бытовой электронике приблизил данный класс техники к возможностям профессиональных научных приборов. Продемонстрирована процедура определения реологических характеристик жидкостей на основе обработки видеоизображений и численного анализа получаемых зависимостей. В целом алгоритм выполнения представляемого реологического теста чрезвычайно прост и доступен любому исследователю. Фактически тест заменяет реологические испытания, проводимые с использованием относительно дорогостоящего прибора CaBER1 [17]. Методика применима для исследований реологических свойств биологических жидкостей, при растяжении капель которых заметно формирование утончающейся капиллярной нити. К данному классу жидкостей относятся, например, мокрота [7, 10, 39, 40], синовиальная жидкость [19, 34], жидкости, участвующие в репродуктивном процессе (см., например, [20]), и др. Все они могут быть исследованы с применением предложенного метода. Представленный алгоритм биореологических исследований может быть реализован в специальном приложении для смартфона для использования в самых разных условиях, необязательно лабораторных.

About the authors

A. N Rozhkov

References

Statistics

Views

Abstract - 197

PDF (Russian) - 55

Refbacks

• There are currently no refbacks.