КИЇВСЬКИЙ НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ
ІМЕНІ ТАРАСА ШЕВЧЕНКА
ГУБА АНТОНІНА ОЛЕКСАНДРІВНА
УДК 532.5
СТАЦІОНАРНІ ТА РІВНОМІРНО-ОБЕРТОВІ КОНФІГУРАЦІЇ
ТОЧКОВИХ ВИХОРІВ
01.02.05 - механіка рідини, газу та плазми
Автореферат дисертації на здобуття наукового ступеня
кандидата фізико-математичних наук
Київ – 2008
Дисертацією є рукопис
Робота виконана у Київському національному університеті імені Тараса Шевченка
Науковий керівник: доктор фізико-математичних наук, професор
МЕЛЕШКО В’ячеслав Володимирович,
Київський національний університет імені Тараса Шевченка,
завідувач кафедри теоретичної та прикладної механіки
Офіційні опоненти: доктор фізико-математичних наук, професор
НІКІШОВ Володимир Іванович,
Інститут гідромеханіки НАН України,
заступник директора
кандидат фізико-математичних наук, доцент
ЧЕРНІЙ Дмитро Іванович,
Київський національний університет імені Тараса Шевченка,
доцент кафедри обчислювальної математики
Захист відбудеться 16 квітня 2008 року о 14.30 годині на засіданні
спеціалізованої вченої ради К26.001.21 у Київському національному університеті імені Тараса Шевченка за адресою: 03127, м. Київ, вул. Глушкова, 2, корпус 7, механіко-математичний факультет, ауд. 41.
З дисертацією можна ознайомитись у бібліотеці Київського національного університету імені Тараса Шевченка (01033, м. Київ, вул. Володимирська, 58).
Автореферат розісланий 5 березня 2008 р.
Вчений секретар
спеціалізованої вченої ради
кандидат фізико-математичних наук А.В. Ловейкін
ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ
Актуальність теми. Теорія вихрових рухів представляє собою один з найкращих підходів до розуміння природи хаотизації потоку та початкового розвитку турбулентності. У випадку нев’язкої рідини вихрова динаміка забезпечує фізичний приклад нелінійних гамільтонових систем нескінченної розмірності та представляє постійний інтерес у зв’язку з дослідженнями хаотичних властивостей динамічних систем. Динаміка вихрових структур являється важливим розділом фізики рідини, газу та плазми, так як всі реальні течії являються вихровими. Особлива увага дослідників традиційно приділена двовимірним вихровим структурам. Зменшення розмірності задачі дозволяє суттєво спростити дослідження, виявити основні властивості та закономірності взаємодії масштабних вихрових структур та з меншими зусиллями досягти необхідних результатів.
Задачі про пошук конфігурацій точкових вихорів, зокрема стаціонарних та рівномірно-обертових конфігурацій вихорів однакової інтенсивності, крім механіки рідини, мають важливе значення в області небесної механіки, фізики надтекучого гелію II. Вивчення руху конфігурацій з невеликою кількістю вихорів поблизу границь найпростішої форми (прямолінійної або колової) дає уявлення про вплив геометрично більш складних границь на природу порядку та хаосу. Дослідження еволюції точкових вихорів, що розміщені по концентричних колах має велике значення для аналізу вихрових доріжок Кармана, що в свою чергу дозволяє вивчати процеси вихроутворення за тілами зі слабкими обтікаючими властивостями. Рівномірно-обертові моделі точкових вихорів застосовуються при пошуку та аналізі стійких конфігурацій вихрових структур, що впливають на формування атмосферних циклонів і океанографічних течій, при проектуванні апаратури в хімічній промисловості.
Задача про формування динамічними системами із скінченою кількістю степенів вільності стаціонарних та рівномірно-обертових конфігурацій належить до найбільш цікавих задач механіки рідини, газу та плазми. Пошук таких конфігурацій та аналіз їх стійкості потребує глибокого розуміння особливостей динаміки руху систем, що розглядаються, й суттєвого звуження діапазонів параметрів багатомірних систем, розробки нових алгоритмічних підходів чисельно-аналітичного розв’язання задачі, що й обумовлює актуальність даної дисертації.
Зв’язок роботи з науковими програмами, планами, темами. Дослідження, що складають дисертаційну роботу, виконані у відповідності до держбюджетної теми “Механіка рухомих деформівних середовищ та експериментальні методи механіки і низькочастотного електромагнітного зв’язку телесистем для похилого і горизонтального буріння нафтогазових свердловин” (№ 01БФ038-02), яка виконувалась у Київському національному університеті імені Тараса Шевченка (2001-2005), а також у рамках проекту INTAS 04-80-7297 “Vortex Dynamics” (2005-2007).
Мета і завдання дослідження. Метою представлених досліджень являється пошук нових стаціонарних та рівномірно-обертових симетричних та несиметричних конфігурацій систем точкових вихорів в ідеальній нестисливій рідині на необмеженій площині.
Для досягнення означеної мети ставились наступні завдання:
- проаналізувати особливості та основні закономірності взаємодії систем точкових вихорів на необмеженій площині;
- розробити ефективний чисельно-аналітичний метод для пошуку нових стаціонарних та рівномірно-обертових конфігурацій точкових вихорів однакової інтенсивності;
- провести класифікацію отриманих симетричних та несиметричних стаціонарних та рівномірно-обертових конфігурацій точкових вихорів однакової інтенсивності відповідно до даних, що містяться в сучасній літературі;
- проаналізувати стійкість знайдених конфігурацій точкових вихорів однакової інтенсивності.
Об’єктом дослідження являються стаціонарні та рівномірно-обертові конфігурації точкових вихорів в ідеальній нестисливій рідині на необмеженій площині.
Предмет дослідження – вплив кількості точкових вихорів однакової інтенсивності на розташування та стійкість стаціонарних та рівномірно-обертових конфігурацій та проведення аналогового експерименту.
Метод дослідження базується на розв’язанні нелінійної алгебраїчної системи рівнянь, що описує рух точкових вихорів однакової інтенсивності в ідеальній нестисливій рідині на необмеженій площині. В якості початкового наближення вибирається стаціонарна конфігурація вихорів порядку та стаціонарна точка потоку рідини. В стаціонарній точці рідини розміщується точковий вихор, інтенсивність якого, по мірі проведення ітерацій, поступово збільшується від нуля до інтенсивності решти вихорів. На кожному ітераційному кроці розв’язується нелінійна система алгебраїчних рівнянь порядку за допомогою чисельного методу Ньютона-Рафсона.
Аналіз стійкості знайдених конфігурацій проводився шляхом чисельного інтегрування (метод Рунге-Кутта 4 порядку) гамільтонової системи руху точкових вихорів.
Наукова новизна одержаних результатів. Спираючись на основні закономірності динаміки точкових вихорів на необмеженій площині та на основні властивості рівнянь руху в гамільтоновій формі, було отримано наступні нові наукові результати:
- представлено новий чисельно-аналітичний метод для пошуку рівномірно-обертових конфігурацій систем точкових вихорів однакової інтенсивності на необмеженій площині, який базується на аналізі стаціонарних точок поля швидкості в рівномірно-обертовій конфігурації системи з меншою кількістю вихорів;
- доповнено відомий, так званий, “Лос-Аламоський каталог” стійких вихрових конфігурацій (структури 83, 91 та 101 в наведених позначеннях), які розміщені по концентричних колах при ;
- знайдено послідовність нових несиметричних рівномірно-обертових конфігурацій систем точкових вихорів однакової інтенсивності. Показано, що такі структури виникають у системах, починаючи з 5 точкових вихорів, а не 8 вихорів;
- показано, що загальна кількість несиметричних конфігурацій збільшується по мірі збільшення кількості точкових вихорів в системах, що розглядаються;
- проведено чисельний аналіз стійкості вихрових конфігурацій, дослідження показали, що більшість з симетричних конфігурацій являються стійкими, тоді як несиметричні конфігурації являються нестійкими відносно малих збурень початкових координат.
Практичне значення одержаних результатів. Практичне значення досліджень визначається широким колом вказаних вище застосувань задачі, наявністю комплексу програм та алгоритмів, які досліджують явище, а також великою кількістю отриманих чисельних результатів. Матеріали дисертації використовуються в спеціальному курсі „Динаміка концентрованих вихорів”, що читається на механіко-математичному факультеті КНУ імені Тараса Шевченка.
Апробація результатів роботи. Основні результати по темі дисертації доповідались та обговорювались на наступних конференціях:
IV Міжнародна конференція “Актуальні проблеми механіки деформівного твердого тіла” (м.Донецьк, червень 2006);
Міжнародна конференція “Dynamical system modeling and stability investigation (DSMSI)” (м.Київ, травень 2007);
Міжнародна конференція “Euler Equations: 250 Years On (EE 250)” (France, Aussois, June 2007);
Міжнародна науково-технічна конференція пам’яті акад. В.І. Моссаковського (1919-2006) “Актуальні проблеми механіки суцільного середовища і міцності конструкцій” (м.Дніпропетровськ, жовтень 2007).
У повному обсязі робота доповідалась на:
семінарі “Сучасні проблеми механіки” Київського національного університету імені Тараса Шевченка під керівництвом член.-кор. НАН України Улітка А.Ф. та проф. Мелешка В.В. (м. Київ, грудень 2007);
Республіканському семінарі при Інституті гідромеханіки НАН України під керівництвом академіка НАН України Грінченка В.Т. (м.Київ, січень 2008).
Публікації та особистий внесок здобувача. Результати дисертаційної роботи викладені у 7 працях, з них 3 опубліковані у рецензованих наукових журналах із переліку фахових видань ВАК України для здобуття наукового ступеня кандидата фізико-математичних наук зі спеціальності 01.02.05 – механіка рідини, газу та плазми [1,2,3] та 4 тез доповідей наукових конференцій [4,5,6,7].
Всі результати роботи отримані автором самостійно. У роботах, виконаних у співавторстві, теоретичні дослідження та чисельні розрахунки виконані здобувачем. Постановка задач належить науковому керівнику Мелешко В.В. та Х. Арефу – співавтору по роботі [3]. Обговорення отриманих результатів виконані спільно з усіма співавторами.
Структура та обсяг дисертації. Дисертаційна робота складається зі змісту, вступу, 5 розділів, висновків, списку використаних джерел. Викладена на 127 сторінках, із них 102 сторінки основного тексту, 20 рисунків, 6 таблиць, бібліографічні посилання складено з 148 джерел.
ОСНОВНИЙ ЗМІСТ РОБОТИ
У вступі шляхом критичного аналізу та порівняння з відомими розв’язаннями наведеної задачі подано інформацію про актуальність досліджень, сформульовано мету та наукову новизну роботи, представлено питання, що виносяться на захист, а також практичне значення одержаних результатів.
У першому розділі представлено узагальнений огляд літератури по темі дисертації, проаналізовано процес виникнення основних напрямків досліджень у даній області. Зазначені невирішені проблеми, аналіз яких дозволив сформулювати основні завдання дисертації.
Відзначено великий внесок у дослідження в області вихрової динаміки робіт таких науковців: Г.Гельмгольц, фундаментальна стаття якого поклали початок теорії вихорів (1858), Е.Г. Борд, А.В. Борисов, Н.С. Васильєв, Д.Н. Горячев, В.Ф. Козлов, Л.Г. Куракін, В.В. Мелешко, Е.А. Новіков, М.А. Соколовський, Л.Г. Хазін, В.І. Юдович, H. Aref, S. Boatto, L.J. Campbell, D. Crowdy, D. Dritschel, W. Grebli, T.H. Havelock, H. Kirchhoff, J.Marshall, P.K. Newton, K.O`Neil, H. Poincare, P.G. Saffman, M.M. Sano, G.J.F.van Heijst, J.J. Thomson, W. Thomson, H. Villat та багатьох інших.
У другому розділі визначено основні поняття теорії вихрових рухів на необмеженій площині, виписано загальні рівняння, що описують рух систем точкових вихорів інтенсивності на необмеженій площині в ідеальній нев’язкій рідині. Вихорі розміщені на комплексній площині та мають координати, при відсутності зовнішнього потенціального потоку рівняння руху точкових вихорів на необмежені рідині.
Виписано перші інтеграли рівнянь руху систем точкових вихорів, що пов’язані безпосередньо з незалежністю функцій від часу та її інваріантністю відносно паралельного переносу та повороту координат:
Представлено основні властивості систем вихорів.
Властивість 1. Якщо для даної конфігурації вихорів одночасно змінити знаки всіх інтенсивностей, то в усі наступні моменти часу система буде проходити через ті ж конфігурації, через які пройшла до цього моменту.
Властивість 2. Нехай в момент часу існує конфігурація, в якій всі вихорі колінеарні, тобто лежать на одній прямій. Тоді конфігурації в момент часу для будь-яких являються відображенням одна одної відносно цієї прямої.
Властивість 3. Система вихорів не може проходити більше, ніж через дві колінеарні конфігурації. Час, необхідний для переходу з однієї колінеарної конфігурації в іншу, завжди однаковий в процесі руху.
Властивість 4. Якщо дві системи та складаються з однакової кількості вихорів, до того ж інтенсивність кожного вихору першої систем пропорційна множнику інтенсивності другої системи , то початкові положення обох конфігурацій подібні. При цьому, масштаби довжин системи рівні масштабам довжин системи, помноженим на множник. Тоді наступна конфігурація другої системи через проміжок часу буде подібна до конфігурації першої системи через час. Час та пов’язаний співвідношенням.
Властивість 5. Якщо інтенсивності вихорів мають однаковий знак, то взаємна відстань між вихорами обмежена під час всього руху.
Властивість 6. При умові фазові потоки системи на площині порядку інтегралів траєкторно-еквівалентні, якщо постійні інтегралів пов’язані співвідношенням.
Наведено означення стаціонарних конфігурацій. До стаціонарних конфігурацій відносяться:
положення рівноваги, при яких швидкість кожного вихору конфігурації дорівнює нулеві;
твердотільні трансляції, при яких вихорі рухаються з постійною швидкістю;
відносні рівноваги, при яких вихорі обертаються навколо деякого центра (центра завихрення) як тверде тіло з постійною кутовою швидкістю;
конфігурації колапсу, при яких всі вихорі прямують до центру завихрення (або віддаляються від нього), а відстані між ними зменшуються таким чином, що конфігурація залишається геометрично подібною до початкової.
В представленій дисертаційній роботі під стаціонарними конфігураціями розуміємо лише відносні рівноваги точкових вихорів, а саме поняття стаціонарного руху розглядаємо лише для гамільтонових систем. Тому стаціонарними вважаємо ті конфігурації, які рівномірно обертаються навколо центра завихрення з постійною кутовою швидкістю.
Функція току рідини, зумовлена конфігурацією точкових вихорів інтенсивності , визначається з виразу:
Поле швидкості потоку рідини пов’язане з полем функції току рівняннями:
Також проведено аналогію з задачами небесної механіки та представлено теореми про кількість деяких конфігурацій, зокрема колінеарних.
Третій розділ дисертації присвячено пошуку стаціонарних конфігурацій рівномірно намагнічених сталевих голок, що рухаються в рідині. Даний експеримент вперше було відтворено Майєром (1878), і саме завдяки цьому експерименту набула розвитку одна з найцікавіших задач вихрової динаміки про рівномірне обертання систем точкових вихорів.
Слід відмітити, що вперше на аналогію між обертанням вихрових структур та електромагнітними явищами звернув увагу Гельмгольц в своїй визначній роботі, яка поклала початок теорії вихрових рухів (1858). І хоча повної аналогії між плаваючими магнітами та точковими вихорами не існує, все ж таки, результати, знайдені експериментально, мають певну аналогію з чисельними та аналітичними розрахунками. В розділі описано експериментальну установку та представлено основні результати експерименту.
В експерименті декілька намагнічених сталевих голок (експеримент проводився з кількістю голок від 2 до 12) розміщувались в корки та занурювались у посудину з водою діаметром 14-16 см та висотою 12 см. Голки розміщувались таким чином, що всі їх додатні полюси знаходились вище поверхні води. До посудини з намагніченими голками зверху над поверхнею площини води підводився сильний магніт протилежним з голками полюсом. Внаслідок чого додатні полюси намагнічених голок відштовхувались одна від одної з силами, величини яких змінювались обернено пропорційно квадрату відстаней між ними, та утворювали через деякий проміжок часу (1 - 1,2 хв.) нерухомі конфігурації. Коли кількість намагнічених голок менше чотирьох, вони розміщувались у вершинах правильних трикутника та чотирикутника.
Отримані конфігурації намагнічених точкових голок при кількості плаваючих голок від 5 до 10 представлені на рис.1. Видно, що 5 та 6 магнітів можуть знаходиться в вершинах правильних многокутників, як без магніту в центрі многокутника, так і з розміщеним в центрі многокутника магнітом, схематично це можна представити у вигляді 4+1 та 5+1. Конфігурація 7 магнітів розміщується, відповідно до запропонованої схеми, у вигляді 7=6+1. Аналогічно, 8=7+1 та 8=6+2; для 9 магнітів 9=8+1 та 9=7+2; 10 магнітів розташовуються у вигляді 10=7+3 та 10=8+2.
Під час відтворення експерименту були також помічені конфігурації намагнічених голок, що розміщувались несиметричним чином, або симетрично відносно вісі, і перебували у даному положенні деякий час нерухомо (приблизно 10-15 сек.). Потім вони продовжували рухатись та приймали стійке симетричне положення. Деякі з цих „нестійких конфігурацій” зображено на рис.2. З рисунку видно, що помічені положення вихорів мають певну аналогію з чисельно отриманими та представленими в роботі несиметричними конфігураціями точкових вихорів. Зокрема, конфігурація „тимчасово стійких” намагнічених голок (а) аналогічна до чисельно отриманої несиметричної конфігурації семи точкових вихорів 71 (лише повернуту на кут ), а конфігурація (с) – схожа на конфігурацію 91 - несиметричну конфігурацію 9 вихорів. Зрозуміло, що отримані конфігурації є нестійкими, але помічені аналогії з чисельними розрахунками ще раз підтверджують важливість експерименту Майєра для вихрової динаміки.
У четвертому розділі дисертації розглянуто рівномірно-обертові конфігурації точкових вихорів однакової інтенсивності. У цьому випадку рівняння руху вихорів мають вигляд.
Розглядається випадок, коли система точкових вихорів рівномірно обертається навколо центру завихрення з постійною кутовою швидкістю , тому розв’язок системи рівнянь.
Функція току рідини зумовлена конфігурацією точкових вихорів, що обертаються з постійною кутовою швидкістю, приймає вигляд:
В даному розділі запропоновано новий метод для побудови рівномірно-обертових конфігурацій точкових вихорів однакової інтенсивності, його чисельна реалізація суттєво спрощує обсяг обчислень при знаходженні нових конфігурацій систем точкових вихорів. Метод базується на розв’язанні нелінійної алгебраїчної системи рівнянь руху точкових вихорів. В якості початкового наближення вибирається стаціонарна конфігурація точкових вихорів та стаціонарна точка потоку рідини в системі координат, що обертається. В стаціонарній точці рідини розміщується точковий вихор, інтенсивність якого, по мірі проведення ітерацій, поступово збільшується від нуля до інтенсивності решти вихорів. При цьому поле швидкостей не зазнає змін. На кожному ітераційному кроці розв’язується нелінійна система алгебраїчних рівнянь порядку , в результаті визначається нова рівномірно-обертова конфігурація з вихором.
Таким чином, для переходу від рівномірно-обертової конфігурації точкових вихорів до системи, що складається з вихору, використовується наступна схема:
1) Вибираємо в якості вихідної системи довільну конфігурацію рівномірно-обертових точкових вихорів.
2) Визначаємо кутову швидкість обертання системи точкових вихорів.
3) Знаходимо в потоці рідини, що розглядається, точки, в яких наведена швидкість збоку вихорів та зовнішнього потоку дорівнює нулеві або має екстремальні значення.
4) В нерухомих точках рідини розміщуємо додатковий точковий вихор інтенсивності та поступово з певним кроком збільшуємо інтенсивність цього вихору до 1. На кожному кроці розв’язується система нелінійних рівнянь відносно невідомих.
5) При розв’язок системи дає нове положення рівномірно-обертової конфігурації системи точкового вихору.
Аналіз системи рівнянь та її розв’язків показує, що нова конфігурація вихорів визначається, в першу чергу, початковим наближенням (конфігурацією точкових вихорів) та вибором точки, в якій розміщується додатковий вихор змінної інтенсивності. Дослідження показали, що стаціонарних точок в рідині може бути декілька, однак не кожна з них може привести до нової конфігурації точкових вихорів. Слід відмітити, що при розв’язанні системи нелінійних рівнянь точковий вихор, спочатку розміщений в різні стаціонарні точки потоку рідини, може потрапляти в одні й ті самі рівномірно-обертові конфігурації систем точкових вихорів.
При чисельній реалізації методу інтенсивність вихору збільшувалась дискретно на кожному кроці розв’язання задачі. При розв’язанні системи рівнянь застосовувався метод Ньютона-Рафсона.
В результаті досліджень сформовано доповнений аналог (фрагмент при), так званого, „Лос-Аламоського каталогу”, який вважається найбільш повним зібранням стійких вихрових конфігурацій для , що розміщуються на вкладених одне в одне концентричних колах (рис.4). Отриманий фрагмент каталогу, принаймні при , відрізняється від вищевказаного наявністю трьох нових конфігурацій 83, 91 та 101. В роботі класифіковано отримані вихрові структури на правильні, полігональні та розміщені по концентричних колах, а також проведено порівняльний аналіз з класами рівномірно-обертових конфігурацій точкових вихорів, наведеними в літературі.
Метод знаходження рівномірно-обертових вихрових конфігурацій дозволив також знайти несиметричні вихрові структури, які виникають. Їх кількість збільшується зі збільшенням кількості вихорів (при та по одній конфігурації, при по три конфігурації, при - сім конфігурацій та при - дев’ять конфігурацій). Побудовано каталог несиметричних конфігурацій. Більшість з представлених несиметричних вихрових структур являються новими. Також в роботі побудовано лінії току всіх знайдених рівномірно-обертових конфігурацій точкових вихорів однакової інтенсивності.
Таблиця 1. Значення інваріантів руху симетричних конфігурацій
№ | Конфігурація | H1 | I1 | P1 | Q1 | H2 | I2 | P2 | Q2 |
1 | 31 | -0,131598 | 2,999912 | 0 | 0 | -0,12812 | 2,91456 | 0 | 0 |
2 | 41 | -0,318666 | 6,0025 | 0 | 0 | -0,31376 | 5,88125 | 0 | 0 |
3 | 51 | -0,598126 | 10,0004 | 0 | 0 | -0,59516 | 9,92763 | 0 | 0 |
4 | 52 | -0,587235 | 9,99824 | 0 | 0 | -0,58089 | 9,84139 | 0 | 0 |
5 | 61 | -0,978449 | 15,0088 | 0 | 0 | -0,97531 | 14,9318 | 0 | 0 |
6 | 62 | -0,979375 | 15,0005 | 0 | 0 | -0,97286 | 14,839 | 0 | 0 |
7 | 71 | -1,47966 | 20,9971 | 0 | 0 | -1,47326 | 20,8384 | 0 | 0 |
8 | 81 | -2,09398 | 28,0026 | 0 | 0 | -2,08775 | 27,848 | 0 | 0 |
9 | 82 | -2,08098 | 27,9975 | 0 | 0 | -2,07392 | 27,822 | 0 | 0 |
10 |