Движение твердой частицы в нелинейных волнах на поверхности жидкости тема автореферата и диссертации по механике, 01.02.01 ВАК РФ

<34-7/ЧЬ'Ч-ъ

МИНИСТЕРСТВО ОБЩЕГО И ПЮФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ

ЧИТИНСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

На правах рукописи

Верхотуров Анатолий Русланович

ДВИЖЕНИЕ ТВЕРДОЙ ЧАСТИЦЫ

В НЕЛИНЕЙНЫХ ВОЛНАХ НА ПОВЕРХНОСТИ ЖИДКОСТИ

Специальность 01.02.01. - "Теоретическая механика"

Научный руководитель - доктор техн. наук, профессор О.АБаландин

Диссертация на соискание ученой степени кандидата физико-математических наук

Чита-1998

СОДЕРЖАНИЕ

ВВЕДЕНИЕ................................................................................................................................................4

1. Силы, действующие на твердую частицу

в волнах на поверхности жидкости..................................................................................14

1.1. Силы, действующие на твердую частицу

со стороны жидкости .................................................................................................................15

1.2. Волны на поверхности жидкости и их основные характеристики, влияющие на движение твердой

частицы........................................................................................................................................................................20

1.3. Определение сил, действующих на твердую

частицу в волнах............................................................................................................................................31

2. Математическая модель движения твердой частицы

в поверхностных волнах..................................................................................................................34

2.1. Дифференциальные уравнения движения частицы................. 34

2.2. Дифференциальные уравнения движения частицы

с учетом поправки Осеена ............................................................................................................52

3. Частные случаи движения частицы в поверхностных

волнах..............................................................................................................................................................................58

3.1. Дифференциальное уравнение движения частицы

без учета затухания волн....................................................................................................................58

3.2. Исследование уравнения движения методом

фазовой плоскости........................................................................................................................................59

3.2.1. Положения равновесия частицы на поверхности

волны..............................................................................................................................................................................59

3.2.2. Фазовые траектории дифференциального уравнения

движения частицы....................................................................................................................................80

3.3. Движение частицы на поверхности

идеальной жидкости..............................................................................................................................91

3.3.1. Движение частицы без учета вязкого сопротивления и собственного циклического

движения жидкости................................................................................................................................91

3.3.2. Движение частицы на поверхности идеальной

жидкости в общем случае............................................................................................................101

3.4. Численное интегрирование дифференциального

уравнения движения частицы ................................................................................................107

3.4.1. Основные цели и последовательность проведения

расчетов....................................................................................................................................................................107

3.4.2. Результаты численного решения дифференциального

уравнения движения частицы..................................................................................................110

3.4.3. Сравнение численных и аналитических решений ................... 129

4. Экспериментальное исследование движения частиц

в поверхностных волнах................................................................................................................134

4.1. Описание экспериментальной установки............................................................134

4.2. Порядок проведения и результаты эксперимента......................................134

4.3. Сравнение экспериментальных и теоретических

результатов ..........................................................................................................................................................139

ЗАКЛЮЧЕНИЕ................................................................................................................................144

СПИСОК ЛИТЕР АТУРЫ......................................................................................148

ВВЕДЕНИЕ

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

Одним из условий успешного применения и развития современных вибротехнологий явились глубокие исследования в области теории колебаний и виброперемещений, результаты которых отражены в работах И.И.Блехмана, Г.Ю.Джанелидзе, И.И.Артоболевского, И.И.Быховского, И.Ф.Гончаревича, В.О.Кононенко, Я.Г.Пановко, Р.Ф.Нагаева,

Р.Ф.Ганиева, Л.Е.Украинского, К.В.Фролова и других известных отечественных и зарубежных ученых.

Среди различных направлений в теории вибрационного движения важное место занимают исследования по динамике твердых тел, частиц и взвесей в жидкой среде в волновых или вибрационных полях. Современные исследования и обзор результатов в этой области наиболее полно были даны в работах [26,30] и других. Эти исследования обусловлены необходимостью решения ряда актуальных задач, связанных с проблемами рационального природопользования , защиты окружающей среды, создания современных комбинированных технологий добычи полезных ископаемых и разработки техногенных месторождений и другими задачами.

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

В представляемой диссертационной работе исследуется движение твердого плавающего тела ( частицы ) в однородной жидкой среде под дей-

ствием бегущих регулярных поверхностных волн, распространяющихся от источника вибрационного типа.

Среди известных моделей задачи наиболее простым является случай движения тела в идеальной несжимаемой жидкости. При этом на частицу, движущуюся в неограниченном объёме жидкости, при условии существования потенциала скорости, действует сила, обусловленная эффектом присоединённых масс жидкости [57].

При движении сферической частицы в вязкой несжимаемой жидкости сила вязкого сопротивления определяется известной формулой Стокса, полученной им при помощи линеаризации уравнений движения тела. Такое упрощение возможно при медленном стационарном обтекании сферы при малых числах Рейнольдса, когда силы вязкого сопротивления превосходят по величине инерционные силы [30,57].

Впоследствии результаты, полученные Стоксом, были уточнены Осееном и обобщены Буссинеском на случай произвольного неравномерного движения сферы. В результате сила сопротивления жидкости была представлена в виде дополнительных поправок к формуле Стокса [30,54,55].

Задача о движении тела в жидкости рассматривалась в трудах Н.Е.Жуковского, С.А.Чаплыгина, А.М.Ляпунова, Н.Е.Кочина, Л.Н.Сретенского и других известных ученых [38,50,51,58,87]. Для их работ характерно применение при различных условиях задачи методов теории функций комплексного переменного. В работах [55,57,86,87] дается вычисление гидродинамических сил и моментов, действующих на полностью погруженное тело.

Новые модели задачи изложены в работах [85,93,112,115]. В статье [115] рассмотрено движение частицы в нестационарном потоке идеальной жидкости. Кроме обычных сил в рассмотрение вводится член, учитывающий неравномерность движения жидкости, выводится уравнение дви-

жения частицы. В работе [112] сила сопротивления сферы рассматривается в зависимости от её ускорения. Принимается, что в идеальной жидкости сила сопротивления пропорциональна квадрату ускорения.

Ряд экспериментальных исследований выполнен в последнее время учеными СССР, России, Украины, Японии и других стран [23, 28, 29, 121, 122, 124]. Силы сопротивления, действующие на круглый цилиндр, колеблющийся в продольном и поперечном направлении в потоке жидкости, были определены в исследовании [124]. Экспериментальное определение силы, действующей на цилиндр со стороны колеблющегося потока жидкости, коэффициента сопротивления и инерционного коэффициента сопротивления проведено в широком диапазоне изменения частоты, амплитуды и скорости потока для цилиндров различных размеров. Для вычисления горизонтальной силы сопротивления применяется формула Морисона [135], представляющая продольную силу в виде суммы инерционной составляющей и силы вязкого сопротивления. В статье [121] приведены результаты экспериментального определения сил сопротивления и инерции полностью погруженной сферы на регулярном и нерегулярном волнении.

Одновременно проводятся и теоретические исследования [39, 79, 85, 93, 122, 136]. В работе [93] рассматриваются колебания цилиндра, плавающего на волнении. Главный вектор и главный момент гидродинамических сил, действующих на тело, определяются интегрированием сил давления по поверхности погруженной части тела.

Динамика вертикального цилиндра, плавающего на регулярном волнении, исследуется в работе [39]. При определении сил, действующих на цилиндр, используется подход, основанный на теории качки корабля А.Н.Крылова. Полные силы сопротивления и возмущающие гидродинамические силы определяются интегралами по погруженной части поверхности цилиндра, которые не выражаются в элементарных функциях и вычисляются с помощью компьютера. Нелинейные уравнения движения ци-

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

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

Задаче о движении твердых плавающих частиц под действием поверхностных волн посвящены исследования [10,43,44]. Для этих работ характерно применение энергетического метода, при котором сила, действующая на частицу со стороны жидкости, определяется через среднюю плотность потока энергии волн [10,43], а также через плотность потока мощности [44]. Полученные уравнения движения при решении приближенными методами показывают, что частица движется в направлении распространения волны, и скорость частицы стремится в пределе к фазовой скорости или к групповой скорости волн.

Определение характеристик волнового поля, в котором движется твердая частица, и задание параметров волн с учетом их нелинейности связано с основными результатами теории волн на поверхности жидкости. Теория поверхностных волн является важным разделом классической гидродинамики и широко представлена в литературе исследованиями Стокса, Рэлея, А.И.Некрасова, Г.Ламба, Н.Е.Кочина, Л.Н.Сретенского и многих других известных ученых [40,53,54,56,86,106]. Волновое движение на поверхности идеальной жидкости, в зависимости от постановки задачи и принимаемых упрощений реальной физической модели, описывается либо линейной, приближенной теорией волн (теория волн малой амплитуды), либо нелинейной теорией (теория волн конечной амплитуды). В обоих случаях исходными являются уравнения Эйлера движения идеальной жидкости совместно с определенными граничными и начальными условиями. Дви-

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

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

Наиболее широко развита приближенная теория волн [50,54,87,106]. Упрощения здесь достигаются с помощью линеаризации уравнений Эйлера в предположении, что скорости движения жидкости невелики и можно пренебречь их квадратами. Линейная теория волн дает бегущую (прогрессивную) волну синусоидального профиля. Из уравнения движения частиц жидкости определяется, что их траектории есть эллипсы, большая полуось которых расположена вдоль горизонтальной оси. Для жидкости бесконечной глубины эти траектории выглядят круговыми. Для всех характеристик волны определяемые результаты являются приближенными, а по вопросу траекторий частиц жидкости и вовсе не соответствуют действительности, как это было показано в нелинейной теории.

Нелинейной теории волн (теории волн конечной амплитуды) посвящено большое количество исследований, обзор которых и основные результаты изложены в [40,53,54,55,56,57,62,63,67,69, 82,86,87,88,92,106] и других работах. Впервые Стоке, решая нелинейную задачу о волнах на поверхности идеальной жидкости [54,87], при условии существования потенциала скорости, показал, что, в отличие от линейной теории, для нелинейных волн характерны такие эффекты, как разомкнутость траекторий частиц жидкости и существование приповерхностного течения жидкости в сторону распространения волн, несимметричность профиля волны, подъем среднего уровня жидкости, увеличение фазовой скорости волны и другие. Усредненная по периоду волны скорость приповерхностного течения, или скорость волнового переноса жидкости, выражается экспоненциальной зависимостью от глубины жидкости и быстро убывает с ростом глубины. Решения этой задачи впоследствии были также даны Рэлеем [54].

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

Нелинейная теория гравитационно-капиллярных волн развивается в работах [82,83,84,119,120]. Определяется профиль волны, фазовая скорость, скорость приповерхностного течения с учетом капиллярных эффектов. Более сложные и неоднозначные зависимости приповерхностного течения от амплитуды волны, вязкости, капиллярности, глубины жидкости выявлены в работах Харрисона, Лонге-Хиггинса, Хогана и других [119,131,132,133]. В работе [131] показано, что возможно совершенно противоположное распределение скорости приповерхностного течения по сравнению с тем, что дает

теория Стокса: при малой глубине жидкости у ее поверхности волновое течение направлено назад, против распространения волн, а у дна, наоборот, по направлению движения волны. Возможен и другой вариант: у поверхности течение направлено в сторону движения волн, в средней части потока против движения волн и у дна - вновь в сторону распространения волн. Величина этого течения может в 2...3 раза превосходить значение скорости, определяемое формулой Стокса, и достигать величины 20...80% от значения фазовой скорости волны.

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

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

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

В последние годы в этом направлении выполнен ряд работ на уровне статей, научных отчетов, экспериментальных исследований. При этом большинство работ выполнены в рамках линейной теории волн (теории волн малой амплитуды) [9,11,43,44 и др.].

Несмотря �