Прикладная математика и механика, 2019, T. 83, № 3, стр. 413-427

ТОНКОСТРУКТУРНЫЕ КОМПОНЕНТЫ ВСПЛЕСКА КАПЛИ

А. Ю. Ильиных 1*

1 Институт проблем механики им. А.Ю. Ишлинского РАН
Москва, Россия

* E-mail: ilynykh@ipmnet.ru

Поступила в редакцию 05.04.2019
После доработки 10.04.2019
Принята к публикации 11.04.2019

Полный текст (PDF)

Аннотация

Тонкие структуры переноса вещества свободно падающей капли, погружающейся в глубокую жидкость, изучаются методом фоторегистрации с большими пространственным и временным разрешениями. В результате более 500 экспериментов, проведенных в широком диапазоне основных размерных параметров, для смешивающихся жидкостей отмечено дискретное распределение вещества капли по деформированной поверхности принимающей жидкости и равномерное для несмешивающихся. Дискретное распределение реализуется каскадом тонких волокон – лигаментов, образованных капельным веществом. Предложена общая схема распределения вещества капли для случая взаимодействия смешивающихся пар жидкостей, определено влияние на структуру схемы таких параметров задачи как числа Вебера и Онезорге, относительная разность коэффициентов поверхностного натяжения, а также топографии дна.

Ключевые слова: удар капли, коэффициенты поверхностного натяжения, брызги, капиллярные волны, лигаменты

Введение. В картине всплеска при падении капли в жидкость разделяют следующие компоненты течения: пелена тонких струй, вылетающих из области контакта при соприкосновении жидкостей [1], каверна [2], венец с краевым шевроном, с вершин зубцов которого вылетают вторичные капли и струи [3]. Всплеск сопровождается генерацией поверхностных капиллярных волн на поверхности принимающей жидкости [4], венца и каверны [8] и даже на поверхности первичной капли при отрыве [7] или контакте с принимающей жидкостью, выбросом тонкой и толстой кумулятивных струй [5, 6]. При погружении капли в принимающую жидкость в жидкость и в газовую среду излучаются высокочастотные акустические пакеты [6].

Высокоразрешающей визуализацией области контакта капли с принимающей жидкостью установлен факт попадания на поверхность капли мелких брызг, на ней генерируются группы капиллярных волн. Направление вылета брызг и, следовательно, вероятность их попадания на поверхность капли определяется относительной разностью коэффициентов поверхностного натяжения взаимодействующих жидкостей [9].

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

1. Параметры задачи. Картины течений, возникающих при свободном падении капли в принимающую жидкость (рис. 1), зависят от физических свойств контактирующих сред, которые характеризуются термодинамическими и кинетическими параметрами, а также условиями опытов. В число влияющих физических величин входят следующие: плотности (воздуха ${{{\rho }}_{a}}$, жидкостей капли ${{{\rho }}_{d}}$ и принимающей ${{{\rho }}_{t}}$), кинематические ${{\nu }_{{a,d,t}}}$ и динамические ${{\mu }_{{a,d,t}}}$ вязкости, коэффициенты поверхностного натяжения ${\sigma }_{d}^{a}$, ${\sigma }_{t}^{a}$, коэффициенты диффузии вещества капли $\kappa _{t}^{d}$ и маркирующей примеси $\kappa _{t}^{{md}}$ в принимающей жидкости, температуры (${{T}_{d}} = {{T}_{t}} = {{T}_{a}}$), условия опытов (диаметр $D$ и величина скорости $U$ капли в момент первичного контакта, высота свободного падения H, глубина слоя принимающей жидкости ${{h}_{l}}$, величина ускорения свободного падения $g$). Выделяются дополнительные параметры: масса капли $M$, площадь свободной поверхности ${{S}_{d}}$, величина импульса $p$, кинетическая ${{E}_{k}}$ и доступная потенциальная поверхностная энергия ${{E}_{\sigma }}$ падающей капли.

Рис. 1.

Постановка и размерные параметры задачи.

В опытах использовался широкий набор рабочих жидкостей (водные растворы ализариновых чернил, спирта, неорганических солей, глицерина, уксусной кислоты, подсолнечное и касторовое масла, нефть, молоко и молочные сливки). Значения размерных параметров изменялись в диапазоне: высота падения H = 40–260 см, диаметр капли D = 0.3–0.5 см, глубина слоя жидкости ${{h}_{l}}$ = 0.2–8 см, плотность ${\rho }$ = 0.87–1.26 г см–3, вязкость ${\mu }$ = 0.01–15 г см–1 с–1, коэффициент поверхностного натяжения на границе с воздухом ${\sigma }$ = 27–75 г с–2.

На основе выделенных размерных параметров для определения требований к методике эксперимента проведен масштабный пространственно-временной анализ системы уравнений. Наборы линейных и временных масштабов задачи разделены на подгруппы.

Первая группа линейных масштабов включает капиллярно-гравитационный параметр $\delta _{g}^{\gamma } = \sqrt {\gamma {\text{/}}g} $, входящий в дисперсионное уравнение коротких поверхностных волн, и диссипативно-капиллярные масштабы $\delta _{\gamma }^{\nu } = {{\nu }^{2}}{\text{/}}\gamma $. Вторая группа линейных масштабов, которые зависят от скорости капли, включает кинематический $\delta _{g}^{U} = U_{d}^{2}{\text{/}}g$, а также тонкие капиллярный $\delta _{U}^{\gamma } = \gamma {\text{/}}U_{d}^{2}$ масштаб и масштаб Прандтля $\delta _{U}^{\nu } = \nu {\text{/}}{{U}_{d}}$.

По аналогии, одна часть собственных временны́х масштабов включает только параметры среды – $\tau _{g}^{\gamma } = \sqrt[4]{{\gamma {\text{/}}{{g}^{3}}}}$, $\tau _{\gamma }^{\nu } = {{\nu }^{3}}{\text{/}}{{\gamma }^{2}}$, другая – размер капли – $\tau _{\gamma }^{d} = \sqrt {{{D}^{3}}{\text{/}}\gamma } $, $\tau _{\kappa }^{\gamma } = \kappa D{\text{/}}\gamma $, $\tau _{\nu }^{\gamma } = \nu D{\text{/}}\gamma $, $\tau _{\nu }^{D} = {{D}^{2}}{\text{/}}\nu $, третья группа – ее скорость – $\tau _{U}^{d} = D{\text{/}}{{U}_{d}}$, $\tau _{g}^{U} = {{U}_{d}}{\text{/}}g$. Значения масштабов использованных в экспериментах жидкостей изменяются в широком диапазоне. Для воды наибольший и наименьший масштабы длины составили ${{{\delta }}_{{\min }}}$ = ${\delta }_{{\gamma }}^{{\nu }}$ = = $1.4 \times {{10}^{{ - 6}}}$ см и ${{{\delta }}_{{\max }}}$ = $\delta _{g}^{\gamma }$ = 0.27 см, соответственно, а временны́е масштабы – ${{{\tau }}_{{\min }}}$ = = ${\tau }_{{\gamma }}^{{\nu }}$ = $1.9 \times {{10}^{{ - 10}}}$ с и ${{{\tau }}_{{\max }}}$ = ${\tau }_{\nu }^{D}$ = 17.6 с.

2. Безразмерные числа. Отношения размерных параметров задачи задают наборы безразмерных комбинаций, включающие традиционные числа: Рейнольдса $\operatorname{Re} $ = = ${{UD} \mathord{\left/ {\vphantom {{UD} {\nu }}} \right. \kern-0em} {\nu }}$ = 150–104, Фруда ${\text{F}}{{{\text{r}}}_{d}}$ = ${{U}^{2}}{\text{/(}}gD{\text{)}}$ = 140–200, Вебера $W{{e}_{d}}$ = ${{U}^{2}}D{\text{/}}{{{\gamma }}_{d}}$ = 400–1000, Бонда ${\text{Bo}}$ = $g{{D}^{2}}{\text{/}}\gamma _{d}^{a}$ = 2.4–3.4, Онезорге ${\text{Oh}}$ = $\nu {\text{/}}\sqrt {\gamma D} $ = $1.8 \times {{10}^{{ - 3}}}{\kern 1pt} - {\kern 1pt} 1.7 \times {{10}^{{ - 1}}}$. При различии в составах вещества капли и принимающей жидкости протекающие процессы характеризуются еще одной группой безразмерных комбинаций, включающей относительную разность плотностей контактирующих сред ${{{\text{R}}}_{\rho }}$ = ${\text{(}}{{\rho }_{t}} - {{\rho }_{d}}{\text{)/(}}{{\rho }_{t}} + {{\rho }_{d}}{\text{)}}$ (число Атвуда), коэффициентов поверхностного натяжения ${{\operatorname{R} }_{\sigma }}$ = ${\text{(}}{{\sigma }_{t}} - {{\sigma }_{d}}{\text{)/(}}{{\sigma }_{t}} + {{\sigma }_{d}}{\text{)}}$, динамических вязкостей ${{\operatorname{R} }_{\mu }}$ = ${\text{(}}{{\mu }_{t}} - {{\mu }_{d}}{\text{)/(}}{{\mu }_{t}} + {{\mu }_{d}}{\text{)}}$, которые в экспериментах принимали как положительные, так и отрицательные значения.

Поверхностная энергия, обусловленная анизотропией действия атомно-молекулярных сил на контактных границах, считается равномерно распределенной в тонком слое толщиной порядка размера молекулярного кластера ${{\delta }_{\sigma }} = {{10}^{{ - 6}}}$ см (рис. 2). Основную долю энергии капли (более 97%) составляет кинетическая энергия, распределенная во всем объеме капли, однако, ввиду тонкости границы, плотность поверхностной энергии на несколько порядков превышает плотность кинетической энергии. В момент контакта за короткое время $\Delta {{t}_{s}}\sim {{10}^{{ - 10}}}$ с (диссипативно-капиллярный масштаб времени) происходит уничтожение поверхностных слоев толщиной порядка размера нескольких размеров молекул ${{\delta }_{s}}\sim {{10}^{{ - 8}}}$ см, который соотносится с диссипативно-капиллярным масштабом, и выброс из области слияния отдельных тонких струй. При этом в выражении для термодинамического потенциала Гиббса

$d{{g}_{s}} = - sdT + VdP + \sigma d{{S}_{b}} + {{\mu }_{n}}d{{N}_{n}}$
исчезают дополнительные члены, зависящие от химического потенциала ${{\mu }_{n}}$ и концентрации соответствующих компонентов $d{{N}_{n}}$ [10].

Рис. 2.

Схема распределения термодинамических потенциалов: (а) – полет капли; (б) – начало процесса слияния; I – потенциал $dg$ в толще жидкости, II – потенциал $d{{g}_{\sigma }}$ в приповерхностном слое, III – потенциал $d{{g}_{s}}$ непосредственно на границае раздела – поверхность молекулярного контакта жидкости с газовой фазой, IV – потенциал $d{{g}_{2}}$ слившихся жидкостей в области соприкоснувшихся приповерхностных слоев.

За более длительное время $\Delta {{t}_{\sigma }} = {{\delta }_{\sigma }}{\text{/}}U\sim {{10}^{{ - 8}}}$ с происходит слияние жидкостей (уничтожение приповерхностных слоев) в пятне первичного контакта погружающейся капли с принимающей жидкостью толщиной ${{\delta }_{\sigma }}\sim {{10}^{{ - 6}}}$ см. Преобразование доступной потенциальной поверхностной энергии в другие формы энергии сопровождается формированием течений, повышением температуры и давления в области слияния. При этом термодинамический потенциал приобретает типичную для толщи жидкости форму $d{{g}_{s}} = - sdT + VdP$. По мере дальнейшего погружения капли в жидкость (продолжительностью $\tau \sim {{10}^{{ - 3}}}$ с) характер передачи энергии изменяется – свободная энергия непрерывно поступает лишь в узком кольце контакта приповерхностного слоя капли с принимающей жидкостью.

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

3. Экспериментальная установка. Опыты выполнены на стенде “Тонкоструктурные быстропротекающие процессы” (ТБП), входящем в состав гидрофизического комплекса уникальных установок для моделирования гидродинамических процессов в окружающей среде “УСУ ГФК ИПМех РАН”, схематическое изображение которого приведено на рис. 3.

Рис. 3.

Стенд ТБП для регистрации мелкомасштабных быстропротекающих процессов на поверхности жидкости: оптическая кювета 1 с принимающей жидкостью глубиной ${{h}_{l}}$, оптическая система регистрации 2, управляемая компьютером 7 и блоком регулируемых задержек 3 с датчиком прохождения капли 4, источники света видимого диапазона 5, 6, защитные и отражающие экраны 8, дозатор капель 9.

Исследования проводились в режиме формирования выраженной центральной кумулятивной струи, который задавался путем выбора высоты свободного падения и диаметра капли. Дозатор 9 устанавливался на высоте $H$ = 30–260 см над поверхностью принимающей жидкости.

Область наблюдения освещалась постоянными источниками света видимого диапазона. Для освещения использовались студийные софиты ReyLab Xenos RH-1000 мощностью 1 кВт 6 и светодиодные источники Optronis MultiLED 5 со световым потоком 7700 лм, которые располагались над кюветами или сбоку от них. Требования к освещению и, следовательно, схема расположения осветительной техники и отражающих поверхностей задавались геометрией и масштабами исследуемых течений. Применялись схемы освещения в прямых и в отраженных лучах.

Регистрация картины течений проводилась видеокамерой Optronis CR3000x2, управляемой с компьютера, с частотой съемки от 4 до 40 тыс. кадров/с или фотоаппаратами Canon EOS350D 2 с объективами Canon EFS 18–55 мм и EF 50 мм. Фотоаппараты и видеокамера размещались на штативах с 5 степенями свободы. Настройка видеокамеры или фотоаппарата осуществляется с компьютера 7. Линия визирования регистрирующей аппаратуры в одних опытах проходила вдоль покоящейся поверхности принимающей жидкости, в других – составляла разные углы: незначительные отклонения (до 5°) в обе стороны, 20°–25° к горизонту при съемке сверху (для исследования формы венца и распределения вещества капли по его поверхности), 70°–75° (регистрация сверху) для исследования распределения вещества капли по дну каверны. Расстояние от передней линзы фотоаппаратов до области регистрируемых течений составляла 7–8 см в большинстве экспериментов.

Для увеличения изображения использовались реверсивное кольцо JJC RR-EOS диаметром 58 мм и набор макроколец Canon толщиной 12, 20 и 36 мм в разных комбинациях, позволившие разрешать элементы с линейными размерами порядка 10 мкм. Состав набора дополнительных колец объектива определялся исходя из масштабов и местоположения регистрируемых компонентов течений. Для разрешения воспроизводимых тонкоструктурных компонентов на фотоаппаратах устанавливалась минимальная выдержка 1/4000 с.

Опыты с подкрашенными каплями (водные растворы чернил, перманганата калия, медного и железного купороса, молоко, спиртовые растворы бриллиантовой зелени) и несмешивающимися жидкостями (окрашенные подсолнечное и касторовое масла, нефть) проводились с использованием ряда квадратных чашек Петри размером $12 \times 12 \times 1.5$ см или глубоких емкостей размером $16 \times 16 \times 8$ см для оперативной смены принимающей жидкости и обеспечения чистоты проведения эксперимента.

Для исследования влияния профиля дна на картину растекания капли по деформированной поверхности неглубокой жидкости (${{h}_{l}} < D$) были изготовлены кюветы с полусферическими и коническими углублениями (диаметры основания составляли 10 и 20 мм, угол при вершине 60°), клиновидным вырезом (с углом раствора 60° и основанием 20 мм) и кювета, на дне которой располагался стеклянный шар диаметром 20 мм.

4. Циклические эффекты. При уничтожении свободных поверхностей в выбросе мелких капель наблюдается ряд циклических эффектов. Длительное действие механизма уничтожения свободных поверхностей, ускоряющего течения, предопределяет геометрию поступления вещества капли в растущий венец, а также обеспечивает продолжительность процесса удлинения тонких струек и выброса капелек с их вершин. По мере погружения капли размер мелких капелек увеличивается при уменьшении нормального угла их вылета. Первые группы капелек вылетают насильно (под углом 10° к горизонту на рис. 4а), зубья шеврона с отделяющимися крупными каплями располагаются под углом 45°. Зональная структура капельных облаков указывает на существование радиальных периодических движений, почти однородных по азимутальной переменной.

Рис. 4.

Циклические эффекты в выбросе мелких брызг в фазе первичного контакта: набор струй, выбрасываемых при падении капли раствора чернил в воду (а), кольцевые структуры брызг на поверхности капли водного раствора медного купороса погружающейся в воду (б), зональная структура капельных облаков (в).

Попадание мелких брызг на поверхность погружающейся капли генерирует группы капиллярных волн, периодически расположенных по линии контакта жидкостей. Форма фронта волны, которая определяется величиной и направлением скорости капельки относительно принимающей поверхности, однозначно определяет направление траектории движения источника (рис. 4б). Здесь в силу свойств дисперсионного уравнения капиллярных волн в картине течения выражены опережающие возмущения – дуги, располагающиеся перед источником (в отличие от гравитационных волн, которые отстают от источника). Вторичные капли отрываются от зубьев шеврона эшелонами в несколько этапов (рис. 4в).

Направление и диапазон углов вылета тонких струй и капель определяется динамикой столкновения и действием поверхностных сил, зависящих от относительной разности коэффициентов поверхностного натяжения ${{R}_{\sigma }}$. Экспериментально установлено, что при значениях ${{R}_{\sigma }} \leqslant 0$ часть вылетающих струй попадает на поверхность капли, генерируя группы коротких волн капиллярного масштаба (рис. 5).

Рис. 5.

Капиллярные волны на поверхностях капель чистой воды, погружающейся в раствор спирта (а) и масло (б), а также капли раствора соли – в чистую воду (в).

В случае положительных значений ${{R}_{\sigma }}$ (например, при погружении капель этанола и подсолнечного масла в воду) угол разлета струй меняется, и капиллярные волны на поверхности первичной капли не обнаруживаются.

5. Общая картина течения (смешивающиеся жидкости). В большом числе экспериментов со смешивающимися жидкостями (а также эмульсиями и суспензиями) отмечено, что вещество капли на деформированной поверхности принимающей жидкости распределяется дискретным образом, собирается в тонких линейчатых элементах – лигаментах. В верхнем ряду венца лигаменты образуют вертикальную линейчатую структуру. На поверхности каверны выражена сетчатая структура, включающая трех-, четырех- и пятиугольные ячейки, на дне каверны сохраняется центральное пятно (рис. 6). Линейчатые структуры устойчиво воспроизводятся, наблюдаются во всех фазах процесса слияния капли: на растущей каверне и венце, а также на заполняющейся каверне и погружающемся венце, сохраняются до распада кумулятивной струи. Положение и пространственная ориентация волокон сетки со временем меняются [11].

Рис. 6.

Линейчатые структуры на поверхностях венца и каверны при растекании капель: чернил (а), перманганата калия (б).

В линейчатой структуре распределения вещества капли по внутренней поверхности каверны выделяется три масштаба. Первый из них определяет длину структур ${{\delta }_{l}} = Ut$, второй – диссипативно капиллярный – определяет толщину слоя $\delta _{\gamma }^{\nu } = {{\nu }^{2}}{\text{/}}\gamma $, третий – линейный масштаб тонких окрашенных нитей и полосчатых структур ${{\delta }_{w}}$ – находится в диапазоне между капиллярным и капиллярно-гравитационным масштабами. Полосчатые структуры в процессе эволюции венца распадаются на более тонкие окрашенные нити, толщиной от 0.005 до 0.04 см.

Линейчатые структуры наблюдаются на нескольких начальных фазах процесса слияния капли, на растущей каверне и венце, а также на заполняющейся каверне и погружающемся венце, сохраняются до распада кумулятивной струи. Образование линейчатых структур объясняется действием больших флуктуаций давления в сравнительно малых по объему областях уничтожения свободной поверхности сливающихся жидкостей. Большие флуктуации давления ускоряют совместное тонкослойное течение соприкасающихся жидкостей и формируют неравномерное распределение вещества капли по поверхности венца и каверны.

Для характеристики влияния основных параметров задачи проведена параметризация картин распределения вещества капли по деформированной поверхности принимающей жидкости. Картина распределения вещества капли по внутренним поверхностям венца и каверны представляет собой сетчатую структуру с центральным пятном, ограниченным окружностью R1 (рис. 7). Далее по мере удаления от центра дна каверны располагаются ярусы сетки, количество которых зависит от условий эксперимента и фазы процесса (на схеме приведено пять ярусов $k = 5$). Каждый ярус сетки ${{R}_{i}}{{R}_{{i + 1}}}$ образован ребрами ${{\operatorname{f} }_{i}}$, в которых концентрируется вещество капли. Ребра образуют ячейки сетчатой структуры (преимущественно трех- и четырехугольные, однако встречаются и пятиугольные), в которых вещество капли диффузно распределено в меньшей концентрации или не обнаруживается вовсе. Внешний ярус сетки образован вертикальными линейчатыми структурами ${{{\text{f}}}_{k}}$ (${{{\text{f}}}_{5}}$ на рис. 7), располагающимися на поверхности венца и, частично, каверны.

Рис. 7.

Схема картины распределения вещества капли по внутренней поверхности короны.

Ярусы сетки ограничиваются условными кольцами ${{R}_{i}}$, на которых располагаются узлы сетки (узлы повышенной концентрации вещества капли). Кольцом ${{R}_{{k + 1}}}$ (${{R}_{6}}$) обозначена кромка шеврона венца – граница картины распределения вещества капли. Исходящие из кольца ${{R}_{{k + 1}}}$ линейчатые структуры ${{{\text{f}}}_{{k + 1}}}$ обозначают зубцы шеврона, с вершин которых отрываются вторичные капли, также содержащие вещество основной капли. Положения капель и вершин зубьев шеврона также образуют кольцо.

Параметрами картины распределения вещества капли служат количество ярусов сетки $k$, изменяющееся при условиях опыта и состава жидкостей, диаметры колец ${{R}_{i}}$ и уровень их целостности (отношение окрашенных участков колец к длине окружности).

6. Влияние высоты (скорости) падения капли на картину распределения вещества. С учетом сил тяжести, Стокса и плавучести, действующих на свободно падающую в воздухе каплю воды диаметром $D = 0.42$ см (в пренебрежении отклонением формы капли от сферической), скорость капли меняется с высотой падения, как показано на рис. 8.

Рис. 8.

Характер изменения скорости свободного падения капли с высотой падения.

Значения скорости капли в момент контакта с принимающей жидкостью, капиллярного масштаба, чисел Вебера и Фруда при разных высотах свободного падения капли, устанавливаемых в экспериментах, приведены в табл. 1.

Таблица 1.

Параметры экспериментов (разные высоты падения)

Высота падения, H, см Скорость, U, м/с Капиллярный масштаб $\delta _{U}^{\gamma } \times {{10}^{3}}$, см We Fr
40 2.34 1.3 320 130
80 2.82 0.92 460 190
100 2.94 0.85 500 210
120 3.02 0.8 530 220
150 3.09 0.76 560 230
200 3.15 0.74 580 240
260 3.17 0.73 586 244

Дискретный характер распределения синих ализариновых чернил в принимающей жидкости сохраняется в широком диапазоне высот падения капли (рис. 9). При падении с высоты $H = 100$ см в картине течения в фазах роста и начала расплывания венца выражены собранные в ярусы треугольные и линейчатые ячейки (рис. 9, строка I). Плотность окраски увеличивается по мере приближения к остатку капли в центре каверны.

Рис. 9.

Распределение вещества капли по деформированной поверхности принимающей жидкости в фазах: выброса пелены брызг (а), формирования короны (б), развитой короны. Высоты свободного падения капли, см: 100 (I), 150 (II), 260 (III).

С увеличением высоты падения нарастает контрастность зональной структуры, состоящей из элементов разной формы: и линейчатых, и более сложной геометрии (рис. 9, строка II). На стенках венца наблюдаются практические свободные от краски ячейки. При большой высоте падения кольцевая структура распределения красителя, образованная элементами сложной формы, становится преобладающей уже на стадии роста венца (рис. 9, строка III). Под действием сил поверхностного натяжения вершина венца стягивается в купол.

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

При погружении капли чернил, свободно падающей с высоты 100 см, из контактной линии жидкостей отходят широкие окрашенные линейчатые структуры шириной порядка 0.03 см, между которыми находятся слабоокрашенные области. На рис. 9а в строке I четко визуализируется кольцо перехода из каверны в венец. Линейчатые структуры при переходе через кольцо дробятся на лигаменты шириной порядка 0.005 см.

Лигаменты формируют один ярус сетки $k = 1$, радиус остатка капли ${{R}_{1}} = 0.33$ см, радиус кромки шеврона составляет ${{R}_{2}} = 0.62$ см, отдельные струи, образующиеся на кромке шеврона, достигают размеров ${{{\text{f}}}_{2}} = 0.27$ см. На следующем кадре (рис. 9б, строка I) картина распределения вещества капли становится сетчатой. Надежно определяются 4 яруса ($k = 4$). Значения радиусов центральных пятен и отношения радиусов колец к радиусу центрального пятна ${{R}_{{i1}}} = {{{{R}_{i}}} \mathord{\left/ {\vphantom {{{{R}_{i}}} {{{R}_{1}}}}} \right. \kern-0em} {{{R}_{1}}}}$ приведены в табл. 2. Целостность (заполненность) кольца верхнего яруса ${{\chi }_{{{{R}_{4}}}}}$ (рис. 10, $H = 100$ см), определенная методом фотометрии как отношение окрашенных областей к длине окружности, составляет 59%. С развитием короны целостность крайнего кольца снижается. Число ярусов сетчатой структуры, относительные радиусы колец и их целостности для разных фаз всплеска и разных высот свободного падения капли (I–III) в соответствии с рис. 9 приведены в табл. 2.

Таблица 2.

Значения параметров распределения вещества капли при разных значениях высоты падения, см: 100 (I), 150 (II), 260 (III)

  (а) (б) (в)
I $k = 1$, ${{R}_{1}} = 0.33$ см, ${{R}_{{2i}}} = 1.88$ $k = 4$, ${{R}_{1}} = 0.45$ см, ${{R}_{{2i}}} = 1.13$, ${{R}_{{3i}}} = 1.27$, ${{R}_{{4i}}} = 1.56$, ${{R}_{{5i}}} = 1.78$, ${{\chi }_{{{{R}_{4}}}}} = 59\% $ $k = 4$, ${{R}_{1}} = 0.44$ см, ${{R}_{{2i}}} = 1.38$, ${{R}_{{3i}}} = 1.88$, ${{R}_{{4i}}} = 2.18$, ${{R}_{{5i}}} = 2.43$, ${{\chi }_{{{{R}_{4}}}}} = 34\% $
II $k = 2$, ${{R}_{1}} = 0.38$ см, ${{R}_{{2i}}} = 1.38$, ${{R}_{{3i}}} = 2.01$ $k = 4$, ${{R}_{1}} = 0.63$ см, ${{R}_{{2i}}} = 1.26$, ${{R}_{{3i}}} = 1.43$, ${{R}_{{4i}}} = 1.57$, ${{R}_{{5i}}} = 1.65$, ${{\chi }_{{{{R}_{3}}}}} = 90\% $ $k = 2$, ${{R}_{1}} = 0.75$ см, ${{R}_{{2i}}} = 1.32$, ${{R}_{{3i}}} = 1.36$, ${{\chi }_{{{{R}_{2}}}}} = 83\% $
III $k = 3$, ${{R}_{1}} = 0.63$ см, ${{R}_{{2i}}} = 1.48$, ${{R}_{{3i}}} = 1.76$, ${{R}_{{4i}}} = 2.2$ $k = 2$, ${{R}_{1}} = 0.66$ см, ${{R}_{{2i}}} = 1.22$, ${{R}_{{3i}}} = 1.46$, ${{\chi }_{{{{R}_{2}}}}} = 64\% $ $k = 2$, ${{\chi }_{{{{R}_{2}}}}} = 64\% $
Рис. 10.

Распределение краски в кольце ${{R}_{4}}$ картины течения, образованного при падении капли чернил с высоты: $H = 100$ см и $H = 150$ см.

Подобная сетчатая картина наблюдалась при высотах свободного падения капли $H < 100$ см, We < 500. Однако уже при высоте $H = 120$ см (We = 530) выраженная сетчатая структура приобретает кольцевой характер (целостность крайнего кольца ${{\chi }_{{{{R}_{3}}}}} = 70\% $ при более плотных нижних), а при высоте $H = 150$ см, We = 560 (рис. 10) целостность кольца достигает ${{\chi }_{{{{R}_{3}}}}} = 90\% $.

Далее с увеличением высоты (скорости, числа Вебера) уменьшается число ярусов картины распределения вещества и целостность кольцевых элементов, образованных вихревыми кольцами, связанных тонкими волокнами (рис. 11).

Рис. 11.

Зависимости числа ярусов (а) и целостности колец (б) в картине распределения вещества капли в фазе развитого венца (рис. 9в) от числа Вебера.

В картине распределения вещества капли по деформированной поверхности принимающей жидкости с увеличением высоты свободного падения капли уменьшается число ярусов сетки (рис. 11а): с пяти ярусов при $H = 80$ см до двух при $H = 200$ см. Линейчатые структуры, образующие ячейки, истончаются; вещество капли концентрируется в кольцах, целостность которых увеличивается и достигает максимального значения при $H = 150$ см, при больших значениях $H$ кольцевые структуры распадаются на вихревые кольца (целостность колец падает, рис. 11б).

7. Влияние вязкости капли на картину распределения вещества. Разные концентрации подкрашенного водного раствора глицерина позволили оценить влияние вязкости на картину растекания вещества капли по деформированной поверхности принимающей жидкости. Вязкости растворов и числа Онезорге в экспериментах приведены в табл. 3

Таблица 3.

Параметры экспериментов (разные концентрации глицерина)

Концентрация Вязкость μ, г см–1 с–1 Oh Относительная разность вязкостей ${{R}_{{\mu }}}$
0 0.01 $18 \times {{10}^{{ - 4}}}$ 0
40 0.038 $66 \times {{10}^{{ - 4}}}$ –0.58
90 2.35 0.42 –1
100 15 2.65 –1

По мере увеличения вязкости (концентрации глицерина) материал капли распределяется по поверхности каверны более равномерно при частичном сохранении линейчатых структур на поверхности венца (рис. 12). Картина распределения капли чернил приведена на рис. 12а. Здесь отмечается четыре яруса сетчатой структуры ($k = 4$) с центральным пятном на дне каверны. Отношение радиуса центрального пятна к радиусу каверны составляет ${{R}_{1}}{\text{/}}R = 0.36$.

Рис. 12.

Картины распределения вещества капли при погружении капель: раствора чернил (а), раствора глицерина: 40% (б) и 90% (в), чистого глицерина (г). Высота падения $H = 80$ см.

На следующем кадре (рис. 10б) на поверхности каверны располагается всего три яруса структуры, отношение радиусов составляет ${{R}_{1}}{\text{/}}R = 0.43$. Зависимость числа ярусов структур и отношения радиусов от числа Онезорге приведены на рис. 13.

Рис. 13.

Изменение числа ярусов (а) и отношения радиусов центрального пятна и каверны (б) с числом Онезорге.

С увеличением вязкости капли число ярусов сетчатой структуры уменьшается, как показано на рис. 13а. При 90% концентрации глицерина в капле в структуре распределения вещества остаются лишь вертикальные структуры на внутренней поверхности венца, далее вещество капли распределяется практически равномерно. Центральное пятно заполняет все бо́льшую площадь каверны (рис. 13б).

8. Влияние топографии дна на картину распределения вещества. Погружение капли чернил в тонкий слой воды (${{h}_{l}}\sim D$) показал, что форма дна бассейна оказывает существенное влияние на картину распределения вещества капли лишь в случае непосредственного контакта (рис. 14). Контакт каверны с гранями углубления приводит к искажению формы каверны и сетчатого рисунка – существенно увеличивается отношение радиуса первого яруса к радиусу центрального пятна.

Рис. 14.

Погружение капли раствора чернил в тонкий слой воды (${{h}_{l}} = 2$ мм) с профильным дном с углублением: сферическим $d = 20$ мм (а), коническим $d = 20$ мм (б), сферическим $d = 10$ мм (в).

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

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

При падении капли нефти на чистую поверхность воды шеврон венца оказывается полностью окрашенным, как и внутренняя часть венца. Вторичные капли, вылетающие с вершин шеврона, также преимущественно состоят из вещества капли (рис. 15а).

Рис. 15.

Картина распределения нефти по поверхности воды: чистой (а), запыленной (б), покрытой тонкой нефтяной пленкой (в).

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

Заключение. В картине распределения вещества капли по деформированной поверхности принимающей жидкости важную роль играют быстрые процессы, связанные с уничтожением контактирующих поверхностей и высвобождением части энергии, которая расходуется, в частности, на энергию трансляционного движения новых компонентов течений и их поверхностную энергию. Движение контактной границы при начальном погружении капли осуществляется по гладкой поверхности принимающей жидкости в направлении от области соприкосновения, а позднее, при погружении половины капли – по возмущенной поверхности ко дну формирующейся каверны.

В зависимости от условий эксперимента и свойств сред картина пространственного распределения вещества капли различна. При погружении капли вещества смешивающего ($\kappa > 0$) с принимающей жидкостью формируется сетчатый рисунок с числом ярусов $k \geqslant 2$. Значение $k$ на разных фазах процесса зависит от высоты свободного падения капли, причем уменьшается по мере увеличения высоты падения (для капли раствора чернил: от $k = 5$ при $H < 100$ см до $k = 2$ при $H = 260$ см). Количество вещества, концентрирующегося в узлах сетчатого, возрастает при снижении концентрации в ребрах сетки.

Погружение вязких несмешивающихся жидкостей (нефть, масло) формирует центральное пятно ${{R}_{1}}$, покрывающее практически всю внутреннюю поверхность венца ($\kappa = 0$, $k = 0$), с шеврона венца ${{R}_{{k + 1}}} = {{R}_{1}}$ выступают зубцы ${{{\text{f}}}_{{k + 1}}} = {{{\text{f}}}_{1}}$, с которых выбрасываются вторичные капли. Положение границы ${{R}_{1}}$ зависит от относительного коэффициента поверхностного натяжения сред ${{R}_{\sigma }}$. При ${{R}_{\sigma }} \leqslant 0$ граница центрального пятна проходит по кромке венца шеврона ${{R}_{{k + 1}}}$ или вблизи нее, вторичные капли содержат незначительное количество вещества капли. При ${{R}_{\sigma }} > 0$ граница центрального пятна не совпадает с кромкой шеврона, а лежит на внешней поверхности венца, количество вещества основной капли во вторичных существенно возрастает.

Из результатов эксперимента следует, что параметр $k$ чувствителен к относительной разности вязкостей взаимодействующих сред (уменьшается с увеличением вязкости капли), при этом существенно возрастает диаметр центрального пятна ${{R}_{1}}$.

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

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

Работа выполнена при финансовой поддержке Российского научного фонда (проект 19-19-00598 “Гидродинамика и энергетика капли и капельных струй: формирование, движение, распад, взаимодействие с контактной поверхностью”). Эксперименты проведены на стенде ТБП (стенд для регистрации тонкоструктурных быстропротекающих процессов), входящем в Гидрофизический комплекс для моделирования гидродинамических процессов (УНУ “ГФК ИПМех РАН”).

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

  1. Thompson D.W. On Growth and Form. Cambridge: Univ. Press, 1992. 368 p.

  2. Engel O.G. Crater depth in fluid impacts // J. Appl. Phys. 1996. V. 37. № 4. P. 1798–1808.

  3. Чашечкин Ю.Д., Прохоров В.Е. Тонкая структура всплеска при падении капли на свободную поверхность покоящейся жидкости // Докл. РАН. 2011. Т. 436. № 6. С. 768–773.

  4. Zhu G.-Z., Li Z.-H., Fu D.-Y. Experiments on ring wave packet generated by water drop // Chinese Sci. Bull. 2008. V. 53. № 11. P. 1634–1638.

  5. Ray B., Biswas G., Sharma A. Regimes during liquid drop impact on a liquid pool // J. Fluid Mech. 2015. V. 768. P. 492–523.

  6. Прохоров В.Е., Чашечкин Ю.Д. Генерация звука при падении капли на поверхность воды // Акустич. ж. 2011. Т. 57. № 6. С. 792–803.

  7. Прохоров В.Е., Чашечкин Ю.Д. Динамика отрыва одиночных капель в воздушной среде // Изв. РАН. МЖГ. 2014. № 4. С. 109–118.

  8. Чашечкин Ю.Д., Прохоров В.Е. Гидродинамика удара капли: короткие волны на поверхности венца // Докл. РАН. 2013. Т. 451. № 1. С. 41–45.

  9. Ильиных А.Ю., Чашечкин Ю.Д. Гидродинамика контакта падающей капли со свободной поверхностью жидкости // Изв. РАН. МЖГ. 2016. № 2. С. 3–12.

  10. Teschke O., de Souza E.F. Water molecule clusters measured at water/air interfaces using atomic force microscopy // Phys. Chem. 2005. V. 7. № 22. P. 3856–3865.

  11. Ильиных A.Ю., Чашечкин Ю.Д. Гидродинамика погружающейся капли: линейчатые структуры на поверхности венца // Изв. РАН. МЖГ. 2017. № 2. С. 152–164.

Дополнительные материалы отсутствуют.