О гиперболических регуляризациях законов сохранения Палин Владимир Владимирович

Введение к работе

Диссертация посвящена исследованию поведения на больших временах решений задачи Копти для гиперболических регуляризации законов сохранения (в терминологии С. Bardos) или систем законов сохранения с релаксацией (в терминологии Gui-Qiang Chen, Levermore С. D) см. например [7]. В одной из последних работ Максвелла была поставлена задача вывода уравнений гидродинамики из кинетической теории газов, т.е. системы уравнений Навье-Стокса сжимаемой жидкости из кинетического уравнения. Л. Больц-ман в статье "О максвелловском методе вывода уравнений гидродинамики из кинетической теории газов" делает следующее предположение: "явствует с очевидностью, что незадолго до смерти он, должно быть, предпринял длинное и детально разработанное исследование этого вопроса, которое, однако, не было опубликовано." Насколько известно, если это исследование и существовало, то так и осталось в виде рукописей. Важность подхода Максвелла была понятна специалистам высокого уровня, каким, безусловно, был Больц-ман, по могла быть не оценена "широкой научной общественостыо". Повторить вычисления, о которых говорит Больцман, физикам, по видимому, в то время не удалось, поскольку первая работа Чепмена, в которой ставилась та же задача и проводился анализ интеграла столкновений появилась спустя 28 лет после указанной статьи Больцмана (и через 43 года после самой работы Максвелла). Задача вывода уравнений гидродинамики из кинетической теории газов определила математические проблемы кинетики для гиперболических регуляризации законов сохранения:

d_tUi + dw_x f\u, v) =0, і = 1,..., тп, (1)

d_tv_k + div-e g^k{u,v) + b^k(u)v = 0, k = m+l,...,N.

Здесь x Є R^d, u{x, t) : R^d xl₊^ R^m, v(x,t) : R^d xl₊^ R^N~^m, b - матрица релаксации порядка (N — m) x (N — m), потоки

f{u, v)eR^d, і = 1,..., m; g^k{u, v) eR^d, к = 1,..., N- m,

и - консервативные переменные, v - неравновесные переменные, т - число консервативных переменных . Главная часть системы (1) - нестрого гиперболична:

Определение. Система называется нестрого гиперболической, если характеристическая матрица

(2)

V 9u{u,v) g_v(u,v)

имеет только вещественные (возмо-жно кратные) корни т = Tj(,,u,v), j = 1,..., N.

Это условие выполнено, если система (1) - симметризуема. Примерами таких систем являются момонтные аппроксимации кинетических уравнений (например, кинетического уравнения Больцмана, описывающего неравновесные процессы гидродинамики, Фоккера-Планка, описывающего динамику броуновских частин,, Больцмана-Пайерлса, опиоывющего процессы теплопе-реноса в кристалах), простейнтая гиперболическая регуляризация трением системы уравнений Эйлера изептропической модели газовой динамики [21], расширение Дирака-Швиндлера системы уравнений Максвелла и т.д.

Изучению гиперболических регуляризации систем законов сохранения (системы законов сохранения с релаксацией) посвящено много работ как российских, так и зарубежных авторов. Прежде всего это касается исследования феномена релаксации [2], в частности вопросов устойчивости и сингулярного предела при стремлении времени релаксации к нулю в работах С. Bardos, С. D. Levermore [1], [2], R. Е. Caffish, G. С. Papanicolaou [4], G. Q. Chen, Н. Frid [8], Е. Yu. Panov [16] и т.д.. С исследованиями Чепмена (см.[5], [6]), порожденными проблемой Максвелла о выводе уравнений гидродинамики из кинетического уравнения, связана задача о существовании "промежуточного аттрактора" для систем вида (1). При выводе уравнений гидродинамики из кинетической теории весьма существенно получить простую функциональную зависимость коэффициентов переноса от вида потенциала взаимодействия и, тем самым, упростить анализ получающихся уравнений. Чепмен и Энског в работе [7] выдвинули гипотезу, что для "физически корректных" моделей механики сплошных сред влияние моментов высшего порядка, т.е. большей части переменных г>, "несущественно". Они предположили существование операторного уравнения состояния

v = Qu_: (3)

выражающее неравновесные переменные через консервативные (проекция в фазовое пространство консервативных переменных), которое замыкает систему законов сохранения в(1)

dtW + d_xf(w,Qu{w)) = 0. (4)

так, что специальные решения UchEns = {и = w, v = Qw} задачи Копій для системы (1), определяемые решениями w системы проекции (4), образуют инвариантное притягивающее многообразие MchEns (промежуточный аттрактор). Т.е. для любого решения U = (u,v) задачи Коти (1) с начальными данными /|г=о ⁼ {u,v) существуют начальные данные wq — T(t/,o,vo) для системы проекции (4), так что в некоторой норме невязка U — UchEns между U и специальным решением UchEns — (w, Оли) стремиться к нулю при t —» оо. Более того, если в фазовом пространстве консервативных переменных

w —» 0, когда I —> оо,

то невязка С/я стремится к нулю быстрее чем UchEns- Существенно то, что мы остаемся в классе гиперболических систем с релаксацией (возможно псев-додифференциальньтх). Целью диссертации является обоснование гипотезы Чепмена-Энскога для некоторого класса систем вида (1).

Много попыток обоснования гипотезы Чепмена-Энскога связано с построением формальных асимптотических разложений для операторного уравнения (3) для систем (1) с жесткой релаксацией, когда матрица B(u)v во втором уравнении в (1) заменяется на -B(u)v, eCl (см. например литературу в монографии Gorban A. N., Karlin I. V. [12]). Однако, при таком подходе возникла проблема устойчивости получаемых приближений. Первое приближение (порядка є¹), названное Навье-Стокс приближением, подтвердило гипотезу -для систем моментов Грэда кинетического уравнения Больцмана оно в точности совпало с системой Навье-Стокса сжимаемой жидкости, более того, необходимым и достаточным условием правильного знака вязкости явилось условие устойчивости линеаризованных на состояниях равновесия систем моментов. Однако, следующие приближения оказались неустойчивыми (см. Бобылев А. В. [23]), несмотря на устойчивость исходной системы моментов, что привело к проблемам в обосновании асимптотических раложений. Впервые обоснование гипотезы Чепмена-Энскога было получено Gui-Qiang Chen, Levermore С. D. and Tai-Ping Luui [2](1994) для одномерных систем размера,

2x2 пида (1) с жесткой релаксацией (как иллюстрация общего результата были рассмотрены р— система законов сохранения с жесткой релаксацией и система уравнений паводковой воды). При выполнении условия устойчивости Лакса на характеристические скорости уравнения проекции и характеристические скорости исходной системы в этой работе было доказано, что невязка между ретцением задачи Копти для исходной системы и соответствующим специальным решением в слабом смысле стремится к нулю при стремлении к нулю времени релаксации. Задача рассматривалась на конечном временном интервале [О.Т]. и предполагалось, что время релаксации стремится к нулю. В работе J-F. Coulombei и Т. Goudon [9)(2004) глобальные гладкие решения многомерной изотермической системы уравнений Эйлера с жесткой релаксацией были построены также на конечном временном интервале [О, Т]. Показано, что при стремлении времени релаксации к нулю плотность сходится к решению уравнения теплопроводности. В терминологии проектора (3) здесь существует проекция в фазовое пространство консервативной переменной - плотности. В этом случае Навье-Стокс приближение совпадает с уравнением теплопроводности, но построения проектора и инвариантного притягивающего многообразия специальных решений UchEns ⁼ (w, Qw) в работе [9] не проведено, хотя при оценке невязки между плотностью и соответствующим решением уравнения теплопроводности неявно используется некоторая аппроксимация оператора (3). Позднее (2005-2007). в работах Е.В. Радкевича [8], [22], [32] для некоторых конкретных моделей, например, систем моментов не выше третьего порядка и систем моментов не выше четвертого порядка газа фононов для глобальных гладких решений задачи Копти (на бесконечном временном полуинтервале [0, со)) построен проектор (3). Проблема существования операторного уравнения (3) сведена к задаче о существовании гладкого решения параметрической системы алгебраических уравнений для символа псевдодифференциального по пространственным переменным оператора Q. Условия разрешимости этой системы алгебраических уравнений определяются соотношением времен релаксации. Оказалось, что коэффициенты символа матричного оператора Q являются функциями с внутренним слоем (типа кинка), плохо приближаемыми полиномами. Последнее объясняет природу неустойчивости пост-Навье-Стокс приближений. Более того, показано, что неприводимые проекции (3) определяют на больших временах основые динамики моделируемого неравновесного процесса.

В диссертации впервые для общей системы (1) проведен линейный анализ

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

Используемые обозначения. Опишем обозначения, которые будут использоваться в дальнейшем. Буквой Е будем обозначать единичную матрицу. Под Єї, ..., Єдг будем понимать единичные векторы-столбцы линейного пространства C^jV, т. е. вк - к-тьт столбец матрицы Е. Линейную оболочку множества векторов г?і, ..., v_m (не обязательно линейно независимых) будем обозначать Linjvi, ...,v_m}. Матрицу алгебраических дополнений к матрице А будем обозначать cof(A). Определитель матрицы А будем обозначать del,(А) либо |Л|. Во всех случаях, когда размеры матриц, участвующих в некотором выражении, не указаны ранее явно, будем предполагать, что они таковы, что операции матричного сложения и умножения в этом выражении определены корректно.

Основные результаты.

В первой главе сформулирована математическая формализация гипотезы Чепмена-Энскога о существовании проекции для линеаризованных задач. В первом параграфе показано, что существование проектора (операторного уравнения состояния (3)) эквивалентно разрешимости квадратного матричного уравнения

Ai (Ли + Л12Р21) = Л₂і + Л₂₂ Ai,

где Л = Ai^+ В, матрица Л - квадратная размера N х iV, матрица Ли - квадратная размера т х т, остальные блоки матрицы Л имеют соответствующий размер.

Основным результатом второго параграфа нерпой главы является следующая теорема:

Теорема. Пусть |Л| ф 0. Пусть, кроме того, найдутся векторы Vi,..,v_m такие, что: 1. V = Lin{vj}"¹ - собственное подпространство .матрицы Л, т.е. AV = V.

2. Вект,оры. vi,.., v_m, e_m и, -, &n образуют базис. Тогда квадратное матричное уравнение

Лп(Лц + Л12Р21) = Л₂і + Л₂оАі

разрешимо и наоборот.

Третий параграф первой главы посвящен обобщению этой теоремы на случай det(A) = 0. Основным результатом этого параграфа является следующая теорема:

Теорема. Пусть система векторов v\,...,v_m задает базис линейного пространства V, инвариантного для лштрицы А и такого, что Lin{vi, ...,iv_n,e_m+h ---5671} - базис W¹. Запишем эти векторы но столбцам,

в виде матрицы I „ j. Тогда соответствующее этому набору векторов

(с учетом порядка,) решение матричного уравнения, записывается в виде

* 21 ⁼ ^21^11

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

В четвертом параграфе первой главы исследуется вопрос о том, при каких условиях построенные в предыдущих параграфах решения квадратного матричного уравнения непрерывно зависят от параметра . Основным результатом этого параграфа является следующая теорема: Теорема. Пусть матрица А непрерывно зависит от, параметра ₇ обратима 'при всех ^ Є Но и ее собственные значения являются, однократными V ф Е*_; где мнооїсество Е* конечно. Тогда для того, чтобы квадратное матричное уравнение, соответствующее матрице А, имело решение, непрерывно зависящее от параметра . необходимо и достаточно, чтобы существовало собственное подпространство V размерности т такое, что Lin{V, e_m+i,..., е^} = C^N при всех (єЕ₀.

В пятом параграфе получено условие разрешимости матричного уравнения Ляпунова АХ — ХВ = С более слабое, чем условие несовпадения множеств собственных значений матриц А и В. С помощью этого результата доказано необходимое и достаточное условие полного разделения динамик, т.е. приведения матрицы Л к блочно-диагональному виду:

Теорема. Пусть матрица Л обратима, и существует базис Vi,...,v_m собственного подпространства V тлкого, что Linjui,..., v_m, е_т+\,..., е^} — ~R^N, V не может быть расширено до собственного подпространства матрицы Л размерности m + 1 путем добавления к базису Vi,...,v_m присоединенного вектора матрицы, Л. Тогда существуют мат,рицы, Ро\_; Q\o такие, чт,о

( Е -Q_n \( Е 0 \ / Е 0\(EQ₁₂\(M_u 0 \ V О Е )\ -Р<21 Е ) \ Р₂₁ Е ) \ О Е ) V 0 Моо )

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

Шестой параграф первой главы содержит пример исследования линеаризованной системы с использованием теоремы о полном разделении динамик.

В седьмом параграфе первой главы с помощью теоремы о полном разделении динамик выписывается представление решения задачи Коши для системы уравнений

д,и + Лд_хи + Ви = 0

в виде суммы трех слагаемых. После замены переменных в образах Фурье U — S~¹u, и если существует полное разделение динамик, то решение может быть представлено в виде

U = U_Ch + Ucor + U_H.

Далее сформулировано основное определение:

Определение. Будем говорить, чгпо проектор Р удовлетворяет условию Ьо-корректности по Чепмену-Энскогу в классе 7i = {(ZYo, Vo)} начальных данных, если для всех начальных данных (Uq, Vo) Є 7Ї найдется Т₀ > О такое, что для всех t > Tq выполнена оценка

Ш)~ ^{е ,г>т<}"

где К, 5 > 0 - некоторые константы.
\ Восьмой параграф первой главы посвящен исследованию вопроса о том,

какие условия на начальные данные матрицу Л необходимо наложить для того, чтобы существовал проектор Р в пространство консолидированных переменных, удовлетворяющий условию Ьг-корректности по Чепмену-Энскогу

в некотором классе ТС начальных данных. Для того, чтобы сформулировать ответ на этот вопрос более кратко, вводится следующее условие: Условие. Будем говорить, что для, пары множ:еств Гі(), Го() выполнено 'жесткое условие щели, если

3₇ > 0 : /о(Г₂) - Ьо(Гі) > 7.

Основным результатом этого параграфа (и одним из центральных результатов диссертации) является следующая теорема:

Теорема. Пусть .матрица Л т.акова, что найдется такое ко > О, что V : || > ко все собственные значения матрицы. Л являются, алгебраически однократными и па них выполнена оценка

|А(оі<а(і + іеіА

где С\, d\ - некоторые конслпанты. Пусть, кроме того, Г і - множество всех собственных значений матрицы Л_; соответствующих задающему разделение динамик собственному тюдпространству V, Го - все остальные собственные значения К, и для Т\, Го выполнено жесткое условие щели. Тогда, если (Uq, Vo) - образы Фурье от начальных данных - принадлежат множеству

H={(U₀,Vo): ||ЫЬ|| ^0, (1 + |ЄІ)⁵"^л|М₂₂|^^л/22ЬіуоЄІ,₂(М)},

т.е. начальные данные для, неравновесных переменных достаточно гладкие, а начальные данные для консервативных переменных не равны, нулю, то соответствующий разделению динамик^ проектор Р удовлетворяет, условию Ьо-корректности по Чепмену-Эискогу в классе ТС начальных данных. При этом для констант, К, 5 из определения Lo-корректности по Чепмепу-Энскогу верно следующее: К зависит от \\Uq\\, ||Vo||; 6 зависит, от 7 и свойств собственных подпространств матрицы Л.

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

Вторая глава посвящена переносу описанной для линейного случая техники на нелинейные уравнения. А именно, для систем вида

где f\(v) и fo(v,w) - нелинейные члены, /(0,0) = 0, пока/зано, что (при условии су и чествования глобального решения из класса С([0, +оо), Н²) П С¹ ([0, +оо), II¹)) для достаточно малых начальных данных ф є Н² из существования проекции для линеаризованной задачи следует существование проекции для нелинейной задачи. Во втором параграфе с помощью метода последовательных приближений показывается существование ретпений задачи (5) для малых начальных данных. Третий и четвертый параграфы этой главы посвящены построению нелинейного оператора проекции Чепмена-Энскога и изучению его свойств. Принципиально в диссертации вопрос существования проекции Чепмена-Энскога рассматривается для глобальных гладких траекторий, что приводит к условию существования достаточно широкой щели. Условие вырожденной щели приводит к необходимости исследовать задачу для глобальных слабых решений с дополнительной гладкостью (например, для глобальных слабых решений в пространствах Орлича).

Третья глава посвящена использованию описанной в первых двух главах техники для исследования ряда конкретных физических задач а также содержит различные дополнительные факты, связанные с рассматриваемым классом систем. Так, в первом параграфе изучается гиперболическая регуляризация системы уравнений Максвелла [27, ?, 14].

Второй параграф посвящен исследованию двумерной 13-моментной аппроксимации для кинетического уравнения Больцмана [8].

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

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

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

Шестой параграф третьей главы посвящен формулировке условий, при ко-

торых элементы матрицы оператора проекции Р в образах Фурье могут быть представлены в виде (*)^аР(І|²)> ^ГДЄ а; Є Z2, р - ветцественнозначная функция, что позволяет упростить линейный анализ проблемы І^-корректнооти по Чепмену-Энскогу. Тем более, такую структуру символов имеют многие задачи механики сплошных сред.

Благодарности. Автор искренне благодарен своему научному руководителю, доктору физико-математических наук профессору Евгению Владимировичу Радкевичу за постановку задач и постоянное внимание к работе, а также доктору физико-математических наук профессору Ирине Викторовне Асташовой за полезные замечания.

Похожие диссертации на О гиперболических регуляризациях законов сохранения

Дискретные уравнения вида dun/dt = F(un-1, un, un+1) (n C Z) с бесконечным набором локальных законов сохранения Ямилов Равиль Исламович

Обобщенные решения задачи Копи для квазилинейных законов сохраненияПанов, Евгений Юрьевич

Системы законов сохранения класса ТемпляАгафонов, Сергей Иванович

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

Асимптотическое по времени поведение решения начальной задачи Коши для закона сохранения с нелинейной дивергентной вязкостьюГасников Александр Владимирович

Неограниченные решения скалярных законов сохраненияЛысухо, Полина Валерьевна