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

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

НИКИФОРОВА Светлана Витальевна

ПРИЛОЖЕНИЯ ГРУПП ЛИ К КОНСТРУИРОВАНИЮ ДИВЕРГЕНТНЫХ ФОРМ УРАВНЕНИЯ ЭЙЛЕРА И МОДЕЛИРОВАНИЮ В ЗАДАЧАХ ЛАМИНАРНОГО ПОГРАНИЧНОГО СЛОЯ

Специальность 05 13 18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ

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

Казань - 2007

003059089

Диссертация выполнена в Казанском государственном техническом университете им А Н Туполева

Научный руководитель

доктор физико-математических наук, профессор Гараев Кавас Гараевич

Официальные оппоненты

доктор физико-математических наук, профессор Чугунов Владимир Аркадьевич

доктор физико-математических наук, профессор Кирпичников Александр Петрович

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

Научно-исследовательский институт математики и механики им Н Г Чеботарева (г Казань)

Защита состоится «29» мая 2007 г в 14 час 00 мин на заседании диссертационного совета Д 212 079 01 в Казанском государственном техническом университете имени А Н Туполева по адресу 420111, г Казань, ул Карла Маркса, 10, в зале заседаний Ученого совета КГТУ

С диссертацией можно ознакомиться в научной библиотеке КГТУ им А Н Туполева

Автореферат разослан « Л » апреля 2007 г

Ученый секретарь диссертационного совета доктор физико-математических наук, ,

профессор п г данилаев

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

Актуальность темы. Групповой анализ дифференциальных уравнений возник как научное направление в работах выдающегося норвежского математика XIX века Софуса Ли (1842-1899 г г) и служил главной составной частью его важнейшего творения - теории непрерывных групп Первоначальная основная задача группового анализа - вопрос о разрешимости в квадратурах обыкновенных дифференциальных уравнений - была решена самим Ли, но не нашла практического применения

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

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

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

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

Одним из основных инструментов построения законов сохранения физических процессов, допускающих вариационную формулировку, является первая теорема Эмми Нетер, которая устанавливает связь между инвариантностью вариационного интеграла относительно конечномерной группы Ли и дивергентными формами уравнения Эйлера Эта теорема дает достаточное условие существования законов сохранения для уравнений Эйлера Н X Ибрагимов дал новое доказательство этой теоремы для п-мерных интегралов на языке инфинитезимальных операторов Ли, что в сочетании с техникой группового анализа, развитого Л В Овсянниковым, дает возможность удобного способа построения законов сохранения Используя понятие слабого лагранжиана, он также установил не только достаточные, но и необходимые условия существования законов сохранения уравнений Эйлера

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

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

Цели работы.

Получение неклассических первых интегралов в простейшей задаче вариационного исчисления

Конкретизация идей С - инвариантности функционалов применительно к вариационным задачам с 2 и 3 независимыми переменными ^ Получение новых автомодельных решений уравнений ламинарного пограничного слоя при сверхзвуковых режимах обтекания в плоском и осесимметричном случаях

Методы исследования. Для решения поставленных задач в работе использовались инфинитезимальный аппарат Ли-Овсянникова, теория инвариантных вариационных задач Э Нетер, теория ламинарного пограничного слоя, методы вычислительной математики

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

Научная новизна. Впервые найдены неклассические первые интегралы уравнения Эйлера в простейшей задаче вариационного исчисления

Дано приложение однопараметрической группы Ли к вариационной задаче с подвижным концом

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

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

соответствующие фактор-системы, которые представляют собой системы

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

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

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

Апробация. Основные результаты диссертационной работы докладывались и обсуждались на-

■S Всероссийской молодежной школе-конференции по математическому

моделированию, алгебре и геометрии (г Казань, 1998г), S Международной летней школе-семинаре по современным проблемам

теоретической и математической физики (г Казань, 1999г), S VIII Четаевской Международной конференции по аналитической механике,

устойчивости и управлению движением (г Казань, 2002г ), S Всероссийской молодежной школе-конференции по математическому

моделированию, алгебре и геометрии (г Казань, 2003г ), ■S XVIII Международной научной конференции «Математические методы в

технике и технологиях ММТТ-18» (г Казань, 2005г ), ■S II Всероссийской научной конференции «ММ-2005» (г Самара, 2005 г ),

XXV Российской школе по проблемам науки и технологий, посвященной 60-летию Победы (г Миасс, 2005г ),

Международной молодежной научной конференции, посвященной 1000-летию города Казани «Туполевские чтения» (г Казань, 2005 г )

Публикации. Основные результаты диссертационной работы опубликованы в журналах (3 статьи), в трудах и материалах Международных и Всероссийских конференций (8 тезисов) В работе [2], основываясь на теории Ь* - инвариантности функционалов, разработанной К Г Гараевым, автором построена дивергентная форма уравнения Эйлера-Остроградского, и задача вычисления экстремального значения функционала по области V сведена к вычислению интеграла по поверхности 5, ограничивающей эту область Автор является дипломантом конкурса научных работ по присуждению именных стипендий Главы Администрации г Казани по итогам II семестра 2001-2002 г г

Положения, выносимые на защиту.

Неклассические первые интегралы уравнения Эйлера в простейшей задаче вариационного исчисления

Конкретизация идей Ь* - инвариантности функционалов применительно к вариационным задачам с 2 и 3 независимыми переменными

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

Структура и объем диссертационной работы. Диссертация состоит из введения, трех глав, заключения, библиографического списка использованной литературы и приложений Работа изложена на 84 страницах основного текста, иллюстративный материал представлен в виде 5 рисунков, приложения содержат 33 таблицы, библиография включает 53 наименования

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В главе 1 рассмотрены приложения групп Ли к вариационным интегралам

В разделе 1.1. получены новые первые интегралы для уравнения Эйлера в простейшей задаче вариационного исчисления

ь

Рассматривается функционал •/[>>(*)] = ^{х,у,у')сЬс, где у(х) и 1.(х,у,у')-

а

достаточно гладкие функции своих переменных

Случай I Лагранжиан задается в виде Ь- /(у')Р(х,у) Требуется построить первые интегралы Необходимое и достаточное условие инвариантности функционала относительно непрерывной группы преобразований с

инфинитезимальным оператором X-<^х{х,у)-~—ь^у(х,у)~—

запишется в виде

(ах+Ь)" ^Чау+Л)3

= 0 При этом

координаты оператора имеют вид =ах+Ь,£у =ау+с1

Первый интеграл для уравнения Эйлера запишется в виде

[ау+с1- /(ах+Ь)]К(х, у) ЩУ'^-+/<У Жх, у)(ах+Ь) = сошг ¿У

Случай 2 Лагранжиан задается в виде 1 = а(х,у)у'2 +Ь{х,у)у'+с(х,у) Первый интеграл для уравнения Эйлера в этом случае запишется в виде

а(х, у)Щу (х, у) у' - р{х)у 2 ]+ Ь(х, у)су (х, >■) + с(х, у)р(х) = сот Г,

где а(х,у),р(х) = £х- произвольные дифференцируемые функции, функции Ь(х,у) ,с(х,у) и координата ¿¡у(х, у) удовлетворяют определенным соотношениям

В разделе 1.2. рассматривается функционал

-Ф(-*.>')]= Я У, "(*> у\их (х, у), и у (х, у) )с1хс1у (1)

Для лагранжиана

Ь = А(х, у,и)и\ +В(х,у,и)ихиу +С(х,у, + /•"(*, у, г*), (2)

где А(х,у,и),В(х,у,и),С(х,у,и),Р(х,у,и)- непрерывно дифференцируемые в области 5 функции переменных х,у,и, построена дивергентная форма уравнения Эйлера - Остроградского вида

1 он.

(?„ - "А - "А )— + Ц,

(3)

где Вх,Оу - операторы полного дифференцирования соответственно по

независимым переменным х и у, 4х(х,у,и\^у(х,у,и\^и(х,у,и) - координаты

инфинитезимального оператора X - и)^+(х, >>, +ди (х, у, и)~

Теорема 1 Если коэффициенты лагранжиана А(х,у,и),В(х,у,и),С(х,у,и), Р(х,у,и) и координаты инфинитезимального оператора удовлетворяют определенным соотношениям, то существует дивергентная форма уравнения Эйлера - Остроградского вида (3)

Там же формулируются прямая и обратная задачи Прямая задача задана структура лагранжиана, требуется найти координаты инфинитезимального оператора, допускаемого данным функционалом Обратная задача задана группа непрерывных преобразований, требуется найти условия на коэффициенты лагранжиана, для которых существуют дивергентные формы для уравненияй Эйлера-Остроградского (обратная задача была использована для интегрирования уравнений в частных производных еще самим Софусом Ли, а для вариационной задачи - сформулирована К Г Гараевым)

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

В разделе 1.3. дается конкретизация идей классической и конформной инвариантности в задаче экстремума для двойного интеграла Рассматривается функционал (1) Экстремальное значение функционала с лагранжианом (2) определяется по формуле

-М -и £ -и £ \2Ли+Ви)+£\ Ли2 + Ви и +

IV н х х у у Х У / •'ХI х х у у у

¡ь-

где координаты инфинитезимального оператора £х(х,у,и\^у{х,у,и\^и(х,у,и)г коэффициенты лагранжиана А(х, у, и), В(х, у, и), С(х, у, и), Р(х,у, и) и параметр а удовлетворяют определенным соотношениям

Полученные результаты обобщены на трехмерный случай Рассматривается функционал

Аи{х,у,2)]= \\\Ь\с,у,2,и(х,у,г\их (:х,у,г\и (х,у,г\и2 {х,у,г))гЫуск (4)

V

Теорема 2 Если функционал (4) конформно-инвариантен относительно непрерывной группы преобразований с оператором

^и экстРемаль и(х>У'2) найдена, то экстремальное

значение функционала можно вычислить по формуле

/V **

(fw ~«x4x+ |cos(n,x)+

8их

5L

Su ~"x4x ~uyty-"ztz -+Цу

huy

-Ux£x-"yZy-Uz4z)^-+L£z cos(n.z)

cos(«, y)+

dS

где соэ^.л^соз^,у), 005(77,2) - направляющие косинусов внешней нормали к поверхности

Если лагранжиан в функционале (4) не зависит явно от х,у,г или и, то

•"iSix

s

s s s

L-u,

"x 8L ' Su,

Su y J I 8uy 8uz

( 8L

eos(n,y)-Uy -cos(/J,x)+-eos(n,z)

V3"* ^z

L-u,

5 L 8u,

\

cos(n, z)-M.

8L , ч Si , ч

cos( n, x)+-COS(í?,J')

8u

du

8L , . 8L , . 8L , ч cos(n, x)+-eos (л, y)+-- COS(íi, z)

У 1иК

dS, dS, dS,

8u

Su,-

8u-

Таким образом, задача вычисления экстремального значения функционала по области V сведена к вычислению интеграла по поверхности 5", ограничивающей эту область

В разделе 1.4. дано приложение теории групп Ли к вариационной задаче с подвижным концом

Рассматривается функционал ^у&^^'^Ы.х^^^сЬс, где функции

х0

у(х) и ¿(д:,у,у') - достаточно гладкие функции своих переменных

Уравнение Эйлера для заданного функционала, как известно, имеет вид

¿у-Ох(Ьу') = 0 Если концы экстремали закреплены, то имеем дело с краевой задачей для у=у{х) Если же правый конец не закреплен, то для отыскания х\

и У = >'(*1) следует использовать условия трансверсальности (А -у'1.у' )&с + 1.у'5у = 0 при х = х, В частности, если у = у{х\) задано, а переменная X] неизвестна, то для ее отыскания имеем уравнение

1-у'1у

х-х

1

= 0

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

•а

переносов по координате хс оператором Х=~ и уравнение Эйлера допускает

первый интеграл И, следовательно, граничная задача сводится к задаче

¿-/¿/=с, у(хо) = уо

Показано, что с-0 в двух случаях

1) Правый (левый) конец экстремали у = у(х) закреплен, если координаты

инфинитезимального оператора = 0, ¿;у = 0 при х = х, (при х = хо )

2) Координаты X] и у = у{х1) свободны

В главе 2 с использованием инфинитезимального аппарата Ли-Овсянникова конструируются новые автомодельные решения уравнений ламинарного пограничного слоя в задаче обтекания непроницаемого цилиндрического тела сверхзвуковым потоком газа

В разделе 2.1. рассматриваются уравнения ламинарного пограничного слоя на непроницаемом цилиндрическом теле, обтекаемом потоком газа под нулевым углом атаки

¿<о»>+|С">-0

р = рЯТ,м = МеОтЬ(т)

с граничными условиями

у = 0,*>0 и = 0,и = 0,Я = #№(х), 3'->» и->ие(х),Н ->#е, х = 0,.у>0 и = ие(0),Н = Не

Здесь приняты обычные в теории пограничного слоя обозначения В переменных Дородницына эта краевая задача записывается в виде

?21

и5 + XV/ = О

иц/ц+у>К~

Ь(г)к,-^±\к2+2а2икя} \--\у}\^Ъ{т)Кг +2Ъ{т)иК,-2и^^КЯ+2а}иЯ2 ]

(5)

^Рг

Здесь Я, = д2и/д12 ,К, = 32у//с'2

/ = 0 « = 0,^=0,^ = 1-^ 5=0 и= 1,У=0

(6)

В разделе 2.2. исследуются групповые свойства системы (5) с граничными условиями (6) Инфинитезимальный оператор, допускаемый этой

системой, имеет вид Х^^+С^+Си-^+С^+С^+Ск^+Ск

д

дК

& ь" ди ' ' ' ** 8Я

Определяющие уравнения для координат инфинитезимального оператора

записываются в виде

5«

дз

Ь(т)

КИ+2а}иРГ

ди 8э ды Зг ди дз

+ +

^ ' •'[¿(¡r) fc(r) Pr dr v ' V lfj

+ (' _ + ^~-^ - "21} + + 2«c?"ÍH

dr2

x(\'2 + + 2cc^uRZk +

+\\—Lya<aX1^Wr)R2 -u^^ipt + laluR^^s +wR~-»2j+

+ 2«2„«2 jJ+J, j2¿(r)J?f ^ + и2^ + 2«2<K)«2 +

+ +2aea>2¿rs +2a}ucjKR + 2a}uR2y

la» 5/ J s/f Uv 5í J

Частное решение этой системы отыскивалось в виде A\s + Л2,

CíM= íwM=Ciw+«Bi(s),^ =0, c*(r)=АЛ, як,

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

переменной s

Подробно рассмотрены два случая Случай 1. Число Прандтля Рг = 1, функция безразмерной температуры 6(г) = 1 (что соответствует линейной зависимости вязкости от температуры), температура стенки постоянна Tw=const

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

у +J_ +1 / i ---ll ае' 2с* / ' с начальным условием ос,(0) = 0 (что

соответствует обтеканию затупленного тела сверхзвуковым потоком газа)

Здесь С[, С* - произвольные постоянные, штрих означает дифференцирование по переменной х = х/1

Важно отметить, что в окрестности плоской критической точки это

уравнение имеетрешение ае = Схт, где C = Q(1-C*)/C*, m = C*/(l-C*), что совпадает с известным классическим решением

Группа, допускаемая системой (5), порождается операторами

v. _,3 , i 3 w д R д КЗ v__р-г-л d . д

Инварианты однопараметрической группы, порождаемой оператором Xi,

имеют вид /](^) = «,/2(^) = ^,/з = w i11 2, где £ = f s-11 2 Соответствующая фактор-система запишется в виде d/, d2/,

2S 1 d? У 2 "

= (7) i=o/,=o,ij=o,;2=i-r, (8)

2

ь / _7

£ «> /[ -» 1, /г -> 0

Здесь С* ,rw- произвольные постоянные

Случай 2 Pr*l, b(r) — COtlSt, Т\у — const

Показано, что ае(х) = const (что соответствует обтеканию клина сверхзвуковым потоком газа) Группа, допускаемая системой (5), порождается операторами

v „ д ,t д w 8 R 8 К 8 _ д _ d^/.s 8 .„ы,^ д

Соответствующая фактор-система в этом случае запишется в виде

Л,

1-/2-Л2

J* dl\ _ dl-j _ 2g d4 d£

0,

de

£ = 0 /1=0,/з = 0,/2 = 1-ги,

f —>co /j —>l,/2 —»0

= 0,

ЛГ

= 0

(10)

Здесь C*,b,C = 2aJ,rw- произвольные постоянные

В разделе 2.3. для решения краевых задач (7), (8) и (9), (10) использовался метод пристрелки

Случай 1. Решение краевой задачи (7), (8) не зависит непосредственно от вида функции ае(х) Последняя определяется значением постоянной С* В результате решения двухточечной нелинейной краевой задачи получены

I <

значения /j(0), ^(0) при различных значениях параметров С*и rw, которые

приведены в Приложении

На рисунке 1 представлены распределения безразмерной скорости ае(х) на внешней границе пограничного слоя и соответствующие им формы

обтекаемого профиля f(x*) при различных значениях параметра С* е (0,0 5],

где ось х* совпадает с осью симметрии плоского профиля, /¿¿с)2

При С* й (0,0 5] задача Коши допускает лишь тривиальное решение, не имеющее физического смысла

Зависимость формы тела от давления определялась по модифицированной формуле Ньютона Форма профиля определялась согласно

дифференциальному уравнению I ^-рю/ро /^-рсо/ро'}.

где отношение давлении до и после скачка уплотнения в точке торможения

РО

2 5 /

потока определяется по формуле 1/мЩ /\66iMl

О 0 1 02 03 04 05 Об 07 08 09

01 02 03 04 05 06 07 08 09

Рис 1 Зависимость ае(х) и зависимость формы профиля от /(х*) 1- С* = 0,1, 2- С* = 0,2, 3 - С* = 0,3, 4- С* = 0,4, 5- С* = 0,5

На рисунке 2 представлены характерные графики функций

пРи Рг=07, = 0 5, С* = 05 Видно, что инварианты носят монотонный характер и удовлетворяют граничным условиям

Получены формулы для аэродинамических характеристик (касательного напряжения трения и локального теплового потока) на обтекаемом профиле Характерные графики при Рг = 0 7, тк =0 5 представлены на рисунке 3

Рис 3 Распределения касательного напряжения трения и локального теплового потока 1- С* = 0,1, 2- С* = 0,25, 3 - С* = 0,3, 4- С* = 0,5

Го

3 / 2

1

С 2у

/(.г-О

/¡(0),

Ох

90 =

Рг^о

г/0--О

г-О

/2(0)

Л

Случай 2. Для краевой задачи (9), (10) расчеты для ае проводились по формуле Аккерета ае = |1 -/?/-у/мД-1 />/1+5/Л/Д , где р - угол полураствора

тонкого клина, Мсо- число Маха Результаты вычислительного эксперимента при Рг = 0 7, г№б(0,1), М«, е [2,5] приведены в Приложении

В главе 3 рассмотрена задача обтекания тела вращения потоком газа под нулевым углом атаки

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

В разделе 3.2. получены формулы для касательного напряжения трения и локального теплового потока

_ МеОЬ(Т

'.1

г/(г-0

Г-0 2 , г <3х

/¡(0),

<70 :

с рМеоЬ(.Т„)Те0\Ут

2у1(У 1) 0^11-а* г

Рг\ие0

: £ 2 V/(^ ~ 0 2 .

/НО)

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

1 Получены новые первые интегралы уравнения Эйлера в простейшей задаче вариационного исчисления

2 Основываясь на теории L*- инвариантности функционалов построена дивергентная форма уравнения Эйлера-Остроградского для вариационной задачи с двумя независимыми переменными Задача вычисления экстремального значения функционала по области V сведена к вычислению интеграла по поверхности S, ограничивающей эту область

3 Дано приложение теории групп Ли к вариационной задаче с подвижным концом

4 Исследованы групповые свойства уравнений ламинарного пограничного слоя в задаче обтекания непроницаемого цилиндрического тела сверхзвуковым потоком газа Получены определяющие уравнения для отыскания координат инфинитезимального оператора Ли В общем случае эти уравнения разрешить не удалось Однако подробно рассмотрены два случая

Случай 1 Рг = 1, Ь(т) = 1, tw = const, безразмерная скорость на внешней границе пограничного слоя подчиняется специальному дифференциальному уравнению, которое имеет решение, отличное от классического

Построена фактор-система, поставлена и решена соответствующая краевая задача Получено новое автомодельное решение Это стало возможным потому, что в качестве математической модели были взяты уравнения пограничного слоя в переменных Дородницына, которые редуцируют исходные уравнения в уравнения, содержащие меньшее число искомых функций, нежели исходные Использование аппарата Ли-Овсянникова позволило свести уравнения в частных производных к обыкновенным дифференциальным уравнениям, для которых формулируется соответствующая краевая задача Построены распределения безразмерной скорости aeix) на внешней границе пограничного слоя и соответствующие им формы обтекаемого профиля при различных значениях параметра С* е (0,0 5] Получены формулы для

аэродинамических характеристик обтекаемого тела (касательного напряжения трения и локального теплового потока) Построены характерные графики

Случай 2 Pr*l, b(r) = const, ru = const

Условие инвариантности уравнений пограничного слоя относительно оператора приводит к постоянной скорости на внешней границе ае = const, что соответствует случаю обтекания клина в сверхзвуковом потоке

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

Получены формулы для касательного напряжения трения и локального теплового потока

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

1 Никифорова С В О существовании неклассических первых интегралов в простейшей задаче вариационного исчисления // Казань Материалы Всероссийской молодежной школы-конференции по математическому моделированию, алгебре и геометрии - 1998 - С 25

2 Гараев К Г , Никифорова С В Использование классической и конформной инвариантности в проблеме экстремума для кратного интеграла - Казань Новейшие проблемы теории поля // Труды XI Международной летней школы-семинара по современным проблемам теоретической и математической физики Под редакцией проф Аминовой А В - 1999 -С 23-27

3 Никифорова С В Теоретико-групповой подход к исследованию одной задачи вариационного исчисления // Казань Вестник КГТУ - 1999 - №4 -С 53-57

4 Никифорова С В О достаточных условиях приводимости краевой задачи для уравнения Эйлера к задаче Коши // Казань Труды VIII Четаевской Международной конференции по аналитической механике, устойчивости и управлении движением - 2002 - С 353

Никифорова С В О достаточных условиях сведения краевой задачи для уравнения Эйлера к задаче Коши // Казань Вестник КГТУ - 2003 - №2 -С 41-42

Никифорова С В О новых фактор-системах уравнений пограничного слоя при сверхзвуковых режимах течения И Казань Материалы Всероссийской молодежной школы-конференции по математическому моделированию, алгебре и геометрии - 2003 - С 29

Никифорова С В О новых фактор-системах уравнений пограничного слоя при сверхзвуковых режимах течения И Казань Вестник КГТУ - 2004 - №3 -С 65-67

Никифорова С В Об одном численном методе решения уравнений пограничного слоя при сверхзвуковых режимах течения // Казань Сборник трудов XVIII Международной научной конференции «Математические методы в технике и технологиях ММТТ-18» — Т 1 -2005 —С 27-28 Никифорова С В Определение формы профиля в задаче обтекания сверхзвуковым потоком газа // Труды II Всероссийской научной конференции «ММ-2005» -Самара -2005 - Часть 2 - С 187-189

0 Никифорова С В Инвариантные решения уравнений пограничного слоя при сверхзвуковых режимах // Екатеринбург. Краткие сообщения XXV Российской школы по проблемам науки и технологий, посвященная 60-летию Победы -2005 -С 45-47

1 Никифорова С В О формах профилей в сверхзвуковом потоке, соответствующих новым автомодельным решениям уравнений пограничного слоя - Туполевские чтения // Материалы Международной молодежной научной конференции, посвященной 1000-летию города Казани - 2005 - Т 2 -С 78-79

Автор бшгодарит доцента В А Овчинникова за помощь оказанную при организации вычис чений

Формат 60x84 1/16 Бумага офсетная Печать офсетная Печл 1,25 Услпечл 1,16 Услкр-отт 1,21 Уч-издл 1,04 Тираж 100 Заказ К82

Типография Издательства Казанского государственного технического университета 420111 Казань, К Маркса, 10

Введение.

Групповой анализ дифференциальных уравнений возник как научное направление в работах выдающегося норвежского математика XIX века Софуса Ли (1842-1899 г.г.) [51,52] и служил главной составной частью его важнейшего творения - теории непрерывных групп. Первоначальная основная задача группового анализа - вопрос о разрешимости в квадратурах обыкновенных дифференциальных уравнений - была решена самим Ли, но не нашла практического применения.

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

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

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

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

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

Одним из основных инструментов построения законов сохранения физических процессов, допускающих вариационную формулировку, является первая теорема Эмми Нетер [53], которая устанавливает связь между инвариантностью вариационного интеграла относительно конечномерной группы Ли и дивергентными формами уравнения Эйлера. Эта теорема дает достаточное условие существования законов сохранения для уравнений Эйлера. Н.Х. Ибрагимов [19] дал новое доказательство этой теоремы для n-мерных интегралов на языке инфинитезимальных операторов Ли, что в сочетании с техникой группового анализа, развитого Л.В. Овсянниковым, дает возможность удобного способа построения законов сохранения. Используя понятие слабого лагранжиана, он также установил не только достаточные, но и необходимые условия существования законов сохранения уравнений Эйлера.

В первые годы развития теории пограничного слоя ученые пытались найти автомодельные решения, которые физически описывали бы некоторые частные течения. Фундаментальностью и глубиной отличаются труды А.А. Дородницына 1942-1948 гг. [16-18] по теории пограничного слоя в сжимаемом газе, определившие развитие этого раздела аэродинамики. Их идеи неразрывно связаны с природой сжимаемости и поэтому стали основой современных аналитических и численных методов расчета пограничного слоя в газе, включая самые сложные случаи с теплопередачей, излучением, протеканием равновесных и неравновесных физико-химических процессов. В работе [17] предложено преобразование переменных, которое теперь стало классическим и носит имя автора. Групповые свойства уравнений ламинарного пограничного слоя получили дальнейшее развитие в работах И.И. Пухначева [46], Ю.Н. Павловского [44,45], B.C. Каплан [22], К.Г. Гараева [1-6,8], В.Г. Павлова [6,8,13], С.А. Дербенева [12-15], В.А. Овчинникова [8,43]. С тех пор новые полученные автомодельные решения, имеющие физическую интерпретацию, автору неизвестны.

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

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

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

Приведем краткое содержание диссертационной работы.

Структура и объем диссертационной работы. Диссертация состоит из введения, трех глав, заключения, библиографического списка использованной литературы и приложений. Работа изложена на 84 страницах основного текста; иллюстративный материал представлен в виде 5 графиков; приложения содержат 33 таблицы; библиография включает 53 наименования.

ЗАКЛЮЧЕНИЕ

1. Получены новые первые интегралы уравнения Эйлера в простейшей задаче вариационного исчисления. построена дивергентная форма уравнения Эйлера-Остроградского для вариационной задачи с двумя независимыми переменными. Задача вычисления экстремального значения функционала по области V сведена к вычислению интеграла по поверхности S, ограничивающей эту область.

3. Дано приложение однопараметрической группы К вариационной задаче с подвижным концом.

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

Случай 1. Pr = 1; Ь(т) = 1 (линейная зависимость вязкости от температуры); т = const; безразмерная скорость на внешней границе пограничного слоя подчиняется специальному дифференциальному уравнению имеет решение, отличное от классического.

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

2. Основываясь на теории L*- инвариантности функционалов у с начальным условием а (0) = 0, которое е редуцируют исходные уравнения в уравнения, содержащие меньшее число искомых функций, нежели исходные. Использование аппарата Ли-Овсянникова позволило свести уравнения в частных производных к обыкновенным дифференциальным уравнениям, для которых формулируется соответствующая краевая задача.

Построены распределения безразмерной скорости ае(х) на внешней границе пограничного слоя и соответствующие им формы обтекаемого профиля при различных значениях параметра С* е (0;0.5].

Получены формулы для аэродинамических характеристик обтекаемого тела (касательного напряжения трения и локального теплового потока). Построены характерные графики. Случай 2. ?гф\; Ь(т) = const; г = const.

Условие инвариантности уравнений пограничного слоя относительно оператора приводит к постоянной скорости на внешней границе a - const, что соответствует случаю обтекания клина в сверхзвуковом потоке.

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

1. Гараев К.Г. Групповые свойства уравнений нестационарного пространственного пограничного слоя несжимаемой жидкости. //Труды КАИ. Казань. - 1970. - Вып.119. - С.47-53.

2. Гараев К.Г. Замечание к теории Нетер // Изв. высш. учеб. заведений, Математика. 1989. - №5. - С. 69-71.

3. Гараев К.Г. Группы Ли и теория Нетер в проблеме управления с приложениями к оптимальным задачам пограничного слоя. Казань: Изд-во КГТУ им. А.Н. Туполева. 1994. 240 с.

4. Гараев К.Г. Приложения непрерывных групп преобразований к дифференциальным уравнениям. М.: Соросовский Образовательный Журнал. - 1998. - №12. - С. 113-118.

5. Гараев К.Г., Дружинин Г.В., Павлов В.Г. Анализ автомодельности и расслоение уравнений нестационарного пограничного слоя на пластине методами теории групп Ли. // (Изв. высш. учеб. заведений). Авиационная техника. 1976. - №4. - С.27-30.

6. Гараев К.Г., Кусюмов А.Н., Павлов В.Г. Об управлении температурой поверхности сферы, обтекаемой высокоскоростным потоком вязкого газа. // (Изв. высш. учеб. заведений). Авиационная техника. 1987.- №2. - С.22-25.

7. Гараев К.Г., Овчинников В.А., Павлов В.Г. К задаче оптимизации теплообмена в ламинарном пограничном слое сжимаемого газа. //

8. Изв. высш. учеб. заведений). Авиационная техника . 1984.- №4. -С. 18-21.

9. Гельфанд И.М., Фомин С.В. Вариационное исчисление. М., Л.: Физматгиз, 1961. - 228 с.

10. Гинзбург И.П. Теория сопротивления и теплопередачи. JL: Изд-во ЛГУ. 1970.568 с.

11. П.Гошек И. Аэродинамика больших скоростей. М.: Изд-во иностр. лит-ры. 1954. 546 с.

12. Дербенев С.А. Групповые свойства уравнений пограничного слоя при наличии магнитного поля и химических реакций. // Казань: Труды КАИ. 1970. - Выпуск 119.

13. Дербенев С.А., Павлов В.Г. Инвариантно-групповые свойства уравнений пограничного слоя электропроводящей жидкости при наличии магнитного поля. // Казань: Труды КАИ. 1970. - Выпуск 125.

14. Дербенев С.А. Некоторые исследования уравнений пограничного слоя. // Казань: Автореферат диссертации. -1971.

15. Дербенев С.А. Групповые свойства уравнений пограничного слоя гиперзвукового потока газа. // Казань: Труды КАИ. 1972. - Выпуск 144.-С.82-86.

16. Дородницын А.А. Пограничный слой в сжимаемом газе // Прикл. матем. и механ. 1942. Т. 6. Вып. 6. С. 449-486.

17. Дородницын А.А. Ламинарный пограничный слой в сжимаемом газе // Сб. теоретических работ по аэродинамике. М.: Оборонгиз. 1957. -С.140-173.

18. Дородницын А.А. Об одном методе решения уравнений ламинарного пограничного слоя. // Прикл. математика техн. физика. 1960.-№3.-С.111-118.

19. Ибрагимов Н.Х. Инвариантные вариационные задачи и их законы сохранения. Теор. и мат. физика. - 1969. - T.I, №3. - С.350-369.

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

21. Ибрагимов Н.Х. Азбука группового анализа // Новое в жизни, науке и технике. Сер. Математика, кибернетика. -М: Знание. 1989. - №8. -48 с.

22. Каплан B.C. Групповые свойства и инвариантные решения уравнений стационарного трехмерного пограничного слоя. // Труды ЦАГИ. 1978. - Вып. 1857. - 79 с.

23. Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. -Т. II. М.: Физматгиз. - 1963.

24. Краснов Н.Ф. Аэродинамика. М.: Высшая школа. 1971. 630 с.

25. Ладыженская О.А. Метод конечных разностей в теории уравнений с частными производными // УМН. 1957. - Т.12, №5. - с.123-149.

26. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука. 1979. 847 с.

27. Лю-Шень-Цюань. Расчет ламинарного пограничного слоя в сжимаемом газе при наличии отсоса или вдува // Журнал вычислительной математики и математической физики. 1962. № 5. С. 868-883.

28. Ляшко И.И., Макаров В.Л., Скоробогатько А.А. Методы вычислений. (Численный анализ. Методы решения задач математической физики). Киев: Вища школа, 1977. 408 с.

29. Миеле А., Халд Д. Трехмерные тела минимального полного сопротивления // Теория оптимальных аэродинамических форм. -М.: Мир.-1969.-с.328-347.

30. Никифорова С.В. О существовании неклассических первых интегралов в простейшей задаче вариационного исчисления. // Казань: Материалы Всероссийской молодежной школы-конференции по математическому моделированию, алгебре и геометрии,- 1998.-С.25.

31. Никифорова С.В. Теоретико-групповой подход к исследованию одной задачи вариационного исчисления. // Казань: Вестник КГТУ. -1999.-№4.-С. 53-57.

32. Никифорова С.В. О достаточных условиях приводимости краевой задачи для уравнения Эйлера к задаче Коши. // Казань: Труды VIII Четаевской Международной конференции по аналитической механике, устойчивости и управлении движением. 2002.- С.353.

33. Никифорова С.В. О достаточных условиях сведения краевой задачи для уравнения Эйлера к задаче Коши // Казань: Вестник КГТУ. -2003.-№2.-С. 41-42.

34. Никифорова С.В. О новых фактор-системах уравнений пограничного слоя при сверхзвуковых режимах течения // Казань: Материалы Всероссийской молодежной школы-конференции по математическому моделированию, алгебре и геометрии. 2003. -С.29.

35. Никифорова С.В. О новых фактор-системах уравнений пограничного слоя при сверхзвуковых режимах течения // Казань: Вестник КГТУ. 2004. - №3. - С. 65-67.

36. Никифорова С.В. Определение формы профиля в задаче обтекания сверхзвуковым потоком газа // Труды II Всероссийской научной конференции. Самара. -2005. - Часть 2. - С. 187-189.

37. Никифорова С.В. Инвариантные решения уравнений пограничного слоя при сверхзвуковых режимах // Екатеринбург: Краткие сообщения XXV Российской школы по проблемам науки и технологий, посвященная 60-летию Победы. 2005. - С.45-47.

38. Овсянников JI.B. Группы и инвариантно-групповые решения дифференциальных уравнений. // Доклады Академии наук СССР. -1958. Т.118, N 3. - С.439-442.

39. Овсянников J1.B. Групповые свойства дифференциальных уравнений. Новосибирск: Изд-во СО АН СССР. - 1962. - 240 с.

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

41. Овчинников В.А. Устойчивость сдвиговых течений при переменных физических свойствах жидкости. // Казань: Автореферат диссертации. 1974.

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

43. Павловский Ю.Н. Групповые свойства управляемых динамических систем и фазовые организационные структуры. I. Группы, характеризующие динамические системы. // Вычислительная математика и физика. 1971. - №4. - С. 862-872.

44. Пухначев И.И. Групповые свойства уравнений Навье-Стокса в плоском случае. // Прикл.механика и техн.физика. 1960. - №1. -С.83-90.

45. Самарский А.А., Гулин А.В. Численные методы. М.: Наука, 1989. 432 с.

46. Тихонов А.И., Васильев А.Б., Свешникова А.Г. Дифференциальные уравнения. -М.: Наука. 1980.

47. ШлихтингГ. Теория пограничного слоя. М.: Наука, 1974. 711 с.

48. Hidman P.G. Generalited coordinate farms of governing fluid equations and associated geometrically juduced errors // AIAA Journ. №5. 1983. -P.47-57.

49. Lie S., Scheffers G. Vorlesungen uber Differentialgleichungen mit bekannten infinitesimalen Transformations gruppen. Leipzig, 1891. -568 s.

50. Lie S., Engel F. Theorie der Transformations gruppen. Bd. 1-3. Leipzig, Teubner, 1883-1893. - 638 s., 554 s., 830 s.