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

• Гарифуллин, Артур Рафаилевич
• кандидат физико-математических наук
• 2009

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

Введение

Глава 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 уравнений газовой динамики и появляются при рассмотрении регулярных частично инвариантных подмоделей.

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

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

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

• групповая классификация гидродинамической системы ранга 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 ~ '

27

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 входят в список ВАК.

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

Заключение диссертации по теме «Механика жидкости, газа и плазмы», Гарифуллин, Артур Рафаилевич

Заключение

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

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

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

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

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

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

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

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.