автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическое моделирование волн на поверхности двухфазной среды

ВВЕДЕНИЕ.

Глава 1. МОДЕЛИРОВАНИЕ ПОВЕРХНОСТНЫХ ВОЛН.

1Л. Основные уравнения.

1.2. Граничные условия.

1.3. Нелинейная краевая задача.

Глава 2. ЛИНЕЙНАЯ ЗАДАЧА О

ПЛОСКИХ ПОВЕРХНОСТНЫХ ВОЛНАХ.

2.1. Постановка и решение задачи.

2.2. Фазовая скорость, декремент затухания и амплитуда волны.

2.3. Расчеты для конкретных сред.

2.4. Траектории частиц несущей и дисперсной фазы.

2.5. Другой способ решения линейной задачи.

Глава 3. НЕЛИНЕЙНЫЕ ВОЛНЫ НА

ПОВЕРХНОСТИ ДВУХФАЗНОЙ СРЕДЫ.

3.1. Постановка нелинейной задачи.

3.2. Решение задачи.

3.3. Возмущение концентрации дисперсной фазы и форма свободной поверхности.

Глава 4. НЕЛИНЕЙНАЯ ЗАДАЧА О ПРОСТРАНСТВЕННЫХ

ВОЛНАХ НА СЛОЕ ДВУХФАЗНОЙ СРЕДЫ.

4.1. Постановка задачи.

4.2. Решение линейной задачи.

4.3. Определение второго приближения по малому параметру.

4.4. Основные параметры волны, возмущение концентрации и форма свободной поверхности.

Глава 5. ВОЛНЫ НА СВОБОДНОЙ ПОВЕРХНОСТИ СМЕСИ С

НЕОДНОРОДНОЙ КОНЦЕНТРАЦИЕЙ ДИСПЕРСНОЙ ФАЗЫ

5.1. Основные уравнения.

5.2. Линейная задача о плоских волнах.

5.3. Асимптотическое решение линейной задачи.

5.4. Расчеты для конкретных сред.

Диссертация посвящена исследованию распространения волн на свободной поверхности двухфазной среды.

Многофазные и, в частности, двухфазные среды широко распространены в природе, встречаются в технологических процессах. Они характеризуются наличием макроскопических (по отношению к молекулярным масштабам) включений. Самыми простыми являются дисперсные смеси, которые состоят из двух фаз. Математическое описание таких сред по сравнению с однофазными осложнено в связи с необходимостью учета структуры фаз, эффектов межфазного взаимодействия. В результате увеличивается число параметров, которые следует учесть для получения представления о движении многофазной системы. Поэтому основной проблемой моделирования смесей является получение замкнутой системы уравнений. К числу первых работ в этой области относятся [44], [48], [22], [55]. В статье X. А. Рахматулина [44] впервые была получена замкнутая модель движения многофазной среды. В качестве замыкающего условия было принято предположение о равенстве давлений в фазах, широко применяемое в дальнейшем. С данной целью также часто используются эмпирические соотношения. Так, в работе [23] предложена двухжидкостная модель одномерного пузырькового потока, использующая для замыкания системы уравнений соотношение Зубера - Финдлея.

Моделирование межфазного взаимодействия требует изучения процессов, происходящих вблизи отдельных частиц, анализа закономерностей их движения, столкновений, деформации. Этот вопрос рассмотрен во многих работах, например, [51], [46], [42], [43], [39], [40].

Описание движения многофазной смеси производится с помощью многоскоростных или диффузионных моделей. При многоскоростном описании каждой фазе соответствуют свои макроскопические параметры (давление, плотность, скорость и т. д.), определенные во всей области, занимаемой смесью. Затем, исходя из этих величин, определяются параметры, характеризующие смесь в целом. При этом вводится понятие объемной концентрации - доли объема, занимаемой отдельной фазой [39]. Уравнения динамики многофазной среды обыкновенно получают путем записи физических законов сохранения массы, импульса и энергии применительно к каждой фазе (феноменологический метод) [39], [40], [45], [46], либо путем последовательного осреднения уравнений, описывающих процессы в микромасштабе [48], [40]. Наиболее полно эти вопросы рассмотрены в монографиях Р. И. Нигматулина [39], [40]. Также ведется разработка новых методов. Например, в работе [24] изложены основные принципы построения уравнений для двухфазной смеси на базе общих законов механики и теории обобщенных функций. Уравнения механики многофазных сред рассмотрены и в ряде других публикаций, например, [21], [31], [50], [26], [27] и др.

Эффекты неоднофазности, существенно осложняющие изучение смесей, наиболее полно проявляются при распространении волн. Если ударные и звуковые волны в многофазных и особенно в двухфазных смесях изучены достаточно полно [26], [40], [37] и др., то исследования волн, распространяющихся по поверхности таких сред, почти не представлены в научной литературе.

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

Дайка [17] и Найфэ [35], [36]. В современной постановке первый метод Стокса разработан Ю. 3. Алешковым [1], применительно к магнитогидродинамическим волнам использовался в работах [4], [15]. Метод разложения по малому амплитудному параметру используется и в данной работе. Сущность второго метода Стокса [61 ] состоит в использовании потенциала скорости и функции тока в качестве независимых переменных, что позволяет сделать область, в которой ищется решение задачи, фиксированной. Дальнейшее развитие теория нелинейных поверхностных волн получила в работах А. И. Некрасова, Леви-Чивита, Н. Е. Кочина, Я. И. Секерж-Зеньковича, JI. Н. Сретенского [38], [57], [47]. Последние разработки в данной области касаются волн, распространяющихся по поверхности неоднородных жидкостей [41], [2], [3], [18], [19], [15], [33]. Значительное внимание уделяется корректности постановок краевых задач, а также вопросам устойчивости [20].

В настоящее время теория поверхностных волн для однофазных жидкостей хорошо разработана, но вопрос о распространении волн по поверхности многофазных сред почти не изучен. Однако исследование этого вопроса представляет как практический, так и теоретический интерес. Полученные результаты могут быть использованы при изучении влияния примесей на волновые параметры. Так, в работе [59] указывается на существенное влияние взвешенных примесей на распространение прибрежных волн Атлантического океана в районе Нормандии. Чтобы ответить на вопрос о том, как влияет наличие в жидкости взвешенных частиц или пузырьков на распространение поверхностных волн, необходимо найти зависимости фазовой скорости, частоты, декремента затухания волны, а также формы свободной поверхности от концентрации дисперсной фазы и характеристик межфазного взаимодействия. Близкая к этой задача была рассмотрена в статье [30]. В этой работе в линейном приближении были исследованы стоячие монохроматические волны на поверхности раздела слоев жидкости и смеси этой жидкости с твердыми частицами, численно проанализирована устойчивость этой поверхности. Однако выше указанных зависимостей получено не было. В работе [52] задача гидродинамики для слоя многофазной жидкости при колебании ограждающих конструкций сводится к задаче для однофазной среды с переменной по глубине плотностью. Однако, этот метод не позволяет определить изменение концентрации примесей при волновом движении смеси. В работе [32] рассмотрены бегущие волны на поверхности газожидкостной среды в рамках односкоростной равновесной по давлениям в фазах модели. В статье [58] рассмотрены гравитационные волны в жидкости, содержащей пузырьки газа, однако акцент делается на изучение влияния прохождения волны на траекторию пузырька.

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

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

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

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

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

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

В диссертации принята тройная нумерация формул. Первая цифра указывает номер главы, вторая - номер параграфа, третья - номер формулы.

Основные результаты и выводы диссертации.

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

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

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

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

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

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

7. Траектории частиц несущей и дисперсной фазы в линейном приближении являются разомкнутыми, что объясняется затуханием волны. При малом влиянии сил трения траектории практически являются эллипсами, что соответствует классическим результатам. Дисперсные частицы с меньшей, по сравнению с несущей жидкостью, плотностью (р2 <р°) движутся по траекториям большего размера, чем при р^ > р°.

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

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

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

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

ЗАКЛЮЧЕНИЕ

1. Алешков Ю. 3. Теория волн на поверхности тяжелой жидкости. Л.: Изд-во. Ленингр. ун-та, 1981. 196 с.

2. Алешков Ю. 3. Теория взаимодействия волн с преградами. J1.: Изд-во Ленингр. ун-та, 1990. 372 с.

3. Алешков Ю. 3. Течение и волны в океане. СПб.: Изд-во С.-Петерб. ун-та, 1996. 228 с.

4. Алешков Ю. 3., Баринов В. А., Тактаров Н. Г. О распространении нелинейных магнитогидродинамических поверхностных волн // Магнитная гидродинамика. 1989. №4. С. 79-86.

5. Баринов В. А. Решение уравнений плоского установившегося движения жидкости // Ж. вычисл. матем. и матем. физ. 2001. Т.41. №1. С. 130-134.

6. Баринов В. А., Бутакова Н. Н. Волны на поверхности слоя дисперсной жидкости // Математическое и информационное моделирование. Тюмень: Изд-во Тюмен. ун-та, 2000. С. 57- 63.

7. Баринов В. А., Бутакова Н. Н. Исследование распространения волн по свободной поверхности двухфазной жидкой смеси // Вестн. ТюмГУ. 2001. №2. С.182-190.

8. Баринов В. А., Бутакова Н. Н. Нелинейные волны на свободной поверхности дисперсной смеси // Тр. Среденеволж. матем. о-ва. 2002. Т.34. №1. С. 47-53.

9. Баринов В. А., Бутакова Н. Н. Нелинейная задача о волнах на свободной поверхности двухфазной смеси // Математическое и информационное моделирование. Тюмень: Изд-во Тюмен. ун-та, 2002. Вып. 4. С. 133-144.

10. Баринов В. А., Бутакова Н. Н. Нелинейные волны на свободной поверхности дисперсной смеси // Материалы Всероссийской конференции «Теория и приложения задач со свободными границами». Барнаул: Изд-во Алт. унта, 2002. С. 10-12.

11. Баринов В. А., Бутакова Н. Н. Волны на свободной поверхности двухфазной среды // Прикл. мех. и техн. физ. 2002. Т.43. №4. С. 1-9.

12. Баринов В. А., Бутакова Н. Н. Пространственная задача о волнах на свободной поверхности двухфазной смеси // Вестн. ТюмГУ. 2002. №3. (в печати)

13. Баринов В. А., Бутакова Н. Н. Волны на свободной поверхности смеси с неоднородной концентрацией дисперсной фазы // Вестн. ТюмГУ. 2002. №3. (в печати)

14. Баринов В. А., Тактаров Н. Г. Математическое моделирование магнито-гидродинамических поверхностных волн. Саранск: Изд-во Мордов. ун-та, 1991.96 с.

15. Бутакова Н. Н. Поверхностные волны на слое дисперсной жидкости // Материалы международной конференции студентов и аспирантов по фундаментальным наукам « Ломоносов 2000». Москва: Изд-во МГУ, 2000. Вып. 4. С. 323-323.

16. Ван-Дайк М. Методы возмущений в механике жидкости. М.: Мир, 1967.

17. Габов С. А. Задачи динамики стратифицированных жидкостей. М.: Наука, 1986. 288 с.

18. Габов С. А. Введение в теорию нелинейных волн. М.: Изд-во МГУ, 1988. 176 с.

19. Габов С. А. Новые задачи математической теории волн. М.: Наука, 1998. 448 с.

20. Дейч М. Е., Филлипов Г. А. Газодинамика двухфазных сред. М.: Энергоиз-дат, 1981. 472 с.

21. Иорданский С. В., Куликовский А. Г. О движении частиц жидкости, содержащей мелкие частицы // Изв. АН СССР. МЖГ. 1977. №4. С. 12-20.

22. Канцырев Б. Л., Ашбаев А. А. Двухжидкостная гидродинамическая модель пузырькового потока // Прикл. мех. и техн. физ. 2001. Т.42. №6. С. 64-72.

23. Криль С. П., Берман В. П. Основные принципы построения уравнений механики двухфазных гетерогенных сред // Вестн. Нац. Техн. ун-та Украины. «Киев, политехи, ин-т». 1999. Т36. №2. С 420-427.

24. Курош А. Г. Курс высшей алгебры. М.: Наука, 1971. 432 с.

25. Кутателадзе С. С., Накоряков В. Е. Тепломассообмен и волны в газожидкостных системах. Новосибирск: Наука, 1984. 302 с.

26. Лабунцов Д. А., Ягов. В. В. Механика двухфазных систем. М.: Изд-во МЭИ, 2000. 373 с.

27. Лаврентьев М. А., Шабат Б. В. Методы теорий функций комплексного переменного. М.: Наука, 1987. 688 с.

28. Ландау Л. Д., Лифшиц Е. М. Гидродинамика. М.: Наука, 1986. 735 с.

29. Лобов Н. И., Любимов Д. В., Любимова Т. П. Поведение двухслойной системы жидкость взвесь в вибрационном поле// Изв. РАН. МЖГ. 1999. №6. С. 55- 62.

30. Лойцянский Л. Г. Механика жидкости и газа. М.: Наука, 1973. 848 с.

31. Ляпидевский В. Ю. Катящиеся волны в газожидкостной среде // Прикл. мех. и техн. физ. 2002. Т.43. №2. С. 39-43.

32. Ляпидевский В. Ю., Тешуков В. М. Математические модели распространения длинных волн в неоднородной жидкости. Новосибирск: Изд-во СО РАН, 2000. 420 с.

33. Матвеев Н. М. Методы интегрирования обыкновенных дифференциальных уравнений. Л.: Изд-во ЛГУ, 1955.

34. Найфэ А. X. Методы возмущений. М.: Мир, 1976. 456 с.

35. Найфэ А. X. Введение в методы возмущений. М.: Мир, 1984. 536 с.

36. Накоряков В. Е., Покусаев Б. Г., Шрейбер И. Г. Волновая динамика газо- и парожидкостных сред. М.: Энергоатомиздат, 1990. 248 с.

37. Некрасов А. И. Точная теория волн установившегося вида на поверхности тяжелой жидкости. М.: Изд-во АН СССР, 1951. 95 с.

38. Нигматулин Р. И. Основы механики гетерогенных сред. М.: Наука, 1978. 336 с.

39. Нигматулин Р. И. Динамика многофазных сред. Ч. 1. М.: Наука, 1987. 464 с.

40. Овсянников JI. В., Макаренко Н. И., Налимов В. И. и др. Нелинейные проблемы теории поверхностных и внутренних волн. Новосибирск: Наука, 1985. 319 с.

41. Протодьяконов И. О., Люблинская И. Е., Рыжков А. Е. Гидродинамика и массообмен в дисперсных системах жидкость твердое тело. Л.: Химия, 1987. 336 с.

42. Протодьяконов И. О., Люблинская И. Е. Гидродинамика и массообмен в системах газ жидкость. Л.: Наука, 1990. 349 с.

43. Рахматулин X. А. Основы газодинамики взаимопроникающих движений сжимаемых сред //Прикладн. матем. и механика. 1956. Т. 20. №2. С. 184195.

44. Рахматулин X. А. Газовая и волновая динамика. М.: Изд-во МГУ, 1983. 200 с.

45. Coy С. Гидродинамика многофазных систем. М.: Мир, 1971. 536 с.

46. Сретенский Л. Н. Теория волновых движений жидкости. М.: Наука, 1977. 816 с.

47. Телетов С. Г. Вопрос гидродинамики двухфазных систем. I. Уравнения гидродинамики и энергии // Вестн. МГУ. Сер. Математика. 1958. №2. С. 15-27.

48. Тихонов А. Н., Самарский А. А. Уравнения математической физики. М.: Наука, 1977. 736 с.

49. Файзуллаев Д. Ф., Умаров А. И., Шакиров А. А. Гидродинамика одно- и двухфазных сред и ее практическое приложение. Ташкент: Фан, 1980.

50. Фортье А. Механика суспензий. М.: Мир, 1971. 264 с.

51. Шейнин И. С., Заболотная В. А. Движение многофазной жидкости с переменной по глубине плотностью при колебаниях ограждающих конструкций // Известия ВНИИГ им. Б. Е. Веденеева. 1969. Т. 91. С. 223-233.

52. Barinov V. A., Butakova N. N. Waves on the free surface of a two-phase medium // J. of Appl. Mech. and Tech. Phys. 2002. Vol. 43. № 4. pp. 512-518.

53. Botsh B. Hydraulishe Kennwerte fur Nachklarbecken: Definitionen und Ver-gleich mit Augaben des Arbeitsblattes ATV-A B1 // Korrespond/ Abvasser. 1998. №7. P. 1289-1300.

54. Green A. E. Naghdi P. M. On basic equations for mixtures // Quart. J. Mech. and Appl. Math. 1969. V.2 2, Pt. 4. P. 427-438.

55. Iooss G., Rossi M. Nonlinear evolution of the two-dimensionnal RayleighTaylor flow //Exp. J. Vtch., B/fluids. 1989. Vol. 8. №1. P. 1 -22.

56. Levi-Civita J. Determination rigoureuse des ondes permanents d'ampleur finite //Math. Ann. 1925. № 93. P. 264-314.

57. Liu Chun-rong, Zhou Xian-chu. The motion of a bubble in waves // J. of Hydrodynamics. 1997. Ser. B.l. P. 56-62.

58. Louaked M., Saidi A. Pointwise control and particle analysis for parabolic equation// Proc. 7th Intern. Symp. Comput. Fluid Dynamics., Beijing, 1997. P. 228234.

59. Stokes G. G. On the theory of oscillatory waves // Math, and Phys. Papers. 1880. Vol. l.P. 197-229.

60. Stokes G. G. Suplement to a paper on the theory of oscillatory waves // Math, and Phys. Papers. 1880. Vol. 1. P. 314-326.