Серия «Математика»
2016. Т. 18. С. 93—109
Онлайн-доступ к журналу: http://isu.ru/izvestia
УДК 517.946 ЫБО 35К05
Исследование совместности переопределенной системы для многомерного уравнения
О О ____/ О и \ Ф
нелинейной теплопроводности (частный случай)*
Г. А. Рудых
Иркутский государственный университет
Э. И. Семенов
Институт динамики систем и теории управления имени В. М. Матросова СО
Аннотация. В статье исследуется многомерное уравнение нелинейной теплопроводности. Это уравнение представлено в виде переопределенной системы дифференциальных уравнений с частными производными (число уравнений больше числа искомых функций). Как известно, переопределенная система дифференциальных уравнений может быть несовместной, у нее может не существовать ни одного решения. Поэтому для установления факта существования решений и степени их произвола проводится анализ данной переопределенной системы дифференциальных уравнений. В итоге проведенного исследования получены не только достаточные, но и необходимые и достаточные условия совместности переопределенной системы дифференциальных уравнений с частными производными. На основе этих результатов с использованием уравнения Лиувилля и теоремы о необходимом и достаточном условии потенциальности векторного поля излагается подход, позволяющий в ряде случаев конструктивно построить точные неотрицательные решения многомерного уравнения нелинейной теплопроводности с конечной скоростью распространения возмущений. Среди построенных точных решений имеются и такие, которые не являются инвариантными с точки зрения групп точечных преобразований и групп Ли-Беклунда. Особое внимание уделено уравнению со степенным коэффициентом нелинейной теплопроводности. Это уравнение является квазилинейным параболическим уравнением с неявным вырождением. Данное уравнение из параболического дифференциального уравнения второго порядка вырождается в нелинейное эволюционное уравнение первого порядка типа Гамильтона - Якоби.
* Работа выполнена при частичной финансовой поддержке РФФИ (проект № 15-08-06680) и Совета по грантам Президента Российской Федерации для государственной поддержки ведущих научных школ Российской Федерации (НШ-8081.2016.9).
РАН
Ключевые слова: многомерное уравнение нелинейной теплопроводности, конечная скорость распространения возмущений, точные неотрицательные решения.
1. Введение
В работе на основе уравнения Лиувилля [27; 30] и теоремы о потенциальных операторах [3; 26] излагается подход, позволяющий конструктивно построить новые точные неотрицательные решения многомерного уравнения нелинейной теплопроводности с конечной скоростью распространения возмущений [10; 5; 19]. При этом некоторые построенные решения не являются инвариантными с точки зрения групп точечных преобразований и групп Ли-Беклунда [14; 7].
2. Многомерное уравнение нелинейной теплопроводности с конечной скоростью распространения возмущений
Рассмотрим многомерное уравнение нелинейной теплопроводности
ut = V- (K(u)Vu) , u = u(x, t) : Q x R+ ^ R+, (2.1)
где x e Йп; Q С Йп — область; R+ = (0, +то); Й+ = [0, +то); u(x,t) > 0 — функция характеризующая температуру среды; K(u) — функция определенная при всех u e Й + ; K(u) > 0 при u > 0; K(0) > 0; K(u) e C2 (Й+) P| C (Й+) — коэффициент нелинейной теплопроводности. Уравнение (2.1) возникает во многих прикладных задачах и принадлежит к классу так называемых неявно вырождающихся параболических уравнений. Строгая математическая теория этих уравнений берет свое начало в сравнительно недавних работах [1; 10; 24; 25; 28; 31]. Всюду в этой работе относительно функции K(u) будем предполагать, что
1
ос. (2.2)
о
Известно [15; 8; 9], что сходимость интеграла (2.2) является необходимым и достаточным условием конечной скорости распространения возмущений в процессах, описываемых уравнением (2.1). Далее, пусть
u(x, 0) = uo(x), x e Йп,
/ n \ 1/2
uo(x) = 0, |x| = ¡Y^ xk) >r e Й+, (2.3)
т.е. u0(x) равна нулю вне замкнутого шара S(r) = {x e Йп : |x| < r < и является финитной [4] функцией с носителем supp uo(x) =
{х £ Мга : ио(х) = 0}. С другой стороны, так как исходя из специфики уравнения (2.1) мы рассматриваем лишь неотрицательные решения, то
Причем для финитности функции uo(x) необходимо и достаточно, чтобы множество suppu0(x) было ограниченным. Иначе mes supp u0(x) < +œ>, что в силу формулы (2.2) обеспечивает финитность по x G Rn решения u(x,t) задачи Коши (2.1), (2.3). Тем самым [13], в задаче Коши (2.1), (2.3) скорость распространения возмущений является конечной, если в любой момент времени t G [0, T] существует замкнутый шар S(r(t)) = {x G Mn : |x| < r(t) < и для всех t G [0,T] вы-
полняется включение suppu0(x,t) С S(r(t)). Кроме того сходимость интеграла (2.2) является [19] необходимым и достаточным условием финитности по x G Mn решения u(x,t) задачи Коши (2.1), (2.3). Итак, если K(u) G C2 (M+) П C (M+) и выполняется неравенство (2.2), то уравнение (2.1) помимо процессов нестационарной фильтрации описывает диффузию и распространение тепла в сплошной среде с большими температурными перепадами [10]. В частности, хорошо известным свойством уравнений типа нестационарной фильтрации является конечная скорость изменения носителей их решений. Первые общие результаты в этом направлении были получены в работах [15; 8; 9].
Теперь представим исследуемое уравнение (2.1) в виде следующей эквивалентной системы
Система уравнений (2.4), (2.5) является переопределенной (число уравнений превосходит число искомых функций, подлежащих определению) относительно и(х,Ь), где {(х,Ь) £ Мга — достаточно гладкая по х £ Мга, Ь £ М вектор-функция. Ниже будем предполагать, что Г(х,Ь) £ С1 (С ^ Мга), где С С Мга+1 — открытое множество; С = П х I; I = {Ь : 0 < Ь0 < +го}; П — проекция С в Мга.
Ясно, что у переопределенной системы уравнений (2.4), (2.5) может вообще не существовать решений. Поэтому для установления факта существования решений и степени их произвола необходимо провести исследование и анализ совместности переопределенной системы (2.4), (2.5). Отметим, что методы исследования совместности переопределенных систем дифференциальных уравнений с частными производными изложены в работах [23; 16; 20; 22]. С другой стороны, уравнение нелинейной теплопроводности (2.1) и линейное многомерное уравнение Лиувилля (2.4) можно рассматривать как систему с дифференциальной связью (2.5). Причем метод дифференциальных связей [16; 20; 22]
suppu0(x) = {x G Mn : u0(x) > 0}.
ut + V- (uf(x,t)) = 0,
(2.4)
(2.5)
существенно использует теорию совместности систем дифференциальных уравнений. Наконец, в работах [11; 12] рассматривалось приложение этого метода для уравнения (2.1) со степенным коэффициентом теплопроводности.
Соотношение (2.4) при заданной {(х,Ь) € С1 (С — Мп) является уравнением Лиувилля [27; 30] относительно температуры и(х, Ь) для системы обыкновенных дифференциальных уравнений (ОДУ)
й
х = £(х,¿), хей", х = —х(£). В силу формулы (2.2) имеет смысл функция
и
Ф(и) = Ф(0) = 0. (2.6)
о
При этом в силу монотонности Ф(и) существует обратная функция Ф-1(п), причем 0 < п < Ф(гс>) <
Ниже при определенных предположениях относительно вектор-функции {(х,Ь) € С1 (С — Мп) и с использованием теоремы о необходимых и достаточных условиях потенциальности векторного поля вектор-функции {(х,Ь) € Мп будет доказано, что скалярная функция и(х,Ь), определяемая из соотношения (2.5) и удовлетворяющая уравнению Лиувилля (2.4) является точным неотрицательным решением многомерного уравнения нелинейной теплопроводности (2.1).
3. Основные результаты
Пусть {(х,Ь) = —А(Ь)х — В(Ь), где А(Ь) = [а^(Ь)] — вещественная симметричная их п матрица; В(Ь) = (Ь1(Ь),... ,Ьп(Ь))' — вектор-столбец; ац(Ь), Ьг(Ь) € С1 (М+); г,] = 1,2,... ,п. Тогда переопределенная система уравнений (2.4), (2.5) примет следующий вид
щ = и ■ ЬгА(Ь) + (Уи, А(Ь)х + В(Ь)), (3.1)
К(и)Уи = (А(Ь)х + В(Ь)) ■ и, (3.2)
п
где ^А(Ь) = ^ ац(Ь) — след матрицы А(Ь); (■, ■) — скалярное произ-
г=1
ведение в Мп. Поскольку матрица А(Ь) является симметричной, т.е. ач (Ь) = а^(Ь), тогда по теореме о потенциальных операторах [3; 26] уравнение (3.2) в силу (2.2), (2.6) обладает решением
и(х,Ь) = Ф"^[г(х,Ь)]+), (3.3)
г(х, *) = ± (х, А(*)х) + (х, В(*)) + С^),
(3.4)
где [•]+ = шах{[-|, 0}; Ф_1(г) — функция обратная к Ф(и), существующая в силу монотонности последней. Матрица А(Ь), вектор-столбец В(Ь) и скалярная функция С(Ь) подлежат определению, причем а^ (Ь), Ъг(Ь) £ С1 (М+). Отметим, что при п = 1 аналогичный случай рассматривался в работах [29; 17; 18].
Особую практическую значимость имеет уравнение (2.1) с коэффициентом нелинейной теплопроводности К (и) степенного вида. Поэтому в данной статье ограничимся важным частным случаем Ф_1(г) = (А[г(х,Ь)]+)1/х, т.е. К (и) = их, А £ М+. В этом случае имеет место следующий результат.
Утверждение 1. Пусть выполнены соотношения (2.2), (2.6) и матрица А(Ь) с элементами а^(Ь) £ С1 (М+), вектор-столбец В(Ь) с компонентами Ъг(Ь) £ С1 (М+), скалярная функция С(Ь) £ С1 (М+) удовлетворяют следующей системе ОДУ:
А(Ь) = 2А2(Ь) + АА(Ь)^А(Ь),
В(Ь) = 2А(Ь)В(Ь) + АB(í)trA(í), С (Ь) = |В(Ь)|2 + АС (Ь)^А(Ь).
(3.5)
Тогда функция
и(х, Ь) = А
2
х,А(Ь)х +(х, В(Ь))+ С(Ь)
+,
(3.6)
является точным неотрицательным решением уравнения пористой среды (нестационарной фильтрации)
иг
= V • ( иЫи
(3.7)
Уравнение (3.7) является квазилинейным параболическим уравнением с неявным вырождением [10; 1], т. е. уравнение (3.7) параболично при и(х, Ь) = 0 и вырождается при и(х, Ь) = 0. Последнее означает, что при и(х, Ь) = 0 это уравнение из параболического дифференциального уравнения второго порядка вырождается в нелинейное эволюционное уравнение первого порядка типа Гамильтона - Якоби. С вырождением уравнения (3.7) связаны некоторые особые свойства его решений, например конечность скорости распространения носителей решений [10]. В свою очередь, с конечностью скорости распространения носителей решений уравнения (3.7) связаны многие другие типичные свойства [1; 10]: наличие режимов с обострением (отсутствие глобальных по времени
л
1
решений), эффекты локализации, инерции (конечной или бесконечной временной задержки) начала распространения носителя решения.
Приступим к исследованию системы ОДУ (3.5). Предположим, что при Ь = 0 заданы: вещественная симметричная матрица А(0) € Ып вектор-столбец В(0) € Мп, 1(М) и скаляр C(0) € М, т.е.
А(Ь)
А(0), В(Ь)
В(0), C(Ь)
= С (0),
(3.8)
г=о г=о г=о
где Мп(М) — множество п х п матриц с элементами из М; Мп,к(М) — множество п х к матриц с элементами из М [21]. Известно [6; 2], что всякая вещественная симметричная матрица А(0) может быть приведена к диагональной форме с помощью некоторого ортогонального преобразования. Иначе, существует вещественная ортогональная матрица 5 € Мп(М) такая, что
А(0) = 5.0(0)5', Б(0) = diag [¿1(0),..., йп(0)] ,
где йк (0) € М — собственные значения матрицы А(0); к = 1, 2,...,п; 55' = 5 5 = I; I — единичная матрица. Покажем, что если функции А(Ь), В(Ь), С(Ь) определены при Ь = 0, то решение задачи Коши для системы ОДУ (3.5) с начальными условиями (3.8) сводится к решению задачи Коши для некоторого скалярного нелинейного ОДУ первого порядка. Итак, покажем справедливость следующего результата.
Теорема 1. Пусть заданы вещественная симметричная матрица А(0) € Мп(М), вектор-столбец В(0) € МпД(М), скаляр С(0) € М и вещественная ортогональная матрица 5 € Мп(М). Пусть кроме того, г(Ь) — вещественное решение задачи Коши
п Л й = Ц [1-24(0)^)]"г(0)=0, т = -г{1), АеМ+. (3.9)
к=1
Тогда функции А(Ь), В(Ь), С(Ь) определяемые формулами
А(Ь) = ¿(Ь)5д(Ь)Б(0)5' = ¿(Ь)5д(Ь)Б(0)5' А(0),
В(Ь) = ¿(Ь)5д(Ь)Б(0)5' В(0), С(Ь) = ¿(Ь)С(0) + г(г) (В(0), В(Ь)) =
т
С(0) + г(Ь)(д(Ь)5 В(0), ¿(Ь)5 В(0)
(3.10)
(3.11)
(3.12)
являются решением задачи Коши (3.5), (3.8). Кроме того А(Ь) — вещественная симметричная матрица для всех Ь € domainA(t)7 где
Я(Ь) = diag
[1 — 2й1 (0)г(Ь)]
1
[1 — 2йп(0)г(Ь)]
1
(3.13)
йк (0) € М, йк (0) = 0, к = 1, 2,...,п.
Доказательство. Прежде всего введем обозначение Н(Ь) = ^А(Ь). Тогда задача Коши (3.5), (3.8) запишется в виде
А(Ь) = 2А2 (Ь) + АН(Ь)А(Ь), А(Ь)
= А(0),
г=о
В(Ь) = 2А(Ь)В(Ь) + АЬ(Ь)В(Ь), В(Ь)
В(0),
г=о
= С (0).
(3.14)
(3.15)
(3.16)
г=о
С(Ь) = |В(Ь)|2 + АЬ(Ь)С(Ь), С(Ь) Решение задачи Коши (3.14) будем отыскивать в виде
А(Ь)= вБ(Ь)в', (3.17)
П(Ь)=(^[с11(Ь),...,с1п(Ь)] , (3.18)
где Б(0) = 0, (1к £ С1 (М+) — вещественные собственные значения матрицы А(Ь). Подставляя (3.17) в уравнение (3.14) и умножая последнее слева на в , а справа на в, приходим к задаче Коши
Б(Ь) = 2Б2(Ь) + щг)п(г), Б(Ь)
= Б(0).
(3.19)
г=о
Исходя из матричного уравнения (3.19) с помощью несложных выкладок, нетрудно проверить справедливость следующих соотношений
П П г „
П Лк (Ь) = П ао ехр (Ап + 2) Н(г) (т к=1 к=1
(
к=1
г ?
(1к{ь)=(^ | — ^ У^ У ^т^ ^т ^) ^ У ^т^ ^т
к о о о
(И
-1 г
1 exp
(3.20)
(3.21) , (3.22)
где Н(Ь) = ^ (к; (к(0) =0; к = 1, 2,...,п. Далее, для упрощения записи к=1
формулы (3.22), введем в рассмотрение функцию
г ? х(Ь) = J exp А J Н(т) (т
¿(0) = 0.
(3.23)
Из формулы (3.23) вытекает равенство
г
¿(¿) = ехр
А У Н(т) йт
, ¿(0) = 1. (3.24)
Выражение (3.22) с учетом соотношений (3.23), (3.24) упрощается и принимает вид
0Ж*) , .
= 1-24(0)^)' (3'25)
причем
п
1 п
п [¿(*)] ^ йк(0)
п ад
п
к=1 П [1 - 2йк(0>(£)] к=1
Помимо этого, из формул (3.20), (3.24) следует, что
пп . . Ап+2 . .
П = [¿(*)] Л
Лп+2 Л
к=1 к=1
Очевидно, что два последних соотношения приводят к справедливости задачи Коши (3.9) и квадратуры
/Л
~[[[1-2<1к(0)г(г)]*<1г = г + Т, (3.26)
к=1
где Т £ М определяется из начального условия ¿(0) = 0. Итак, в силу (3.25) заключаем, что вещественные собственные значения йк(¿) матрицы А(Ь) запишутся в виде
где ¿(¿) — вещественное решение задачи Коши (3.9). Таким образом, с одной стороны, формула (3.24) приводит к зависимости
"<<>4§- (з-28>
О другой стороны, из соотношений (3.9), (3.27) получим, что
Итак, согласно (3.28), (3.29) заключаем, что функция ¿(Ь) помимо (3.9) удовлетворяет задаче Коши
■¿(г) = А
Е-
(к (0)
¿2(Ь), г(0) = 0, ¿(0) = 1.
1 - 2(к(0)г(Ь)
Очевидно также (см. (3.25)), что имеет место формула
Б(Ь) = ¿(№(1)0(0),
(3.30)
(3.31)
где Q(t) — матрица вида (3.13). Далее, так как А(0) = вБ(0)в', то Б(0) = в А(0)в и легко устанавливается, что соотношения (3.17), (3.31) приводят к справедливости цепочки равенств (3.10).
Теперь перейдем к построению решения задачи Коши (3.15) при условии, что матрица А(Ь) определяется согласно (3.10). Итак, подставляя (3.10) в уравнение (3.15), умножая последнее слева на в и вводя в рассмотрение вектор- столбец Х(Ь) = в В(Ь), приходим к задаче Коши
Х(Ь) = [2Б(Ь) + АЬ(Ь)Т] Х(Ь), Х(0) = в В(0),
(3.32)
причем
В(Ь) = вХ(Ь), Х(Ь) = со1оп(Х1(Ь),.. .,Хп(Ь)). (3.33)
В силу того что диагональная матрица Б(Ь) и скалярная функция Н(Ь) нами определены через решение задачи Коши (3.9) (см. (3.31), (3.28), (3.29)), то из (3.32) имеем
Хк(Ь) = [2йк(Ь) + ЩЬ)]Хк(Ь), Хк(0) = 7к,
где 7к £ М; к = 1, 2,...,п. В итоге приходим к зависимости
(3.34)
Нетрудно убедиться, что соотношение (3.34) с учетом формул (3.24) (3.25) упрощается
Г * 1 Г * 1
Хк (Ь) = Чк ехр 2 [ (к (т) (т ехр А [ Н(т) (т
о о
1к
1 - 2(к(0^(Ь)
¿(Ь).
(3.35)
Итак, используя формулы (3.13), (3.31), а также зависимости 7 = Х(0), Х(0) = вВ(0) из (3.35) получим, что
Х(Ь) = Б-1 (0)Б(Ь)Х(0) = ■ (Ь^(Ь)Х(0) = ¿(Ь^(Ь)в' В(0).
Тем самым выражение (3.33) приводит к справедливости цепочки равенств
В(Ь) = в0-1(0)П(Ь)Х(0) = в0-1(0)П(Ь)в В(0) = ¿(Ь^(Ь)В(0).
Наконец, принимая во внимание (3.28) и учитывая, что вектор-столбец В(Ь) нами определен (см. (3.11)), найдем решение задачи Коши (3.16). Несложно проверить, что решение задачи Коши (3.16) имеет вид
C (t) =
C (0) +
|В(т)|'
¿(г)
dr
z(t),
(3.36)
где |Б(т)|2 = B(0),SD-2(0)D2(т)S B(0) ; Z(j) = 0. Отсюда получаем
о
|В(т)|'
z(T)
г t
dT = B(0), SD-2 (0)
о
D2ir) z{t)
dT
S B(0)
Далее нетрудно убедиться, что для любого к = 1, 2,... ,п выполняется соотношение
г
J Z(T ) о
Тем самым, имеет место представление
z(t)"
|В(т)|' ¿(т)
dT =
B(0),BWi|
Окончательно, из (3.36) следует справедливость формулы (3.12). Легко видеть, что найденные функции (3.10)—(3.12) удовлетворяют начальным условиям (3.8), а следовательно и задаче Коши (3.14)—(3.16).
Итак, для завершения доказательства осталось показать симметричность матрицы A(t) (см. (3.10)) для всех t из области её определения, т. е. для всех t Е domain A(t). С этой целью введем в рассмотрение матрицу G(t) = SQ(t)S , где Q(t) — диагональная матрица, определяемая посредством (3.13). Ясно, что G(t) — невырожденная симметричная матрица. Тогда A(t) = Z(t)G(t)A(0), где A(0) — вещественная симметричная матрица. В первую очередь, убедимся, что матрицы A(0) и G(t) коммутируют. Поскольку выполняется цепочка равенств
Г 1 — 1
A(0)G-1 (t) = A(0) SQ(t)S' = A(0)SQ-1 (t)S' =
A(0)S [I - 2z(t)D(0)] S' = A(0) Iss' - 2z(t)SD(0)S'
A(0) [I - 2z(t)A(0)] = A(0) - 2z(t)A2(0) = [I - 2z(t)A(0)] A(0) = G-1A(0),
t
t
t
то справедливо соотношение G(t)A(0) = A(0)G(t), т. е. матрицы A(0), G(t) коммутируют. Кроме того, имеем
[G(t)A(0)]' = Al(0)Gl (t) = A(0)G(t).
Отсюда следует, что G(t)A(0) — симметричная матрица. Поэтому таковой является и матрица A(t). Наконец заметим, что формулами (3.10)-(3.12) определяются все решения задачи Коши (3.14)—(3.16) с вещественной симметричной матрицей А(0) G Mn(R). Теорема доказана. □
Суммируя результаты этого раздела заключаем, что имеет место
Теорема 2. Пусть выполнены условия (2.2)-(2.6). Пусть вещественная симметричная матрица A(t), вектор-столбец B(t) и скалярная функция C(t) определяются соответственно формулами (3.10)-(3.12). Тогда уравнение пористой среды (нестационарной фильтрации) (3.7) обладает точным неотрицательным решением (3.6).
Пример 1. Пусть вещественная симметричная A(t) является диагональной, т. е. A(t) = diag [ai(t),..., an(t)]. В этом случае система ОДУ (3.5) принимает вид
ài(t) =
2ai(t) + \Y,àk (t) k=1
ài(t), bi(t) =
2ài (t) + X^àk (t) k=1
bi(t),
c(t) = Xc(t)Y àk(t) + Y b2k(t),
(3.37)
k=1
k=1
причем ài(t), bi(t), c(t) G C1 (R+) для i = 1, 2,...,n, где X G R+. Функ-
ции ài (t) будем отыскивать в следующем виде ài (t) =
Si(t)
ГЬ (t)
U=1
Тогда первые n уравнений системы ОДУ (3.37) запишутся
si(t) = -2
ГЬ (t)
j=1
(3.38)
Следовательно, справедлива цепочка равенств вг^) = = ... =
зп(Ь) и для ] = 1, 2,...,п — 1 имеют место соотношения
в](1) = сз + вп(1), ] е М. (3.39)
В итоге, система ОДУ (3.38) с учетом (3.39) сводится к квадратурам
n— 1
П (sn + cj) j=1
dsn = T - 2t, T G R.
(3.40)
1
A
Если все Cj = 0, то из (3.39), (3.40) имеем Sj(t) =
T - (Xn+2)t
2
An+2
+
X G
М+. Таким образом, в этом частном, уравнение (3.7) обладает решением
|х + Ь|2
г An An+2
+c T - (Xn + 2)t
- +
2 [Г - (An + 2)i]
где b G Rra — произвольный постоянный вектор; c G R
+
Теперь рассмотрим случай, когда Cj = 0 и квадратура (3.40) при фиксированном n G N интегрируется в явном виде. Пусть X = 2, n = 2. Тогда из формул (3.39), (3.40) получим
si(t) = r(t) + j, s2(t) = r(t)-j, ci G R \ {0},
где
1 1 r2 1
r(t) = -[e(t)]î+^[e(t)]-î,
6(t) = 12(T - 2t) - ci + 2
36(T - 2t)2 - 6(T - 2t)ci
+
Отсюда следует, что уравнение нелинейной теплопроводности щ = V- (и2Уи) , и = и(х, г) : П х М+ ^ R+, х с R2, имеет точное неотрицательное решение вида
u (xi,x2,t) =
-i--1__i--1. (r(t) + f ) (r2(t) - f) (r(t) - f ) (r2(t) - f)
2c2 exp bi(t)+2c3 exp b2(t) + 2Ç(t)
где bi(t) = 6p(t) - ci^(t), b2(t) = 6p(t) + ciф(t),
<p(t) =
r(t) dt
2
Ф(^ =
+
dt
(r4t) - c4)
c2 y v(t) dt + cl j w(t) dt + c4
(r4t)-î)'
exp[4p(t)],
v(t) = exp[8<p(t) - 2аф(Щ , w(t) = exp[8^(t) + 2аф(Щ , cl G R, l = 1,2,3,4.
2
2
Отметим, что при определенных предположениях полученные решения могут уходить на бесконечность за конечное время, т. е. имеет место так называемый режим с обострением [5; 19].
Список литературы
1. Антонцев С. Н. Локализация решений вырождающихся уравнений механики сплошной среды / С. Н. Антонцев. - Новосибирск : Ин-т гидродинамики СО АН СССР, 1986. - 108 с.
2. Беллман Р. Введение в теорию матриц / Р. Беллман. - М. : Наука, 1976.
3. Вайнберг М. М. Вариационные методы исследования нелинейных операторов / М. М. Вайнберг. - М. : Гостехтеориздат, 1956.
4. Владимиров В. С. Уравнения математической физики / В. С. Владимиров. -М. : Наука, 1971.
5. Квазилинейное уравнение теплопроводности: обострение, локализация, симметрия, точные решения, асимптотика, структуры / В. А. Галактионов, В. А. Дородницын, Г. Г. Еленин, С. П. Курдюмов, А. А. Самарский // Современные проблемы математики. Новейшие достижения. Итоги науки и техники. - М. : ВИНИТИ АН СССР, 1987. - Т. 28. - С. 95-205.
6. Гантмахер Ф. Р. Теория матриц / Ф. Р. Гантмахер. - М. : Наука, 1966.
7. Ибрагимов Н. Х. Группы преобразований в математической физике / Н. Х. Ибрагимов. - М. : Наука, 1983.
8. Калашников А. С. О возникновении особенностей у решений уравнения нестационарной фильтрации / А. С. Калашников // Журн. вычисл. математики и мат. физики. - 1967. - Т. 7, № 2. - C. 440-443.
9. Калашников А. С. Об уравнениях типа нестационарной фильтрации с бесконечной скоростью распространения возмущений / А. С. Калашников // Вестн. МГУ. Сер. мат. механики. - 1972. - № 6. - С. 45-49.
10. Калашников А. С. Некоторые вопросы качественной теории нелинейных вырождающихся параболических уравнений второго порядка / А. С. Калашников // Успехи мат. наук. - 1987. - Т. 42, № 2. - С. 135-176.
11. Капцов О. В. Линейные определяющие уравнения для дифференциальных связей / О. В. Капцов // Мат. сб. - 1998. - Т. 189, № 12. - С. 103-118.
12. Капцов О. В. Методы интегрирования уравнений с частными производными / О. В. Капцов. - М. : Физматлит, 2009. - 184 с.
13. Мартинсон Л. К. Исследование математической модели переноса нелинейной теплопроводности в средах с объемным поглощением / Л. К. Мартинсон // Математическое моделирование. Процессы в нелинейных средах. - М. : Наука, 1986. - С. 279-309.
14. Овсянников Л. В. Групповой анализ дифференциальных уравнений / Л. В. Овсянников. - М. : Наука, 1978.
15. Олейник О. А. Задача Коши и краевые задачи для уравнений типа нестационарной фильтрации / О. А. Олейник, А. С. Калашников, Чжоу-Юй-Линь // Изв. АН СССР. Сер. мат. - 1958. - Т. 22, № 5. - С. 667-704.
16. Рождественский Б. Л. Системы квазилинейных уравнений / Б. Л. Рождественский, Н. Н. Яненко. - М. : Наука, 1978.
17. Рудых Г. А. Об одном подходе построения частных точных решений квазилинейного уравнения теплопроводности с N -пространственными переменными / Г. А. Рудых, Э. И. Семенов // Препринт № 6 ИрВЦ СО АН СССР. Иркутск. — 1991. — 21 с.
18. Рудых Г. А. Построение точных решений многомерного квазилинейного уравнения теплопроводности / Г. А. Рудых, Э. И. Семенов // Журн. вычисл. математики и мат. физики. - 1993. - Т. 33, № 8. - С. 1228-1239.
19. Режимы с обострением в задачах для квазилинейных параболических уравнений / А. А. Самарский, В. А. Галактионов, С. П. Курдюмов, А. П. Михайлов. -М. : Наука, 1987. - 480 с.
20. Сидоров А. Ф. Метод дифференциальных связей и его приложения в газовой динамике / А. Ф. Сидоров, В. П. Шапеев, Н. Н. Яненко. — Новосибирск: Наука. — 1984.
21. Хорн Р. Матричный анализ /Р. Хорн, Ч. Джонсон. - М. : Мир, 1989.
22. Шапеев В. П. Метод дифференциальных связей и его приложение к уравнениям механики сплошной среды : дис. . . . д-ра физ.-мат. наук / В. П. Шапеев. -Новосибирск, 1987.
23. Яненко Н. Н. Теория совместности и методы интегрирования систем нелинейных уравнений в частных производных / Н. Н. Яненко // Тр. IV Всесоюз. мат. съезда. Т. 2. - Л. : Наука, 1964. - С. 613-621.
24. Aronson D. G. The porous medium equation / D.G. Aronson // Some problems in nonlinear diffusion: Lecture Notes in Math. - Springer Verlag, 1986. - N 1224.
25. Aronson D. G. Regularity of flows in porous medium: a survey / D. G. Aronson // Nonlinear Diffusion Equations and Their Equilibrium States. - N. Y. : Springer, 1988. - Vol. 1, N 1. - P. 35-49.
26. Berger M. S. Perspectives in nonlinearity / M. S. Berger. - N.-Y. ; Amsterdam, 1968.
27. Kaplan W. Some methods for analysis of the flow in phase space. Proc. of the symposium on nonlinear circuit analysis / W. Kaplan. - N. Y., 1953. - P. 99-106.
28. Meirmanov A. M. Evolution Equations and Lagrangian Coordinates / A. M. Meirmanov, V. V. Pukhnachov, S. I. Shmarev. - Walter de Gruyter. Berlin, N. Y., 1997.
29. Rudykh G. A. Application of Liouville's equation to construction of special exact solutions for the quasilinear heat equation / G. A. Rudykh, E. I. Semenov // IMACS Ann. Comput. and Appl. Math. - 1990. - Vol. 8. - P. 193-196.
30. Steeb W. H. Generalized Liouville equation, entropy and dynamic systems containing limit cycles / W. H. Steeb // Physica A. - 1979. - Vol. 95, N 1. -P. 181-190.
31. Vazquez J.L. The Porous Medium Equation: Mathematical Theory. Oxford Mathematical Monographs. Clarendon Press, Oxford, 2007.
Рудых Геннадий Алексеевич, доктор физико-математических наук, профессор, Институт математики, экономики и информатики, Иркутский государственный университет, 664003, г. Иркутск, ул. К. Маркса, 1, тел.: (3952)242214, (e-mail: [email protected])
Семенов Эдуард Иванович, кандидат физико-математических наук, старший научный сотрудник, Институт динамики систем и теории управления СО РАН, 664033, Иркутск, ул. Лермонтова, 134, тел.: (3952) 453099, (e-mail: [email protected])
G. A. Rudykh, E. I. Semenov
Research of Compatibility of the Redefined System for the Multidimensional Nonlinear Heat Equation (Special Case)
Abstract. In paper the multidimensional equation of nonlinear heat conductivity is investigated. This equation is presented in the form of the overdetermined system of the differential equations with partial derivatives (the number of the equations are more than number of required functions). It is known that the overdetermined system of the differential equations can be not compatible, at it can not exist any solution. Therefore, for establishment of the fact of existence of solutions and degree of their arbitrariness the analysis of this overdetermined system of the differential equations is carried out. As a result of the conducted research not only sufficient, but also necessary and sufficient conditions of compatibility of the overdetermined system of the differential equations with partial derivatives are received. On the basis of these results with use of the equation of Liouville and the theorem of a necessary and sufficient condition of potentiality of the vector field the approach allowing to construct in some cases exact non-negative solutions of the multidimensional equation of nonlinear heat conductivity with a final velocity of propagation of perturbations is stated. Among the constructed exact decisions are available also such which are not invariant from the point of view of groups of pointed transformations and Lie-Backlund's groups. The special attention is paid to the equation with degree-like coefficient of nonlinear heat conductivity. This equation is the quasilinear parabolic equation with implicit degeneration. This equation from the parabolic differential equation of the second order degenerates in the nonlinear evolutionary equation of the first order like Hamilton-Jacobi.
Keywords: multidimensional nonlinear heat equation, finite velocity of propagation of perturbation, exact nonnegative solutions.
References
1. Antontsev S.N. Localization of solutions of degenerate equations of continuum mechanics (in Russian). Novosibirsk, Institut gidrodinamiki SO AN SSSR, 1986. 108 p.
2. Bellman R. Introduction to the theory of matrices (in Russian). Moscow, Nauka, 1976.
3. Vajnberg M.M. Variational methods for the study of nonlinear operators (in Russian). Moscow, Gostekhizdat, 1956.
4. Vladimirov V.S. Equations of mathematical physics (in Russian). Moscow, Nauka, 1971.
5. Galaktionov V.A., Dorodnicyn V.A., Elenin G.G., Kurdyumov S.P., Samarskij A.A. The quasilinear heat equation: peaking, localization, symmetry, exact solutions, asymptotic behavior, structures (in Russian). Sovrem. probl. matem. Novejshie dostizheniya. Itogi nauki i tekhniki, Moscow, VINITI AN SSSR, 1987, vol. 28, pp. 95-205.
6. Gantmaher F.R. The theory of matrices (in Russian). Moscow, Nauka, 1966.
7. Ibragimov N.H. Groups of transformations in mathematical physics (in Russian). Moscow, Nauka, 1983.
8. Kalashnikov A.S. On the occurrence of singularities in the solutions of the equation unsteady filtration (in Russian). Zhurn. vychis. matem. i matem. fiziki, 1967, vol. 7, no 2, pp. 440-443.
9. Kalashnikov A.S. On the equations of unsteady filtration type with infinite perturbation propagation velocity (in Russian). Vestn. MGU. Ser. mat. mekh., 1972, no 6, pp. 45-49.
10. Kalashnikov A.S. Some questions in the qualitative theory of nonlinear degenerate parabolic equations of second order (in Russian). UMN, 1987, vol. 42, no 2, pp. 135176.
11. Kaptsov O.V. Linear determining equations for differential constraints (in Russian). Matem. sbornik, 1998, vol. 189, no 12, pp. 103-118.
12. Kaptsov O.V. Methods of integration of partial differential equations (in Russian). Moscow, Fizmatlit, 2009. 184 p.
13. Martinson L.K. A study of mathematical model of the non-linear transfer thermal conductivity with volume absorption media (in Russian). Matematicheskoe modelirovanie. Processy v nelinejnyh sredah, Moscow, Nauka, 1986, pp. 279-309.
14. Ovsyannikov L.V. Group analysis of differential equations (in Russian). Moscow, Nauka, 1978.
15. Olejnik O.A., Kalashnikov A.S., Chzhou-Yuj-Lin' Cauchy problem and boundary value problems for equations of unsteady filtration (in Russian). Izv. AN SSSR. Ser. mat, 1958, vol. 22, no 5, pp. 667-704.
16. Rozhdestvenskij B.L., YAnenko N.N. Systems of quasilinear equations (in Russian). Moscow, Nauka, 1978.
17. Rudykh G.A., Semenov E.I. An approach of constructing exact solutions of partial quasilinear heat equation with N — spatial variables (in Russian). Preprint № 6 IrVC SO AN SSSR. Irkutsk, 1991. 21 p.
18. Rudykh G.A., Semenov E.I. Construction of exact solutions of the multidimensional quasilinear heat equation. Computational Mathematics and Mathematical Physics, 1993, vol. 33, no 8, pp. 1087-1097.
19. Samarskij A.A., Galaktionov V.A., Kurdyumov S.P., Mihajlov A.P. Modes with peaking in problems for quasi-linear parabolic equations (in Russian). Moscow, Nauka, 1987. 480 p.
20. Sidorov A.F., Shapeev V.P., Yanenko N.N. The method of differential constraints and its applications in gas dynamics (in Russian). Novosibirsk: Nauka, 1984.
21. Horn R., Dzhonson CH. Matrix analysis (in Russian). Moscow, Mir, 1989.
22. Shapeev V.P. The method of differential constraints and its application to the equations Continuum Mechanics (in Russian). Diss. dokt. fiz.-mat. nauk. Novosibirsk, 1987.
23. Yanenko N.N. The theory of compatibility and integration methods of nonlinear systems PDEs (in Russian). Trudy IV Vsesoyuznogo matematicheskogo s"ezda. Leningrad, Nauka, 1964, vol. 2, pp. 613-621.
24. Aronson D.G. The porous medium equation. Some problems in nonlinear diffusion. Lecture Notes in Math. Springer Verlag, 1986, no 1224.
25. Aronson D.G. Regularity of flows in porous medium: a survey Nonlinear Diffusion Equations and Their Equilibrium States. New York, Springer, 1988, vol. 1, no 1, pp. 35-49.
26. Berger M.S. Perspectives in nonlinearity. New York, Amsterdam, 1968.
27. Kaplan W. Some methods for analysis of the flow in phase space. Proc. of the symposium on nonlinear circuit analysis, New York, 1953, pp. 99-106.
28. Meirmanov A.M., Pukhnachov V.V., Shmarev S.I. Evolution Equations and Lagrangian Coordinates. Walter de Gruyter, Berlin, New York, 1997.
29. Rudykh G.A., Semenov E.I. Application of Liouville's equation to construction of special exact solutions for the quasilinear heat equation. IMACS Ann. Comput. and Appl. Math, 1990, vol. 8, pp. 193-196.
30. Steeb W.H. Generalized Liouville equation, entropy and dynamic systems containing limit cycles. Physica A, 1979, vol. 95, no 1, pp. 181-190.
31. Vazquez J.L. The Porous Medium Equation: Mathematical Theory. Oxford Mathematical Monographs. Clarendon Press, Oxford, 2007.
Rudykh Gennadii Alekseevich, Professor, Institute of Mathematics, Economics and Information Science, Irkutsk State University, 1, K. Marx, Irkutsk, 664003, tel.: (3952)242214, (e-mail: [email protected])
Semenov Edward Ivanovich, Senior Scientist; Matrosov Institute for System Dynamics and Control Theory SB RAS, Post Box 292, 134, Lermontov st., Irkutsk, 664033, tel.: (3952) 453099 (e-mail: [email protected])