Потенциалы сопряжённых пар гиперболических уравнений и их прикладные аспекты Текст научной статьи по специальности «Математика»

Ю. И. Фёдоров

ПОТЕНЦИАЛЫ СОПРЯЖЕННЫХ ПАР ГИПЕРБОЛИЧЕСКИХ УРАВНЕНИЙ И ИХ ПРИКЛАДНЫЕ АСПЕКТЫ

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

Ключевые слова: потенциалы сопряжённых пар гиперболических уравнений второго порядка.

Yu. Fedorov

POTENTIALS OF CONJUGATE PAIRS OF HYPERBOLIC EQUATIONS AND THEIR APPLICATIONS

The article introduces the concept of a potential conjugate pair for second order hyperbolic equations and an associated with it concept of dual equation. The applications of the potential for integral representations of solutions of hyperbolic equations and the parameter continuation are regarded, and the potential practical value is discussed.

Keywords: potentials conjugate pairs of second order hyperbolic equations.

В односвязной области D плоскости XOYрассматривается линейное гиперболическое дифференциальное уравнение с частными производными второго порядка с двумя переменными

L(u) ° uy + a(X y)ux + b(x y)uy + c(x y)u = 0 (1)

с коэффициентами, удовлетворяющими условию

a(x, y), b(x, y), c(x, y), ax(x, y), by(x, y) e C(D) . (2)

Условия гладкости решения u (x, y): u(x, y) e C1 (D) n C2 (D).

Введём основное понятие рассматриваемого метода, то есть понятие потенциала сопряжённой пары гиперболического уравнения, и рассмотрим построение потенциалов. Пусть L*(u) — дифференциальный оператор, действующий на функцию v(x,y), формально сопряжённый с оператором L(u) [3]:

L (v) = vxy - (av)x - (bv)y + cv.

Если выполняются условия (2), то "u(x, y), v(x, y) e Cl(D) n C2(D) справедливо тождество Грина:

2(vL(u) - uL4 v)) = (vuy - uvy + 2auv)x - (uvx - vux - 2buv)y. (3)

Из этого тождества следует, что если Ь(и) = 0, Ь(у) = 0 в области Б, то

(иух - уих - 2ЬиУ)у = (уиу - иуу + 2аиУ)х, то есть Ру (х, у) = 0х(х, у), (4)

где Р(х, у) = иух - уих - 2Ьиу, ((х, у) = уиу - иуу + 2аиу.

Для выражения Р(х,, у)Сх + ((х, у) Су выполнены условия известного признака полного дифференциала в односвязной области Б [2]. Поэтому выражение

Р(х, у)Сх + ((х, у)Су = (иух - уих - 2ЬиУ)Сх + (уиу - иуу + 2аиу)Су

является полным дифференциалом в области Б, то есть существует функция Z (х, у), определённая в области Б такая, что

<С1 = Р(х, у)Сх + ((х, у) Су ° (иух - уих - 2ЬиУ)Сх + (уиу - иу + 2аиу)Су.

В векторном анализе условие Ру (х, у) = ((х(х, у) означает, что векторное поле

Г(х, у) = (Р(х, у); (((х, у)) — безвихревое в Б [2]. При этом функцию и(х, у), у которой Р(х,у)с<х + ((х,у)Су =-Си(х,у), называют скалярным потенциалом безвихревого векторного поля Г. Из этих соображений далее рассматривается и(х, у), а не Z(x, у):

Си (х, у) = - Р ■ Сх - ((■ Су ° (уих - иух + 2ЬиУ)Сх + (иу - уиу - 2аиу)Су. (5)

Сравнивая форму первого дифференциала функции двух переменных

Си = их ■ Сх + иу ■ Су

и формулу (5), заключаем, что это дифференциальное равенство равносильно системе уравнений

их = уих - иух + 2Ьиу,

(6)

иу = иУу - уиу - 2аиу.

Здесь и(х, у) и У(х, у) — любые решения из О(Б) п С2(Б) уравнений Ь(и) = 0 и

*

Ь (у) = 0 соответственно.

Если М((х0, уо) — фиксированная начальная точка линии интегрирования, М(х, у) — точка области Б, М0М с Б, то и(х, у) восстанавливается по своему полному

дифференциалу с точностью до вещественного постоянного слагаемого и0= и(х0, уо) по формуле

м

и(х, у) = ио + | (уи5 - иу5 + 2Ьиу)ёэ + (иу( - щ - 2аиу)&,

м„

(7)

где 5, t — переменные интегрирования, криволинейный интеграл не зависит от формы дуги М0М [2]. Рассматриваемые дуги кусочно-гладкие.

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

Р(х, у) = (Рх, у); 0(х, у)) = (иУх - уих - 2Ьиу уиу - иУу + 2аиу),

потенциалом которого является и(х, у), формируется с помощью решений и(х, у) исходного и у(х, у) сопряжённого уравнений [5].

Отметим, что в книге [2] в сходной ситуации, но для обыкновенных дифференциальных

^ 2 и с1и ^ 2 ё

операторов Ь(и) = а0(х)—- + а1(х)--+ а2(х)и и Ь*(и) = —-[а0(х)у\--[а1(х)у\ + а2(х)у

ёх -1-- -1 -1--

ёх

ёх2

ёх

в тождестве уЬ(и) - иЬ(у) = Р(х), где Р(х) = а0(уи - иу) + (а1 - а0')иу, функция Р(х) называется конъюнкцией и(ху) и у(х) относительно оператора Ь.

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

Представление (6) потенциала позволяет записать уравнение (1) в дивергентной форме:

(уих - иух + 2Ьиу)у - (иуу - уиу - 2аиу)х = 0.

Теперь определим понятие двойственного уравнения, связанное с понятием потенциала сопряжённой пары. Для этого выведем уравнение, которому удовлетворяет потенциал и(х, у). С этой целью исключим функцию и(х, у) из системы уравнений (6). Это можно сделать, например, с помощью следующих элементарных преобразований системы (6):

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

иу = -уиу + (уу - 2ау)и.

^ = их - ^- 2Ь^ • и у

= -иу +

V

и

Полагая — - 2Ь,— - 2а е С1(П), первое уравнение системы (8) продифференцируем

V V

по у, а второе — по —

и )у = ^ - IV - 2Ь) ' "У - IV - 2Ь)у ' ^

(9)

( Чу.)

V v J

= -"ух +

V

- 2а

V

) ( ■ и„ +

V

- 2а

V

Л

■и.

J -

Затем выразим из системы уравнений (8) производные и- и иу и подставим их выражения в систему (9):

V

( Чу )

V v J

-У

V -2Ь

и (

у V

+

V

- 2а

V

)

= "у- +

V

( V.. )(Ч (V

^ + | - 2Ь | ■и

V V V

V.,

V

- 2а

V

J

и

J J

) (V.

- 2Ь | ■ и,

у

+

V

J

Л

V

- 2а

V

и.

J -

Выполним элементарные преобразования и умножим обе части каждого уравнения на у(х, у):

Ч„.—Ч-I —-2Ь |■П, = V",,,,-| -2Ь

-У

V

Ч —Х-Ч ■

у V у

У -У

( V )

- 2а

V

иV -1 — I ■ иV + 2Ь„ ■ Vи,

V

J

(

V

- 2а

V

Л

Ч- = -™у- + I V - 2Ь

V

- 2а

V

Л

"V+

J

( V )

V v J

■ иу - 2а, ■ Vи.

Сложив почленно уравнения этой системы, получим следующее уравнение:

Ч-у -

-а

Ч- - (7 - Ь) ЧУ = ( ЬУ - а-)

а ■ V и.

(10)

Уравнение, которому удовлетворяет потенциал Ч(-, у), в данном случае (10), будем называть двойственным уравнению (1), а решения двойственного уравнения — двойственными решениями к и(-, у). Таким образом, потенциал Ч(-, у) — двойственное решение к и(-, у).

Если в области В выполняется условие аХ = Ьу, то двойственное уравнение становится линейным гиперболическим уравнением. Таким образом, справедлива теорема 1.

V V

Теорема 1. Если выполняются условия (2), ах = Ьу, — - 2Ь,— - 2ае 0\П), то каковы

V V

бы ни были решения и(х, у) и \(х, у) исходного и сопряжённого уравнений класса С1(0) п С2 (В), соответствующий потенциал Щх, у), определяемый системой уравнений (6) и построенный на сопряжённой паре и(х, у), у(х, у), является решением двойственного линейного гиперболического уравнения

+

Г Vy ^

а —-V Г

■ их+I Ь - V IЩ = 0. (11)

Потенциал Щ(х, у), как двойственное решение, имеет интегральное представление (7).

Вопрос о виде двойственного уравнения, если условие ах = Ьу в области В не выполняется, обсуждается в других работах автора.

Так как коэффициенты двойственного уравнения зависят от у(х, у), то двойственное уравнение — это не одно уравнение, а класс уравнений. Выбирая различные у(х, у), будем получать различные уравнения этого класса и различные двойственные решения, то есть различные потенциалы. Объединяет эти уравнения то, что потенциал Щ(х, у) является их решением и имеет интегральное представление (7). Эта точка зрения позволяет рассматривать некоторые уравнения как двойственные к другим, исходным, более изученным (модельным, базовым) и найти двойственное решение (общее решение или решение начальной, краевой задачи) в виде потенциала и(х, у) с помощью интегрального представления (7), если известны и(х, у) и у(х, у). Таким образом, предлагаемая теория потенциалов предоставляет ещё один метод интегрирования уравнений.

Простейшей реализацией метода интегрирования с помощью потенциалов сопряжённых пар является продолжение решений гиперболических уравнений по параметрам. В качестве примера рассмотрим уравнение Эйлера — Пуассона — Дарбу (Э-П-Д) с равными параметрами:

А А

Е( А, А; и) ° иу +-их--иу = 0 (12)

у - х у - х

в области С = {(х, у)е Я2 : у > х}, где А — действительный параметр [3]. Зная решение

уравнения Э-П-Д с положительными параметрами, построим его продолжение для отрицательных параметров.

Хотя уравнение Э-П-Д исследовано достаточно подробно и найдено общее решение при различных значениях параметров А, В [4], нам важно показать возможность построения продолжения решений уравнения (12) с помощью нового метода потенциалов сопряжённых пар. При 0 < А < 1 решения уравнения (12) хорошо известны [3]; решения этого

уравнения будем обозначать и(А, А, х, у). Функция у(х, у) = (у - х)24 будет одним из решений сопряжённого уравнения

Е» ° Vx7 -

Г А ^ Г А ^

+

у- х

■ V

0.

у-х I

A B

Полагая в (1) a(х, у) =-, ^х, у) =--, с(х, у) = 0, убеждаемся, что условия

у-х у-х

теоремы 1 выполняются (в том числе и условие ах = Ьу), находим коэффициенты двойственного уравнения вида (11):

а^ = А. ((у-х)2А)у = ь-Л = . А ((у-х)2А)х А

V у - х (у - х)2 А у - х V у - х (у - х)2 А у - х и записываем двойственное уравнение:

Е(- А, - А; Ц) ° и„ + Ь4) их - Цу = 0. (13)

у - х у - х

Двойственное уравнение (13), так же как и исходное (12), является уравнением Э-П-Д с равными, но отрицательными параметрами, и потенциал Ц(х, у) является решением двойственного уравнения.

Система уравнений (6) и формула (7), определяющие потенциал Ц(х, у) сопряжённой пары и(А, А, х, у), ^(х, у) = (у - х)24, в данном примере принимают вид

их = ( у - х)2 Ч,

(14)

иу = -(у - х)2 ч,

( х, у)

и (х, у) = Ц(0,0) + | (Г - з)2 Аизёз - (г - з)2 Аи(Л. (15)

(0,0)

Итак, двойственное уравнение (13) является уравнением Э-П-Д с равными отрицательными параметрами. Потенциал Ц(х, у), вычисленный по формуле (15) и являющийся решением двойственного уравнения, будет продолжением решения уравнения Э-П-Д для отрицательных параметров -1 < А < 0. Преобразование (14), (15), в результате которого получается потенциал Ц(х, у), меняет знак каждого параметра уравнения Э-П-Д на противоположный: Ц(х, у) = и(-А, -А, х, у).

До сих пор мы рассматривали гиперболическое уравнение (1), записанное в первом каноническом виде [3]. Пусть теперь в односвязной области В плоскости ХОУдано гиперболическое уравнение во втором каноническом виде

А (и) ° их - иуу + а(х У)их + ь(х У)иу + с(х У)и = ^ с( х, у) е О(П), а( х, у), Ь(х, у) е С>( В).

Сопряжённое уравнение имеет вид

А* (^ - ^ - (av)x - (Ьv)у + ^ = 0.

Как и для уравнения (1), "и(х,у), у(х,у)е 0(П)П 02(П) верно тождество Грина [6]

(и) - и1*\(V) = (шх - тх + аиУ)х -(шу - жу - Ьш)у.

Если теперь Ц(и) = 0, £\(у) = 0, то (уих - иух + аиу)х = (уиу - иуу - Ьиу) . Поэтому выражение (уиу - иуу - Ьиу)ёх + (уих - иух + аиу)ёу является полным дифференциалом в области О. Функция и(х, у), определённая в области О таким образом, что ёи = -(уи - иуу - Ьиу)ёх- (уих - иух + аиу)ёу, также названа нами потенциалом сопряжённой пары и(х, у) и у(х, у). Потенциал и(х, у) определяется системой уравнений

их = иУу - уиу + Ьиу, иу = иУх - уих - аиу

или криволинейным интегралом, не зависящим от формы дуги интегрирования

(х, у)

и(х, у) = и(х0, у0) + | (иу - уи( + Ьиу)ёв + (иу - уи5 - аиу)А.

(xo, у>)

Если а + Ьх = 0 в области О, то такой потенциал и(х, у) является решением двойст-

ух

венного уравнения

их - иу +1 а - 2 V и +

Г V > Ь+2

V у У

иу = 0.

Потенциалы сопряжённых пар имеют важное прикладное значение в математических моделях, которые описываются задачами для линейных гиперболических уравнений второго порядка с двумя переменными [7]. В качестве простого примера рассмотрим математическую модель процесса распространения плоских звуковых волн (малых возмущений) вдоль оси хв покоящейся среде [1]. Обозначим через у время, и(х, у) — скорость возмущённой среды, р(х, у) — давление в этой среде. Постоянные г и с характеризуют свойства покоящейся среды: г — плотность, с — скорость звука.

В рамках этой модели задача нахождения функции и(х, у) сводится к решению одномерного волнового уравнения [2]

4(и) ° с2 • их - иу = 0. (16)

Сопряжённое уравнение имеет вид

Ц» ° с2 • - Ууу = 0. (17)

Волновой оператор Ь0(и) в (16) является самосопряжённым: Е0(у) = Ь0(у) .

Выберем в области В решение сопряжённого уравнения (17) ^х, у) ° г, где г — постоянная: плотность покоящейся среды. В соответствии с рассматриваемым методом потенциалов сопряжённых пар составим тождество Грина (полагая ^х, у) ° г)

V • А(и) - и • =( с 2 • г • их ) х-( г • иу ) У .

Если Ц,(и) = 0, С0(V) = 0 , то (с2 • г• их)х = (г• иу) . Поэтому выражение (г• иу)ёх + (с2 • г • их)ёу является полным дифференциалом в области В. Потенциал Ц(х, у) сопряжённой пары и(х, у) и ^ х, у) ° г в области В здесь определяется полным дифференциалом <Ш = -(г• иу)ёх — (с2 • г• их)ёу и является скалярным потенциалом векторного поля

Р(х у) = (г • иу; с2 •г • их).

Потенциал Ц(х, у) можно определить либо системой уравнений

'Ц- = -г • иу,

и У (18)

иу =-г • с • их ,

либо криволинейным интегралом, не зависящим от формы дуги интегрирования.

( х, У)

I

(xb; У0)

и(X; у) = и(х^ У))- I г• и • с2 •г• из • &.

Система уравнений (18), которой мы определили потенциал сопряжённых пар ц(х, у), описывает распространение плоских звуковых волн (малых возмущений) в покоящейся среде и называется системой уравнений акустики [1]. В литературе по математической физике её принято записывать в виде

иу +1 •их = 0,

г

иу + г• с2 • их = 0,

причем потенциал и(х, у) в этой системе выполняет роль давления р(х, у) в этой среде. Первое уравнение системы (18) продифференцируем по х, а второе — по у

ихх = -г • иух, , ч

и Г (!9)

иуу =-г • с • иху.

Из системы уравнений (19) исключим иху(х, у) и получим уравнение, которому удовлетворяет потенциал Щх, у):

Потенциал Щ(х, у) (давление р(х, у)) является решением двойственного уравнения (20). Двойственное уравнение совпадает с исходным уравнением (16), поэтому исходное уравнение (16) является самодвойственным.

Следовательно, в рассматриваемой математической модели, описываемой уравнениями акустики, потенциал сопряжённых пар Щ(х, у) имеет важный прикладной смысл: потенциал Щ(х, у) равен давлению р(х, у) в среде.

Итак, введённые автором понятия потенциала сопряжённой пары и двойственного уравнения позволили обосновать новый метод нахождения интегральных представлений решений и продолжения решений по параметрам для гиперболических уравнений второго порядка. Рассмотрена математическая модель, описываемая уравнениями акустики, в которой раскрыт важный прикладной смысл потенциала сопряжённых пар.

1. Годунов С. К. Уравнения математической физики. М.: ГИФ-МЛ, 1979. 392 с.

2. Корн Г., Корн Т. Справочник по математике. М.: Наука, 1973. 720 с.

3. Кошляков Н. С., Глинер Э. Б., Смирнов М.М. Дифференциальные уравнения математической физики. М.: ГИФ-МЛ, 1962. 767 с.

4. Нахушев А.М. Уравнения математической биологии. М.: Высшая школа, 1995. 304 с.

5. О свойствах линейных дифференциальных операторов, связанных с полными дифференциалами / Ю. И. Фёдоров // Тез. докл. на 3-й Международной конференции, посвященной 85-летию чл.-корр. РАН, профессора Л. Д. Кудрявцева. М.: МФТИ, 2008. С. 341-344.

6. Тихонов А. Н, Самарский А. А. Уравнения математической физики. М.: Наука, 1972. 736 с.

7. Фёдоров Ю. И. Математические аспекты динамики сорбции газов // Известия Оренбургского государственного аграрного университета. 2014. № 3 (47). С. 46-48.

1. Godunov S. K. Uravnenija matematicheskoj fiziki. M.: GIF-ML, 1979. 392 s.

2. Korn G., Korn T. Spravochnik po matematike. M.: Nauka, 1973. 720 s.

3. Koshljakov N. S., Gliner E. B, Smirnov M. M. Differentsial'nye uravnenija matematicheskoj fiziki. M.: GIF-ML, 1962. 767 s.

4. Nahushev A. M. Uravnenija matematicheskoj biologii. M.: Vysshaja shkola, 1995. 304 s.

5. O svojstvah linejnyh differentsial'nyh operatorov, svjazannyh s polnymi differentsialami / Fjodorov Ju. I. // Tez. dokl. na 3-j mezhdunarodnoj konferentsii, posvjashchjonnoj 85-letiju chl.-korr. RAN, professora L. D. Kudrjavtseva. M.: MFTI, 2008. S. 341-344.

6. Tihonov A. N., Samarskij A. A. Uravnenija matematicheskoj fiziki. M.: Nauka, 1972. 736 s.

7. Fjodorov Ju. I. Matematicheskie aspekty dinamiki sorbtsii gazov // Izvestija Orenburgskogo gosu-darstvennogo agrarnogo universiteta. 2014. № 3 (47). S. 46-48.

c2-U - U = 0.

** jy

(20)

СПИСОК ЛИТЕРАТУРЫ

REFERENCES