Пространство начальных данных второй краевой задачи для параболического дифференциально-разностного уравнения9) Текст научной статьи по специальности «Математика»

УДК 517.955.4

ПРОСТРАНСТВО НАЧАЛЬНЫХ ДАННЫХ ВТОРОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ ПАРАБОЛИЧЕСКОГО ДИФФЕРЕНЦИАЛЬНО-РАЗНОСТНОГО УРАВНЕНИЯ9

А.М. Селицкий

Вычислительный центр им. А.А.Дородницына РАН, ул. Вавилова, 40, Москва, 119333, Россия, e-mail: selitsky@mail.ru

Аннотация. Рассматривается вторая краевая задача для параболического дифференциально-разностного уравнения (п > 1). Получено описание пространства начальных данных в терминах пространств Соболева при предположении гладкости решений в некоторых подобластях исходной области, получено, что сильное решение будет существовать при начальных функциях из W21(Q). Тем самым выявлен новый класс операторов, для которых верна гипотеза Т. Като о корне квадратном из оператора. Используемые методы основаны на теории эллиптических функционально-дифференциальных уравнений и теории интерполяции пространств Соболева с краевыми условиями.

Ключевые слова: функционально-дифференциальные уравнения, теория полугрупп, теория интерполяции.

1. Введение

Первая краевая задача для параболического дифференциально-разностного уравнения со сдвигами по пространственным переменным рассматривалась в работах [?]-[?]. В работе [?] была установлена связь такого рода задач с нелокальными задачами. Вопрос сильной разрешимости в терминах интерполяционных пространств изучался в работах [?], [?] и др., точное описание этих пространств в терминах пространств Соболева получено в [?]. Оказывается, что возникающие при этом пространства начальных данных связаны с известной проблемой Т. Като о корне квадратном из оператора, построенного по секториальной форме (замечание 2.29 в [?], гл. VI, §2). В работе [?] были указаны широкие классы функционально-дифференциальных операторов, удовлетворяющих гипотезе Т. Като.

Вторая и третья краевые задачи для параболического дифференциально-разностного уравнения рассматривались в работах [?], [?]. В этих работах рассматривались области размерности п > 1. Но вопрос о явном описании пространств начальных данных оставался открытым. Третья краевая задача в случае п = 1 изучалась в работе [?], там вопрос о прастранстве начальных данных решался путем сведения к системе уравнений.

9Работа выполнена при финансовой поддержке РФФИ (проект №10-01-00837), Программы ОМН РАН «Современные проблемы теоретической математики», проект «Оптимальные алгоритмы решения задач математической физики» и Программы №3 фундаментальных исследований ОМН РАН.

В перечисленных работах было показано, что в отличии от параболических уравнений гладкость решений может сохраняться, а может и нарушаться внутри области даже при бесконечно дифференцируемой начальной функции, что и создает сложности в описании интерполяционных пространств.

В данной работе получено описание пространства начальных данных для второй краевой задачи для параболического дифференциально-разностного уравнения (п > 1) в терминах пространств Соболева при предположении гладкости решений в некоторых подобластях исходной области. Используемые методы основаны на теории эллиптических функционально-дифференциальных уравнений [?], [?] и теории интерполяции пространств Соболева с краевыми условиями [?].

Заметим, что краевые задачи для параболических функционально-дифференциальных уравнений возникают в теории нелинейных оптических систем с двумерной обратной связью (см., например, [?] - [?]).

Параболические функционально-дифференциальные уравнения с запаздываниями по времени изучались многими авторами, наиболее общий случай таких уравнений, содержащих переменные запаздывания в старших производных, рассматривался в работах [?], [?].

Автор выражает благодарность А. Л. Скубачевскому за постановку задачи, а также М. С. Аграновичу за указание на работы Р. Сили.

2. Постановка задачи

1. Пусть Q - ограниченная область в Кп (п > 2) с кусочно гладкой границей

N0 __

дQ = М|,

1=1

где М| - (п - 1)-мерные многообразия класса С“ открытые и связные в топологии дQ. Предположим, что в окрестности каждой точки х е К = дQ \ и1Ц область Q

I

удовлетворяет условию конуса, в частности, она диффеоморфна плоскому углу, если п = 2.

Введем ограниченные разностные операторы Кц, ^ : L2(Rn) L2(Rn) по формулам (^и)(х) = а^и(х + И) (И = 1,..., п),

Ием (2.1)

(^и)(х) = а|Ии(х + И) (I = 0, 1,...,п).

ИеМ

Здесь аиИ, а|И - комплексные числа, множество М состоит из конечного числа векторов с целочисленными координатами.

Введем линейные операторы Р,^, и К|(Э. Оператор ^: L2(Q) L2(Rn) является оператором продолжения функций нулем вне Q; оператор Р(з: L2(Rn) L2(Q) — оператор сужения функций на Q; операторы RIjQ, ^ : L2(Q) L2(Q) определены по формулам R|jQ = PQR|jІQ, = PQR|ІQ.

Рассмотрим дифференциально-разностное уравнение

п ( У п

Щ(хд)- Киаи^хЛ) х,+ R|QUx¡(x,t)+ RoQи(х, t) = ^хд), (х, t) е Qт , (2.2)

у=1 |=1

с краевым условием

п

R|jQUxj ообСу, Х|) = 0, (хД) е Гт, (2.3)

и=1

и начальным условием

ик=0 = ф(х), хец, (2.4)

где Qт = Qx (0, Т), 0 < Т < оо, Гт = (ЭQ \ К) х (0, Т), V - единичный вектор внешней

нормали к Гт, f е L2(QT), ф е L2(Q).

2. Пусть W2(Q) - пространство Соболева комплекснозначных функций из L2(Q), имеющих все обобщенные производные вплоть до к-го порядка из L2(Q) с нормой

/2

|a|<k Q

Введем полуторалинейную форму aR[v, w] в L2(Q) с областью определения W2(Q) по формуле

n ( У n

aR[v, w]— RijQVxj ,wxi L2(Q) + (RiQVxi, w)L2(Q) + (RoqV, w)L2(Q) . (2.5)

i,j —1 i—1

Разностные операторы RUQ, RiQ, R0Q: L2(Q) L2(Q) - ограниченные. Поэтому существует число c0 > 0 такое, что

|aR[v, w]| < Со v w2(Q) w w2(q) (v-we W2(Q))" (26)

Поскольку полуторалинейная форма aR[v, w] непрерывна по w в W 1(Q), существует линейный ограниченный оператор AR: W2(Q) [W2(Q)] такой, что

(CARv, w) — aR[v, w] (v, w e W2(Q)), (2.7)

где [W2(Q)] -сопряженное пространство к W2(Q).

Определение 2.1. Оператор AR называется сильно эллиптическим, если существуют числа с1 > 0 и с2 > 0 такие, что

Re(CArv, v) > c v W21(q) - С2 v L2(Q) (v e W2(Q)). (2.8)

В дальнейшем мы будем предполагать, что оператор AR, соответствующий уравнению (2.2), сильно эллиптический. Необходимые и достаточные условия сильной эллиптичности оператора AR в алгебраической форме сформулированы в работе [?] (леммы 2.2 и 2.3). В этом случае естественно назвать задачу (2.2)- (2.4) второй краевой задачей для параболического дифференциально-разностного уравнения.

Далее мы будем предполагать, что в неравенстве (2.8) с2 = 0. В противном случае положим и = z еС2(.

3. Сильные решения

1. Обозначим через W2к’°(От) пространство функций и е L2(QT), имеющих все обобщенные производные по х вплоть до к-го порядка из L2(QT ) с нормой

{ Г1'2

^5'0№)= ^ Ю*\(х,ч12 с1хс1^ .

|а|<к ^

Определение 3.1. Функция и е W2’°^т) называется слабым решением задачи (2.1)-(2.4), если

-ЫУ( + + ^их> + (^ и)у dxdt =

Qт І’J =1 1=1

= ^ dxdt + Ф^=о dx (3.1)

Qт Q

для всех V е W21(QT) таких, что v|t=T = 0.

Тождество (3.1) можно получить из уравнения (2.2), используя формальное интегрирование по частям.

Введем гильбертово пространство

W = V е L2 (0,Т; W2^)У : vtеL2 0,Т; 52^)

с нормой

1 Т Т 1 '2

'' W = 0 V 5,2 (Q)dt + 0 Vt ^ДО)]1 ^ ■

Здесь производные понимаются в смысле распределений на QT.

Теорема 3.1. Пусть оператор AR сильно эллиптическ1(й. Тогда задачу* (2.1) - (2.4) имеет единственное слабое решение ие Wдлявсех f е L2 0, Т; [5 2^)] и феL2(Q).

П

П

Доказательство см. в [?] (теорема 3.1).

2. Для того, чтобы сформулировать определение сильного решения задачи (2.2) -(2.4), введем понятие т-секториального оператора.

Определение 3.2. Полуторалинейная форма а^, w] называется секториальной формой с вершиной Y, если ее числовая область значений является подмножеством угла вида

п

{^ є С : | агд(^-у)! < 0}, 0 <0 < ^ , YЄR.

Лемма 3.1. Пусть оператор АК сильно эллиптический и с2 = 0. Тогда форма а^, w] является секториальной с вершиной 0.

Доказательство см. в [?] (лемма 4.1).

Определение 3.3. Линейный оператор А называется т-секториальным с вершиной Y, если оператор В = А + а1 для некоторого а є К удовлетворяет условиям:

1) существует ограниченный оператор (В + XI )-1 при ReX > 0,

2) (В + XI )-1 < ^еХ)-1, при ReX > 0,

и числовая область значений оператора А является подмножеством некоторого угла вида

п

| агд(£ -Y)! < 0, 0 < 0 <- .

Поскольку W2I(Q) - полное пространство, форма аК[V, w] замкнута. Из леммы 3.1 данной работы и первой теоремы о представлении (см. [?], гл. VI, §2, теорема 2.1) следует существование т-секториального оператора Дк : D(AR) L2(Q), D(AR) с L2(Q), с вершиной 0 такого, что

аК[v> w] = (^^ w)L2(Q) (v є D(ДR)> w є W2^)); (3.2)

более того, D(AR) плотно в W21(Q). Заметим, что в силу (2.7) соотношение (3.2) можно записать в виде

АКУ = ДRV (V Є D(ДR)). (3.3)

В пространстве D(ДR) введем скалярное произведение по формуле

w)D(AR) = (ДRV> ДRw)L2(Q) + W)L2(Q) ■

Поскольку оператор ДК замкнутый, пространство D(ДR) является гильбертовым. Введем гильбертово пространство

W(ДR) = 'м є L2(0, Т; D(AR)) : wt є L2(Qт)

со скалярным произведением

Т Т Т

мМдр) = (ДRV, ДRW)L2(Q) dt + (V, + (vt, м^2(0) dt.

0 0 0

Определение 3.4. Слабое решение u(x, t) задачи (2.2) - (2.4) называется сильным решением, если u е W(Ar).

3. Для доказательства существования сильного решения мы используем теорию полугрупп.

Определение 3.5. Сильно непрерывная полугруппа операторов {Tt} (t > 0) в гильбертовом пространстве H называется сжимающей, если Tt < 1 (t > 0).

Обозначим Лш = {z е C : | argz| < ш}, где 0 < ш < п.

Определение 3.6. Семейство линейных ограниченных операторов {Tz} (z е Лш) в H называется аналитической полугруппой вЛш, если 1) функция z Tz аналитическая в Лш,2) T0 = I и lim Tzx = x (x е H), 3) Tzi+z2 = Tzi Tz2 (z1, z2 е Лш). Полугруппа

z^0, 7бЛш

{Tt} (t > 0) называется аналитической, если существует аналитическое продолжение Tz оператор-функции Tt в некоторый сектор Лш.

Теорема 3.2. Оператор (-Ar) является генератором аналитической сжимающей полугруппы.

П Оператор Ar - m-секториальный с вершиной 0. Из теоремы 1.24 в работе [?], гл. IX, §1, следует, что оператор (-Ar) является генератором аналитической сжимающей полугруппы. И

Пусть Х0 и X1 - банаховы пространства и X1 с X0 линейно и непрерывно. Обозначим через S полосу

S = {z : 0 < Re z < 1}.

Через H (X1, X0) обозначим пространство непрерывных ограниченных функций f :""5 Х0 аналитических в 5, таких что f (j + iy) е Xj и f (j + iy) непрерывна как Xj-значная функция от y, j =0, 1, -то <y< +то, с нормой

f H(X1, X0) = n=0JX1 suP f (j + iy) Xj .

j 0, 1 — то <y <+ oo

В дальнейшем нам понадобится интерполяционное пространство

[X1, X0]e = {а е X0 : 3f (z) е H (X1, X0) с f (0) = a}

с нормой

а [X1, Xü]0 = f(¡J=a f H(X1, X0) .

Также нам понадобится неравенство

f(0 + iy) [X1,X0]0 < f h(X1,X0) , (3-4)

справедливое при всех действительных у (см. неравенство (1.4) в [?]).

Теорема 3.3. Задача (2.2) - (2.4) для любых f є L2(QT) и ф є ^(Д^, L2(Q)]1/2 имеет единственное решение, которое определяется по формуле

t

и(х, ^ = ^ф(х) + Т^(х, s) ds, (3 5)

0

где {Tt} ^ > 0) - аналитическая полугруппа с генератором (-Д^.

Доказательство см. в [?] (теорема 4.2).

4. Пространство начальных данных

1. В теореме 3.3 принадлежность начальной функции пространству ^(Д^, L2 ^^/2 является необходимым и достаточным условием существования сильного решения. Таким образом, возникает проблема описания этого интерполяционного пространства. В случае первой краевой задачи для параболических дифференциально-разностных уравнений в работе [?] при некоторых предположениях относительно области Q и сильно эллиптического оператора ДR было доказано, что выполняется равенство

^(д„), ыо)],/2 = \W2cq).

В [?] этот результат был обобщен на более широкий класс сильно эллиптических функционально-дифференциальных операторов, в частности, были сняты предположения, наложенные в [?], при этом использовался тот факт, что разностный оператор отображает \^^) в W2^), что не верно в случае второго краевого условия. В данной работе получен результат аналогичный [?], только для описания интерполяционного пространства использовался не только метод следов, как в [?], но и комплексный метод.

2. Для того, чтобы сформулировать условия, накладываемые на область Q и сильно эллиптический оператор Д^ введем некоторые обозначения. Обозначим через G аддитивную группу, порожденную множеством М. Обозначим через Qr открытые связные компоненты множества Q\ и (ЭО + її).

ИєЄ

Определение 4.1. Множество Qr называется подобластью. Множество К всех подобластей Qr называется разбиением области Q.

Разбиение Р естественным образом распадается на непересекающиеся классы. А именно, подобласти Qr1, Qr2 є R принадлежат одному и тому же классу, если существует вектор h є G такой, что Qr2 = Qr1 + К Обозначим подобласти через Qsl, где s -номер класса ^ = 1, 2,...), а I - номер подобласти в s-ом классе. В силу ограниченности области Q каждый класс состоит из конечного числа N = N ^) подобластей Qsl и N ^) < (<3іат Q + 1)п. Вообще говоря, множество классов счетно. Обозначим через Іі^ вектор h є G такой, что Qs1 + h = 05|. Очевидно, ^1 = 0.

Рассмотрим следующее множества:

К = Qn (Эо + h1) п [(ао + И2)\(ЭО + М .

И1, єG

Для простоты доказательств будем предполагать, что К с К.

Приведем вспомогательные результаты лемм 7.3-7.5 работы [?], §7.

Лемма 4.1. А) Пусть х0 є дQsI nдQ. Предположим, что существует последовательность точек {хк} таких, что хк х0 при к ^ то и хк є Qsк1 к, ^к, 1к) = (s, I). Тогда

х0 є К.

Б) Пусть х0 є Цп ЭQpIn 0Qqm, (р, I) = (q, т). Предположим, что существует последовательность точек {хк} таких, что хк х0 при к ^ то и хк є Qsкк, (sk, lk) = (р, I), т). Тогда х0 є К.

Обозначим через Гр открытые связные (в топологии дQ) компоненты множества ЭQ \ К. Очевидно, Гр є Сто.

Лемма 4.2. Пусть (Гр + И) п Q = 0 для некоторого И є G. Тогда либо Гр + И с Q, либо существует Гг с дQ \ К такое, что Гр + И = Гг.

В силу леммы 4.1 мы можем разбить множество {Гр+И : Гр+И с^, р = 1, 2,...; И є G} на классы следующим образом. Множества Гр1 + И1 и Гр2 + И2 принадлежат одному и тому же классу, если 1) существует И є G такой, что Гр1 + И1 = Гр2 + И2 + И и 2) в случае ГР1 + И1, Гр2 + И2 с дQ направления внутренних нормалей к дQ в точках х є ГР1 + И1 и х - И є Гр2 + И2 совпадают. Множество Гр с дQ может принадлежать только одному классу, а множество Гр + И с Q — не более, чем двум классам. Обозначим множество Гр + И через Г^, где г — номер класса, а j — номер элемента в данном классе (1 =J (г)). Без ограничения общности будем предполагать, что

Гг1, ■ ■ ■ , ГrJo с Q, ГГ\_І0+1 , ■ ■ ■ , Г^ с ^ (0 < ^ = -І0(г) <(г)).

Теорема 4.1. Пусть мп-1(К п дQ) = 0 и К с К, числа S0 и г0 различных классов подобластей Ssl и поверхностей Ггт конечны, и каждая подобласть Qsl (s = 1, ..., S0;

I = 1, ..., N ^)) липшицева. Пусть также оператор Др сильно эллиптический с с2 = 0 и решение задачи ДR = f принадлежит W22(QsI). Тогда

р(до, w^ (3.6)

П 1. Опишем пространство D(AR). По условию теоремы О(Ар) с W 12^) nW2^0. Очевидно, что ^их е W21(Qsl). Проинтегрируем по частям следующее выражение по

области 0^:

П

(RijQUxj)Vx¡ + (^аих> + (Коаи}\7 dx =

¡=1

о*

(^иоих_|)х* + ^¡О^ + (^0Ои) V dx +

¡,j =1 ¡=1

+ (RuQUxi)V С0Б(У, x¡} dx (3.7)

^ j ЗОз! \К

ЭОз:\К ¡'■=1

для всех V е W,1(0). После суммирования по всем s и I получаем

П

(R¡jQUxj)Vx¡ + (К^> + (Коои)у dx =

0 и=1 ¡=1

П

= ^ dx + (R¡jQUxj}V ^ ^ dx, (3.8)

, ■ j дОз! \К

о ^ ЭОз!\к и 1

п

где f е L2(Q} определена по формуле f (x} = - (R¡jQUxj)x¡ + (R¡QUx¡} + (RоQU)

¡,■=1 ¡=1

при x е Qsl. Тогда, так как u удовлетворяет тождеству (ARU, v)L2(0) = v)L2(0), то

п

(R¡ соб(у, x¡} dx = 0. (3.9)

Sl | ¡=1 д Оз! \К

s’l ЗОз!\К и '

Таким образом, функция u е W 2(0) пW 2(0^) принадлежит D(AR) тогда и только тогда, когда выполняются условия (4.4). Условия (4.4) определяют замкнутое линейное подпространство в W2(0) |(0^) (более подробно см. [?]). Обозначим через W2 р(0)

пространство функций, принадлежащих W2(0) гW|(0^), и удовлетворяющих условиям (4.4). Так как 0 е с(А^, и оператор AR : W|^(0) L2(0) ограниченный, по теореме Банаха об обратном операторе нормы U D(AR) и U w2,R(Q) эквивалентны в

Так как необходимым и достаточным условием принадлежности функции u е W2(0^) пространству W21(0) является совпадение следов на смежных границах, то пространство W2^(0) — это подпространство пространства W|(0^) с условиями (4.4) и условием совпадения следов на смежных границах подобластей 0^.

2. Докажем теперь, что W2’R(0}’ L2(0) 1/2 = W21’R(0), где W2^(0) — это пространство W21(0sl} с условием совпадения следов на смежных границах подобластей 05|, т. е. пространство W2(0). Для доказательства воспользуемся комплексным методом интерполяции, см. [?]. Пусть u е W^(О), L2(0) 1/2 и f е Н W|’R(0}’ L2(0)

п

п

п

п

п

такая, что f(1/2) = и. Так как W|,к(0) с W22(Qsl). Из неравенства (3.4) и того, что W22(Qsl)I L2(Q) 0 = W220(Qsl) (в силу интерполяционных теорем 7.1, 9.1 и 9.2 в [?] и теоремы 5 в [2], гл. VI, §3 о продолжении функций для областей с липшицевой границей), следует, что сужение функции f на Ре ї < 1 —непрерывное отображение в W2(Qsl) аналитическое в 1 < Ре z < 1. Если обозначить Т оператор взятия следа, то отображение z Tfнепрерывно в L2(0Qsl) при 0 < Реї < 1 и аналитическое при 0 < Ре z < 1. Таким образом, функция f пинадлежит W2^0 и удовлетворяет условию совпадения следов на смежных границах подобластей Qsl при Реї >0, а значит и є W21’R(Q). И

Литература

1. Лионе Ж.-Л., Мадженес Э. Неоднородные граничные задачи и их приложения / М.: Мир, 1971.

2. Stein E.M. Singular integrals and differentiability properties of functions // Math. Series.

№30 f Princeton: Princeton Univ. Press, 1970.

INITIAL DATE SPACE FOR SECOND BOUNDARY VALUE

PROBLEM OF PARABOLIC DIFFERENTIAL-DIFFERENCE EQUATION

A. M.

Selitskii Dorodnitcyn Computing Centre of

RAS,

Vavilova St., 40, Moscow, 119333, Russia, e-mail: selitsky@mail.ru

Abstract. The second mixed problem for parabolic differential-difference equation (n > 1) is studied. The description of initial data space is obtained in the terms of Sobolev's spaces in the case when the solution is smooths in some subdomains of original domain. It is obtained also that the strong solution exists if the initial function belongs to W 1(Q). So, the new class of operators satisfying the Kato hypothesis relative tc? the operator square root is found. We have used the theory of functional differential equations and the theory of interpolation of Sobolev's spaces with boundary conditions.

Key words: functional differential equations, theory of semigroups, interpolation theory.