Подмодели сжимаемой жидкости и инвариантно-групповые решения

Гарифуллин, Артур Рафаилевич

Подмодели сжимаемой жидкости и инвариантно-групповые решения тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Гарифуллин, Артур Рафаилевич АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Уфа МЕСТО ЗАЩИТЫ

2009 ГОД ЗАЩИТЫ

01.02.05 КОД ВАК РФ

Диссертация по механике на тему «Подмодели сжимаемой жидкости и инвариантно-групповые решения»

Автореферат диссертации на тему "Подмодели сжимаемой жидкости и инвариантно-групповые решения"

003487Б62

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

ГАРИФУЛЛИН АРТУР РАФАИЛЕВИЧ

ПОДМОДЕЛИ СЖИМАЕМОЙ ЖИДКОСТИ И ИНВАРИАНТНО-ГРУППОВЫЕ РЕШЕНИЯ

01.02.05 — Механика жидкости, газа и плазмы

АВТОРЕФЕРАТ

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

Уфа - 2009

1 О ДЕК 2009

003487662

Работа выполнена в Учреждении Российской академии наук Институте механики Уфимского научного центра РАН

Научный руководитель: Официальные оппоненты:

Ведущая организация:

доктор физико-математических наук, профессор Хабиров Салават Валеевич доктор физико-математических наук, профессор Газизов Рафаил Кавыевич доктор физико-математических наук Рылов Анатолий Игоревич

Учреждение Российской академии наук Институт математики с вычислительным центром Уфимского научного центра РАН

Защита состоится «24 » й£кдфя2009 года в /2С^часов на заседании диссертационного совета Д 212.013.09 в Башкирском государственном университете по адресу: 450074, г. Уфа, ул. Заки Валиди, 32 в аудитории 216 физико-математического корпуса.

С диссертацией можно ознакомиться в библиотеке Башкирского государственного университета.

Автореферат разослан « 20 » МояУ^а 2009 года.

Ученый секретарь диссертационного совета, д.т.н., профессор

Ковалева Л.А.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Отыскание и классификация точных решений уравнений газовой динамики (УГД), которые описывают широкий круг явлений в механике сплошной среды, являются актуальными проблемами. Используя подгруппы преобразований, допускаемые УГД, можно находить подмодели (системы уравнений с меньшим количеством независимых переменных) движений газа и параметрические семейства точных решений. Полный список подмоделей еще не завершен. Важной задачей является построение и исследование новых подмоделей с целью получения точных решений и их интерпретации с точки зрения физики. Основой этих построений служит оптимальная система подалгебр, допускаемых УГД. На подалгебрах можно построить множество подмоделей: инвариантных, частично инвариантных, дифференциально-инвариантных. Построение проводится с помощью инвариантов. Первая группа инвариантов выражается через независимые переменные. Вторая группа инвариантов содержит газодинамические функции. Если часть инвариантов назначить функциями других, то получается представление решения, с помощью которого определяется подмодель. Представление решения характеризуется рангом и дефектом. Ранг — это число независимых переменных в подмодели. Дефект — это число функций общего вида. Когда дефект равен нулю, получаются инвариантные подмодели, когда дефект не равен нулю, получаются частично инвариантные подмодели. Инвариантная подмодель, у которой время есть одна из независимых переменных, приводится к эволюционному каноническому типу. В противном случае инвариантная подмодель приводится к стационарному каноническому типу.

В работе рассмотрен класс новых подмоделей УГД, описывающих движение сжимаемой жидкости, когда давление и температура могут достигать больших значений. Рассмотрены 12 двумерных подалгебр и одна четырехмерная подалгебра. Последняя подалгебра порождает 5 подмоделей, для одной из которых найдены все решения и даны их физические интерпретации. При исследовании подмоделей важно вычислить допускаемую алгебру. В диссертации ставится и частично решена задача групповой классификации гидродинамической системы ранга 2 стационарного типа. Такие системы более

общие, чем инвариантные подмодели ранга 2 УГД и появляются при рассмотрении регулярных частично инвариантных подмоделей.

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

• построение и исследование новых подмоделей сжимаемой жидкости на подалгебрах, взятых из оптимальной системы подалгебр 13-мерной алгебры Ли, которая допускается моделью сжимаемой жидкости;

• нахождение всех точных решений одной переопределенной подмодели, описание движений сжимаемой жидкости, соответствующих этим решениям;

• групповая классификация гидродинамической системы ранга 2 стационарного типа.

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

Научная новизна.

• построены новые подмодели сжимаемой жидкости;

• найдены новые частично инвариантные решения уравнений сжимаемой жидкости и дана их физическая интерпретация;

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

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

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

Разработан новый метод интегрирования переопределенных систем, основанный на использовании новых переменных лагранжево-го типа.

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

Практическая значимость. Полученные точные решения можно использовать в качестве тестовых задач для численных методов.

Результаты работы использовались при составлении отчета о НИР Института механики "УНЦ РАН № 01.200.211712 инв. № 02.2.006 07719 за 2008 г.

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

Апробация работы. Основные результаты, представленные в диссертации, докладывались на следующих конференциях и семинарах:

• 34-ой Региональной молодежной конференции «Проблемы теоретической и прикладной математики» (Екатеринбург, 2003);

• II конкурсе научных работ молодых ученых и аспирантов Уфимского научного центра и АН РБ (Уфа, 2003);

• 35-ой Региональной молодежной конференции «Проблемы теоретической и прикладной математики» (Екатеринбург, 2004);

• Всероссийской молодежной научной школе-конференции «Че-ботаревские чтения по проблемам современного группового анализа и его приложениям в нелинейной механике» (Казань, 2004);

• XX Всероссийской школе-семинаре «Аналитические методы и оптимизация процессов в механике жидкости и газа» (Абрау-Дюрсо, 2004);

• II Всероссийской конференции «Актуальные проблемы прикладной математики и механики» (Абрау-Дюрсо, 2004);

• Всероссийской научной конференции «Современные проблемы физики и математики» (Стерлитамак, 2004);

• III конкурсе научных работ молодых ученых и аспирантов Уфимского научного центра и АН РБ (Уфа, 2004);

• 36-ой Региональной молодежной конференции «Проблемы теоретической и прикладной математики» (Екатеринбург, 2005);

• Российской конференции «Механика и химическая физика сплошных сред» (Бирск, 2007):

• Международной научной конференции «Дифференциальные уравнения и смежные проблемы», посвященной юбилею академиков Ильина В.А. и Моисеева E.H. (Стерлитамак, 2008);

• Научных семинарах Института механики Уфимского научного центра РАН под руководством д.ф.-м.н. С. Ф. Урманчеева, профессора В. Ш. Шагапова.

В 2005 г. работа связанная с результатами второй главы диссертации была удостоена диплома за II место на III конкурсе научных работ молодых ученых и аспирантов Уфимского научного центра и АН РБ.

Публикации. Основные результаты исследований опубликованы в 10 печатных работах, из них 4 входят в список ВАК. Список представлен в конце автореферата.

Объем и структура работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы. Объем диссертации составляет 140 страниц, в том числе 13 рисунков и 18 таблиц. Список литературы состоит из 78 наименований.

КРАТКОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Во введении отражена актуальность проблем, рассматриваемых в диссертации, сформулированы цели исследования, отмечены научная новизна и теоретическая значимость работы. Приведены основные итоги по выполнению научных исследований в рамках программы ПОДМОДЕЛИ школой академика Л.В. Овсянникова. Кратко изложена структура диссертации.

В п.1.1 первой главы выписана модель движения сжимаемой жидкости:

Ш + р-1Ур = 0, . + р <1^5=0, 1)5 = 0, (1)

где I) = дь + и ■ V, V = (дх, ду, д2), и = (и, и, ги) — вектор скорости, р — плотность, р — давление, 5 — энтропия.

Уравнение состояния имеет вид:

р = = + (2)

где В ^ 0, 7 ^ 0,1- постоянные, В-( > 0, д(Б) — произвольная функция энтропии. Последнее уравнение в (1) для энтропии иногда заменяется уравнением для давления Бр+ра2(Иуй = 0, где о2 = /р — квадрат скорости звука.

Показано, что уравнение состояния (2) для специальных функций д(5) согласуется с уравнением состояния, описывающим поведение реальных плотных сред (жидкостей и твердых тел), когда давление и температура могут достигать больших значений. Модель (1), (2) допускает алгебру Ли £13 размерности 13 с базисом, который в декартовой системе координат имеет вид:

Х\ = дх, Х2 = ду, Хз = д2, Хю = дг, Х\$ = др,

Х4 = 1дх + ди, Х5=Ьду + дг„ Xв = tдz + д■ш,

Х7 = удг - гду + иди - тд„, Х8 = гдх - хдг + и)ди - идт, (3)

Хд = хду - удх + иду - ьди, Хц = + хдх + уду + гд2,

2 27 Х\2 = - иди - удь - и)дш--~рдр---рдр,

7—1 7~1

где 7 ф 0,1.

Рассмотрены 12 двумерных подалгебр из оптимальной системы подалгебр алгебры Ли ¿13 когда параметр 7 Ф 0,1 произвольный. В п.1.2 приведены основные формулы построения инвариантных и частично инвариантных подмоделей через инварианты подалгебр. В п. 1.3 вычислены инварианты рассматриваемых двумерных подалгебр и сведены в таблицы. На трех подалгебрах получены инвариантные подмодели ранга 2, которые приведены к стационарному каноническому типу. В п.1.4 показано, что для оставшихся девяти подалгебр получаются регулярные частично инвариантные подмодели ранга 3 дефекта 1. Показано, что частично инвариантные подмодели редуцируются к инвариантным подмоделям на одномерных подалгебрах. Полученные инвариантные подмодели ранга 3 приведены к эволюционному и стационарному каноническим типам.

Во второй главе рассмотрена 4-мерная подалгебра из оптимальной системы подалгебр алгебры Ли Ь13. Операторы базиса подалгебры следующие Х\, Х4, Х12, -Х13• Разыскиваются регулярные частично инвариантные решения ранга 2 дефекта 2 с представлением:

Теорема 1 Регулярные частично инвариантные решения (4) ранга 2 дефекта 2 системы (1), (2) определяются только пятью подмоделями. Для представления решения Ц) с уточнением

получается редуцируемая подмодель с представлением решения ранга 2 дефекта 1 на подалгебре {Х[, Х4, Х12}, состоящая из инвариантной подсистемы

и = и(Ъх), и = г 11)1(2/, г), т = Ь 1и)1(у,г), 2

р = Г^р1(у,г), р = р(г,х).

(4)

и = хи!{у,г)+ и0(Ь,у,г), р = рх(у,2) (5)

Бгщ = щ - и\, Лхг)! + рх гр1у = VI, Ди)1 + р11ри = и>1,

Dipi + 7BpJ{viy + wls) = yPi ~ jBplui, и уравнения для функции uq:

tuot + DIUÛ — -щщ,

где D\ = v\dv + w^dz-

Для представления решения (4) с уточнением

u = t-1{x + K-l{ui(y,z)-e\nt)),

Р = i_1(Pi(2/, z)+Kx- e(ln t + 1)), 7 = -1,

(7)

получается редуцируемая подмодель с представлением инвариантного решения ранга 2 на подалгебре {еХг + КХи, Хц + КХ13}:

Бгщ = е - К2ргХ, - + РГ*^ = «ь

+Р\1Ри =101, £>1Р1 +Рх(щу + ил2) = -2р1, (8)

- -Зр^1^ + ади) =р\ + ВрГ1 ~ ы1>

где е, К ^ 0 — постоянные.

Для представления решения (4) с уточнениями

и^К'1

t t-'Ui

(у, z)p0 - Г1 (Â(t)epo + B(t)K2pô1)

+ХГ1, p = t ■r-1, p = i ^ (p1(i/,z)po + BpJ) + i<ri-1a;+

[7 (S(t) + l) ^pô1 - (Â(t) + l) epo] ,

(9)

7_- 1

и = К-1 (г^и1(у,г)р0 - А{г)ер0 - В(1)К2р, р = г, р = Кх + Г^ (р1(у,г)р0 +Вр1) +

(Ю)

(И)

где е, К ф 0, ро ф 0 — постоянные, получаются соответственно переопределенные подмодели:

7+1

£>1«1 =-тЩ+£, + Р1у = VI, £>1101 +Р1г = мь

7-1

3-7 ^ 27

+ И>1г = -£>1Р1 = -~Р1 - иХ\

7-1 7-1

(12)

З7 — 1

в\и\ = -—111 -ь е, £>1«! +р1„ = vi, £>1ш1 +р1г = к71,

7-1

(13)

2 27

+ V)и =--, £>хР1 =--Р1 - иг,

7—1 7—1

£>И>1 + Пу = "1, £>1101 + Р1г = 101,

2-7 п 27 <14>

Щу + ти =-г, £>1Р1 =--Р1

7—1 7—1

7-1 7-

Замечание 1 В Теореме 1 при £ < 0 в формулах для представлений р, р переменную Ь нужно брать по модулю.

Теорема 2 Подмодель (14) имеет решения только для четырех значений показателя у уравнения состояния (2)

7 = 2: р1 = ы = Ш1 — 0;

7 = 3/2: р! = 0, К^^Асоб^,

где А = у бш 6 4- г сое 5 — х(0> X = х(0 ~~ произвольная функция, т(1) = —усоэ^-Ь гвтй, 5 — постоянная;

7 = 4/3: Р1 = 0, и1=у-Ф(0, -ш^г-Ф^),

где функции Ф, Ф связаны равенством

(г - Ф)Ф' = (у - Ф)Ф';

7-2/3: Р1 = -3(у2 + .г2), VI — —2у, ил = -2г.

В п.3.1 третьей главы рассмотрена гидродинамическая система ранга 2 стационарного типа:

= гшх + уиу + + дРу) = аи,

= + + ВГ1^ + д-»Ру) = а

^ = иДх + ^ + Д(иг + = аЛ = Яа^*, (15)

^ = иРх + уРу + А(Н, Р){их + г»у) = ар,

Рк = июк + ти* =ак,к = 1,..., п,

где фиксированные функции ди[х,у) > 0, дх'{х,у) > 0, Д2 = диду -д2(х,у) > 0; аУ{х,у,и,и,Я,Р,и}1 ,...,юп), V = и,и,Я,Р,1,...,п, — линейные или квадратичные функции по переменным и, V, и>к. Функция А считается произвольной.

Системы вида (15) появляются при рассмотрении регулярных частично инвариантных подмоделей ранга 2 УГД. Решается задача групповой классификации этой системы по функции А. Операторы алгебры Ли, допускаемой системой (15), имеют вид:

х = ?дх + $уду + £иди + + ?°Лд„н + £ндц + £рдр,

где = С(х, у, и, и, Л, Р, го1,.., гоп). В качестве независимых производных выбираются величины ыж, иу, ги^, Ру, Рх, Ру, г = 1,..., п. Условия инвариантности системы (15) относительно оператора X записываются с неопределенными множителями:

Х^-а") = (\»+\1их+\уу+ХЪКу+\»Р1Рх+\»Р2Ру+\^1)(ар-Рр),

где Х1 — продолжение оператора X. Неопределенные множители х", Лу, ] = и, v, Р, Р1, Р2,1,..., п, зависят от переменных я, у, и, и, и/, Л, Р. Зависимые производные выражены через независимые из уравнений системы (15).

Условия инвариантности содержат не более чем квадратичные слагаемые по независимым производным. Приравнивание нулю коэффициентов при различных квадратичных и линейных слагаемых по независимым производным приводит соответственно к 2-уравнениям и 1-уравнениям. Эти уравнения проинтегрированы в п.3.2 и п.3.3 соответственно. Оставшиеся уравнения не содержат независимых производных и образуют О-уравнения, они приведены в п.3.3.

Лемма 1 Общее решение 2-уравнений имеет вид:

Р = А? = 0. (16)

Лемма 2 Интегрирование 1-уравнений дает представление для координат оператора X и для неопределенных множителей А" при Аф 0:

= у, 1,Р)), = - 2а; - 2п),

е = Р£{х,у) + &(*>ь), А" = -иц, А» = -«/х, Ай-Д^,(17)

Хк — 2 А2 еля-1 - тр) - Ар = —ш - £р — Ад,

где ¡1 = 2 Д2 (ш/> - R~l)шJ + шзБр, а> = ш(х,у^\^3(Я,Р)), J = - 2£?ии + 2 Д2 га(Д, Р),

А9Р + Д5д = 0, АтР + Птн = АВ.'1. (18)

Кроме того, 1-уравнения включают равенства для определения координат £,х,

п = $ +9{9у)~1й - 1(9"Г1дд" = £ +9(диГ1% ~ ¿(дТ^д",

А2 ((9иГ^ + (Л-^) =дд-\д {{9иГ1дди + (Л" ЭД, (19) где д = £хдх + £уду, и классификационное соотношение:

Ар(Р£[ + + - 2п) = (20)

При А = 0 в равенствах (17), (18) нужно считать (р = £р{х,у,Р), т = О, Лй = -Щрр + И/г, а 5(Д,Р) заменяется па Р.

О-уравнения имеют вид:

и2£х + 2ич£у + + аи(£ + + 2ш)+

+ ьи3Ъ) + Я~\ди£,р + + и{ишх + ушу)+ +и(иЗх + уЗу)щ + - Хаи - (ар - Аа^)А" = 0;

еще одно уравнение получается из предыдущего уравнения заменой и <—> и, х <—> у\ + шак + + + {иЗх + уЗу + аиЗи + а?

-Хак - (ар - Ааf)Лfc = 0; (21)

и(£р - 2п + £ + фх + «(£? - 2п + й + Щ)у + а?ш-—.Ха^* — ишх — уо1у — (иЗх + V Зу + аи Зи + avЗv)u}J— - (ар - Аа?)Я~1ХК = 0; <р + у£р + А[и{£ + + + +ии;х + +{иЗх + ьЗу + аи Зи + а" J„)шJ + а'ш,^ + а^ (ы + ££)]--Хар - (ар - Аа?)ХР = 0.

В п.3.4 с помощью классификационного соотношения (20) с точностью до преобразований эквивалентности для системы (15) (растяжение по Я, Р, перенос по Р) получено 13 случаев специализации функции А, что согласуется с ранее полученными результатами по задаче групповой классификации для инвариантных подмоделей УГД.

В п.3.4 с помощью полученных формул вычислены алгебры допускаемые подмоделями (6) и (8) из Теоремы 1.

В четвертой главе в силу Теорем 1 и 2 формулируется

Утверждение 1 Система (1), (2) имеет решения вида (11) только для четырех значений показателя 7

7 = 2 : и = 2у, г; = ги = О, р = Ь 2, р = (Вр0 - х2) Ь 5 = (В(р0 - 1) - х2) Г4, а = у/2В\Ц~1-,

(22)

7 = 3/2 : и = 3у, г/ = О, ю = р = Г4,

р = (Вр^2 - Зх2) Г6, 5 = [В(р1/2 - 1) - 3*2) 1

-6 (23)

а = ^/ЗВ/ЭД-1;

,,, .х у - Ф(0 г-Ф(0 7 = 4/3: « = 4-, и= " , го =-

р = Гв, р=(Вр£/3-6аф-8, (24)

5=(в(р^/3- 1) - б®2)) *-8, а — -у/^В/З^!-1, где функции Фг Ф связаны равенством (г — Ф)Ф' = (у — Ф)Ф';

7 = 2/3: и = Д р = 46, р=(вро1/3-Зг2)г4,

Г (25)

я = (¿?(р0-1/3 -1) - зг2), а = у'ВДГ \

где г = |я|, и = —2~; и, V, из — компоненты вектора скорости, 5 — функция энтропии, а — скорость звука.

Замечание 2 Далее в формулах (22)-(25) положено ро = 1; е положено В — 2, а в (23)-(25) В = 7""1 в силу преобразования подобия.

Для каждого из решений (22)-(25) движение рассматривается в области, в которой давление положительно. Предполагается, что на поверхности с нулевым давлением жидкость переходит в другую фазу, которая не может быть описана уравнением состояния (2). Такие поверхности далее называются поверхностями фазового перехода (фазовыми границами). Звуковой поверхностью называется поверхность, на которой скорость частиц жидкости равна скорости звука.

В п.4.1 рассмотрено решение (22). Оно описывает коллапс частиц в точку или мгновенный точечный источник при одномерном неизоэнтропическом движении (рис. 1). Решение обрезается характеристиками и склеивается в непрерывное решение одномерной подмодели в конечной области (рис. 2). Кривая С02 — мировая линия частицы в сжимающейся под действием поршня жидкости. Кривая Coi — мировая линия той же частицы в расширяющейся жидкости. Звуковые характеристики ОВ и О А обрезают области определения точных решений, описывающих коллапс и мгновенный источник. В области BOA доказано существование и единственность непрерывного одномерного движения. Показана возможность непрерывного периодического неизоэнтропического движения сжимаемой жидкости в ограниченной области под действием поршня.

Рис. 1. Характеристики одномерной подмодели сжимаемой жидкости на решении (22)

Рис. 2. Сопряжение сходящегося и расходящегося движений жидкости

В п.4.2 рассмотрено решение (23). Мировые линии частиц в ре-

шении имеют вид х — М3, У = Уо, г — гщ1 + где Ь0, у0, и;0 —

постоянные. Решение рассматривается в области х € (—\/2/3, \/2/3). Плоскости а: = ±\/2/3 есть плоскости фазового перехода. При £ > О решение описывает мгновенный источник, когда частицы при Ь = О вылетают с плоской кривой Ь (рис. 3) заданной уравнениями х = О, У = У о, -г = х{Уо)- Частицы разлетаются либо по кубическим параболам

х = к(г~х{у)) , У = Уо, (26)

где к = 6о/г0д, Ьо Ф 0, №о 0, либо по прямой параллельной оси х при Ьо ф 0, иио = 0, либо по оси параллельной оси г при 6о ^ 0, №о = 0, либо остаются на месте в точке (0, уо, х{Уо)) при Ьо = Юо = 0. Сечение звуковой поверхности Эят2 + (г —

Х(у)) = 1 плоскостью у = уо образует звуковую линию (эллипс), внутри которой частицы жидкости двигаются с дозвуковыми скоростями, вне ее с сверхзвуковыми скоростями. Рассматривается движение в области х € [0, \/2/3), г > х(у) (рис. 3). В области

х 6 (—\/2/3,0] траектории получаются отражением относительно плоскости х = 0. Движение ограничивается стенками: плоскостями у = а < 0, у = Ь > 0, х - 0, где а, Ь — постоянные, и цилиндром Су1 с направляющей кривой Ь и образующей параллельной оси х. Траектории

частиц в области г < х(у) получаются отражением траекторий из области г > х(у) относительно образующей цилиндра. На рис. 3

1(0) +1 г

Рис. 3. Звуковая поверхность и траектории частиц

представлены траектории частиц, вылетающих из точки с параметром Уй — 0. Звуковая линия изображена штрих-пунктирной линией, уровень х — \/2/3 изображен штриховой линией. Частицы, вылетающие из точки с параметром уо ф 0, двигаются по траекториям полученным параллельным переносом картины траекторий с центром (0,0, х(0)) в центр (0, уо,х(уо))- Таким образом, решение при £ > 0 в данной постановке можно интерпретировать как движение сжимаемой жидкости в канале расслоенное по плоскостям параллельным плоскости хОг из мгновенного дозвукового источника через звуковую поверхность к фазовой границе.

В рамках полученного решения можно рассмотреть другие виды каналов. Например, канал, у которого стенка х — 0 заменяется стенкой, образованной параболами (26) с одинаковым значением к.

В п.4.3 рассмотрено решение (25), которое интерпретируется как гладкое расширение (движение при £ < 0) или сжатие (движение при Ь > 0) жидкости под действием сферического поршня. Траектория поршня есть одна из мировых линий решения. Решение рассматривается в области г < 2-1/2. Сферический слой радиуса г = есть поверхность фазового перехода. Внутри шара радиуса г = 2-1 жидкость движется с дозвуковыми скоростями, в сферическом слое 2"1 < г < 2~1^2 со сверхзвуковыми скоростями. Если на поршне возникло сферическое возмущение, то часть возмущения останется на поршне (характеристика Со), другая часть возмущения попадет в жидкость по характеристике С+ или С-. В случае расширяющегося поршня возмущение, попавшее в жидкость, распространяется следующим образом. Внутри шара радиуса г = 2~г возмущение распространяется в направлении к центру шара. В сферическом слое 2-1 < г < 2™1^2 возмущение распространяется в направлении от центра слоя. В слое радиуса г = 2"1 часть возмущения задерживается, другая часть возмущения доходит до слоя радиуса г 2~~1/"2. В случае сжимающегося поршня сферическое возмущение, возникшее на поршне, распространяется по характеристике С_ в направлении к центру. Возмущение при отражении по характеристикам С+ от центра и при отражении по характеристикам С_ от поршня со временем сходится к центру (рис. 4).

В п.4.4 рассмотрено решение (24). Оно рассматривается в обла-

Рис. 4. Распространение сферического возмущения в жидкости сжимаемой сферическим поршнем; 1 — траектория поршня

Рис. 5. Звуковая поверхность и траектории частиц в области х 6 [0, \/2/4)

сти х € (-л/2/4, л/2/4). Плоскости х = ±л/2/4 есть плоскости фазового перехода. При t > 0 решение описывает движение сжимаемой жидкости из мгновенного дозвукового источника через звуковую поверхность к фазовой границе. Возможно два случая многообразия источника: кривая в плоскости х = 0 (рассмотрена окружность единичного радиуса с центром в начале координат) или точка в начале координат. Уравнения мировых линий в первом случае имеют вид х = bot4, у — —Cot sin/о + cos ¿о, г = Cot cos Iq + sin ¿о, где bo, Со ф 0, lo € [0,2тг) — постоянные. Звуковая поверхность есть эллипсоид с полуосями \/2/4, ч/2, \/2. Рассмотрено решение с Со одного знака, чтобы не рассматривать столкновения частиц. Траектории частиц, вылетевших из одной точки, лежат в плоскости у cos Iq + z sin Iq = 1 касательной к цилиндру у2 + z2 — 1 и представляют собой параболы четвертой степени или прямую в плоскости х = 0. На рис. 5 показано движение частиц, вылетевших из точки с параметром ¿о = 0. Образующая d цилиндра есть предельное положение, к которому стремятся траектории частиц, для которых Ь^С^4 —+ оо. Плоскость Р — плоскость фазового перехода. Картина движения частиц, вылетевших из точек с параметром 1о € (0,27г), получается поворотом вокруг оси х на угол ¿о картины движения частиц, вылетевших из точки с параметром 1о = 0. Звуковая поверхность в области х 6 (—v/2/4,0] полу-

чается отражением относительно плоскости х = 0 звуковой поверхности из области х € [0,\/2/4). При ¿о < 0 траектории получаются отражением относительно плоскости х = 0 траекторий построенных при 6о > 0, а при Со < 0 — отражением относительно образующей цилиндра траекторий построенных при Со > 0.

Уравнения мировых линий во втором случае имеют вид х — 6o¿4, у ~ v0t, z = wat, где £>о, vq, Wq — постоянные. Звуковая поверхность есть эллипсоид с полуосями 1/4, 1, 1. Частицы движутся либо по параболам четвертой степени, либо по оси х, либо по радиальным прямым в плоскости х = 0 или остаются в точечном источнике (начале координат). Проекции парабол на плоскость yz есть радиальные линии. Частицы, двигающиеся по параболам или вдоль оси х, достигают плоскостей х — ±>/2/4, на которых происходит фазовый переход. Получены уравнения характеристического коноида, который является геометрическим местом всех бихарактеристик, выходящих из данной точки K(to,xo). На рис. 6 показано сечение плоскостью,

а) б) в)

Рис. 6. Поверхности уровня характеристического коноида при возникновении возмущения в области дозвуковых — а), звуковых — б), сверхзвуковых — в) течений жидкости. 1 — звуковая поверхность, 2 — траектория частиц, 3,4 — поверхности уровня коноида, 5 — плоскости фазового перехода

содержащей траекторию частицы и точку возникновения возмущения (xo,sq)j где 5о — уУо + z'u- ДВУХ поверхностей уровня t = const характеристического коноида. Выяснено, что если возмущение возникло внутри звуковой поверхности, то поверхность уровня коноида полученная в данный момент времени содержит в себе поверх-

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Получены все инвариантные подмодели с 2-мя и 3-мя независимыми переменными в каноническом виде и некоторые регулярные частично инвариантные подмодели сжимаемой жидкости.

2. Найдены все решения одной частично инвариантной переопределенной подмодели сжимаемой жидкости методом обобщенных лагранжевых переменных.

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

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

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

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Подмодели сжимаемой жидкости на двумерных подалгебрах / А. Р. Гарифуллин // Сибирский журнал индустриальной математики. - 2003. — Т. 6. — № 1(13). — с. 16 - 26.

2. Подмодели сжимаемой жидкости на двумерных подалгебрах / А. Р. Гарифуллин // Труды 34-й Региональной молодежной конференции «Проблемы теоретической и прикладной математики». — Екатеринбург, 2003. — С. 101-105.

3. Подмодели ранга 2 стационарного типа уравнений газовой динамики для одной 4-мерной подалгебры / А. Р. Гарифуллин // Труды 35-й Региональной молодежной конференции «Проблемы теоретической и прикладной математики». — Екатеринбург, 2004. - С. 119 - 123.

4. Решения с линейным полем скоростей регулярной переопределенной подмодели сжимаемой жидкости ранга 2 дефекта 1 / А. Р. Гарифуллин // Труды Всероссийской научной конференции "Современные проблемы физики и математики". Т. 1. — Уфа: Изд-во «Гилем», 2004. — С. 32 - 35.

5. Решения регулярной переопределенной подмодели сжимаемой жидкости ранга 2 дефекта 1 / А. Р. Гарифуллин // Материалы всероссийской молодежной научной школы-конференции «Чеботаревские чтения по проблемам современного группового анализа и его приложениям в нелинейной механике». — Казань, 2004. - С. 25 - 31.

6. Групповая классификация гидродинамической системы ранга два стационарного типа / А. Р. Гарифуллин // Сибирский журнал индустриальной математики. — 2004. — Т. 7, — Л'5 3(19). — С. 66 - 75.

7. Общее решение переопределенной подмодели сжимаемой жидкости ранга 2 дефекта 1 / А. Р. Гарифуллин // Труды 36-й Региональной молодежной конференции «Проблемы теоретической и прикладной математики», — Екатеринбург, 2005. — С. 122-125.

8. Пример сферически симметричного движения сжимаемой жидкости / А. Р. Гарифуллин // Сибирский журнал индустриальной математики. — 2007. — Т. 10, — № 2(30). — С. 66-75.

9. Расслоенное по параллельным плоскостям движение сжимаемой жидкости в канале из мгновенного дозвукового источника через звуковую поверхность к фазовой границе / А. Р. Гарифуллин // Труды международной научной конференции «Дифференциальные уравнения и смежные проблемы». Т. 1. — Уфа: Изд-во «Гилем», 2008. — С. 54-58.

10. Непрерывное сопряжение специальных неизэнтропических одномерных движений газа / А. Р. Гарифуллин, С. В. Хабиров // Труды института математики и механики УрО РАН. — Т. 14, — № 1. - 2008. - С. 22 - 30.

Гарнфуллин Артур Рафаилович

ПОДМОДЕЛИ СЖИМАЕМОЙ ЖИДКОСТИ И ИНВАРИАНТНО-ГРУППОВЫЕ РЕШЕНИЯ

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

Лицензия на издательскую деятельность ЛР№ 021319 от 05.01.99г.

Подписано в печать 18.11.2009 г. Формат 60x84/16. Усл. печ. л. 1,32. Уч.-изд. л. 1,38. Тираж 100 экз. Заказ 761.

Редакциопно-издательский центр Башкирского государственного университета 450074, РБ, г. Уфа, ул. Заки Валиди, 32.

Отпечатано на множительном участке Башкирского государственного университета 450074, РБ, г. Уфа, ул. Заки Встиди, 32.

Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Гарифуллин, Артур Рафаилевич

Введение

Глава 1. УРАВНЕНИЯ ПОДМОДЕЛЕЙ СЖИМАЕМОЙ ЖИДКОСТИ НА ДВУМЕРНЫХ ПОДАЛГЕБРАХ

1.1. Уравнения движения сжимаемой жидкости и допускаемая ими алгебра.

1.2. Инвариантные и частично инвариантные подмодели уравнений движения сжимаемой жидкости (общая теория).

1.3. Инвариантные подмодели на двумерных подалгебрах.

1.4. Редукция регулярных частично инвариантных подмоделей к инвариантным

Глава 2. РЕШЕНИЕ ПЕРЕОПРЕДЕЛЕННОЙ ПОДМОДЕЛИ РАНГА 2 ДЕФЕКТА 2 СЖИМАЕМОЙ ЖИДКОСТИ ДЛЯ ОДНОЙ ЧЕТЫРЕХМЕРНОЙ ПОДАЛГЕБРЫ

2.1. Определенные и переопределенные подмодели ранга 2 дефекта 2 сжимаемой жидкости для одной четырехмерной подалгебры и их редукция.

2.2. Решения переопределенной подмодели ранга 2 дефекта 2 в случае линейного поля скоростей

2.3. Решения переопределенной подмодели ранга 2 дефекта 2 в ортогональной системе координат.

2.4. Общее решение переопределенной подмодели ранга 2 дефекта

Глава 3. ГРУППОВАЯ КЛАССИФИКАЦИЯ ГИДРОДИНАМИЧЕСКОЙ СИСТЕМЫ РАНГА 2 СТАЦИОНАРНОГО ТИПА

3.1. Система определяющих уравнений допускаемой алгебры

3.2. Интегрирование 2-уравнений системы определяющих уравнений

3.3. Интегрирование 1- и 0-уравнений системы определяющих уравнений

3.4. Групповая классификация по уравнению состояния.

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

Глава 4. ФИЗИЧЕСКАЯ ИНТЕРПРЕТАЦИЯ РЕШЕНИЙ ПЕРЕОПРЕДЕЛЕННОЙ ПОДМОДЕЛИ РАГА 2 ДЕФЕКТА 2 СЖИ

МАЕМОЙ ЖИДКОСТИ

4.1. Непрерывное сопряжение инвариантного схождения и расширения через характеристики одномерной нестационарной подмодели сжимаемой жидкости (7 = 2).

4.1.1. Постановка задачи о сопряжении сходящегося и расходящегося движений сжимаемой жидкости

4.1.2. Построение решения задачи Гурса.

4.1.3. Доказательство сходимости.

4.1.4. Сопряжение расходящегося и сходящегося движений сжимаемой жидкости

4.2. Расслоенное по параллельным плоскостям трансзвуковое движение сжимаемой жидкости в канале из мгновенного дозвукового источника (7 = 3/2).

4.3. Задача о трансзвуковом движении сжимаемой жидкости под действием сферического поршня с границей перехода в другую фазу (7 = 2/3).

4.3.1. Точное решение и характеристики.

4.3.2. Физическая интерпретация решения, распространение сферического возмущения при t < 0.

4.3.3. Физическая интерпретация решения, распространение сферического возмущения при t > 0.

4.4. Задача о движении сжимаемой жидкости из мгновенного дозвукового источника через звуковую поверхность к фазовой границе (7 = 4/3).

Введение диссертация по механике, на тему "Подмодели сжимаемой жидкости и инвариантно-групповые решения"

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

В современной теоретической физике при исследовании строения материи (микромир) или Вселенной (макрокосмос) теоретико-групповые методы играют основополагающую роль. В механике сплошных сред в этом вопросе наблюдался пробел. Систематические исследования по применению методов группового анализа к моделям механики сплошной среды были начаты JI.B. Овсянниковым и его школой в конце 50-х годов прошлого столетия [21, 22]. В работах JI.B. Овсянникова, Н.Х. Ибрагимова, В.В. Пухначева, А.А. Никольского, J1.B. Капитанского, Ю.Н. Павловского, А.А. Бучнева, В.О. Бытева и других авторов впервые были изучены групповые свойства дифференциальных уравнений механики жидкости и газа и показано как эти свойства можно использовать для решения физически важных задач [2, 3, 7, 8, 9, 18, 19, 39, 40, 41]. В работе [23] JI.B. Овсянниковым была выдвинута программа ПОДМОДЕЛИ. Более детальное описание можно найти в [26]. Здесь же приведены базовые сведения, необходимые для выполнения программы ПОДМОДЕЛИ применительно к уравнениям газовой динамики (УГД) pDu + Vp — 0, Dp + pdivu = 0, DS = 0, (0.1) где D — dt + й • V, V = (<9Ж, ду, dz), и — (и: v, w) — вектор скорости, р — плотность, р — давление, S — энтропия, р = f(p, 3) — уравнение состояния, да = д/да. Последнее уравнение для энтропии можно заменить уравнением для давления Dp + pa2divu = 0, где а2 — fp — квадрат скорости звука. Энтропия S определена с точностью до замены S —> g(S) с произвольной непостоянной функцией д. Это преобразование энтропии является преобразованием эквивалентности для уравнения состояния и всегда допускается УГД. Уравнения (0.1) представляют собой модель невязкого нетеплопроводного газа, который движется в отстутствии внешних источников энергии и силовых полей.

Целью этой программы является исчерпывающее использование симметрии модели для построения классов точных решений (подмоделей) этих уравнений. Промежуточные итоги по выполнению программы для уравнений газовой динамики были опубликованы в работе [34].

Кратко остановимся на основных понятиях группового анализа [22]. Говорят, что система дифференциальных уравнений Е допускает группу G непрерывных преобразований, если Е остается неизменной под действием любого преобразования из G. Фундаментальное свойство допускаемой группы состоит в том, что группа G действует на множестве решений системы Е, переводя любое решение системы снова в решение. Это справедливо и для любой подгруппы Н группы G.

В рамках программы ПОДМОДЕЛИ точные решения УГД строятся с помощью подгрупп Н допускаемых уравнениями группы Галилея, расширенной растяжением. Каждая подгруппа Н допускаемой группы G также допускается исходной моделью Е. Подгруппа Н имеет инварианты, конечные и (или) дифференциальные. Установление дополнительных соотношений между инвариантами подгруппы Н выделяет из множества решений модели Е определенный класс точных решений. Такие решения выражаются через новые искомые функции (инварианты), удовлетворяющие выводимой из Е системе дифференциальных уравнений, которая называется факторсистемой Е/Н. Поэтому фактор система Е/Н называется подмоделью исходной модели Е. Число сг независимых переменных в факторси-стеме называется рангом подмодели. В стандартном случае четырехмерного пространства событий (время и три координаты), на котором определена система Е, ранг подмодели может принимать значения 3, 2, 1, 0.

Существенно различаются два типа подмоделей: инвариантные (ИП) и частично инвариантные (ЧИП). Для ИП все искомые функции выражаются через инварианты, поэтому система уравнений называется определенной. В ЧИП только часть искомых функций имеет инвариантное представление. На оставшиеся «лишние» функции при построении подмодели не накладывается никаких дополнительных ограничений, они считаются зависящими от всех независимых переменных исходной модели. Число S таких функций определяет дефект ЧИП. Уравнения ЧИП содержат инвариантную подсистему для функций, имеющих инвариантное представление, и переопределенную подсистему для «лишних» функций. Исследование переопределенной подсистемы для «лишних» функций на совместность (приведение в инволюцию) часто является не простой задачей, что определяет сложность исследования ЧИП. Таким образом, ЧИП характеризуется парой чисел: сг и 5, задающих соответственно число независимых переменных в подмодели и число искомых функций, не имеющих инвариантное представление. Иногда в процессе приведения переопределенной системы в инволюцию лишние функции уточняются, и получается ЧИП с тем же рангом а, но меньшим дефектом 6' < 5 на подгруппе Н' С Н меньшей размерности. Говорят, что ЧИП (сг, на подгруппе Н редуцировалась к ЧИП (<т, 5') на подгруппе II' [22]. Если при приведении в инволюцию системы все лишние функции принимают инвариантное представление, то есть $ = 0, то данная подмодель совпадает с ИП ранга сг, построенной на подгруппе меньшей размерности Н' С Н. Тогда говорят, что ЧИП (сг, 5) на подгруппе Н редуцировалась к ИП (сг, 0) на подгруппе Н'. Установление редукции позволяет избежать большого количества фактически ненужной работы.

Два решения системы Е называются несущественно различными относительно группы G, если одно из них можно перевести в другое некоторым преобразованием этой группы. Если такого преобразования нет, то два решения называются существенно различными относительно группы G. Рассмотрим множество Н — решений, получаемых относительно всевозможных подгрупп Н С G. Любые два решения из этого множества несущественно различны, если соответствующие им подгруппы сопряжены (подобны) отностельно внутренних автоморфизмов допускаемой группы G. Действие внутренних автоморфизмов разбивает множество подгрупп группы G на классы подобных. Существенно различные решения получаются относительно различных классов подобных подгрупп. Таким образом, задача перечисления всех существенно различных решений сводится к разбиению подгрупп группы G на классы подобных и определению представителей этих классов. Совокупность таких представителей называется оптимальной системой подгрупп и обозначается символом QG. Решения, соответствующие подгруппам из QG, образуют оптимальную систему решений. Все остальные решения можно получить из этой системы с помощью действия группы G.

В теории Ли каждой группе преобразований G ставится в соответствие некоторая алгебра дифференциальных операторов L. Это соответствие является взаимооднозначным и справедливо для подгруппы и подалгебры. Внутренним автоморфизмам группы G соответствуют внутренние автоморфизмы алгебры L, действие которых разбивает множество подалгебр алгебры L на классы подобных. Совокупность представителей этих классов называется оптимальной системой подалгебр и обозначается символом ©L. Оказывается, что задачу о нахождении оптимальной системы подгрупп QG удобнее решать как задачу построения оптимальной системы подалгебр ©L. По этой системе однозначно восстанавливается оптимальная система решений, которая дает совокупность подмоделей исходной модели

УГД (0.1) в случае общего уравнения состояния допускают 11-мерную алгебру Ли Ьц, операторы базиса которой в декартовой системе координат имеют вид:

Хх = дх, Х2 = ду, Хг = dz, Х4 = tdx + ди, Х5 = tdy + dv,

Х6 = tdz + dW: Х7 = ydz - zdy + vdw - wdv,

X8 = zdx - xdz + wdu - udw, X9 = хду - ydx + udv — vdu,

Хю = dt, Хц = tdt + xdx + уду + zdz. В цилиндрических координатах x — (x,r,9), u — (U,V,W), y = rcos9, z = r sin 9, и = U, v = V cos6 — W sin9, w — Vsm9 + Wcos9 эти операторы примут следующий вид [55]:

Е.

0.2)

Xi=dx, Х2 = cos 9дг - — (дв + Wdv ~ Vdw), г cos 9

Х3 = sin 9дг +-(дв + Wdv ~ Vdw), Х4 = tdx + ди, г

Хъ = cos 9{tdr -dv)- ^t {дв + Wdv - (v - jj dw^j , sin9{tdr + dv) + ^-t (de + Wdv - (V - Q dw) , Х7 = дв, X8 = sin в(гдх - xdr + Vdu ~ Udv)+ cos 9

Хд — — cos в(гдх — хдг -f- Vdu ~ Udv) + + sin в (wdu - Udw + Wdv - Vdw)) ,

Xiq = dt, Хц = tdt + xdx + rdr.

Дифференциальные уравнения, описывающие физические процессы, часто содержат параметры или функции, которые находятся экспериментально и поэтому пе строго фиксированы. В этом случае говорят, что система содержит произвольный элемент А в виде указанных параметров и функций. Группа преобразований, допускаемая системой независимо от имеющегося произвольного элемента, называется основной группой. При конкретизации элемента А допускаемая системой группа может только расширяться. При этом уравнение состояние так же конкретизируется (элемент А связан с уравнением состояния). Возникает задача групповой классификации: для данной системы дифференциальных уравнений Е найти основную группу и все ее спецификации элемента А (спецификации уравнения состояния), дающие расширение основной группы.

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

Как говорилось выше, массив решений при реализации программы ПОДМОДЕЛИ определяется оптимальной системой подалгебр. О его широте в случае УГД можно судить по следующим цифрам: для общего уравнений состояния газа оптимальная система подалгебр содержит 221 представитель [26], для УГД с уравнением состояния p = Bp< + g{S), (0.3) где В ф 0,77^0,1 — постоянные, Bj > 0, g(S) — произвольная функция энтропии, 606 представителей [50], для политропного газа 1342 представителей [4], для одноатомного газа 1817 представителей [61].

На основе программы ПОДМОДЕЛИ к настоящему времени получено большое количество конкретных примеров точных решений УГД. Качественное исследование подмоделей двумерных, винтовых, вращательных и других движений газа приведено в работах [15, 16, 17, 27, 47, 48, 51].

Примеры нетривиальных частично инвариантных решений можно найти в [13, 14, 29, 30, 38, 66]. Класс движений газа, в которых давление зависит только от времени (барохронные движения) изучен в работах [64, 65, 67]. Введено понятие «простого» (с групповой точки зрения, аналога постоянного) решения [35]. В этой же работе описаны все простые решения, не относящиеся к специальным типам движения газа. Построено решение, описывающее двумерное стационарное течение газа с замкнутыми линиями тока [36]. Это решение является аналогом двумерного течения идеальной жидкости из [63]. В работе [11] приведены все инвариантные подмодели с двумя независимыми переменными для УГД с общим уравнением состояния.

В процессе реализации программы ПОДМОДЕЛИ были получены ответы на вопросы, имеющие общетеоретическое значение. Создан эффективный алгоритм построения оптимальных систем подалгебр конечномерных групп Ли [25], дан критерий .х-автономии группы, допускаемой некоторой системой дифференциальных уравнений [24]. В работе [28] введен новый классификационный признак для частично инвариантных подмоделей. Замечено, что все инвариантные подмодели газовой динамики могут быть записаны в некоторой канонической форме [46], а в [32] доказано это. Дан алгоритм получения подмоделей в канонической и дивергентной формах, пригодный для реализации на компьюторе. Реализации этого алгоритма в виде действующей программы посвящена работа [62]. Установлено, что все инвариантные подмодели имеют эквивалентное представление в виде интегральных законов сохранения [31]. Получена теорема, задающая иерархию инвариантных подмоделей [33], которая нужна для упорядочивания и структуирования множества подмоделей, которое может быть достаточно большим и трудно обозримым.

Настоящая диссертация основана на материале, полученным автором при исследовании в рамках программы ПОДМОДЕЛИ для уравнений сжимаемой жидкости (0.1), (0.3). Этим объясняется тематическая направленность изложенного ниже материала.

Актуальность работы. Отыскание и классификация точных решений уравнений газовой динамики, которые описывают широкий круг явлений в механике сплошной среды, являются актуальными проблемами. Используя подгруппы преобразований, допускаемые УГД, можно находить подмодели (системы уравнений с меньшим количеством независимых переменных) движений газа и параметрические семейства точных решений. Полный список подмоделей ещё не завершен. Важной задачей является построение и исследование новых подмоделей с целью получения точных решений и их интерпретации с точки зрения физики.

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

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

• групповая классификация гидродинамической системы ранга 2 стационарного типа.

Перейдем к описанию содержания диссертации.

В первой главе рассматривается модель движения сжимаемой жидкости (0.1), (0.3). Показано, что уравнение состояния (0.3) для специальных функций g(S) согласуется с уравнением состояния, описывающим поведение реальных плотных сред (жидкостей и твердых тел), когда давление и температура могут достигать больших значений. Выписывается базис операторов 13-мерной алгебры Ли L13, допускаемой моделью движения ежимаемой жидкости. Он вычислен в [26] и состоит из базиса операторов (0.2) алгебры Ли Ьц, допускаемой УГД в случае общего уравнения состояния, дополненного операторами растяжения и переноса

Х12 = tdt - иди - vdv - wdw - (7 - 2)рдр - 7рдр, Хгз = др: (0.4) которые в цилиндрических координатах принимают вид

Х12 = tdt - иди - Уду - Wdw - (7 - 2)рдр - 7рдр, Х13 = др, где 7 = 27/(7 — 1): 7 2. Рассмотрены 12 двумерных подалгебр из оптимальной системы алгебры Ли L13, возникающих когда параметр 7^0,1 произвольный. Кратко приведена известная теория построения инвариантных и частично инвариантных подмоделей через инварианты подалгебр. Вычислены инварианты рассматриваемых двумерных подалгебр и сведены в таблицы. На трех подалгебрах получены инвариантные подмодели ранга 2, которые приведены к стационарному каноническому типу (время не есть инвариант подалгебры) по правилу, приведенному в работе [53]. Для оставшихся девяти подалгебр получаются регулярные частично инвариантные подмодели ранга 3 дефекта 1. Регулярные означает, что в независимых инвариантах нет искомых газодинамических функций. Показано, что частично инвариантные подмодели редуцируются к инвариантным подмоделям на одномерных подалгебрах. Из полученных девяти инвариантных подмоделей ранга 3 две подмодели приведены к эволюционному (время есть инвариант подалгебры) каноническому типу, семь подмоделей приведены к стационарному каноническому типу по правилу, описанному в работе [53].

Во второй главе рассмотрена 4-мерная подалгебра из оптимальной системы подалгебр алгебры Ли L13, допускаемой уравнениями движения сжимаемой жидкости. Операторы базиса подалгебры следующие: дх, Х4 = tdx + ди, Ххз = др,

0.5)

2 27

Х12 = tdt - иди - vdv - wdw---рдр---рдр.

7—1 7 — 1

Разыскиваются регулярные частично инвариантные решения ранга 2 дефекта 2. Представление регулярного частично инвариантного решения ранга 2 дефекта 2 таково u = u(t,x), V = t~lvi(y,z), w = t~lwi(y,z), p = f^FTp! (0, г), p = pit, ж).

Доказана

Теорема 0.1. Регулярные частично инвариантные решения ранга 2 дефекта 2 уравнений движения сжимаемой жидкости (0.1), (0.3) построенные на подалгебре с базисом (0.5) определяются только пятью подмоделями. Для представления решения (0.6) с уточнением

1 -27 u = t~ хщ(у, z) +u0(t,y, z), p = t-r~1pi(y,z) (0.7) получается редуцируемая подмодель с представлением решения ранга 2 дефекта 1 на подалгебре {Xi, Х4, Х12}, состоящая из инвариантной подсистемы

A«i = mi-«i, Dxvi + pi'lpiy = vi, Diwi + p^pu = wi,

Dm + Pi(viy + wu) = pi ~ '

Dipi + 7BpJ(vly + wiz) =--pi - jBplui,

7-1 и уравнения для функции щ: tuot + D1U0 = —itiuo, где D\ = г>1<Эу + w\dz.

Для представления решения (0.6) с уточнением

0.8) и = Гг(х + К'1 (щ(у, z) - elnt))

0.9)

V = «^(Pifo, + ДГж - e(lnt + 1)), 7 = -1, получается редуцируемая подмодель с представлением инвариантного ре

13 шения ранга 2 на подалгебре {eXi 4- КХ12, Х4 + КХ\3}:

Dim = е — К2р11, D1V1 + p^piy = vi,

D\W\ + рг Vz = wi, Dipi + P\(viy + wu) = -2/01, (0.10)

Dipi - - В Pi l(viy + wlz) = pi + Bpi1 - ub где e, К ф 0 — постоянные, D\ — v\dy + W\dz.

Для представления решения (0.6) с уточнениями и = хГ1 + К 2)po - r1 (A(t)epo + B{t)K2pu l)] , р = гФт, p = (pife, *)A> + S/oJ) +

7-1 , 27 7 (B(t) + l) K2Pq1 ~ + l) epo] , (0.11) и = if"1 (^г/1(з/, z)po - A(t)epo - B(t)K2pq , p = t"^,

P = Kx + Г^т (pi(3/, 2)/90 4- Bpl) +

0.12) p = t 7 , 2 U = --ret , p = t -У-1, 7-1 (0.13) где e, К ф 0, /?о ф 0 — постоянные, получаются соответственно пере- определенные подмодели: 7+1

D\Ui =--щ + е, Div I + Ply = Vi, DiWi + pu = wh

7-1

3-7 27 viy + w\z =--, Dipi =--pi - ui]

7 — 1 7 — 1

З7 — 1

D1U1 =-+ +piv = vi, Diwi+piz = wi,

7-1

0.15)

2 27 =-D1P1 =-~P\-Ui\

7 — 1 7 — 1

D1V1 + piy = vi, D1W1 + pu = wi,

0.16)

2-7 27

Viy + W lz =-D1P1 =--pi,

7 — 1 7 — 1 где D\ = v\dy + w\dz.

Замечание 0.1. В Теореме 0.1 при t < 0 в формулах для представлений р, р переменную t нужно брать по модулю.

Для подмодели (0.16) найдены решения в случае линейного поля скоростей, в случае ортогональной системы координат и в общем случае. Результаты сформулированы в виде теорем.

Теорема 0.2. Подмодель (0.16) в случае линейного поля скоростей

Щ = ацу 4- a\2z + а\, w\ = а2\у + а22z + а2, где a,i, aij — постоянные, допускает 3 типа решений, которые классифицируются по кратным волнам с точностью до преобразований подобия УГД [22, §23]:

1. 0-волна v\ — w\ — pi — 0 при 7 = 2;

2. 1-волна vi = ацу + auz, wi = a2iy + a22z: pi = 0, Oil + «22 = 1, «11^22 — «12«21 = 0 при 7 = 3/2;

3. 2-волна

Vi — У, Wi — z, pi — 0 при 7 -- 4/3; vi = —2y, wi - —2z, pi = -3(y2 + z2) при 7 = 2/3.

Теорема 0.3. Подмодель (0.16) в ортогональной системе координат с линиями тока в качестве координатных линий имеет решения только с линейным полем скоростей. Решения классифицируются по кратным волнам [22, § 23] с точностью до преобразований подобия УГД и образуют три типа решений:

1. 0-волна vi = w\ = р\ = 0 при 7 = 2;

2. 1-волна г>1 = (у sin во + z cos во) sin во, w\ = {у sin во + z cos во) cos во,

Pi — 0 при 7 = 3/2, где во — постоянная;

3. 2-волна vi — у, wi = z, pi = 0 при 7 - 4/3; vi = -2 у, wi = -2 z, pi = -3 {у2 + z2) при j = 2/3.

Теорема 0.4. Подмодель (0.16) имеет решения только для четырех значений показателя 7 уравнения состояния (0.3)

7 = 2: pi = v\ = wi = 0;

7 = 3/2: pi = 0, г/i = isinJ, w\ — A cos 5, где A = у sin 6 + z cos 5 — x(0> X — x(0 ~~ произвольная функция, т(1) = —у cos J + z sin 6: 5 — постоянная;

7 = 4/3: pi=0, Уг=у-Ф(1), wi = z-^f{l), где функции Ф, Ф связаны равенством z - ф)ф' = (у - ф)ф';

7 = 2/3: Pl = -3(i/2 + z2), = «л = -2г.

В третьей главе рассмотрена гидродинамическая система ранга 2 стационарного типа

Fu = иих + vuy + R~\guPx + gPy) = аи,

Fv = uvx + vvy + R-\gPx + gvPy) = a\

Fr = uRx 4- vRy + R(ux + vy) = aR = ifaf, (0.17)

Fp = uPx + vPy + A(R, P)(ux + vy) = ap,

Fk = uw* + vwy = afc, где к — 1,., n, фиксированные функции ^(ж, ?/), ди{х, ?/) > 0, ^(яг, у) > 0, Д2 = gugv д2 > Q. a?{x,y,u,v,R,PiWl,.,wn), v = u,v,R,P, l,.,n, — линейные или квадратичные функции по переменным и, и, wk. Функция Л считается произвольной (произвольным элементом) и определяется уравнением состояния.

Системы вида (0.17) появляются при рассмотрении регулярных частично инвариантных подмоделей ранга 2 УГД. Решается задача групповой классификации этой системы, включающая вычисление допускаемой ею алгебры и возможных специализаций элемента А. Операторы алгебры Ли, допускаемой системой (0.17), представляются в виде: х = сдх+еду++?dv+ckdwk+tRdR+fdP, где = У> u-> vj Р> Р> wl> ••> wn). В качестве независимых производных выбираются величины их, vyi wly, Ry, Рх, Ру, г = 1, .,п. Далее записываются условия инвариантности системы (0.17) относительно оператора X с неопределенными множителями:

X\{Fy - av) = (Л" + Ких + АХ + ЛuRRy + ЛuPlPx + \vP2Py + А>*)(</ - Fp), где Xi — продолженный на первые производные оператор X. Неопределенные множители A", AJ, j — и, v, R, PI, Р2,1,., п, зависят от переменных ж, и, и, u;\ R, Р. Зависимые производные выражены через независимые из уравнений системы (0.17).

Условия инвариантности содержат не более чем квадратичные слагаемые по независимым производным. Искомые функции ^ не зависят от производных. В результате расщепления уравнений по степеням независимых производных получается переопределенная система уравнений для нахождения координат допускаемого оператора. Приравнивание нулю коэффициентов при различных квадратичных и линейных слагаемых по независимым производным приводит соответственно к 2-уравнениям и 1-уравнениям [24]. Оставшиеся уравнения не содержат независимых производных и образуют 0-уравнения [24].

Интегрирование 2- и 1-уравнений приводит к двум Леммам. Лемма 0.1. Общее решение 2-уравнений имеет вид:

Z" = e(x,y), f = №,n) лу = 0. (0.18)

Замечание 0.2. Утверждение Леммы 0.1 согласуется с работой [49]. Лемма 0.2. Интегрирование 1-уравнений дает представление для координат оператора X и для неопределенных множителей Хи при А ^ 0: с = иси+иех+<, с = ™++ к = У, < Р)), = - 2W - 2п),

Ptffoi/J+^Or.y), Xu = -Ufi, Xv = А* = Яд, (0.19)

Afc — 2 Л2 ^(Яг1 - шр) - Хр = -ш-еР~ Ац, где 11 = 2 Л2 (rap — R~l)uj 4- w = uj(x1y,wt1JJS(R,P))>

J = gVy? 2guv + flfV+ 2 Л2 т(Л, P),

ASp + RSr = 0, Amp + йтд = AR~l. (0.20)

Кроме того, 1-уравнения включают равенства для определения координат £у: п = ^ + g{gvT4yx - \{gvYldgv = g + ^ - \{guYldgu,

1 - 2 (0.21)

А2 (G/T1^ + (sT1® = dg - -g [(guYldgu + (</T , г^е 5 = £хдх + iydy, и классификационное соотношение: 6f) + ArR(£i - 2w - 2n) = (0.22)

При A = 0 в равенствах (0.19), (0.20) нужно считать £p — £p(x, у, P), m — 0; Xя — —R^pp + R/i, a S(R,P) заменяется на P.

Замечание 0.3. В Лемме 0.2 выражения для Хр, ц уточняют формулы из работы [Щ.

О-уравнения имеют вид: и2ехх + 2 uv£y + v2eyy + аи{£ + uujjJu + 2и)+ av{ixy + uujjJv) + R~l{gu£ + + u{uux + vuy)+ +u{uJx + vJy)uj J + aluujwi — Xau — (ap — Aaf)\u = 0; еще одно уравнение получается из предыдущего уравнения заменой и <—> г/, х <—>■ у;

4- шак + + + +(и Jx 4 vJy + auJu + avJv)£kj

-Xak - (ap - Aaf )A* = 0; (0.23) u(g - 2n + g + фх + v(g -2n + g + фу + afo;

Xaf — wwx — vuy — (uJx 4- vJy 4- auJu + avJv)uj— —tfujyji - (ap - Aaf)R-1XR = 0; v& + + + + + фу + U0JX + VU>y+ (uJx + vJy + auJu -I- avJv)u)j 4- alujwi + af (w 4-~Xap - (ap - Aaf)Xp = 0.

Эти уравнения более точные, чем в [49].

С помощью классификационного соотношения (0.22) получено 13 случаев специализации элемента А (уравнения состояния). Десять случаев совпадает с результатами групповой классификации УГД в случае общего уравнения состояния [26, табл. 1]. В оставшихся трех случаях получаются более общие уравнения состояния, чем в [26, табл. 1]. Они появляются при групповой классификации подмоделей, построенных по подалгебрам 2.6 и 2.19 из [26, табл. 6]. Групповая классификация для первой подмодели приведена в [68], а для второй в [12].

С помощью полученных формул вычислены алгебры допускаемые подмоделями (0.8) и (0.10) из Теоремы 0.1.

В четвертой главе в силу Теорем 0.1 и 0.4 формулируется Утверждение 0.1. Уравнения движения сжимаемой жидкости (0.1), (0.3) на подалгебре с базисом (0.5) имеют решения вида (0.13) только для четырех значений показателя 7: о ож п 1 Вр0-х2 7 = 2: и = v = w = P=~t2> р =--' В(р0-1)-х2 V2B (0-24) г4 ' а~

X „ 2 — Y(IJ) 1

7 = 3/2: и = 3—, i; = 0, w =-p =

BpT - За;2 B(pS/2 - 1) - За2 /3~1

P =-i5-' 5 =-*-' a=sV2BW ж 2/-Ф(0 1

7 = 4/3: и = 4—, w =-p =

BpT ~ „ B{pl'Z - 1) - 6x2 /41 1 p =-P-' ^ =-P-' V3 Й' где функции Ф; Ф связаны равенством (z — Ф)Ф' = (у — Ф)Ф';

7 = 2/3: г/= р = р = (я,"1'3 - 3r2)

0.26)

0.27)

S = (ад1'3 - 1) - Зг2) Л г где г = \х\, U = -2-; и, v, w — компоненты вектора скорости, р — плотность, р — давление, S — функция энтропии, а — скорость звука.

В решениях (0.24)-(0.27) берется ро = 1. В решении (0.24) взято В = 2 и рассматриваются одномерные движения сжимаемой жидкости. В решениях (0.25)-(0.27) меняются масштабы координат, скорости, давления, функции энтропии и скорости звука х —>• л/уВж, й —► у/^Вй, р —» -S1 —> 7-05, a —> у/^Ва, масштабы времени и плотности не меняются t —» t, р р. В решении (0.26) масштабы функций Ф, Ф меняются как и масштаб координат. Изменение масштабов сделано для того, чтобы в решениях не было постоянной В. Уточнения для решений (0.24)-(0.27) примут соответственно вид:

7 = 2 : р = (2 - х2) Г\ S = -хН~\ а = 2I*!"1; (0.28)

7 = 3/2: р = 3(2/9 - ж2) Г6, 5 = -Зж2Г6, a = (0.29)

7 = 4/3 : р = 3/4 (1 - 8ж2) Г8, 5 = -6x2t~\ a = М"1; (0.30) 7 = 2/3: р = 3/2(1 -2r2)i4, 5 = -3r¥, a = И"1. (0.31)

Для каждого из решений (0.24)-(0.27) с уточнениями (0.28)-(0.31) движение жидкости рассматривается в области, в которой давление положительно. Предполагается, что на поверхности с нулевым давлением жидкость переходит в другую фазу, которая не может быть описана уравнением состояния (0.3). Такие поверхности далее называются поверхностями фазового перехода (фазовыми границами). Звуковой поверхностью называется поверхность, на которой скорость частиц жидкости равна скорости звука.

Решение (0.24) с уточнением (0.28) описывает коллапс частиц в точку или мгновенный точечный источник при одномерном изоэнтропическом движении. Решение обрезается характеристиками и склеивается в непрерывное решение одномерной подмодели в конечной области. Показана возможность непрерывного периодического неизоэнтропического движения сжимаемой жидкости в ограниченной области под действием поршня.

Решение (0.25) с уточнением (0.29) описывает расслоенное по параллельным плоскостям трансзвуковое движение сжимаемой жидкости из мгновенного дозвукового источника. Это движение можно интерпретировать как движение жидкости в канале, если ограничить движение стенками. Для частиц, вылетевших из данной точки, звуковая линия есть эллипс, лежащий в плоскости движения частиц и с центром в точке вылета частиц. Под звуковой линией понимается линия, на которой скорость частиц равна скорости звука.

Хорошо известно решение задачи о расширяющемся поршне для одномерных движений со сферическими волнами (см., например, [37, §21]). Это решение автомодельное (инвариантно по группе растяжений) и существует, когда поршень расширяется в покоящийся политропный газ из точки (начала координат) с постоянной скоростью. По покоящемуся газу идет сферическая ударная волна. Есть инвариантные решения, для которых мировые линии движения частиц являются ветвями парабол г = где а > 0 — постоянный показатель автомодельности, Rq характеризует скорость движения частиц [1, §20].

Решение (0.27) с уточнением (0.31) неавтомодельное, но частично инвариантное решение уравнений движения сжимаемой жидкости (0.1) с уравнением состояния (0.3) вида р = -р2/3 + S, где S — функция энтропии. Мировые линии для этого решения задаются равенством г = Rot , где

Ro — постоянная, характеризующая скорость частиц, а показатель отрицательный. Решение описывает расширение при t < 0 и сжатие при t > О сжимаемой жидкости под действием сферического поршня. Жидкость движется с трансзвуковыми скоростями. В области движения жидкости не возникает ударных волн. Описан характер распространения сферического возмущения, возникшего на поршне, в одномерной постановке задачи. При сжатии жидкости возмущение со временем сходится к центру. При расширении жидкости часть возмущения задерживается в сферическом слое данного радиуса, другая часть возмущения доходит до поверхности фазового перехода.

Решение (0.26) с уточнением (0.30) описывает движение сжимаемой жидкости из мгновенного дозвукового источника через звуковую поверхность к фазовой границе. Возможно два случая многообразия источника: кривая в плоскости х = 0 или точка в начале координат. В первом случае, когда многообразие источника есть окружность единичного радиуса, звуковая поверхность представляет собой эллипсоид. Во втором случае звуковая поверхность так же является эллипсоидом, построены поверхности уровня t = const характеристического коноида. Выяснено, что если возмущение возникло в точке, расположенной внутри звуковой поверхности, то поверхности уровня коноида вложены друг в друга. Если возмущение возникло в точке на звуковой поверхности, то поверхности уровня коноида касаются друг друга в одной точке. Если возмущение возникло в точке за звуковой поверхностью, то поверхности уровня коноида имеют общие сечения.

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

На защиту выносятся следующие положения:

- групповое свойство гидродинамической системы стационарного типа с двумя независимыми переменными, которая обобщает частично инвариантные подмодели ранга 2 уравнений газовой динамики;

- новые инвариантные подмодели сжимаемой жидкости ранга 2 и 3 на 12-и двумерных подалгебрах и подмодели ранга 2 на одной 4-мерной подалгебре из оптимальной системы 13-мерной алгебры Ли, допускаемой уравнениями движения сжимаемой жидкости;

- способ получения всех решений переопределенной подмодели сжимаемой жидкости ранга 2 дефекта 2, полученной на 4-мерной подалгебре;

- физическая интерпретация решений подмодели ранга 2 дефекта 2.

Научная новизна работы состоит в следующем:

- построены новые подмодели сжимаемой жидкости;

- найдены новые частично инвариантные решения уравнений сжимаемой жидкости и дана их физическая интерпретация;

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

Теоретическая значимость.

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

Разработан новый метод интегрирования переопределенных систем, основанный на использовании новых переменных лагранжевого типа .

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

Результаты работы использовались при составлении отчета о НИР Института механики УНЦ РАН № 01.200.211712 инв. № 02.2.006 07719 за 2008 г.

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

- 34-ой Региональной молодежной конференции «Проблемы теоретической и прикладной математики» (Екатеринбург, 2003);

- II конкурсе научных работ молодых ученых и аспирантов Уфимского научного центра и АН РБ (Уфа, 2003);

- 35-ой Региональной молодежной конференции «Проблемы теоретической и прикладной математики» (Екатеринбург, 2004);

- Всероссийской молодежной научной школе-конференции «Чеботарев-ские чтения по проблемам современного группового анализа и его приложениям в нелинейной механике» (Казань, 2004);

- XX Всероссийской школе-семинаре «Аналитические методы и оптимизация процессов в механике жидкости и газа» (Абрау-Дюрсо, 2004);

- II Всероссийской конференции «Актуальные проблемы прикладной математики и механики» (Абрау-Дюрсо, 2004);

- Всероссийской научной конференции «Современные проблемы физики и математики» (Стерлитамак, 2004);

- III конкурсе научных работ молодых ученых и аспирантов Уфимского научного центра и АН РБ (Уфа, 2004);

- 36-ой Региональной молодежной конференции «Проблемы теоретической и прикладной математики» (Екатеринбург, 2005);

- Российской конференции «Механика и химическая физика сплошных сред» (Вирск, 2007)

- Международной научной конференции «Дифференциальные уравнения и смежные проблемы», посвященной юбилею академиков Ильина В.А. и Моисеева Е.Н. (Стерлитамак, 2008);

- Научных семинарах Института механики Уфимского научного центра РАН под руководством д.ф.-м.н. С. Ф. Урманчеева, профессора В. Ш. Ша-гапова.

В 2005 г. работа связанная с результатамрг второй главы диссертации была удостоена диплома за II место на III конкурсе научных работ молодых ученых и аспирантов Уфимского научного центра и АН РБ.

Публикации. Основные результаты диссертации опубликованы в 10 работах [69]-[77], из них 4 входят в список ВАК.

Заключение диссертации по теме "Механика жидкости, газа и плазмы"

Заключение

На основании полученных в работе результатов можно сделать следующие выводы:

2. Найдены все решения одной частично инвариантной переопределенной подмодели сжимаемой жидкости методом обобщенных лагранже-вых переменных.

Список источников диссертации и автореферата по механике, кандидата физико-математических наук, Гарифуллин, Артур Рафаилевич, Уфа

1. Баутин С.П., Дерябин С.Л. Математическое моделирование истечения идеального газа в вакуум. Новосибирск: Наука. 2005.

2. Бучнев А.А. Группа Ли, допускаемая уравнениями движения идеальной несжимаемой жидкости // Динамика сплошной среды. ИГ СО АН СССР, 1971. Вып. 7. С. 212-214.

3. Бытыев В.О. Групповые свойства уравнений Навье-Стокса // Численные методы механики сплошной среды. 1972. Т. 3, № 5. С. 13-17.

4. Головин С.В. Оптимальная система подалгебр для алгебры Ли операторов, допускаемых уравнениями газовой динамики в случае полит-ропного газа // Новосибирск, 1996. Препринт/Ин-т гидродинамики. Сиб. отделение РАН. № 5-96.

5. Головин С.В. Точные решения для эволюционных подмоделей газовой динамики // ПМТФ. 2002, т. 43, № 4. С. 3-14.

6. Ибрагимов Н.Х. Классификация инвариантных решений уравнений двумерного нестационарного движения газа // ПМТФ. 1966. № 4. С. 19-22.

7. Ибрагимов Н.Х. Групповые свойства некоторых дифференциальных уравнений. Новосибирск, 1967. 60 с.

8. Ибрагимов Н.Х. Группы преобразований в математической физике. М.: Наука, 1983. 280с.

9. Капитанский JI.B. Групповые анализ уравнений Навье-Стокса и Эйлера при наличии вращательной симметрии и новые точные решения этих уравнений // ДАН СССР. 1978. Т. 243, № 4. С. 901-904.

10. Лапко Б.В. Построение оптимальных систем подгрупп группы Ли преобразований, допускаемой уравнениями газовой динамики // ДСС. 1973. Вып. 14. С. 112-119.

11. Мамонтов Е.В. Инвариантные подмодели ранга два уравнений газовой динамики // ПМТФ. 1999. Т. 40. № 2. С. 50-55.

12. Мамонтов Е.В. Групповые свойства 2-подмоделей класса S уравнений газовой динамики // Вестник НГУ. 2007. Том 7, вып. 1. С. 72-85.

13. Мелешко С.В. О неизэнтропических стационарных пространственных и плоских нестационарных двойных волнах // ПММ. 1989. Т. 53. Вып. 2. С. 255-260.

14. Мелешко С.В. Об одном классе частично инвариантных решений плоских течений газа // Дифференц. уравнения. 1994. Т. 30. № 10. С. 1825-1827/

15. Мелешко С.В. Групповая классификация уравнений двумерных движений газа // ПММ. 1994. Т 57, Вып. 4. С. 56-62.

16. Мелешко С.В. Групповая классификация уравнений движений газа в постоянном поле сил // ПМТФ. 1996. Т 37, № 1. С. 42-47.

17. Мустаев А.Ф. Подмодели винтовых галилеево-инвариантных течений в газовой динамике. Уфа. 1999. (Препр. / БашГПУ)

18. Никольский А.А. Инвариантное преобразование уравнений движения идеального одноатомного газа и новые классы их точных решений // ПММ. 1963. Т 27. № 3. С. 496-508.

19. Никольский А.А. Инвариантные преобразования уравнений движения идеального газа для специальных случаев // Инженерный журнал. 1963. Т 3. № 1. С. 140-142.

20. Овсянников JI.B. Групповые свойства уравнения нелинейной теплопроводности // Доклады АН СССР. 1959. Т. 125, № 3. С. 83-90.

21. Овсянников JI.B. Групповые свойства дифференциальных уравнений. Новосибирск, 1962. 240 с.

22. Овсянников JI.B. Групповой анализ дифференциальных уравнений. М.: Наука, 1978. 399 с.

23. Овсянников Л. В. Программа ПОДМОДЕЛИ. Новосибирск: Ин-т гидродинамики СО РАН, 1992. 11 с.

24. Овсянников Л. В. О свойстве х-автономии // Доклады РАН. 1993. Т. 330, № 5. С. 559-561.

25. Овсянников Л.В. Об оптимальных системах подалгебр // Доклады РАН. 1993. Т. 333, № 6. С. 702-704.

26. Овсянников Л. В. Программа ПОДМОДЕЛИ. Газовая динамика // Прикл. математика и механика. 1994. Т. 58, вып. 4. С. 30-55.

27. Овсянников Л.В. Изобарические движения газа // Дифференциальные уравнения. 1994. Т. 30, № 10. С. 1792-1798.

28. Овсянников Л.В. Регулярные и нерегулярные частично инвариантные решения // Доклады РАН. 1995. Т. 343, № 2. С. 156-159.

29. Овсянников Л.В. Особый вихрь // ПМТФ. 1995. Т. 36, № 3. С. 45-52.136

30. Овсянников JI.В. Регулярные типа (2,1) подмодели уравнений газовойдинамики // ПМТФ. 1996. Т. 37, № 2. С. 3-13.

31. Овсянников Л.В. Инвариантные интегральные законы сохранения //

32. Докл. РАН. 1996. Т. 351. № 5. С. 599-602.

33. Овсянников Л.В. Каноническая форма инвариантных подмоделей га-зовой динамики. Новосибирск, 1997 (Препр. / СО РАН, Ин-т гидродинамики; № 3 97). 41 с.

34. Овсянников Л.В. Об иерархии инвариантных подмоделей дифференциальных уравнений // Докл. РАН. 1998. Т. 361. Ш 6. С. 740-742.

35. Овсянников Л.В. Некоторые итоги выполнения программы «ПОД-

36. МОДЕЛИ» для уравнений газовой динамики // Прикл. математикаи механика. 1999.Т. 63, вып. 3. С. 355-372.

37. Овсянников Л.В. О «простых» решениях уравнений динамики политропного газа // ПМТФ. 1999. Т. 40. № 2.Прикл. математика имеханика. 1999.Т. 63, вып. 3. С. 355-372.

38. Овсянников Л.В. Плоские течения газа с замкнутыми линиями тока // Докл. РАН. 1998. Т. 361. № 1. С. 51-53.

39. Овсянников Л.В. Лекции по основам газовой динамики. Москва,

40. Ижевск: Ин-т компьютеных исследований, 2003.

41. Овсянников Л.В., Чупахин А.П. Регулярные частично инвариантныеподмодели уравнений газовой динамики // Прикладная математикаи механика. 1996. Т. 60, вып. 6. С. 990-999.

42. Павловский Ю.Н. Исследование некоторых инвариантных решенийуравнений пограничного слоя // Журнал вычисл. матем. и матем.физики. 1961. № 2. С. 280-294.

43. Пухначев В. В. Групповые свойства уравнений Навье-Стокса в плоском случае // ПМТФ. 1960. № 1. С. 83-90.

44. Пухначев В.В. Инвариантные решения уравнений Навье-Стокса, описывающие движения со свободной границей // ДАН СССР. 1972. Т.202, № 2. С. 302-305.

45. Сидоров А.Ф. Избранные труды. Математика. Механика. М.: Физ-матлит., 2001.

46. Станюкович К.П. Неустановившиеся движения сплошной среды. М.:1. Гостехиздат, 1955.

47. Тешуков В.М. Пространственный аналог центрированных волн Ри-мана и Праидтля-Майера // Журнал прикладной математики и технической физики. 1982. № 4. С. 98-106.

48. Хабиров С. В. Нестационарное инвариантное решение уравнений газовой динамики, описывающее растекание газа до вакуума // Прикл. математика и механика. 1988, том 52, вып. 6. С. 967-975.

49. Хабиров С.В. К анализу инвариантных подмоделей ранга три уравнений газовой динамики // ДАН. 1995, т. 341, № 6. С. 764-766.

50. Хабиров С.В. Подмодель винтовых движений в газовой динамике // Прикл. математика и механика. 1996, т. 60, вып. 1. С. 53-65.

51. Хабиров С.В. Винтовые движения в газовой динамике с давлением и плотностью, зависящими только от времени // Матем. заметки. 1996. Т. 59. Вып. 1. С. 133-141.

52. Хабиров С.В. Об общих свойствах инвариантных подмоделей ранга два в газовой динамике // Проблемы механики и управления. Институт механики УНЦ РАН. Уфа. Гилем, 1996. С. 102-111.

53. Хабиров С.В. Оптимальные системы подалгебр, допускаемых уравнениями газовой динамики. Уфа: 1998. (Препр. / УНЦ РАН. Ин-т механики).

54. Хабиров С.В. Подмодель вращательных движений газа в однородном поле сил // Прикл. математика и механика. 1998, т. 62, вып. 2. С. 263271.

55. Хабиров С.В. Подмодель вращательных движений в газовой динамике // ПМТФ. 1998, т. 39, № 6. С. 37-45.

56. Хабиров С.В. Приведение инвариантной подмодели газовой динамики к каноническому виду // Мат. заметки. 1999. Т. 66, вып. 3. С. 439444.

57. Хабиров С. В. Течения газа со спиральными поверхностями уровня // ПМТФ. 1999. Т. 40, № 2. С. 34-39.

58. Хабиров С.В. Инвариантные решения ранга 1 в газовой динамике // Тр. Междунар. науч. конф. « Моделирование, вычисления, проектирование в условиях неопределенности — 2000». Уфа: изд. УГАТУ, 2000. С. 104-115.138

59. Хабиров С.В. Нерегулярные частично инвариантные решения (2.1)эволюционного типа // Тр. Междунар. науч. конф. « Комплексныйанализ, дифференциальные уравнения и смежные вопросы». Т. 2.

60. Уфа: Ин-т математики с ВЦ УНЦ РАН, 2000. С. 151-155.

61. Хабиров С.В. Непрерывное радиальное ограниченное движение газапод действием поршня // ПМТФ. 2004, т. 45, № 2. С. 124-135.

62. Хабиров С.В. Инвариантные решения уравнений газовой динами-ки // Вестник УГАТУ. 2001, № 1(3) С. 47-52.

63. Хабиров С.В. Задача Гурса о непрерывном сопряжении радиальныхпрямолинейных движений газа // Мат. заметки. 2006, вып. 4. С. 601606.

64. Хабиров С.В. Частично инвариантные решения для подмодели радиальных движений газа // ПМТФ. 2007. Т. 48, № 5. С. 26-34.

65. Черевко А.А. Оптимальная система подалгебр для алгебры Ли операторов, допускаемых системой уравнений газовой динамики с уравнением состояния р = /(5)р5у/3. Новосибирск. 1996 (Препр. / СО РАН,

66. ИН-т гидродинамики; №4-96).

67. Черевко А.А. Построение канонических систем дифференциальныхуравнений для инвариантных подмоделей газовой динамики // Выч.технологии. 1998. Т. 3. № 6. С. 92-96.

68. Черный Г.Г. Плоские установившиеся автомодельные вихревые течения идеальной жидкости // Изв. РАН. МЖГ. 1997. № 4. С. 39-53.

69. Чупахин А.П. О барохронных движениях газа. 1997. // ДАН. 1997,т. 352, Ш 5. С. 624-626.

70. Чупахин А.П. Барохронные движения газа: общие свойства и подмодели типов (1,2) и (1,1). Новосибирск. 1998. (Препринт / Институтгидродинамики СО РАН; № 4-98). 79с.

71. Чупахин А.П. Небарохронные подмодели типов (1,2) и (1,1) уравнений газовой динамики. Новосибирск. 1999. (Препринт / Институтгидродинамики СО РАН; № 1-99). 40с.

72. Чупахин А.П. О регулярных подмоделях типа (1,2) и (1,1) уравненийгазовой динамики // ПМТФ. 1999. Т. 40. № 2. С. 40-49.

73. Khabirov S. V. Submodel of the Spiral Stationary Motion in Gas

74. Dynamics. Modern Group Analysis VII. Thondheim: Mars Publishers.1999. Norway. P. 181-187.

75. Гарифуллин А.Р. Подмодели сжимаемой жидкости на двумерных подалгебрах // Сибирский журнал индустриальной математики. 2003. Т. 6, № 1(13). С. 16-26.

76. Гарифуллин А.Р. Подмодели сжимаемой жидкости на двумерных подалгебрах // Труды 34-й Региональной молодежной конференции «Проблемы теоретической и прикладной математики». Екатеринбург, 2003. С. 101-105.

77. Гарифуллин А.Р. Подмодели ранга 2 стационарного типа уравнений газовой динамики для одной 4-мерной подалебры // Труды 35-й Региональной молодежной конференции «Проблемы теоретической и прикладной математики». Екатеринбург, 2004. С. 119-123.

78. Гарифуллин А.Р. Решения с линейным полем скоростей регулярной переопределенной подмодели сжимаемой жидкости ранга 2 дефекта 1 // Труды Всероссийской научной конференции «Современные проблемы физики и математики». Т. 1. Уфа: Гилем, 2004. С. 32-35.

79. Гарифуллин А.Р. Групповая классификация гидродинамической системы ранга два стационарного типа // Сибирский журнал индустриальной математики. 2004. Т. 7, № 3(19). С. 66-75.

80. Гарифуллин А.Р. Общее решение переопределенной подмодели сжимаемой жидкости ранга 2 дефекта 1 // Труды 36-й Региональной молодежной конференции «Проблемы теоретической и прикладной математики». Екатеринбург, 2005. С. 122-125.

81. Гарифуллин А.Р. Пример сферически симметричного движения сжимаемой жидкости // Сибирский журнал индустриальной математики. 2007. Т. 10, № 2(30). С. 45-52.

82. Гарифуллин А.Р., Хабиров С.В. Непрерывное сопряжение специальных неизэнтропических одномерных движений газа // Труды института математики и механики УрО РАН. Т. 14, № 1. 2008. С. 22-30.