УДК 517.97
Задачи управляемости для модифицированного уравнения переноса в случае переопределения на
выходящем потоке
Исмаил Ахмед Абдел Басет
Кафедра дифференциальных уравнений и математической физики Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198
Доказана локальная разрешимость обратных задач в случае переопределения на выходящем потоке для линейного и нелинейного модифицированного уравнения переноса, где управлением является индикатриса рассеяния.
Ключевые слова: модифицированное уравнение переноса, обратные задачи, дифференциальные операторы, функционал Минковского.
1. Введение
В работе доказана локальная разрешимость обратной задачи для линейного и нелинейного модифицированного уравнения переноса, которую можно трактовать как задачу управляемости. Аналогичный результат линейной задачи для обычного уравнения переноса получен ранее в работе [1]. В основе доказательства лежат известные свойства решений линейных (прямой и обратной) задач, а также двукратное применение теоремы об обратной функции в соответствующих функциональных пространствах (в одном из двух случаев пространство прообразов определяется как множество решений соответствующей задачи). Ранее эта методика в обратной задаче для квазилинейного уравнения теплопроводности была проведена в [2]. Исследованию линейных уравнений переноса посвящено большое количество работ (см., например, [1,3-10] и их библиографии).
2. Предварительные сведения. Обозначения
Введём следующие функциональные пространства и обозначения.
1) Пусть область О С К" строго выпуклая, а ограниченное замкнутое множество V содержится в шаровом слое {у Е К" : 0 < у0 ^ |у| ^ у1 < то}.
2) Ь^ {Б,ЬР(У)) —пространство классов существенно ограниченных функций на Б со значениями в Ь Р(у), 1 < р < то, где Б = О х V х (0, Т).
3) НР(Б) = {и Е ЬР(Б) : щ, (у, V)« Е ЬР(Б), и|г_ Е ЬР(Т-)}, где Г- = 7_ х [0, Т] и 7_ = {(х,у) Е дО х V : (у,пх) < 0}, а пх —внешняя нормаль к границе дО области О в точке х. Это пространство является банаховым относительно нормы
|и\\яр(Л) =
Р II Р,В
+\\щ
|р,в+\\(у, V)«'Р
р,В
+ \\и|г
II р
\\Р,Г-
1/Р
1 < р < то.
и
Статья поступила в редакцию 28 февраля 2008 г.
Работа выполнена по направлению «Функциональные методы в дифференциальных уравнениях и междисциплинарных исследованиях» Инновационной образовательной программы РУДН.
4) Wtp(D) = |е(х,у,г) е Ьр(Б) : ^ е Ьр(В)} с нормой
\\Е\\1(П) +
дЕ
р
ЬР(В)\
1/р
1 < р < ж.
дЬ
5) Нр(О х V) = {< е Ьр(О х V) : (V, У)<р е Ьр(О х V),<11_ е Ьр(^-)}, с нормой
1/р
к„ =
Ьр(ОхУ) р(Ох У) + \М7- \\ьр(7_)
Сх^у — наименьшая константа вложения X в У, т.е. \\-\\у^ Сх^у\\'\\х, С(Х, У) — множество линейных непрерывных операторов из X в У.
3. Обобщённая разрешимость обратной задачи для линейного модифицированного уравнения переноса определения пары функций и и ]
В этом пункте перейдём к рассмотрению коэффициентной обратной задачи щ(х, V, Ь) + (у, V) и(х, V, Ь) + Т(х, V, Ь) и(х, V, Ь) =
= ! J(х,у',Ь,у) а(х,у',Ь)Лу' + Е(х,у,Ь), где (х,у,Ь) е D, (1)
У
и(х, V, Ь) = ц, (х, V, Ь) е 7_ х [0, Т], (2)
и(х, V, 0) = <<, (х, V) е О х V, (3)
где область О = х такая, что каждый х е О можно записать сле-
дующим образом у + VII = х — характеристическая прямая вдоль вектора V, при фиксированном V; причём £ — параметр из [£_, (+], где = (-(у, V), (+ = (+(у, V) : 0 < (у, V) < (+(у^) при всех (у,^)) е х V и (-(у= (+(у^) при (у, V) е дпъ х V, т.е. у + V (± е дО± = {х е дО : их) > 0}.
По условию переопределения на выходящем потоке, точнее на его части,
и|г+ = V(у + v(+,v,Ь), (у + v(+,v,Ь) е Г+ = 7+ х [0,(+ - (-), (4)
т.е. задача одновременного определения функции и и коэффициента J. При этом коэффициент J будем искать в виде
J(х, V, Ь, V) = ](х, у)д(х, V, Ь, V1) + Н(х, V, Ь, V),
где ] — неизвестная, а д и Н — заданные функции. Под обобщённым решением обратной задачи (1)—(4) будем понимать пару функций {и,]}: и е Ир^), ] е Ьр(О х V), почти всюду удовлетворяющих условиям (1)—(4). Пусть данные обратной задачи (1)—(4) удовлетворяют следующим условиям: Т, Тt е Ь^Ф); д, gt е
Ь^; Ь2(у)), | / д(х^, 0,vl)a(x,vl, 0)dv/| > д2 > 0; к,Нг е Ьр(^); е Ьр(D);
У
<<, (V, V)<р е Ьp(ОхV), <р1ч- е Ьр(у-); е Ьр(Г-); е Ьр(Г+). Исследование разрешимости рассмотренной задачи проведём на множестве К = {{и,]} : и е Ир^), J е Ьр(О х V), \\и\\рр>в< к!, \\з\\рр^ху< к2, \Ы\рр>в< кз, У^, V)u\\Pp!D< кз} в предположение, что выполнены условия согласования
<(х,у) = ^х^, 0) при почти всех (х^) е (5)
<(х, V) = V(х, V, 0) при почти всех (х, V) е (6)
Tv- 1diamО < а (постоянные к\,К2,Кз и а ввиду их громоздкости будут выписаны ниже).
Теорема 1. Если выполнены вышеуказанные предположения, то существует единственное обобщённое решение {и,]} обратной задачи (1)-(4), принадлежащее множеству К.
Доказательство. Обратную задачу (1)-(4) сведём к эквивалентной ей в классе НР = НР(Б) х ЬР(О х V) системе интегральных уравнений второго рода
и(у + у£,у,г) = ъ(£,у,у,г) + У
(Ри)(у + ув,у,в - £4 г)+
+](у + ув,у) /9(у + ув, у,в - £ + г, у')а(у + ув,у',в - £ + г)
V
+ ! Н(у + ув,у,в - £ + г, у')а(у + ув, у', в - £ + г) ¿у'
V
¿в, (7)
з(у + у(С+ - г),у) =
19(у + у(С+ - г), у, 0, у')а(у + у(С+ - г), у') ¿у'
1
(у+ус+,у,г)+
-(Ри)(у + ув,у,в - С+ + г)+
Г дЕ
+ (у, V)p(y + у(с+ - г),у) -у — (у + ув,у,в - С+ + г)йв - Е (у + у(С+ - г),у, 0)-
с+
- (Рр)(у + у(с+ - г), у, 0) - I г д
+ ](у + ув, у) — (9(у + ув, у, в - С+ + г, у')а(у + ув, у', в - С+ + г)) ¿у'+
V
Г д
+ дг № + ув, у,в - (+ + г, у')а(у + ув, у', в - С+ + г)) ¿у' ¿в-
V
- /Ну + у(С+ - г),у,0,у')а(у + у(С+ - г),у',0)ёу'|, (8)
где
ъ(^,у,у,г) = р(у + у($ - г),у)+ у е(у + ув,у,в - ^ + г)&в, ?) = £ - г, когда 0 £ - С-, С_ < £ < С+;
£
Ъ(£,у,у,г) = »(у + у(_,у,г - £ + С_) + !Е(у + ув,у,в -£+г)<1в, п = (_,
С-
когдаё-С_ Т, с_ <ё< С+;
(Ри)(у + у£, у, г) = -Т,(у + у£, у, г)и(у + у£, у, г).
Обозначив операторы правых частей уравнений (7), (8) соответственно через А1, Л2, всю систему для краткости можно записать в виде {и,]} = Л{и,3} =
{Л1{и,3 },Л2 {и,3 }}.
Введём в рассмотрение ограниченное замкнутое множество
К = {{и,.П е Ир : \\и\\рр^< к1, \трр>а*у< к2, \ЫГр^< кз, \\(v, V)u\\Pp!D< кз} ,
где
к2 =
т(V)\д\\а\
1
к1 = УМ + УЬЛ+к,
IV\\+(1 + \\Р\\)\М\ + \\Е\Ы т(У )\М\а\\+к
кз = \\<\\ + \\1\\+2\\Е\\+2\\Р\\к1 + \\д\\\\а\\к2 + 2^ЫУ)\М\а\\+к,
\\Р\\= \\Т\\, к — произвольная положительная постоянная. Взяв операторы Л1, Л2 в явном виде, легко убедиться, что множество значений оператора Л лежит в пространстве Ир. Более того, если
Tv0 ^ат О <С1 = к2
\\Е\\+^ т^)\\1г\\\\а\\ + \\Р\\к1 + Л т(У)\\д\\\\а\\к
2
то оператор Л будет непрерывно отображать множество К в себя.
Оценим теперь расстояние между образами двух произвольных элементов {и1,31},{и2,32} из множества К в метрике пространства Ир:
Л1{и1,]1}-А1 {и2,]2} < (Tv0 diamО)С2 {и,]}
р^
Л2 {и1, 31} - Л2 {и2, 32 }
р,ОхУ
< ^Tv01diam О)' С2д2р {и, 3}
1/2
и„
2
р
р
р
Л{и1,31} - Л{и2, 32 }
<
ир
5 (Tv0 1diam О^ + д2 р (Tv0 ^ат О^
1/2
х С2
{и3,3 }
ир
+ ЗС2 (\\и\\рр^ +
х
р
р,ОхУ
где С2 =max{yP\\р, (т^))р/2УgУpУaУp}, и = щ - щ, 3 = 31 - 32. На основании полученных неравенств имеем следующую оценку
А {и 1,31} — А {и2,з'2}11 ^ < 5(Ту0 diam О) + д2р(Ту0 <11ат О) 2
+
+ 6
5 (Tv-1di&m О) + д2-Р (^-^ат О) 2 }О^У{и,3}\\* < (а2 + 6а) С|у{и,^}у
ИЯр
где а = 5(Tv- 1diam О) + д2 'p(Tv22ldiam О)1/2 Таким образом, если
Tv2 1diam О < а =
1 + 20д*р[ у/9 + С-2 - 3 ) - 1
2
р
р
2
р
то вторая степень оператора А будет сжимающей, и, основываясь на принципе сжимающих отображений, можем утверждать, что существует единственный элемент из множества К, являющийся решением системы интегральных уравнений (7), (8), а следовательно, и обобщённым решением обратной задачи (1)-(4). Теорема доказана. □
Замечание 1. По аналогии с теоремой 1 можем доказать, что существует единственное обобщенное решение {и]} из класса Нж(Б) х ЬЖ(О х V).
4. Обратная задача для нелинейного модифицированного уравнения переноса в случае
переопределения на выходящем потоке, где управлением является индикатриса рассеяния 3
Пусть данные обратной задачи (1)-(4) удовлетворяют следующим условиям:
Е,Ег Е Ьж(Б); Е,Ег Е Ьж(Б); а,аг Е Ьж(Б,Ьр(У)); 9,9г Е Ь00{О х У,Ьр(0,Т)),
| / 9(х, у, 0, у')а(х, у', 0)ёу'| > 92 > 0 при (х, у) Е О х V; Ь,Ы Е Ь р(Б); р, (у, V)p Е V
Ьр(О х V); ц Е Шр;(Г_); V Е Шр,(Г+); Ту_1ё1ашО < а и выполнены условия согласования
р(х,у) = ц(х,у, 0) при почти всех (х,у) Е
р(х, у) = V(х, у, 0) при почти всех (х, у) Е 7+.
Теорема 2. Если выполнены вышеуказанные предложения, то исследуемая обратная задача (1)-(4) имеет единственное обобщённое решение {и]} класса Нр(Б) х Ьр(О^).
Поскольку оператор
А : Ьр(О х V) Ш'р(Г+), 3 ^ "
линейный непрерывный и биективный (так как решение единственно), то по теореме Банаха А — изоморфизм, т.е.
А-1 : Ш1р(Г+) Ьр(О х V)
непрерывен и, таким образом, V V 3 ! решение ] линейной обратной задачи, причём
\\J1\Lp(GxV)^ С\\V ^¿(Г^ (9)
при некотором С > 0 (С зависит от Е, Е, а, 9, Ь, ц, О, V, а также от норм дифференциальных и интегральных операторов, введённых в доказательстве теоремы). Пусть ц(х,у,г) = 0, (х,у,г) Е 7_ х [0, (+ - (_); р(х,у) = 0, (х,у) Е О х V;
Ь(х,у,г) = 0 на Б. Положим Ьи = щ+(у, V) и+Е и, Н = {и Е Hp(Б)|3J Е Ш%(Б):
п
и — решение задачи (1)-(4)} с нормой \\и\\я= \\Ьи\\^(о). Тогда ЬН = Ш£(Б) Ь : Н ^ Ш£(Б) — изометрический изоморфизм Н на Ш£(Б) [4, с. 124], Н — полно
и непрерывно вкладывается в Нр(Б). Оттуда же следует, что 3С > 0, Vг1 Е
]0,Т] , Vи Е Н : \\и|(=41 \\ьр{ох^)^ С\\и\\н. Отметим также, что по теореме о
(10)
следах для Нр(Б) оператор
Ли : Н ^ Ьр (О х V),
и ^ и^г=11
непрерывен в силу непрерывности вложения И в Ир Рассмотрим теперь нелинейную обратную задачу
Ut(x, v, t) + (v, V)u(x, v, t) + v, t)u(x, v, t) + [S(u)] (x, v, t) =
= J J(x,v,t,v')a(x,v',t)dv', где (x,v,t) E D,
V
где нелинейная часть
i
[S(u)](x,v,t) = J J Q(x,v,x',v',t)a (u(x' ,v',t'))dx' dv' dt' (11)
0 GxV
с условиями
u(x,v, 0) = 0, (x,v) E G x V, (12)
u(x,v,t) = 0, (x,v,t) E Y- x [0,Z+ - Z-), (13)
u|r+ (x,v,t) = v (y + v(+ ,v,t), (y + v(+,v,t) E y+ x [0, Z+ - Z-), (14)
J(x,v,t)= j(x,v) g(x,v,t). (15)
Будем для краткости обозначать эту задачу обратной задачей II'. Положим x(x, v) = j(x, v) g(y + vZ+ ,v, t).
Определение 1. Пусть X,Y — нормированные пространства. Отображение f : X ^ Y называется строго дифференцируемым в точке xo, если
f (x + h) - f (x) = f'(xo)h + R(x, h),
где f'(xo) E L(X,Y), а \\R(x,h)\\/\\h\\^ 0 при x ^ xo, h ^ 0. Напомним теорему об обратной функции.
Теорема 3. Пусть X, Y —банаховы пространства, U —открытое множество в X, а отображение f : U ^ Y строго дифференцируемо на U и f'(xo) : X ^ Y —изоморфизм X на Y для некоторой точки x0 E U. Тогда существует такая окрестность U' точки xo, что f осуществляет гомеоморфизм U' на открытое множество f (U'), f '(x) — изоморфизм X на Y при x E U', f-1 : f (U') ^ X строго дифференцируемо на f (U') и (f-1)' f (x)) = [f '(x)}-1 при x E U' (т. е. f — диффеоморфизм U' на f (U') класса C1).
4.1. Формулировка основного результата
Всюду в дальнейшем будут использованы принятые ранее обозначения. Пусть а : rot ^ rot дважды непрерывно дифференцируемо на rot, причём
|a(u)| < Cilul,
la'(u)l < Ci, а'(0) = 0, (16)
la''(u)l < C2
(примерами таких функций могут быть a(u) = \/1 + u2 — 1, a(u) =
u2
vi , u2
,3
Vl+u2
a(u) = —2 и многие другие).
Далее, Q(x, V, х', V,Ь) непрерывна на О х V х О х V х [0, Т] и 2 ^ р < ж. Введём оператор
т
Б : ИР(Р) м Wtp(D), [Б(и)](х^,г) = ! [ Q(x,v,x'а(и(х',£)) ¿х'^'М'.
0 ОхУ
Теорема 4. Пусть £,^ е Ь^(Р), \^(у + v(+,v,t)\ > £о > 0, (у + v(+,v,t) е Г+; е ЬХ(Р); аа е Ь(Р х V); д,дь е Ь(О х У,Ьр(0,Т)), \д(у + v(+,v,t)\ > до > 0; <р, ^, V)? е Ьр(О х V), » е Wtp(Г-); V е Wtp(Г+), Г± = 7± х [0, С+ — С-); Tv-хё1аш О < 0, и выполнено условие V(у + vZ+ ^^) = 0 при (х^) е Тогда .задача II', где у —искомая, а д —априори заданная функция, имеет единственное решение (и, у) в окрестности нуля в Н(Р) х Ьр(О х V) при достаточно малых по норме V е Нр(О х V).
Замечание 2. Доказательство теоремы 4 сводится к проверке выполнимости условия теоремы об обратной функции, при этом главными моментами являются проверка строгой дифференцируемости соответствующих отображений и выбор функциональных пространств, связанных с задачей II'.
4.2. Строгая дифференцируемость оператора Б Предложение 1. Оператор
т
Б : Нр(Р) м Wtp(D), [Б(и)](х^^) = ! [ Q(x,v,x'а(и(х'^',£)) Ах'ММ
0 Ох У
строго дифференцируем на Нр(Р) (доказательство изложено подробно в работе [11]), где
1а : ИР(Б) ^ Wf(Б),
и(х,у,~Ь) ^ а(и(х,у,~Ь)^, (17)
: Wf(D) ^ Wp(Б), т
и,(х,ю,£) ^ J J Я(х,у,х',у'^)и,(х'дх'&и'dí/, (18)
0 ОхУ
Б(и) = [1д ◦ 1а](и). (19)
4.3. Завершение доказательства основного результата
По теореме об обратной функции 3 существуют такие открытые окрестности нуля и и VI в Н и Wt(D) соответственно, что оператор £ : и м Ьи + Б (и) осуществляет диффеоморфизм и на V1 класса С1. При этом оператор ц = £-1 : Ух м и строго дифференцируем на VI как отображение в Н и, следовательно,
в Нр(Р) и ц'(,1) = [£'(ц(,1))] при J е VI. В этом случае, в силу линейности и непрерывности оператора х(х^) м х(х^)д(х^^)/д(у + vZ+,v,t) из Ьр(О х V) в W£(D) следует, что оператор Р : Ьр(О х V) м Н(Р), х(х^) м Р(х) = П (х(х^)д(х^^)/д(у + vZ+,v,t)) строго дифференцируем по Фреше в окрестности нуля VI из Ь р (О х V), а операторы
01 : У2 - ЬР(С х V),
X - (РЫ), ,
г+
02 : V2 - ЬР(С х V),
X - (V, V)(P(х))
03 : V2 - ЬР(С х V),
X - Б(Р(х))
(20) (21) (22)
также строго дифференцируемы на У2. Далее будем искать решение нелинейной задачи II' в виде и = Р(х), где и — решение нелинейной задачи II' с 3(х,ь^) = Х(х,у)д(х,у,г)/д(у + у(+,у,г). Введём оператор
г
+
+
М (х) =
X а ^у'
- (Р(X)), - (V, V)(P(X)) - ^ (Р(X))
М/(0)х = |У X а - (Р'(О^), Ц - (V, V) (Р'(0)х) |
1
Е(у + Ч+^'О
(23)
(24)
где Р'(0)х = п'(0) (х(х,у)д(х,у,£)/д(у + т.е. решение прямой линейной
задачи с 3(х,у,{) = х(х,у)д(х,у,£)/д(у + (решение существует, как по-
казано в [12], и оператор ц'(0) непрерывен).
Подставив в (11) и1х=у+ъс++ из (14) и 3(х,у,£) = х(х,у)д(х,у,£)/д(у + приходим к уравнению
М (х) = у(у + (25)
к которому в силу непрерывности М'(0), а также существования и ограниченности [М'(0)]-1 (это следует из корректной разрешимости линейной обратной задачи), можно применить теорему 3. Согласно теореме 3 существуют такие открытые окрестности нуля и* в Ьр(О х У) и У* в что М : х ^ V осуществляет
диффеоморфизм и* на У* класса С1. Поэтому V V £ У* 3! х £ и* : М(х) = V.
Теперь, сопоставляя определение оператора М (см. (23) и исходное уравнение, на основе полученного результата заключаем, что при и = Р(х) пара (и,]) — единственное решение задачи II'. Заметим,что принято считать решением пару (и]), тогда как решением является ], поскольку и однозначно определяется через ]. Теорема доказана.
5. Заключение
X
г
+
+
+
г
+
X
X
Аналогичные результаты можно получить для случая, когда управлением является множитель в правой части (т. е. в функции источников) с тем же условием переопределения.
Литература
1. Волков Н. П. Обратные задачи для математических моделей физических процессов. — М.: МИФИ, 1991. — С. 16.
2. Сухинин М. Ф., Акпата Э. О задаче управляемости для квазилинейного уравнения теплопроводности // Вестник РУДН. серия Математика. — № 3(1). — 1996. — С. 119.
3. Орловский Д. Г. Решение одной обратной задачи для уравнения переноса с интегральным переопределением. — М.: МИФИ, 1991. — 71 с.
4. Прилепко А. И., Волков Н. П. Обратные задачи определения параметров нестационарного уравнения переноса по переопределениям интегрального типа // Диф. уравнения. — Т. 23, № 1. — 1987. — С. 124-136.
5. Kaper H. G., Lekkerkerker G. G., Heitmanek J. Spectral in Linear Transport Theory. — Basel, 1982.
6. Greiner G. Spectral Properties and Asymptotic Behavior of the Linear Transport Equation // Math. Z. — Vol. 185, No 2. — 1984. — Pp. 167-177.
7. Voigt J. Spectral Properties of the Neutron Transport Equation // J. Math. Anal. Appl. — Vol. 106, No 1. — 1985. — Pp. 140-153.
8. Волков Н. П. Анализ математических моделей физических процессов. — М.: Энергоатомиздат, 1988. — 11-19 с.
9. Волков Н. П. Теоретико-функциональные методы в задачах математической физики. — М.: Энергоатомиздат, 1986. — 22-26 с.
10. Тихонов И. В. Корректность обратной задачи с финальным переопределением для нестационарного уравнения переноса // Вестн. Моск. Ун-та. Сер. 15, Вычисл. Матем. и Киберн. — № 1. — 1995. — С. 56.
11. Хамди Н. Обратная задача для нелинейного уравнения переноса с финальным переопределением. — Рус. — Деп. в ВИНИТИ, 22.07.2003, № 1435-В2003. — 2003.
12. Волков Н. П. Обратные задачи для нестационарного кинетического уравнения переноса с разрывными переменными: Дис. канд физ.-мат. наук / МГУ. — М., 1986.
13. Сухинин М. Ф. Избранные главы нелинейного анализа. — М.: РУДН, 1992.
14. Канторович Л. В., Акилов Г. П. Функциональный анализ. — М.: Наука, 1977.
15. Прилепко А. И, Волков Н. П. // Диф. уравнения. — Т. 24, № 1. — 1988. — С. 136-146.
16. Владимиров В. С. // Труды Математического института им. В.А. Стеклова АН СССР. — № 61. — 1961.
17. Ладыженская О. А. Краевые задачи математической физики. — М.: Наука, 1973.
18. Сухинин М. Ф. // Теоретическая и математическая физика. — Т. 103, № 1. — 1995. — С. 23.
19. Кириллов А. А., Гвишиани А. Д. Функциональный анализ. — М.: Наука, 1988.
20. Красносельский М. А. Интегральные операторы в пространствах суммируемых функций. — М.: Наука, 1966.
UDC 517.97
Problems of Controllability for Modified Nonlinear Transport Equation in Case of Redetermination in Output Flow
Ismail Ahmed Abdel Baset
Department of Differential Equations and Mathematical Physics Peoples' Friendship University of Russia Miklukho-Maklaya str., 6, Moscow, Russia, 117198
The theorem on local controllability for a linear and nonlinear system described by the modified transport equation in the case of redetermination in output flow is proved.