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

Вычислительные технологии

Том 7, № 2, 2002

ГРУППОВОЙ АНАЛИЗ ИНТЕГРОДИФФЕРЕНЦИАЛЬНЫХ КИНЕТИЧЕСКИХ УРАВНЕНИЙ. РЕЗУЛЬТАТЫ И ПЕРСПЕКТИВА

Ю. Н. Григорьев Институт вычислительных технологий СО РАН Новосибирск, Россия С. В. МЕЛЕШКО Институт теоретической и прикладной механики СО РАН

Новосибирск, Россия e-mail: grigor@ict.nsc.ru

The paper contains a survey of original authors' studies on group theory analysis of integrodifferential kinetic equations of Boltzmann, Smolukhovsiy, Vlasov and others. Scheme of direct method for calculation of widest Lie groups (algebras) admitted by these equations is given. The groups, optimal systems of their subgroups and new invariant solutions obtained by authors for concrete equations are considered.

Введение

Кинетические уравнения (КУ) составляют основу математического аппарата кинетических теорий разреженного газа, плазмы, переноса излучений, коагуляции дисперсных систем. Известны и стали классическими кинетические уравнения Больцмана, Власова, Ландау, Смолуховского [1-3].

Характерной особенностью КУ является наличие в них нелокальных (интегральных) операторов. В кинетические модели, где главную роль играют коллективные (осреднен-ные) взаимодействия частиц, нелокальность обычно входит в коэффициенты дифференциальных операторов через некоторые интегральные функционалы от решения. Примерами могут служить уравнения Власова для бесстолкновительной плазмы или уравнение Бхатнагара — Гросса — Крука в динамике разреженного газа.

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

© Ю.Н. Григорьев, С. В. Мелешко, 2002.

В математической теории КУ особое место всегда занимали исследования их инвариантных решений, которые напрямую связаны с фундаментальными свойствами симметрии этих уравнений. Как известно, для изучения инвариантных решений дифференциальных уравнений широко используется аппарат группового анализа [4, 5]. Выработанный здесь подход включает следующие этапы:

1. Нахождение тем или иным способом группы Ли, допускаемой уравнением.

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

3. Поиск и вычисление в явном виде отдельных инвариантных решений.

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

1. Нахождение допустимой группы Ли из определяющих уравнений

Для дифференциальных уравнений с локальными дифференциальными операторами допустимые группы Ли вычисляются явно путем построения и решения определяющих уравнений [4, 5]. Эту схему до работ авторов [6, 7] не удавалось непосредственно использовать для КУ из-за входящих в них нелокальных операторов.

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

Ф

где u — общий вектор переменных.

Постановка задачи предполагает, что в уравнение (1.1) входят операторы, действующие в бесконечномерных функциональных пространствах. Допускаемая уравнением (1.1) однопараметрическая группа преобразований G\(X) состоит из таких преобразований Ta, которые любое решение u уравнения (1.1) преобразуют в решение u' того же уравнения:

u' = TaU,

Ф

Каждое из преобразований Та £ 0\(Х) определяется значением группового параметра а. Через X обозначен инфинитезимальный оператор группы 0\, который в символической записи выглядит как

X = П

Здесь

ад = ^

— вектор проекций оператора X на координатные операторы ди. С помощью П(и) любое из преобразований (1.2) в малой окрестности а = 0 можно представить как инфинитези-мальное:

и' = и + П

Нелокальным операторам, входящим в уравнение (1.1), среди координатных векторов, входящих в ди, соответствуют производные Фреше. По определению [4] уравнение (1.1), допускающее группу С^Х), инвариантно относительно ее инфинитезимального оператора X. Это условие записывается в форме

X Ф

Если вычислить производную уравнения (1.3) по групповому параметру а в точке а = 0, то из соотношений (1.4)-(1.6) следует

^ = X Ф(и) = 0.

Таким образом, выражение

аФ(и') , ч

—Г1 =0

является эквивалентным (1.7) условием инвариантности. Рассматривая (1.8) на многообразии всех решений уравнения (1.1), получаем определяющие уравнения для однопарамет-рической группы С1^). Из них в принципе можно вычислить координаты П(и) оператора X. Если П(и) найдено, то допускаемая группа ) выстраивается интегрированием системы уравнений Ли [4]:

аи'

аа

Построенное таким образом множество всех различных однопараметрических подгрупп, которые допускает уравнение (1.1), образует искомую основную группу Ли.

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

— = П

2. Полные группы Ли и оптимальные системы подгрупп для частных кинетических уравнений

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

2.1. Кинетическое уравнение Больцмана с приближенным асимптотическим интегралом столкновений

Рассматривалось однопараметрическое семейство кинетических уравнений

1

Ф

0

Уравнение (2.1) приближенно описывает релаксацию однородного газа в области больших молекулярных энергий. Здесь х = г>2/2; /(х,*) — функция распределения (ФР) молекул по энергиям; 0(/) — некоторый функционал от /. Используется изотропная модель рассеяния со степенным потенциалом взаимодействия и = -1 [12]. Показатель а, который является параметром семейства, определен как а = (5 — V— 1)-1/2.

Уравнение (2.1) имеет смысл приближенной промежуточной асимптотики кинетического уравнения Больцмана при х ^ то и больших, но конечных Внимание к нему было привлечено в связи с проблемой расчета кинетических процессов при больших надтеп-ловых энергиях молекул. Несмотря на простоту (2.1), оно сохраняет характерные черты полного уравнения Больцмана, в том числе нелокальный интегральный оператор столкновений.

Приложение общей схемы к уравнению (2.1) выглядит следующим образом. Пусть и = (х, /) — вектор переменных. Инфинитезимальный оператор группы записывается в виде

X = £(х, /)дл + п(ж, /)д + С(х, /)д,. Определяющее уравнение (1.8) получается в форме

1

д Ф Г

— = хаА^(х, *) + ^(0, *)/(х, *) + /(х, *Шх, *) - 2 ^ /((1 - в)х, ¿Швх, *) = 0. (2.2) да |а=о ,/

0

Здесь /(х,*) — произвольное решение (2.1); Dt — оператор полного дифференцирования по

^(х, *) = С(х, /(х, *)) - £(х, /(х, / - п(х, /(х, *))/.

Решение полученного определяющего уравнения (2.2) строилось методом расщепления в классе формальных степенных рядов. Функции ^(х,*),/(х,*) представлялись рядами вида

^(х,*) = ^ г„(*)хп, /(х,*) = ^ £„(*)хга

п=0 п=0

Этот подход использовался нами и для других кинетических уравнений, рассматриваемых в данной работе. Характерные технические приемы построения решений определяющих уравнений можно найти в [9, 10]. Здесь мы опускаем достаточно длинные выкладки и ограничиваемся формулировкой полученных результатов.

Общее решение определяющего уравнения в данном случае имеет вид

С = (с1 + С2х + сза)/, £ = сзх, п = с4 - С1*,

где с^ с2, c3, c4 — произвольные постоянные. Часто вместо рассмотрения в терминах групп Ли удобно использовать эквивалентный язык алгебр Ли [4, 5]. Основные результаты группового анализа уравнения (2.1) можно сформулировать здесь следующим образом.

оо

оо

Предложение 2.1. Основная (наиболее широкая) алгебра Ли Ь4, допускаемая интег-родифференциальным уравнением (2.1), однозначно определяется базисными операторами

X! = д, Х2 = х/д/, Хз = а/д/ + хдж, Х4 = /д/ -

Оптимальная система одномерных подалгебр алгебры Ь4 порождается операторами

Хь Х2, Хз, Х4 + вХз, Х2 + X!, Хз + X!, Х4 ± Х2,

где в — произвольная постоянная.

Система (2.4) позволяет выделить все классы инвариантных решений уравнения (2.1), существенно различных относительно группы С4 в том смысле, что ни одно из них нельзя получить из другого с помощью преобразований Та £ С4.

Каждому из операторов (2.4) соответствует определенное представление инвариантных решений (2.1), вид которых приводится ниже:

Х! : / = £(у^ У = х

Хз : / = ха#(у), у =

Х4 + вХз : / = Х^(у), У = ,

Х2 + Х! : / = еХ£(у), у = х,

Хз + Х! : / = ха#(у), у = хе-4,

Х4 ± Х2 : / = Г(!±х)£(у), у = х.

Для построения фактор-уравнений необходимо вид этого решения подставить в уравнение (2.1) и учесть представление для независимых переменных. В данном случае это делается элементарно, и мы не будем на этом останавливаться.

2.2. Кинетическое уравнение Смолуховского

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

V оо

/^ = 2/^в(V - г^)/(г - гМ)/(гМ) - /М)^^(^)/(гь*), (2.5) 0 0

/ (г, 0) = /о (г).

Здесь /(г,£) — ФР частиц по объемам. Симметричная неотрицательная функция

в (г, г!) = в (г!,г) > 0

называется ядром уравнения коагуляции [3]. При коагуляции двух частиц образуется одна частица суммарного объема. Такой характер взаимодействия отражается в структуре аргументов в интегральном операторе (2.5).

Нами были найдены преобразования ФР f (г,£) и временной переменной позволяющие свести уравнение (2.5) к уравнению вида

1

т) = ж7 J ^

0

В частности, это возможно сделать для ядер вида

в

при некоторых значениях параметров а, Ь, более удовлетворительно описывающих взаимодействие коагулирующих частиц [3]. Основной результат, полученный нами для (2.5), состоит в следующем.

Предложение 2.2. Основная (наиболее широкая) алгебра Ли Ь4, допускаемая интег-родифференциальным уравнением (2.5), подобна алгебре ¿4 (2.3) для уравнения (2.1) и в переменных (2.6) определяется базисными операторами

Х1 = дТ, Х2 = жбд, Х3 = хдх — , Х4 = — тдТ.

Оптимальная система подалгебр, допускаемых уравнением (2.5), совпадает с системой (2.4).

Вид существенно различных инвариантных решений также совпадает с найденными для уравнения (2.1).

2.3. Кинетическое уравнение Власова

Кинетическое уравнение Власова рассматривалось в рамках задачи об одномерных высокочастотных колебаниях в плазме. Задача описывается следующей системой уравнений:

ft + v/x - E/v = 0, /t=t0 = /ü

Ev = 0, Et = J v/dv, Ex = 1 -J /dv,

где ФР /(x,v,t) и напряженность электрического поля E зависят от пространственной координаты x и скорости v, определенных в R1. Здесь кинетическое уравнение в приближении самосогласованного поля формально не содержит такого нелокального оператора, какими были интегралы столкновений в предыдущих случаях. В последние два уравнения нелокальные операторы входят в виде моментов ФР, взятых по скоростному пространству R1.

В [13] для системы (2.7) была найдена допускаемая ей алгебра Ли L5 с базисными операторами

Xi = dt, X2 = 5x, X3 = x5x + v5v + E&e - fdf,

X4 = cos tdx — sin + cos tdE , X5 = sin tdx + cos + sin tdE.

Она была построена путем перехода от (2.7) к бесконечной системе дифференциальных уравнений для моментов ФР /. При этом вопрос о полноте L5 и ее возможных расширениях остался нерешенным. Построив определяющие уравнения для алгебры Ли, допускаемой системой (2.7), и найдя их общее решение описанным выше способом, мы получили следующий результат.

Предложение 2.3. Алгебра Ли L5 (2.8) системы (2.7) полна и не допускает расширений.

2.4. Система кинетических уравнений Больцмана релаксирующего многокомпонентного газа

Система кинетических уравнений Больцмана, описывающих пространственно однородную эволюцию к равновесию N-компонентной смеси атомарных газов, рассматривалась для случая изотропной максвелловской модели молекулярного рассеяния [2]. При этом в терминах характеристических функций

для а = 1 , ..., N система кинетических уравнений записывается в виде

N 1

У. ^«в /

в=1

0

Инфинитезимальный оператор наиболее широкой группы Ли, допускаемой системой (2.9), разыскивался в форме

N

X = Г(и)34 + (и)дха + С ),

а=1

где и = (¿, XI, ... , зд, ^1, • • • , ) — общий вектор переменных. К определяющим уравнениям для коэффициентов оператора X присоединялись дифференциальные связи = 0, а = в, а, в = 1, ... , N. Эти связи давали дополнительную подсистему определяющих уравнений

Се+ Се - 4>«л$,р+ 4) - ¥>а,*в(Се+ Сз) = 0, а = в, а, в = 1, ..., N.

Анализ решений этой подсистемы в классе формальных рядов по степеням компонент вектора и позволил привести функциональную зависимость коэффициентов X к виду

е4 = т с = с(^,х«,£), с = ,*), а = 1,

С учетом этого система определяющих уравнений записывается так:

N 1

А^а + ^ ^«в (0) - ^ ^а(жа(1 - (Ха5)-

1 в=1 0

- ^

Здесь = (а - еа - А — оператор полного дифференцирования по ¿. Полу-

ченные нами результаты [7-9] для более простых случаев позволили предположить, что система (2.9) допускает алгебру Ли А, определяемую базисом операторов

N NN

Х1 = Х2 = ^ Ха, Хз = ^ , Х4 = ^ -

а=1 а=1 а=1

Инвариантность (2.9) относительно Х1, ..., Х4 проверяется непосредственно. Факторизация по ним существенно упрощает решение определяющих уравнений (2.10), которое также строилось в классе степенных рядов. Результаты можно сформулировать следующим образом.

Предложение 2.4. Наиболее широкая алгебра Ли, допускаемая системой (2.9), есть ¿4 с базисом операторов Х1, Х2, Х3, Х4.

Предложение 2.5. Классы инвариантных решений (2.9), существенно различных относительно Ь4, определяются оптимальной системой одномерных подалгебр с инфи-нитезимальными операторами

X, Х4 + вХз, X — X, Х4 ± Х2, X + Х3.

На этой основе было дано исчерпывающее описание класса инвариантных решений, обобщающего на случай многокомпонентной смеси известное точное решение Бобылева — Крука — Ву (БКВ) уравнения Больцмана для простого газа. Строгая формулировка этих результатов, данная в работах [10, 11], потребовала бы много места. Здесь же мы ограничимся следующим утверждением.

Предложение 2.6. При некоторых необходимых и достаточных условиях БКВ-ре-шения имеют вид нестационарных максвелловских ФР

/а(Иа,*) = Па[2пТ(¿)]-3/2 вхр ( — ^^ ,

а = 1, ..., /, I < г,

соответственно для компонент а = I + 1, ... , г

/„(^) = П„[2ПГ М]-3/2

1 + ^6<"> (-

3

Г (¿) у 2Г (¿) 2

а для остальных компонент а = г + 1, ..., N

ехр

2Г (^

/а К,£) = Па[2пГ (¿)]-3/2

1 +

1 — Ь(а)Па ,

1 — Г (¿) а=+1 а

Г(¿) " "

3

N

Е Па

а=г+1

\2Г (¿) 2

еХр

2Г(^ ,

где Ь(а), а = I + 1, ...,г — произвольно выбираемые коэффициенты, такие что У6(а) > 0,

г

Е Ь(а)Па < 1.

а=г+1

Найденные БКВ-решения с произвольными коэффициентами Ь(а) существенно расширяют множество БКВ-решений в элементарных функциях. Особый интерес представляют нестационарные максвелловские ФР, которые ранее в рамках БКВ-класса не рассматривались. Подбирая специальным образом параметры компонент и используя найденные решения, можно моделировать различные режимы релаксации. Интересные с физической точки зрения возможности такого моделирования обсуждались нами в работах [14, 15].

2

V

а

2

V

а

3. Классификация инвариантных решений полного уравнения Больцмана

Для полного уравнения Больцмана (УБ), в котором ФР зависит от семи переменных, а интеграл столкновений содержит пятикратное интегрирование, групповые свойства и

инвариантные решения до последнего времени изучались исключительно на основе эвристических ad Ное методов. При этом, как правило, форма допустимых преобразований постулировалась априори. В настоящее время в работах [6, 16, 17] и некоторых других найдено, что полное УБ допускает 11-параметрическую группу Ли Оц. Известны также некоторые расширения Оц для частных случаев межмолекулярных потенциалов. Естественно, что в силу эвристического подхода вопрос о полноте найденных групп остается до сих пор открытым. Тем не менее было бы очень полезно построить для Оц оптимальную систему подгрупп, что позволило бы выполнить классификацию множества инвариантных решений УБ и начать их систематическое исследование, начиная с простейших. В теории групп известны общие алгоритмы построения оптимальных систем, однако их практическая реализация для групп большой размерности, таких как Оц, связана с большим объемом рутинных вычислений. В качестве примера можно привести ситуацию с системой уравнений газовой динамики, которая для разных уравнений состояния допускает группы Ли подобной размерности. Несмотря на то, что эти полные группы найдены еще в 70-е годы, соответствующие оптимальные системы подгрупп были вычислены лишь в последнее время с использованием методов компьютерной алгебры [18, 19]. Оптимальные системы подгрупп известных групп Ли, допускаемых полным УБ, рассматривались нами в работах [20, 21].

3.1. Допустимая группа Ли полного уравнения Больцмана

Полное уравнение Больцмана, описывающее эволюцию ФР / (£, х, V) в прямом произведении пространств Я+ х Я х Я, имеет вид [2]

1 ^ I =1

}

V* = + w + дп), w* = + w — дп), и = V — w, д = IV — w| = IV* — w*|, |п| = 1. Здесь £ € Я+ — время; х = (ж, у, г) € — пространственные координаты; V = (и,^,^) €

?э _______________________________Л ёп

Я — молекулярные скорости; а д, — — дифференциальное сечение рассеяния. Для

V д / _ _

степенных потенциалов межмолекулярного взаимодействия и (г) ж г_(^-1) (V > 2) сечение

имеет вид а = д71 ( ), 7 = (V — 5)/(^ — 1). Значение V = 5 (7 = 0) соответствует

д

максвелловским молекулам, а предел V ^ то — молекулам-твердым сферам.

В работе [6] мы строили группу Ли О(Та) точечных преобразований Та для УБ (3.1) в виде

/ = х, V; а)/', £ = т

где а — групповой параметр, В — некоторая 3 х 3-матрица. Дополнительно предполагалось, что преобразование нелинейного интеграла столкновений обладает следующим свойством обобщенной "однородности":

J(/',/')= х', V')/(/,/).

Неизвестные функции в (3.1) находились из основного свойства допустимой группы Ли: для каждого значения параметра a преобразование Ta переводит произвольное решение УБ вновь в некоторое его решение. Подстановка (3.3) в уравнение (3.1) и использование метода расщепления [4] позволили вычислить допустимую группу G(Ta) в явном виде. Для дальнейшего нам необходима соответствующая найденной группе алгебра Ли инфи-нитезимальных операторов. Для случая произвольного сечения а имеет место алгебра Ли Ln(X) с базисом операторов:

Xi = dx, X2 = dy, X3 = dz, X4 = tdx + du, X5 = tdy + 0V, Хб = tdz + dw,

X7 = ydz — zdy + vdw — wdv, X8 = zdx — xdz + wdu — udw, X9 = xdy — y5x + u5v — vdu, X10 = dt, Xn = tdt + xdx + ydy + zd* — fdf.

Для степенных потенциалов алгебра L11(X) расширяется до алгебры L12(X) оператором

X12 = tdt — udu — vdv — w5w + (7 + 2) fcf. Наконец, для случая 7 = — 1 имеется дополнительный оператор

X13 = t2dt + txdx + (x — tv)dv,

который отвечает проективному преобразованию [4]. Действие производной df в операторах X11 и X12 надо рассматривать как дифференцирование по Фреше. Заметим, что в [17] группы Ли с приведенными здесь операторами были заявлены как полные группы, допускаемые полным уравнением Больцмана (3.1), хотя проведенные вычисления были выполнены практически тем же ad hoc методом, что ив [6]. В этой связи необходимо еще раз подчеркнуть, что доказательство полноты допустимой группы возможно только путем построения общего решения определяющих уравнений, как это было сделано в упоминаемых выше работах. Подобное доказательство для многомерной группы G11 еще предстоит провести.

3.2. Классификация подалгебр

Классификация подалгебр позволяет подразделить множество Н-решений УБ, инвариантных относительно группы С11, на классы эквивалентности. Два Н-решения /1 и /2 принадлежат одному классу эквивалентности, если существует преобразование Га £ О такое, что /2 = Га /1. В противном случае /1, /2 принадлежат разным классам и называются существенно различными Н-решениями. Список всех существенно различных Н-решений (по одному представителю от каждого класса) образует оптимальную систему инвариантных решений, определяющих искомую классификацию. Для ее нахождения строится оптимальная система подалгебр допустимой алгебры Ли Ь [4, 5]. По определению есть максимальное множество подалгебр N £ Ь, произвольная пара которых неподобна относительно внутренних автоморфизмов алгебры Ь. Для алгебр малой размерности вычисление достаточно просто, но с повышением размерности Ь вычислительные трудности возрастают многократно.

Однако для алгебры Ь11 (X), допускаемой УБ (3.1), имеет место обнаруженное нами замечательное обстоятельство, позволяющее избежать огромного объема вычислений. Справедливо следующее утверждение [20].

Предложение 3.1. Алгебра Ли Ьц(Х), допустимая полным уравнением Больцмана, изоморфна алгебре Ли ¿и (У), допускаемой системой уравнений газовой динамики

+ рУи = 0, р1Г + = 0, ^ + Жр.р^

0.

Здесь р, р — плотность и давление газа, как и выше; г € Я+

+ '

х

(ж, у, г) € Д; V

(м, € Д, только у — теперь вектор макроскопической скорости газа.

Как известно [4], для произвольной функции состояния А(р,р) система (3.4) допускает 11-параметрическую группу Ли с инфинитезимальными операторами

У = дх, У> = ду, Уз = д^, У4 = гдх + д«

У7 = удг - гду + - ид Уд = жду - удх + мд^ - -уд«, У1

Ув =

= д4,

, У5 = гду + д^, Уб = ¿д* + дад, гдх - жд^ + идад - мдад, У11 = + ждх + уду + гд^.

Доказательство изоморфизма следует из инвариантности коммутаторов обеих алгебр относительно отображения алгебры Ь11(Х) на алгебру Ь11(У), которое проверяется непосредственно.

Для политропного газа, когда А(р, р) = кр, существует расширение алгебры Ь11(У) двумя дополнительными операторами

У12 = гд4 - мд„ - -д^ - идад + 2рдр, У13 = рдр + рдр

до алгебры Ь13(У) = Ь12(У) Ф {У13}. Аналогично доказывается, что алгебра Ь12(Х) изоморфна подалгебре Ь12(У).

Так как оптимальная система подалгебр полностью определяется таблицей коммутаторов, которые инвариантны относительно изоморфизма, то для классификации и нахождения существенно различных Н-решений полного уравнения Больцмана (3.4) можно непосредственно использовать оптимальную систему подалгебр для уравнений газовой динамики, построенную в недавней работе [23]. Правда, есть определенные различия при ее использовании для нахождения существенно различных инвариантных решений уравнения Больцмана и системы уравнений газодинамики из-за разного количества неизвестных функций и независимых переменных. Наибольший интерес представляют решения с одной и двумя независимыми инвариантными переменными. Для уравнения (3.1) с произвольным сечением рассеяния оптимальная система подалгебр Ь11(У) [23] дает 11 различных классов инвариантных решений от одной независимой инвариантной переменной и 38 классов с двумя независимыми переменными [20, 21]. Найденные представления инвариантных решений с одной независимой переменной приведены в табл. 1.

Таблица 1

№ Представление ФР Индекс

1 С

2 Г1^2 +

3 С

4

5 Ь-1^

6 ¥>

7 у

8

9 е£Уф

10

11 ж-1^

Индекс (т, п) в последней колонке идентифицирует представителя оптимальной системы подалгебр: т есть размерность соответствующей подалгебры, а п — номер подалгебры в табл. 6 работы [23]. Буква С означает, что данное представление надо рассматривать в цилиндрической системе координат (х,г, V, Ж), q = VV2 + . Остальные представления соответствуют декартовым координатам. Можно видеть, что для большинства представлений инвариантные решения либо не существуют, либо не имеют физического смысла, хотя в последнем случае могут быть полезны для численных тестов.

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

3.3. Допустимая группа Ли полного уравнения Больцмана в Фурье-представлении

Анализ полученных представлений инвариантных решений полного УБ весьма затруднен сложной структурой интеграла столкновений. Это заставило нас в поисках возможного упрощения обратиться к групповому анализу уравнения Больцмана в Фурье-представлении. Для степенных потенциалов взаимодействия последнее имеет вид [24]

дф . д2ф д£ дхду

СЫ У 07 ^ф (ф (у -

м

--п

1 2

-

У

ТО^--

1 2

-

где и = у — ; С

27 п-3/2Г

Г

С использованием уже описанного метода для уравнения (3.5) была найдена допустимая алгебра Ли Ьд с базисными операторами

У = дх, У> = ду, У3 = дх,

У = удг — хду + — ид, У5 = гдх — хдг + ид — мдад,

Уб = хду — удх + мд^ — -уд«, У7 = д4, Уз = ¿д4 + хдж + уду + гд^ — ^д^,

Уд = ¿д4 + ид« + -д^ + + (7 — 1)^д,

Меньшая размерность найденной алгебры позволила непосредственно вычислить оптимальную систему подалгебр, применяя развитый для этого в [23] двухшаговый алгоритм. Как и выше, таким образом была получена классификация инвариантных решений уравнения (3.5). В частности, найдены представления решений с одной и двумя независимыми инвариантными переменными. Представления инвариантных решений с одной независимой переменной приведены в табл. 2.

Таблица 2

№ Представление ФР Индекс

1 7.2

2 7.3

3 з дав-4) 7.4

4 7.5

5 7.6

6 6.1

7 к7-1г-7 ^(г^/Л) 6.2

8 и7 ¿-1,0(д/-и) 6.3

9 6.4

Заметим, что хорошо известное БКВ-решение для простого газа принадлежит классу 7.4 при 7 = 0. Но теперь этот класс расширен на произвольные степенные потенциалы.

Заключение

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

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

Для таких уравнений представляется также возможным выполнить их групповую классификацию по отношению к "произвольному элементу"[4], т.е. найти все возможные расширения основной допустимой алгебры для конкретных спецификаций этого элемента. Например, это вполне достижимо и представляет интерес для уравнения Смолуховского, где таким элементом является ядро коагуляции.

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

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

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

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

Список литературы

[1] БАЛЕСКУ Р. Статистическая механика заряженных частиц. М.: Мир, 1984.

[2] ФЕРЦИГЕР Дж., КАПЕР Г. Математическая теория процессов переноса в газах. М.: Мир, 1976.

[3] ВОЛОЩУК В. М., СЕДУНОВ Ю. С. Процессы коагуляции в дисперсных системах. Л.: Гидрометеоиздат, 1975.

[4] ОВСЯННИКОВ Л. В. Групповой анализ дифференциальных уравнений. М.: Наука, 1978.

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

[6] ГРИГОРЬЕВ Ю.Н., Мелешко С. В. Исследование инвариантных решений кинетического уравнения Больцмана и его моделей. Новосибирск, 1986 (Препр. / АН СССР. Сиб. отд-ние. ИТПМ; №18-86).

[7] ГРИГОРЬЕВ Ю.Н., Мелешко С. В. Групповой анализ интегродифференциального уравнения Больцмана // Докл. АН СССР. 1987. Т. 297. C. 323-327.

[8] Grigoryev Yu.N., Meleshkü S.V. Group theoretical analysis of kinetic Boltzmann equation and its models // Arch. Mech. 1990. No. 6. P. 693-701.

[9] Grigoryev Yu.N., Meleshkü S.V. Group analysis of kinetic equations // Russ. J. Numer. Anal. Math. Modell. 1995. Vol. 10, No. 5. P. 425-447.

[10] ГРИГОРЬЕВ Ю.Н., Мелешко С. В. Полная группа Ли и инвариантные решения системы уравнений Больцмана многокомпонентной смеси газов // Сиб. мат. журн. 1997. Т. 38, №3. C. 511-525.

[11] ГРИГОРЬЕВ Ю. Н., Мелешко С. В. Групповой анализ уравнений Больцмана релак-сирующего многокомпонентного газа // Докл. РАН. 1998. Т. 360, №4. C. 459-462.

[12] Ernst M. H. Nonlinear model-Boltzmann equations and exact solutions // Phys. Rep. 1981. No. 78. P. 1-171.

[13] ТАРАНОВ В. Б. О симметрии одномерных высокочастотных движений бесстолкнови-тельной плазмы // ЖТФ. 1982. №46. C. 1271-1277.

[14] Grigoryev Yu. N., Meleshko S. V. Bobylev — Krook — Wu-modes for multicomponent gas mixtures // Phys. Rev. Letters. 1998. Vol. 81, No. 1. P. 93-95.

[15] Grigoryev Yu.N. Bobylev — Krook — Wu-modes: explicit form and applications for modeling // Book of Abstracts of 21st Intern. Symp. on Rarefied Gas Dynamics, Marseille, France, July 26-31, 1998. Pt II. P. 33-34.

[16] Бунимович А. И., КРАСНОСЛОБОДЦЕВ А. В. Инвариантно-групповые решения кинетических уравнений // Изв. АН СССР. Механика жидкости и газа. 1982. №4. C. 135140.

[17] Bobylev A. V. Boltzmann equation and group transformations // Math. Models Methods Appl. Sci. 1993. No. 3. P. 443-476.

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

[19] ЧЕРЕВКО А. А. Оптимальная система подалгебр для алгебры Ли операторов, допустимых уравнениями газовой динамики с уравнением состояния p = f (S)р5/3. Новосибирск, 1996. (Препр. / РАН. Сиб. отд-ние. Институт гидродинамики им. М. А. Лаврентьева; №4-96).

[20] Grigoryev Yu. N., Meleshko S. v., Sattayatham P. Classification of the Invariant Solutions of the Boltzmann Equation // Phys. A: Math. and General.1999. No. 32. P. L337-L343.

[21] Grigoryev Y. N., Meleshko S.V. Classification of invariant solutions of the full Boltzmann equation // Proc. Xth Intern. Conf. on Waves and Stability in Continuous Media (WASCOM 99), 7-12th June, Italy, 1999.

[22] Grigoryev Y. N., Meleshko S.V. Group Classification and Representations of Invariant Solutions of the full Boltzmann Equation // Proc. 22nd Intern. Symp. on Rarefied Gas Dynamics, Sydney, Australia, 9-14 July 2000.

[23] Овсянников Л. В. Программа "ПОДМОДЕЛИ". Газовая динамика // Прикладная математика и механика. 1994. №58. C. 30-55.

[24] Григорьев Ю.Н., МИХАЛИЦЫН А. Н. Спектральный метод численного решения кинетического уравнения Больцмана // Журн. вычисл. математики и мат. физики. 1983. Т. 23, №6. C. 1454-1463.

Поступила в редакцию 23 ноября 2001 г.