О существовании решения параболической задачи на графе с краевыми условиями, содержащими производные по времени Текст научной статьи по специальности «Математика»

Владикавказский математический журнал 2011, Том 13, Выпуск 3, С. 42-52

УДК 517.956

О СУЩЕСТВОВАНИИ РЕШЕНИЯ ПАРАБОЛИЧЕСКОЙ ЗАДАЧИ НА ГРАФЕ С КРАЕВЫМИ УСЛОВИЯМИ, СОДЕРЖАЩИМИ ПРОИЗВОДНЫЕ ПО ВРЕМЕНИ

Р. Ч. Кулаев

Рассматривается начально-краевая задача для параболического уравнения, заданного на геометрическом графе, с краевыми условиями содержащими производную по времени. Устанавливается теорема существования решения краевой задачи, дающая представление решения в виде контурного интеграла.

Ключевые слова: краевая задача на графе, метод контурного интеграла.

В настоящей работе рассматривается смешанная задача для уравнения д д ( д \

т и = д^ + + (М)еГх[0,Т] = Гг, (1)

заданного на геометрическом графе Г.

В каждой граничной вершине а графа Г решение уравнения должно удовлетворять условию

дд

а(а) —— и(а, ¿) + (3{а) — и(а, ¿) + 5(а) и(а, ¿) = 0, а е <9Г, £ > 0. (2)

д X д и

А в каждой внутренней вершине а на решение уравнения (1) накладываются условия непрерывности и одно условие согласования

ик(а, и) — ико (а,и) = 0, к,к0 € I(а),

дд ак(а)—ик(а,1) + рк(а)—ик(а,1) + 5к(а)ик(а,1)= 0, а € V.

дх кХ ' 7 ' дЬ

кет (а)

В начальный момент времени Ь = 0 ставится условие

и(х, 0) = 0. (4)

В уравнении (1) полагаем р € С 2[Г], т^ег р(х) > 0, д € С [Г], / € С2 [Г] П С [Гт]. В условиях (2), (3) {ак(а)}ке1(а), {вк(а)}¿6/(а), [$к(а)}ке1(а) — наборы вещественных чисел, свои для каждой вершины графа. При этом считаем, что ак (а) > 0, вк (а) ^ 0, к € I(а), для любой вершины. В условиях согласования (3) считаем, что все производные подсчитаны в направлении к вершине а.

2011 Кулаев Р. Ч.

Рассматриваемая задача имеет естественную физическую интерпретацию. Она моделирует процесс распространения тепла в системе стержней, копирующей граф Г. В узловых точках системы стержни могут быть спаяны или соединены муфтами, каждая из которых имеет свою теплоемкость, что выражается условиями (2), (3).

Отметим, что в работе автора [1] установлена разрешимость параболической задачи для случая однородного уравнения (1) при неоднородных краевых и начальном условиях (2)-(4).

1. Предварительные сведения

Под геометрическим графом в настоящей работе понимается одномерное стратифицированное многообразие, вложенное в Ж" и обозначаемое через Г [2]. Ребра графа — это пространственные гладкие кривые, не имеющие самопересечений. Вершина графа — точка, являющаяся концом одного или нескольких ребер. Считая ребра графа занумерованными, обозначаем их через 7и, а вершины — через а^ или bj (при этом предполагается, что нумерация вершин независима от нумерации ребер). Через V обозначаем множество вершин графа, которые являются концевыми точками двух и более ребер. Такие вершины мы называем внутренними. Вершины графа, не принадлежащие V, называем граничными и обозначаем их через дГ. Если вершина а является концевой точкой ребра 7и, то будем говорить, что ребро 7и примыкает к вершине а. Множество индексов всех ребер, примыкающих к внутренней вершине а, обозначим I(а). Всюду далее полагаем, что граф Г является конечным и связным множеством в Ж".

Под функцией на графе понимается отображение и : Г — С. Через и и будем обозначать сужение функции и на ребро Yi. Если а — произвольная вершина (граничная или внутренняя) графа Г, то под ик (а) понимается ик (ж), ж £ 7к. Через С [Г] будем

обозначать пространство функций, равномерно непрерывных на каждом ребре графа. Мы не требуем непрерывности функций пространства С [Г] на всем графе Г, поэтому величины ии(а) и иц(а) не обязаны совпадать при к = г. Через С (Г) — подпространство из С [Г] функций непрерывных на всем графе.

Введем понятие производной функции. Для этого введем в рассмотрение функцию ф(ж), взаимно однозначно отображающую каждое ребро ^и на некоторый интервал (0,1и) при 1и > 0. Функцию ф(ж) будем называть метрической, а величину 1и — длиной ребра 7и. Для сужения функции ф(ж) на ребро 7и существует обратное отображение ж(ф) интервала (0,1и) на ребро 7и. Функцию и(ж) назовем дифференцируемой на графе Г, если для каждого ребра 7к дифференцируема по переменной ф суперпозиция и(ж(ф)). При этом полагаем и'{ х) = Обозначим через Сп[ Г] пространство функций и(х) £ С [Г],

производные которых до порядка п включительно существуют и принадлежат пространству С [Г], а через Си (Г) обозначаем Си [Г] Р| С (Г). Далее мы будем рассматривать ком-плекснозначные функции нескольких переменных, у которых одна (пространственная) переменная имеет областью своего изменения граф Г. Области изменения других переменных могут быть различны. Везде ниже полагаем, что все рассматриваемые функции равномерно непрерывны по пространственной переменной на каждом ребре графа. Рассматривая отображения вида и : Гт — С , через Сга>5[Гт] обозначаем пространство всех функций, определенных на Гу, имеющих равномерно непрерывные производные д^п^ на каждом множестве 7и х [0,Т], 7и £ Г.

Под интегралом функции и £ С [Г], взятым по графу Г, понимаем сумму криволинейных интегралов по всем ребрам графа.

В дальнейшем нам понадобятся следующие результаты работы [3].

Спектр соответствующей (1)—(4) спектральной задачи состоит из последовательности собственных значений, не имеющей конечной предельной точки. Все собственные значения расположены в полосе П конечной ширины с центром в начале координат комплексной плоскости, границы которой параллельны вещественной оси. Обозначим через Д(р) характеристический определитель спектральной задачи, где р = \/—Х. Выбросим из р-плоскости внутренности малых кругов радиуса г с центрами в нулях рк функции Д(р) и обозначим через Сг, оставшуюся часть р-плоскости. Пусть р £ Сг. Обозначим через С(ж,£,р) функцию Грина [2, 4] интегрального оператора, обращающего спектральную задачу. Функция Грина удовлетворяет следующему асимптотическому представлению [3]

ы т ы т Г«(1)(г)Г

dx

k=1

л/Рк{х)

x ( 4 (!)

exp [ip j QV(z)<b) exp ( - ip | e<2>(*)ds) [cf(6]

ak bk

ak

+ exP (ip I QM(z)dz ) + (¿P)"-1 £

V'Pkix)

bk , 4

exp - ip J ©(2)(z)d^ J exp - ip J ©(2) (z)dz^ [zk2 (0]

42)(p)

A(p)

ak

bk

+ exP(ipj ©(1) (z)d^ [Zk!) (e)]

4f(p)

A(p)

(5)

где функции E(S)(p)/A(p) ограничены некоторым числом в Cr,

dx!

■Kk (x,e,p) = <

0, e G Yn, x е Yj, n = j;

(ip)v-1 [41 (x, e)] exp (ip f ©(1) (z)dz) ,

Yk = [ak,bk], ak ^ e ^ x ^ bk; -(ip)v-1 [pki (x,e)] exp (ip f ©(2) (z)dz^ ,

Yk = [ak,bk], ak ^ x ^ e ^ bk,

(6)

e^Oz) = 1/ y/p(z), e®(z) = -l/y/p(z) при Imp )0и e^Oz) = -l/y/p(z), ©(2)(z) = l/^Kz) при Imp < 0; = (2©^ Ю)"1, = (©^(аОГЧ^(6/

\/pk{x), s = 1, 2, v = 0,1, 2, а квадратные скобки [p] обозначают сумму р + 0(1/р).

х

2. Существование решения задачи (1)—(4)

Через П мы по-прежнему обозначаем полосу комплексной р-плоскости, содержащую все собственные значения соответствующей спектральной задачи. Пусть гр £ (0, Обозначим через Ь разомкнутый контур, лежащий в верхней полуплоскости комплексной

р-плоскости, не имеющий общих точек с множеством П, и такой, что прямые arg р = ф и arg р = п — ф являются асимптотами. Обозначим через On последовательность окружностей радиусов rn с центром в начале координат р-плоскости, целиком лежащих в Cr, причем rn — при n — Пусть, далее, Ln есть часть контура L, расположен-

ная внутри окружности On. Полагаем, что L\Ln совпадает с достаточно отдаленными частями прямых arg р = ф и arg р = п — ф.

Теорема. Если f £ C2'2[Гу], то решение задачи (1)-(4) существует. Это решение можно представить в виде

t

u(x,t) = -¿/р/ G(x,Z,p) j e-^-^fiMdr^dß. (7)

L Г 0

< Для доказательства теоремы достаточно показать возможность дифференцирования под знаком интеграла (7).

Проинтегрируем по частям интеграл

t

1 , 1 ^ 1 d

J ехр ( - p\t - г)) /(£, r)dr = j2 /(£, i) - 1 e~p2t /(£, 0) - 11 / (£, t)

t

1д

(8)

р4 dt

Рассмотрим интеграл

t

I j j ^ m-r)*rMi>

L Г 0

= ¡"wwh^ {7mt) - 7e"hmo)

L Г

t=0

(9)

t

Р

Как видно из асимптотических представлений (5) и (6) и формулы (8) при 11т р| главной частью подынтегральной функции (9) является

1 ( (£ 1 С(х,£,р)П£,1) *

р dts dxv J ' ' ' p(£)

Г

k=1 ^k £

bk £ ч

-Jexp^-ipJ e^(z)dzy^(x,Ofk(^t)^\

—>

k-1 w i

bk i 4

- j exp (-ip j (10)

x x '

Интегрированием по частям можно убедиться в том, что (10) при v ^ 1 на контуре L

убывает как р-2, равномерно относительно t ^ 0 и x е Г. Следовательно, интеграл (9)

при v ^ 1, s ^ 2 сходится равномерно относительно t ^ 0 и x е Г. Это дает нам право

дифференцировать под знаком интеграла (7) один раз по x и два раза по t и переходить

к пределу при t ^ 0, x ^ а, а е 0Г U V. Поэтому lim u3(x,t) = 0. Далее, принимая во

t^ü

внимание (9) и возможность дифференцирования под знаком интеграла один раз по x и один раз по t, получим:

t

u(x,t) = S p S G(x,Z,p) i e-p2(t-r)f{^T)dT^-

dtsdxv K ' J J ' dtsdxv J v J J

dp

/■ 1 dv i ds df = С(ж, p) T-— /(£, i) -77- dp

J p dxv J v 7 <9is J vs' 7 F

L Г

L Г

Г 1 / ds+1 df

L Г

L Г

L

Ь Г 0

В силу (9) у первого интеграла по контуру Ь в правой части (11) подынтегральная функция убывает как р-2, следовательно, представляя его в виде предела интеграла по контуру 0+, убеждаемся, что он равен нулю. Точно так же убеждаемся, что равны нулю третий и пятый интегралы в правой части (11).

Подставим функцию и(ж,£) из (7) в краевые условия (2), (3). Пользуясь линейностью функционалов, определяющих условия (2), (3), разобьем выражение, получающееся после подстановки функции и(ж, ¿) в краевые условия, на шесть частей, соответствующих каждому слагаемому в (11). Все нечетные слагаемые равны нулю. Равенство нулю четных слагаемых следует из того, что функция Грина С(-,£,р) удовлетворяет однородным краевым условиям (2), (3).

Таким образом, функция, определяемая формулой (7), удовлетворяет всем условиям (2), (3). Поэтому для доказательства теоремы нужно показать, что (7) удовлетворяет уравнению (1).

x

x

Поскольку мы не имеем права дифференцировать два раза под знаком интеграла (7), выделим в нем главную часть. Для этой цели введем обозначение

G(1) (x,С,р) =

0, х £ Yn, С £ Yj, n = j;

(гр)-1рк£1 (х, С) exp ^¿р / в(1) (z) dz

х,£ е = [ak,bk], ak ^ ж ^ bk; (12)

— (гр)-1 (х,С) exp (VIв(2) (z) ^ ,

G 7fc = [a>k,bk], ak ^ ж ^ bk,

G(2) (x, С, р) = G(x, С, р) — G(1) (x, С, р). (13)

Обозначим через Dx — дифференциальный оператор Dx = р{х) + -§^.р{х) + q{x). Тогда

d

—жг Dx u(х, t) +7ri — щ (ж, t)

dt

Dx —

L Г 0 (14) t

+ / p J j e-p2^f(t,t)dt^dp.

L Г

Согласно теореме 3 подынтегральная функция второго слагаемого в (14) убывает как р-2 на контуре Ь. Следовательно, этот интеграл допускает дифференцирование два раза по х и один раз по Ь. Тогда, очевидно,

-тгг (dx - u3(x,t) = (dx-^jpj gV(x,z,p) j e-^-^fiMdr

L Г

t

+ j p(dx + p^j j g^(x,^,p) j e-^-^fimdr^dp

L Г 0

- jpj g^(x,c,p)f(c,t)dt^dp.

(15)

L Г

Подынтегральная функция третьего слагаемого на контуре Ь убывает как р-2, поэтому этот интеграл равен нулю. Тогда, принимая во внимание свойства функции Грина и (13), из (15) получаем

u(x,t) = (dx -j р j GV(x,£,p)j

L Г

- jp(dx + p2)jg^(x,^p)j e-^-^fimdr^dp jpj e-^-^f{i,t)dtdp.

L0

t

t

t

t

t

Ввиду того, что на контуре Ь функция

I С(1) (ж,£,р)| е-р2(4-т)/(£,т) ^т

убывает как р 3, во втором слагаемом (16) дифференциальный оператор можно вынести за знак интеграла по Ь.

-7П - п3 (ж, ¿) = ^ + <?(ж)) j Р j °{1) (ж> Р)

Ь Г

//С«(ж,е,р) (17)

0 Ь 4 Г

00

Теперь вычислим выражение в фигурных скобках из (17). Для сокращения записи введем обозначение

^^у1, (18)

Как видно из определения (см. конец § 1)

Ф? (ж,е,т ) = -Фк5)(ж,е,т). Принимая во внимание (12) и (13), непосредственным вычислением находим

í

р{х)&+р2) / / е~рч'~т)Шт)(1тШ

1к 0

{х х

1|ехр (гр I

«к 4

4 Ьк х

х i 1 ехр (гр | е^*) ^р^ж^)

0 х 4

00

(х х

j ехр^гр^ ©(1) (г) ^ фЦ (ж,£,т) ^

"«к

4

Ьк х ч

+ У ехр^грУ ©(2) (г) Ф^ (х,^,т) (Л (т

х £ '

* Г

= 1 е-р2(1-тН Рк (х){©к)(х) Ф^ (х,х,т)+©к2) (х) фк2) (х,х,т о

(х х

¿е^Ог) I ехр (гр I ©^(^^

ак £

х х ч

/ ехр (¿р I в®^)

Ьк £ '

(х х

©^(ж) I ехр ^р I ©«(г)^

ак £ х х ч

+в!2)(ж) I ехр (гр I

Ьк £ '

(х х

/ехр (^/©«(^¿Ф^т)^

ак £ х х ч ч

ехр ^р I ©(2) (г) ^ Ф12) (ж, т) 1 I йт. (19)

Ьк

Из определения и (18) следует, что

Рк(х^©к1)(х)фк1) (х,х,т)+©к2)(х) Фк2) (х,х,т^ = /к(х,т). Подставляя (19) в (16) и принимая во внимание (12), получим при х € 7к

* , х

{о, -|) ифл) = £ - |) / /

о ак

х х

х У ехр ^ — р2(Ь — т) + ©(1) (г) + У фк2) (х, т)

Ь £ Ьк

х ^ *

х у ехр р2(1-т)+гр ! в^^й^йрйЛйт - ! |№(ж) ^ ©^ (ж)

Ь £ ^ 0

х х

^У Ф(1) (хЛ,т) ! гр ехр ^ — р2(Ь — т) + 1р ^ ©(1)(г) й^йрй^

ак

+^©12)(ж) I Ф^2)(ж ,£,т) 1гре^[-р2{1-т) + гр I Ьк Ь 4

х х

+2рк{х) (©11}(ж) У I ípexp(-p2(t-т) + íp ^

ак Ь 4

х х

+©12)(ж) I АФ(2)(Ж)^Г) I гРехр(-р2^-т)+гр ^ ©(2)(г) ¿г) йрд^

Ьк Ь 4

х х

+рк(х)[1 ^ ф£°(а:,£,т) I ехр (-р2(;£-т) + гр I ©(1)(,г)(1г)ф^

ак Ь 4

х х \

+ I (ж, 6 т) I ехр(-р2^-г)+гр ! {г) \йр. (20)

Ьк Ь 4

Заметим, что во всех интегралах, встречавшихся в (20), вещественные части степеней показательных функций отрицательны при р £ Ь, что дает хорошую сходимость этих интегралов, которые в общем виде можно записать следующим образом:

1к = J(гр)к ехр ( — р2а + грЪ) йр = ^{%р)к ехр (гра + <^р,

ЬЬ

где а > 0, Ь ^ 0, к = 0,1.

Представляя интеграл в виде предела интегралов по контурам Ьп и вводя новую

Ь2

переменную г = гра + можно показать, что

22

т V71" Г V71" - /щ\

/0 = ^е 4а, Д =--4а 21

л/а л/а 2 у/ а

Тогда, вычисляя по формуле (21) интегралы, входящие в (20), получим

, , , , & о

-Ш п3(ж, ¿) = ^ [Р*(х) ^ - ^

0 7к 4 4

_ í ^ ж

г/^ /»ьц*?^^/

0 у 7 Чк 4

х х

+2©« (ж) (I вМ (г) ± Ф« (ж, г)) ехр ( - (| ©« (г) ^

4

-2(* - г) I рк (ж) ^ ф£> (ж, г) ехр ^ - (I в? (г) йг^ ^ | дт.

7к 4 4

х

х

Вычисляя оператор (рк(ж) — -щ) от первого интеграла, получим

д

-т - — ^и3(х,Ь)

1к 4 £ 7 = —л/ттг Ит

1к £

= -ЛИ Ит / + ( / _N 2

1к £

1к £

Здесь мы воспользовались определением функции Ф^ и условием / € С2[Гу]. Покажем, что предел в правой части последнего равенства равен —-^/(ж, ¿). Для этого разобьем интеграл по 7к на три части

х—е х+е Ьк /=/ + / + /' 0 <£« 1

1к ак х—е х+е

В силу того, что показатели экспонент при £ € (ак,х — е) и (х + е,Ък) отрицательны, а функция Ф^(х,£,Ь) имеет непрерывные производные по х и £ на [ак,Ък], то первый и третий интегралы стремятся при ¿1 ^ 0 к нулю. Рассмотрим второй интеграл:

х+е (1) х „

х—е £

х+е (0) х 0

= Г №*.е)М(.г) 1 С 1 (Г <Ь \у 1 ^иК) V «1 V/ v5W/ )

х— е £

рЙ^х,х)/к(х,г)\ ( ( } йг У^ _

бхр

л/ртЩь1 л/Рк(х)Ь ) V уШ/ / л/РкШ

(0) х+е х 2

¿к!(х,х) /к(х,и) [ 1 ( 1(1 (г ^ ^ (£

+2 ——-■==- / —— ехр — -— 11 11

л/рк{х) У 2^ \ V\/рк(г); ) л/Щ)

х—е £

Из определения функции и условия /к € С2[^к] следует, что при е ^ 0 и Ь1 ^ 0 первый интеграл стремится к нулю. А второй интеграл путем замены переменной сводится к интегралу Эйлера — Пуассона и стремится к —\рк при ¿1 —0.

х—е

Таким образом, получаем при x G Yk

/4i(x,x)fk(x,t)

(0)

2

«

/k(x,i). >

Литература

1. Кулаев Р. Ч. Теорема существования для параболической задачи на графе с краевыми условиями содержащими производную по времени // Тр. ин-та математики и механики УрОРАН.—2010.— Т. 16, вып. 2.—2010.—С. 139-148.

2. Покорный Ю. В. и др. Дифференциальные уравнения на геометрических графах.—М.: Физмат-лит, 2004.—272 с.

3. Кулаев Р. Ч. Асимптотика решения краевой задачи на графе // Изв. вузов Сев.-Кавк. регион. Естеств. науки.—2010.—Вып. 4.—С. 17-21.

4. Покорный Ю. В., Карелина И. Г. О функции Грина задачи Дирихле на графе // Докл. АН СССР.—1991.—Т. 318, № 3.—С. 942-944.

Статья поступила 7 сентября 2010 г. кулаев руслан черменович

Южный математический институт ВНЦ РАН и РСО-А, старший научный сотрудник лаб. теории операторов РОССИЯ, 362027, Владикавказ, ул. Маркуса, 22 E-mail: kulaevrch@mail.ru

The mixed problem for the parabolic equation set on geometrical graph, with boundary conditions containing a derivative on time is considered. The theorem of existence of the solution of the boundary problem giving representation of the solution in the form of circuit integral is established.

ON EXISTENCE OF THE SOLUTION OF PARABOLIC PROBLEM ON THE GRAPH WITH A BOUNDARY CONDITIONS CONTAINING DERIVATIVES ON TIME

Kulaev R. Ch.

Key words: differential equation on the graph, the circuit integral method.