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

Вычисление операторов, разрешающих задачи теплопроводности в прямых цилиндрах, с использованием полугрупповой симметрии

к.ф.-м.н. доц. Иванов Д.Ю.

Университет машиностроения 8-916-219-22-64, ivanovdyu@yandex.ru

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

Ключевые слова: метод граничных элементов, метод разделения переменных, операторный метод, теплопроводность.

Введение

В настоящей работе рассматривается операторный метод численного решения начально-краевых задач для однородного уравнения теплопроводности с постоянными коэффициентами в цилиндрической области Qx IY на конечном временном промежутке IT (здесь

IT = (0, T], IY = (0, Y); Q = Q+ или Q = Q , Q+ - плоская открытая ограниченная односвяз-

ная область, Q— = R2 \ Q+ , R = ( —да, ). Задачи решаются с нулевыми начальными условиями, неоднородными граничными условия первого или второго рода на боковой поверхности цилиндра и нулевыми граничными условия первого, второго или третьего рода на основаниях цилиндра. Для численного решения таких задач применимы три известных подхода: метод граничных элементов (МГЭ), метод конечных элементов и метод конечных разностей. Для получения сеточного оператора, разрешающего задачу в фиксированной точке цилиндрической области, более предпочтительным с точки зрения экономии вычислений является МГЭ, в рамках которого приходится исключать промежуточные величины только в узлах на поверхности цилиндра, а не в некоторой прилегающей к ней трехмерной области, как в остальных двух методах. Заметим, что использование готовых разрешающих операторов дает ощутимую экономию при получении решения в небольшом количестве точек по сравнению, например, с методом Гаусса, применяемого в рамках МГЭ. Кроме того, наличие готовых разрешающих операторов необходимо для реализации различных методов численного решения обратных граничных задач теплопроводности, например, градиентных [1, с. 150].

Вычисление сеточного оператора, разрешающего трехмерное по пространственным координатам граничное интегральное уравнение (ГИУ) в рамках МГЭ, требует значительных затрат. В работе [2] предложен экономичный метод вычисления операторов, разрешающих рассматриваемые здесь задачи в цилиндре. Экономия достигается благодаря разделению вычислений по переменным x е 5Q (5Q - граница области Q ) и y е IY и проведению вычислений в алгебре полиномов, образованных степенями некоторого полугруппового оператора. При этом коэффициенты полиномов - операторы в пространстве функций, заданных на множестве 5Q, тогда как полугрупповой оператор действует в пространстве функций на

множестве 1У х 1Т. Численный метод строится на основе аналитического метода, являющегося комбинацией метода ГИУ и метода Фурье [3]. При этом ГИУ являются двумерными в том смысле, что определяемые ими векторные потенциалы заданы на плоском множестве О, так что сами ГИУ решаются на его границе дО . Коэффициенты полиномов получаются при решении соответствующих "плоских" задач на основе МГЭ. Заметим также, что кроме экономичности, преимуществом предложенного в [2] подхода является использование одномерных граничных элементов вместо двумерных. Наконец, полученные разрешающие операторы, имеющие вид полиномов, могут быть использованы для экономичного и устойчивого численного решения обратных граничных задач теплопроводности на основе канонической факторизации таких полиномов [4].

Метод [2] использует простейшую кусочно-постоянную аппроксимацию по параметру полугруппы и вдоль кривой дО. В настоящей работе исследуется возможность увеличения точности данного метода за счет использования вдоль названных множеств аппрокимаций более высокого порядка, использующих квадратичные функции. Для полученного алгоритма доказаны сходимость и устойчивость, установлены более высокие порядки аппроксимации относительно шагов дискретизации вдоль параметра полугруппы и кривой дО , сделан анализ экономичности метода, приведены результаты вычислительных экспериментов по сравнению точностей различных реализаций метода. Поскольку увеличение точности практически не влечет увеличения объема вычислений, в заключение работы сделан вывод о предпочтительности предлагаемой здесь реализации метода по сравнению с описанной в работе [2].

1. Постановка задачи и ее аналитическое решение на основе векторных потенциалов

Пусть дО - кривая класса С2. Рассмотрим четыре краевые задачи (/ = 1,2):

а2Аи±= Ви± (хе°±), и± ^^^, д„и± 1^ = ^, (1.1)

где: и± = и± (х) и w± = w± (х) (/ = 1,2) - векторные функции со значениями в пространстве Ь2 = Ь2 (1У х 1Т), заданные на множествах О и дО, соответственно (все вводимые здесь пространства функций считаем комплексными); п - нормаль к кривой дО в точке х, направленная внутрь области О+; А = д2 + д^ (непрерывность и диффе-

ренцируемость векторных функций предполагается здесь в норме пространства их значений); а > 0 - коэффициент температуропроводности. Замкнутый оператор В определен в пространстве Ь2 на множестве 0(В() 0(Ву) как В = В{ + Ву + р .

Здесь - замкнутый оператор в Ь2(1Т): (В( ¥)(/) = /'(/), заданный на абсолютно непрерывных на 1Т функциях /(¿) е Ь2(1Т) таких, что /(0) = 0; Ву - замкнутый оператор в Ь2( I,) : (Ву ¥)(у) = -а2/"(у) , заданный на абсолютно непрерывно дифференцируемых на 1У функциях /(у)еЬ2(17) таких, что /'(0)-/) = /'(У) + /У) = 0 (0<ку <¥); р + т > 0, где т - наименьшее собственное значение оператора Ву. Заметим, что в настоящей работе рассматриваются только линейные пространства и операторы.

В работе [3] доказана однозначная разрешимость задач (1.1) в классе С (О) С 2(О) при любых w± е С(дО) (С(О" ) и Ст (О" ) - пространства непрерывных и т раз непрерывно дифференцируемых векторных функций со значениями в Ь2). Решения имеют вид векторных потенциалов — криволинейных интегралов первого рода:

и± (х) = ОДу± = |дп1 К (г) у± (х') М , и± (х) = а2(х)у± = | К(г) у± (х') Ж", ( х еО ), (12а)

дО дО

где функция у± (х) находится из соответствующего ГИУ:

С? = + (-1)г2-1+С,, (0,Т)(х)= \ дпК(гЩх')Ж' (*ейП). (1 2с)

эп

Для уравнений (1.2Ь) доказана устойчивая однозначная разрешимость в пространствe Е2 = ¿2(ЭПх¡г х 1Т) для любой правой части (см. теорему 11 [3]). Здесь г = |х-х'|, у±(х) -векторная функция со значениями в Ь2, определенная на ЭП; п1 и п2 = п - нормали к кривой ЭП, проходящие через точки X их, соответственно, и направленные внутрь области П+; дифференцирование Эп и Э осуществляется по точкам х' и х, соответственно; К (г) ( г > 0) - ограниченные операторы в пространстве Ь2, определяемые равенствами:

К (г) Г = | а(г,т)и (т) Щ йт (Г е Ь2),

1т

где: а(г,т) = (4^т)-1ехр[-г2/(4а2т)], и(т) = и,(т)иру (т), ир (т) = ертиу(т).

Операторы и (г) образуют С0 -полугруппу, порождаемую оператором Б{: (и (г) Г) =

/(/ -т) при г< ?. Операторы иу (г) образуют С0-полугруппу, порождаемую оператором Ву

(иу(т)Щ)(у) = |Л(у,ут/(у')йу' (ГеЗД), т>0), (13)

ь

где: Л(у, у \т) = ^ .=1 ¥] (у) ¥] (у')ехр (-т,т), при этом т - собственные значения оператора

Ву (0 < т <т2 <. ., т1 = 0 лишь при Н0 = Л = 0), (у) - соответствующие собственные функции, нормированные в Ь2 (1Т ). Заметим, что \\и( (т)|| = 1 при т< Т, и (т) = 0 при т> Т и ||иу (т)|| < ехр (-цт) . Поэтому ||и(т)| < ехр [-(р + т)т]< 1, и(т) = 0 при т> Т. Таким образом, все рассматриваемые С0-полугруппы сжимающие, а и{ (т) и и (т) - нильпотентные.

Наряду с задачами в прямом цилиндре Пх 1Г здесь рассматриваются и соответствующие "плоские" задачи в области П. Последние допускают аналогичную формулировку на основании уравнений (1.1), где В = В( + р (р > 0) — операторы в пространстве Ь2 = Ь2(1Т) ,

обладающие сходными свойствами с операторами В = В( + Ву + р ( р > -т ). Решения "плоских" задач также однозначно определяются равенствами (1.2) в соответствующих пространствах, при этом полугруппу (т)ир(т) заменяем на полугруппу е-рти( (т) , также сжимающую и нильпотентную. Вследствие этого основные результаты, полученные ниже (теоремы 1 - 6, следствия 1 - 3), справедливы для задач (1.1) как в цилиндре, так и в плоской области.

Введем в пространстве С(ЭП) норму: ГЦ ,эп = тах|Г(х)|| , что делает это простран-

II НС(Эп) хеЭП ^2

ство банаховым. Условимся оператор А, отображающий банахово пространство В в банахово пространство С, обозначать как А (В ® С). Обозначим через Ск(ЭП) (к = 1,2,...) банаховы пространства, состоящие из элементов Г е С(ЭП), таких, что Г(х) е Б(В!) при х е ЭП

и ВЩ е С(ЭП) (I =1, к X с нормой Щсв (ЭП) =^к=ЛВ1 Щ|С(ЭП). Будем считать, что

С (ЭП) = С(ЭП). В работе [2] доказана теорема:

Теорема 1. Пусть ЭПе С2 . Тогда операторы О± (Свк (ЭП) ® Свк (ЭП)) (к = 0,1,...) огра-

ниченно обратимы.

Далее введем в рассмотрение параметрические уравнения кривой ЗО: х1 = х1(^), х2 = х2(¿) , где 5 - длина дуги, откладываемой от некоторой фиксированной точки в определенном направлении и заканчивающейся в точке х = (х1, х2). Функции х1(5), х2(5) периодические с периодом 28 (£ - половина длины ЗО) осуществляют взаимнооднозначное отображение множества 18 = (-£, £] на множество ЗО . Условимся обозначать через 5 и 5 ' значения параметра, соответствующие точкам х и х '.

Обозначим через Ск(ЗО) (к = 1,2,...) банаховы пространства, состоящие из функций ¥ е С(ЗО), имеющих равномерные на множестве ЗО производные ¥(/): (¥(/))(х(5)) = С11(х(5))/(/ = 1,к), с нормой \\f\C(ЗО) =Х/к=Л¥(/ЦС(ЗО)' Будем считать, что

С °(ЗО) = С (ЗО).

В работе [2] доказана следующая теорема.

Теорема 2. Пусть ЗО е Ск+2. Тогда операторы 0± (Ск(ЗО) ® Ск(ЗО)) ограниченно обратимы (к = 1,2,...).

Согласно теоремам 1 и 2 уравнения (1.2а) устойчиво однозначно разрешимы в банахо-

вых пространствах св (зо)=С (зо) п св (зо) с нормой ц/|(зо) /||С (зо) /|с(ЗО) при

любых /, п = 0,1,....

2. Построение и обоснование вычислительного алгоритма, оценка его эффективности

Так как ЗО - кривая класса С2, то она не имеет точек самопересечения, а функции (я,.?') = Зп 1пг"1 могут быть доопределены при я = я' до непрерывных на множестве 18х18.

На основании этого получаем при ^е 1° = 18 \{0} (а = 5 '-5) и т >0 оценки:

|Зпа(г,т)| £ с1 (а/т)2 exp [-с21/(4а2)], (2.1)

где: Я = а2/т, с1=(&7га2>) шах я')!, с2 = тт(г/|<т|).

С учетом оценок (2.1) операторы Ог могут быть представлены в виде:

вг ¥ = | А(т)и(т)¥сСт (¥ е С(ЗО)),

1Т

где: Д.(т) (т > 0) - ограниченные операторы (С(ЗО) ® С(ЗО)):

( Аг (т) ¥ )(х)= |ЗП а(г,т) Г(х') (¥ е С(ЗО), х еЗО).

ЗО

При т > 0 функции А1 (т) непрерывны в равномерной операторной топологии, причем справедливы оценки:

+СО

||4(г)||<сг-1/2, с = с, |г2ехр[-с2У/(4а2)]^. (2.2)

-со

Пусть ТУ/2 е N = {1,2,...} . Введем в рассмотрение операторы <Э1 (С(ЗО) —» С(дС2)):

. N-1 „

С, = {тп=пкт,кт =Т/Ю,

1Т п=О

где: им(т) - квадратичная функция на каждом интервале [т„,тп+2]: ик(тп+1) = 11(тп+])

(п = 0, N 2 -1, у = 0,2).

Операторы Gi ограничены в силу оценок (2.2) и < 3 . При условии ¥ е С"(дО.)

имеет место оценка:

С (<Ю)

< 2

-1

БЩ

С (5П)

я.

(2.3)

На основании оценок (2.2) и (2.2) получаем оценку:

сл-ол

С(дП)

<с

л/т7!^!

С(дП)

из которой вытекает следующее утверждение:

Теорема 3. Пусть ЭО е С2. Операторы в, (С"(сП) ^С(еО)) сходятся при N ® ¥ по операторной норме к соответствующим операторам 01 (СБ (ЭО) ® С(ЭО)) с порядком аппроксимации О (Н3).

Пусть Ь/2 е N. Разобьем кривую ЭО на дуги Т1 = [£, £+1 ] (£ = ¡Н, Н = 2£/Ь, I = 0,Ь -1, Ь е N). Введем в рассмотрение пространства УЬ сеточных функций ЧЬ со значе-

ниями Щ в пространстве Ь2, заданных в узлах £ . Определим в УЬ норму:

= шах| Щ

и зададим проекционные операторы РЬ (С(ЭО) ® УЬ): = РЬЧ, Щ = Ч(£), причем на каждом промежутке [ £21, £21+2] (I = 0, Ь/ 2 -1) функции (£) описываются квадратными трехчленами. Очевидно, что \\РГ\\ <3.

Введем в пространстве Уь ограниченные операторы (УЬ—>УЬ):

ы-1

и=0

Заметим, что операторы ^ Л имеют вид скалярных квадратных матриц порядка Ь . При ст е справедливы оценки:

I ои

||ЭП а(г,т) с1т <с , с =^|ехр -с2АД4а2) ¿//I,

а при Ч е С3(ЭО) имеют место оценки:

Р Ч - Щ

С (ЭО)

< 2-1

КЗ)

С (<Ю)

НЗ.

(2.4)

(2.5)

Определение 1. Пусть АЬ - ограниченные операторы в соответствующих пространствах УЬ . Будем говорить, что операторы АЬ сходятся при Ь ® ¥ к ограниченному оператору А (Ск (ЭО) ® С(ЭО)) по операторной норме, если ||АЬР Ч - РЬА Щ\у ® 0 при Ь ®¥ равномерно в шаре Щ |\Ск ) < 1.

На основании оценок (2.4), (2.5) и ||с7^(г)|| < 3 получаем оценку:

бЛ-вЛ <£с||г(3)|| /г3,

1 1 с(дп) II Нс(ап) Л из которой вытекает следующее утверждение:

Теорема 4. Пусть дС1 е С2. Тогда операторы сходятся при по операторной

норме к соответствующим операторам Сг1 (С3(ЭСХ) ~^С(дС2)) равномерно по N с порядком

аппроксимации О (Н£).

Введем в рассмотрение ограниченные операторы = +(—I)12-1 + (5, и

Ь

2

I

т

0

0

= +(-1)' 2 1 + о (Уь —> Уь). Поскольку ||(5? 0| < 6сфт~, что при условии:

12с^<1, (2.6)

в силу теоремы Банаха оператор (5~0 ограниченно обратим, следовательно, ограниченно обратим оператор и имеют место равенства:

И , (п = \,Ы-\). (2.7)

п=0 т=1

Множество значений N е N, обеспечивающих выполнение условия (2.6), обозначим через N '. Оценка (2.6) аналогична оценке (4.3) в работе [2], поэтому справедливо следующее утверждение.

Теорема 5. Пусть дС1еС2. Тогда операторы (X е IV, ТУ е №) ограничены в со-

вокупности.

Определение 2. Будем говорить, что ограниченные операторы А^) (УЬ ® УЬ) сходятся при Ь, N к ограниченному оператору А (С^п (ЗО) ® С(ЗО)) по операторной норме, ес-

ли

АЬN)Рь ¥ - РьА ¥

® 0 при Ь, N ® ¥ равномерно в шаре £ 1.

11 11С/ п ( ЗО )

Из теорем 1-5 вытекает следующее утверждение.

Следствие 1. Пусть дС1 е С5. Тогда операторы ) сходятся при Ь, N —> ос

по операторной норме к соответствующим операторам (О~) (С3В3(ЗО) ® С(ЗО)) с порядком аппроксимации О (И3 + И3).

Пусть х еО. Тогда функция а(г, т) непрерывна по (х, х ', т) еОх ЗО х [0, ¥) (а(г ,т = 0) = 0). При фиксированном х еО введем в рассмотрение ограниченные операторы О1(х) (С(ЗО) ® Ь2) (I = 1,2) (т >0):

О1 (х) ¥ = | Аг (х, т) и(т) ¥ сСт ,

1т

(А1(х,т)¥)(х)= |ЗП1 а(г,т)¥(х ')сСз', (А2(х,т)¥)(х)= | а(г,т)¥(х ')'.

ЗО ЗО

Операторы О. (х) (х еО') ограничены в совокупности на любом замкнутом подмножестве О 'сО в силу непрерывности функции А1 (х,т) по (х,т)еЗОх[0,¥) в равномерной операторной топологии (А\к) (х, т) = 0).

Введем в рассмотрение ограниченные операторы Сг(х) (С(ЗГ2)

1Т т=0

где: N = МЫ, Ме ]Ч, тт=ткт9 \ = йг/М; - квадратичная функция на каждом ин-

тервале [тт,тт+2]: О^т^Щт^) (т = 0,Й/2-1= ОД). Пусть ЗОе С2. Тогда интеграл I Зп 1пгСз' ограничен в области О [5, с. 361].

J ЗОI .

Пусть тах£^|Зп 1пг-1 ¿¿у' и с2 = (4;г) 12£. При хеО и т> 0 справедливы

оценки:

Ь

(, (х, т)I < Спт~], где с" = с. max Л2~г ехр[- Л/(4я2)].

Теорема 6. Пусть 5П е С2. Тогда операторы (х) (—>Ь2) сходятся при N —>оо и фиксированном М по операторной норме к соответствующим операторам 01 (х) (

—>1^) равномерно по хеО с порядком аппроксимации .

Данное утверждение вытекает из следующих оценок:

1

< 2~1cn P?3f hi 1п(Г//гг),

' II Нс(Ш) г v ' г'

C (ЭП)

j4(x,r)[c/(r)-c7^(r)]f£/r

<£Q|lB3f| 7z2.

' II IIC(SQ) г

с (an)

Ко второй оценке приходим в результате следующих выкладок:

т h

||U(t)f-ий(г)f|| < 2-1 JU(t)B3f (t-t)2dt + 2-1 (r2/h2)JU(t)B3f(h-t)2dt

2 h h 4-1 (t2/h2) J U(t) B3f (2h -1)2 dt - (t/h) JU(t) B3f (h -1)2 dt

+

2 h

4-1 (t/ h ) J U (t) B3f (2h -1 )2 dt

< 2 B 3f th;5

с (an)

осуществленных при условии t e [0, ht ].

Пусть L/2 e N. Введем в рассмотрение ограниченные операторы G}(х) (Vl^L2)\

G,(x)fL = { 4(х.t)Un(t)Pf f ¿r = XGim(x)U(fm)fr .

j m=0

iT

Заметим, что операторы Gim(x) суть скалярные матрицы размерности 1 xL.

Далее, пусть dWd - кривая в области W, параллельная границе ЭП и находящаяся от нее на расстоянии d [5, с. 263]. Область между dQd и ЭП обозначим через Wd, а через Q" обозначим ее дополнение в области W. Если х e and, то r > d. Если х e Wd, то г > 2~V0 > 2_1с|сг0|, при этом г0 = |х'-х0|, сг0 = s' — s0. Здесь х0 е 5Q - основание перпендикуляра, опущенного из х на ЭП; s0 - значение параметра я, соответствующее точке х0. Имеем следующие оценки:

Т S _т

J J|a(r,t)|dtdsf < (4p)-1 matx {l/4exp[-l/(4a2)]} J da/J2-1 c\a\ Jt~3/4dt (х eWd),

8Q 0 Ле[0,оо) ^ o

T

| ^\a{r,r)\drdsf <c2 max r"1 exp -d2/[Aa2r^ (xeQ"),

an о

т

11 max Л,exp4a2) (x e Q).

(2.8)

an о

По аналогии с определением 1 будем понимать сходимость по операторной норме ограниченных операторов АЬ (х) (Уь ® Ь2) к соответствующим ограниченным операторам

А(х) (Ск(д€1)^Ь2). На основании оценок (2.5), (2.8) и <3 получаем следующее

утверждение:

h

0

0

С

0

0

0

Теорема 7. Пусть дС1 е С2. Тогда операторы Сг{ (х) сходятся при Ь —» оо по операторной норме к соответствующим операторам СДх) (С3(дС1) —>С(дС2)) равномерно по N и х е О (М фиксировано) с порядком аппроксимации О(И3) •

Введем обозначения: /?~(х) = , = С1{х){о~) . По аналогии с опреде-

лением 2 будем понимать сходимость по операторной норме ограниченных операторов А[м)(х) (¥ь ® Ь2) к соответствующим ограниченным операторам А(х) (С^ (ЭО) ® Ь2) при ->оо.

В силу оценок (2.8), ||с7^(г)||<3 и ||р||<3 операторы <5(х) (¿е!Ч, хеО) при

фиксированном М ограничены в сокупности. Поэтому на основании теорем 6, 7 и следствия 1 делаем вывод о равномерной сходимости операторов К~ (х) в области .

Следствие 2. Пусть дС1 е С5. Тогда операторы Д±(х) сходятся при —> ос

и фиксированном М по операторной норме к соответствующим операторам (х)

(С3В3(ЭО) ® Ь2) равномерно по х е О с порядком аппроксимации О(Н2т + И3).

Операторы Я~(х) имеют вид:

т=0

где: Д т(х) - матрицы размерности \хЬ, вычисляемые с помощью уравнений:

п ___

= (к = 0,М-1, т = 0,М-1). (2.9)

к=О

При вычислении матриц <5,. п, Сг1т (х) интегрирование по г осуществляется точно, а по л' с помощью элементарных формул Гаусса с четным количеством узлов. Учитывая ограниченность в совокупности операторов

О (х) (х еО') на любом замкнутом подмножестве О'сО несложно получить следующее утверждение.

Следствие 3. Пусть ЭОеС5 и О 'с^ - замкнутое подмножество. Тогда операторы Щ(х) сходятся при —> ос и фиксированном М по операторной норме к соот-

ветствующим операторам Е± (х) (С3В3 (ЭО) ® Ь2) равномерно по х еО ' с порядком аппроксимации О (И + Из).

В силу следствий 2 и 3 точность решения задачи (1.1), полученного с помощью оператора (х), может снижаться при приближении к границе области . С этой целью введен

параметр М, который позволяет в силу тех же утверждений скомпенсировать уменьшение точности вычисления при приближении точки х к границе дС2, не увеличивая объема вычислений, необходимых для получения оператора СДх) (последнее обстоятельство более подробно обсуждается в конце параграфа).

Операторы (х) являются дискретными лишь по координатной переменной я. Для завершения построения вычислительного алгоритма остается осуществить их дискретизацию по координатным переменным ^ и у в случае В = В1 + Ву + р или только по / в случае

В = В, + р.

Рассмотрим вначале случай В = В( + В + р , Ь2 = Ь2(1У х 1Т) . Введем в рассмотрение

пространства Ст'" (т,п = 0,1,...), состоящие из функций /(х,у,t), имеющих непрерывные на множестве дОх 17 х 1Т производные д/(х,у,£) (3 = 0,т ) и д)/(х,у,^ (к = 0,п). В про-

странстве C0'0 введем норму: ||f||с0,0 = max |f (x,y,t)|, и будем считать, что

(x, y,t Iy xIj,

с0,0 с C(ao).

Пусть yj = jhy ( hy = y/J, j = 0, J , J/2 e N). В пространстве Cопределим проекционные операторы PJ : fJ = PJ f, при этом fJ (x, y,, t) = f (x, y,, t) , и на каждом промежутке

[У23, У23+2] ( 3 = 0, J|2 -1) функции ^ (х, у, t) при фиксированных (х, {) описываются квадратными трехчленами. Кроме того, пусть = гп = п\ (К = Нт, п = 0, N, Ы/2 е N ). Аналогично с помощью кусочно-квадратичной интерполяции вдоль переменной I определим в пространстве Соператоры РЫ: = РЫГ (/Ы(х,у,tn) = /(х,у,tn)). Операторы PJ и РЫ ограничены: ЦЛ £ 3, \\РЫ\\ < 3 . Если Г е С3'3, то имеют место оценки:

11« Г - 1с(дО) (Ц/ У, 4с0.0 К + 1/ У, ОЦ^ К ) . (2.10)

При фиксированном Л/eN операторы Р~(х) (£е1Ч, ограничены в

совокупности в силу такой же ограниченности операторов ¿гДх) (ДЛ^е]Ч,.хеГ2) и

(Ь е N, Nе N ', х еО). Операторы РЬ, , РЫ (Ь, J,Nе N) также ограничены в совокупности. Поэтому на основании следствий 2, 3 и оценки (2.10) имеем следующее утверждение.

Следствие 4. Пусть дО е С5 и О' с О - замкнутое подмножество. Если w± е С£3 (дО) п С3'3, то функции й± (х) = Кт(х) РЬ Р,Р^± ( N е N ') сходятся при Ь, J, N и фиксированном М к решению задачи (1.1) в норме пространства Ь2 равномерно по х еО' с порядком аппроксимации О(К + К + Ку + К) и равномерно по хеО с порядком аппроксимации О (К + К3 + Ку + К ) . Если же w± е С0 0, то существует число с > 0, не зависящее от

Ь, 3 е N и ТУ е IV, такое, что шах й. (х) <с т, I .

х-со II ' Нд II Не00

Для сравнения: аналогичные функции йг(лг), полученные в работе [2], сходятся равномерно по х еО' с порядком аппроксимации О (к32 + И21п К + Ку + К ) и равномерно по х еО с порядком аппроксимации О (Нт + К + Ку + К ) .

Интегральные выражения (1.3) позволяют вычислить полугрупповые операторы Ту (г) (г >0) на интерполянтах PJГ в виде конечных сумм:

(Ту (г)Р^ Г ) (х, у, t) = £ (г, у)/(х, у, t), (2.11)

3=0

причем значения функций г]з (г, у) вычисляются точно. Действие тождественного оператора Ту (0) на интерполянтах также описывается равенствами (2.11), где

Ъ3+к (г у) = (рк ((у - у23 У Ку ) , при этом у е [у23 , у23+2 К 3 = ^ J|2 -1,

^к (х) = П (х - к')/ П (к - к) ( х е [0,2]), к = 0,2.

к'=0,2(к'Ак) / к'=0,2(к'^к)

Окончательно равенство Й7(х) = может быть записано в следующем

скалярном виде:

Ь 3 N¡+2

й; (х, и у) = Е Е Е ви,,,п (*> У> 0 < (*/ > У у > ) • (2.12а)

I=1 ]=0 п=0

Здесь Щ/2 - целая часть дроби I / (2^ ); коэффициенты 0~1 ]к (х, у, I) вычисляются с помощью формул:

Мл

ви,у,о (Х> У>'0 = Е (Х' (а) А ) >

т=0

Кп+I"1 _

т=Кп

_ _

£=0,1 т=Кп_к

в?ш+2(х>У>0 = Е М*иАх>У)<Ро А) • (2.12Ь)

т=М,-2М+\

В свою очередь, здесь //;/;;;,(х.у) = г.щ,(х)<? /,г'";;Дгт. у). где /;.„,, (х) - элементы матриц (х); ЛГ, - целая часть дроби ///?г, ]У,0 = ТУ, -МЛГ,; Кп = 2М(п-1) + +1 (и = 1,^/2 + 1

Заметим, что увеличение параметра М не приводит к измельчению временной сетки, а лишь к увеличению объема вычислений разрешающего оператора, определяемого коэффициентами ]п (х, у, ^). Заметим также, что хотя в данной реализации шаг по времени Н1 совпадает с шагом по параметру полугруппы НТ (но не совпадает с Ьт!), но ничто не мешает сделать его не зависящим от Нт. Если, например, /гí//íг=M1GN и кт/кт =М2 е]Ч, то будут иметь место формулы (2.12), где М = М1М2.

Итак, формула (2.12а) позволяет вычислить приближенное решение й~(х^,у) задачи (1.1) в цилиндре в любой точке (х,у,еОх 1у х 1Т. При наличии готовых коэффициентов ]п (х, У, ^) для этого требуется »LJN умножений и делений (УД). Коэффициенты .п (х, у, вычисляются с помощью коэффициентов г*т1{х) и у), которые получаем

независимо друг от друга, при этом основные затраты 0( L3 N2) УД приходятся на вычисление элементов обратной матрицы по формуле (2.7), т.е. получение сеточного оператора, разрешающего двумерное ГИУ. Так, например, непосредственно на вычисление ] п(х>У>0 по формулам (2.12Ь) требуется 0(ЫЫ) УД, на ^/(х) по формулам (2.9) -

0(1}И2М) УД. Если решение й~ (х, у) получаем при Ь0, */0, значениях переменных х, у, /, соответственно, то затраты на вычисление остаются прежними, тогда как,

например, объем вычислений для й~{х^,у) и .п (х, у, возрастает в раз, для

КтМ) "В К Ра3-

Таким образом, в рамках данной задачи используется возможность разделения вычислений по переменным х и у. Это позволяет значительно сократить расходы на получение разрешающих ее сеточных операторов, которые, в свою очередь, обеспечивают выигрыш при получении решения в случае L0 J0N << LJN . Действительно, если не использовать раз-

Серия «Естественные науки»

деления переменных х и у, то для нахождения сеточного оператора, разрешающего трехмерное ГИУ на полной поверхности цилиндра ЭО х У и равномерной временной сетке, требуется О(И:33 N2) УД. Если же решать трехмерное ГИУ методом Гаусса, то каждый раз при изменении значений граничной функции мт (х, у, /) необходимо О(1}J2N2) УД, поэтому при условии L0J0N << LJN использование готовых разрешающих операторов более эффективно.

Реализации схем, описанных в работе [2] и в настоящей статье, отличаются лишь на этапе вычислений матриц (7. п, (7 т (х), требующих, соответственно, ()(/?М ), О(ЬЫМ) УД.

Поскольку LN << LN2 и LNM << LN, то можно считать, что в целом объем вычислений в обоих случаях одинаков.

В случае В = В( + р, L2 = L2(IT ) по аналогии получаем утверждение:

Следствие 5. Пусть дС2еС5. Если е С^3(дС1)г)С3, то функции щ(х) = Кг(х)РьРлг\\т ( N е К') сходятся при L, N ® да и фиксированном М к решению задачи (1.1) в норме пространства L2 равномерно по х е О' с порядком аппроксимации О (И + И + И3) и равномерно по х е О с порядком аппроксимации О(И + И3 + И ). Если же w± е С0, то существует число

с > 0, не зависящее от Ь е N и N е 14', такое, что шах й7(х) < с рут „.

ХС_П II ' Нд II ' Пс°

Для сравнения: аналогичные функции йг (х), полученные в работе [2], сходятся равномерно по х еО' с порядком аппроксимации О (И32 + И21п И + И ) и равномерно по х еО с порядком аппроксимации О (Ит + И + И] ).

Удаляя переменную у и индекс ], получаем формулы, аналогичные (2.12), где ,ЛХ) = ^±тАх)е~РТ'" ■ Заметим, что при М = 1 формулы (2.12) определяют приближенные решения задач (1.1) в плоской области О, полученные в рамках МГЭ в непрямой формулировке на основе кусочно-квадратичной аппроксимации граничной функции м, (х, /) [6, с. 168].

3. Вычислительные эксперименты

Рассмотрим численное решение внутренних трехмерных и "плоских" задач (1.1) в случае, когда граница ЭО представляет собой окружность радиуса Я = 1, при этом Т = У = 1, а = 1, р = 0. Решения "плоских" задач вычисляются на окружностях ЭWd, отстоящих от границы дС1 на расстоянии ¿7, в узлах (точки 5, получаются из в результате подобного отображения ЭО на ЭWd). Решения трехмерных задач вычисляются на круговых цилиндрических поверхностях оОс/ /1, в узлах ). Таким образом, /_,, = Л, ,/п =./ , Лг0 = N . Мелкость пространственно-временной сетки определяется значениями Ь=Ы = 30, ./ = 16. При вычислении матриц (7. я, С/ ш (х) интегрирование по 8' осуществляется с помощью 12-

точечной формулы Гаусса. Граничная функция на цилиндрической поверхности задается формулой м+(р,у,/) = /(у)/(¿)Бтр (р - полярный угол), а на границе круга - формулой м+(р,/) = /(Оэтр, при этом /(х) = 16х2 (1 -х)2. "Точные" решения и, (, = 1,2) "плоских" и трехмерных задач (1.1), заданные в узлах (£,,?„) и соответственно, находятся с

помощью функций Грина, причем интегрирование функции Грина по временной переменной

на интервале [0; 5х10-4] осуществляется численно с помощью 12-точечной формулы Гаусса, а все остальные интегралы от функции Грина вычисляются с помощью точного почленного интегрирования ее как ряда.

Приближенные решения "плоских" и трехмерных задач (1.1), полученные с помощью формулы (2.12) в точках (£,,?„) и , соответственно, обозначим через й1.. Близость

сеточных функций й* и й;+ определяется величиной относительного среднеквадратичного отклонения: ёи{ = ||Ли;||х 100%, где Ли. -П;+, ||-|| - среднеквадратичная норма. Будем обозначать через <377д/ и ¿ш(А/ значения ¿¡¡г для "плоских" и трехмерных задач, соответственно, при заданном значении М .

Как уже отмечалось, решения йг "плоской" задачи при М = 1 можно считать полученными в рамках классического "двумерного" МГЭ на основе кусочно-квадратичной аппроксимации. Можно предположить, что его точность не ниже точности соответствующего "трехмерного" МГЭ, примененного к задачам (1.1) в цилиндре Ох 17 и реализованного на Ь х У х N граничных элементов (при равных в обоих случаях Ь и N ). Поэтому будем сравнивать точность предлагаемого метода с точностью Зй1. Кроме того, в приведенной ниже

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

Таблица

Относительные среднеквадратичные отклонения Зи1 (в процентах)

М "плоские" к0 = кг = да И0 = кг = 0 И0 = Иг = 1

5их Зи2 Зи1 Зи2 Зи1 Зи2 Зи1 Зи2

1 0.013 0.0067 0.54 0.22 0.31 0.059 0.39 0.090

(0.21) (0.28) (0.27) (0.63) (0.29) (0.28) (0.31) (0.25)

0.9 2 0.0028 0.0038 0.028 0.014 0.016 0.0029 0.017 0.0027

(0.17) (0.24) (0.26) (0.24) (0.23) (0.24) (0.21) (0.20)

4 0.0035 0.0042 0.022 0.014 0.0087 0.0027 0.0081 0.0031

(0.16) (0.22) (0.25) (0.23) (0.22) (0.23) (0.20) (0.20)

1 0.012 0.0069 0.22 0.29 0.25 0.11 0.29 0.16

(0.23) (0.27) (0.41) (0.53) (0.27) (0.29) (0.32) (0.27)

0.5 2 0.012 0.0069 0.13 0.043 0.072 0.0078 0.086 0.011

(0.17) (0.24) (0.28) (0.25) (0.18) (0.25) (0.21) (0.21)

4 0.0051 0.0042 0.026 0.015 0.22 0.0035 0.22 0.0044

(0.16) (0.22) (0.23) (0.22) (0.16) (0.23) (0.20) (0.20)

1 0.26 0.093 2.31 0.98 1.91 0.35 2.22 0.50

(0.81) (0.44) (5.31) (2.77) (7.19) (1.85) (7.31) (2.25)

0.1 2 0.018 0.019 0.83 0.23 0.44 0.077 0.58 0.11

(0.34) (0.26) (2.52) (0.91) (3.13) (0.61) (3.28) (0.71)

4 0.020 0.0056 0.23 0.049 0.095 0.019 0.13 0.028

(0.13) (0.24) (0.77) (0.27) (0.74) (0.30) (0.88) (0.29)

1 0.29 0.21 1.67 1.32 1.28 0.61 1.44 0.94

(0.63) (0.41) (3.84) (3.52) (5.35) (2.34) (5.76) (2.96)

0.05 2 0.29 0.085 0.80 0.37 0.48 0.12 0.72 0.27

(0.33) (0.35) (2.22) (111) (2.87) (0.84) (3.02) (1.12)

4 0.29 0.050 0.44 0.14 0.35 0.050 0.43 0.18

(0.15) (0.27) (1.08) (0.42) (1.24) (0.40) (1.36) (0.56)

Согласно данным, приведенным в таблице, точность 8и][ меньше 8и-!, но

Серия «Естественные науки»

ôïtf1 « StL-"2 уже при небольших значениях Ml/M2 (Sir « 8Ti ). Напомним, что увеличение M влечет увеличение объема вычисления сеточных операторов, но не измельчение временной сетки. Это позволяет быстро вычислять высокоточные разрешающие операторы, допускающие также экономную реализацию. Действительно, вычисление аналогичных сеточных операторов на основе классического "трехмерного" МГЭ требует в J3/M » 500 раз больших затрат, поэтому преимущества предлагаемого метода, основанного на разделении переменных, очевидны. Заметим также, что согласно данным таблицы схема, описанная в настоящей работе и основанная на квадратичной интерполяции вдоль параметра полугруппы t и вдоль кривой 5Q, дает ощутимое увеличение точности по сравнению со схемой [2], хотя при этом объем вычислений, необходимых для получения разрешающих операторов, практически не возрастает.

Заключение

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

Литература

1. Алифанов О.М., Артюхин Е.А., Румянцев С.В. Экстремальные методы решения некорректных задач и их приложение к обратным задачам теплообмена. М.: Наука, 1988. 287 с.

2. Иванов Д.Ю. Экономичный метод вычисления операторов, разрешающих некоторые задачи теплопроводности в прямых цилиндрах // Актуальные проблемы гуманитарных и естественных наук. 2014. № 9(68). С. 16 - 32.

3. Иванов Д.Ю. Решение двумерных краевых задач, соответствующих начально-краевым задачам на прямом цилиндре // Дифференц. уравнения. 2010. Т. 46. № 8. С. 1094 - 1103.

4. Иванов Д.Ю. Обоснование одного алгоритма численного решения обратных граничных задач теплопроводности, построенного с учетом полугрупповой симметрии таких задач // Ж. вычисл. матем. и матем. физ. 1998. Т. 38. № 12. С. 2028 - 2042.

5. Михлин С.Г. Курс математической физики. М.: Наука, 1968. - 575 с.

6. Бреббия К., Теллес Ж., Вроубел Л. Методы граничных элементов. М.: Мир, 1987. - 524 с.