CC BY
10
ВЕСТНИК ЮГОРСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
_2018 г. Выпуск 3 (50). С. 57-66_
DOI: 10.17816/byusu2018057-66 УДК 517.956
В. В. Ротко
ОБРАТНЫЕ ЗАДАЧИ ДЛЯ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ КОНВЕКЦИИ-ДИФФУЗИИ С ТОЧЕЧНЫМИ УСЛОВИЯМИ ПЕРЕОПРЕДЕЛЕНИЯ
Работа поддержана грантом РФФИ и Правительством ХМАО-Югры № 18-41-860003 р_урал_а.
В работе рассматривается вопрос о корректности в пространствах Соболева обратной задачи об определении функции источников в квазилинейной параболической системе второго порядка. Проблемы подобного вида возникают при описании процессов тепломас-сопереноса, диффузионных процессов, процессов фильтрации и во многих других областях. Главная часть оператора линейна. Неизвестные функции, зависящие от времени, входят в нелинейную правую часть. В том числе, в этот класс задач входят и коэффициентные обратные задачи об определении младших коэффициентов в параболическом уравнении или системе. В качестве условий переопределения рассматриваются значения решения в некотором наборе внутренних точек. В качестве краевых условий берутся условия Дирихле или условия задачи с косой производной. Задача рассматривается в ограниченной области с гладкой границей. Однако результаты допускают обобщения и на случай неограниченных областей таких, в которых соответствующие теоремы о разрешимости прямой задачи имеют место. Приведены условия, гарантирующие локальную по времени корректность задачи в классах Соболева. Условия на данные задачи минимальны. Полученные результаты являются точными. Задача сводится к операторному уравнению, существование решения которого доказывается при помощи априорных оценок и теоремы о неподвижной точке. Полученное решение обладает всеми обобщенными производными, входящими в уравнение и принадлежащими пространству Lp с р > п + 2, и обладает необходимой дополнительной гладкостью в некоторой окрестности точек переопределения.
Ключевые слова: параболическое уравнение, обратная задача, тепломассоперенос, краевая задача, функция источников.
V. V. Rotko
INVERSE PROBLEMS FOR MATHEMATICAL MODELS WITH THE POINTWISE OVERDETERMINATION
In the article we examine well-posedness questions in the Sobolev spaces of an inverse source problem in the case of a quasilinear parabolic system of the second order. These problem arise when describing heat and mass transfer, diffusion, filtration, and in many other fields. The main part of the operator is linear. The unknowns occur in the nonlinear right-hand side. In particular, this class of problems includes the coefficient inverse problems on determinations of the lower order coefficients in a parabolic equation or a system. The overdetermination conditions are the values of a solution at some collection ofpoints lying inside the spacial domain. The Dirichlet and oblique derivative problems under consideration. The problems are studied in a bounded domain with smooth boundary. However, the results can be generalized to the case of unbounded domains as wellfor which the corresponding solvability theorems hold. The conditions ensuring local (in time) well-posedness of the problem in the Sobolev classes are exposed. The conditions on the data are minimal. The results are sharp. The problem is reduced to an operator equation whose solvability is proven with the use of a priori bounds and the fixed point theorem. A solution
possesses all generalize derivatives occurring in the system which belong to the space Lp with p > n + 2 and some additional necessary smoothness in some neighborhood about the overdetermination points.
Keywords: parabolic equation, inverse problem, heat-and-mass transfer, boundary value problem, source function.
Введение
Мы рассматриваем вопрос об определении вместе с решением функции источника и коэффициентов, характеризующих параметры среды в квазилинейных математических моделях конвекции-диффузии. Пусть в - область в Еп с границей Г класса С2 и (( = (0,Т) X С. Соответствующая параболическая система имеет вид:
Ьи = щ + А( Ь,х,0)и = й (х,Ь,и^и) 1 [ ( Ь,х)чЖЬ),(Ь,х) е (, (1)
где А - матричный эллиптический оператор вида А(Ь.х.Б)и = — у=1ац (Ь,х)их.х.+
1 а I их. + а 0 и, а а I - Их И матрицы, и - вектор длины Л и Ч(Ь) = ( Ч!.( Ь),Ч2(0, ■ ■ -,Чт(0) - неизвестные функции, подлежащие определению вместе с решением , которые входят как в правую часть, так и в сам оператор как коэффициенты. Структуру оператора А мы опишем ниже. Система (1) дополняется начальными и граничными условиями:
и \ с=о=и0,Ви \ 5 = д,Б = ( 0,Т)хГ, (2)
где В и = ]_У1 (Ь,х)их + а( Ь,х)и или В и = и, щ - ¿-я координата внешней единичной нормали к и ( ) ( ) матрицы-функции размерности , принадлежащие
классу (Б). Условия переопределения записываются в виде:
и(х[^) = р[(= 1,2, ■.,,5,т = Ия. (3)
Таким образом, дополнительные условия являются данными замеров решения (например, концентрации переносимого вещества) в определенных точках области. Обратная задача состоит в нахождении решения и уравнения (1) и функций ч^Ь) (I = 1,2, ■.., т) по данным (2), (3).
Проблемы подобного вида возникают при описании процессов тепломассопереноса, диффузионных процессов, процессов фильтрации и во многих других областях (см. [1]). Большое количество результатов было получено в случае линейной по своим аргументам функции . Можно отметить работу [2], где получена теорема существования и единственности решений задачи (1)-(3) в пространствах Гельдера в случае И = 1 , г = 1. В случае п = 1, г = 1 и С = Е некоторые результаты получены в работах [3, 4]. В работах [5, 6] была рассмотрена также и задача об определении младшего коэффициента в параболическом уравнении. Общие теоремы о разрешимости абстрактных задач такого вида в квазилинейном случае получены в монографии [7]. Однако стоит отметить, что результаты [7] применимы к задаче (1)-(3) в случае некоторых специальных функций д и в случае краевых условий, таких, что область определения оператора Ь не зависит от времени. Однако даже в этом случае в [7] используются излишним условиям гладкости и согласования на данные. Ряд задач, входящих в класс (1), был рассмотрен в работе [8]. Численные методы решения различных модельных задач, входящих в класс (1)-(3), рассматривались, например, в работах [9, 10] и многих других. В частности, в работе [10] численно решалась обратная задача (1)-(3) об определении функции источника для квазилинейной системы параболических уравнений. Наши результаты обобщают результаты работы [14] (см. также [15]). В данной работе мы получим локальную теорему о разрешимости задачи.
Опишем содержание работы. В первом параграфе описаны условия на данные задачи и сформулированы основные результаты. Во втором параграфе приведено их доказательство. Обозначения функциональных пространств стандартные (см., например, [12]).
Обратные задачи для математических моделей конвекции-диффузии
с точечными условиями переопределения Определения, обозначения и вспомогательные результаты
Пусть Е - банахово пространство. Символом Lр( G; £") ( G - область в Мп) обозначаем пространство сильно измеримых функций, определенных на G, со значениями в Е, наделенное нормой | | | | и(х) | | Е| | L (С) [12]. Мы также используем пространства Гельдера Са( G). Обозначения Wp( G;E), Wp( Q ;Е) пространств Соболева являются стандартными (см. определения в [12, 13]). Если Е = С (Е = М) или Е = С п (Е = Мп), тогда последнее пространство обозначается через Wp(Q ). Аналогично используем обозначения Wp (G) или Са( G ) вместо Wf( G^) или Са( G^). Таким образом, включение u G Wf( G) (или u G Са( G)) для данной вектор-функции ( ) означает, что каждая из ее компонент принадлежит
Wps( G) (или Са( G)). В этом случае норма вектора есть просто сумма норм координат. То же самое соглашение принимаем для матриц-функций. Для интервала / = ( О, Г) положим Wp5, r( Q) = Wpr(/ ; Lp( G )) n L р(/; WpS( G)). Соответственно, r(S) = Wpr(/ ; L p(Г )) n
L p(/ ; WS( Г)). Аналогично определяем пространство Гельдера Сг' s( Q ).
Пусть ßs(Xj) - шар радиуса 5 с центром в точке Xj. Далее мы используем следующие обозначения: QT = ( О, т) X G, S7 = ( 0,7) X Г. Дан набор точек {ху} из (3), параметр 5 > О
назовем допустимым, если ( ) , ( ) ( ) для , . Пусть
G5 = и j Ва(х), Q 5 = ( О, Г) X G5, Q| = ( 0,т) X G5.
Рассмотрим задачу (1), (2) и сформулируем один вспомогательный результат. Мы будем предполагать, что у нас выполнены условия. Условия на коэффициенты:
а и G C(Q),afc G L p( Q ), 7, ö G С 1/2 2(S),p > n + 2 ; (4) а ¿у G L oo ( О, Г ; WpS( G5)), afc G Lp(О, Г ; WpS( G5)), i ,7 = 1,2.....n, k = 0, 1.....n, (5)
для некоторого допустимого 5 > О и s G (""+", 1 ]. Опишем условия параболичности оператора
L. Рассмотрим матрицу А0( t,x, f) = — £Пу=1 а ¿у ( t,x)f£fy ( f G Мп) и предположим, что найдется постоянная , такая, что корни полинома ( ( ) ) ( - единичная матрица) удовлетворяют условию
Re р< — 51 | f | 2 VfGMn V(x,t)GQ. (6)
Пусть В0u = и в случае условий Дирихле в (2) и В0и = £П= 17у d^.u в противном случае.
Условие Лопатинского может быть записано в виде: для любой точки ( ) , и операторов ( ) и ( ), записанных в локальной системе координат в этой точке (ось направлена по нормали к S и оси ,...,уп _ 1 лежат в касательной плоскости в точке (х0, t0 )), система
(ЯЕ + А 0(X0 , t0 , if ' , öyn))v(z) = О, В 0(х0 , ^ , if ', 5Уг>( О ) = Лу, (7)
где f ' = ( fi,. .. ,fn _ 1 ), yn G М+, имеет единственное решение из с( М ), ограниченное на бесконечности при всех , , и таких, что . Мы
также предполагаем, что
2--
U0(x)GWp р( G),^GWp2fco' k0(S), (8)
( ) ( ) ( )
где в случае условия Дирихле и , в противном случае
2+s-—
u 0(х) G Wp p( G5) д л я h e к о то р о г о д о пу ст и м о г о 5 > О и (9)
s G ((n + 2 )/р, 1] .
Даны постоянные 81 < 82 < 8. Построим вспомогательную функцию р Е С^ (, такую, что р = 1 в области и р = 0 в Сд2.
Теорема 1. Пусть выполнены условия (4)-(9) для некоторого допустимого 8 > 0 и 5 Е ((п + 2 )/р, 1] , / Е Ьр( ( т), / р Е Ьр( 0,т; Шр( в3)) и т Е ( 0, Т] . Тогда существует единственное решение и Е Шр'( т) задачи:
Ьи = и£ + Аи = / ((х, Ь) Е ( ),и\£=0 = и0(х),Ви\5 = д. (10)
Причем ри£ Е Ьр(0,т;Шр(С8)), ри Е Ьр(0,т;Шр +5( С8)). Если д = 0, и0 = 0, то справедлива оценка:
I I иI I + I I ри£I I ьр(0,х-п*(с8)) + | I ри| I ьр(о+*(в 8)) ^
С [ | | /11 1р,(д*) + 11 р! | | Ьр(0,г;\¥$(С6))], (11)
где постоянная с не зависит от /, решения и и т Е (0, Т].
Доказательство. Доказательство этой теоремы может быть найдено в работе [17, Теоремы 1, 2]. Основное утверждение теоремы известно, см., например, [13]. Дополнительная гладкость по существу вытекает из известных результатов о внутренней гладкости решений параболических и эллиптических задач.
Основные результаты
Рассматривается уравнение
г
и £ + А и = /(х, Ь, и, Уи) + (х, Ь) Ч £( Ь) + !0(х, Ь). (12)
¿=1
Мы рассматриваем задачу (1)-(3) о восстановлении правой части уравнения вида / = ^= 1 й(х,Ь)ц£(Ь) + /0(х,Ь) и коэффициентов, в частности входящих в главную часть уравнения (1 ). Предположим, что оператор имеет вид:
А
= Ь0~ ^ Чк (Ь)Ьk,
к=г+1
п п
Ьки = — ^ а^ (Ь,х)их.х. + ^ ак (Ь,х)их. + а¡к(Ь,х)и,
¿,У=1 ¿ = 1
где Поскольку неизвестные могут встречаться в главной части урав-
нения, мы будем искать их в классе ( ). Построим матрицу ( ) размерности
тХт, чьи строки с числами с (] — 1 )!г + 1 по }Л занимают вектор-столбцы: ( ( ) ( ) ( ) ( ) ( ))
Мы предполагаем, что
р 1 Е С[0,Т]),и= р¡(0)Ц = 1,2.....б,),
_ (13)
у£,а Е С1 /г,2 (5),р > п + 2,
а £ Е С (() П Ь т (0, Т; Шр5( С5)), а к Е Ьр( ( ) П Ь р ( 0, Т; Шр5( С5)) ( )
/ Е Ьр(()ПЬр(0,Т;Шр5(С5)) (I = 1.....т), (15)
Обратные задачи для математических моделей конвекции-диффузии
с точечными условиями переопределения для некоторого допустимого 8 > 0, 5 > ( 2 + п)/р , и к = 0 , г + 1 ,. . ., т;
а [ 0,Г] ) (16)
для всевозможных значений /,_/', к, /, а. Также необходимо выполнение следующего условия: Условие Существует число 80 > 0, такое, что
| с1 е t В( 0 | > 8 0 н а(0 ,Г). Заметим, что элементы матрицы В принадлежат С( [ 0 , Г] ). Рассмотрим систему
V/с( 0 ) + Iо" о(ху, 0 ) — / ( 0, ху, и о(ху), Рио(ху)) =
г т
= ^Ц о д/д(ху,0)+ ^ Цо ^ Iди о (ху) /с = 1 /с=г+1
по отношению к вектору ¿7 о = ( Цо х , Цо 2,..., Ц от1). При выполнении условия (А) система имеет единственное решение. Обозначим Л о = I о — г+х Ц о д I Пусть Вк - шар радиуса с центром в нуле в Е(п+1 ) п.
Условие (В). Функция /(х, и, р ) непрерывна по совокупности переменных (и, р £ ( ) ), для любого найдется постоянная такая, что
| /(Г,х,и 1 ,р 1) —/(*;,х,и2 ,р 2 ) | < Мо( | "1 — и2 | + | Р1 — Р2 | ), для всех (и 1 ,р 1 ), (и2 ,р 2 ) £ВК, /(Г,х,и,р ) £ С([0, Г] X В^; С5)).
Условие (C). Функция /(х, и,р ) дифференцируема по параметрам (и,р ) £ Е(п+1 при
п.в. ( Г,х)£<? 5 и _
/м(Г,х,и,р ),/р.(Г,х,и,р ) £ С([0,Г] X Вк;Жр5( С5)) V/? > 0,
для любого найдутся функции Ф 1(^х),Ф 2(^х) £ с([0,Г];1 р(С)), такие, что
| | /м(х, и 1 ,р ^ — /и(х, и2 ,р 2 ) | | £(мл) <
< | Ф 1( ^х) | | и 1 и2 | + | Ф 2 ( £,х) | | Р 1 — Р 2 | ,
1 1 /р ,( ^х,и ^р ^ —/р ,( ^х,и2 ^ 2 ) 1 1 <
< | Ф 1( ^х) | | и 1 и 2 | + | Ф 2 ( £,х) | | Р 1 — Р 2 |
для всех ( ) ( ) . Здесь величины ( ) - соответствующие
матрицы Якоби и | | • | | цд^ - норма в пространстве линейных непрерывных отображений из в
При указанных условиях теорема существования примет следующий вид. Теорема 2. Пусть условия (А)-(^, (13)-(16) выполнены. Предположим также, что оператор Мо = д с + Л о параболический и выполнено условие Лопатинского, т. е. условия (6), (7) выполнены. Тогда найдется число ( такое, что на промежутке существует
единственное решение (и,ц 1 ,ц2 , ...,цш) задачи (1)-(3), такое, что и £ 1р( 0,т о;Ир2( С)), и с £ I р( @т°), Ц ¿( 0 £ С( [0,т о ] ), / = 1,. ..,т. Кроме того, <р и £ 1р( 0,т о;И^2+5( С)), <р и с £ ( ( )). Доказательство. Найдем решение задачи.
г
Ф с + ЛоФ = ^ Ц о ¿/ ((х,0 £ (?),Ф | с=о = и о(х), ВФ | 5 = д. (17)
к=1
По теореме 1 Ф £ Жр2, ^ <?), <рФс £ I р( 0,Г; Жр5( С)), <рФ £ I р( 0,Г; Жр2+5( С)). Из теоремы III 4.10.2 в [11] вытекает, что <р Ф £ С([ 0,Г]; Жр2+5 -2/р( С )) с С([ 0,Г]; С2+5 -2 /р -п/р(С)). Сле-
довательно, <р Ф £ С( [О,Г] ;С2( С )). Из уравнения (17) вытекает, что Ф [(хуД) £С([ О,Г] ). Сделаем замену переменных (Ц = (Ц о+ 9 1 и и = V + Ф в (1). Получим новую задачу:
ЬV = V[ + Ло^- ^ ц 1 /Ь /V =
/с=г+1
/ //I
(18)
1=1 1=Г+1
1?|£=0 = 0,В1?|5 = О,
1(ху,0 = $( О = !/>;( О - Ф(ху,г) £ С\[ О, Г] ),#(О) = ф[( О) = О . (19)
Стандартные теоремы о разрешимости параболических задач гарантируют локальную разрешимость задачи (18) при фиксированном векторе ц ^^. Однако у нас используются специальные классы данных и решений. Поэтому мы наметим доказательство. Определим пространство Яо , т = {и £ И^1,2( @ т): <р и £ Ь р( О,т ; С)), <ри £ Ь р( О ,т ; И^2+5( С))} . В качестве нормы в этом пространстве возьмем величину
I I и I I Яот = I I и I I м/р1 , 2(<гт) + I I и I I Ьр(0,т; И$(С)) + I I <Ри I I 1р(0,т; Мр2 +5(С)).
Используя теорему 1, можем сказать, что решение и £ Н0, т задачи (18) есть решение операторного уравнения:
V = М0" 1 /( г, х, V + Ф , 171 + 17Ф ) + т т (20)
^ МО/сЬл V + МОг,х)ц^(г),
+
/с=г+1 к=1
где М0 и = и [ + Л 0и, /д. = // при к < г и /к = Ь кФ при к > г. Зафиксируем параметр ( ( ( ) ))) введем обозначения
( ) ( ( )) ( 0, т; И/р(С )),
| | V | | 0,т = | | V | | с( [ 0, т; С ЧС)) + | | V | | С( [ 0 , т] ; С 1+51(С)).
Теорема 1, условия на функции / и неравенство Гельдера гарантируют оценку
т
I I МО ( г,х)(г) - ц2(г)) I I Яо,т < ^ I I Ц 1 - Ц2 I I с([0,т]), (21)
к=1
где с0 - постоянная, не зависящая от параметра т. Пусть | | ЦЦ | | С([0,т] ) < г0. Аналогично, используя условия на коэффициенты, получим:
т
I I МО 1 ^ Ц 1 к Ь/VI I Яо,т<С1Г0 | | V I I Яод. (22)
/с=г+1
Уравнение (20) можно переписать в виде:
т
V = Й( V) + /0 ,/0 = МО 1 ^ /к ( г,х)ц 1 к( г). (23)
/с=1
Покажем, что если параметр достаточно мал, то на малом промежутке времени уравнение (23) разрешимо в классе Я0, т. Фиксируем параметр г0 < 1 / 4 с 1 и шар Вг° = { ц £
Обратные задачи для математических моделей конвекции-диффузии
с точечными условиями переопределения С([0,т]): | | ц| | С([0,т] )< г0 }. Фиксируем также шар ВЯо = {V е : | | V| | Яот < И0} с И0 = 2 ( | | И( 0 ) | | Яот + сх г0 ). Получим оценки. Пусть е 5Ко (/ = 1, 2 ). Имеем, используя теорему 1, условия на f и (22), что
| | ) - И(у2) | | Яо,т < с3 | | V! - | | 0,т + сх Г^ | | - | | ^ . (24)
При оценке нормы | | f( г, х, V 1 + Ф, V+ 17Ф ) — Я ^х, + Ф, V+ 17Ф ) | | Уот используем неравенства:
| | f ( г, х, V! + Ф, V V! + V Ф ) — Я г, х, + Ф , V + V Ф ) | | ,( 0, Т) <
<С4 | | 1?2 | | с([0,Г] ; с 1(ё)),
вытекающее из условия (B), и неравенство
| | <р(Я г, XV! + Ф , V V ! + VФ ) — я г,х^2 + Ф, V V2 + VФ )) | | ^ 0т; ^^
< С51| -Р2 II о,-т>
(25)
(26)
где постоянные с ^ не зависят от т. Чтобы получить последнее неравенство, используем представление:
( ( ) ( ) ( ( ) ( ))) ( ( ) ( ) ( ( ) ( )))
I ( г, х, V !( г, х + Л) + тД V, V(V !( г, х + Л) + тД v( г, х + Л))) • д v( г, х + Л) + 'о
п
^ /р. ( г, х, V ! + тД V, V( V ! + тД V)) • д г, х + Л) йт, Д V = v2 — V ! , 1=1
f( г,х^ 2 (г,х) + Ф( г,х)^( V 2( г,х) + Ф ( г,х))) —Я г, х, V!(г,х) + Ф( г,х), V( V !( г,х) + Ф( г,х)))
= I /и ( г, х,
п
V !( г, х) + тД v( г,х), V(V ! + тД v)( г, х)) • Д V
+ ^ /р . ( г, х, V ! + тД V, V(V ! + тД v))Д V;,;. йт 1=1
и определение нормы в пространстве /1р( С ). Оценим правые части в (25), (26). Теорема III
2_2 /т) 1 2
4.10.2 в [11] гарантирует, что V е С( [0,т] ;/ 'Р(С)) (класс / ' (@т) вложен в С( [ 0, т] ;/2 _2/р( С))) и что <р V е С( [ 0, т] ;/2 _ 2/р+5( С)). Можем считать, что постоянные вложения в обоих случаях не зависят от т (поскольку четное продолжение (относительно точки т) на промежуток [0, 2 т] и далее нулем сохраняет норму и класс Н0 т для функций V, которые обращаются в ноль при г = 0). Используя теоремы вложения (см. теорему 4.6.2 в [12]) и интерполяционные неравенства, а также формулу Ньютона-Лейбница, имеем:
| | 17 | | С([0,т];СЧС)) < с6 | | 17 | | с([ 0 , т]; <+п/р(С)) << С? | | 17 | | с( [0,т] ; Ж2 -2 /р(С)) | | 17| | С([0,т];!р(С)) <
С( [ 0, т] ; ' " "(С))--"'""с( [0,т] ; Жр2 - 2 / Р(С))'
< С» | | " И С([0,т];<-2/"(<Ч)1И11ит)т(1-<')(!-1/'') < (27)
< С9 |V | | ^[0,];<-'.'(С))||^||и.)т(1-',Х1-1/!" << Сют!!«,..^.
где в = (р + п)/( 2 р — 2 ), // = ( 1 — б )( 1 — 1/р ) и постоянные с не зависят от т и V. Оценим второе слагаемое, входящее в норму . Аналогично имеем:
| | ф V | | С( [0,т]; С ^(С)) < С1 1 | | ф V | | с( [ 0, т]1(С))
< с- I I ФV I I с^т]м - 2/^(С)) I I V I I ^)Т(1о01)(1о1/^ < (28)
< С1 3 I I Ф V I I £[0,т]м - 2/Р-(С)) I I V I I ЦЧТ(1о01)(1о1/р) << С14 I I V I I ,
где в 1 = ( 1 + п /р + 5 1 )/( 2 + 5 — 2 /р ), /^ = ( 1 — в 1)( 1 — 1 /р ) и постоянные с не зависят от т и V. Таким образом, используя (27), (28), оценку (24) можно переписать в виде
I I ЯМ — Й0?2) | | Яо,т < С1 5 Т0 | | V! — V2 | | Яо,т + С1 Г | | Vl — V2 | | Но,т, (29)
где постоянные не зависят от и - некоторая постоянная.
Выберем такое, что
Т о С1 5 < 1 / 4, С 1 5 = С15(Я 0 ), тогда оператор удовлетворяет оценке
I I ) — Я( V 2) I I Яо,т < I I Vl — V2 I I Яо,т/ 2 VV¿ £ ВКо ,Т < Т0 ,ЦЦ £ В Ч (30)
В частности, имеем неравенство
I I Я00 — Я( О) I I Яо,т < I I V | | Яо,т/2 Vv £ В^т < Т0,ЦЦ £ В Ч (31)
Из (23) вытекает неравенство
I I Я^) + /0 I I Яо,т < I I V | | Яо,т/2 + Я 0/ 2 < Я 0, Vv £ ВКо, (32)
которое гарантирует, что отображение ( ) переводит шар в себя и является в нем сжимающим. По теореме о неподвижной точке получим, что уравнение (23) разрешимо. Получим необходимую оценку. Пусть два решения, отвечающие двум различным векторам ЦI = ( ц..., ц^) ( / = 1, 2 ), из шара В г°. Таким образом,
т т
V; = МО 1(/(г,х,т,- + Ф,7 Ту + 7Ф)— X цуЬ ^у+^/К г,х)ц/( г)). (33)
/с=г+1 1=1
Вычитая два уравнения (33) и используя полученные выше оценки, мы выводим
I I И1 — V 2 I I Яо,т < С1 6 I I Ц 1 — Ц 2 I I с([ 0 , т ], (34)
где постоянная не зависит от .
Пусть V £ Н0т. Полагая х = Ху в (18), мы придем к системе равенств
^ + Мг,ху)— X Ц"Ь^Сгху) =
/с=г+1
(35)
т.
= /( г, ху, v( г, ху) + Ф( г, ху), 7( V + Ф )( г, ху)) + X Л ( г, ху)ц 1 ¿( г),
¿=1
где ] = 1 ,. . .,5. Эта система может быть переписана в виде ВсЦ = 1/> + /?( д), где координаты векторов и Я( д) с номерами от (/ — 1 )Л + 1 до _/' Л совпадают с векторами ^ — /( г, ху, Ф( г, ху), 7 Ф( г, ху)) и л оv( г, ху) — г+1 ц 1 к Ь кV( г, ху) — /( г, ху, v( г, ху) +
Обратные задачи для математических моделей конвекции-диффузии
с точечными условиями переопределения Ф( t, Ху), 7v(t,Xy) + 7 Ф( t, Ху)) +/(t, Ху, Ф( t, Ху), 7Ф ( t, Ху)), соответственно. В оператор ft входит функция v = v( q), которая определяется через функцию q как решение задачи (18). Таким образом, приходим к системе
¿| = В- 1^ + В- q), (36)
где оператор В- 1i?( q): С( [0,т] — С( [ 0, т] ограничен. Более того, он удовлетворяет условиям теоремы о неподвижной точке на некотором малом промежутке времени. Получим оценки.
Отметим, что по построению t | t=0 = 0. Следовательно, найдется параметр т 1 < т 0 , такой, что
I I Ч> I I с( [о,т] )<ro/2 VT<To- (37)
Пусть q7 £ С( [ 0, т] ) два вектора из шара Вг° и v 1 , v 1 - соответствующие решению задачи (18). Далее имеем, что (с = v — v2)
г
I I В- ^(qT) — В- ^(qT) I I C([0,T]) < n(£( £ I I Дао(ху,0I I c( [o,r] +
j=l \a\<2
( )( )) ( )( )) ( ))
Далее, используя теоремы вложения и оценки (25), (26), оценим правую часть этого неравенства через
( ( )) ( )
Используя рассуждения из доказательства оценок (29) и вышеприведенные неравенства, аналогично получим оценку
I I В- ^(qT) — В- ^(qT) I I С([о,Т] ) < I I V1 — V2 I I н°,т, (38)
где постоянная не зависит от и - некоторая положительная постоянная. Выберем т2 < т 1 , такое, что у3 т2 < 1 / 2. При таком выборе т2 оператор В 1ft( q) будет сжимающим и будет переводить шар В г° в себя. Следовательно, уравнение (36) разрешимо.
Пусть v есть соответствующее решение задачи (18). Покажем, что полученные решения ( v, q) есть решение обратной задачи (18), (19) эквивалентной нашей. Действительно, полагая х = Ху в (18), придем к равенствам
т1
*(ЬХу)+МЬХу)— £ ,1^,Ху) =
k=m+1
т1
= /(t, Ху, v(t, Ху) + Ф( t, Ху), 7(v + Ф )( t, Ху)) + £ /г ( t,x,0 q 1 г( t) .
¿=1
Вычитая эти равенства из (35), получим v( t, Ху) = t). Единственность решений задачи (18), (19) вытекает из полученных в процессе доказательства оценок.
Литература
1. Marchuk, G. I. Mathematical Models in Environmental Problems [Text] : V. 16: Studies in Mathematics and its Applications / G. I. Marchuk. -Amsterdam : Elsevier Science Publishers, 1986. - 217 p.
2. Прилепко, А. И. Теоремы разрешимости и метод Ротэ в обратных задачах для уравнения параболического типа [Текст] / А. И. Прилепко, В. В. Соловьев // Дифференциальные уравнения. - 1987. - Т. 23. - № 10. - C. 1791-1799.
3. Afinogenova, O. A. Stabilization of the solution to the identification problem of the source function for a one-dimensional parabolic equation [Text] / O. A. Afinogenova, Yu. Ya. Belov, I. V. Frolenkov // Doklady Mathematics. - 2009. - Vol. 79, 1. - P. 70-72.
4. Белов, Ю.Я. О задаче и дентификации функциии сточника для уравнения типа Бюр-герса [Текст] / Ю. Я. Белов, К. В. Коршун // J. of Siberian Federal University. Mathematics & Physics. - 2012. - 5(4). - P. 497-506.
5. Кулиев, М. А. Многомерная обратная задача для параболического уравнения в ограниченной области [Текст] / М. А. Кулиев // Нелинейные граничные задачи. - 2004. - № 14. -C.138-145.
6. Прилепко, А. И. О разрешимости обратных краевых задач определения коэффициента перед младшей производной в параболическом уравнении [Текст] / А. И. Прилепко, В. В. Соловьев // Дифференциальные уравнения. - 1987. - Т. 23, 1. - С. 136-143.
7. Prilepko, A. I. Methods for solving inverse problems in Mathematical Physics [Text] / A. I. Prilepko, D. G. Orlovsky, I. A. Vasin. - New-York Marcel Dekker Inc, 1999. - 709 p.
8. Пятков, С. Г. О некоторых классах коэффициентных обратных задач для параболических систем уравнений [Текст] / С. Г. Пятков, М. Л. Самков // Математические труды. -2012. - Т. 15, 1. - C. 155-177.
9. Ozisik, M. N. Inverse Heat Transfer [Text] / M. N. Ozisik, H. R. B. Orlande. - New York : Taylor & Francis, 2000. - 314 p.
10. Mamonov, A. V. Point source identification in nonlinear advection-diffusion-reaction systems [Text] / A. V. Mamonov, Y-H. R. Tsai // Inverse Problems. - 2013. - V. 29, 3. - 26 p.
11. Amann, H. Linear and quasilinear parabolic problems [Text] : V. 1 Abstract Linear Theory / H. Amann. - Berlin ; Boston ; New-York : Birkhauser, 1995. - 342 p.
12. Трибель, Х. Теория интерполяции. Функциональные пространства. Дифференциальные операторы [Текст] / Х. Трибель. - Москва : Мир, 1980. - 664 c.
13. Denk, R. Optimal Lp-Lq-estimates for parabolic boundary value problems with inhomoge-neous data [Text] / R. Denk, M. Hieber, J. Pruss // Math. Z. - 2007. - V. 257, □ 1. - P. 193-224.
14. Ладыженская, О. А. Линейные и квазилинейные уравнения параболического типа [Текст] / О. А. Ладыженская, В. А. Солонников, Н. Н. Уральцева. - Москва : Наука, 1967. - 736 c.
15. Пятков, С. Г. Об определении функции источника в квазилинейных параболических задачах с точечными условиями переопределения [Текст] / С. Г. Пятков, В. В. Ротко // Вестник ЮУрГУ. Серия«Математика. Механика. Физика». - 2017.- Т. 9. - № 4. - С. 19-26.
16. Pyatkov, S. G. On some parabolic inverse problems with the pointwise over determination [Electronic resource] / S. G. Pyatkov, V. V. Rotko // AIP Conference Proceedings. - 2017. -V. 1907. - URL: https://aip.scitation.org/doi/abs/10.1063/1.5012619.
17. Pyatkov, S. G. On some parabolic inverse problems with the pointwise overdetermination [Text] / S. G. Pyatkov, V. V. Rotko // AIP Conference Proceedings. - 2017. - V. 1907. -P. 020008.
