В настоящее время понятие устойчивости в физике представляется одним из центральных и ключевых при описании картины мира в целом и частных явлений [1]. Наблюдаемые пространственные структуры, воспроизводимое временное поведение физических систем, так или иначе, характеризуются свойством устойчивости. Наверное, остановиться на каком-то узком определении устойчивости при таком разнообразии сторон физического мира просто невозможно. При изучении различных явлений, процессов, устойчивых структур приходится сталкиваться с переходными стадиями, а также с физическими системами, находящимися в метаста-бильных состояниях, но их описание по-прежнему является неотъемлемой частью науки об общем свойстве системг характеризоваться, некоторой степенью устойчивости. Тем не менее, в первую очередь описанию подлежат воспроизводимые, т. е. в определенном смысле устойчивые состояния и структуры, т.к. именно они отражают картину мира на каждом промежутке времени. В книге [2] устойчивость определяется как «свойство системы быть невосприимчивой к малым возмущениям». Это значит, что основополагающее значение имеет вопрос об эволюции возмущений: система устойчива, если возмущения затухают со временем [3,4]. В некотором смысле данная диссертация тоже посвящена обсуждению вопросов устойчивости различных физических процессов и явлений, а. именно, динамике однородных жидкостей и смесей, пространственным структурам, теплои массопереносу в конвективных системахпри наличии дополнительного воздействия, со стороны переменных силовых полей различной природы. К их числу можно отнести различные вибрационные воздействия, переменные или вращающиеся электрические и магнитные поля, переменный нагрев границ полости и т. д. Несмотря на разнообразие ситуаций, наложение переменного силового поля на конвективную систему, чаще всего приводит в высокочастотном пределе к некоторому отличному от нуля осредненному действию [5]. В частности высокочастотные колебания полости, заполненной жидкостью, при наличии пространственной неоднородности температуры могут вызывать регулярные осредненные движения даже в невесомости. В настоящее время этот результат действия вибраций принято называть явлением термовибрационной конвекции [6, 7]. Возникающее при этом конвективное течение состоит из двух компонент — колебаний с частотой вибраций и осредненного течения. Если период колебаний много меньше всех гидродинамических и тепловых времен, а амплитуда смещения в некотором смысле мала, то может быть применен метод осреднения [8], который позволяет получить замкнутую систему дифференциальных уравнений для осредненных полей скорости, температуры и давления. В теории? тепловой? конвекции этот метод был впервые развит в работе С. М ЗеньковскошшИЖСимоненко [9] дляизучения влияния высокочастотных вибраций на конвективную устойчивость равновесия горизонтального подогреваемого снизу слоя жидкости. При определенных условиях, когда «медленная» составляющая скорости равна нулю возможно состояние механического квазиравновесия. Находясь в состоянии квазиравновесия, жидкость в среднем неподвижна, однако, вследствие вибрационного воздействия имеют место пульсации скорости, температуры и давления. Если неоднородность температуры достаточно велика, то квазиравновесие становится, неустойчивым и в жидкости пороговым образом может возникать некоторое осредненное конвективное течение.
Вне зависимости от природы высокочастотного переменного силового поля его осредненное действие, как правило, нельзя игнорировать в ходе реализации различных технологических процессов [10], а также, при проектировании технических устройств [11], в которых присутствует конвективный тепломассоперенос.
Таким образом, в данной диссертации предпринята попытка рассмотрения с единых позиций влияния различных высокочастотных переменных силовых полей на конвективные процессы в широко используемых гидродинамических системах. Диссертация состоит из восьми глав, включая литературный обзор, 96 иллюстраций и списка литературы из 202 наименований, который, конечно, не исчерпывает всего разнообразия различных ответвлений рассматриваемой тематики.
В первой главе производится обзор теоретических и экспериментальных работ по влиянию переменных силовых полей на различные гидродинамические системы, в том числе, рассматриваются работы по вибрационной конвекции. Во второй части I главы методом осреднения выводятся уравнения термовибрационной конвекции для однородных жидкостей.
Результаты исследования, устойчивости наклонного слоя жидкости относительно произвольных пространственных возмущений, полученные с помощью уравнений термовибрационной конвекции, представлены в главе II. Рассматриваются два варианта высокочастотного вибрационного воздействия: вертикальные и поперечные слою вибрации. Производится расчет валиковой конвекции для разных углов наклона слоя и анализируется смена конвективных режимов при увеличении надкритичности. Обсуждаются результаты численного исследования устойчивости плоскопараллельного течения в наклонном слое при подогреве сбоку относительно плоских и спиральных возмущений.
В главе III разработана численная методика и в невязком приближении выполнен расчет спектра собственных частот капиллярного моста. Найденные значения резонансных частот использовались позднее при попытке разработки методики управления процессом выращивания кристаллов по методу расплавленной зоны при помощи высокочастотных вибраций.
В главе IV изучается влияние вращающегося магнитного поля на конвективные течения в цилиндрической жидкой зоне. В первом параграфе получены аналитические выражения для компонент вектора магнитной индукции и силы Лоренца для произвольных частот вращения магнитного поля. Найдены ограничительные рамки диапазона параметров, когда магнитное поле в жидкой зоне можно считать однородным в каждый момент времени. В пределе относительно высоких частот получены выражения для осредненной силы Лоренца. Для расчета конвективных течений внутри расплава выведена формула, значительно упрощающая вычисление силы Лоренца в приближении однородного поля.
Во втором параграфе четвертой главы излагаются результаты расчета термокапиллярного конвективного течения в цилиндрической жидкой зоне, находящейся во вращающемся магнитном поле в условиях невесомости. Расчеты, выполненные численно методом сеток, показали, что вращающееся магнитное поле приводит к появлению азимутального движения в расплаве, при этом интенсивность конвективного движения в меридиональной плоскости уменьшается. Выяснены особенности влияния вращающегося магнитного поля на распределение примеси в жидкой зоне, в том числе, вблизи движущегося фронта кристаллизации.
В пятой главе представлены результаты исследования свободной и термовибрационной конвекции в ячейке Хеле — Шоу при подогреве снизу. Изучены устойчивые и переходные режимы тепловой конвекции в полости с разным соотношением сторон. В широком диапазоне значений управляющих параметров задачи построены карты устойчивости конвективных режимов. Теоретически подтверждено существование новых типов устойчивых конвективных режимов, которые были названы пульса-ционными течениями. Изучены различные причины нарушения симметрии типичного четырехвихревого режима с перезамыканием вихрей.
Большое внимание уделено переходным течениям. Рассмотрены прикладные аспекты теории вибрационной конвекции.
В главе VI теоретически исследовано влияние на тепловую конвекцию в связанных каналах, оказываемое наличием примеси в жидкости. В широком диапазоне управляющих параметров численно методом конечных разностей изучены переходные режимы и установившиеся конвективные течения. Показано, что в зависимости от величины надкритично-сти в бинарной смеси возможны как стационарные, так и колебательные установившиеся режимы тепловой конвекции. В рамках предложенной модели удалось теоретически описать «перебросовые» колебания, наличие которых демонстрирует достаточно «сложное» поведение гидродинамической системы при малых надкритичностях. Особое внимание было уделено расчету распределения концентрациипримеси" поперек каналов. В-ходе расчетов было показано, что именно явление термодиффузии ответственно за колебательный характер конвекции вблизи-порога.
Исследование влияния высокочастотных вибраций на порог конвекции и форму надкритических движений бинарной смеси в связанных каналах при подогреве снизу проводится в VII главе. В первую очередь рассматривается влияние вертикальных вибраций на колебательные течения вблизи порога. В VIII главе численно и аналитически рассмотрено влияние вертикальных высокочастотных вибраций на тепловую конвекцию бинарных смесей в ячейке Хеле — Шоу, имеющей широкие грани высокой теплопроводности. Проведена аналогия между явлениями, наблюдавшимися в связанных каналах и ячейке Хеле — Шоу. В заключении перечислены наиболее ценные с точки зрения автора научные результаты, полученные в ходе работы над диссертацией.
Основные результаты исследований изложены в 93 различных печатных работах, в том числе, в 16 статьях, опубликованных в реферируемых журналах, учитываемых ВАК при защитах докторских диссертаций.
В подавляющем числе публикаций [141−143], [153], [155], [157], [167], [168], [175], [176], [183−187], [189], [190], [195−200], [202] представлены результаты совместных экспериментальных и теоретических исследований. Все эксперименты выполнялись соавторами: И. А. Бабушкиным или А. Ф. Глуховым. В постановке некоторых экспериментов принимал участие Г. Ф. Путин. Вся теоретическая часть в этих работах принадлежит автору диссертации. Работы [90−92], [113], [154] были сделаны без соавторов- [120], [151], [152], [156], [166], [170] подготовлены совместно со студентами, выполнявшими дипломные и курсовые работы. В теоретических работах [86], [89], [104], [112], [122], [132] автор участвовал в постановке задач, выполнял расчеты и проводил обобщение результатов. Материалы диссертации докладывались на многочисленных научных конференциях, конгрессах и симпозиумах. Кроме того, автор неоднократно выступал с докладами, на Пермском городском, гидродинамическом семинаре под руководством проф. Г. З: Гершуни (1994, 1995, 1996, 1997), и Д. В. Любимова (2002, 2004, 2005, 2006, 2008, 2009), а также в Институте механики сплошных сред УрО РАН на семинаре под руководством проф. А. Ф. Пшеничникова (2005).
Автор глубоко благодарен своему учителю — Г. З. Гершуни за вовлечение в науку и незабываемые минуты общения. Также особую признательность хотелось бы выразить заведующему кафедрой теоретической физики Д. В. Любимову, благодаря чуткому руководству которого у автора сложились условия, способствовавшие работе над диссертацией. Автор выражает признательность И. А. Бабушкину, А. Ф. Глухову и Г. Ф. Путину за плодотворное сотрудничество в плане предоставления возможности теоретику соприкоснуться с условиями реального эксперимента и огромное спасибо всем членам Пермской гидродинамической школы за неоценимую моральную поддержку, благодаря которой смогла состояться данная диссертация.
Заключение
.
В диссертации предпринята попытка рассмотрения с единых позиций влияния различных переменных силовых полей на конвективные процессы в широко используемых гидродинамических системах. Автор стремился показать, что, несмотря на разнообразие конвективных систем (наклонный слой жидкости, ячейка Хеле — Шоу, жидкая зона, конвективная петля), в высокочастотном пределе воздействие переменных силовых полей различной природы приводит чаще всего к некоторому существенному осредненному действию. Вне зависимости от природы переменного силового поля это осредненное воздействие нельзя игнорировать при проектировании различных технических устройств [202] и оптимизации технологических процессов. Перечислим теперь основные результаты диссертационной работы.
1. Проведено исследование устойчивости механического квазиравновесия наклонного слоя жидкости, находящегося в статическом поле тяжести и подвергающегося действию высокочастотных вибраций относительно нормальных пространственных возмущений. Показано, что при подогреве снизу в случае поперечных слою вибраций наиболее опасными всегда являются спиральные возмущения, в то время как для вертикальных вибраций помимо плоских и спиральных возмущений более сложные пространственные моды могут отвечать за порог устойчивости квазиравновесия в достаточно широком диапазоне углов наклона слоя. В нелинейной постановке исследованы конвективные режимы в плоском наклонном слое конечной длины при подогреве снизу и поперечных вибрациях. Усиление, нагрева при определенном вибрационном воздействии всегда приводит к перестроению ячеистых режимов в течение с большей длиной волны, так что, в конечном счете, в слое устанавливается одновихревое течение.
2. Изучены свободные колебания капиллярного моста, равновесная форма которого определяется силами поверхностного натяжения и статическим полем тяжести. В соответствии с условиями эксперимента в невязком приближении найдены значения 25 «нижних» уровней спектра собственных частот колебаний капиллярного моста для различных значений управляющих параметров. Результаты данных расчетов, были использованы впоследствии в ходе масштабного экспериментального и теоретического исследования, связанного с попыткой управления с помощью высокочастотных вибраций процессом выращивания кристаллов по методу жидкой зоны.
3. Выполнен расчет магнитного поля и силы Лоренца в расплавленной полупроводниковой жидкой зоне при воздействии внешнего вращающегося магнитного поля. Решение получено для ограниченной твердыми цилиндрическими массивами жидкой* зоны, находящейся в условиях невесомости. Для разных частот вращения произведен расчет термокапиллярной конвекции и индуцируемого вращающимся, магнитным полем течения в азимутальном направлении. Рассмотрены процессы перераспределения примеси в расплаве при наличии движения фронта кристаллизации. В ходе численного моделирования выяснилось, что вращающееся магнитное поле не сильно, но, тем не менее, сглаживает радиальные неоднородности примеси в жидкой зоне и является мощным средствомдля бесконтактного аксиального перемешивания расплавов.
4. Теоретически исследованы надкритические режимы свободной тепловой конвекции в ячейке Хеле — Шоу при подогреве снизу. Построены карты устойчивости конвективных режимов в случае идеально теплопроводных и теплоизолированных широких граней. Для полостей с разным соотношением сторон теоретически обнаружены новые устойчивые и переходные конвективные течения пульсационного типа. Проведено сравнение результатов численного моделирования и экспериментальных данных для полости с теплоизолированными вертикальными гранями. Выявлено хорошее согласие теории и эксперимента.
5. Отдельно исследованы характеристики типичного для ячейки Хеле — Шоу автоколебательного течения в виде четырех перезамыкающихся с течением времени диагональных вихрей. Отметим, что до настоящего времени этот режим специально в деталях не изучался, несмотря на то, что с ним часто приходилось сталкиваться при исследовании тепловой конвекции в ячейке Хеле — Шоу. Показано, что четырехвихревое течение с перезамыканием диагональных вихрей является регулярным в широком диапазоне значений надкритичности. В силу необычайной устойчивости пространственная картина и временное поведение этого режима практически не меняются, даже когда он становится стохастическим. Вследствие необычайной вариативности оно наблюдается в ячейке Хеле — Шоу в самых разных условиях: в полостях с разным соотношением сторон широких граней, в однородных жидкостях и смесях, в случае как высоко теплопроводных, так и теплоизолированных широких граней. Оно может возникать в результате воздействия, как высокочастотных вибраций, так и вибраций конечной частоты, и продолжает существовать при «включении» других осложняющих факторов. В случае сильно неоднородного нагрева в связи с зависимостью вязкости от температуры у данного колебательного четырехвихревого течения нарушается симметрия между верхними и нижними вихрями. Также, при определенных условиях у этого течения может наблюдаться спонтанное нарушение «лево — правой» симметрии.
6. Рассмотрено влияние, оказываемое высокочастотными горизонтальными вибрациями на конвекцию в ячейке Хеле — Шоу, находящейся в однородном гравитационном поле. Численно и аналитически изучены нелинейные режимы вибрационной конвекции в надкритической области. Показано, что горизонтальные вибрации, ориентированные продольно широким граням полости, понижают порог устойчивости квазиравновесия. Для ячейки 2×20×40 найдены области существования однои двухвихре-вых стационарных течений, изучены нестационарные регулярные и хаотические режимы термовибрационной конвекции. Обнаружены новые регулярные и хаотические режимы в ячейке Хеле — Шоу, которые были названы пульсационными течениями. Расчеты показывают, что пульсацион-ные течения являются результатом нелинейного взаимодействия «нижних» мод, которые отвечают за реализацию регулярных надкритических конвективных режимов. Проведено сопоставление результатов расчетов с экспериментальными данными. Выполнен расчет термовибрационной конвекции в ячейке Хеле — Шоу при воздействии вертикальных вибраций. Показано, что вертикальные вибрации повышают порог устойчивости механического равновесия, стабилизируют все наблюдавшиеся в ходе расчетов течения и отодвигают переход к нерегулярным, течениям* в область больших надкритичностей.
7. С помощью прямого численного моделирования проанализировано влияние коротких инерционных внешних сигналов на конвективные течения в ячейке Хеле — Шоу, подогреваемой снизу пульсирующим точечным источником тепла. Произведен расчет температурных полей и формы течений в случае произвольно ориентированных относительно вертикали коротких периодических сигналов. Разработана методика, позволяющая определять частоту, амплитуду и направление внешнего инерционного сигнала. Результаты расчетов предложено использовать при проектировании прибора, способного регистрировать сильные инерционные воздействия. На основе результатов моделирования выполнен подбор рабочей жидкости и найдены оптимальные значения геометрических параметров полости конвективного датчика.
8. Теоретически изучены надкритические конвективные движения бинарной смеси в связанных каналах конечной высоты. Подтверждено, что в отличие от однородных жидкостей, в смеси наблюдается жесткое возбуждение конвекции вне зависимости от знака термодиффузии, имеют место специфические переходные течения и колебательные режимы конвекции. Предложен механизм, объясняющий наблюдаемые явления, подтвержденный теоретическим решением задачи. Разными методами получены амплитудные кривые, а также поля скорости, температуры и концентрации примеси в поперечном сечении каналов, иллюстрирующие конкуренцию термодиффузионного и термогравитационного механизмов конвекции. Проведено сравнение с экспериментальными данными. В частности, из расчетов стационарных режимов следует, что там, где возникает подъемное течение, в среднем имеет место избыток тяжелой примеси, а в канале с опускным течением — недостаток. Таким образом, конвективная петля может использоваться как установка для разделения смесей на компоненты. С одной-стороны, процесс разделения* смесей происходит в динамике, на фоне конвективного переноса, т. е. быстрее, чем в статике, с другой стороны, в отличие от термодиффузионной колонны в связанных каналах отсутствуют встречные гидродинамические потоки, т. е. в области высоких надкритичностей течение дольше остается ламинарным. 9. Численно и аналитически исследовано влияние вертикальных высокочастотных вибраций на тепловую конвекцию бинарных смесей в связанных каналах конечной высоты при подогреве снизу. Показано, что вертикальные вибрации стабилизируют конвективные течения и повышают порог устойчивости механического равновесия. Рассмотрена свободная и термовибрационнаяконвекция бинарной смеси в ячейке Хеле — Шоу. Проведена аналогия между явлениями, наблюдавшимися в связанных каналах и ячейке Хеле — Шоу.