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

• Степанова, Ирина Владимировна
• кандидат физико-математических наук
• 2008

Оглавление диссертации кандидат физико-математических наук Степанова, Ирина Владимировна

Введение

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

1.1 Описание системы уравнений

1.2 Определяющие уравнения.

1.3 Групповая классификация системы (1.4).

1.4 Групповая классификация уравнений стационарного турбулентного пограничного слоя

1.5 Построение фактор-систем и инвариантных решений

1.6 Об одном стационарном автомодельном решении.

Введение диссертации (часть автореферата) на тему «Групповая классификация и точные решения уравнений двух моделей гидродинамики»

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

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

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

В настоящее время изучение симметрий дифференциальных уравнений проводится в рамках современного группового анализа, объединяющего в себе три научных направления — метод первого интеграла, классический групповой анализ С. Ли и дискретно-групповой анализ. Обыкновенные дифференциальные уравнения были первым объектом приложения групп Ли к дифференциальным уравнениям. Позднее Биркгоф привлек внимание к приложениям групп Ли к дифференциальным уравнениям механики жидкости [8]. Однако в данном случае использование групповых свойств дифференциальных уравнений для поиска частных решений было использовано в единичных случаях и носило, в основном, иллюстративный характер.

Систематические исследования по применению методов группового анализа для широкого круга физически важных задач были начаты Л.В. Овсянниковым и его учениками в конце 50-х годов прошлого столетия [33,34]. Применительно к уравнениям механики жидкостей и газа, эти работы продолжаются Л.В. Овсянниковым, его научной школой по настоящее время. В аналитической механике теоретико-групповые методы были использованы еще в работах А. Пуанкаре и Н.Г. Четаева. Одни из первых работ по исследованию групповых свойств уравнений механики жидкости и газа в нашей стране были выполнены Ю.Н. Павловским [37], В.В. Пухначевым [39], С.В. Хабировым [57]. В работах этих авторов показано, каким образом групповые свойства можно использовать для решения физически важных задач. В настоящее время наряду с указанными авторами исследование уравнений механики сплошной среды продолжается В.К. Андреевым, О.В. Бытевым, О.В. Капцовым, С.В. Мелешко, С.И. Сенашовым, А.П. Чупахиным, С.В. Головиным и др. [1,2,10,31,46,58,65]. Из работ зарубежных авторов следует отметить монографию П. Олвера [36].

В 1991 г. на Седьмом Всесоюзном съезде по теоретической и прикладной механике JI.B. Овсянниковым была предложена концепция программы ПОДМОДЕЛИ, направленная на полное и систематическое изучение групповых свойств различных моделей механики сплошной среды [35]. Под руководством JI.B. Овсянникова группой исследователей ведется активная работа по реализации этой программы для уравнений газовой динамики. Работа над данным проектом привела к оформлению более четких алгоритмов, используемых в групповом анализе, а также к расширению его теоретической базы. Что касается исследования других моделей, то здесь следует отметить работы [65,66], которые посвящены систематическому изучению уравнений Навье-Стокса с помощью теоретико-групповых методов.

Определим основные понятия, использующиеся далее.

Рассматривается пространство Z — W1 х Rm = X х U. Переменные разделяются на два типа: х = (а;1,., хп) — независимые переменные, и = (it1,., ит) — функции от х.

Система дифференциальных уравнений Е допускает группу G преобразований всех участвующих в Е величин (независимых и зависимых переменных), если система Е остается неизменной при всех преобразованиях, принадлежащих группе G. Здесь группа G есть локальная од-нопараметрическая группа Ли с каноническим параметром а е А с 1, которая задана преобразованиями Та : х — /(ж,а), х £ Мп. Фиксируя точку х и изменяя параметр а, получаем кривую в пространстве Шп, касательный вектор £ к которой имеет компоненты df(x,a) ем = а=О да

Инфинитезимальным оператором группы G называется дифференциальный оператор

Пространство „ дх1' г=1

Z = X х U х U х . х U, и = ( n °Sua . I к Iks I. OXJl.X^ J называется к-м продолжением пространства Z. Независимыми переменными в продолженном пространстве Z являются переменные х, функк ции и и все производные и^ ^ до к-го порядка включительно.

Пусть в пространстве Z действует локальная группа Ли G с инфинитезимальным оператором X = u)dxi -\-rja(x, и)диа. Действие группы

G естественным образом распространяется на пространство Z. При этом к преобразования производных вычисляются по обычным правилам математического анализа. Продолженной группе G отвечает продолженный к оператор

X = ?дх< + 7fdu* + Cfdu? + . + Cl.iAt.ik> а = 1,га; i, ii, .Ли — 1,., n, по повторяющемуся индексу производится суммирование. Координаты продолженного оператора X вычисляютк ся по формуле - + «Ч-агде Di — оператор полного дифференцирования по переменной %\

Di = dxi + ufdua + . 4- ufj^jdu?^, + ••••

Важно помнить, что во всех этих формулах координаты х,и и их производные и% • должны рассматриваться как независимые переменные.

J1 • • 4 J s

В случае, когда система дифференциальных уравнений Е содержит произвольный элемент А в виде неопределенных параметров и функций, возможно подчиненных некоторым условиям, естественным образом возникает задача групповой классификации. Такую систему обозначим через Е(А), а через GE{A) — группу, допускаемую системой Е(А). Тогда ядром основных групп называется группа GEq, равная пересечению всех групп GE(A). Задача групповой классификации заключается в следующем: для системы Е(А) найти ядро основных групп GEq и все специализации произвольного элемента А, дающие расширения группы

GE0.

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

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

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

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

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

Очень многие технически важные течения являются турбулентными. В таких течениях на главное движение налагается пульсационное движение, в результате чего возникает перемешивание отдельных частей жидкости. Это пульсационное движение в своих деталях настолько сложно, что возможность его теоретического расчета является затруднительной. Между тем перемешивание жидкости, вызываемое пульсационным движением, придает турбулентному течению особенности, резко отличающие его от ламинарного течения. Наблюдения позволяют обнаружить, что при турбулентном течении скорость и давление в фиксированной точке пространства не остаются постоянными во времени, а очень часто и очень неравномерно изменяются. Такие изменения скорости и давления, называемые пульсациями, являются наиболее характерным признаком турбулентности. Так как вследствие сложности пульсационного движения чисто теоретический расчет турбулентного течения до настоящего времени невозможен, то закономерности развитого течения приходится искать лишь для осредненных по времени величин, характеризующих это движение. Появление в уравнениях осредненного турбулентного течения членов, содержащих произведение пульсациоииых компонент скорости, входящих в функцию турбулентного трения, делает их незамкнутыми. На этой почве возникли полуэмпирические теории турбулентности, содержащие как правило две эмпирически определяемые "константы турбулентности". Плодотворное развитие этого направления связано с работами Л. Прандтля (1925) [38], К. Тейлора (1915, 1932), А.Н. Колмогорова (1941). Фундаментальное значение в дальнейшем развитии теории турбулентности имели работы JI. Лойцянского [28], М.Д. Миллионщико-ва [32], Г. Шлихтинга [63]. Классические опыты И. Никурадзе [26] дали обоснование полуэмпирическим теориям и способствовали их развитию. Но тем не менее до последнего времени теория турбулентности не могла обойтись при количественном анализе без опытных коэффициентов, далеко не всегда ясной природы. В настоящей работе сделана попытка исследования основных закономерностей осредненного пристенного турбулентного течения жидкости без привлечения эмпирических коэффициентов. А именно, предлагается использовать метод группового анализа для исследования уравнений модели плоского турбулентного пограничного слоя.

В последнее время расширился круг задач, связанных с так называемыми естественно- или свободноконвективиыми течениями. Естественная конвекция — один из видов макроскопического движения, который интенсивно изучается в современной фундаментальной механике жидких сред. Это связано не только с развитием математического аппарата и средств вычислительной техники, но и, возможно, с тем, что этот характерный вид течений, возникающий под действием сравнительно малых сил плавучести, вызванных малыми градиентами плотности жидкости, находящейся в поле массовых сил, часто встречается в природных условиях и технологических процессах, имеющих широкое применение и познавательную ценность. Именно термоконцентрационной конвекцией вызывается образование пространственно-регулярной по глубине тонкой структуры океана, атмосферы планет и звезд, мантии Земли и других геофизических систем [24]. Конвективные течения часто упорядочены и в горизонтальном направлении. Развитие экспериментальных и теоретических исследований естественной конвекции привело к ее выделению в самостоятельный раздел механики жидкости и газа. Изучению моделей конвекции посвящено множество работ. В качестве примера следует привести исследование модели микроконвекции [41], конвекции с учетом эффектов термодиффузии и диффузионной теплопроводности [12], а также стоит отметить монографии [1, 2], в которых наряду с классическими моделями исследуются уравнения термокапиллярного движения, уравнения конвекции с коэффициентами переноса, зависящими от температуры. Что касается исследования явления термодиффузии или, так называемого эффекта Соре, то его изучению посвящены работы [17,42,60-62,64]. Результаты исследований в этой области применяются в теплоэнергетике, металлургии, метеорологии, химии, кристаллофизике и. т. д. Действительно, такие усложненные модели с большей точностью (по сравнению с классическими) описывают реальные физические процессы и в последнее время активно используются в вычислительной гидродинамике. В связи с этим является актуальной задача качественного исследования таких моделей с помощью методов группового анализа.

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

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

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

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

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

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

Первая глава посвящена исследованию групповых свойств уравнений плоского турбулентного пограничного слоя. Исследованы случаи нестационарного и установившегося течения. Появление в модели турбулентного течения в пограничном слое членов, отвечающих за дополнительное касательное напряжение, делает их незамкнутыми, что, естественно, затрудняет исследование. В работе предлагается считать, что дополнительное касательное напряжение д есть функция продольной пространственной координаты у и производной от продольной компоненты скорости по у. иу. Задача групповой классификации уравнений турбулентного пограничного слоя решалась относительно этой функции: д(у,иу). Была найдена основная алгебра Ли, которая оказалась бесконечномерной в случае нестационарного течения и конечномерной в случае установившегося течения. Выделены ядра основной алгебры Ли в случаях, когда классифицируемая функция д зависит только от иу или д = — иу + а(у), где а(у) ф const — произвольная нелинейная гладкая функция. Результат групповой классификации представлен в виде таблиц, где указаны специализации классифицируемой функции и операторы, расширяющие основную алгебру Ли для найденных специализаций. Нужно особо отметить, что при составлении этих таблиц существенно использовались преобразования эквивалентности, также найденные в работе. Заметим также, что когда классифицируемая функция д = О, исследуемая система описывает движение в ламинарном пограничном слое. Групповой анализ такой системы уже проводился Л.В. Овсянниковым в зависимости от того, является функция давления произвольной или заданной. Результаты, полученные JI.B. Овсянниковым, согласуются с результатами, полученными в данной главе, как в случае нестационарного, так и в случае установившегося течения [34].

Далее построено несколько примеров инвариантных систем на вычисленных операторах для соответствующих специализаций классифицируемой функции. К сожалению, далеко не все полученные факторсисте-мы удается проинтегрировать аналитически. Для случая, когда исходная система описывает стационарное течение и функция турбулентного напряжения зависит лишь от одной переменной, а именно д = д(иу), вычислена оптимальная система подалгебр первого порядка. На операторе растяжения, вошедшем в оптимальную систему подалгебр, построена факторсистема, которая, в отличие от исходной системы, уже является системой обыкновенных дифференциальных уравнений, описывает автомодельное решение и становится проще для интегрирования. Для численного решения этой задачи с использованием метода Рунге-Кутта удалось поставить следующие граничные условия: условия равенства нулю скоростей на стенке (условия прилипания) и равенства нулю турбулентного напряжения вне пограничного слоя. Поскольку граничные условия сформулированы в таком виде, то помимо метода Рунге-Кутта использовался метод пристрелки для поиска начального условия на стенке для турбулентного трения. Следует особо отметить, что в силу произвольности функции д(иу) она была выбрана согласно эмпирической формуле Прандтля: д = ozoUy, где ао = pi2 (р — плотность среды, I — длина пути смешения) — некоторая положительная постоянная, введенная Прандт-лем. Для сравнения найденных характеристик турбулентного течения с ламинарным, на том же операторе растяжения построено автомодельное решение для ламинарного пограничного слоя и сделаны некоторые выводы относительно поведения скоростей, касательного напряжения и толщины пограничного слоя в случае турбулентного и ламинарного течения.

Во второй главе рассматривается другая не менее значимая модель современной гидродинамики: модель конвекции бинарной смеси с учетом эффекта термодиффузии и сил плавучести. После введения новых переменных для температуры (Т = (x~D)T/(aD), где х ф D ф 0 — коэффициент температуропроводности, D ф 0 — коэффициент диффузии, а ф 0 — коэффициент Соре) и концентрации (с = Т+с) система упростилась: уравнения для температуры и концентрации стали иметь одинаковую дифференциальную структуру. Для данной системы была решена задача групповой классификации относительно функции Я, определяющей силу плавучести, с учетом, что она зависит от трех параметров: давления, температуры и концентрации. После построения классифицирующего уравнения стало ясно, что нужно исследовать отдельно три случая: 1) функция R не зависит от давления, 2) функция R зависит от давления линейно, 3) R зависит от давления нелинейно. Результат исследования представлен в виде таблиц, где выписаны специализации классифицируемой функции и операторы, допускаемые исследуемой системой в зависимости от найденных специализаций. Следует отметить, что при составлении этих таблиц существенно использовались преобразования эквивалентности, найденные в работе.

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

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

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

1. Результаты групповой классификации уравнений плоского нестационарного турбулентного пограничного слоя.

2. Результаты групповой классификации уравнений плоского стационарного турбулентного пограничного слоя.

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

4. Результаты групповой классификации трехмерных уравнений конвекции с учетом эффекта термодиффузии и сил плавучести.

5. Точные решения инвариантных подмоделей уравнений термодиффузии и их физическая интерпретация.

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

XLV Международной конференции "Студент и научно-технический прогресс"(Новосибирск, 2007),

VII Всероссийской конференции по математическому моделированию и информационным технологиям (Красноярск, 2006),

VIII Всероссийской конференции по математическому моделированию и информационным технологиям (Новосибирск, 2007),

Конференции молодых ученых Института вычислительного моделирования СО РАН (Красноярск, 2007, 2008 гг.),

Конференции молодых ученых Красноярского научного центра (Красноярск, 2007, 2008 гг.),

Международной научной конференции "Современные проблемы математического моделирования и вычислительных технологий - 2008" (Красноярск, 2008 г.),

Семинаре Института вычислительного моделирования СО РАН "Математическое моделирование в механике" под руководством профессора В.К. Андреева.

Семинаре Института гидродинамики им. М.А. Лаврентьева СО РАН под руководством академика Л.В. Овсянникова.

Публикации. Основные результаты диссертационной работы опубликованы в работах [43,48-55].

Структура и объем работы. Диссертация состоит из введения, двух глав, заключения и списка литературы, который содержит 67 наименований. Общий объем диссертации 98 страниц, включая 9 рисунков и 8 таблиц, которые вынесены в конец соответствующих глав.

Заключение диссертации по теме «Дифференциальные уравнения», Степанова, Ирина Владимировна

Заключение

Сформулируем основные результаты диссертационной работы:

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

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

3. Решена задача групповой классификации трехмерных уравнений термодиффузии относительно функции R, определяющей силу плавучести и произвольно зависящей от трех параметров: температуры, давления и концентрации. Указана замена переменных для температуры и концентрации, с помощью которой уравнение диффузии можно привести к однородному, в результате чего исходная система существенно упрощается. При решении классифицирующих уравнений выделяются три случая зависимости функции R от давления: a) R зависит от давления линейно, b) R не зависит от давления, c) R зависит от давления нелинейно.

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

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

Список литературы диссертационного исследования кандидат физико-математических наук Степанова, Ирина Владимировна, 2008 год

1. Андреев В.К., Капцов О.В., Пухначев В.В., Родионов А.А. Применение теоретико-групповых методов в гидродинамике. Новосибирск: Наука, 1994. 319 с.

2. Андреев В.К., Бублик В.В., Бытев В.О. Симметрии неклассических моделей гидродинамики. Новосибирск: Наука, 2003. 352 с.

3. Барренблатт Г.И. Турбулентные пограничные слои при очень больших числах Рейнольдса// Успехи математических наук, 2004, Т. 59, № 1. С. 45-62.

4. Бекежанова В.Б. Об устойчивости плоского слоя в модели микроконвекции: Дис. . канд. физ.-мат. наук. Красноясрк, 2003.8. биркгоф Г. Гидродинамика. М.: ИЛ, 1963. 244 с.

5. Бочаров О.Б., Васильев О.Ф., Овчинникова Т.Э. Приближенное уравнение состояния пресной воды вблизи температуры максимальной плотности// Известия АН, 1999, Т. 35, №4, С. 556-558.

6. БЫТЕВ В.О. Групповые свойства уравнений Навье-Стокса// Численные методы механики сплошной среды. 1972. Т. 3, N2 5. С. 13-17.и. бэтчелор Дж. Введение в динамику жидкости. м.: Мир, 1973. 760 с.

7. Гершуни Г.З., Жуховицкий Е. М., Непомнящий А. А. Устойчивость конвективных течений. М.: Наука, 1972. 392 с.

8. Гершуни Г.З., Жуховицкий е. М. Конвективная устойчивость несжимаемой жидкости. М.: Наука, 1989. 320 с.

9. ГИЛЛ А. Динамика атмосферы и океана. Т.2, М.: Мир, 1986. с.

10. ГРЮ К.Э., ИВБС T.JL Термическая диффузия в газах. М.: Государственное издательство технико-теоретической литературы, 1956. 184 с.

11. Дейч М.Е., Зарянкин А.Е. Гидрогазодинамика. М.:Энергоиздат, 1984. 384 с.

12. ЗАЙЦЕВ В.Ф., Полянин А.Д. Справочник по обыкновенным дифференциальным уравнениям. М.: Физматлит, 2001. 576 с.

13. ИБРАГИМОВ Н.Х. Опыт группового анализа обыкновенных дифференциальных уравнений. М.-.Знание, 1991. 47 с.21. ибрагимов н.х. Группы преобразований в математической физике. М.: Наука, 1983. 280 с.

14. ИБРАГИМОВ Н.Х. Азбука группового анализа. М.:3нание, 1989. 45 с.

15. Компаниец ji.a., Якубайлик Т.В., Гаврилова Л.В., Гуре-вич К.Ю. Модели экмановского типа в задачах гидродинамики. Новосибирск:Наука, 2007. 156 с.26. кутателадзе С.с. Пристенная турбулентность. Новоси-бирск:Наука, 1973. 228 с.

16. Ландау Л.Д., Лившиц Е.М. Гидродинамика. М.: Наука, 1986. 734 с.28. лойцянский Л.г. Механика жидкости и газа. М.: Дрофа, 2003. 840 с.

17. ЛоЙЦЯНСКИЙ Л.Г. Ламинарный пограничный слой. М.: Физматгиз, 1962. 480 с.

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

19. ОлВЕР П. Приложения групп Ли к дифференциальным уравнениям. М.: Мир, 1989. 639 с.

20. ПУХНАЧЕВ В.В. Модель конвективного движения при пониженной гравитации// Моделирование в механике. 1992.Т. 6(23),№ 4. С. 4756.42. рабинович Г.Д. Разделение изотопов и других смесей термодиффузией. М.:Атомиздат, 1981. 142 с.

21. Родионов А.А., степанова И.В. Групповая классификация уравнений модели конвекции с учетом сил плавучести// Вычислительные технологии, Новосибирск, 2008, Т. 13, № 5,С. 61-69.

22. РЫЖКОВ И.И. Инвариантные подмодели и точные решения уравнений термодиффузии: Дис. . канд. физ.-мат. наук. Красноярск, 2005.

23. Рыжков И.И., Андреев В.К. Групповая классификация и точные решения уравнений термодиффузии // Дифференциальные уравнения. 2005. Т. 41, № 4. С. 508-517.

24. Сенашов С.И., Киряков П.П., Яхно А.Н. Приложение симмет-рий и законов сохранения к решению дифференциальных уравнений. Новосибирск, 2001. 192 с.

25. СЛЕЗКИН Н.А. Динамика вязкой несжимаемой жидкости. М.: Государственное издательство технико-теоретической литературы, 1995. 520 с.

26. СТЕПАНОВА И.В. Групповая классификация уравнений стационарного плоского турбулентного пограничного слоя// Вестник КрасГУ, 2006, № 1. С. 114-119.

27. СТЕПАНОВА И.В. Инвариантная подмодель уравнений плоского стационарного турбулентного слоя// Материалы конф. молодых ученых ИВМ СО РАН, Красноярск: ИВМ СО РАН, 2007, С. 30-36.

28. СТЕПАНОВА И.В. Групповая классификация уравнений плоского нестационарного пограничного слоя// Вычислительные технологии, Новосибирск, 2007, Т. 12,№ 6. С. 101-108.

29. СТЕПАНОВА И.В. Групповые свойства уравнений турбулентного пограничного слоя// Тезисы докладов VII Всероссийской конференции по математическому моделированию и информационным технологиям. Красноярск: ИВМ СО РАН, 2006. С. 70.

30. СТЕПАНОВА И.В. Групповая классификация уравнений одной модели конвекции// Тезисы докладов VIII Всероссийской конференции по математическому моделированию и информационным технологиям. Новосибирск: ИВТ СО РАН, 2007. С. 75.

31. СТЕПАНОВА И.В. Об инвариантных моделях термодиффузионных течений// Сборник трудов конференции молодых ученых Красноярского научного центра СО РАН. Красноярск: КНЦ СО РАН, 2008. С. 44-45.

32. СТЕПАНОВА И.В. О некоторых точных решениях модели конвекции с учетом сил плавучести// Материалы конф. молодых ученых ИВМ СО РАН, Красноярск: ИВМ СО РАН, 2008, С. 10-17.

33. СТЕПАНОВА И.В. Об инвариантных свойствах уравнений модели конвекции// Тезисы докладов XXXIII Дальновосточной математической школы-семинара им. академика Е.В.Золотова, Владивосток: Изд-во Дальнаука, 2008, С. 99-100.

34. УИТТЕКЕР Э.Т., ВАТСОН Дж.Н. Курс современного анализа. Т.2. М.:Физматлит, 1963. 516 с.

35. ХАБИРОВ С.В. Классификация нелинейных волновых уравнений// Труды института метематики и механики УрО РАН, 2007, Т. 13, № 2. С.1-15.

36. ЧЕРЕВКО А.А, Чупахин А.П. Стационарный вихрь Овсянникова. Препринт №1, Новосибирск: ИГиЛ, 2005, 52 с.

37. ШАШИН В.М. Гидромеханика. М.:"Высшая школа", 1990. 384 с.

38. Шерстянкин П.П., Куимова Л.Н., Толоконцева О.Ю., Хохлов В.В. Потенциальная температура в и плотность р* глубоких персных вод на примере озера Байкал// ДАН, 2001, Т. 378, № 1. С. 98-102.

39. Шерстянкин П.П., Куимова Л.Н., Потемкин В.Л. Основные закономерности термохалинного режима глубинной зоны озера Байкал на основе Т, S-анализа// ДАН, 1997, Т. 355, № 5. С. 638-687.

40. ШЕРСТЯНКИН П.П. Равновесные процессы в природных водах// ДАН, 1998, Т. 360, № 6. С. 819-822.63. шлихтинг Г. Теория пограничного слоя. М.: Издательство иностранной литературы, 1956. 528с.

41. Gebhart В., Mollendorf J.С. Buoyancy-induced flows in water under cjnditions in which density extrema may arise // J. Fluid Mechanics, 1978, V. 89, part 4, P. 13-17.

42. Fushchych W., Popowych R. Symmetry reduction and exact solutions of the Navier-Stokes equations. I // Nonlinear mathematical physics. 1994. V. 1, № 1. P. 75-113.

43. Fushchych W., Popowych R. Symmetry reduction and exact solutions of the Navier-Stokes equations. II // Nonlinear mathematical physics. 1994. V. 1, № 2. P. 158-188.

44. TRITTON, D.J. Physical Fluid Dynamics. Oxford University Press, 1988. 519 p.