Аппроксимационные свойства метода конечных суперэлементов Федоренко

Вычислительные технологии

Том 13, Специальный выпуск 4, 2008

Аппрокеимационные свойства метода конечных суперэлементов Федоренко*

*

С. А. Лазарева Московский государственный технический университет им. Н.Э. Баумана, Россия e-mail: Lazarevas@gmail.com

This paper addresses a development of the Fedorenko finite super element method (FSEM). A priori error estimates are obtained; saturability of the method in Sobolev spaces is determined. The conditions leading to the derivatives convergence are indicated, corresponding error estimates are given. Regularity and accuracy of the numerical solution on the Sobolev class of functions are shown. The investigation deals with the Dirichlet problem, but can be extended over a general class of linear elliptic problems.

Введение

Метод конечных еуперэлементов (МКСЭ) впервые предложен в работах Федоренко и его коллег [1, 2], Он входит в класс численных методов, основанных на декомпозиции области в сочетании с выбором особой аппроксимации решения. Приближенное решение МКСЭ заведомо содержит в себе некоторые из свойств, присущих рассматриваемой задаче. Это позволяет проводить расчеты ряда "сложных" вычислительных задач, В работах [1-7] эффективность МКСЭ подтверждена примерами решения разнообразных физических проблем. Для их авторов преимущественный интерес представляют задачи, которые характеризуются наличием множества резких особенностей, проявляющихся на малых пространственных областях. Такие особенности могут представлять собой как "сингулярности" самого решения, так и резкие неоднородности области расчета.

Рассмотрим простую задачу Дирихле для дифференциального уравнения Лапласа в двумерной области П С Я2. Расчетная область представляет собой квадрат с исключенными из него кругами, радиус которых мал по сравнению с размерами П. Полагаем, что в окрестностях таких мелких отверстий сосредоточены все резкие "сингулярности" решения:

где и — искомая функция; дП — граница расчетной облаети; д — некоторая известная функция на дП,

* Работа выполнена при частичной финансовой поддержке Российского фонда фундаментальных исследований (грант № 06-01-00421).

© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2008.

-An = 0 В Q u = g на 0Q,

Подобно обычному методу конечных элементов (МКЭ) расчетная область разбита на некоторое число подобластей, называемых суперэлементами. Каждое место сосредоточения особенности (отверстие, неоднородность и т, п.) заключено строго в одном из суперэлементов. Каждая базисная функция МКСЭ Фг(х) единообразно задается в каждом из суперэлементов 0к и является решением задачи Дирихле следующего вида [2-5]:

-ДФг = 0вПк, Фг = ^ на д0к,

где граничные базисные функции <г, заданные на д0к, принимают значения

<г

в узлах суперэлемента Р,-, ] = 1,... ,п, включая возможную границу отверстия в данном суперэлементе, обозначенную через Р0; — символ Кронекера, Узлы Р^ с индексами ] = 1,... ,п расположены только па границе суперэлемента: на его ребрах, в углах. Функции <г(х) удовлетворяют граничному условию пер во го рода на д0 согласно (2):

<г

С узлов па ребра границы суперэлемента функции <г продолжаются некоторым "стандартным" интерполянтом: полиномиальным, кусочно-линейным, сплайн-интерполянтом и т. д.

Заметим, что сингулярности задачи в окрестностях отверстий учтены посредством базисной функции с нулевым индексом Ф0(х) в каждом из суперэлементов. Остальные функции Фг(х), г = 0, при надпчнп отверстия в суперэлементе 0к обращаются в нуль на его границе согласно (3), Если в суперэлементе отверстия нет, то Ф0(х) = 0,

Решение задачи внутри каждого отдельного суперэлемента будем искать при помощи построенного базиса:

п

и

г=0

Таким образом определено приближенное решение МКСЭ и(х) во всей расчетной области 0 = ик0к. Неизвестные значения аг найдены с помощью схемы метода Бубнова— Галеркина при выборе функций Фг(х) в качестве базисных и пробных [3-7],

Для исследования МКСЭ в работах [3, 4] предложен теоретический алгоритм, позволяющий строить его аппроксимации по отношению к широкому кругу задач математической физики. Исследования опираются на общую запись формулы Грина и охватывают слабые решения класса задач, описываемых линейными эллиптическими уравнениями, Установлена связь аппроксимаций МКСЭ с проекционными методами и показано, что для сходимости метода на пространстве слабых решений необходимо и достаточно сходимости аппроксимаций, задаваемых в пространстве их следов. Помимо качественного анализа метода проведено и его детальное численное исследование. Предложены различные варианты МКСЭ, Подтверждена расчетная эффективность для задач, содержащих резкие особенности в расчетной области как в пространственно-двумерном, так и в трехмерном случае [5, 6],

Дальнейший интерес представляет получение априорных оценок погрешностей МКСЭ, Явным образом выделяется вопрос о влиянии изначального выбора варианта аппроксимации на границе разбиения на результирующую точность расчетов метода. Невыясненным является также факт сходимости либо расходимости ошибок производных гладкого решения произвольного порядка. Данная работа посвящена результатам, полученным в ходе дальнейшей разработки и исследования МКСЭ Федоренко,

1. Аппроксимация МКСЭ в пространстве H1^

Качественный анализ различных вариантов МКСЭ осуществляем на примере модельной задачи Дирихле для уравнения Лапласа (1), (2) в пространственно-двумерном случае, Исследование проводим в предположении достаточной гладкости функции граничного условия g и границы области для того, чтобы искомое решение принадлежало пространству Соболева HД(П), R > 1 [8],

Как правило, суперэлемент Qk является многоугольником, и необходимо учитывать тот факт, что принадлежит клaccv C0 гладкости. Мы вправе рассматривать лишь тот случай, когда — многоугольная граница либо граница, состоящая из конечного числа гладких кривых, т. е, = Uilki, l = 1,..., L. Введем обозначения: |11 — длина hi] v — порядок граничного сплайна построенного при разбиении границы S = UkOQk на (N — 1) отрезок; U — приближенное решение; п^u — его сплайн-ннтерпо-лянт [9].

Будем говорить [9], что метод имеет насыщение (насыщаем) на классах HR(Q), R > 1, если 3R0, для которого:

1) lim sup llu — nNu||„= 0 при R < R0;

N — «ей*

2) sup Hu — п^u|H 1

3) при R0 < R существует u G HД(П) и те зависящая от v константа c > 0, такая, что 11 u — п^(u)HH 1(0) > c sup hu — nNuhH 1(0) .

Следующее утверждение базируется на теореме о разрешимости задачи Дирихле в пространстве Соболева H^П) со следами класса H1/2 (S), Оно устанавливает свойство насыщаемости МКСЭ и априорные оценки в H 1(П), Его получение использует теорию интерполяции пространств [10], методы МКЭ [11], результаты насыщаемости при полиномиальной интерполяции в пространстве L2 [12], Исследование МКСЭ в пространстве Соболева H 1(П) сводится к работе со следами решений из H1/2 (S), где мы можем пользоваться и свойствами стандартных оценок при их соответствующем корректном преобразовании,

Утверждение 1. Пусть МКСЭ соответствует интерполяции лагранжевыми сплай-

SS

сеткп |11. Тогда МКСЭ имеет насыщение по гладкости в пространстве H 1(П) на классах Hs(H), s > 1, s G R. Классом насыщения является Hv+3/2(Q), порядком насыщения — O(1/Ns). Справедливы следующие априорные оценки погрешностей метода в пространстве H 1(П).

При v > R — 3/2:

IIм ~~ й11я!(П) - CMR-l/2MR-3/2 (j, 1)Д-1 1/2 1М1ял(П) >

где константа С зависит от параметров исходного оператора и постоянных в неравенствах вложения, а М^-!/2 и М^-3/2 зависят только от К. При V < К — 3/2

\\и — и||Я1(П) < С  111 + НЯ^+3/2(П) ,

где константа С зависит только от V, параметров исходного оператора и постоянных неравенств вложения.

Повышение порядка полиномов на суперэлементных границах повлечет за собой ап-прокеимант более высокой точности МКСЭ в пространстве Н!(0) (для V > К—3/2). Это подтверждено полученными оценками, В случае аппроксимации производных решения при помощи МКСЭ в пространстве Н!(0) такая ситуация изменится, а исследование потребует привлечения дополнительных фактов. Проблема регулярности ограничения приближенного решения МКСЭ на отдельный суперэлемент представляет как отдельный интерес, так и служит для получения оценок погрешностей производных. Далее приведем асимптотику приближенного решения МКСЭ в окрестностях углов суперэлементного разбиения и выясним свойства его гладкости. Результаты, приведенные ниже для угла суперэлемента, можно обобщить на всю расчетную область.

2. Гладкость решения. Аппроксимация производных

Утверждение 2. Пусть Л — один из углов суперэлемента с границей дЛ раствора а, 0 < а < п; Р — его вершина; Ртв — полярная система координат, связанная с Л,

Р

дЛ 1 и дЛ^, представляют собой полиномы порядка не выше V. Тогда приближенное и

и

т9 (ад соъ(дв) + Ъ9 ъ\п(дв)) ,п/а / {Ъ Е < V} ,

9=0

[т9 (а9 соъ(дв) + сдв соъ(дв)) + т91п тс9 ът(дв)] +

9=]'к/а 0<Ч<У

+ Е т9 (ад соя(дв) + Ъя яп(дв)) ]п/а Е {Ъ Е Z,Ъ < V} У]

9=]п/а 0<д<у

1, 2,...,

где ад

А2; ЬдБт(да) = (А* - А2 сов^а)); сд = ^ ~

1

А2

соъ(да) 9

0,

А\, А2 - коэффициенты решения (граничного полинома <) па границах дЛ!, дЛ2 перед т9.

и

в соболевских пространствах. Если п/а Е Z и п/а < V, то решение и нерегулярно по отношению к граничному условию так, что и Е Нп/а(0к) для произвольных а9, Ъ9, вд. В противном случае решение в угле Л обладает бесконечной гладкостью, а приближенное решение в суперэлементе 0к имеет максимальную гладкость по отношению к граничному условию па д0к.

С использованием асимптотики приближенного решения в углах суперэлементов могут быть получены априорные оценки погрешностей производных решения МКСЭ

локально — в данных углах. Для этого использовано разложение произвольной гармонической функции в области суперэлемента П^ со следами клаееа ПНд-1/2(!ы).

I

Разложение может быть найдено из результатов для весовых пространств Соболева "с неоднородной нормой" [13].

Сведем сначала результаты для априорных оценок погрешностей приближения производных решения первого порядка в МКСЭ в норме пространства Н 1(Л) (пли, что одно и то же, оценок самого решения в норме Н2(Л) в угле Л),

Утверждение 3. Пусть вариант МКСЭ соответствует полиномиальной аппроксимации порядка не выше V на границах дЛ в окрестности угла Л суперэлементного разбиения раствора 0 < а < п. Тогда априорные оценки погрешностей МКСЭ в пространстве Н2 (Л) в этом угле имеют следующий вид.

1. Случай п/а / {Ь € Ь < V} . Для многоугольного суперэлемента это эквивалентно выполнению: а/п не принадлежит множеству рациональных приведенных дробей, или же для рационального дробного раствора а/п выполнено V = 0,V = 2 при а < п/2. Тогда при V < Я — 3/2

IIм — м||2,Л < С 111 Ыя>+1

где С зависит только от V, при V > Я — 3/2

R-3/21"R-5/2 UR-5/2

H R-1/2

HR-3/2

где константы Мд_3/2, Мд_5/2 зависят только от Я.

2, Случай п/а € {Ь € Ь < V} . Для многоугольного суперэлемента это эквивалентно выполнению: величина а/п предетавима в виде приведенной рациональной дроби и V > 3, либо V = 2и а = п/2. Тогда при V < Я — 3/2

2,Л

дЛi q=1

где константы зависят только от v, cq соответствуют (6). Погрешность МКСЭ на классе HR(Q) в пространстве H 2(П), т. е. величин a sup ||u — п^ u||20 не обладает сходимо-

«eHR(0) ,

стью к нулю при |11 ^ 0 . Существует некоторая константная погрешность, и справед-

ливо

При v > R — 3/2

II |R-3/2 II" \\2,K — jjR—5/2

Cv,a < sup ||u — п,u12 0 «ей R(0) ,

R—3/2

CR Ыня-1/2^) + CR|du/дп|нл-з/2(оЛ) + CR' X/ |cq1

6Л; q=1

и константы зависят от показателя гладкости R, cq соответствуют (6). Погрешность МКСЭ на классе HД(П) в пространетве H2(П) не обладает сходимостью к нулю при |11 ^ 0. Существует некоторая константная погрешность, и справедливо:

«eHR

Утверждение 4. Для сходимости погрешности метода к нулю в пространстве H2 (Л) при n/a Е {b Е Z,b < V} необходимо и достаточно выполнения условий

AJ cos(qa) — A2 = 0, Уд = jn/a, q = 1,..., min{v, R — 3/2}, j = 1, 2,...

Соответствующие оценки совпадают с п, 1 предыдущего утверждения.

Утверждение 5. Пусть Pxy — система координат такая, что ось Px направлена по одной стороне угла Л и ось Py — по другой его стороне. Для сходимости погрешности метода к нулю в пространстве H2 (Л) при п/a Е {b Е Z,b < v} необходимо и достаточно, чтобы след граничного интерполянта ф на сторонах угла ЗЛ^, i = 1, 2, был следом некоторого гармонического полинома pH(x, y) в координатах Oxy в этом угле.

Можно показать, что при M Е Z, M > 3 норма ошибки ||u — u\\Mл расходится. Метод конечных суперэлементов не аппроксимирует производные решения порядка больше единицы в норме H 1(П),

Заключение

Работа посвящена качественному анализу МКСЭ Федоренко, Получены априорные оценки погрешностей МКСЭ в шкале пространств Соболева на примере задачи Дирихле, Данный вопрос связан с получением аппроксимантов повышенного порядка точности, Определена гладкость приближенного решения. Разрешен вопрос о сходимости производных любого порядка, при наличии такой сходимости результат подтвержден характерными априорными оценками. Метод конечных суперэлементов позволяет разрешать задачи, содержащие мелкие "сингулярности" в расчетной области, обладая при этом погрешностями приближения решения, оцениваемыми только относительно ограничений этого решения на гладких частях суперэлементных границ. Это не требует использования сеток, сгущающихся в окрестностях "сингулярностей", и связано с выбором особых аппроксимирующих пространств. Исследование позволяет установить особенности, свойственные специфике аппроксимаций МКСЭ.

Список литературы

[1] Метод конечных суперэлементов в задачах конвекции-диффузии / В.Т. Жуков, Н.Д. Новикова, Л.Г. Страховская, Р.П. Федоренко и др. \!.. 2001 (Препр. 1111\! им. М.В. Келдыша РАН. № 8).

[2] Федоренко Р.П. Введение в вычислительную физику. М.: МФТИ, 1994.

[3] Galanin \!.. Savenkov е. Fedorenko finite superelement method as special Galerkin approximation // Math. Modelling and Analysis. 2002. Vol. 7, N 1. P. 41-50.

[4] Галанин М.П., Савенков Е.Б. К обоснованию метода конечных суперэлементов Федоренко // Журн. вычисл. математики и мат. физики. 2003. Т. 43, № 5. С. 711-727.

[5] Galanin \!.. Lazareva S., Savenkov Е. Fedorenko finite superelement method and its applications // Comp. Methods in Appl. Math. 2007. Vol. 7, N 1. P. 3-24.

[6] Galanin M.. Lazareva S., Savenkov E. Numerical investigation of the finite superelement method for the 3D elasticity problems // Math. Modelling and Analysis. 2007. Vol. 12, N 1. P. 39-50.

[71 Galanin XL. Savenkov Е., Temis J. Finite Ssperelements Method for elasticity problems // Math. Modelling and Analysis. 2005. Vol. 10, N 3. P. 237-246.

[8] Лионе Ж.-Л., Мл.т/к к и не Э. Неоднородные граничные задачи и их приложения. М.: Мир, 1971.

[9] Локуциевский О.В., Гавриков М.Б. Начала численного анализа. М.: Янус, 1995.

[10] Трибель X. Теория интерполяции, функциональные пространства, дифференциальные операторы. М.: Мир, 1980.

[11] Сьярле Ф. Метод конечных элементов для эллиптических задач. М.: Мир, 1980.

[12] Бавенко К.И. Теоретические основы и конструирование численных алгоритмов задач математической физики. М.: Наука, 1979.

[13] Kozlov V.A., Maz'ya V.G., Rossmann J. Elliptic boundary value problems in domains with point singularities. Providence: AMS, 1997.

Поступим в редакцию 20 февраля 2008 г.