Резонансные колебания цилиндрической жидкой капли в вибрационном поле

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

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

В теоретических работах обычно рассматриваются системы с простой геометрией поверхности раздела. Начало этой тематике положили работы Юнга [1] и Лапласа [2]. Одними из первых работ, где исследовались малые свободные капиллярные колебания систем с поверхностями раздела сред относительно их равновесной сферической формы, были работы лорда Рэлея [3−4] и лорда Кельвина [5], изучавших каплю невязкой несжимаемой жидкости. Было показано, что решение можно представить в виде линейной суперпозиции мод с нормальной зависимостью от времени, в основе которых лежат сферические гармоники. Для собственной частоты колебаний была получена формула Рэлея, связывающая частоту с меридиональным числом сферической гармоники и параметрами системы.

Механическое равновесие жидкого столба (жидкой зоны) и струи относительно малых свободных капиллярных колебаний теоретически и экспериментально исследовались в работах Платэ [6,7] и лорда Рэлея [8,9]. Именно в этих работах было найдено предельное отношение между высотой жидкого столба Ь и его радиусом Я: Ь = 2 к Я. При больших высотах столб становится неустойчивым и разрушается. Это явление получило название неустойчивости Рэлея. Так же в работе Рэлея [8] рассматривается неустойчивость цилиндрической струи, окруженной другой жидкостью. Собственные частоты свободного жидкого столба можно найти в книге Дамба [10].

В последнее время интерес к подобным конфигурациям вызван их использованием в различных технологических процессах. Например, жидкая зона рассматривается при изучении процессов роста полупроводниковых кристаллов [11,12]. Основное внимание в таких работах уделяется течению внутри жидкой зоны при наличии нагрева, вертикальных вибраций, магнитного поля. Также рассматривалась фильтрация через пористую среду [13,14]. Отметим, что подавляющее большинство работ посвящено рассмотрению цилиндрического жидкого столба (жидкой зоны) окруженного газом, влияние которого не учитывается. Таким образом, боковая поверхность столба рассматривается как свободная.

Жидкость с малой вязкостью рассмотрена в [10], где найдено приближенное выражение для декремента затухания свободных колебаний сферической капли. В качестве первого приближения им использовалось решение для невязкого случая, на основе которого определялась вязкая диссипация в объеме жидкости. Случай произвольной вязкости жидкости был исследован в работе [15]- подробное описание может быть найдено в книге [16].

Задачу о колебаниях сферической невязкой несжимаемой жидкости, взвешенной в невязкой несжимаемой жидкости другой плотности, можно найти в [10]. В работе [17] получено дисперсионное соотношение для случая произвольной вязкости жидкостей и поверхности раздела жидкостей, обладающей вязко-эластическими свойствами типа Буссинеска — СкривенаБупара, а также рассмотрен ряд важных предельных случаев.

Та часть работы [17], которая касается влияния вязко-эластических свойств поверхности раздела жидкостей, предварялась исследованиями, приведенными в [10], который рассмотрел колебания сферической капли маловязкой несжимаемой жидкости, покрытой нерастяжимой пленкой, и аналогичной работой [18] о пузыре, отделенном нерастяжимой пленкой от неограниченного объема несжимаемой жидкости малой вязкости. В обоих случаях было найдено, что нерастяжимая пленка на поверхности жидкости существенно увеличивает декремент затухания колебаний. В работе [17] показано, что результаты для нерастяжимой пленки могут быть получены, если считать, что дилатационная вязкость или дилатационная эластичность поверхности раздела жидкостей велики.

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

Частоты осесимметричной моды и коэффициенты затухания собственных малоамплитудных колебаний найдены в работе [20]. Рассматривались большие, но конечные числа Рейнольдса. Затухание было вызвано учетом вязкого пограничного слоя на твердых поверхностях. Частоты и коэффициенты затухания осесимметричной и трансляционной моды собственных колебаний вязкой жидкой зоны в поле тяжести были найдены в работе [21]. В этих работах было показано, что значения частот и коэффициентов затухания увеличиваются с увеличением радиуса капли и уменьшаются при увеличении числа Бонда.

Во всех вышеупомянутых работах рассматривались решения с отделенной временной переменной, так называемые нормальные моды, собственные частоты которых находились из проблемы собственных значений. Оставлялись без внимания переходные решения, которые формально можно получить как линейную суперпозицию нормальных мод, но которые, как целое, обнаруживают иное поведение во времени, чем каждая нормальная мода в отдельности, лишь асимптотически приближаясь к наименее затухающей нормальной моде. Кроме того, анализ нормальных мод не решает вопроса о начальной генерации завихренности поля скорости жидкости, например для системы, которая приводится в движение из состояния покоя. Задача Коши для малых свободных капиллярных колебаний сферической в равновесии капли несжимаемой жидкости, взвешенной в несжимаемой жидкости другой плотности и вязкости, была рассмотрена в [22]. В работе [23] были рассмотрены частные случаи капли жидкости в вакууме и пузыря в объеме жидкости, преимущество которых состоит в том, что они допускают аналитическое решение в виде дифференциального уравнения, содержащего обратное преобразование Лапласа [23], или интегро-дифференциального уравнения [24].

В настоящее время большое распространение получили и активно развиваются методы определения физических параметров жидкости (коэффициент поверхностного натяжения, плотность, коэффициент вязкости и т. п.) по характеристикам свободных колебаний их систем. Например, в [25,26] проводились измерения поверхностного натяжения и вязкости капли жидкости, которая левитировала в другой жидкости с помощью звуковых волн. Акустические волны служат не только для «подвешивания» капли в несмешивающейся жидкости, но и для возбуждения колебаний. Благодаря хорошему согласию с линейной теорией, такая методика используется довольно часто [21].

Так, одним из методов определения коэффициента поверхностного натяжения является измерение частот собственных капиллярных колебаний капли жидкости, см., например, [27,28]. Если капля имеет в равновесии сферическую форму, амплитуда колебаний много меньше размера капли, и жидкость невязкая и несжимаемая, то коэффициент поверхностного натяжения может быть вычислен по формуле Рэлея [3]. Различные факторы могут привести к отклонению частот, рассчитанных по формуле Релея при заданном коэффициенте поверхностного натяжения, от реально наблюдаемых собственных частот. В литературе рассмотрены следующие причины изменения спектра частот: вращение капли жидкости [29]- несферичность формы капли [30]- влияние поля, подвешивающего каплю в поле силы тяжести [31,32]- конечность амплитуды колебаний [33,34]- влияние вязкой диссипации [15]. Требование как можно более точного определения коэффициента поверхностного натяжения обуславливает необходимость уточнения формулы Релея. Часто подобное уточнение не может быть получено в замкнутой аналитической форме, что обусловливает применение методов теории возмущений и численных методов.

С помощью измерения частот свободных капиллярных колебаний может быть определено также поверхностное натяжение на поверхности раздела двух несмешивающихся жидкостей. Для этого можно рассмотреть двухслойную систему, которая в равновесии имеет вид жидкой сферической капли, окруженной концентрическим сферическим слоем другой жидкости, что проделано в работах [35,36]. Теоретический анализ нормальных мод малых свободных капиллярных колебаний в такой системе для случая невязких несжимаемых жидкостей, а также экспериментальные наблюдения были выполнены в работе [37]. Для каждого меридионального числа они нашли две моды свободных колебаний с двумя значениями собственных частот, которые являются двумя решениями некоторого биквадратного уравнения. Рассматривались предельные случаи: капли жидкости в вакууметвердой внутренней каплиповерхности внешней жидкости, граничащей с твердой оболочкойтонкого или толстого слоя внешней жидкости. Такая же система, но для случая несжимаемых жидкостей высокой вязкости была рассмотрена в работе [38], где было найдено приближенное выражение для декремента апериодического затухания колебаний.

Подобные методы сложно развить на основе исследований жидкой зоны, т.к. в этой конфигурации не возможно рассматривать каплю изолированной. С другой стороны, существующие методы определения физических параметров не работают при высоких температурах, а имеющиеся данные очень неточные. Например, при процессе роста полупроводникового кристалла методом плавающей зоны, температура расплава составляет 1000 -3000°С [39,40]. Поэтому в таких случаях, параметры жидкости (расплава) пытаются определить, изучая поведение капиллярных волн на свободной поверхности жидкой зоны [21]. Есть надежда, что, получив в эксперименте частоты и время затухания таких волн, можно определить вязкость и поверхностное натяжение. В работе [41] рассматривалась возможность применения жидкой зоны для определения микро-ускорений при проведении экспериментов по выращиванию кристаллов на борту спутников и космических станций.

Вынужденные колебания. При исследовании систем с поверхностями раздела сред при наличии вибрационного воздействия особое внимание уделяется резонансным явлениям.

Первые работы по этой тематике принадлежат Фарадею, который при рассмотрении влияния одночастотных вертикальных вибраций на горизонтальную границу раздела жидкость — газ экспериментально обнаружил и описал явления возбуждения параметрических колебаний границы раздела (так называемая «рябь Фарадея» или «волны Фарадея») [42]. Теоретическое объяснение этого явления для данной системы было дано в [43] (см. также [44]) с использованием линейного анализа устойчивости для невязкого потенциального течения. Было показано, что уравнения, определяющие поведение системы, эквивалентны системе уравнений Матье, которые исследовались с помощью теории Флоке [45]. Учет малой вязкости производился на основе невязкого течения, то есть путем добавления диссипативного слагаемого в уравнение Матье или вязкого слагаемого в уравнение Бернулли. Полученное дисперсионное соотношение приводило к результатам, согласующимся с экспериментами для маловязких жидкостей [46]. Однако, как показано в работе [47], где на основе полных уравнений выполнен анализ Флоке устойчивости плоской горизонтальной поверхности раздела двух вязких несжимаемых жидкостей, метод Ламба в данном случае не применим, так как существенной становится вязкая диссипация в пограничных слоях, окружающих поверхность раздела, которую этот метод не учитывает. В работе [48], в которой сравниваются и анализируются уравнения Матье для невязкой, слабовязкой и вязкой жидкости. Нелинейная теория возбуждения параметрической неустойчивости для случая малой вязкости жидкости рассмотрена в работе [49].

Горизонтальная граница раздела при вертикальных вибрациях многократно исследовалась в лабораторных экспериментах [50,51] в которых образовывались рельефы на поверхности жидкости как схожие с конвективными структурами (полосы, квадраты, гексагоны, спирали), так и необычные структуры типа треугольников и пр. [52]. В работе [53] было найдено, что при высоких частотах вибраций основным эффектом является сжатие капли вдоль направления оси вибраций, причем этот эффект имеет место как для тяжелой, так и для легкой капли. При двухчастотном (частоты были близки между собой) акустическом взвешивании сферической капли [54,55] наблюдались похожие явления: капля сжималась вдоль вертикальной оси и окончательно принимала форму яблока (dog-bone form). Если разность частот акустического поля была близка к частоте квадрупольной моды собственных колебаний, то происходило резонансное возбуждение квадрупольных колебаний капли. Как показано в [54], это резонанс вынужденных колебаний, а не параметрический резонанс. В работе [56], как и в предыдущей работе, экспериментально исследовалось поведение капель жидкости в одночастотном акустическом поле. Но в отличие от указанных выше работ, исследовался также процесс разрушения капли. Особое внимание уделялось изменению равновесной формы капли при изменении величины акустического давления и механизма разрушения кали в интенсивном акустическом поле. Когда отношение экваториального и равновесного радиуса составляло примерно 1.3, то капля представляла собой фактически плоский диск. На поверхности этого диска параметрическим образом возбуждалась рябь, вызванная взаимодействием двух основных волн, одна из которых движется в радиальном направлении, а другая в азимутальном. Однако эта неустойчивость имеет природу отличную от параметрической неустойчивости Фарадея.

При исследовании поведения цилиндрической капли несжимаемой невязкой жидкости между двумя твердыми поверхностями в одночастотном вибрационном поле [57] был получен для средней формы эффект аналогичный [53]: капля сжималась вдоль оси вибраций независимо от соотношения плотностей. Необычное поведение капель было обнаружено в экспериментах, описанных в работе [58]: тяжелая капля вытягивалась вдоль оси вибраций, а легкая сжималась.

Параметрический резонанс, обнаруженный в [42], наступал в случае равенства внешней частоты удвоенной частоте собственных колебаний б) = 2С1п, со — частота вибраций, Ол — частота п — моды собственных колебаний. Это результат линейного (по амплитуде вибраций) анализа неустойчивости [43,59,60]. При исследовании нелинейного поведения и анализа устойчивости во втором порядке вводят в рассмотрение еще одну моду собственных колебаний. Можно классифицировать соотношением частот 2-х мод форму поверхности, образуемую параметрическим возбуждением волн Фарадея (так называемые внутренние резонансы) [46]. Свойства системы подбирают так, чтобы отношение частот С1п и От этих двух мод составляло определенную величину. Например, в [61] рассматривался случай, когда соотношение частот 1:1, а в [62] 1:2 как наиболее опасные в случае полубесконечного слоя жидкости. В работе [63] в слое конечной глубины рассматривалась неустойчивость в третьем порядке при соотношении частот 1:3, которые для такого случая становиться более опасными, чем рассмотренные в предыдущих работах.

В [64], в отличие от большинства работ, рассматривалась амплитуда вибраций поверхности и объема сравнимая с размерами газового пузыря при одночастотном воздействии. Основная цель работы заключалась в исследовании нелинейного взаимодействия между осцилляциями объема и формы поверхности пузыря, которые возникали или из-за деформации формы или при неравновесном состоянии объема. Для пузырей сферической формы рассматривалось условие резонанса 2:1, для имеющих несферическую форму 1:1 и 2:1, т. е. или О.0−2С1т (результат малоамплитудного приближения), где О0 — частота радиальной моды собственных колебаний, — частота т — ой моды собственных колебаний формы.

В работе [65] исследована параметрическая неустойчивость вынужденных колебаний формы почти сферической каплипоказано, что для главных резонансов неустойчивость появляется при выполнении условия синхронизма со = Ои + Оя+, со — частота вибраций,, Ол+1 — частоты двух соседних мод собственных колебаний. Такое условие резонанса связано с тем, что моды собственных колебаний взаимодействуют через трансляционную моду.

В книгах [66−68] описано явление супергармонического (или нелинейного, в терминах [69]) резонанса в механических системах, который возникает при условии 2а> = &euro-1т, т. е. при равенстве частоты вибраций половине частоты собственных колебаний. В работе [57] было обнаружено подобное явление: возникал нелинейный резонанс при равенстве внешней частоты половине частоты основной моды собственных колебаний цилиндрической капли. В [70] рассматривался линейный со =, супергармонический 2со = С1 и параметрический резонанс со = 2П при исследовании колебаний в нелинейной механической системе с одной степенью свободы. Рассматривалось влияние различных параметров на устойчивость системы и появление резонансных колебаний.

Задачам с многочастотным внешним вибрационным воздействием в последнее время уделяется пристальное внимание. При экспериментальном изучении волн Фарадея в [71] при использовании двух кратных частот (тсо и псо) вибраций были получены новые необычные рельефы поверхности при наличии внутренних резонансов т: п=4:5 и 6:7. Теоретическое объяснение этих экспериментов дано в [72].

В [73] использовалось рассеивание акустической модулированной волны с несущей частотой со, и модуляционной частотой сор на поверхности газового пузырька для измерения его размеров и плотности. Если постепенно приближать частоту сор к частоте собственных колебаний, то амплитуда колебаний стенок пузырька нарастает и достигает максимума при частотах колебаний стенок щ±со (линейный резонанс). В [74] рассматривался параметрическии резонанс со, ±-сор/ 2 в подобной системе.

В [75] на примере уравнения Дуффинга рассматривалась устойчивость системы с двумя положениями равновесия при двухчастотном вибрационном воздействии с большой и маленькой частотой. Проводились сравнения с поведением такой же системы под действием низкочастотных вибраций и наложением высокочастотного шума.

В [76] изучалось поведение системы с двумя степенями свободы и с квадратичной и кубической нелинейностью при двухчастотном внешнем воздействии. Рассматривались два случая: в отсутствии внутреннего резонанса и с различными возможностями внутреннего резонанса (1:1, 2:1, 3:1, 1:2, 1:3).

В [77] исследовался линейный резонанс нелинейного пружинного маятника при внешнем двухчастотном воздействии. Такая система, например, может моделировать поперечную качку судна. Приведена зависимость роста резонансной амплитуды от параметров системы.

Работа [78] посвящена поведению нелинейного осциллятора при многочастотном воздействии. Рассматривались три случая: (1) внешние частоты далеки друг от друга, но близки к собственной частоте осциллятора,.

2) внешние частоты близки друг к другу, но далеки от собственной частоты,.

3) внешние частоты близки друг к другу и собственной частоте. Исследовалась устойчивость вынужденных колебаний, рассматривалась влияние вязкости и параметра расстройки частоты.

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

Влияние малой вязкости жидкостей исследовалось в работе [65,80] феноменологически, искусственным добавлением в амплитудные уравнения членов, отвечающих за диссипацию энергии, что приводит к появлению конечного порога по амплитуде вибраций и сдвигу резонансной частоты. Малая вязкость оказывает дестабилизирующее влияние, сдвигая область неустойчивости в область устойчивости для невязкого случая. Было обнаружено, что в случаях нулевой вязкости и вязкости, стремящейся к нулю, области неустойчивости в пространстве параметров «амплитуда вибраций — расстройка частоты вибраций» не совпадают. Подобный эффект был обнаружен и при исследовании параметрической неустойчивости полусферической капли маловязкой жидкости на вибрирующей подложке.

81]. В последней работе учесть малой вязкость проводился строго, аналогично [49].

Экспериментально стоячие волны на свободной поверхности вибрирующего жидкого моста впервые рассмотрел [82]. Он же впервые попробовал определить частоты собственных колебаний осесимметричной моды. Подобные эксперименты, выполненные позднее описаны в работах [83,84]. В работах [85,86] экспериментально определялись частоты собственных колебаний жидкого моста окруженного другой жидкостью близкой плотности.

В работе [87] рассматриваются нелинейные вынужденные колебания конечной амплитуды вязкого жидкого моста в поле тяжести, верхний торец которого совершает гармонические колебания. Особое внимание уделялось линейному резонансу на собственных частотах осесимметричной моды. Резонансная частота определялась по максимальному значению отношения коэффициента деформации и амплитуды вибраций. Коэффициент деформации определялся как интеграл по все высоте капли от абсолютного значения отклонения боковой поверхности от равновесной формы. Наблюдается сдвиг частот собственных колебаний в сторону уменьшения по сравнению с линейной теорией [21]. В то же время увеличиваются значения коэффициентов затухания. С ростом амплитуды вибраций уменьшается значение резонансной частоты. Слабо нелинейный анализ осесимметричных вынужденных колебаний моста маловязкой жидкости проведен в работе [88].

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

Динамика контактной линии. Первой работой, посвященной поведению контактного угла, была работа [1], в которой удалось связать значение равновесного краевого угла со значениями свободной энергии поверхностей раздела фаз (фактически, с коэффициентами поверхностного натяжения), не зная ничего о структуре переходной зоны. В этой работе было получено так называемое уравнение Юнга (или уравнение Юнга — Лапласа), которое связывает разность давлений на поверхности раздела и кривизну поверхности. Подобное уравнение независимо получил Лаплас [90,91].

Работы, посвященные поведению контактного угла и движению линии контакта, можно разделить на два типа. В первом рассматривается установившееся движение контактной линии по подложке. Растекание жидкости в этом случае происходит или за счет межчастичного взаимодействия (например, ван-дер-ваальсовое притяжение молекул жидкости к подложке) или за счет статических внешних воздействий (например, сила тяжести). Подробный обзор основных методов и приближений, применяемых в теории смачивания, приведен в [90]. Ранние экспериментальные [92] и теоретические работы [93] посвящены растеканию жидкости при малых капиллярных числах Са. К настоящему времени имеется довольно много экспериментальных работ [94], в которых рассматриваются процессы при более интенсивном движении и при конечных значениях краевого угла.

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

При прямом равновесном краевом угле обычно используют эффективное граничное условие, предложенное в [95] при рассмотрении стоячих капиллярно-гравитационных волн между двумя вертикальными стенками. Это условие предполагается линейную зависимость между скоростью движения контактной линии и отклонением краевого угла от его равновесного значения: о) д (дх к ' где ^ - отклонение поверхности раздела от равновесного состояния, х — нормальная к твердой стенке координата, увеличивающаяся вглубь жидкости, Л — феноменологическая постоянная (так называемый параметр смачивания или капиллярный параметр), имеющая размерность скорости. В трехмерном случае правая часть (1) принимает вид ЛЯУ^, где п — внешняя нормаль к твердой поверхности. Частными случаями граничного условия (1) являются условия фиксированной контактной линии (Л = 0) и постоянного краевого угла (Л = оо). Как показано в [95], за исключение этих двух случаев, граничное условие (1) приводит к затуханию колебаний.

В указанной работе также проводится качественное сравнение с экспериментальными работами разных авторов. В частности, показано, что заметные расхождения между экспериментальными измеренными и теоретическими рассчитанными (учитывалось только вязкое трение в пограничном слое) декрементами затухания колебаний можно объяснить диссипацией энергии в окрестности контактной линии. При сравнении с результатами экспериментов [96] показано, что шлифовка поверхности увеличивает параметр смачивания от значений порядка единицы до нескольких десятков. Таким образом, параметр Л характеризует не только химические свойства пары материалов жидкость — подложка, но и степень обработки последней.

В работах [97−98] предполагается, что изменение контактной линии может происходить не в фазе с изменением краевого угла, т. е. капиллярный параметр может быть комплексным.

Частоты и коэффициенты затухания осесимметричной моды малоамплитудных собственных колебаний жидкого моста рассматривались в работе [20]. Динамика контактной линии учитывалась с помощью эффективного граничного условия (1). В тоже время, принималось во внимание существование тонкого вязкого пограничного слоя на твердой поверхности. Контактная линия движется по поверхности погранслоя. Значения частот собственных колебаний уменьшались с ростом капиллярного параметра. Наибольшие частоты имеет жидкий мост с закрепленной контактной линией. С увеличением отношения равновесного радиуса и высоты капли значения частот возрастали. Не до конца остается выясненным вопрос правомерности учета вязкого пограничного слоя на твердой поверхности. Эффективные граничные условия как раз и предназначаются для того, чтобы не рассматривать процессы диссипации в самом пограничном слое.

В работах [99,100] использовано более сложное граничное условие, допускающее неоднозначную зависимость краевого угла от скорости движения контактной линии:

Здесь в = д?/дх — отклонение краевого угла от равновесного значения (в [100] также предполагается, что равновесный краевой угол является прямым). Ниже полагаем вх~вг=вс. Как видно из условия (2), в этом случае малые отклонения краевого угла не приводят к движению контактной линии. Отметим, что граничное условие (2) соответствует неподвижной контактной линии, как в случае, А = 0, так и при достаточно большом вс, и условию (1.1), при Ос —^ 0. Результаты экспериментальных работ [101,102] хорошо описываются формулой (2) при малых отклонениях краевого угла.

В [103] при теоретическом рассмотрении затухания стоячих волн в закрытом сосуде, было сделано предположение о том, что контактный угол а/.

0,-02<<9<<9, — Л (<9 + <92),<9<-02.

2) может не изменяться в процессе движения, т. е. контактная линия скользит по стенке. При сравнении результатов с экспериментами [101] было выделено два режима движения: i) контактный угол пропорционален гидродинамической скорости контактной линии v при v < vc cos0e (v)-cos0c=-?— (3) ii) контактный угол не зависит от модуля скорости при v > vc: cos Qc (v) — cos (4) M.

Т.о. изменение краевого угла не зависят от скорости движения контактной линии при больших амплитудах колебаний. Отметим, что при вс=я/2 и 0е (у) -вс<^, условие (3) в точности совпадет с условием (1).

В работе [104] рассматривались аналогичные режимы поведения капиллярно-гравитационных волн. В малоамплитудном режиме результаты эксперимента хорошо описывались условием (1). Результаты при больших амплитудах объяснялись в рамках моделей [95,98,103] с введением некоторого эффективного капиллярного коэффициента Xeff.

В [105] обнаружено три режима движения контактной линии: при малой амплитуде колебаний контактная линия фиксирована. При увеличении амплитуды наблюдается согласованное по фазе (однако, не удовлетворяющее условию (2)) изменение скорости контактной линии и краевого угла. При большой амплитуде колебаний наблюдается свободное движение контактной линии, описываемое формулой (4). Расхождение экспериментальных данных и формулы (2) может быть, по-видимому, объяснено тем, что пороговое значение отклонения краевого угла вс лежит за пределами применимости линейной теории.

Начиная с работ [6−9], выполненных еще в XIX веке, интерес к жидкому мосту не ослабевает. Однако, как следует из приведенного обзора, основное внимание уделяется капле окруженной газом, влияние которого на движение капли не учитывается. В подавляющем большинстве работ при рассмотрении собственных колебаний исследуется только осесимметричная мода. При изучении вынужденных колебаний обычно используется осесимметричное вибрационное воздействие. Также автору диссертации не известны работы по изучению параметрической неустойчивости вынужденных колебаний жидкого столба.

Общая характеристика работы.

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

117 Выводы.

В первой части этой главы рассматривалась модельная задача. Изучалось поведение цилиндрической капли окруженной другой жидкость между двумя параллельными плоскостями. На поверхности раздела учитывался тонкий вязкий пограничный слой. Контактная линия была свободной.

Рассмотрены собственные колебания капли. Полученные частоты собственных колебаний совпадают с частотами цилиндрического жидкого столба [10]. В первом порядке разложения найдена поправка к частоте, которая вызвана диссипацией энергии в вязком пограничном слое.

Исследована устойчивость вынужденных колебаний по отношению к малым возмущениям. Параметрический резонанс наступал при выполнении условия синхронизма: частота вибраций равняется сумме частот двух соседних мод собственных колебаний. Найдено выражение, описывающие резонансные области. Показана, что малая вязкость приводит к появлению порога амплитуды вибраций и сдвигу области неустойчивости при сравнении с нулевой вязкостью. В результате этого, малая вязкость, как видно из рис. 3.2, приводит к дестабилизации вынужденных колебаний.

Область неустойчивости в предельно малой вязкости не совпадает с областью для невязкого случая (см. рис. 3.3). Подобный результат был также обнаружен и для сферической [80] и полусферической капли [81].

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

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

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

Напомним, что m — азимутальное волновое число, к, I — осевое волновое число. В случае продольных вибраций параметрический резонанс возникал при условии равенства частоты вибраций сумме частот одной моды: со = Q, +Q ,.

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

Заключение

.

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

Движение цилиндрической капли идеальной жидкости со свободной контактной линией под воздействием многочастотных вибраций рассматривалось в первой главе.

1. Показано, что в главном порядке разложения по амплитуде вибраций форма капли не изменяется, и она смещается, как целое (трансляционная мода).

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

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

4. Исследована устойчивость вынужденных колебаний относительно малых возмущений. Показано, что параметрический резонанс появляется при выполнении условия синхронизма: одна из внешних частот равна сумме частот соседних мод собственных колебаний. Найден параметр расстройки и построены области неустойчивости.

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

1. В первой части главы рассмотрены колебания и осредненная деформация цилиндрической капли жидкости, окруженной другой жидкостью и находящейся между двумя твердыми поверхностями. Равновесный краевой угол равенг/2.

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

3. Обнаружено, что действительные части частот (т = 2,3,4,.) могут обращаться в ноль в некотором интервале значений капиллярного параметра. Ширина этого интервала зависит от пространственного параметра Ь.

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

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

6. Показано, что капля может в среднем растягиваться или сжиматься вдоль оси вибраций в зависимости от соотношения плотностей жидкости, внешней частоты и капиллярного параметра. Так, тяжелая капля может растягиваться, а легкая сжиматься, что качественно согласуется с результатами экспериментов в работе [57].

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

8. Показано, что набольшими собственными частотами обладает капля с прямым равновесным краевым углом.

В третьей главе диссертации рассматривалось влияние вязких пограничных слоев на устойчивость вынужденных колебаний капли. В первой части рассматривалась модельная задача. Изучалось поведение цилиндрической капли окруженной другой жидкость между двумя параллельными плоскостями. На поверхности раздела учитывался тонкий вязкий пограничный слой. Контактная линия была свободной.

1. Исследована устойчивость вынужденных колебаний по отношению к малым возмущениям. Параметрический резонанс наступал при выполнении условия синхронизма: частота вибраций равняется сумме частот двух соседних мод собственных колебаний. Найдено выражение, описывающие резонансные области. Показана, что малая вязкость приводит к увеличению области неустойчивости при сравнении с нулевой вязкостью.

2. Область неустойчивости в предельно малой вязкости не совпадает с областью для невязкого случая. Подобный результат был также обнаружен для сферической [79] и полусферической капли [80].

3. Рассмотрено поведение капли вязкой жидкости со свободной поверхностью, зажатой между двумя твердыми параллельными поверхностями с закрепленной контактной линией под воздействием поперечных продольных вибраций. Учитывался вязкий пограничный слой на поверхности раздела твердая поверхность — капля.

4. Рассмотрены собственные колебания. Получены частоты собственных колебаний. В первом порядке разложения найдена поправка к собственной частоте, вызванная наличием диссипации в вязком пограничном слое.

5. Исследована неустойчивость вынужденных колебаний капли по отношения к малым возмущениям. При поперечных вибрациях, как и в первой части, параметрическая неустойчивость возникала при выполнении условия синхронизма: а = Птк+Пт+11.

В случае продольных вибраций параметрический резонанс возникал при условии равенства частоты вибраций сумме частот одной моды:

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