Сираев Р.Р., Брацун Д.А. «Возбуждение релаксационных колебаний на искривленной межфазной границе в условиях внутренней задачи» Вычислительная механика сплошных сред, 17, № 1, с. 87-97 (2024)

Численно исследуется колебательный режим концентрационной конвекции Марангони при абсорбции поверхностно-активного вещества (ПАВ) из однородного внешнего раствора внутрь капли воды. Эффект обусловлен действием постоянной силы тяжести, которая способствует оседанию молекул ПАВ в водной среде. Такой вариант колебательной конвекции, возникающей в условиях внутренней задачи, недавно обнаружен экспериментально. В настоящей работе рассматривается случай химически инертной системы, в которой реакции отсутствуют. Эффекты деформации межфазной поверхности считаются незначительными и во внимание не принимаются. Математическая модель включает уравнения Навье–Стокса, записанные в приближении Хеле–Шоу и Буссинеска, уравнения переноса ПАВ в системе. Полагается, что характерное время адсорбции ПАВ мало по сравнению со временами его диффузии в обоих растворах;это позволяет игнорировать формирование поверхностной фазы. Краевая задача включает условие равновесия системы, учитывающее разные значения химического потенциала в фазах. Показано, что капля воды является аккумулятором ПАВ, который диффундирует из органической фазы. Задача реализована в размерном виде с помощью пакета COMSOL Multiphysics с использованием набора физических констант для уксусной кислоты, которая как многие представители карбоновых кислот обладает свойствами ПАВ в отношении воды. Обнаружено, что прямое численное моделирование системы способно воспроизвести релаксационные колебания, наблюдаемые в эксперименте, только при дополнительном феноменологическом предположении о неньютоновской реологии межфазной границы, которое ранее предложено для внешней задачи. Обсуждается возможный физический механизм запаздывания в срабатывании неустойчивости Марангони. Показывается, что периодические колебания генерируются внутри капли за счет конкуренции эффекта Марангони и гравитационно-зависимой конвективной неустойчивости раствора. В результате численного моделирования определены структуры конвективного движения на межфазной поверхности и около нее, получена зависимость от времени интенсивности течения, найден диапазон изменения периода колебаний.