Построение фундаментального решения для одного уравнения нечетного порядка Текст научной статьи по специальности «Математика»

УДК 517.951.5 Вестник СПбГУ. Математика. Механика. Астрономия. 2018. Т. 5 (63). Вып. 2 MSC 35A08

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

Б. Ю. Иргашев

Намангангский инженерно-строительный институт, Узбекистан, 160103, Наманган, пр. Дустлик, 12

Для цитирования: Иргашев Б. Ю. Построение фундаментального решения для одного уравнения нечетного порядка // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2018. Т. 5(63). Вып. 2. С. 244-255. https://doi.org/10.21638/11701/spbu01.2018.205

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

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

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

Dpu (x, y) + aDyu (x, y) = F (x, y, D1yu, D1yu,..., Dp-1u, Dy-1u) , (1)

где Dfu = ^-pr, P,q€ N,p > q, a = const.

Это уравнение имеет только лишь одну систему характеристик — y = const, которые являются р-кратными характеристиками. Уравнение (1) при р = 3, q = 1 включает в себя известное уравнение Кортевега—де Вриза, а при р = 3, q = 2 описывает стационарное течение трансзвукового газа [1].

Метод построения фундаментального решения для уравнения (1) дан в [2]. Используя эту методику, для уравнений

D2nu - (-1)n Diu = 0, D 2xn+1u + (-1)n D\u = 0, D Xn+1u + (-1)n D1yu = 0,

© Санкт-Петербургский государственный университет, 2018 244 https://doi.org/10.21638/11701/spbu01.2018.205

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

03хи(х,у) + 02уи+-01уи = /(х,у) У У у

для случая а = 0, используя метод подобия и автомодельного решения, в работах [6, 7] построены и изучены некоторые свойства фундаментальных решений, выраженные через специальные функции. Случай произвольного а исследован в [8].

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

В данной работе для области П = {(х, у) : 0 < х < р, 0 < у < д} построено фундаментальное решение и (х, у, С, п) уравнения

Ь [и] = (-1)" и (х, у) - (-1)т (х, у) = 0, (2)

где (2п + 1, 2т) = 1 (т. е. взаимно просты), п, т € N п > т.

Для фундаментального решения и (х, у, С, п) уравнения (2) доказано следующее

р 9

свойство: функция ^ (х, у) = / / и (х, у, С, п)д (С, п) ¿С^п удовлетворяет в области П

0 0 _ уравнению £[«] = <? (ж, у), где д (х, у) непрерывно дифференцируема в О,.

Основная часть. В работе [9] для уравнения (2) построены частные решения ия (х, у, С, п), в = {0,1,..., т — 1}. Напомним вкратце это построение. Сначала методом подобия строятся дельтообразные решения и* (х, у, С, п) уравнения (2), которые имеют вид

/ (*)

и.*(х У, С, п) =

7Г — Г}\ 2п+1

где

, * /* / 2п+1 \ / 2п+1 \

/;(*)= / ехр -Л— /ЗПсов + ил, -

,! V /V / |у_^|2п + 1

х-С

4в + 1

/З3= /3{- г/3|, = /3| = соэ^, = —-тг,

4т

и доказывается их основное свойство.

Теорема 1. Для € С [а; 6] справедливо 'равенство

о

Иш / и**(х,у,С,п)^(С)^С = ^(хо).

Далее строятся частные решения ия (х, у, С, п), такие что

(х У, С, п) = |у — п|0 Л И, * = (х — С) |у — п|а,

Вестник СПбГУ. Математика. Механика. Астрономия. 2018. Т. 5 (63). Вып. 2 245

где fs (£), Ь выбираются так, чтобы выполнялось равенство

(х, у, е, п) = sgn (у - п) и; (х, у, е, п).

Справедлива теорема [9].

Теорема 2. Имеют место следующие формулы:

1 2т—2 _ / _ \

и = а (2т -2)! £ ^ К- " С"1)' / г"^"1/; (г) ¿г 1 , < > О,

1 2т—2 / Г \

\ /

а = —2т/(2п + 1), Ь = —2па — 1, с+,,с—, —определенным образом вычисленные постоянные, и для них справедливо 'равенство

Приступим теперь к построению фундаментального решения. Для дальнейших выкладок нам понадобятся оценки функции /; (£) на бесконечности. Справедлива следующая теорема.

Теорема 3. Имеет место соотношение

причем:

1) если [22Г^1] четное, то

ит (/; (*))<*> = (-1)* с^ г + 1 + к),к = о,1,2...,

у 2т у

¿¿т, (/; (*))<*> = с^Г /■•)./■• п. 1.2..,

^ если Р^Л ] нечетное, то

(/; (*))<*> = (-1)к +1 + а) , к = о, 1,2...,

¿щ^ (/; (*))<*> = с^г +1 + л) , л = о, 1, 2...,

где

Г2п+11

[ 2ттг ]

(-1)^ сое да-[*§£])-г.), *->-«>,

Г2и±11 + 1

(-1)^- -п (I да _ да» _оо.

сь =

Доказательство. Имеем равенство (для упрощения записи индекс s будем пропускать)

f * (t) = J exp(-Ac,0i) (cos(At)cos(e2Ac) - sin (At) sin (e2Ac)) dA =/i - /2,

0

где

2n + 1

с = -,

2m

I1 (t) = J exp (-Ace1) cos (At) cos (^2Ac) dA =

1 ^ U-- I = ..... ^

e л cos Ac/32dsin At = - / Ac_ie-A л sin (Ac/32 + zs) sin At dA. tt

00

Пусть теперь [с] (целая часть числа) четно, * > 0 (остальные случаи доказываются

аналогично). Тогда можем записать

+^

Ь (¿) = (-1)^ I втА* (А^е-*"* 8ш(Ас/32 + г«))"0" dA. о

Сделаем замену А* = г, отсюда получим

М (с)[с] + 1 [ (гу-\-с\-1 ^ а \ ■ . \ ■ ¿г

D-^T- / Ы exp (-у ft) sin /?2 + ¿sj sinz —.

0

Имеем равенство

Ит (í) = (-l)^ (c)[c]+1sinzs J zc-[c]-lsmzdz =

0

/2n + 1 Л . . / тг / 2n + 1

= ("!) 2 Г + 1 sin sin'

2m s 2 2m

Аналогичным методом получим асимптотику для /2 (t):

2n+ 1 2m

/2(i) = (-#^ j cosXtfr-ie-*" freos (Xе p2+zaf) ° dA,

0

сделаем замену At = z, тогда

'.M-f-D*^ / (f)^"'«p(-(f) A)«.((f) A + =

0

/ \

M (с)_и+1 f rhl f (Z\° ffz^o

= ("!)" J ^^coszexp^-J ftjeos^-j fo + zsJ dz,

0

отсюда

Ит (¿) = (-1)^ (С)[с] + 1С08 2Я J ХС-[С]~1С0ЪХ<1Х

[с] + 1 Сов

0

(-1) * Г —- + 1 \C0SZSC0S

2т 2 2т

Значит, при четном Р^1] имеем

2п+ 1 2т

= + (з)

Теперь покажем, что соотношения (3) можно формально дифференцировать. Итак, нужно показать, что выполняется соотношение

(/; (*))' ^(а + 2) а—2.

Действительно, имеем

(/; (*))' = 11 (*) —12 (*),

1[{г) = - ! Аехр(-Ас/^)со8Ас/3|8тАЫА =! \ ехр(-Xе/3{) сов Xе 1362(1 сов М = 0 0

J ехр (-Xесов Xесов \tdX- ! ф{Хс ех р (-Xе¡в{) сов Xе¡в62 сов Хг ¿X-

0 0

^ У ев!Ас ехр (—Лсв!) вт Лсв| сов Лг ¿Л

1 I 7 \

=-- Д (г)-с / Лс ехр (—Лс/3®) вт (Xе (И2 + га) сое ¿А I .

t

Отсюда получаем

V с ! ЛС/Э^^^/ЛС/:

У ЛС ехр (—ЛСв!) вт (ЛСв| + ¿^сов ЛЫЛ =

1+[с] +00

~|е]+12 / (Асехр(-Ас/^)8т(Ас/328 + ^))([с]+1)8тАЫА~ -А (4):

0

значит,

Л (*) ~ "-у""71 (*) •

248 Вестник СПбГУ. Математика. Механика. Астрономия. 2018. Т. 5 (63). Вып. 2

Аналогично можем записать отсюда

(/*)' (*) с*8Г(а + 2) а—2.

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

Фундаментальное решение и (х, у, С, п) будем искать в виде

т— 1

и (х, у, С, п) = К ^ и« (х, у, С, п)

в 1

я=0

где К — некоторая постоянная, подлежащая определению. Предварительно вычислим следующую сумму:

т—1 ^^ т—1

у^ = (Л \0 _ 2т—2__,+

2^ з* ^ ^ (0„ _ 2п±1 п , л\ 2^

" ' (2П-М(1+,))5

- __(2т~2У-_Г (^±1 (1 + 4)) х

^ /тг ( (2п + 1) (1 + Я + (2т - ^ - 1) (2т - 4в - 1)' X / СОЭ I I

^ 12 V 2т

я=0 у у

^-Л1 /тг / (2п + 1) (1 + з) + (2т -3-1) (2т - 4в - 1)

Ь С°8 ^ V ^

2 2 2т

я=0 у

Далее рассмотрим два случая. Первый случай: пусть з четно, тогда 2т— 1 — з нечетно и выполняется равенство

т 1 '7Г . ^ , 7Г (2п +!)(!+ з) — (2т — 3 — 1) (4й + 1) „ 2 _ ^ _ ^-2т~

= у^ /тг(2п + 2)(1+Я _ 7Г (2т-;?•-!) _ тг

I ; 8111 I 2 2т т в 2

я=0 у

т—1

2 т — з | -I ^-ч

= ^^ соэ (а — ¿в),

я=0

Гдеа = |12±Ш±а ¿= 1(22^11,

^2 т ' т

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

г_1

+ _ t ^n+m+1+i 2m__(2m - 2)!

E c+s = (-i)n

=o 2n+lj!(2n-2^1(1+^-

X oin (("+2)(l+j)„

V 2m v ,JJJ Sin(^)

4 7 V 2m /

В случае нечетного j аналогичные действия приведут к равенству

m— 1

Ec+s =0.

=0

m— 1

Для Y1 c—s имеем следующий результат:

s=0 j

1

Е

cjs

c—, j = 0, 2,4,..., 2l < m - 1, 0, j = 1, 3, .. ., 2l + 1 < m - 1,

где

^n+m+i+i 2m (2m-2)!

c

— = (-1)n

2n

2m v 7 7 sin(^)

2m

Учитывая вышеизложенное, отметим некоторые свойства функции U (ж, y, п):

1) lim D2s—1U (x,y,^,n)= lim D2s—1U (ж, y, n) = 0, s = 1,...,m - 1;

y^n+0 y y^r—0 y

2) lim D2sU (ж, y, £, n)= lim D2sU (ж, y, £, n), s = 0,..., m — 1;

y^n+0 y y^n — 0 y

3) lim U (x,y,^,n)= lim U (x,y,^,n), s = 0,..., 2n.

Пусть теперь в области П = {(ж, y) : 0 < ж < p, 0 < y < q} нам дано неоднородное уравнение

L [u] = g (ж, y), (4)

где функция g (ж, у) непрерывно дифференцируема в Покажем, что функция

p q

F (ж, y) = J j U (ж, y, £, n)g (£, n) d^dn

00

при некотором подборе константы K удовлетворяет неоднородному уравнению (4).

Предварительно сделаем разрез области П прямыми п — е, П + е (е достаточно мало), в результате получим две области П1 = {(£, п) : 0 <C<P, 0 <n<y — е} и П2 = {(£, п) : 0 < C < p, y + е < n < q}. Далее вычислим частные производные функции F (ж, y) по каждой области и устремим е к нулю. Имеем

х

дР , —— = кт

ду е^+0

у—£ р

9 Р

^ ¡¿п/щх, у, V) 9 (£, п) + I dvJ и (х, У> V) 9 (£, п)

0 0 у+е 0

ур

^ / (и (х, у, С, у — е) д (С, у — е) — и (х, у, с, у + е) д (С, у + е)) ¿С+

0

^у—е р 9 Р

00

У+е 0

Теперь перейдем к пределу. Если учесть свойство 2 функции и (х, у, С, п), то первое слагаемое обнулится. В результате получим

Р 9

дР Г Г д

ду

00

Аналогично,

д2Р

тгт = Ит

ду2 е^+0

у —е р 9 Р

д ди д ди

^./1111.1 ^ у ^ *

00

У+е 0

/Уди (х, у, С, у — е) ди (х,у,С,у + е) Л л ,

0

у—е Р

9 ур

+

Дт I <&, |

00

У+е 0

Переходим к пределу. Если учесть свойство 3, то первый предел обнулится и в итоге получим

Р 9

д2Р /" /"д2и

ду2 ду2

00

Таким образом, будем иметь формулу

Р 9

д2т—1 р г п д2т—1и

0у2т-1 =У _/ дуЪт-^Ы^--

00 т—1

Яд 2и

-т^-дЫ^т).

= ( /¿п /и;(х,у,с,п)д(С,п)¿С — /¿п /и*(х,у,с,п)д(С,п)¿С1 .

«=0 \0 0 У 0 /

У0

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

9 Р

Пусть

y p

1 (ж y) = J dn J u: (ж y, C n) g (C n) dC

00

тогда

y+Ay p y0

y p

+ ¿; / dr]J ^y + Ay' ^ ^ " ^Vi L ^ ^ =

00

Для определенности будем считать, что Ay > 0 (случай Ay < 0 рассматривается аналогично). Учитывая теорему 1, можем записать

lim /1 = g (ж, y), Ay^0

y p

r [du;(x,y + eAy,z,ri) ...

J J --9

00

Покажем, что справедливо равенство

y p

= J = к,,,* ,5,

00

Имеем

p y

J J --9 ( ^ =

0 0

p ( y ^

-J J Us: (ж, y + 0Ay,e,n) g (C, n)l0 -/ Us: (ж, y + 0Ay,e,n) g^ (C,n) dn| dC =

00 p

= - у (0Ay)° / ((ж - C) (0Ay)°) g (C, y) dC+

0

p

+ J (y + 0Ay)° / ((ж - C) (y + 0Ay)°) g (C, 0) dC+

0

p y

+ / dtf и: (ж, y + 0Ay, C, n) g'n (C, n) dn.

00

Так как функция / ограничена и a +1 > 0, то во всех трех интегралах можно переходить к пределу под знаком интеграла, что и означает верность равенства (5).

Итак, получаем

р 9 Р 9

д 2т /* /* /* /* д 2т

¿V™ У У ^ ^ ??) ^Л? = 2тойг (ж, у) + У У

0 0 0 0

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

Р 9 Р 9

д2п+1 I- > I- > д 2п+1

У У и(х>У^)9{£,1П)№1П = ] У

0 0 0 0

отсюда следует формула

Ь [и] = (-1)т+1 2т3 (х,у). Таким образом, фундаментальное решение уравнения (2) имеет вид

( 1)т+1 т—1

^ (ж, у, ^ V) = —5--(х>

т я=0

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

Литература

1. Наполитано Л. Д., Рыжов О. С. Об аналогии между неравновесными и вязкими инертными течениями при околозвуковых скоростях // Журнал вычислительной математики и математической физики. 1971. Т. 11, №5. С. 1229-1261.

2. Block H. Sur les equations lineaires aux derives parielles a carateristiques multiples // Ark. Mat. Astron. Fys. Note 1. 1912. Vol. 7, No. 13. P. 1-34; Note 2. 1912. Vol. 7, No. 21. P. 1-30; Note 3, 1912-1913. Vol. 8, No. 23. P. 1-51.

3. Cattabriga L. Una generalizationi del problema fundamentale di valori al contorno per eguationi paraboliche lincari // Annali di matematica pura ed applicata. 1958. T. 46. P. 215-247.

4. Cattabriga L. Potenziali di linea e di dominio per equazioni non paraboliche in due variabili a caratteristiche multiple // Rendiconti del seminario matimatico della univ. di Padava. 1961. Vol. 31. P. 1-45.

5. Абдиназаров С. Краевые задачи для уравнений с кратными характеристиками: дис. .. .док. физ.-матем. наук. Ташкент: Институт матем. АН РУз, 1992. 239 с.

6. Джураев Т.Д., Апаков Ю.П. Об автомодельном решении одного уравнения третьего порядка с кратными характеристиками // Вестник Самарского государственного технического университета. Сер. Физ.-матем. науки. 2007. Вып. 2(15). C. 18-26. https://doi.org/10.14498/vsgtu525

7. Джураев Т.Д., Апаков Ю.П. К теории уравнения третьего порядка с кратными характеристиками, содержащего вторую производную по времени // Украинский матем. журнал. 2010. Т. 62, №1. С. 40-51.

8. Засорин Ю. В. Точные решения сингулярных уравнений вязких трансзвуковых течений // Доклады АН СССР. 1986. Т. 287, №6. С. 1347-1351.

9. Irgashev B. Y. On Partial Solutions of One Equation with Multiple Characteristics and Some Properties of the Fundamental Solution // Ukrainian Mathematical Journal. 2016. Т. 68, № 6. С. 868—893.

Статья поступила в редакцию 15 июля 2017 г.; рекомендована в печать 21 сентября 2017 г. Контактная информация:

Иргашев Бахром Юсупханович — канд. физ.-мат. наук; [email protected]

Composition of fundamental solution for an odd order equation

B. Yu. Irgashev

Namangan Engineering and Construction Institute, pr. Dustlik, 12, Namangan, 160103, Uzbekistan

For citation: Irgashev B. Yu. Composition of fundamental solution for an odd order equation. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2018, vol. 5(63), issue 2, pp. 244-255. https://doi.org/10.21638/11701/spbu01.2018.205

In recent years an enormous attention has been given to examining non-classical equations in quotient derivatives. On the one hand, this is explained by the lack of studying of these equations in theoretical way. On the other hand, often their enclosure to various problems of mechanics, physics and technics has been found out. To less studied non classical equations of high orders belongs the equation in this form:

Dpx u (x, y) + aDq u (x, y) = F (x,y, Dl u, D^ u,..., D%-1u, Dq-1u) ,

^ dsu

l)tu = ——, p, q e N ,p > q, a = const. dts

This equation has only one system of characteristics y = const which is a p-multiple characteristics. This equation with p = 3, q = 1 involves a famous equation of Cortevega — de Vrize, and when p = 3, q = 2 it describes stationary flow of transonic gas.

In this paper in the area of О = {(x,y): 0 < x < p, 0 <y<q} fundamental solution U (x,y,£,n) for an equation

L [u] = (-1)n D2xn+1u (x, y) - (-1)m D2ymu (x, y) = 0,

(2n + 1, 2m) = 1, n,m £ N, n > m

was built.

Keywords: odd-order equation, multiple characteristics, fundamental solution.

References

1. Napolitano L. G., Ryzhov O.S., "An analogy between non-equilibrium and viscous inert flows at transonic velocities", USSR Computational Mathematics and Mathematical Physics 11(5), 166—207 (1971).

2. Block H., "Sur les equations lineaires aux derives parielles a carateristiques multiples", Ark. Mat. Astron. Fys.; note 1 7(13), 1-34 (1912); note 2 7(21), 1-30 (1912); note 3 8(23), 1-51 (1912-1913).

3. Cattabriga L., "Una generalizationi del problema fundamentale di valori al contorno per eguationi paraboliche lincari", Annali di matematica pura ed applicata 46, 215-247 (1958).

4. Cattabriga L., "Potenziali di linea e di dominio per equazioni non paraboliche in due variabili a caratteristiche multiple", Rendiconti del seminario matimatico della univ. di Padava 31, 1-45 (1961).

5. Abdinazarov S., Boundary-Value Problems for Equations with Multiple Characteristics (Doctoral-Degree Thesis (Physics and Mathematics), Tashkent, 1992) [in Russian].

6. Dzhuraev T. D., Apakov Y. P., "On the self-similar solution of a third-order equation with multiple characteristics", Vestn. Samar. Gos. Tekh. Univ., Ser. Fiz.-Mat. Nauk, issue2(15), 18-26 (2007) [in Russian]. https://doi.org/10.14498/vsgtu525

7. Dzhuraev T. D., Apakov Y. P., "On the theory of the third-order equation with multiple characteristics containing the second time derivative", Ukrainian Mathematical Journal 62(1), 43—55 (2010) [in Russian].

8. Zasorin Y. V., "Exact solutions of singular equations for viscous transsonic flows", Doklady Akademii Nauk SSSR 287(6), 1347-1351 (1986) [in Russian].

9. Irgashev B.Y., "On Partial Solutions of One Equation with Multiple Characteristics and Some Properties of the Fundamental Solution", Ukrainian Mathematical Journal 68(6), 868-893 (2016). https://doi.org/10.1007/s11253-016-1263-9

Author's information:

Bakhrom Yu. Irgashev — [email protected]