Двумерное эллиптическое уравнение с нелинейным источником, содержащим степени первых производных Текст научной статьи по специальности «Математика»

УДК 517.952

ДВУМЕРНОЕ ЭЛЛИПТИЧЕСКОЕ УРАВНЕНИЕ С НЕЛИНЕЙНЫМ ИСТОЧНИКОМ, СОДЕРЖАЩИМ СТЕПЕНИ ПЕРВЫХ ПРОИЗВОДНЫХ

TWO-DIMENSIONAL ELLIPTIC EQUATION WITH NONLINEAR SOURCE CONTAINING POWERS OF FIRST DERIVATIVES

И.В. Рахмелевич I.V. Rakhmelevich

Нижегородский государственный университет им. Н.И. Лобачевского, Россия, 603950, г. Нижний Новгород, пр. Гагарина, 23

Nizhny Novgorod State University, 23 Gagarin Ave, Nizhny Novgorod, 603950, Russia

E-mail: igor-kitpd@yandex.ru

Аннотация

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

Abstract

There is investigated two-dimensional non-autonomous elliptic equation, right side of which contains the arbitrary nonlinearity on unknown function and power-law nonlinearities on its first derivatives. There are received the solutions of this equation with additive and multiplicative separation of variables, and also solutions of travelling wave type, self-similar and spherically symmetric solutions. The conditions of existence of founded solutions are formulated. The dependence of solutions on parameters of equation is investigated.

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

Keywords: partial differential equation, elliptic equation, ordinary differential equation, separation of variables, power-law nonlinearity, solution of travelling wave type, self-similar solution.

Введение

Исследование дифференциальных уравнений в частных производных со степенными нелинейностями является важной частью современной нелинейной математической физики. Этой проблеме уделяется большое внимание как в справочниках и учебных пособиях [1,2], так и в оригинальных работах [3-7]. При этом особый интерес представляют нелинейные неавтономные уравнения, так как они описывают нелинейные явления в нестационарных и неоднородных системах. Настоящая работа посвящена изучению решений эллиптического уравнения с правой частью, явно зависящей от координат и содержащей произвольную нелинейность от искомой функции и произведение степенных функций от ее первых производных. Анализ решений данного уравнения проводится с помощью метода разделения переменных (РП), который широко используется при исследовании нелинейных дифференциальных уравнений в частных производных [1-9]. При этом находятся условия существования полученных решений, которым должна удовлетворять правая часть уравнения.

Постановка задачи. Простейшие решения

Рассмотрим двумерное уравнение второго порядка эллиптического типа относительно

неизвестной функции и(х'у) , с правой частью, включающей произвольную нелинейность по искомой функции и степенные нелинейности по ее первым производным:

д2и д2и

дх2 дУ2

= ё (и )/( х У)

ди дх

Р1

ди

дУ

\Р2

(1.1)

Здесь Р1'- параметры уравнения, ё (и)'/ (х'у) - заданные функции (их вид будет уточняться в дальнейшем).

Теорема 1.1.

Пусть выполнены условия: ё(и) = ёо ехР(Уи) / (х, у) = /1( х)/2(у)

1. Тогда при условии, что

/2(У) = /20(2а2У + Ь2)-Р2 ехР ~У(а2У2 + Ь2У)

уравнение (1.1) имеет следующее решение:

(1.2)

(1.3а)

и (х, у) = X (х) + а2 у + Ь2 у + с2

(1.3б)

причём Х (х) удовлетворяет обыкновенному дифференциальному уравнению (ОДУ):

X"(х) = gJl(х)[Х'(х^1 ехр[уХ(х)\-2а2, где коэффициент gl = g0/20 ехР(Ус2 } . 2. Аналогично, при условии, что

У1(х) = /ю(2а'Х + ¿1)-р1 ехр^уЦх2 + Ь'Х)\

(13в)

(1.3г)

(14а)

(1.4б)

уравнение (1.1) имеет следующее решение: и( х, у) = ахх2 + Ь1х + с1 + У (у)

причём У (у) удовлетворяет следующему ОДУ: у '(у) = g 2/2 ( у) [у '(у)]р 2 ехр[уУ(у)\-2а1, (^

где коэффициент

g2 = g0f10 ехР(УС1). (14г)

Доказательство.

1. Согласно известной схеме аддитивного РП [1,2], решение уравнения (1.1) ищем в виде:

и( х у) = Х (х) + У (у). (15)

Подставляя (1.5) в уравнение (1.1) и учитывая условия (1.2), получаем:

X "(х) + У "(у) = х)Ф2(у), (1.6)

где

Ф1(х) = ^(х)[Х'(х^1 ехр[уХ(х)], Ф2(у) = Ь(у)\?'(у)]32 ехр[уУ(у)]. (16а)

Дифференцируя уравнение (1.6) по х, можно привести его к виду: X""( х)

—гтт = g 0Ф2( у)

Ф1(х) . (1.7)

Так как левая часть (1.7) зависит только от х, а правая часть - только от у, то это уравнение можно удовлетворить, только если обе его части равны постоянной, откуда

следует, что Ф2(у) = а = соп^ . Подставляя это выражение в (1.6), получаем:

У" (у) = Ш х) - X '' (х)

(1.8)

где ^ . Аналогично (1.7), левая и правая части уравнения (1.8) зависят от

разныхпеременных, и поэтому должны быть равны некоторой постоянной. Отсюда, с

учетом (1.6а), следует, что Х (х) удовлетворяет обыкновенному дифференциальному уравнению:

X'(х) - ё/( х)[Х'(х)Р' ехр[уХ (х)] = -2а2,

(1.9)

У( у)

а - переопределенной системе из двух уравнений:

Иу)Г2 ехр[у7(у)] = ~2//2(у), П>) = 2а2.

Из второго уравнения системы (1.10) находим, что У(у) = а у + ¿2 у + С2. Подставляя это выражение в первое уравнение системы (1.10), получаем функцию /г(у) :

—о 2

/2 (У) = /20 (2а2У + Ь2 У ехР — У(а2У + Ь2У)

(111)

У20 а2ехр( Выряжяа отсюда ~

где -720 2 1 2 . Выражая отсюда 2, получаем формулу (1.3г) для коэффициента ё1. Из (1.9) и (1.11) следует, что первая часть теоремы доказана.

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

у(у) и переопределенную систему уравнений для Х(х) :

У "(У) — ё2/2(У)[Г'(У)]Р 2 ехр[уУ (у)]= —2а1

(1.12)

[Х '(х)]р1 ехр[уХ (х)] = /1( х), Х "(х) = 2а1. (113)

Из второго уравнения системы (1.13) находим, что Х (х) = а1х + ¿1х + С1. Подставляя это выражение в первое уравнение системы (1.13), получаем функцию /1(х) :

/1(х) = /ю(2а1х + —р1 ехр[—у(а1х2 + ^х)]

+ ¿^п (1.14)

/10 = ~1 ехр(—Ус1)

где 710 1 ' 1. Выражая отсюда постоянную 1, получаем формулу (1.4г) для коэффициента ё2 . Из (1.12) и (1.14) следует, что вторая часть теоремы доказана. Пример.

Пусть у = а2 = 0, /20 = ¿2Р2, /1(х) - 1. Тогда, из (1.2) и (1.3а) ё(и) - ё0, /(х,У) - 1, т.е. в условиях данного примера уравнение (1.1) автономное и явно не зависит от искомой

функции. Решение, определяемое теоремой 1.1, имеет вид и( х'у) = Х (х) + ¿2 у + С2, где

Х ( х) определяется выражениями:

X (х) =

11 ч2-^

Ы1 -р1)х + С1 )1-р1 + Xо при р1 * 1,р1 * 2

^(2 -Р1)

а1ехр( glx) + Xо

1 g1

--1п| х - хо + X0

при при

в = 1

Р1 = 2

Здесь а1'С1'хо'Xо - произвольные постоянные. Теорема 1.2.

Пусть выполнены условия:

g(и) = gоu т, f (х, у) = /1( Х)f2(y) 1. Тогда, при условии что

f2(у) = f2о(у -уо)1-в1 -у уравнение (1.1) имеет решение вида

и( х у) = Ь2X (х)( у - уо)

(115)

(1.16а)

(1.16б)

где

X (х)

является решением ОДУ:

(1.16в)

X" (х) = glfl( х)[х' (х^1 [X (х)]т+р2 Здесь ¿2' -произвольные постоянные, а Й1 определяется выражением:

Й = g^•f20Ь2

(1.16г)

Здесь и всюду далее используется обозначение в + Р2. 2. При условии что

f2 (у) = f20 (С ^ (Яу) + С2 ск(Ху))-Р2 (С1 сЬ (Яу) + С2 З^Яу))1-^1

(1.17а)

уравнение (1.1) имеет решение вида и(х, у) = X(х)(с1 сЬ (Яу) + с2 8Ь(Яу))

где X (х) является решением ОДУ: X" (х) = x)[X'(х)]р1 [X(х)]т+р2 - Я2X(х),

а коэффициент Й1 определяется выражением:

Й1 = ЙоЛЛ"2

3. При условии что

(1.17б)

(117в)

(117г)

/2(У) = /20(с2cos(V) — ^ЧАу^^Ш^^ + с2 йп(АУ))1—т—р1. (118а) уравнение (1.1) имеет решение вида

и(x, У) = Х(х)(с1cos (ХУ) + с2 , (1.18б)

где Х (х) является решением ОДУ: Х"(х) = ёь/1(х)[Х'(х)]в1 [Х (х)Г2 + X2 Х (х)

(1.18в)

причем коэффициент ё1 определяется выражением (1.17г). Доказательство.

1. Используя мультипликативное РП [1,2], решение уравнения (1.1) ищем в виде: и( x, У) = Х (х)У (У). (119)

Подставляя (1.19) в уравнение (1.1) и учитывая условия (1.15), получаем: Х "(х) У "(у) , , , ,

—— + —— = ё0 VI (х)у 2 (у)

Х(х) У(у) (1.20)

где

х) = /1(х)[Х '(х)]р1 [Х (х)рр2 -1, V 2( у) = /2(у)[У '(у )]р 2 [У (у)Гр1 —1

(1.20а)

Аналогично доказательству теоремы 1.1, в результате дифференцирования (1.20) по х, получаем:

1 й ( Х"(х)^

у1(х) йх ^ Х(х) ) ёоV2(у). (121)

Правая и левая части уравнения (1.21) зависят от разных переменных, откуда следует, что V 2(У) = ~2 = ю^. (1.21а)

Подставляя выражение (1.21а) в (1.20), получаем:

УМ = х)—

У (у) Х (х)

Здесь ё1 ё0а. Разделяя переменные в (1.22) и учитывая (1.20а), находим уравнение для Х(х) , и переопределенную систему из двух уравнений для У (у) :

(122)

- й/1( x)[X'(х)]р1 [X (х)Г+р 2-1 = -а А (х) , (1.23)

= а [У ' (у)]р2 [У (уХГ -1 =-0^-

у (у) , /2(у). (1.24)

При анализе системы (1.24) рассмотрим возможные частные случаи. Случай а): а = о. Тогда из первого уравнения системы (1.24) находим, что

У (у) = Ь2( у - уо). (125)

Подставляя (1.25) во второе уравнение системы (1.24), получаем выражение для /2(у) : /2(у) = (у - уо)1"^1. (126)

При этом уравнение (1.23) упрощается следующим образом: X" (х) = й/{х)[Г(х)Г1 [X(х)]т+р2. (1.27)

Для нахождения коэффициента Й1 сравним соотношения (1.26) и (1.16а), откуда а2 ЛЛ . Далее, учитывая, что Й1 Йо<~2, получаем выражение (1.16г).

Случай б): а = Я > о. Тогда из первого уравнения системы (1.24) получаем, что

У (у) = с сЬ (Яу) + С2 Б^Яу)

(128)

Аналогично предыдущему случаю, находим /г( у), подставив (1.28) во второе уравнение (124):

(1.29)

/2 (у) = ~2Я-Р2 (с «Ь (Яу) + С2 сЬ(Яу))-р2 (с1 сЬ (Яу) + С2 8Ь(Яу))1-т-Р1 . Случай в): а = -Я < о . Тогда из первого уравнения системы (1.24)

У (у) = с1 с°в (Яу) + с2 «1п(Яу) , (1.3о)

а функция /г( у) определяется выражением:

/2(у) = ~2Я-Р2(с2со«(Яу)-с1 8т(Яу))-р2(с1сов(Яу) + с2 8т(Яу))1-у-р1 . (1.31) Из (1.20) получаем уравнение относительно X (х) для случаев б), в) :

X"(х) = x)[X '(х)]р1 [X (х)]7+р2 + Я2 X (х)

(132)

где верхний знак относится к случаю б), а нижний - к случаю в). Выражение (1.17г) для Й1 находится путем рассуждений, аналогичных проведенным при анализе случая а).

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

получены выражения (1.25), (1.28), (1.30) для У(у) , уравнения (1.27), (1.32) для Х(х) , и

выражения (1.26), (1.29), (1.31) для /г( у), при которых найденное решение существует. Теорема доказана.

Замечание. Аналогично теореме 1.2, можно сформулировать утверждение о решениях

-у* /_V 11 % /_V у

уравнения (1.1), которые получаются из (1.16), (1.17), (1.18) путем замены ^ -У> ^ >

/1 ^ /в1 ^ Р2, а также соответствующие условия существования этих решений. Доказательство этого утверждения полностью аналогично доказательству теоремы 1.2.

2. Функциональное разделение переменных

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

Теорема 2.1.

Пусть /(х'у) = ^(2) 2 = с1х + с2у, где ^ - произвольная функция, с1' с2 - некоторые постоянные. Тогда уравнение (1.1) имеет решение типа бегущей волны и( х,у) = и (2), причем функция и (2) является решением следующего ОДУ:

и " ( 2 ) = Сё (и ( 2) )[и' ( 2^ ^ ( 2)

(2.1)

с Р1 с в 2 С = Ч 2

где

с1 + с2

Доказательство.

Используя известный метод функционального РП [1,2], решение уравнения (1.1) ищем в виде:

и( х, у) = и (2) 2 = Х (х) + У (у)

(2.2)

Подставляя (2.2) в (1.1), после дифференцирования и элементарных преобразований получаем:

и"(2) + Х -- (х) + У" (у) = ё (и )[и - 1 / (х, у)[Х'(х)]в1 [У'(у)]32 и"(2) [Х(х)]2 + [У'(у)]2 " ё1 ( )А ( П [Х"(х)]2 + [У(у)]2 (2 3)

Положим Х(х) = СlX, У(у) = с2у . Подставив эти выражения в (2.3) и учитывая, что по

условию теоремы /(у) = ^(2), приходим к уравнению (2.1) относительно и(2) . Теорема доказана.

Следствие 2.1.

Если /(x, у) = 1, то решение уравнения (1.1) типа бегущей волны может быть представлено в неявном виде следующим образом:

1.Если Рб * 1, Рб * 2, то

с1х + с2у - z0 = ^\у(СО(и) + А)] V ёи

2.Если Рь = 2, то

с1х + с2у - z0 = А| ехр(-С0(и)) ёи

3.Если вь =1, то

ёи

(2.4а)

(2.4б)

С1х + С2 у - Zо = |

СО (и) + А

(2.4в)

С,,С2,z0,А V = 2-Ву О(и)

где 1 2 0 - произвольные постоянные, Рь; к ' определяется выражением:

О(и ) = | Й(и )ёи (24г)

Доказательство.

Поскольку по условию р(z) = 1, то уравнение (2.1) можно записать в виде:

и"( z )[и'( z )]1-рь- С*0^ = 0

^ . (2.5)

При условиях вь * 1 Рб * 2 уравнение (2.5) в результате интегрирования сводится к уравнению первого порядка:

[и'(^^(СО (и (+ А). (2.6)

В результате решения уравнения (2.6) получаем:

z - z0 = | \у(С0(и) + А)]- V1 ёи

(2.7)

Подставив в (2.7) выражение С1х + С2у и заменив и ^ и, получаем для данного случая решение (2.4а). Проводя рассуждения, аналогичные представленным выше, для случаев

вь = 2 и вь = 1, получаем решения (2.4б) и (2.4в) соответственно. Доказательство закончено.

Теорема 2.2.

/ (х у) = Ь.х у "р2 р (z) z = х 2 + у2 77

Пусть ' 0 ^ \ ;> где р - произвольная функция. Тогда

уравнение (1.1) имеет решение и( x, у) = и (z), причем функция и (z) является решением ОДУ:

ё (zU'( z)) = Ьй (и ( z ) )[и'( ^ Р ( z ) ^ , (2.8)

Ь = Ьо • 2Рь-* р(z) = 1 Вт = 1

где 0 . В частности, если выполнены условия 4 ' , Рь , то решение

уравнения (1.1) можно представить в неявном виде:

2 2 х + у = с0 ехр

^ йи ^

1

ЪС(и) + а0

(2.9)

с, ап

где 0 0 - произвольные постоянные.

Доказательство.

Аналогично доказательству теоремы 2.1, используем метод функционального РП.

Положим Х(х) = х , У(у) = у . Подставляя эти выражения в (2.3), и учитывая выражение

/ (х у) для , получаем:

+ 1 = Ъ-ё (и ( 2))[и " ( 2)13б—1 Р ( 2) и (2) 22 . (2.10)

— (2и' ( 2 ) )= 2и " ( 2) + и' ( 2)

Умножив (2.10) на 2 и учитывая, что —2 , приходим к уравнению

(2.8).

Предположим теперь, что выполнены условия р (2) =1, Рб = 1. Тогда (2.8) можно представить так:

—(2и"(2) — ЪО(и (2))) = 0 —2 . (2.11)

Понижая порядок уравнения (2.11) и решая полученное уравнение первого порядка, находим:

2 = с0 ехр

' —и ^

I

ЪО(и) + а0

(2.12)

„2 , ,.2

где с0,а0 - произвольные постоянные. Подставляя в (2.12) 2 х + у и заменяя и ^ и. получаем решение в форме (2.9). Теорема доказана.

Теорема 2.3.

Пусть

/(ЛУ) = (с1 + 2ах)—р1 (с2 — 2аУ)—р2 [с1 + 2ах)2 + (с2 — 2аУ)2 ]р(2), (2.13а)

2 = а(х2 — у2)+ с1х + с2У + 20 , (2.13б)

где р - произвольная функция. Тогда уравнение (1.1) имеет решение ^ ,у' ^ причем функция и(2) является решением ОДУ:

и"(2) = ё (и (2))[и' (2 )]Рб р (2). (2.14)

Доказательство.

Аналогично доказательству теорем 2.1, 2.2, используем метод функционального РП.

2 2 т-т X(x) = ах + С>х + х0, У (у) =-ау + С2у + у0 ^

Положим у 7 1 0 ^' у ^ У0. Подставляя эти выражения в (2.3),

и учитывая выражение для / (х'у), получаем уравнение (2.14). Теорема доказана.

Теорема 2.4.

1. Пусть

/(х у) = х_р2у"р1 (х2 + у2 р(zX z = ху, (215)

где р - произвольная функция. Тогда уравнение (1.1) имеет решение и(x, у) = и(z), причем функция и (z) является решением уравнения (2.14).

Доказательство.

Будем искать автомодельное решение уравнения (1.1), которое имеет вид: и( х, у) = хаи (z), z = уха, (2.16)

где а'а - параметры, которые должны быть определены ниже. Подставляя (2.16) в (1.1), получаем:

а+2а

х

22

и"(z) + а(2а + а -1)уха+а-2и'(z) + а(а -1)ха-2и(z) =

1 а у 1 + ^

. х2 ,

= g(xаU(z))f(x, у) {зуха+а-ипх> + а*а-1и(^)|!1 {с"ЧГ'( z)f(217)

Уравнение (2.17) допускает разделение переменных только в том случае, если а = 0 а = 1; тогда z определяется второй формулой (2.15). При указанных значениях параметров (2.17) сводится к виду:

(х 2 + у 2 )и ) = Й (и ^))[и'(z )]Рь / (х, у)ув1 хв2

(2.18)

Из (2.18) с учетом (2.15) следует, что функция и (^ должна удовлетворять уравнению (2.14). Теорема доказана.

Заключение

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

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

Полянин А.Д., Зайцев В.Ф. 2002. Справочник по нелинейным уравнениям математической физики: точные решения. М., Физматлит, 432.

Polyanin A. D. and Zaitsev V. F. 2012. Handbook of Nonlinear Partial Differential Equations, 2nd Edition, Chapman & Hall/CRC Press, Boca Raton.

Полянин А.Д., Зайцев В.Ф., Журов А.И. 2005. Методы решения нелинейных уравнений математической физики и механики. М., Физматлит, 256.

Polyanin, A. D., Zaitsev, V. F., Zhurov A. I. 2005. Metody resheniya nelineynyh uravneniy matematicheskoy fiziki i mehaniki. M., Fizmatlit (In Russian).

Рахмелевич И.В. 2015. О двумерных гиперболических уравнениях со степенной нелинейностью по производным. Вестник Томского государственного университета. Математика и механика, № 1: 12-19.

Rakhmelevich I.V. 2015. O dvumernyh hyperbolicheskih uravneniyah so stepennoy nelineynostiu po proizvodnym. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika, №1: 12-19. (In Russian).

Рахмелевич И.В. 2015. О некоторых новых решениях многомерного уравнения в частных производных первого порядка со степенными нелинейностями // Вестник Томского государственного университета. Математика и механика, № 3: 18-25.

Rakhmelevich I.V. 2015. O nekotorykh novykh resheniyakh mnogomernogo uravneniya v chastnykh proizvodnykh pervogo poryadka so stepennymi nelineynost'yami. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika, No 3: 18-25 (in Russian).

Рахмелевич И.В. 2016. О решениях многомерного дифференциального уравнения произвольного порядка со смешанной старшей частной производной и степенными нелинейностями. Владикавказский математический журнал., 18 (4): 41-49.

Rakhmelevich I.V. 2016. O resheniyakh mnogomernogo differentsialnogo uravneniya so smeshannoy starshey chastnoy proizvodnoy i stepennymi nelineynost'yami. Vladikavkazskiy matematicheskiy zhurnal. 18(4): 41-49 (in Russian).

Рахмелевич И.В. 2016. О редукции многомерных уравнений первого порядка с мультиоднородной функцией от производных. Известия вузов. Математика, № 4: 57-67.

Rakhmelevich I.V. 2016. Reduction of multidimensional first order equations with multi-homogeneous function of derivatives. Russian Mathematics, 60(4): 47-55.

Рахмелевич И.В. 2017. О многомерных уравнениях в частных производных со степенными нелинейностями по первым производным. Уфимский математический журнал, 9 (1): 98-108.

Rakhmelevich I.V. 2016. On multidimensional partial differential equations with power nonlinearities in first derivatives. Ufa Mathematical Journal, 8(4): 98-108.

Miller J. (Jr.), Rubel L.A. 1993. Functional separation of variables for Laplace equations in two dimensions. Journal of Physics A., 26: 1901-1913.

Zhdanov R.Z. 1994. Separation of variables in the non-linear wave equation. Journal of Physics A., 27: L291-L297.