УДК 517.956
В.А. Киричек, Л.С. Пулькина1
ЗАДАЧА С ДИНАМИЧЕСКИМИ ГРАНИЧНЫМИ УСЛОВИЯМИ ДЛЯ ГИПЕРБОЛИЧЕСКОГО УРАВНЕНИЯ
Рассматривается начально-краевая задача с динамическим граничным условием для гиперболического уравнения в прямоугольнике. Динамическое граничное условие представляет собой соотношение, в которое помимо значений производных искомого решения по пространственным переменным входят производные первого порядка по переменной времени. Основной результат статьи состоит в обосновании разрешимости поставленной задачи. Доказано существование единственного обобщенного решения. Доказательство базируется на полученных в работе априорных оценках, методе Галёркина и свойствах пространств Соболева.
Ключевые слова: динамические граничные условия, гиперболическое уравнение, обобщенное решение.
Введение
Математическое моделирование колебательных процессов различной природы обычно приводит к краевым задачам для гиперболический уравнений, граничные условия в которых отражают способы закрепления. Если требуется учесть сопротивление среды, например, при наличии некоторого демпфирующего устройства, то граничные условия должны содержать производную по переменной времени. Такие условия принято называть динамическими граничными условиями. Постановка простейшей задачи с динамическим условием приведена в [1]. Эта задача состоит в нахождении решения уравнения колебаний струны, конец которой испытывает сопротивление среды, пропорциональной скорости ее движения [1]. При этом возникает динамическое граничное условие ких(1,Ь) = —ащ(1^). Задачи с динамическими граничными условиями позволяют резрешить вопросы в таких областях, как генетика, медицина, физика, особенно в задачах, свзяанных с акустикой [2]. Задачи с динамическими граничными условиями возникают при изучении различных процессов: процесса переноса газа сквозь мембраны с учетом физико-химических процессов на поверхности [2], распространения волн в стратифицированных и вращающихся жидкостях [3], различных колебательных процессов [4-6]
В работе рассмотрена задача для гиперболического уравнения с динамическими условиями, содержащими производную по времени первого порядка, и доказано существование единственного обобщенного решения. Ранее эта задача изучалась в [8], но в предлагаемой работе удалось существенно ослабить условия на входные данные, обеспечивающие однозначную разрешимость задачи.
1. Постановка задачи
Рассмотрим уравнение
иа - (а(х,г)их)х + о(х,г)и = /(х,г), (1.1)
где а(х,Ь) > 0 всюду в (т = (0,1) х (0,Т), и поставим задачу: найти в области (т решение уравнения (1.1), удовлетворяющее начальным данным
и(х, 0) = ф(х), щ(х, 0) = ф(х) (1.2)
и граничным условиям
их (0,г) = а(г)щ(0,г), их(1,г) = ¡з(г)щ(1,г). (1.3)
Киричек В.А., Пулькина Л.С., 2017 Киричек Виталия Александровна ([email protected]), Пулькина Людмила Степановна ([email protected]), кафедра уравнений математической физики, Самарский национальный исследовательский университет имени академика С.П. Королева, 443086, Российская Федерация, г. Самара, Московское шоссе, 34.
Обозначим
г0 = {(х,г)-. х = 0,г е [0,т]}, г1 = {(х,г): х = 1, 4 е [0,т]},
Г = Го и гь Ш(От) = {и : и е ), Щ е Ь2(Г)},
№(От) = {V : V е Ш(Ют), у(х, Т) = 0}. Введем понятие обобщенного решения, для чего, следуя известной процедуре [7], сначала выведем равенство
т I I
J !\-UtVt — аитух + cuv)dxdt = ^ ф(х^(х, 0)dx+
0 0 о
т т
^У a(0,t)a(t)ut(0,t)v(0,t)dt — ^ а(1,г)/3 (t)ut(l,t)v(l,t)dt+
о о
т I
+П/(х^ (1-4>
оо
которое является результатом интегрирования соотношения гиЬи = /V по области Ют.
Определение. Обобщенным решением задачи (1.1)—(1.3) будем называть функцию и(х,€) е Ш(От), удовлетворяющую условию и(х, 0) = ф(х) и интегральному тождеству (1.4) для любой функции v(x,t) е Ш(От)
2. Разрешимость задачи
Теорема. Пусть выполняются следующие условия:
с е с (От), а^ е с (О т), а е с (О т), / е Ьх(От), а, в е с2[0,т],а(^ > 0,№) < 0 Ш е [0, Т].
Тогда существует единственное обобщенное решение задачи (1.1)—(1.3). Доказательство
Для доказательства единственности покажем, что функция и(х,Ь), которая представляет собой разность двух обобщенных решений задачи и(х,Ь) = и\(х,Ь) — и2(х,Ь), всюду в области Qт равна нулю. Очевидно, что и(х, 0) = 0 и выполняется тождество
т I
(—utvt — auxvx + cuv)dxdt =
оо
т т
= ! a(0,t)а(t)ut(0,t)v(0,t)dt — ¡^ a(l,t)в(t)ut(l,t)v(l,t)dt. оо Выберем в качестве v(x,t) в этом тождестве функцию
v{x,t)= \ I и(х,п)^; 0 < t < т, (2.1)
1 0; т < г < Т,
где т е [0, Т] и произвольно. Легко видно, что эта функция принадлежит пространству ) и vt = и.
В результате итегрирования по частям получим
I
J(у^(х, т) + a(x, 0^Х(х, 0)^х+
о
т т
a(0,t)а(t)vХ(0,^(0,^ — 2J a(l,t)в(t)vХ(l,t)dt =
оо
т т
= ! (a(0,t)а(t))ttvХ(0,t)dt — ^ (a(l,t)в(t))ttvХ(l,t)dt—
т I т I
а^х + 2 /
0 0 0 0 +ю2(0, 0)(а(0,1)а(1))г\г=0 - ю2(1, 0)(а(1,г)в(г))г\г=0. (2.2)
Сделаем некоторые оценки. Заметим, что в силу условий теоремы существует такое положительное число
С0, что тах\с(х,Ь)\ ^ С0 Применяя неравенство Коши, получим
Ят
т I т I
2\ ! J ^ С0 J J(V2 + vх)dxdt.
0 0 0 0
Рассмотрим теперь два первых слагаемых правой части равенства (2.2). Для того, чтобы их оценить, воспользуемся неравенствами
I I
,2/„ I 2 I „2/
ю2(0,Ь) < 2^ VI(х, Ь^х v2(x,t)dx,
' I
00 I I
,2/„ I 2 I „2/
ю2(1,Ь) ^ V12,(x,t)dx + 2 J v2(x,t)dx,
I
00 которые легко следуют из представлений [6]
0 I
= / „(е,т+ФЛ = / „+ФЛ
хх
С помощью этих неравенств, учитывая также условия теоремы, гарантирующие существование положительных чисел А0, В0 таких, что \(а(0, Ь)а(Ь))м\ ^ А0 , \(а(1 ,Ь)в (Ь))и\ ^ В0 , получим
т т I т I
\ ! v2(0,t)(a(0,t)a(t))ttdt\ < А0(2^ J v2xdxdt + 2 ^ J v2 (x,t)dxdt), 0 0 0 0 0 т т I т I
у! ю2^^)^^)^))^^ < В0(2^ ! v2dxdt + 2 ! J v2(x,t)dxdt). 0 0 0 0 0 Используя условия теоремы и неравенства [7]
I
ю2(0, 0) (ею2х(х, 0) + с(е)ю2(х, 0))dx,
0
I
ю2(1, 0) (ею2х(х, 0) + с(е)ю2(х, 0)^х,
0
получим оценки внеинтегральных слагаемых
I I
\ю2(0, 0)(а(0,Ь)а(Ь))4\4=0\ < А(е^ ю2х (х, 0)dx + с(е)^ ю2(х, 0^х),
00 I I
\ю2(1, 0)(а(1,Ь)/3(Ь))^=0\ < В(е ! ю2х(х, 0)dx + с(е)^ ю2(х, 0)dx).
х
00 Благодаря полученным оценкам придем к неравенству
I т I
Х(х т) , а(х 0)ю1 (х, 0))с1х < сЛ (и2 + ю2х)
J(ux(x ,т) + а(х, 0)юхх (х, 0))dx ^ с\ J ^(и2 + ю2х )dxdt+ 0 0 0 I I
+С2еjv2x(x, 0)^х+сз/^ ^
где С1 = max{2l(Aо + Во), 2(Ао+Во) ,со}, сх = (А + В), сз = схс(е).
Будем считать, не ограничивая общности, что a(x,t) ^ aо > 0, и выберем е = . Тогда aо — схе = аХг > 0.
Перенеся теперь второе слагаемое правой части в левую, получим
I т I I
J(u2(x,т) + ооvX(х, 0))dx * с1 J J(и2 + vX)dxdt + сз ^ у2(х, 0)dx. о о о о
В силу представления (2.1)
т
"2(х 0) * т1иъж
мы приходим к неравенству
I т I
1(и2(х,т) + у vX(х, 0))dx * сА ! J(и2 + у2х)dxdt, (2.3)
2
о о о
Х
где с4 = с1 + сзт2.
Введем вспомогательную функцию
t
w(x,t) = J их(х,
о
Из этого представления следует, что
ух(х, ^ = ,ш(х, ^ — ,ш(х, т), Ух(х, 0) = —,ш(х, т).
Тогда
уХ(х^) * 2т2(х, ^ + 2и)х(х,т).
Пользуясь произволом, выберем т так, чтобы аХ0 — 2с4т > 0 Пусть для определенности аХ0 — 2с4т ^ а0. Тогда для любого т е [0, ]
I т I
Х Х Х Х
1(и+*с5/1(и2+
оо
где 5 = 2с4/ шш{1, }. Применив лемму Гронуолла, получим, что и(х,т) = 0 для любого т е [0; ]. Если рассмотреть теперь задачу с начальными данными на t = , то в результате тех же рассуждений, что и выше, докажем, что и(х,т) = 0 при т е [0; ]. Продолжив этот процесс, за конечное число шагов получим, что и = 0 во всем цилиндре От.
Для доказательства существования обобщенного решения применим метод компактности. Начнем с построения последовательности приближенных решений поставленной задачи методом Галеркина. Рассмотрим последовательность {тп(х)}, где тп е С2(П) П С({}), которая является линейно независимой и образует полную систему в Ш^П). Будем искать приближенные решения задачи в виде
т
ит(х^) = У1 dk (^юк(х) (2.4)
к=1
из соотношений
I
J (итт Wj + + сит^з )dx + a(0,t)а(t)umí(0,t)wj —
I
—о(1, ^в^и1^,^ = ! ^dx, (2.5)
о
которые представляют собой систему дифференциальных уравнений второго порядка относительно с!к ^) :
т т т
XX (Мко +^)Вкз ^к + dk (^Скз =
к=1 к=1 к = 1
= /з ^), (2.6)
I
где Акз = / wk (х^з (х^х,
о
Ву = а(Ь)а(0,Ь)тк(0)т^ (0) - в(Ь)а(1,Ь)тк(1)т^ (I) I
Ск0 = Ла(х,Ь)тк (х)т' (х) + с(х,Ь)тк (х)т^ (х)]Сх,
0
I
и /] (Ь) = / /т^ (x)dx.
0
Добавим начальные условия
dk(0)= , С'к (0) = 5к,
т т
где 7к и 5к — коэффициенты сумм фт(х) = ^ 7ктк(х), фт(х) = ^ 5ктк(х), аппроксимирующих при
к=1 к=1
т эти суммы аппроксимируют ф(х) и ф(х) в нормах Ш2(0,1) и Ьх(0,1) соответственно: фт(х) ^ ф(х)
и фт(х) ^ ф(х). Таким образом, мы получили задачу Коши. Так как матрица Ау является обратимой в силу выбора тк (х), то, учитывая условия теоремы, гарантирующие ограниченность коэффициентов системы, задача Коши однозначно разрешима и С/', € Ь1 (0,Т). Следовательно, последовательность приближенных решений {ит} построена.
Теперь докажем ограниченность этой последовательности в ^^((т). Для этого нужно получить априорные оценки. Умножим (2.5) на (Ь), просуммируем по ] от 1 до т и проинтегрируем по Ь от 0 до т. Получим
т I т I т I
ит ит dxdt + ! J аит^^А + ! J ситит dxdt+
0 0 0 0 0 0
т т т I
+ ! аа(ит(0,Ь))2СЬ - J ва(ит(1,Ь))2СЬ = ^ J /umdt. (2.7)
Н "" J \u't V;")) — J J J u't
0 0 0 0 Проинтегрировав по частям последнее равенство, получим
I
2!((иТ(х,т))2 + а(х,т)(ит(х,т))2)Сх + ^ аа(ит(0,Ь))2Л-00 т I т I
- ! ва(и'1п(1, Ь))2 СЬ =2>1(ит(х, 0))2 Сх +^ ^ а^ит)2СхСЬ+ 0 0 0 0
I т I т I
+ Ц а(ит(х, 0))2Сх - J J сититСхСЬ + ^ J /и^СхЛ. (2.8)
0 0 0 0 0
Применяя представление
(х,т) = J ит& + ит(х, 0),
и
0
получим
I т I I
0 0 0 Так как в силу условий теоремы
J (ит(х, т))2Сх < 2т ^ J(ит)2СхСЬ + 2 J(ит(х, 0))2Сх.
т I т I
<ит? СхСь ^ а111(ит2схж,
0 0 0 0
то приходим к неравенству
I
J((ит(х, т))2 + (ит(х, т))2 + а(х, т)(ит(х, т))2)Сх + 2 J а(0, Ь)а(Ь)(ит(0, Ь))2СЬ-00 т I I
-2^ а(1,Ь)в(Ь)(ит(1,Ь))2СЬ (ит(х, 0))2Сх + ^ а(ит(х, 0))2Сх+
т
I т I т I
+2 J(ит(х, 0))2dx + М ! 1[(ит)2 + (ит)2 + (um)2]dxdt + ! J /2dxdt. о о о о о
Учитывая свойства норм в пространствах Ь2 и ШХ1 и выбор начальных условий задачи Коши, получим
I т
J((um)2 + К™)2 + a(um)2)dx + 2 J aa(t)(um'0,t))2dt-
0
t l
-2 j aß(t)(um(l,t))2dt < M J J''um)2 + 'um)2+ 0 0 0
,m\2nj__n I II-PII2 I лгл|„/,1|2 I ™ll,„l|2
+(их У)сЪсИ + \\/\Ц2 + N (\Щ\12 + т\\ф\\^1).
Так как / является ограниченной функцией в Ьх(От), ф в Ш2 (От), а ф в Ьх(От),то \\/\\2Ь2 + + \\ф\?ш 1 + + \\Ф\\ь2 * К, где К положительная константа. Тогда последнее неравенство примет вид
I т I
у ((ит)2 + ит2 + ю2^ dx * м I у ((ит)2 + ит2 + (и^2)^+к.
о о о
К этому неравенству применима лемма Гронуола, после ее применения и интегрирования по t от 0 до Т получим
т I
((ит)2 + (ит)2 + (<')2Ь=тdxdt * МК(емт — 1). оо
Пусть Мм К(емт — 1) = N, тогда \\^х^^И'ъу1 * N. Значит, последовательность функций ограничена в Шх1, следовательно, можно выделить подпоследовательность, слабо сходящуюся к элементу из этого же пространства, то есть к и е Ш2(От). Покажем, что этот предел и есть обобщенное решение задачи (1.1)—(1.3). Умножим каждое из соотношений (2.5) на свою функцию В^^), В^ (Т) = 0, полученные равенства про-
т
суммируем по ] от 1 до т, а затем проинтегрируем от 0 до Т. Обозначим ц(х,Ь) = ^ В^ (^л^ (х). После
з=1
интегрирования по частям получим тождество
T l
t
00 T
(mm, mm, m ,
'ut nt + aux nx + cu n )dxdt+
+ J aaum'0,t)nmdt -J ßaum'l,t)nmdt
0
T l
m
/nmdxdt,
оо
т
справедливое для любых ц(х,Ь) = ^ В^ (^л^ (х). Переходя к пределу при т ^ то, получим тождество (1.4).
з=1
Заметим, что это тождество выполняется пока лишь для функций ц(х,Ь). Однако известно [7], что множество всех функций такого вида плотно в Ш(От), тождество для предельной функции и(х,Ь) выполняется для любых функций у(х,Ь) е Ш(От).
Следовательно, теорема полностью доказана.
т
Литература
[1] Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука, 1977.
[2] Tobias Louw, Scott Whitney, Anu Subramanian, and Hendrik Viljoen. Forced wave mation with internal and boundary damping // Journal of applied physics, 111,014702 (2012).
[3] Корпусов М.О. Разрушение в неклассических волновых уравнениях. М.: URSS, 2010.
[4] Doronin G.G., Lar'kin N.A., Souza A.J. A hyperbolic problem with nonlinear second-order boundary damping // EJDE, 28, 1-10 (1998).
[5] Бейлин А.Б., Пулькина Л.С. Задача о продольных колебаниях стержня с динамическими граничными условиями // Вестник СамГУ. 2014. № 3(114). С. 9-19.
[6] Пулькина Л.С. Задача с динамическим нелокальным условием для псевдогиперболического уравнения // Известия вузов. Математика. 2016. № 9. С. 42-50.
[7] Ладыженская О.А. Краевые задачи математической физики. М.: Наука, 1973.
[8] Бейлин С.А. Об одной краевой задаче для волнового уравнения // Вестник СамГУ. 2011. № 5(86). С. 12-17.
References
[1] Tikhonov A.N., Samarsky A.A. Uravneniia matematicheskoi fiziki [Equations of mathematical physics]. M.: Nauka, 1977 [in Russian].
[2] Tobias Louw, Scott Whitney, Anu Subramanian, and Hendrik Viljoen. Forced wave motion with internal and boundary damping. Journal of applied physics, 111,014702 (2012) [in English].
[3] Korpusov M.O. Razrushenie v neklassicheskikh volnovykh uravneniiakh [Destruction in nonclassical wave equations]. М.: URSS, 2010 [in Russian].
[4] Doronin G.G., Lar'kin N.A., Souza A.J. A hyperbolic problem with nonlinear second-order boundary damping. EJDE, 28, 1-10 (1998) [in English].
[5] Beylin A.B., Pulkina L.S. Zadacha o prodol'nykh kolebaniiakh sterzhnia s dinamicheskimi granichnymi usloviiami [A problem on longitudinal vibrations of a rod with dynamic boundary conditions]. Vestnik SamGU [Vestnik of Samara State University], 2014, no. 3(114), pp. 9-19 [in Russian].
[6] Pulkina L.S. Zadacha s dinamicheskim nelokal'nym usloviem dlia psevdogiperbolicheskogo uravneniia [A problem with dynamic nonlocal condition for pseudohyperbolic equation]. Izvestiia vuzov. Matematika [Russian Mathematics. (Iz. VUZ)], 2016, no. 9, pp. 42-50 [in Russian].
[7] Ladyzhenskaya O.A. Kraevye zadachi matematicheskoi fiziki [Boundary problems of mathematical physics]. M.: Nauka, 1973 [in Russian].
[8] Beylin S.A. Ob odnoi kraevoi zadache dlia volnovogo uravneniia [On a certain boundary problem for a wave equation]. Vestnik SamGU [Vestnik of Samara State University], 2011, no. 5(86), pp. 12-17 [in Russian].
V.A. Kirichek, L.S. Pulkina2
PROBLEM WITH DYNAMIC BOUNDARY CONDITIONS FOR A
HYPERBOLIC EQUATION
We consider an initial-boundary problem with dynamic boundary condition for a hyperbolic equation in a rectangle. Dynamic boundary condition represents a relation between values of derivatives with respect of spacial variables of a required solution and first-order derivatives with respect to time variable. The main result lies in substantiation of solvability of this problem. We prove the existence and uniqueness of a generalized solution. The proof is based on the a priori estimates obtained in this paper, Galyorkin's procedure and the properties of Sobolev spaces.
Key words: dynamic boundary conditions, hyperbolic equation, generalized solution.
Статья поступила в редакцию 22/1/2017. The article received 22/7/2017.
2Kirichek Vitaliia Alexandrovna ([email protected]), Pulkina Ludmila Stepanovna ([email protected]), Department of Equations of Mathematical Physics, Samara National Research University, Samara, 34, Moskovskoye shosse, 443086, Russian Federation.