CC BY
112
УДК 519.63
О.П. Филатов1
ГЛОБАЛЬНАЯ ТЕОРЕМА СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ РЕШЕНИЯ ПЕРВОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ НЕЛИНЕЙНОГО ИНТЕГРОДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ ПАРАБОЛИЧЕСКОГО ТИПА
Доказана глобальная теорема существования и единственности обобщенного решения первой краевой задачи для нелинейного интегродифференци-ального уравнения параболического типа.
Если правая часть уравнения интегрально ограничена, то имеет место оценка нормы разности двух решений, из которой следуют непрерывная зависимость решения от начальной функции и единственность решения первой краевой задачи.
Рассматриваемая задача обобщает реальные модели измерения уровня несжимаемой жидкости в топливных баках ракет. Поэтому такие задачи имеют актуальные приложения.
Ключевые слова: интегродифференциальное уравнение, параболический тип, первая краевая задача, априорная оценка, обобщенное решение, существование, единственность, непрерывная зависимость.
1. Постановка задачи
Рассматривается первая краевая задача
щ — ¿гу(р(х)Уи) = /(х,г; и), и\р0, и^Т =0 (1-1)
в цилиндре От = Сх[0,Т], где ограниченная область О С Ят, дО £ С2, Ят = дОх х [0, Т] — боковая поверхность цилиндра, у которого верхнее и нижнее основания обозначаются соответственно символами
Бт = {(х,г) : х £ О, г = Т}, Б0 = {(х,г) : х £ О, г = 0},
правая часть уравнения
/(х,г; и) = Г(х,г, и(х,т) ¿хйт),
¿Оь
!© Филатов О.П., 2015
Филатов Олег Павлович (filatov_oleg@samaradom.ru), кафедра уравнений математической физики, Самарский государственный университет, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
функция р £ С1 (С), р(х) ^ р0 = еоивЬ > 0. Измеримая функция Г : С х [0, Т] х х К ^ К удовлетворяет условию Липшица
\Г(х,г,у1) - Г(х,г,у2)\ < Цу1 - у2\
по переменной у при любых допустимых значениях переменных Ь,х,у1,у2 с общей постоянной Ь ^ 0 и интегрально ограничена: \Г(х,г,у)\ ^ Го(х,г) для любых г,х,у, где Го £ Ь2(Ст).
Во всех леммах и теоремах, которые приводятся ниже, предполагается выполнение указанных свойств.
В качестве примера приведем частный случай задачи (1.1), которая связана с измерением уровня несжимаемой жидкости в топливном баке ракеты по уровню жидкости в вертикальной измерительной трубке, которая сообщается с баком в его нижней части и в верхней. Если давление наддува в верхней полости бака обеспечивает постоянство давления жидкости в его нижней части при изменении высоты жидкости в баке Ь(г), то в соответствующей системе координат получим задачу
ди
— - иАи = I (г; и), и\По = ф, и\зт =0
в цилиндре Ст, где
I(г; и) = я(г)(! - щт^, 1(р; и) = 1о - Б J ^ и(х,т)ЗлЗ,т.
1
1(г; и)Г "и Б Jо
Здесь 1(г; и) — высота столба жидкости в измерительной трубке, и = «з = и(х,г) — скорость жидкости в момент времени г £ [0,Т] в точке х = (х1,х2) £ С = = {х £ К2 : х\ + х\ < т2}, Б — площадь поперечного сечения трубки, V — коэффициент кинематической вязкости жидкости. Для вектора скорости жидкости V = («1,«2,Уз) выполняется условие «1 = «2 = 0. Функция д(г) имеет смысл величины вектора ускорения свободного падения, направленного вертикально вниз. В этой модели фигурирует одномерное уравнение Навье — Стокса [1; 2].
Далее используется гильбертово пространство Ш^^т) функций и £ Ь2(Ст), которые имеют обобщенные производные их, £ Ь2 (Ст ),3 = 1, 2,...,т. Скалярное произведение и норма в этом пространстве определяются соотношениями
,, т
(и,«)т = (иV + их, Vх, )йх&, \\и\\т = V7(и,и)т. 3 От j=l
Пространство (Ст) отличается от пространства Ш^^т) тем, что его элементы дополнительно имеют обобщенную производную по г из Ь2(Ст).
Подпространство функций из Ш2^(Ст), след которых на поверхностях Вт и Бт равен 0, обозначим через О0(Ст).
Определение обобщенного решения задачи (1.1) следует схеме, принятой для линейных задач [3; 4]. Функция и £ Ш^1'0(Ст^и^ = 0 называется обобщенным решением задачи, если для любой пробной функции п £ &о(Ст) выполняется равенство
/ игц ¿х ¿г + р(Уи, Уп) ¿х& = ф(х)п(х, 0) ¿х + I(х,г; и)щ ¿х&, (1.2)
и От ^ Оь и О «/ От
где
т
Уп) =53 их, Пх, . j=1
К настоящему времени, начиная с классической книги [5] (первое издание на испанском языке появилось в 1927 г.), опубликовано необозримое количество работ по интегродифференциальным уравнениям. В книге [6] приводятся модели реакции-диффузии, учитывающие эффекты обратных связей, для описания которых используются интегродифференциальные уравнения Вольтерра. Подобные модели рассматривались в [7; 8]. Из последних, наиболее близких работ, отметим [9], где рассматривалась первая краевая задача для параболического интегродифферен-циального уравнения с нелинейностью специального вида. Понятие обобщенного решения в этой работе опиралось на статью [10].
2. Теорема единственности
Зафиксируем число b £ (0,T]. Из [4, теоремы 2, 3, гл. VI, п. 2] следует, что для любой функции u £ W2l'0(Gb) и соответствующей правой части f(x,t;u) существует единственное обобщенное решение w = A(u) £ W21,0(Gb) задачи (1.1), которое удовлетворяет неравенству
Mb < C(M\L2{G) + \\f \\ыоь)1 (2.1)
где постоянная C > 0 не зависит от b £ (0,T], f и ф. Если Uj £ W21,0(Gb) и Wj = = A(uj),j = 1, 2, то функция w = wi — W2 является обобщенным решением задачи (1.1) с правой частью
g(x,t) = f (x,t; ui) — f (x,t; u2) и начальной функцией ф = 0. Из (2.1) получим
\\w\\b < C\\g\\L2{ob). (2.2)
Так как
\g(x,t)\ < L / u(x,T) dr dx , Gt = G x [0,t],
JGt
где u = ui — u2, то с учетом неравенства Коши отсюда следует
\g(x,t)\ < Ly^GlML^Gb),
где \G\— мера Лебега области G. Возведем левую и правую часть этого неравенства в квадрат, а затем проинтегрируем на множестве Gb. После простых преобразований получим оценку нормы
L\ G\ b
\\g\\L2(Gb) \\u\\L2(Gb).
Так как \\u\\L2(Gb) ^ \\u\\b, то отсюда и (2.2) следует неравенство
\\w\\b < q\\u\\b, (2.3)
где
CL\G\b
q=q(b) = ^r.
За счет выбора b > 0 всегда можно обеспечить неравенство
q(b) < 1. (2.4)
Теорема 2.1. Пусть для числа b £ (0, T] выполняется неравенство (2.4). Тогда существует не более одного обобщенного решения задачи (1.1) из пространства W 1'0(Gb).
Доказательство. Пусть п\ и и2 — обобщенные решения. В неравенстве (2.3) и> = и = и\ — и^. Следовательно, (1 — д)||и||ь ^ 0. По условию д < 1, что влечет ||и||ь = 0, а значит, и равенство и\ = и2. Теорема доказана.
Замечание 2.1. Теорема остается справедливой и в случае, если начальный момент времени 0 заменяется на число То € (0,Т), функция
F(x,t, / u(x,T) dxdr)
J о Jg
заменяется на функцию
F(x,t,I + / / u(x,T) dxdT), Jt0J G
где I = const, а цилиндр Gb — на цилиндр G x [То, To + b].
3. Теорема существования
Пусть задано некоторое b £ (0, T], для которого выполняется неравенство (2.4). Построим последовательность функций un £ W1,0(Gb),n = 1, 2,... по следующему правилу. Примем ui =0 и допустим, что функция un уже определена, тогда для правой части f (x,t; un) найдем обобщенное решение un+1 задачи (1.1) из пространства W 1,0(Gb). Последовательность построена. Из (2.3), (2.4) получим
||u3 - u2||b < q\\u2 - ui\\b = q\\u2\\b, ||u4 - u3^b < q||u3 - u2Уь < q2\\u2\\b.
В общем случае, очевидно,
||un+i - un\\b < qn-1\u2\b, n = 2, 3,....
Отсюда следует
hn+lb < ||u2 - ui ||b + ||u3 - u2||b +-----+ ||un+1 - un\\b <
< |u2 - u1||b(1 + q + q2 + ■■■ + qn—1) < .
1-q
Таким образом, получена оценка
HunHb < ^, n =1, 2,.... (3.1)
1 - q
Так как для любого к = 1, 2,... и произвольного n = 1, 2,... выполняется неравенство
qn-1
\\un+k - un\\b < --, (3.2)
1-q
то последовательность u1,u2,... фундаментальна в пространстве W 1,0(Gb). Следовательно, существует предельная функция
u = lim un
п^ж
из пространства W 1,0(Gb). Переходя к пределу при n ^ <ж в равенстве (см. (1.2))
/ un+ir/t dx dt + p(Vun+i, dx dt = J Gb J Gb
= / ф(x)n(x, 0) dx + f (x,t; un)r/t dxdt, Jg J Gb
убеждаемся, что функция и является обобщенным решением задачи (1.1). Следовательно, доказана
Теорема 3.1. Для данной начальной функции ф £ при достаточно ма-
лом Ь > 0 существует обобщенное решение и £ Ш1,0(Оь) задачи (1.1).
Замечание 3.1. Теорема остается справедливой в рамках замечания 2.1. Из неравенства (3.2), после перехода к пределу при к ^ ж, следует оценка
дп-1
\\и - ип\\ъ -, п =1, 2,...
1 - 9
для приближенного решения ип задачи (1.1) и точного и.
4. Глобальная теорема
Сначала докажем лемму о продолжении обобщенного решения. Лемма 4.1. Пусть Ь £ (0,Т) и и1 £ Ш 1,0(Оь) — обобщенное решение задачи (1.1). Тогда существует число 5 > 0 и обобщенное решение и £ Ш1,0(Оь+б) этой же задачи, которое является продолжением решения и1 . Доказательство. Рассмотрим задачу
ди
— - Згу(р(х)Уи) = /(х,г; и1), у\По = ф, у\зъ =0
в цилиндре Оь. Согласно [4, теоремы 2, 3, гл. VI, п. 2] обобщенное решение этой задачи и* существует, единственно и определяется рядом Фурье
= ^ W(1\t)vk(x),
k=i
где vk,к = 1, 2,..., — ортонормированные собственные функции задачи
div(p(x)Vv) = Xv, v\da = 0, х £ G,
отвечающие собственным числам Xk, 0 > Xi ^ X2 ^ ••• ,Xk ^ —œ при к ^ œ. Здесь функция
W(i)k(t) = ФкeXkt + / fk1)(r)eXk(t-T> dT,
J 0
где f (i)(x,t) = f (x, t; ui), ф £ L2(G),fi £ L2 (Gb),
ТО
ф(х) = Фk vk (x), fi(x,t) = ^2 fkl](t)vk (x). k=i k=i
С другой стороны, функция и1 — решение этой же задачи. Следовательно, и* = = и1. Обозначим
ф(2) = и1\^ь, 11 = / / и1(х,т) ¿хЗт. ./о За
Теперь рассмотрим задачу
дг I"г (
— - Злу(рУг) = Е(х,г,Т1 + у у г(х,т) 3,х3,т), г\Вь = ф(2), г^ь+а =0
в цилиндре О х [Ь,Ь + 5], где Бь — нижнее основание, а — боковая сторона
этого цилиндра. По теоремам 2.1, 3.1 и соответствующим замечаниям 2.1, 3.1 при
u
достаточно малом 5 > 0 задача имеет единственное обобщенное решение г = щ, которое совпадает с обобщенным решением задачи
^ - Згу(рУг) = /(2\х,г), г\Въ = ф(2), г\зь,ь+6 = 0,
где
'Ь -¡о
и, следовательно, определяется рядом Фурье
/(2\х,Ь) = Г(х,Ь,1\ + / / щ(х,т) ЗхЗт), ./ь ■! о
тся рядом Фурье
ж
и2(х,г) = ^ (г)ук(х), к=1
где
/г ж
/(2)(т)вХк(г-т) Зт, ф(2)(х) =]Т фк^к(х),
к=1
е(2)(
к = 1
/ (2)(х,г) = ^2 /к2)(*Н (х).
Определим теперь функции
, ч ( и^х,Ь), если Ь € \0,Ъ], х € О, из(x,— Л
и2(х,Ь), если £ € [Ъ,Ъ + 5], х € О, / (3)(х,£) = / (х,Р, и3)
и рассмотрим задачу
ди
— - Зм(рУи) = /{3)(x,t), и\0о = ф, и\яъ+а =0
в цилиндре 0х[0, Ъ+5]. Из [4, теорема 2, 3, гл. VI, п. 2] следует, что ее обобщенное решение и существует, единственно и определяется рядом Фурье
ж
и(х,г) = №к(фк(х), к = 1
где
№ф) = фкеХкг + ( /(3 (т)вХк(г-т) Зт, ио
ж ж
ф(х) = ^2 фк Ук (х), / (3)(х,г) = ^2 /13)№к (х).
(1),
Легко проверяется, что функция №к на отрезке [0, Ъ] совпадает с функцией
(2)
а на отрезке [Ъ,Ъ + 5] — с функцией №к . Следовательно, функция и на множестве О х [0, Ъ] совпадает с функцией и1, а на множестве О х [Ъ,Ъ + 5] — с и2. Это означает, что и = из. Отсюда следует, что функция и является обобщенным решением задачи (1.1) в цилиндре О х [0, Ъ + 5]. Лемма доказана.
Пусть Ъе € (0, Т] — точная верхняя грань тех Ъ > 0, для которых в пространстве № 1,0(ОЬ) обобщенное решение задачи (1.1) существует. По теореме 3.1 выполняется неравенство Ъе > 0.
Лемма 4.2. Обобщенное решение задачи (1.1) существует.
Доказательство. Допустим, что Ье < Т. Тогда по лемме 4.1 решение задачи (1.1) можно продолжить на больший цилиндр, что противоречит определению числа Ье. Лемма доказана.
В следующей теореме постоянная С из неравенства (2.1), функция
к (у) = Су/у(у + 1о), У е я,
где
10 = 4\т\ь2{от )■
Здесь функция Го берется из условия интегральной ограниченности правой части уравнения (1.1).
Теорема 4.1. Если ф е ^(О), то существует единственное обобщенное решение и е Ш1,0(От) задачи (1.1). Если п\ и и2 — решения задачи (1.1), отвечающие начальным функциям ф1 и ф2 соответственно, то выполняется неравенство
\\и1 — и2\\т < к (у) ехр(аТ 2/4), (4.1)
где у = \\ф1 — Ф2\\L2(G), а = (СЬ°\)2.
Доказательство. Существование решения доказано в лемме 4.2. Докажем единственность решения. Пусть й — точная верхняя грань тех Ь е (0, Т], для которых решение задачи (1.1) единственно из пространства Ш1,0(Оь). Покажем, что й = Т. Предположим противное. Тогда ! <Т и существует единственное решение и е Ш 1,0(Ол), при этом начальная функция и\т>Л по теоремам 2.1, 3.1 и лемме 5.1 однозначно определяет продолжение решения и, и оно принадлежит пространству Ш 1,0(О^+г ) при некотором 6 > 0, что противоречит определению числа й.
Докажем неравенство (4.1), используя условие интегральной ограниченности функции Г. Заметим, что функция V = и1 — и2 является обобщенным решением задачи (1.1) с правой частью
д(х,г) = /(х,г; и1) — /(х,г; и2)
и начальной функцией ф = ф1 — ф2. Отсюда
\д(х,г)\2 < Ь2( I \и1 — и2\йх!^2 < Ь2[ I 1 йхйг) \\v\2 = ЬЧ\О\М2.
Проинтегрируем левую и правую части этого неравенства по множеству О1 и
у™^ что ^ \М1, получим
\\д\\< ь2\о\2 0 тм\1 йт. (4.2)
Так как из неравенства (2.1) следует оценка
\\v\t < С\\ф\и2^) + С\\gWL2G)
при любом Ь е [0,Т], то, с учетом (4.2), получим
\П2 < с2(\ф\\^) + 2WФWL2(G)WgWL2(Gt) + \\д\^2^)) < Q(Ф) + с2Ы^2^) <
< Q(ф) + а( т\\vW2 йт,
0
где
Q(ф) = С\WфWL2(G) + То) \\ф\^у
Здесь использовано условие |F(x,t,y)\ < F0(x,t). Обозначим z(t) = IH^. В результате получим интегральное неравенство
z(t) < Q(p) + а ( rz(r) dr.
J о
Отсюда
z(t) < Q(^) exp(at2/2), t € [0, T]
и
Mt ^vQ¥)exp(at2/4), t € [0,T],
что и требовалось. В частности, из неравенства (4.1) следует единственность решения задачи (1.1), но при дополнительном условии интегральной ограниченности правой части уравнения (1.1). Теорема доказана.
Литература
[1] Валландер С.В. Лекции по гидроаэромеханике. СПб.: Изд-во С.-Петерб. ун-та, 2005. 304 с.
[2] Лойцянский Л.Г. Механика жидкости и газа. М.: Наука, 1973. 848 с.
[3] Ладыженская О.А. Краевые задачи математической физики. M: Наука, 1973. 408 с.
[4] Михайлов В.П. Дифференциальные уравнения в частных производных. M: Наука, 1976. 392 с.
[5] Вольтерра В. Теория функционалов, интегральных и интегродифференциальных уравнений. М.: Наука, 1982. 304 с.
[6] Pao C.V. Nonlinear parabolic and elliptic equations. New York; London: Plenum Press, 1992. 777 p.
[7] Нефедов Н.Н., Никитин А.Г., Уразгильдина Т.А. Задача Коши для интегродиффе-ренциального уравнения Вольтерра // Журнал вычислительной математики и математической физики. 2006. Т. 46. № 5. С. 805-812.
[8] Нефедов Н.Н., Никитин А.Г. Начально-краевая задача для нелокального сингулярно возмущенного уравнения "реакция-диффузия" // Журнал вычислительной математики и математической физики. 2012. Т. 52. № 6. С. 1042-1047.
[9] Боголюбов А.Н., Малых М.Д. Об одном классе нелокальных нелинейных уравнений параболического типа // Журнал вычислительной математики и математической физики. 2011. Т. 51. № 6. С. 1056-1063.
[10] Ладыженская О.А., Уральцева Н.Н. Краевая задача для линейных и квазилинейных параболических уравнений. I // Изв. АН ССР. Сер.: Математика. 1962. Т. 26. № 1. С. 5-52.
References
[1] Wallander S.V. Lectures on hydroaeromechanics. Izd-vo S.-Peterb. un-ta, 2005, 304 p. [in Russian].
[2] Loitsiansky L.G. Fluid Mechanics. M., Nauka, 1973, 848 p. [in Russian].
[3] Ladyzhenskaya O.A. Boundary-value problems of mathematical physics. M., Nauka, 1973, 408 p. [in Russian].
[4] Mikhailov V.P. Partial differential equations. M., Nauka, 1976, 392 p. [in Russian].
[5] Volterra V. Functional theory, integral and integro-differential equations. M., Nauka, 1982, 304 p. [in Russian].
[6] Pao C.V. Nonlinear parabolic and elliptic equations. New York-London, Plenum Press, 1992, 777 p.
[7] Nefedov N.N., Nikitin A.G., Urazgil'dina T.A. The Cauchy problem for integro-differential equation of Volterra. Zhurnal vychislitel'noi matematiki i matematicheskoi fiziki [Journal of Computational Mathematics and Mathematical Physics], 2006, Vol. 46, no. 5, pp. 805-812 [in Russian].
[8] Nefedov N.N., Nikitin A.G. Initial-boundary value problem for a nonlocal singularly perturbed reaction-diffusion equation. Zhurnal vychislitel'noi matematiki i matematicheskoi fiziki [Journal of Computational Mathematics and Mathematical Physics], 2012, Vol. 52, no. 6, pp. 1042-1047 [in Russian].
[9] Bogolyubov A.N., Malykh M.D. On a class of nonlocal nonlinear parabolic equations. Zhurnal vychislitel'noi matematiki i matematicheskoi fiziki [Journal of Computational Mathematics and Mathematical Physics], 2011, Vol. 51, no. 6, pp. 1056-1063 [in Russian].
[10] Ladyzhenskaya O.A., Ural'tseva N.N. Boundary value problem for linear and quasi-linear parabolic equations. I. Izv. AN SSR. Ser. Matem. [Proceedings of the Academy of Sciences of the USSR. Series Mathematics], 1962, Vol. 26, no. 1, pp. 5-52 [in Russian].
O.P. Filatov2
GLOBAL THEOREM OF EXISTENCE AND UNIQUENESS
OF THE FIRST BOUNDARY VALUE PROBLEM FOR NONLINEAR INTEGRODIFFERENTIAL EQUATIONS
OF PARABOLIC TYPE
Global theorem of existence and uniqueness of solution of the first boundary value problem for nonlinear integrodifferential equation of parabolic type is proved. If the right-hand side of the equation is integrally bounded, then we have estimate of the norm of the difference of two solutions, which implies continuous dependence of solution on the initial function and uniqueness of solution of the first boundary value problem. The problem under consideration generalizes the real model for measuring the level of incompressible fluid in the fuel tanks missiles. Therefore, such problem have a current application.
Key words: integrodifferential equation, parabolic type, first boundary value problem, generalized solution, existence, uniqueness, continuous dependence.
Статья поступила в редакцию 11/77/2015. The article received 11/77/2015.
2Filatov Oleg Pavlovich (filatov_oleg@samaradom.ru), Department of Mathematics and Mechanics, Samara State University, 1, Acad. Pavlov Street, Samara, 443011, Russian Federation.
