CC BY
43
2015
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Математика и механика
№ 6(38)
УДК 517.968.25
DOI 10.17223/19988621/38/4
Д.Ю. Иванов
УСТОЙЧИВАЯ РАЗРЕШИМОСТЬ В ПРОСТРАНСТВАХ ДИФФЕРЕНЦИРУЕМЫХ ФУНКЦИЙ НЕКОТОРЫХ ДВУМЕРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ТЕПЛОПРОВОДНОСТИ С ОПЕРАТОРНО-ПОЛУГРУППОВЫМ ЯДРОМ
Исследуются возникающие в задачах теплопроводности двумерные граничные интегральные уравнения, операторные ядра которых выражаются через
пространственно-временную С0 -полугруппу. При условии дО е Ск +2 доказана устойчивая разрешимость интегральных уравнений в пространствах к раз непрерывно дифференцируемых на границе Sfl векторных функций со значениями в пространствах типа Соболева, определяемых степенями генератора С0 -полугруппы .
Ключевые слова: граничное интегральное уравнение, теплопроводность, существование, единственность,устойчивость.
Одним из методов, применяемых для аналитического решения достаточно широкого класса краевых задач нестационарной теплопроводности, является метод граничных интегральных уравнений (ГИУ) [1, с. 175]. На его основе для решения начально-краевых задач теплопроводности был разработан ряд численных методов, традиционно объединяемых под общим названием «методы граничных элементов» [2, с. 156; 3-8]. В связи с необходимостью строгого обоснования таких численных методов, связанного с вопросами аппроксимации и устойчивости, были проведены исследования, посвященные устойчивой разрешимости соответствующих ГИУ в пространствах дифференцируемых функций [5, 6, 9, 10].
В настоящей работе исследуются двумерные векторные ГИУ типа Фредгольма второго рода с операторным ядром, выраженным через С0 -полугруппу U(т). Как
показано в работе [11], такие ГИУ позволяют решить векторные краевые задачи первого, второго и третьего рода для линейных дифференциально-операторных уравнений Д2 и = B и в плоской ограниченной односвязной открытой области О+
или ее внешности О- = R2 \ О+ . Операторный коэффициент B является генератором С0-полугруппы U(т) в пространстве L2 = L2(IY хIT) (IY = [0,Y],
IT = [0,T]): Bf = lim т-1(f - U(т) f). В свою очередь, указанные краевые зада-
т^+0
чи суть возможные постановки начально-краевых задач теплопроводности на временном промежутке IT в однородном цилиндре О+ х IY или О- х IY с неоднородными граничными условиями первого, второго и третьего рода на боковой поверхности цилиндра, нулевыми граничными условиями первого, второго или третьего рода (в зависимости от оператора B) на основаниях цилиндра и нулевым
34
Д.Ю. Иванов
начальным условием. Двумерные векторные ГИУ не совпадают с обычными скалярными ГИУ типа Вольтерра - Фредгольма для таких задач теплопроводности, но представляют собой по существу прямые суммы скалярных ГИУ типа Воль-терра-Фредгольма в плоских областях D± , порождаемые спектральным разложением пространственной составляющей С0 -полугруппы.
Преимущество таких двумерных ГИУ по сравнению с обычными состоит в возможности экономного вычисления разрешающих их сеточных операторов в алгебре полиномов, образованных степенями полугруппового оператора. В работах [12, 13] построены соответствующие вычислительные схемы, исследованы их аппроксимирующие свойства и вопросы устойчивости. При обосновании этих положений возникает необходимость установления инвариантности пространств Ck (dD, Hn (IY х IT)) относительно прямых и обратных операторов ГИУ, а также ограниченности таких операторов в указанных пространствах (здесь Ck (dD, HB (IY х IT)) - пространство к раз непрерывно дифференцируемых на SQ векторных функций со значениями в H’B (IY х IT); SQ - граница области D+ ; HB (IY х IT) - пространство типа Соболева, определяемое n+1 степенью оператора B). Основную задачу настоящей статьи составляет доказательство наличия у таких операторов перечисленных свойств.
Хотя исследуемые здесь интегральные уравнения и отличаются сами по себе от традиционно изучаемых, приведем еще несколько отличий полученных здесь результатов от аналогичных результатов других авторов. Разрешимость ГИУ в пространствах дифференцируемых функций, насколько мне известно, исследовалась или на границах типа Сш [6, 9, 10], или липшицевых [5], причем в последнем случае, как и в [6, 9], лишь в пространствах, являющихся соболевскими по всем переменным. В работе [10] доказана устойчивая разрешимость ГИУ второго рода в пространствах типа Ск (IT, Cn (dD)) или Ск (IT, Hn (dD)) (Hn (dD) - пространство Соболева), причем в областях Q произвольной размерности, но в предположении неограниченно гладкой границы. Поэтому стоит отметить, что в настоящей работе устойчивая разрешимость ГИУ в пространстве Ск (dD, HB (IY х IT)) доказана при менее обременительном условии dD е Ск+2.
Работа [14] представляет собой более ранний этап настоящего исследования: результаты получены лишь в пространствах С (dD, HB (IY х IT)) и Ск (dD, L2);
достаточное условие dD е Ск+2, по существу, только объявлено, но отсутствует часть доказательства, где оно используется (здесь теорема 1); исследовались лишь первая и вторая краевые задачи.
В конце настоящей работы сделано указание на справедливость утверждений, аналогичных полученным, для стационарных и нестационарных задач теплопроводности в плоской области и стационарных задач в цилиндре, а также целого класса абстрактных краевых задач для уравнения Д2и = Bu .
Устойчивая разрешимость в пространствах дифференцируемых функций
35
Постановки задач и предварительные сведения
Пусть дО е С2. Рассмотрим четыре краевые задачи (i = 1,2):
а2Д2 u± = Bu± (х = (xj, х2) е О±), (1±)
u± = (х е дО), (2a)
дпи±-П u± = w± (х едО), (2b)
решения которых - функции u± (х) со значениями в L2, определенные на О± . Здесь w± (x) - функции со значениями в L2, заданные на SQ; п - нормаль к кривой дО в точке х, направленная внутрь области О+; Д2 =д2Х + д2 , дп -
сильные производные векторных функций; a > 0 - коэффициент температуропроводности, ц> 0 - коэффициент теплообмена на боковой поверхности цилиндра. Оператор B определен в пространстве L2 как B = Bt + By + p E на пересечении областей определения операторов Bt и By . Здесь E - тождественный оператор; оператор Bt: (Bt f) (y, t) = д{/(y, t), задан на абсолютно непрерывных по t функциях f (y, t) е L2, таких, что дг/ е L2 и f |t=0 = 0 при почти всех y е 1Y ; оператор By : (By f) (y, t) = -a2 д2yf (y, t), задан на абсолютно непрерывно дифференцируемых по y функциях f (y, t) е L2, таких, что д2yf е L2 и
(дyf-Я,о f) | y=о =^yf + Kf) I y=y = 0 (0 <V XY <ж) при почти всех t е 1T ; p > -ц1, где ц1 > 0 - наименьшее собственное значение оператора By (ц1 = 0 лишь при Х0 =XY = 0). Оператор B замкнут как сумма двух замкнутых операторов в пространстве L2, порождающих С0 -полугруппы и действующих вдоль различных переменных [15]. Оператор B порождает экспоненциально убывающую С0-полугруппу U(т): ||U(т)|| < exp[-(p + p.j )т], причем U(т) - нулевые операторы при т> T .
Будем считать, что если по условию значения векторной функции принадлежат банахову пространству, то предельные операции над этими значениями по умолчанию осуществляются в норме этого пространства. Обозначим через С (О')
и Ck (О') пространства непрерывных и k раз непрерывно дифференцируемых функций со значениями в L2, определенных на множестве О' с R2.
Определение 1. Решением уравнения (1 ±) будем называть функцию ы± (х) е С2 (О±) со значениями в D(B) (области определения оператора B), обращающую уравнение (1±) в истинное равенство.
Определение 2. Решением задачи |^±| будем называть функцию
и± е С(О±), являющуюся решением уравнения (1 ±) и удовлетворяющую гра-
36
Д.Ю. Иванов
ничному условию (2а). В случае задачи {р } будем требовать также выполнения
условия: и- ^0 при |х| ^-да.
II IIL2
Определение 3. Решением задачи {P2±} будем называть функцию
и± e C(Q±), являющуюся решением уравнения (1±) и имеющую с внутренней (внешней) стороны Ж правильную нормальную производную dinu:± (дпи± (х±|n) ^ дПи± (х) при равномерно относительно х edQ), опреде-
ляемую равенством (2b): дПщ = w± + nи± . В случае задачи {P2-} будем требо-
^ 0 при |х| ^ да
вать также выполнения условия x и
\\Vu
(|MlL2 =|дxj И||L2 +|дХ2 и|L2).
Пусть n1 и n2 - нормали к кривой 5^, проходящие через точки х' и х соответственно и направленные внутрь области Q+ ; дифференцирование д и dn
осуществляется по точкам Х и х соответственно. Зададим на множествах Q± векторные функции р± (х) (i = 1,2) со значениями в пространстве L2 с помощью криволинейных интегралов первого рода:
Р± (х) = j дnj K(r) v± (х') ds', p± (х) = J K(r) v± (х') ds',
dQ dQ
где r = \х - х'\, v± (х) - векторные функции со значениями в L2, заданные на dQ; K(r) (r > 0 ) - функция со значениями в пространстве ограниченных операторов, действующих в L2 , определяемая равенствами
K(r) f = (4п )-1 J т-1 exp[-rУ(4а2т)] U(т) f dт (f e L2 ).
IT
Согласно работе [11], при условиях dQ e C2 и v± e C(dQ), функции p± и p± являются векторными аналогами потенциалов двойного и простого слоев соответственно. Если dQe C2 и w± e C(dQ), то задачи {p1} однозначно разрешимы и
их решения представимы в виде соответствующих функций pi± с неизвестными v± e C(dQ), однозначно определяемыми ГИУ:
G±v±, + (-1) 2-1 v±+ Gtv±= w± ,
(G f )(х) = J dn1 K(r)f (х')ds',
dQ
(G f)(х) = J [dn2K(r)-nK(r)] f (х')ds'
()
dQ
(х e dQ).
Устойчивая разрешимость в пространствах дифференцируемых функций
37
Анализ интегральных операторов в пространствах дифференцируемых функций
Введем в рассмотрение параметрические уравнения кривой 5Q: x1 = x1 (s), x2 = x2 (s), где s - длина дуги, откладываемой от некоторой фиксированной точки в определенном направлении и заканчивающейся в точке x = (xl, x2). Функции x1 (s), x2 (s), периодические с периодом 2S (S - половина длины 5Q), осуществляют взаимнооднозначное отображение множества IS = (-S, S] на множество 3Q. Очевидно, что если функции xj(s), x2(s) принадлежат классу Ck (IS), т.е. имеют непрерывные производные на замкнутом множестве IS до порядка k включительно, то 9Q e Ck. Условимся далее писать 5Q e Ck, если функции x (s), x2 (s) принадлежат классу Ck (IS). Кроме того, пусть s и s' - значения параметра, соответствующие точкам x и x', и пусть ст = s' - s.
Введем в рассмотрение функции уm (s, s') (m = 0,1,2 ), заданные на множестве
IS х IS при s' Ф s равенствами у m = pm /ст2 , где
Р0 (s, s') = (4a2 ) r2 = (4a2 ) {x (s') - {(s)]2 + [x2 (s ') - x2 (s)]2 } , pi(s,s') = (8na2) r 3^r = (8na2) {-[x(s') -x1(s)]x2(s' ) + [x2(s’) - x2(s)]xj(s')} ,
p2(s,s') = (8na2) r dn2r = (na2) {-[x1(s) -x(s')]x2(s) + [x2(s)- x2(s')]x{ (s) } ,
а при s' = s равенствами
у0(s,s) = (4a2 ) , yj (s,s) = y2(s,s) = (na2) [(s) x2(s) -x"2(s) x'(s)].
Теорема 1. Пусть 5Q e C”+2 (n e Z+ = {0,1,...}). Тогда существуют непрерывные на множестве IS х IS производные 3'sdJs,уm (i = 0,n - j , j = 0,n , m = 0,1,2).
Лемма. Пусть I - замкнутый интервал на вещественной оси. Предположим, что некоторая вещественная функция f (z, Z) имеет в области I х I непрерывные
производные dZdjf (i = 0,n , j = 0,n’), причем n < n’ и djf |^=z = 0 при z e I, j = 0, q -1, где q = n’ - n . Тогда функция h(z, Z), заданная при Z Ф z равенством h(z, Z) - f/(Z- z)q , а при Z = z равенством h(z, z) = dqf lZ=z/q!, имеет в области I х I непрерывные производные ffzdjh при i = 0, n - j , j = 0, n .
Доказательство леммы. Зададим на множестве I х I функции gk { m (z, Z) так, что при Z Ф z
gk,l,m (Z, Z)
dmdZ+lf l.
k!
Z=z
dm dZ+l+1fl.
■ (k +1)!
Z=z
(Z-z ) +... +
dm dn'-1f |
Z=z
(n'-1 -1)!
(Z- z )n'-1-k-l +
_________1________
(n'-1 -1 )!(Z - z )k
Z
Js m sn ' f Iz= (z- t )n ,-1-1dt
z
(l = 0, n’ - k , k = q, n’, m = 0,n)
+
38
Д.Ю. Иванов
(если k +1 = n', то функция gklm (z, Z) определяется только последним интегральным слагаемым), а при Z = z gkj,m (Z, Z) = dm % +f |с= zjk !.Непосредственно проверяется существование непрерывных на I х I производных dzgklm и
gk l m и выполнение равенств
dzgk,l,m = gk,l,m+1 + kgk+1,l,m , dZgk,l,m = gk,l+1,m kgk+1,l,m
при k < n’, m < n и k < n’, l < n' соответственно. Отсюда легко видеть, что производные &z djgq00 представляют собой линейные комбинации функций gktm при k +1 < q + i + j < n', m < i. В силу формулы Тейлора с дополнительным членом в виде определенного интеграла [16, с. 146] имеем h = gq 0 0. Следовательно, если
i + j < n, то производные д\djh существуют и непрерывны. Лемма доказана.
Доказательство теоремы 1. Можно убедиться, что условия леммы выполняются, если f = pm и q = 2, при этом для функций р0 и р1 полагаем z = s', Z = s , а для функции р2 полагаем z = s, Z = s'. Тогда получаем справедливость утверждения леммы для функций h = ym . Теорема доказана.
Определим в пространстве C(dQ) норму: ||/||с(dQ) = sup ||/ (x)|| , что делает
это пространство банаховым. Введем также в рассмотрение банаховы пространства Ck(dQ) (k = 1,2,...), состоящие из функций / е C(dQ), имеющих непрерывные на множестве IS производные /(l): /(l)(s) = dl/(x(s))/ds1 (l = 1,k), с нормой |\/\\Ck (dQ) = max|| /(l )| |c (dQ). Будем считать, что C 0(dQ) = C (dQ).
Оператор A, отображающий банахово пространство S само в себя, условимся обозначать как A[S].
Теорема 2. Пусть dQ е Ck+2. Тогда операторы Gi |^Ck (dQ)J (k е Z+, i = 1,2) всюду определены и ограничены.
Доказательство. Введем в рассмотрение вещественные функции
Xm (s, ст) = ^m (s, s + ct) (m = 0,1,2 ). Пусть / е Ck (dQ). Запишем выражение Gi / в следующем виде:
S
(Gj/)(s)= j |ct|-“ Z(s,ct)/(x(s + ct))dct ,
-S
T
Zi (s, ct)/ = j та/2-1 zt (s, ct, t) U(t) / dт ,
0
z (s, ct, t) = X'+a 2 Xj (s, ct) exp [ (s, ст) X], z2(s,ct,t) = Xa/2 [(s,ct)X-in]exp[-Xo(s,ст)X], (4)
где a е (0,1) - некоторое фиксированное число; X = ст2/т, г] = 2 а2ц ; Zi - функ-
Устойчивая разрешимость в пространствах дифференцируемых функций
39
ции со значениями в пространстве ограниченных операторов, действующих в L2.
В силу теоремы 1 существуют непрерывные при (s, ст) е IS х IS производные
дsхm (l = 0,k). Отсюда с учетом неравенства х0 > 0 получаем существование непрерывных и ограниченных на множестве IS х IS х (0, да) производных dlszi (l = 0, k). В результате, принимая во внимание ограниченность операторов U(т) (т>0) в совокупности, приходим к существованию в операторной норме непрерывных на множестве IS х IS производных dls Zt (l = 0, k). Тогда, используя представление (4) и учитывая, что f е Ck (дО), имеем Gi f е Ck (дО) и оценки:
G f )(l)(s
£ C } |ст| ad1sZl (s, ст) ds 1 f (x(s + ст))dст
l'=0
L2
< 2(1 -a)-1 S1-a cu max_ f(l)(s)l £C\ < 2k+1 (1 -a) 0<l '<l,seI^ L2 l=0
- S1 c
где cit = max lid, Zi (s, ст)||, s е IS , Cl, = l !/[(l -1’ )!l'!] (l = 0, k). Из полу’ 0<l '<l ,(s,CT)eIs xI ^_11
ченных оценок вытекает ограниченность операторов Gt ]Ck (дО) ] . Теорема доказана.
Обозначим через HB пространство функций f е L2, таких, что Bm f е L2
(m = 1, n), с нормой
£
m=01|
вm
-i1/2
. В силу замкнутости оператора B
пространства HB банаховы. Введем в рассмотрение пространства Ck (дО)
(k = 0,1,..., n = 1,2,...), состоящие из элементов f е Ck(дО), таких, что f (x) е HB при x едО и Bmf е Ck (дО) (m = 1, n), с нормой
л ,д^ = max sup || f(l)(s)||нп . Пространства Ckn (дО) также банаховы. Легко
\\rk (дО) _ ААШЛ °иР
"Cn (дО) 0<l<k sеT
2
видеть, что если оператор A [Ck (дО) ] всюду определен и ограничен и коммутирует с операторами Bm (m = 1, n) на множестве CП (дО), то оператор A [Ckn (дО)] также всюду определен и ограничен: ||A||k n < ||A||k (здесь ||A||k n, ||A||k - нормы операторов A [c^ (дО)], A [Ck (дО)] соответственно).
Следствие 1. Пусть дОе Ck+2. Тогда операторы Gi [c^ (дО) ] (k е Z+, n е N, i = 1,2) всюду определены и ограничены.
Доказательство. Используя замкнутость оператора B и выполнение равенств U(т)B = BU(т) на множестве D(B), несложно убедиться в справедливости равенств Gt Bm = BmGt (m = 1, n) на множестве Ckn (дО). Поэтому с учетом теоремы 2 операторы Gi [Ck (дО)] всюду определены и ограничены. Утверждение доказано.
40
Д.Ю. Иванов
Теорема 3. Пусть 3Q e Ck+2, w± e Ck (3Q) (k e Z+, i = 1,2) и v± e L2, где L2 = L2 (3Q х IY x IT). Тогда v± e Ck (3Q).
Доказательство. Представим операторы Gi в виде G{s+ G"e, где
S S
(G'e f)(s) = I K'1S (s, s') f (x(s')) ds', (G"e f )(s) = I K"e (s, s') f (x(s')) ds',
-S - S
K't e (s, s') = Ki (s, s') Фе (ст), K”e (s, s') = Ki (s, s') [1 - Фе (ст)],
Ki(s,s') f = dniK(r) f = P(s,s')IT 2exp[-p0(s,s')/t] U(t) f dт,
It
K2(s,s') = dn2K(r)-nK(r) = I [ft(s,s')t 2-TIT 1 ]exp[-po(s,s')/t]U(t) f dt,
It
Фе (ст) - k раз непрерывно дифференцируемая вещественная функция: фе (ст) = 1 при |ст|<е/ 2, 0 <фе (ст) <1 при е/ 2 <|ст|<е, фе (ст) = 0 при |ст| > е (0<e<S ); K'iz (s, s’), K"e (s, s’), Ki (s, s’) - функции со значениями в пространстве ограниченных операторов, действующих в L2.
Так как 5Q е Ck+2, то на множестве IS х IS при s ' Ф s существуют непрерывные производные dlsp0 (l = 0,k + 2) и ^spl, (l = 0,k +1). Кроме того, р0 >0
при s' Ф s . Поэтому при s' Ф s в операторной норме существуют непрерывные производные dls Kt (s, s ') (l = 0, k +1). Следовательно, если функцию K"e (s, s ') доопределить при s’ = s нулевыми операторами, то на множестве IS х IS существуют непрерывные производные dlsK”e (s, s') (l = 0, k). Тогда, учитывая, что v± e IL2 и w± e Ck (5Q), имеем h± = +2(-1) (w± - G"ev±) e Ck (5Q).
Далее заметим, что согласно теореме 6 [11] функции dn K(r) ограничены на 9Q, а согласно теореме 3 [11] функция rа K(r) ограничена на 9Q при любом
а > 0. Кроме того, кривая 5Q e C2 не имеет точек самопересечения, следовательно, существует c0 = max _ (ст/r).
(s,s')eIS xIS
Зафиксируем некоторое число а e (0,1). Обозначим ci = max _ ||rа2 Ki , где
(s,s ')eIS xIS
функции rа2Ki (s, s' ) доопределены при s ' = s по непрерывности. Полагая е1-а <(1 -а)(8С0аc2S) , имеем неравенства:
S f S
ds <
К| Г < | l| ||ra'2 K'e (s, s ' )|| r-а ds'
- s V-s
S ( e 2 ^
< Са I l I ||rа/2K1 (s,s + ст)|| ст-аdст ds < 4 (1 - а)-1 с0аc2 Sе1-а < 2-1:
-S V-e
Устойчивая разрешимость в пространствах дифференцируемых функций
41
в силу которых к уравнениям (l + 2(-1)г G'ie) v± = h±, эквивалентным соответствующим уравнениям (3± ), применима теорема Банаха: функции v± представимы в виде рядов У"=0 j±n (j±n = (+2(-1)iGi'e^ ) ), сходящихся в норме J. Используя рассуждения предыдущей теоремы и учитывая, что
S
(,е f )))= { |СТГ“ Z (S,СТ) Фе (ст) f (*(s + СТ))СТ
и |фе (ст)| < 1, получаем существование непрерывных производных (+±п) (l = 0, к) на множестве IS и справедливость оценок:
(i)
([±2 g;,] n f )
<( *2 (1 -а)-1 ,1-а с,,)
< 2l+2 (1 -а)1 е1-а сп max_
v ' ,l 0<l<i,SGIS
([±2 Gi'f f )) V)||L s(2‘*2 (1 -а)- 1
max_
0<l '<l, seI?
Ck (3Q).
b
2
Тогда, налагая дополнительное условие е1-а < 2-к-2 (1 -а)с, к , получаем оценки:
(n )(l) (S)
< qn\ h к
J 1 II i ICk (3Q) L?
(q = 2k+2 (1 -а) 1 е1 а cik < 1, l = 0,к, s e IS ), вследствие которых ряды (±n) ’ (s) (l = 0, к) равномерно сходятся на IS в норме L2. Тогда в силу j±n e Ск (5Q) имеем v± e Ск (5Q). Теорема доказана.
Следствие 2. Пусть 5Qe Ск+2. Тогда операторы G± |^Ск (5Q)] (к e Z+, i = 1,2) ограниченно обратимы.
Доказательство. Согласно теореме 11 [11], уравнения (3±) имеют единственные решения v± e J при условии w± e J . Следовательно, в силу теоремы 3 уравнения (3± ) имеют единственные решения v± e Ск (5Q), если w± e Ск (5Q). С учетом теоремы 2 это означает, что операторы Gi± [Ск (5Q) ] являются биекциями, и поскольку они при этом замкнуты, то в силу теоремы о замкнутом графике обратные операторы (g,* ) 1 [Ск (5Q) ] ограничены. Утверждение доказано.
Следствие 3. Пусть 5Qe Ск+2. Тогда операторы G± ^С^ (SQ)] (к e Z+, n e N, i = 1,2) ограниченно обратимы.
Доказательство. Оператор B ограниченно обратим как порождающий экспо-нениально убывающую С0 -полугруппу [17, с. 136, 149]. С учетом коммутативности операторов G± и Bm (m = 1, n) на множестве Сn (SQ) операторы (g,* ) и
42
Д.Ю. Иванов
B т коммутируют на множестве Ck (dQ), следовательно, множество Ckn (dQ) инвариантно относительно (g±- . В силу следствия 1 и теоремы о замкнутом графике операторы (g±- |^Ck (dQ)] ограничены. Утверждение доказано.
Заключение
Следствия 2 и 3 являются обоснованием устойчивой разрешимости ГИУ (3± ) в пространствах СП (dQ) (к = 0,1,..., n = 0,1,..., C (dQ) = Ck(dQ)): при любой правой части w± e C<k (dQ) существует единственное решение v± e Ck (dQ), непрерывно зависящее от w± в норме Ck (dQ).
Класс задач, для которых справедливы полученные результаты, можно расширить. С учетом работы [11] для их получения к оператору B достаточно предъявить следующие требования: (A) оператор B порождает экспоненциально убывающую C0 -полугруппу и может быть расширен до оператора, допускающего
спектральное разложение и также порождающего экспоненциально убывающую C0 -полугруппу; (B) операторы Gi должны быть компактны. Таким требованиям,
кроме оператора B = Bt + By + p E , удовлетворяют также операторы B = p E и
B = Bt + p E (p > 0), в случае которых задачи {Pk } представляют собой соответственно стационарные и нестационарные, , задачи теплопроводности в плоской области Q± , а также оператор B = By + p E (p > -|Tj ), в случае которого
задачи {р1}
представляют собой стационарные задачи теплопроводности в ци-
линдре Q± х IY . Вместе с тем стоит отметить, что компактность операторов Gt использовалась в работе [11] лишь для доказательства существования обратного оператора (Gf)-1 [ L ] на основе теории Фредгольма. В работе [18] доказана ограниченная обратимость операторов Gt [Lx2 ], когда оператор B фактически
удовлетворяет одному условию (A); при этом получены аппроксимации в виде рядов по неотрицательным степеням полугруппового оператора U(H) (H > 0),
сходящиеся к оператору (/Г при H —— +0 в норме LX,. Поэтому представлен-
ные в настоящей статье и работе [11] результаты относительно задач
{P}
и со-
ответствующих им ГИУ могут быть распространены на достаточно большой класс абстрактных двумерных краевых задач для уравнений (1±), определяемый условием (A), внутри эллиптического случая [17, с. 304].
ЛИТЕРАТУРА
1. МазьяВ.Г. Граничные интегральные уравнения. Анализ-4 // Итоги науки и техники. Серия «Современные проблемы математики». Фундаментальные направления. Т. 27. М.: ВИНИТИ, 1988. С. 131-228.
2. Бреббия К., ТеллесЖ., Вроубел Л. Методы граничных элементов. М.: Мир, 1987. 524 с.
Устойчивая разрешимость в пространствах дифференцируемых функций
43
3. Onishi K. Convergence in the boundary element method for heat equation // Teaching for Robust Understanding of Mathematics. 1981. V. 17. P. 213-225.
4. Mcintire E.A. Jr. Boundary integral solutions for the heat equation // Mathematics of computation. 1986. V. 46. No. 173. P. 71-79.
5. Costabel M. Bounndary integral operators for the heat equation // Integral Equation Operator Theory. 1990. V. 13. No. 4. P. 498-552.
6. Noon P.J. The single layer heat potential and Galerkin boundary element methods for the heat equation. Ph. D. Thesis. University of Maryland. 1988.
7. Shirota K., Onishi K. A boundary element Galerkim method for the Dirichlet problem of the heat equation in non-smooth domain // Scientiae Mathematicae. 1998. V. 1. No. 1. P. 107-123.
8. Hongtao Y. On the convergence of boundary element methods for initial-Neumann problems for the heat equation // Mathematics of computation. 1999. V. 68. No. 226. P. 547-557.
9. ArnoldD.N., Noon P.J. Coercivity of the single layer heat potential // Journal of Computational Mathematics. 1989. V. 7. No. 2. P. 100-104.
10. Hongtao Y. A new analysis of Volterra-Fredholm boundary integral equations of second kind // Northeastern Mathematical Journal. 1997. V. 13. No. 3. P. 325-334.
11. Иванов Д.Ю. Решение двумерных краевых задач, соответствующих начально-краевым задачам диффузии на прямом цилиндре // Дифференц. уравнения. 2010. Т. 46. № 8.
С. 1094-1103.
12. Иванов Д.Ю. Экономичный метод вычисления операторов, разрешающих некоторые задачи теплопроводности в прямых цидиндрах // Актуальные проблемы гуманитарных и естественных наук. 2014. № 9. С. 16-32.
13. Иванов Д.Ю. Вычисление операторов, разрешающих задачи теплопроводности в прямых цилиндрах, с использованием полугрупповой симметрии // Известия Московского государственного технического университетата МАМИ. 2014. Т. 4. № 4(22). С. 26-38.
14. Иванов Д.Ю. Анализ двумерных граничных интегральных уравнений, определяющих решения задач теплопроводности в прямых цилиндрах // Перспективы науки. 2014. № 12(63). С. 103-109.
15. Иванов Д.Ю. Замкнутость сумм дифференциальных операторов, возникающих в задачах теплопроводности в пространствах L2 // Тр. Ин-та системного анализа РАН. Динамика неоднородных систем. М.: КомКнига, 2005. Вып. 9(1). С. 111-123.
16. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. Т.2. М.: ФИЗМАТЛИТ, 2001. 810 с.
17. Крейн С.Г. Линейные дифференциальные уравнения в банаховом пространстве. М.: Наука, 1967. 464 с.
18. Иванов Д.Ю. Решение в пространстве L2 интегрального уравнения, соответствующего задаче теплопроводности в однородном прямом цилиндре на временной полуоси // Актуальные проблемы гуманитарных и естественных наук. 2013. № 11-1. С. 20-25.
Статья поступила 26.09.2015 г.
Ivanov D. Yu. STABLE SOLVABILITY IN SPACES OF DIFFERENTIABLE FUNCTIONS OF SOME TWO-DIMENSIONAL INTEGRAL EQUATIONS OF HEAT CONDUCTION WITH AN OPERATOR-SEMIGROUP KERNEL
DOI 10.17223/19988621/38/4
In this paper, we study two-dimensional vector boundary Fredholm integral equations of the second kind with an operator kernel expressed in terms of a spatial-temporal C0 -semigroup. Such two-dimensional integral equations allow one to obtain solutions of vector boundary value problems of the first, second, and third kind for linear differential-operator equations Д2и = Bu in a
planar bounded simply connected domain Q+ or its exterior Q- = R2 \ Q+ . The operator coefficient B is a generator of the C0 -semigroup in space L2(IY x IT) (IY = [0,Y], IT = [0,T]). In
44
Д.Ю. Иванов
turn, these boundary value problems are possible formulations of initial boundary value problems of heat conduction on the time interval IT in a homogeneous cylinder Q+ x IY or x IY with
inhomogeneous boundary conditions of the first, second, and third kind on the lateral surface of the cylinder, zero boundary conditions of the first, second, or third kind (depending on the operator B) on the cylinder bases and zero initial conditions. The main result of this paper is as follows: under condition 5Q e Ck+2 , the spaces Ck(5Q, HnB (IY x IT)) are invariant with respect to direct and inverse operators of the integral equations, and such operators are bounded in these spaces. Here, Ck(5Q, HnB (IY x IT)) is the space of vector functions, k times continuously differentiable
on the border dO with values in the Sobolev type space HnB (IY x IT) defined by powers n + 1 of the operator B.
Keywords: boundary integral equation, heat conduction, existence, uniqueness, regularity.
IVANOVDmitrii Yurievich (Candidate of Physics and Mathematics, Moscow State Academy of Water Transport, Moscow, Russian Federation)
REFERENCES
1. Maz'ya V.G. Granichnye integral'nye uravneniya. Analiz-4. Itogi nauki i tekhniki. Seriya «Sovremennye problemy matematiki». Fundamental'nye napravleniya. Moskow, VINITI Publ., 1988, vol. 27, pp. 131-228. (in Russian)
2. Brebbia C.A., Telles J.C.F. Wrobel L.C. Boundary element techniques. Springer-Verlag,
1984. 464 p.
3. Onishi K. Convergence in the boundary element method for heat equation. Teaching for Robust Understanding of Mathematics, 1981, vol. 17, pp. 213-225.
4. Mcintire E.A. Jr. Boundary integral solutions for the heat equation. Mathematics of computation, 1986, vol. 46, no. 173, pp. 71-79.
5. Costabel M. Bounndary integral operators for the heat equation. Integral Equation Operator Theory, 1990, vol. 13, no. 4, pp. 498-552.
6. Noon P.J. The single layer heat potential and Galerkin boundary element methods for the heat equation. Ph. D. Thesis. University of Maryland, 1988.
7. Shirota K., Onishi K. A boundary element Galerkim method for the Dirichlet problem of the heat equation in non-smooth domain. Scientiae Mathematicae, 1998, vol. 1, no. 1, pp. 107-123.
8. Hongtao Y. On the convergence of boundary element methods for initial-Neumann problems for the heat equation. Mathematics of computation, 1999, vol. 68, no. 226, pp. 547-557.
9. Arnold D.N., Noon P.J. Coercivity of the single layer heat potential. Journal of Computational Mathematics, 1989, vol. 7, no. 2, pp. 100-104.
10. Hongtao Y. A new analysis of Volterra-Fredholm boundary integral equations of second kind. Northeastern Mathematical Journal, 1997, vol. 13, no. 3, pp. 325-334.
11. Ivanov D.Yu. Solution of two-dimensional boundary-value problems corresponding to initialboundary value problems of diffusion on a right cylinder. Differential equations, 2010, vol. 46, no. 8. P. 1104-1113.
12. Ivanov D.Yu. Ekonomichnyy metod vychisleniya operatorov, razreshayushchikh nekotorye zadachi teploprovodnosti v pryamykh tsidindrakh. Aktual'nye problemy gumanitarnykh i est-estvennykh nauk, 2014, no. 9, pp. 16-32. (in Russian)
13. Ivanov D.Yu. Vychislenie operatorov, razreshayushchikh zadachi teploprovodnosti v prya-mykh tsilindrakh, s ispol'zovaniem polugruppovoy simmetrii. Izvestiya Moskovskogo gosu-darstvennogo tekhnicheskogo universitetata MAMI, 2014, vol. 4, no. 4(22), pp. 26-38. (in Russian)
14. Ivanov D.Yu. Analiz dvumernykh granichnykh integral'nykh uravneniy, opredelyayushchikh resheniya zadach teploprovodnosti v pryamykh tsilindrakh. Perspektivy nauki, 2014, no. 12(63), pp. 103-109. (in Russian)
Устойчивая разрешимость в пространствах дифференцируемых функций
45
15. Ivanov D.Yu. Zamknutost' summ differentsial'nykh operatorov, voznikayushchikh v zada-chakh teploprovodnosti v prostranstvakh L2. Tr. In-ta sistemnogo analiza RAN. Dinamika neodnorodnykh sistem. Moskow, KomKniga Publ., 2005. Vyp. 9(1), pp. 111-123. (in Russian)
16. Fikhtengol'ts G.M. Kurs differentsial'nogo i integral'nogo ischisleniya, vol. 2. Moskow, Fiz-matlit Publ., 2001. 810 p. (in Russian)
17. Kreyn S.G. Lineynye differentsial'nye uravneniya v banakhovom prostranstve. Moskow, Nauka Publ., 1967. 464 p. (in Russian)
18. Ivanov D.Yu. Reshenie v prostranstve L2 integral'nogo uravneniya, sootvetstvuyushchego zadache teploprovodnosti v odnorodnom pryamom tsilindre na vremennoy poluosi. Aktu-al'nyeproblemy gumanitarnykh i estestvennykh nauk, 2013, no. 11-1, pp. 20-25. (in Russian)
