Математика к Математическое
моделирование
Ссылка на статью: // Математика и Математическое моделирование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2016. № 02. С. 1-8.
Б01: 10.7463/шаШш.0216.0843737
Представлена в редакцию: 12.03.2016 Исправлена: 28.03.2016
© МГТУ им. Н.Э. Баумана
УДК 517.95
К задаче о наклонной производной для уравнения Лаврентьева - Бицадзе в полуплоскости
Алгазин О. Д.1'", Копаев А. В.1
1МГТУ им. Н.Э. Баумана, Москва, Россия
Уравнение Лаврентьева - Бицадзе решается в полуплоскости. При этом областью эллиптичности также является полуплоскость, а областью гиперболичности - полоса. На одной из прямых, ограничивающих полосу, задана наклонная производная в направлении характеристики, а на другой прямой - границе раздела полосы и полуплоскости - решения сопрягаются краевыми условиями четвертого рода.
Ключевые слова: уравнение Лаврентьева - Бицадзе, наклонная производная, краевые условия четвёртого рода
Введение
Многие проблемы прикладного характера (например, проблемы околозвукового течения сжимаемой среды) сводятся к решению дифференциальных уравнений смешанного типа (уравнению Трикоми, уравнению Чаплыгина и др.). Важной моделью уравнения смешанного типа является уравнение Лаврентьева-Бицадзе [1]. Решению различных краевых задач для этого уравнения посвящено огромное количество статей (см., например, [2], [3], [4], [5], [6]). Поэтому и решение задачи о наклонной производной для уравнения Лаврентьева Бицадзе представляется весьма актуальным.
1. Обозначения
Введём следующие обозначения:
прямая, параллельная оси ; и - полуплоскости, на которые прямая
разбивает плоскость .
В полуплоскости Б (Л > 0 ) рассмотрим уравнение Лаврентьева - Бицадзе
(1)
В полуплоскости Б + - области эллиптичности уравнения (1) - положим
и (х ;у) = и+(х;у) ,
а в полосе С = Б_ П Б = { ( (х; у) Е Е2 : — к < у < 0 ) } - области гиперболичности уравнения (1) - положим
и(х;у) = и~(х; у)
2. Задача 1 о наклонной производной
Пусть q (х) - функция, удовлетворяющая на числовой прямой условию Гёльдера (включая бесконечно удалённую точку) [7, с. 20], при этом пусть q ( оо ) = 0 , а её производная - функция р (х) = q ' (х) - также удовлетворяет на числовой прямой условию Гёль-
дера, включая бесконечно удалённую точку, и р ( х) < —— ( а > 1 , | х | > Е) .
\х\
Рассмотрим следующую задачу.
Найти функцию и+(х ;у) , гармоническую и ограниченную в полуплоскости Б + и функцию , ограниченную и удовлетворяющую в полосе уравнению
и;* — и;у = о (2)
по краевому условию
ди~ ди~
(х; —Л)---— (х; —Л) = р(х)
дх ду
(по заданной на прямой I _ ^ наклонной производной в направлении характеристик уравнения (2), заданных уравнениями х+у = С ) и по условиям сопряжения функций и+( х ;у) и на оси
и_ (х; 0) = и+ (х; 0) (х; 0) = к^ (х ; 0 ) (3)
где - положительное число.
Решим эту задачу. Представим функцию и _( х ;у) по формуле Даламбера
и~(х;у) = /(х + у) + д(х - у) где и - дважды непрерывно дифференцируемые функции. Преобразуем функцию и _( х ;у) :
и~ (х;у) = /(х + у) + д(х + у) - д(х + у) + д(х - у) =
гХ-у гХ-у
= /(х + у) + д(х + у) + = т(х + у) +
^х+у х+у
где введены обозначения
т(х) = /(х) + д(х),/л(х) = д'(х)
Имеем
ди~
—- (х; у) = т (х + у) + ц(х - у) - ц(х + у)
ди~
—— (х; у) = т (х + у) - ц(х - у) - ц(х + у) ду
ди~ ди~
(х;у) ~ "аТ"(х;у) = 2ц(х ~ у)
ди~ ди~
—— (х; -К)---— (х; -К) = 2ц(х + К) = р(х)
ох Зу
1 1 Отсюда получаем, что ц (х + Л) = - р ( х) , а ц ( х) = - р ( х — Л) . Далее,
1 г*-^
{Г(х;у) = т(х + у) + - р(^-Л)^ и~(х;0) = т(х)
Но тогда в силу условий сопряжения и+( х ; 0 ) = т ( х) . А так как функция и+( х ;у) является гармонической и ограниченной в полуплоскости , то она представима в этой полуплоскости интегралом Пуассона
{/+(*; у) = - [
- X)2 + у2
Далее,
ди~ 1 1
(х; у) = т (х + у) - - р(х - у - Л) - - р(х + у - Л)
ди~
—— (х; 0) = т (х) - р(х - Л) ду
ди+ ± [+сот^т_х)2_у2щ
1 Г
.(Х;у) = -
Зу ' к ]_т [^-х)2+у2]2
Применим к последнему интегралу формулу интегрирования «по частям»:
ди+ 1 е-х . 1 Г+0°т'(О(Т-хЖ - (х-у) =--т(Л--- + - -—
Зу ' 7Г - х)2 + у2 °° 7Г }_т (( - х)2 + у2
Так как функция ограничена в полуплоскости , то функция
ограничена на числовой прямой, и первое слагаемое в правой части последнего равенства равно 0. Тогда
ди+ 1
1 Г
■(х;у) = -
Зу ' 7Г }_т (( - х)2 + у2
—— (х; 0) = - —-
ду "к ¿-со К %
Здесь мы должны предположить, что функция т ' ( х) удовлетворяет на числовой прямой условию Гёльдера (включая бесконечно удалённую точку). А так как ещё ограниченная функция т (х) = /^т ' ( й^ + С , где С - действительное число, то т ' ( оо ) = 0 . Тогда для нахождения неизвестной функции получаем сингулярное интегральное уравне-
ние
к Г+0°т'(0^
к Г
т (х) — р(х — К) = —
п 1-е
Рассмотрим интеграл типа Коши
М
1т
2т }_т % - г {Ф (г),/шг<0 Применяя формулы Сохоцкого для интеграла типа Коши [7, с.47 ]
Ф+(х) - ф-(х) = т'(х),
Ф+(х) + Ф~(х)
_ 1 Г
-оэ % Х
получим задачу о скачке на действительной оси:
Ф+(х) — Ф_(х) — р (х — к) = к / (Ф+(х) + Ф_(х) ) (4)
Перепишем краевое условие (4) в виде
Ф+(х)(1 - Ы) - ф-(х)(1 + Ы) = р(х - К) Рассмотрим ещё один интеграл типа Коши
- К)^ _ 1т г > 0;
1т
+ 00 ,
2тп }_т [У (г), 1т г < 0.
Теперь с учётом формул Сохоцкого получаем
Ф+(х) (1 - Ы) - ф-(х)(1 + Ы) = Ч*+(х) - ^-(х)
или
Ч*+(х) - Ф+(х)(1 - Ы) = 4>~(х) - ф-(х)(1 + /а) А отсюда следует, что
Ч*+(г) - Ф+(г) (1 - Ы) = С (1т г > 0) ¥~(г) - Ф~(г)( 1 + Ы) = С (1т г < 0) где - произвольная (комплексная) постоянная. Но тогда
¥+(г) — С ¥~(г) — С
ф+(г) = ТГПИГ(/т2 > 0)'Ф_(г) = ТТТйГ(/т2 <0)
Теперь можем найти неизвестную функцию т'(х) :
Ч*+(х) — С Ч*~(х) — С
т'(х) = Ф+(х)-Ф-(х) =
(1 - /а) (1 + /а) 1 1 1
11 1 Г+сор((-к)с1( 2 ЫС
)
1 + /с/ 2 27П }_т % — х 1 + к2
р(х-}I) к [+сор(^ -К)й$ 2ЫС
[
1 + к2 (1 + к2)л ]_т %-х 1+ к2 Это возможно тогда и только тогда, когда - число чисто мнимое, т. е. когда , где - действительное число. А учитывая, что и , получаем, что
, а
р(х — К) к Г+0°р(^ -К)й$
*'(*) = , , +
1"
1+/С2 (1 + /С2)ТГ ^
Наконец, мы можем найти неизвестную функцию т (х) :
а(х-К) к Г+со
ю = ^ЩГ - (ТТР^ ]_„ р(? " Л) 1п|?"х| ^ + С
где - произвольная (действительная) постоянная. Применим к последнему интегралу формулу интегрирования «по частям»:
^.+00 /-+00 - h) lnlf -x№ = qrC-Г - h) lnlf ~x\\±%-
J — CQ J — СО
{-X
Так как функция q ( х) удовлетворяет в бесконечно удалённой точке условию Гёль-дера и так как q ( оо ) = 0 , то первое слагаемое в правой части последнего равенства равно 0, и
т (х) = ^ + -V /++ C
4 у 1+к2 (1 + к2)п •'-оо ^—х
Рассмотрим пример для
p (х) = (а2 + х 2 ) 2 = (_ 2 (а 2 + х2 ) ) .
Имеем:
(а2+ х2)2 V 2(а2+ х2)у 1
2(1 + /с2)(а2 + (х- К)2) к f+co
т(х) =
2(1+ к2)п }_т (а2 + a-/i)2)(^-x)
+ С =
1 к(х - h)
+ -. ^ п .-гтттг + С =
2(1 + к2){а2 + (х - /i)2) 2а(1 + к2){а2 + (х - /i)2)
k(x-h)-a
2а(1+ fe2)(a2+ (x-h)2)
+C
2a(l+ fe2)((x-h)2+ (y+a)2)
jj-r л = 2k(x+y-h)+a(k2-l)___1 , r
( 'у ) 4 а ( 1 + fe2 ) ( а 2 + (х+у - h) 2 ) 4 ( а 2 + (х - у - h) 2 )
3. Задача 2 о наклонной производной
Аналогично решается следующая задача. Найти функцию tt+(х;у) , гармоническую и ограниченную в полуплоскости D + и функцию tt _ ( х 'у) , удовлетворяющую в полосе G уравнению (2) по краевому условию
dU~ dU~
- -Г— (х; -Л)--— (х; -Л) = р(х)
ох оу
(по заданной на прямой I_ h наклонной производной в направлении характеристик уравнения (2), заданных уравнениями х — у = С ) и по условиям (3) сопряжения функций tt+(х ' у) и Í/_ (х ' у) . Только при решении задачи 2 используется представление
tr (X' у) = т (х — у) + /;_+уУ ц( О df , где т (х) = f (х) + # (х) , а ц(х) = /' (х) = —1 р (х + Л) . В результате получаем
q(x + h) к í+coq^ + h)d£
TW" 1 + fc2 nil+k2))^ f-x +
n+f л у f+0° *(№
U (x;y) = - —-
n Loo (f - x)2 + У2
I гХ+у
U~(x;y) = т(х -у) - - I + Л.) df
Заключение
Таким образом, в работе решена задача о наклонной производной для уравнения Лаврентьева-Бицадзе в полуплоскости, при этом область эллиптичности также является полуплоскостью, а область гиперболичности - полосой. Наклонная производная задана на одной из прямых, ограничивающих полосу, в направлении характеристики гиперболического уравнения, а на другой прямой, ограничивающей полосу и разделяющей ее с полуплоскостью эллиптичности, выполняются краевые условия четвертого рода. Показано, что при сделанных предположениях относительно заданных функций задача имеет единственное решение (с точностью до произвольной действительной постоянной).
Список литературы
1. Бицадзе А.В. К проблеме уравнений смешанного типа // Тр. МИАН СССР. 1953. Т.41. С. 3-59.
2. Сабитов К.Б., Новикова В.А. Нелокальная задача А.А. Дезина для уравнения Лаврентьева-Бицадзе // Известия вузов. Математика. 2016. № 6. С. 61-72.
3. Moiseev E.I., Moiseev T.E., Vafodorova G.O. On an Integral Representation of the Neu-mann-Tricomi Problem for the Lavrent'ev-Bitsadze Equation // Differential Equations. 2015. Vol. 51. No. 8, pp. 1065-1071. DOI: 0.1134/S0012266115080108
4. Солдатов А.П. О задачах типа Дирихле для уравнения Лаврентьева-Бицадзе. Дифференциальные уравнения и динамические системы. Сборник статей // Тр. МИАН. 2012. Т. 278. С. 242-249.
5. Сабитов К.Б., Хаджи И.А. Краевая задача для уравнения Лаврентьева-Бицадзе с неизвестной правой частью // Известия вузов. Математика. 2011. № 5. С. 44-52.
6. Сербина Л.И. Решение одной начально-краевой задачи теории фильтрации с нелокальными краевыми условиями // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2003. Вып. 19. С.16-21. DOI: 10.14498/vsgtu133
7. Гахов Ф.Д. Краевые задачи. М.: Наука. 1977. 640 с.
Mathematics à Mathematical Modelling
Mathematics and Mathematical Madelling of the Bauman MSTU, 2016, no. 02, pp. 1-8.
DOI: 10.7463/mathm.0216.0843737
Received: Revised:
12.03.2016 28.03.2016
© Bauman Moscow State Technical Unversity
Towards the Oblique Derivative Problem for the Lavrentyev - Bitsadze Equation in the Half-Plane
O.D. Algazin1'*, A.V. Kopaev1
mopi66igyandex.ju
1Bauman Moscow State Technical University, Moscow, Russia
Keywords: the Lavrent'ev - Bitsadze equation, oblique derivative, boundary conditions of the fourth kind
The article solves the boundary value problem with an oblique derivative for the Lavrentyev - Bitsadze equation, which is an equation of the mixed elliptic-hyperbolic type. The mixed type equations are used in transonic gas dynamics. A restricted solution to the problem is sought in the half plane, consisting of the upper half plane (where the equation is elliptic) and adjacent band (where the equation is hyperbolic). At the boundary of the domain an oblique derivative towards one of the characteristics is specified. At the boundary of the upper half plane, which is a line of equation type alteration, the matching conditions of the fourth kind are set. In the hyperbolicity band the solution is represented by d'Alembert formula, and in the upper half plane, where the equation is elliptic, the restricted solution is represented by Poisson integral of unknown density. For unknown Poisson integral density a singular integral equation is obtained, which is reduced to the Riemann boundary value problem for holomorphic functions. The solution of this problem is obtained in an explicit form. Thus, the solution to the problem with an oblique derivative for the Lavrentyev - Bitsadze equation was obtained in an explicit form for the case of the half plane accurate to a constant summand. An example of solving the problem to prove the theoretical calculations is provided at the end of the article.
References
1. Bitsadze A.V. On the problem of equations of mixed type. Trudy MIAN SSSR, 1953, vol.41, pp. 3-59. (in Russian).
2. Sabitov K.B., Novikova V.A. Dezin's problem for Lavrent'ev-Bitsadze equation. Izvestiya vuzov. Matematika, 2016, no. 6, pp. 61-72. (In Russian).
3. Moiseev E.I., Moiseev T.E., Vafodorova G.O. On an Integral Representation of the Neu-mann-Tricomi Problem for the Lavrent'ev-Bitsadze Equation. Differential Equations, 2015, vol. 51, no. 8, pp. 1065-1071. DOI: 10.1134/S0012266115080108
4. Soldatov A.P.On Dirichlet-type problems for the Lavrent'ev-Bitsadze equation. Trudy MIAN = Trudy Matematicheskogo Instituta Steklova. 2012, vol., 278, pp. 242-249. (In Rus-
sian). (English version of journal: Proceedings of the Steklov Institute of Mathematics, 2012, vol. 278, no. 1, pp. 233-240. DOI: 10.1134/S0081543812060223).
5. Sabitov K.B., Khadzhi I.A. The boundary-value problem for the Lavrent'ev-Bitsadze equation with unknown right-hand side. Izvestiya vuzov. Matematika. 2011, no. 5, pp. 44-52. (In Russian). (English version of journal: Russian Mathematics, 2011, vol. 55, no. 5, pp 35-42. DOI: 10.3103/S1066369X11050069).
6. Serbina L.I. Saving one initial boundary value problem of filter theory with non-local boundary conditions. Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta. Ser. Fiziko-matematicheskie nauki, 2003, vol. 19, pp.16-21. (In Russian). DOI: 10.14498/vsgtu133
7. Gakhov F.D. Kraevye zadachi [Boundary problems]. Moscow, Nauka Publ., 1977. 640 p.(In Russian).