Построение и исследование неконформного метода конечных элементов для решения задачи Стокса с разрывным коэффициентом тема автореферата и диссертации по математике, 01.01.07 ВАК РФ

Рукавишников, Алексей Викторович АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Хабаровск МЕСТО ЗАЩИТЫ
2005 ГОД ЗАЩИТЫ
   
01.01.07 КОД ВАК РФ
Диссертация по математике на тему «Построение и исследование неконформного метода конечных элементов для решения задачи Стокса с разрывным коэффициентом»
 
Автореферат диссертации на тему "Построение и исследование неконформного метода конечных элементов для решения задачи Стокса с разрывным коэффициентом"

На правах рукописи

Рукавишников Алексей Викторович

ПОСТРОЕНИЕ И ИССЛЕДОВАНИЕ НЕКОНФОРМНОГО МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ ДЛЯ РЕШЕНИЯ ЗАДАЧИ СТОКСА С РАЗРЫВНЫМ КОЭФФИЦИЕНТОМ

01.01.07 - вычислительная математика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Хабщюпск 2005

Работа выполнена в Хабаровском отделении Института прикладной математики ДВО РАН

Научный руководитель: доктор физико-математических наук,

профессор Кобельков Георгий Михайлович

Официальные оппоненты: доктор физико-математических наук,

профессор Катрахов Валерий Вячеславович

Защита состоится 8 декабря 2005 г. в 15 часов на заседании диссертационного совета К 212.294.02 в Тихоокеанском государственном университете по адресу:

680035. г. Хабаровск, ул. Тихоокеанская, 136 ауд. 315л.

С диссертацией можно ознакомиться в библиотеке Тихоокеанского государственного университета.

Автореферат разослан 21 октября 2005 г.

Ученый секретарь

кандидат физико-математических наук, Власенко Виктор Дмитриевич

Ведущая организация: Институт математического моделирования

РАН. г. Москва

диссертационного совета

Вихтенко Э.М

2Ш671

I. ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ1

Актуальность темы. Задача Стокса привлекает внимание исследователей в связи с ее многочисленными приложениями, а также тем фактом, что эффективное решение этой задачи открывает путь к решению нелинейных уравнений Навье-Стокса.

Теория применения нестыкующихся сеток для приближенного решения краевых задач в настоящее время является интенсивно развивающимся направлением вычислительной математики. Использование таких сеток неразрывно связано с разбиением, или декомпозицией, расчетной области. В предположении, что область разбита на несколько непересекающихся или пересекающихся подобластей, исходная краевая задача переформулируется на подобластях так, чтобы выполнялись необходимые условия согласования решения на границе между подобластями (интерфейсе). С этой целью, для решения краевых задач эллиптического типа на декомпозиционной области с нестыкующими-ся сетками на интерфейсе, французскими математиками К. Бернарди, И. Мадей и А. Патера в конце 80-х - начале 90-х годов был предложен и изучен новый подход — метод мортарных конечных элементов.

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

Цель работы — построение схемы неконформного метода конечных элементов на декомпозиционной области с использованием мортарных склеек на границе между подобластями для двумерной стационарной задачи Стокса с разрывным (кусочно-постоянным) коэффициентом кинематической вязкости в эллиптической части уравнения; получение оценок скорости сходимости приближенного решения

'Работа выполнена при финансовой поддержке фонда РФФИ (грант 01-01-00375), фонда РФФИ и Администрации Хабаровского края (грант 04-01-97004) и фонда Содействия отечественной науке

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

Общая методика исследований. В диссертации используется понятийный и математический аппарат метода конечных элементов (М К Э), развитый в работах Ф. Брецци, М. Фортэна, Ф. Сьярле, Г. Стрэнга и других авторов. При этом применяются методы и результаты теории дифференциальных уравнений, вычислительной математики и функционального анализа, в частности, теория пространств С.Л. Соболева.

Научная новизна. Новыми в работе являются следующие основные результаты:

1. Для задачи Стокса с разрывным (кусочно-постоянным) коэффициентом кинематической вязкости предложена новая вариационная постановка, учитывающая условия согласования решения на линии его разрыва (интерфейсе): (1) условие слабой непрерывности компонент вектор-функции скорости; (2) равенство потоков скоростей с давлением на функционалах. Исследован вопрос существования и единственности обобщенного решения.

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

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

Теоретическая и практическая ценность. Результаты и методы ис-

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

Апробация работы. Результаты диссертации докладывались на Второй Международной конференции по вычислительной математике (МКВМ-2004; Новосибирск, 2004 г.), на Дальневосточной школе-семинаре по математическому моделированию и численному анализу (РЕ38-МК*ША'01; Находка, 2001 г.), на Конкурсах-конференциях молодых ученых и аспирантов Хабаровского края (Хабаровск, 2004. 2005 г.г), на Дальневосточных математических школах-семинарах имени академика Е. В. Золотова (2002-2004 г.г.), на совместных семинарах лаборатории математического моделирования в физике и технике ВЦ ДВО РАН, кафедр математического анализа Хабаровского государственного педагогического университета и систем автоматизированного проектирования Дальневосточного государственного университета путей сообщения, а также на семинаре Института математического моделирования (Москва).

Публикации. Основные результаты диссертации опубликованы в шести работах, указанных в конце автореферата.

Структура и объем диссертации. Диссертация состоит из введения, трех глав (с делением на пункты и подпункты) и списка использованной литературы — 110 наименований, общий объем работы 146 стр.

II. СОДЕРЖАНИЕ РАБОТЫ

Во введении дается обзор литературы по тематике и излагается содержание работы.

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

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

В п. 1.1 приведена постановка задачи Стокса в прямоугольной области Q с R?, Q = {(я, у) : 0 < х < 2тг. О < у < 7г} с границей Г, П = ft U Г.

Требуется найти вектор-функцию скорости w = (v,v) и скалярную функцию давления р, удовлетворяющие следующей системе уравнений и граничных условий:

-div(ao(x, у) Vw) + Vp = f, divw = 0ei2,

w = 0 на Г, Jpdxdy = 0, С1)

я

где f = (/1, /2) - известная вектор-функция массовых сил; положительная функция ао(х, у) имеет вид

"»frW-U, (х,у)бП2,

где П = ПхиПг и fii П = 0- Разрыв функции ао(я, у) происходит по отрезку Г12 = {(я, у) : х = 7г, 0 < у < 7г}, который будем называть интерфейсом между соседними подобластями Qi и ÍI2, hi = У)"- 0 < х < 7г, 0 < у < ir}, a íí2 =

В п. 1.2 введены пространства C.J1. Соболева на произвольной выпуклой ограниченной области По с Я2, По = По U Го, с достаточно гладкой границей Г0, 70 - часть Го- Для каждой функции г из Я1 (По) определена функция следа г|7о на 70 из Я^2(70); введены также пространства функций Яда2(70) и (#¿o2(7o))'¡ второе из них является дуальным к первому (относительно Ь2). Определены их аналоги для случая вектор-функций. Кроме этого, введены пространства с нормами: U(fi) = {w = (и, и); V е Ь2{П), v € Ь2(П), ик = и\п„ € Н1{Пк), ик = -yin, е Н\ак). w = О на Г, к = 1,2}, 1Н|и,п = (\М\Л + ||w||

V(fi) = {w = (i/, v) e U(íl); [ti]|r„ e о/2(Г12), [«]|r„ 6 #№2)}, l'w||v,n = ПМ1и,я + 11М11н^(Г12))1/2' где [w]|r„ - вектор-функция разности двух следов на Г12, т.е. [w]|r,3 = wllr12n«i ~~ ш21г,>№:

Y(il) = {w = («, v) € V(fi); l [t/]|r„ ■ udy = 0, / jwj|,-,. • >н\у 0.

F 2 I*1J

1 /2

Vf e (Яоо (Г12))'} с нормой пространства V(fi);

ВД = {р е L2(fi);p* =р|Пк 6 ^(пошкп = (11р11о,п, + 11р11одъ)1/2;

X/R(iY) = {p € X(Cl)\ Jpdxdy = 0} с нормой пространства X(Q)\ n

W(fi) = {w e Y(Ji); £ / div w* • 4Mdxdy - О, Щ e Х(П), fc=in*

= 0|nt} с нормой пространства V(Q); Я*(П) = {z e L2(fi); г* = ¿|fit 6 Я'(П*), A = 1, 2, в 6 N},

ll2ll*,fi = (IM|2,n, + IMI2,nJ1/2;

Н'(П) = {w = (и, v);u£ H'{Cl),v 6 Hs{Cl)},

imi:* = ((Ni;. n)2+(Ni;,n)2)1/2;

М(Г12) — пространство таких вектор-функций F, что каждая их компонента принадлежит (Я^2(Г12))', |И1м,гп = sup > где

< Т7,-ц >,,пз= / Т?-~Р<1У при Р € Ь2(Г12) и р е Ьа(Ги). ги

Рассмотрены их свойства и введены вспомогательные обозначения

V' = Ь2(Я), Х' = Ьг{П), М' = Ь2(Г12).

Обобщенное решение системы (1) определим (п 1.3) как решение вариационной задачи: при заданной вектор-функции ( е V' найти А) £ V х X х М, удовлетворяющие системе уравнений

(2)

а(w, <р) + b(tp, р) + d(tp, А) =< f, <р > V<p е V, 6(w, V) = О Щ е X, d(w, р) = О Vp е М и условию

Jpdxdy = 0. (3)

п

В равенствах (2)

2 t it-rt(w. <р) = Z ak(w, p), iu(w, v?) = J ak Vw* • V<pkdr<U) k= 1 fi* 2 и и

Ы^. p) = £ 6jfc(p, p), &*(<£. p) = / -p* • d\vtf*drdy. k= 1 nk

<f.*>= Е4Ы, lk(e)= ffk <pkdxdy, (f* = fk). t-l n*

Л) - Е </*(<л А). А) = | -1)*'1 А • Лу.

л I I ]>

Для согласования решения в (2), (3) на Г12 использованы условия:

1) равенства потоков скоростей с давлением на функционалах

/ + Р'п0 ' <№у = / + Р2пг) -фЛуЧфе Ь2(Г12),

Гц ' Г12

(4)

п* - внешняя нормаль к Г12 — куску границы области к — 1,2;

2) слабой непрерывности вектора скоростей

/ (^ГипП, - ^1г12лп2) • V йу = О УР € м (5)

Г,з

Для того чтобы замкнуть с помощью равенства (4), полученные на подобластях и ГЬ уравнения, в (2) введена вспомогательная вектор-функция А. определенная из соотношения

]Ъ-фйу= ] + Ргпг) • ФЛу УФ € Ь2(Г12).

Г12 Г|2 ^

Следует отметить, что здесь, в отличие от стандартной обобщенной постановки задачи, расширено пространство вектор-функций и <р с Нд(П) до У(П), что позволяет независимо находить решения на подобластях затем склеивать их с помощью условий (4), (5) на Г12.

В п. 1.4 установлено, что при выполнении определенных условий рассматриваемой задачи Стокса существует единственное обобщенное решение. А именно, имеет место

Теорема 1.1. Пусть а(-, •) - непрерывная билинейная форма на V х V, /)('• •) - непрерывная билинейная форма на V х X, а (1(-, ■) - билинейная форма на V х М; 1т В и 1т Б замкнуты в X' и М' соответственно

(В:У->Х':<Вц>,ф >Х'хХ= &(<?, ф)4ч> X,

£> : V —> М' : < Ю<р, /5 >М'*м= ¿(^,/5) 6 е М), т.е. существует постоянная ко > 0 тагая, что

llVllv.il

и существует постоянная к\ > 0 такая, что

¿{ф, р) .. ..

ЭиР тгтп- - А-1 ? М/АегОг,Г„.

»«V Hvllv.il

где ЦрЦм/КегО'-.Гп = II/5 + Д»11м.г„. £>г : М — V' :

р0€Кети' (Iм

< <Р, вт-р ¿(<р, р)У<реУ,Уре М.

Кроме этого, форма а(-, •) коэрцитивна на V/, т.е. существует постоянная ао > 0 такая, что

а(<р, ф) > а0 IMIv.ii Ч? е

Тогда существует (лу, р, X) € V х X х М. удовлетворяющая системе (2) и условию (3) для произвольной Г б V'. Компоненты и р единственны, а А определяется с точностью до элемента КегЮТ. При этом справедливы оценки

1Мкп < ± iifiio.fi, мкп < г (1 + \m\ofl,

«о «о «о

РИм/л^г» < р(1 + Щ(1 + ^)Р11о.П,

М /Со "О

где ||а|| и ||6)| - положительные постоянные в определениях непрерывности форм а(-, ■) и &(•, •) соответственно.

В п. 1.5 построена схема М К Э. С этой целью в пп. 1.5.1 выполняется триангуляция Тд области Г2. Сначала каждую из подобластей Пк вертикальными и горизонтальными линиями разбивается на элементарные квадраты со стороной 2 Ль к = 1,2. Затем каждый из них диагоналями делится на четыре треугольника, которые обозначаются символом е и называются конечными элементами. Далее, вводится Т^ - триангуляция подобласти П^. а через Пи — 11 е, = II е

обозначаются разбиения О, П*- А" = 1.2, соответственно. В качестве узлов аппроксимации (а^'} для каждой компоненты вектор-функции скорости на П*;, выбираются середины сторон конечных элементов с, совокупность которых обозначается через в(кК Подмножество узлов принадлежащих П^иГ^. определяется как

Я'*'. к = 1.2 Выбор узлов аппроксимации для скалярной функции давления не принципиален В случае /г = =4.112 подмножества узлоп аппроксимации 5'"'

и 51"' для компонент на Г|2Г)П| и Г12ПП2 не совпадают, т.е. сетки не стыкуются на интерфейсе Г^.

В пп. 1.5.2 введены конечно-элементные пространства на

1. ^^'(Пал) - подпространство кусочно-линейных непрерывных на е конечно-элементных функций из пространства 1/2(Г2м), в котором г-я базисная функция принимает значение единицы в г-м узле аппроксимации € и значение нуль во всех остальных узлах. При этом базисная функция $ тождественно равна нулю вне конечных элементов, содержащих узел а) .

{к)

2. Хк (Пки) ~ подпространство кусочно-постоянных конечно - элементных функций из пространства /^(Пм), в котором г-я базисная функция (гр,) принимает значение единицы на г-м конечном элементе и нуль на всех остальных элементах.

Конечно-элементное решение системы (2), (3), на подобластях ищем в виде

У) = Е ' У), у) = Е • («. Е у),* =1,2.

{¿.г.етГ}

Определим (пп. 1.5.3) мортарное конечно-элементное пространство на Г12. Отрезки

До = [(тт, 0), 0г, Лж», = ((тг, Лх + 2 0 - 1)/г!)- (», /ц + 2^Лг)],; = 1, 1,

Д* = [(*, /и + 2 (ЛГ - 1) Л,), (*, 7г)], N е N : ЛГ -

N

образующие разбиение Гш = и А», называются мортарными эле-

1=0

ментами. В качестве узлов аппроксимации выбирается подмножество на Г12; обозначим через Л^Гш) подпространство линейных на Д3, з = 1,..., N - 1. и постоянных на До, Дуу непрерывных функций из пространства ^(Гуц,)-

Функции склейки в пространстве Л/ДГ^ь) ищем в виде

а£М = Е А?• = 5? #•

|=г0 «=0

где у), I = 0,..., N — 1 — базисные функции пространства Лл(Г|й). На функции А£ и к = 1,2, наложены условия А/( = А), = -А/„ вк = в^ = -01, которые в свою очередь являются сеточными аналогами компонент А, определяющие конечно-элементную вектор-функцию склейки А„(Х* = (А£,/?£), к = 1,2).

В пп. 1.5.4 сначала введены сеточные пространства на П/, и на Г12Л, затем исследованы их свойства: Ук(Пк) = Ы = (ьи, ьк)\ икк = «Л|П1 €

4 = € Л^)}, 1Ык = (\Ы\1и„ + 11К1111/а,А1г„),/2. где

1К||1А1 = ( Е = (Д Е ММ?.,)1",

НК]||1/аАГ„ = Л-1/211МПо,г12;

х^ь) = {Р= Рл|Пк е Ык = (1ЫЗ,П,+1Ы1^)1/2;

М,(гт) = {17Л; УН е ЛЬ(Г12Л) X Ль(Г12Л)},

1Ы1~1/2АГ„ = № ||Рл||о,Г1з;

Ул(Пй) = € УЛ(ПЛ);/ К1Г„ПЙ1 - ^11Г12ПП2) • = О,

Гц

Уг7л б М/,} с нормой пространства Ул(^л);

\¥Л(ПЛ) = {лу„ € УЛ(ПА); Е Е / • ^ЗД = О, €

Фи — € с нормой пространства V

Доказано (лемма 1.6), что отображение V//, —» является

нормой на Ул, а согласно введенным конечно-элементным пространствам, имеет место соотношения: Ун V, У/, £ У, УУ^ £ Хь С X, М/, С М. Следовательно, рассмотренный метод является неконформным методом конечных элементов.

Для произвольной вектор-функции € У^(Пь). имеет место оценка ||(™Л]||а/2ЛГ1,, < СД^лН^п^ + ^И^) из которой, в частности, вытекает эквивалентность норм || • |{а,ггЛ и II р !|у„ (лемма 1.9).

В пп. 1.5.5 определено приближенное решение задачи (1): требуется найти рд, \к) € V/, х Хи х М&, удовлетворяющее системе уравнений

ан^н, ¥">) + Ьн&ь Рн) + 4(¥>/.- А/,) -< f. V/. >/.• М™^ ^л) = 0, ик) -= О

ДЛЯ любых (yl,, %'h■ I'll) € V/, X Xh х М/, и условию

I phdxdy = 0. (7)

li

В равенствах (6)

2 2 «/»(wh, ¿h)= Т. akh(wh, ¥>h) = Е Е f ajt VwJ ■ Vipfcdxdy, к=i *=i еб7ч*>«

2 2 Ph) = ,E М¥>л> Ph) = E E J" -p£ • div<pfrdxdy,

к=1 Jr=legj-№e

2 2 £ 4(<Pft, ^h) = E 4л(<£Л, \)= £ ; \ •

= lГц

2 2 < f. >h= E /fch(v>h) = E / f* • dx dy, *=1 fc=lil» 2

4(wfcl Рл) = E dkh(v/h, vh) = / (w£|rjjniIl + (-wJ|riannj)) ■ Vhdy.

к— 1 I i2

Равенством rfh(w/,, Х/Л) = 0 определяется конечно-элементный аналог условия слабой непрерывности (5) вектор-функции w/, на Гхг-

Во второй главе установлены оценки (w—w^) в нормах пространств Vи Ъг(Пл) и (р - Ph) в норме пространства Xh(flh)-

В п. 2.1 на основе вспомогательных утверждений получен аналог второй леммы Стрэнга.

Вначале отмечено, что форма dh(-, •) на интерфейсе Г12 удовлетворяет условию Ладыженской - Бабушки - Брецци (утверждение 2.5) и установлен факт существования единственного обобщенного решения задачи Стокса на декомпозиционной области по М К Э (лемма 2.1).

В пп. 2.1.2 сформулирован и доказан аналог второй леммы Стрэнга:

Теорема 2.1. Пусть w up - компоненты решенья задачи (1), удовлетворяющие системе уравнений (2) и условию (3), и w е у(п) п H2(fi). р е Х/ВЩ п Hl{Sl), f 6 ь2(п). Пусть (w/i, Ph- А/,) - приближенное решение, удовлетворяющее системе уравнений (6) и условию (7), wh 6 Yh(&h)-Ph € Xh{Clh), Xh € М/^Гщ,). Тогда имеет место оценка

'lw - wfcl|v„ + \\p~ph\\xh < <c2(.inf ||w-wk||vfc+ inf \\p-qh\\xh +

|n„(w. + b,M, p) + d,M. Xh)- < f. vl > к

. sup-----------¡Г—oil-)'

-IISY, ll^llv,

где С2 - положительная константа.

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

Далее доказываются

Лемма 2.2. Пусть и р - компоненты решения задачи (1), удовлетворяющие системе уравнений (2) и условию (3), V/ € ЛН2(П), р € Х/Я(П)ПЯ'(П). и триангуляция па каждой подобласти к = 1,2, квазиравномерна. Тогда существует положительная постоянная не зависящая от компонент решения V/ и р задачи (1) и к такая, что имеет место оценка

.ы м ||р-?л|к<СзЛ(И|;1П + |ьн:1П).

Лемма 2.4. Пусть XV и р - компоненты решения задачи (1), удовлетворяющие системе уравнений (2) и' условию (3), и XV 6 п Н2(П), р € Х/ЩЛ) п f е Ь2(П). Пусть

Рн, Ад) - приближенное решение, удовлетворяющее системе уравнений (6) и условию (7), £ £ ^(П^). А/, € М^Г^,)

Тогда существует такая положительная постоянная Сц, не зависящая от р, {, пространств У/^Пл), М^Г^ь) и Ь, что имеет место оценка

8ир М*, 'Л) + ЬнЬЪ р) + 4(А- < Л > 1 ^ <с4л(Н|51П + ||р|Ц1П).

При доказательстве леммы 2 4 использована

Лемма 2.3. Пусть А) - решение задачи (1) в вариационной постановке (2), а А/, - компонента приближенного решения в вариационной постановке (6), и при этом е У(П)ПН2(Я) р £ Х/П(П)Г\Н1(0). Тогда существует положительная постоянная Сч, не зависящая оты. р и Ь., такая, что справедлива оценка

||А-Ал||_1/2АГп<С5/г(|Н|^п+||р||:,п).

В пп. 2.2 2 получена оценка ошибок (\v-Wh) в норме пространства Vи {р - Рн) в норме пространства

Теорема 2.2. Пусть w up- компоненты решения задачи (1), удовлетворяющие системе уравнений (2) и условию (3), и w е Y(Q) П H2(ii), р G X/R(fl) П H\Sl), f € L2(il). Пусть (wh- рь■ Ah) - приблиоюенное решение, удовлетворяющее системе уравнений (б)иусловию (7), wh е Yh(0,h),ph G Xh(Qh),Jh G Мл(Гт). Тогда существует такая положительная постоянная Се, не зависящая от w, р, w)n ph, Xh и h, что для проведенной триангуляции Th области Q и разбиения {Д,};^ интерфейса Г12 имеет место оценка скорости сходимости

I|W - Whllv, + IIP - Pfclk <Cth (||w||^ + Ib||i_n).

В п. 2.3 установлена оценка для (w — w/J в норме пространства

L2(Qh):

Теорема 2.3. Пусть w и р - компоненты решения задачи (1), удовлетворяющие системе уравнений (2) и условию (S), и w е У(П) Л H2^), р е Х/ЩП)ПЙ1(П), f G L2(ii). Пусть (wfc. ph, Ah) - приближенное решение, удовлетворяющее системе уравнений (6) и условию (7), wh е Yh{Clh),ph G Xh(ilh),\ е МЛ(ГШ). Тогда существует такая положительная постоянная Ст, не зависящая от w, р, wл, р^, Ад и h, что для проведенной триангуляции Th области Q и разбиения {Д,}^0 интерфейса Ti2 имеет место оценка скорости сходимости

||w-wA||0,nft<C7/i2(||w||*,n + ||p||i>n).

В третьей главе приведено описание численной реализации некон-форного метода конечных элементов для двумерной задачи Стокса в прямоугольнике с кусочно-постоянным коэффициентом в дивергентно

- градиентной части уравнения.

В п. 3.1 дана постановка дифференциальной задачи, а в п. 3.2 рассмотрена ее конечно-элементная аппроксимация.

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

Множество узлов S**', k = 1,2. разбивается на три типа:

- середины вертикальных сторон элементарных квадратов;

- середины горизонтальных сторон элементарных квадратов;

- середины сторон элементарных треугольников, внутренних по отношению к сторонам элементарных квадратов

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

Из полученной системы исключаются узлы (и соответствующие переменные) третьего типа:

(1) для уравнений, имеющих дивергентную структуру, берется линейная комбинация на конечных элементах двух смежных элементарных квадратов (8 треугольников);

(2) остальные уравнения преобразуются следующим образом: из уравнений, соответствующих узлам третьего типа, выражаются переменные того же типа и подставляются в уравнения, соответствующие соседним узлам первого и второго типов.

П. 3.3 посвящен построению и исследованию итерационного процесса для решения преобразованной системы уравнений, матрица которой имеет седловую структуру

' л в' Т ш

ст 0 г

Для нахождения решения системы (8) построен итерационный процесс с переобуславливанием матрицы системы, состоящий из четырех этапов. На первом этапе находится вектор Со, решение системы .4(о = Я по формуле

= % + (9)

На втором этапе - 17. решение системы ¿>77 = С7'(о —

г}п+] = т}п + а13~1(Ст(0 - г- 6г}п). (10)

Здесь 5: 5 = СТА~1В - дополнение по Шуру матрицы системы (8). Далее, на третьем этапе, ищется поправка т к(ос помощью найденного вектора г} как решение системы Лт = Вт}:

тп+1 = тп + а3пА-1(Вг]-Лтп). (11)

На последнем этапе вычисляется вектор

С = (о ~ т-

Здесь пара векторов ff) - решение системы (8); а,1,, а„. nft - параметры процессов (9)-(11), Со"1- ', fn" и т7"- т" ~ значения векторов на (п + 1)-й и п-й итерации (9)-(11) соответственно; ?/". т° - задаваемые начальные приближения процессов (9)-(11). Л и S -персобуславливающие матрицы для Л и <S соответственно.

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

III. ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ

1. Для задачи Стокса с разрывным (кусочно-постоянным) коэффициентом кинематической вязкости предложена новая вариационная постановка, учитывающая условия согласования решения на линии его разрыва (интерфейсе): (1) условие слабой непрерывности компонент вектор-функции скорости; (2) равенство потоков скоростей с давлением на функционалах. Исследован вопрос существования и единственности обобщенного решения.

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

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

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

СПИСОК РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Рукавишников А.В. Численный метод решения задачи Стокса с разрывными коэффициентами в прямоугольнике, часть I/ Алексей Рукавишников// Препринт N 10/ Хабар, отдел. Ин-та прикл. матем. ДВО РАН. - Владивосток, 2002. - 38 с. (2.37 п.л.).

2. Рукавишников А.В. Численный метод решения задачи Стокса с разрывными коэффициентами в прямоугольнике, часть II/ Алексей Рукавишников// Препринт N 04/ Хабар, отдел. Ин-та прикл. матем. ДВО РАН. - Владивосток, 2003. - 20 с. (1.25 п.л.).

3. Рукавишников А.В. Оценка скорости сходимости неконформного метода конечных элементов для задачи Стокса с разрывными коэффициентами/ А.В. Рукавишников, В.А. Рукавишников// Препринт М 71/ Вычисл. центр ДВО РАН. - Хабаровск, 2003. -57 с. (3.5 п.л.).

4. Рукавишников А.В. О дифференциальных свойствах решения задачи Стокса с разрывным коэффициентом/ Алексей Рукавишников// Препринт N 09/ Хабар, отдел. Ин-та прикл. матем. ДВО РАН. - Владивосток, 2004. - 36 с. (2.25 п.л.).

5. Rukavishnikov A.V. Method of numerical analysis for the Stokes problem with discontinuous coefficient/A.V. Rukavishnikov, V.A. Rukavishnikov// Proceedings of the International Conference on Computational Mathematics. ICCM-2004. Part II/ Eds. G.A. Mikhailov, V.P. Il'in, Yu.M. Laevsky. - Novosibirsk, 2004. - P. 907-915 (1 п.л.).

6. Рукавишников А.В. Численный метод решения задачи Стокса с разрывным коэффициентом/Алексей Рукавишников// Вычислительные методы и программирование. - 2005. - Т.б, М1. - С 21-30 (1.12 п.л.).

9

Рукавишников Алексей Викторович

ПОСТРОЕНИЕ И ИССЛЕДОВАНИЕ НЕКОНФОРМНОГО МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ ДЛЯ РЕШЕНИЯ ЗАДАЧИ СТОКСА С РАЗРЫВНЫМ КОЭФФИЦИЕНТОМ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Лицензия ЛР N 040118 от 15.10.96 г. Подписано в печать 14.10.2005 г. Формат 60 х 84/16. Усл. п.л. 1.0. Уч.-изд.л. 0,92. Тираж 100 экз.

Издательство "Дальнаука"ДВО РАН 690041, г. Владивосток, Радио, 7

Отпечатано в Хабаровском отделении ИПМ ДВО РАН 680000, г. Хабаровск, Запари на, 92

»1953g

РНБ Русский фонд

2006-4 21268

 
Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Рукавишников, Алексей Викторович

Введение .".

Глава 1. Построение неконформного метода конечных элементов для двумерной задачи Стокса с разрывным коэффициентом в эллиптическом операторе

1.1. Классическая постановка задачи Стокса.

1.2. Основные обозначения и определения пространств

1.3. Обобщенная постановка задачи Стокса

1.4. Существование и единственность обобщенного решения задачи Стокса

1.5. Схема метода конечных элементов

1.5.1. Триангуляция исходной области.

1.5.2. Определение конечно-элементных пространств на подобластях.

1.5.3. Определение мортарного конечно-элементного пространства на интерфейсе между подобластями

1.5.4. Свойства конечно-элементных пространств

1.5.5. Определение приближенного решения задачи Стокса

Глава 2. Получение оценок скорости сходимости задачи Стокса с разрывным коэффициентом

2.1. Аналог второй леммы Стрэнга

2.1.1. Существование и единственность приближенного решения задачи Стокса

2.1.2. Формулировка и доказательство аналога второй леммы Стрэнга.

2.2. Оценка скорости сходимости (w — w^) в норме Vh{^h) и (p-ph) в норме Xh(Qh).

2.2.1. Вспомогательные утверждения

2.2.2. Получение оценки скорости сходимости

2.3. Оценка нормы (w — wh) в Ь2(П/г).

Глава 3. Численная реализация неконформного метода конечных элементов для задачи Стокса с разрывным коэффициентом.

3.1. Постановка дифференциальной задачи

3.2. Конечно-элементное представление

3.2.1. Функциональный вид.

3.2.2. Алгебраический вид

3.2.3. Конечно-элементная аппроксимация.

3.2.4. Преобразование полученной системы линейных уравнений.

3.3. Построение итерационного процесса.

3.4. Численный эксперимент и анализ результатов . . . 127 Литература.

 
Введение диссертация по математике, на тему "Построение и исследование неконформного метода конечных элементов для решения задачи Стокса с разрывным коэффициентом"

Настоящая диссертация посвящена разработке и исследованию неконформного метода конечных элементов для решения двумерной стационарной задачи Стокса с кусочно-постоянным коэффициентом в эллиптическом операторе с использованием мортарных склеек на линии его разрыва.

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

Первоначально задача Стокса, в ее классической постановке, была решена методами теории потенциалов. Независимо друг от друга Ф. Одквист [91] и Л. Лихтенштейн [86], благодаря построению гидродинамических потенциалов, нашли классическое решение задачи.

Позднее, в работах Д. Лерэ [81]-[83] была дана вариационная формулировка для уравнений Стокса (в общем случае нелинейных) при изучении слабых решений задачи. Существование решения в такой интерпретации вытекает из классической проекционной теоремы (см., например, [28]), при этом само решение называется обобщенным, а вариационная формулировка — обобщенной постановкой задачи Стокса.

В работах Л. Катабриги [55], С. Агмона, А. Даглиса и Л. Нирен-берга [29], В.А. Солонникова [22]-[24], а также в [25], [2], [6], [76], [110], [87] и [7] разработаны различные подходы к исследованию регулярности как классического, так и обобщенного решений стационарной задачи.

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

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

Созданию и исследованию численных методов решения краевых задач посвящено большое количество работ, которые могут быть условно поделены на две части. К первой из них относятся публикации по аппроксимации задач методом конечных разностей, а ко второй — методом конечных элементов (М К Э).

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

Впервые М К Э был предложен Р. Курантом в [61], но это важное исследование тогда осталось незамечанным. Затем в начале 50-х годов прошлого столетия метод независимо был переоткрыт инженерами. Наиболее ранние ссылки, широко встречающиеся в литературе, относятся к работам Д. Аргириса [31], [32], М. Тернера, Р. Клафа, X. Мартина и Л. Топпа [100]. Название метода было предложено Р. Клафом [60].

Позднее, в 60-ые годы математиками, начиная с работ С.Г. Мих-лина [10], [11], была показана значимость анализа метода Галеркина и метода Ритца при построении подхода нахождения приближенного обобщенного решения краевых задач. Хотя математикам в то время не были известны достижения инженеров, изучавшиеся ими приближенные методы, как показывают работы Р. Варги [101], Ж. Сеа [56] в одномерном случае и Д. Биркгофа, М. Шульца и Р. Варги [45] в многомерном случае, напоминали М К Э.

Более подробно теоретические результаты М К Э изложены в монографиях Ж. Деклу [3], Г. Стрэнга и Г. Фикса [26], Ф. Сьярле [27], Э. Митчелла и Р. Уэйта [9], Г.И. Марчука и В.И. Агошкова [8], в работах И. Бабушки и Э. Азиза [36], П. Равьяра [93], М. Зламала [106]. В [4], [20], [21] описаны применения М К Э к решению задач, возникающих в различных областях физики и техники.

По своей природе методы конечных элементов разделяются на два типа. Первый тип — согласованные (конформные) М К Э (см. [3], [27], [9], [8], [36], [93]) в том смысле, что

1) сеточное пространство Sh является подпространством исходного пространства 5;

2) используемые в определении дискретной задачи билинейные и линейные формы тождественно соответствуют формам исходной задачи.

Второй тип М К Э включает в себя методы, нарушающие хотя бы одно из условий согласованности. Из практики вычислений можно выделить три основных способа таких нарушений:

1) Численное интегрирование (численная квадратура, приближенное интегрирование, приближенная квадратура). При этом подходе все еще имеет место вложение Sh С S, но для вычисления коэффициентов результирующей системы линейных алгебраических уравнений используются квадратурные формулы. (Общее введение в предмет численного интегрирования представлено в работах С. Хабера [77], Ф. Дэ-виса и Ф. Рабиновича [64], дальнейшее его изучение проводилось Л. Мэнсфилдом [89], [90], П. Равьяром [93], Ф. Сьярле [58]).

2) Метод криволинейных конечных элементов. При этом подходе сеточное пространство Sh не содержится в исходном пространстве S. Такое нарушение встречается при аппроксимации краевой задачи, поставленной в области с криволинейной границей. Внутри области аппроксимация проводится с помощью прямолинейных конечных элементов (элементов с плоскими гранями), а вблизи границы — с помощью конечных элементов, имеющих по крайней мере одну криволинейную грань. Криволинейные конечные элементы используются для того, чтобы как можно лучше аппроксимировать границу области. Преимущественно применяются изопараметрические конечные элементы, впервые предложенные Д. Аргирисом и И. Фридом в [33], И. Ерга-тудисом, Б. Айроном и О. Зенкевичем в [66] и получившие развитие в работе [59]. Наряду с ними, рассматриваются и другие, отличные от изопараметрических, криволинейные конечные элементы. В связи с этим, прежде всего, необходимо отметить работы М. Зламала [107] -[109].

3) Неконформный метод конечных элементов. При этом подходе также не выполняется вложение Sh С 5, такое нарушение включения, например, для задач второго порядка вызвано тем, что конечно - элементные функции из пространства Sh при переходе через границы конечных элементов терпят разрыв, то есть не являются непрерывными. Метод был разработан и проанализирован Г. Стрэнгом в работах [98], [99].

Что касается задачи Стокса, то при ее конечно-элементной аппроксимации возникает дополнительная трудность — правильный учет условия несжимаемости. Существуют два способа для ее преодоления. Первый подход состоит в том, чтобы использовать такое конечно - элементное пространство, в котором предполагается, что дискретный аналог уравнения несжимаемости выполнен точно. Однако данный процесс, чаще всего, приводит к усложнению элементов. Методы этого подхода были изучены М. Фортэном в [70]. Второй подход, предложенный М. Крузеем и П. Равьяром, основан на аппроксимации условия несжимаемости. В работе [62] этими авторами были построены и изучены как конформные, так и неконформные конечно-элементные пространства. Дальнейшему исследованию методов конечно-элементной аппроксимации задачи Стокса посвящены также публикации [34], [43], [67] - [69], [72], [73], [102]. Одними из первых вопросами получения оценок погрешности аппроксимации, зависящих от гладкости обобщенного решения, занимались Р. Келлог и Д. Осборн [78], Д. Осборн [92] и Р. Темам [28]. Общий обзор по данной тематике изложен в монографии Ф. Брецци и М. Фортэна [53].

Первыми численными методами для решения стационарной задачи Стокса в переменных "скорость-давление"были итерационные методы К. Эрроу, Л. Гурвица и X. Удзавы, предложенные без обоснования на дифференциальном уровне в работе [35]. Несмотря на то, что с момента появления этих подходов прошло уже более сорока лет, они остаются основой для создания новых методов. Вопросы их применения для задач гидродинамики изучали Ж. Сеа, Р. Гловински и Д. Неделек [57], М. Фортэн [70], М. Фортэн, Р. Пейрэ и Р. Темам [71]. Одни из первых экспериментальных исследований проблемы оптимального выбора итерационных параметров р (или риа) были проведены Д. Бежи в [38] и М. Фортэном в [70]. Теоретическое решение этой проблемы для одного частного случая дано М. Крузеем в работе [63]. Отметим ряд публикаций [37], [50], [51], [65], [74], [95], [96], [103], [88] и [79], в которых предложены и изучены модернизированные подходы к построению итерационных методов решения системы линейных алгебраических уравнений, возникающей в результате дискретизации исходной задачи, как с использованием переобусловливающих матриц, так и без них.

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

Тематика применения нестыкующихся сеток для приближенного решения краевых задач привлекает внимание многих исследователей с конца 80-х - начала 90-х годов минувшего столетия. Использование таких сеток неразрывно связано с разбиением, или декомпозицией, области и учетом условий согласования на линиях ее раздела. В связи с этим К. Бернарди, И. Мадей и А. Патерой [44] был предложен новый подход - метод мортарных конечных элементов, который значительно расширяет область допустимых аппроксимаций.

Во-первых, в разных подобластях возможно применение различных сеток, не связанных друг с другом, что позволяет:

- использовать разные формы сеточных ячеек;

- строить "скользящие"сетки, когда аппроксимация задачи с движущимися границами осуществляется за счет сдвига, или скольжения вдоль интерфейса, одной подобласти относительно другой, а не за счет трансформации всей сетки;

- генерировать сетки в подобластях независимо от генерации сеток в соседних подобластях;

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

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

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

В-третьих, в разных подобластях возможно использование различных типов аппроксимаций, например, метода конечных элементов и метода конечных объемов.

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

Метод мортарных элементов, основанный на выполнении условий слабой непрерывности функций на границах между соседними подобластями, оказывается неконформным. Важно отметить, что в первой версии метода, кроме этого, требовалось удовлетворение функций условиям сильной непрерывности в точках на концах интерфейса между подобластями. Во второй же версии мортарного метода, которая уже стала классической, Ф. Бен Велгасемом в [39] было предложено отказаться от этого требования. В работе [40] Ф. Бен Велгасемом и И. Мадей было показано, что использование второй версии метода дает ощутимое преимущество при исследовании трехмерных краевых задач эллиптического типа. Следует подчеркнуть, что на дискретном уровне подобные неконформные методы независимо от них изучали П. Ле Таллек [84], П. Ле Таллек и Т. Соши [85].

Существуют два подхода для построения метода мортарных элементов.

Первый подход (см., например, [48], [39]) основан на использовании: множителей Лагранжа в интегральном тождестве; дополнительного уравнения в обобщенной постановке задачи — условия слабой непрерывности решения на интерфейсе между подобластями.

При построении приближенного решения задачи получается система линейных алгебраических уравнений, матрица которой обладает седловой структурой и имеет аналогичный алгебраический вид, что и матрица системы при дискретизации задачи типа Стокса. В качестве итерационных процедур для решения вышеупомянутой системы уравнений используют неточный алгоритм Удзавы [50], метод сопряженных градиентов в подпространстве, где седловая матрица положительно определена [52], метод минимальных итераций Ланцоша с блочно-диагональным переобусловливателем [80], [1]. Кроме этого, возможно применение многосеточных методов непосредственно к задачам с сед-ловым оператором [105], [48].

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

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

Методика получения оценок скорости сходимости приближенных по методу мортарных конечных элементов решений к точным решениям краевых задач эллиптического типа в нормах различных пространств разработана Д. Браессом в [49], Б. Волмус в [104], Ф. Бен Белгасемом в [39].

Решению задачи Стокса с постоянным коэффициентом кинематической вязкости на декомпозиционной области с нестыкующимися сетками на интерфейсе между подобластями посвящены публикации Ф. Бен Белгасема [41], [42], в которых применен второй подход, в качестве условия согласования на общей границе подобластей использовано условие слабой непрерывности вектора скоростей. Основная направленность первой статьи состоит в том, чтобы установить inf — sup неравенство (глобальную div - стабилизацию), в случае произвольного разбиения исходной области на непересекающиеся подобласти; во второй статье автором на границах между подобластями предложено использовать (для тетраэдральных сеток) метод с добавлением интерфейсных bubble функций [54].

В настоящей диссертации для двумерной стационарной задачи Стокса рассмотрен случай разрывного (кусочно-постоянного) коэффициента кинематической вязкости в эллиптической части уравнения. Исходная область разбивается на подобласти так, чтобы на каждой из них коэффициент вязкости был постоянен. На декомпозиционной области предложена новая вариационная постановка задачи, учитывающая условия согласования решения на интерфейсе между подобластями. Решение поставленной задачи определяется как обобщенное, что позволяет установить его существование и единственность в специальных пространствах. Вследствие специфики рассматриваемой задачи, а именно, разрыва коэффициента кинематической вязкости, неконформный метод конечных элементов характеризуется следующими особенностями: схема М К Э строится независимо на подобластях на основе определения обобщенного решения соответствующей задачи; при этом конечно-элементное пространство для вектор-функции скорости, не являясь подпространством исходного пространства, содержит непрерывные на конечных элементах и кусочно-линейные на подобластях (разрывные при переходе через границу соседних конечных элементов) вектор-функции, а конечно-элементное пространство для скалярной функции давления, являясь подпространством исходного пространства, состоит из непрерывных на конечных элементах и кусочно-постоянных на подобластях функций; сетки на подобластях имеют постоянные, но отличные друг от друга шаги, вследствие чего они не стыкуются на интерфейсе между подобластями; на общей границе двух подобластей введено мортарное конечно-элементное пространство для выполнения сеточных аналогов условий согласования решения.

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

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

 
Список источников диссертации и автореферата по математике, кандидата физико-математических наук, Рукавишников, Алексей Викторович, Хабаровск

1. Василевский Ю.В. Методы решения краевых задач с использованием нестыкующихся сеток// Труды Математического центра им. Н.И. Лобачевского. - Казань: Унипресс, 1999. - С. 94-121.

2. Ворович И.И., Юдович В.И. Стационарные течения вязкой несжимаемой жидкости// Матем. сб. 1961. - Т.53. - С. 393-428.

3. Деклу Ж. Метод конечных элементов. М.: Мир, 1976. - 92 с.

4. Зенкевич О. Метод конечных элементов в технике. М.: Мир, 1975. - 541 с.

5. Ильин В.П. Методы конечных разностей и конечных объемов для эллиптических уравнений. Новосибирск: Изд - во Ин -та математики, 2000. - 345 с.

6. Кошелев А.И. Об ограниченности в Lp производных от решений эллиптических дифференциальных уравнений// Матем. сб. -1956. Т.38. - С. 359-372.

7. Ладыженская О.А. Математические вопросы динамики вязкой несжимаемой жидкости. М.: Наука, 1970. - 288 с.

8. Марчук Г.И., Агошков В.И. Введение в проекционно-сеточные методы. М.: Наука, 1981. - 416 с.

9. Митчелл Э., Уэйт Р. Метод конечных элементов для уравнений с частными производными. М.: Мир, 1981. - 216 с.

10. Михлин С.Г. Вариационные методы в математической физике. -М.: Наука, 1957. 476 с.

11. Михлин С.Г. Численная реализация вариационных методов. М.: Наука, 1966. - 432 с.

12. Рукавишников А.В. Численный метод решения задачи Стокса с разрывным коэффициентом// Вычислительные методы и программирование. 2005. - Т.6, Я 1. - С. 21-30.

13. Рукавишников А.В., Рукавишников В.А. Оценка скорости сходимости неконформного метода конечных элементов для задачи Стокса с разрывными коэффициентами. Хабаровск, 2003. -57 с. (Препринт N 71/ Вычислительный центр ДВО РАН).

14. Рукавишников А.В. Численный метод решения задачи Стокса с разрывными коэффициентами в прямоугольнике. Владивосток: Дальнаука, 2002, часть I, 38 с. (Препринт N 10/ Хабаровское отделение Ин-та прикладной математики ДВО РАН).

15. Рукавишников А.В. Численный метод решения задачи Стокса с разрывными коэффициентами в прямоугольнике Владивосток: Дальнаука, 2003, часть II, 20 с. (Препринт N 04/ Хабаровское отделение Ин-та прикладной математики ДВО РАН).

16. Рукавишников А.В. О дифференциальных свойствах решения задачи Стокса с разрывным коэффициентом. Владивосток: Дальнаука, 2004. - 36 с. (Препринт Я 09/ Хабаровское отделение Ин-та прикладной математики ДВО РАН).

17. Сегерлинд JI. Применение метода конечных элементов. М.: Мир, 1979. - 392 с.

18. Сильвестер П. Феррари Р. Метод конечных элементов для радиоинженеров и инженеров-электриков. М.: Мир, 1986. - 229 с.

19. Солонников В.А. Оценки тензоров Грина для некоторых граничных задач// Докл. АН СССР. 1960. - Т.130. - С. 988-991.

20. Солонников В.А. Априорные оценки для некоторых граничных задач// Докл. АН СССР. 1961. - Т.138. - С. 781-784.

21. Солонников В.А. Об общих краевых задачах для систем эллиптических в смысле Даглиса Ниренберга, ч. 1 // Изв. АН СССР, сер. матем. - 1964. - Т.28. - С. 665-706.

22. Солонников В.А., Скадилов В.Е. О краевой задаче для стационарной системы Навье-Стокса// Тр. Мат. ин та им. В.А. Стек-лова. - 1973. - Т.125. - С. 196-210.

23. Стренг Г., Фикс Г. Теория метода конечных элементов. ~М.:Мир, 1977. 349 с.

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

25. Темам Р. Уравнения Навье-Стокса. М.: Мир, 1981. - 408 с.

26. Agmon S., Douglis A., Nirenberg L. Estimates near the boundary for solutions of elliptic differential equations satisfying general boundary conditions. Part II// Comm. Pure Appl. Math. 1964. - V.17. -P. 35-92.

27. Agouzal A., Thomas J.-P. Une methode d'element finis hybrides en decomposition de domaines// RAIRO Math. Model. Anal. Numer. -1995. V.29. - P. 749-764.

28. Argyris J. Energy theorems and structural analysis. Part I: General theory// Aircraft engineering. 1954. - V.26. - P. 347-356.

29. Argyris J. Energy theorems and structural analysis. Part II: General theory// Aircraft engineering. 1955. - V.27. - P. 42-58.

30. Argyris J., Fried I. The LUMINA element for the matrix displacement method (Lagrangian interpolation)// The Aeron. J. the Royal aeronautical society. 1968. - V.72. - P. 514-517.

31. Arnold D., Brezzi F. Mixed and non-conforming finite element methods: implementation, postprocessing and error estimates// RAIRO Math. Model. Anal. Numer. 1985. - V.19. - P. 7-35.

32. Arrow K., Hurwicz L., Uzawa H. Studies in linear and non-linear programming. Stanford: Stanford Univ. Press, 1958. - 334 p.

33. Bank R., Welfert В., Yserentant H. A class of iterative methods for solving saddle point problems// Numer. Math. 1990. - V.55. -P. 645-666.

34. Begis D. Analyse numerique de l'ecoulement d'un fluide de bingham// These, Universite Paris VI. 1972.

35. Ben Belgacem F. The mortar finite element method with Lagrange multipliers// Numer. Math. 1999. - V.84. - P. 173-197.

36. Ben Belgacem F., Maday Y. The mortar element method for three dimensional finite elements// RAIRO Math. Model. Anal. Numer. -1997. V.31. - P. 289-302.

37. Ben Belgacem F. The mixed mortar finite element method for the incompressible Stokes problem: convergence analysis// SIAM J. Numer. Anal. 2000. - V.37, N 4. - P. 1085-1100.

38. Ben Belgacem F. A stabilized domain decomposition method with nonmatching grids for the Stokes problem in three dimensions// SIAM J. Numer. Anal. 2004. - V.42, N 2. - P. 667-685.

39. Bercovier M., Pironneau O. Error estimates for finite element method of the Stokes problem in the primative variables// Numer. Math. 1977. - V.33. - P. 211-224.

40. Bernardi C., Maday Y., Patera A. A New nonconforming approach to domain decomposition: the mortar element method// In: Nonlinear Partial Differential Equations and their Applications. Paris: Pitman, 1989. - P. 13-51.

41. Birkhoff G., Schultz M.H., Varga R.S. Piecewise Hermite interpolation in one and two variables with applications to partial differential equati- ons// Numer. Math. 1968. - V.ll. - P. 232-256.

42. Bjorstad P.E., Widlund O.B. Iterative methods for the solutions of elliptic problems on regions partitioned in substructured// SI AM J. Numer. Anal. 1986. - V.23, N 6. - P. 1097-1120.

43. Boland J.M., Nicolaides R.A. Stability of finite elements under divergence constraints// SIAM J. Numer. Anal. -1983. V.20, N 4. - P. 722731.

44. Braess D., Dahmen W., Wieners C. A multigrid algorithm for the mortar finite element methods// SIAM J. Numer. Anal. 1999. -V.37, N 1. - P. 48-69.

45. Braess D. Finite elements: Theory, fast solvers, and applications in solid mechanics. Cambridge: Cambridge University Press, 1997. -322 p.

46. Bramble J., Pasciak J., Vassilev A. Analysis of the inexact Uzawa algorithm for saddle point problems// SIAM J. Numer. Anal. 1997. -V.34. - P. 1072-1092.

47. Bramble J., Pasciak J., Vassilev A. Uzawa type algorithms for nonsy-mmetric saddle point problems// Math. Сотр. 2000. - V.69. -P. 667-689.

48. Bramble J., Pasciak J. A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems// Math. Сотр. 1988. - V.50. - P. 1-17.

49. Brezzi F., Fortin M. Mixed and hybrid finite element methods. -N.Y.: Springer Verlag, 1991. - 350 p.

50. Brezzi F., Marini D. Error estimates for the three-field formulation with bubble stabilization// Math. Сотр. 2001. - V.70. - P. 911-934.

51. Cattabriga L. Su un problema al contorno relativo al sistema di equazioni di Stokes// Rend. Mat. Sem. Univ. Padova. -1961. V.31. -P. 308-340.

52. Cea J. Approximation variationnele des problemes aux limites// Ann. Inst. Fourier. 1964. - V.14. - P. 345-444.

53. Cea J., Glowinski R., Nedelec J. Minimization de fonctionelles non differentiables// Lecture Notes in Math. V.228. - Berlin: Springer -Verlag, 1971. - P. 19-38.

54. Ciarlet P., Raviart P. Interpolation theory over curved elements, with applications to finite element methods// Comput. J. Methods Appl. Mech. Engrg. 1972. - V.l. - P. 217-249.

55. Clough R.W. The finite element method in plane stress analysis// Proceedings of the Second ASCE Conference on Electronic Computation. Pittsburg, 1960. - P. 345-378.

56. Courant R. Variational methods for the solution of equilibrium and vibrations// Bull. Amer. Math. Soc. 1943. - V.49. - P. 1-23.

57. Crouzeix M., Raviart P.A. Comforming and non-conforming finite element methods for solving the stationary Stokes equations // RAIRO Anal. Numer. 1973. - V.7. - P. 33-76.

58. Crouzeix M. Sur 1'approximation des equations differentielles operati-onnelles lineaires par des methodes de Runge-Kutta// These, Universi-te Paris VI. 1975.

59. Davis P.J., Rabinowitz P. Methods of numerical integration. N.Y.: Academic Press, 1975. - 459 p.

60. Dyn N., Ferguson W. The numerical solution of equality — constrained quadratic programming problems// Math. Сотр. 1983. -V.41. - P. 165-170.

61. Ergatoudis I., Irons В., Zienkiewicz О. Curved, isoparametric, "quadrilateral "elements for finite element analysis// Internat. J. Solids and Structures. 1968. - V.4. - P. 31-42.

62. Palk R. An analysis of the finite element method using Lagrange multipliers for the stationary Stokes equations// Math. Сотр. 1976. - V.30. - P. 241-249.

63. Falk R. A finite element method for the stationary Stokes equation using trial functions which do not have to satisfy div v = 0// Math. Сотр. 1976. - V.30. - P. 698-702.

64. Falk R., King J. A penalty and extrapolation method for the stationary Stokes equations// SIAM J. Numer. Anal. 1976. - V.13. -P. 814-829.

65. Fortin M. Calcul numerique des ecoulements des fluides de Bingham et des fluides Newtoniens incompressible par la methode des elements finis// These, Universite Paris VI. 1972.

66. Fortin M., Peyret R., Temam R. Resolution numeriques des equations de Navier-Stokes pour un fluide in compressible// J. de Mecani-que. 1971. - V.18, N 3. - P. 357-390.

67. Girault V., Raviart P. Finite element approximation of Navier-Stokes equations// Lecture Notes in Math. V.749. - Berlin: Springer -Verlag, 1971. - 202 p.

68. Girault V., Raviart P. Finite element method for Navier-Stokes equations. Theory and Algorithms. Berlin: Springer - Verlag, 1986. -376 p.

69. Golub G., Wathen A. An iteration for indefinite systems and its application to the Navier-Stokes equations// SIAM J. Sci. Comput. -1998. V.19. - P. 530-539.

70. Gopalakrishnan J., Pasciak J. Multigrid for the mortar finite element method// SIAM J. Numer. Anal. 2000. - V.37, N 3. - P. 1029-1053.

71. Greco D. Nuove formola integrali di maggiorazione per le soluzioni di un'equazione lineare di tipo elliptico ed applicazioni alia teoria del potenziale// Ric. Mat. 1956. - V.5. - P. 126-149.

72. Haber S. Numerical evaluation of multiple integrals// SIAM Rev. -1970. V.12. - P. 481-526.

73. Kellog R., Osborn J. A regularity result for the Stokes problem in a convex polygon// J. Functional Analysis. 1976. - V.21. - R 397-431.

74. Kobelkov G.M. On numerical methods of solving the Navier-Stokes equations in "velocity-pressure"variables// In: Numerical Methods and Applications. Amsterdam: CRC Press Inc., 1994. - P. 81-115.

75. Kuznetsov Yu.A. Efficient iterative solvers for elliptic finite elements problems on nonmatching grids// Russ. J. Numer. Anal. Math. Modelling. 1995. - V.10, N 3. - P. 187-211.

76. Leray J. Etude de diverses equations integrales nonlineaires et de quelques problemes que pose l'hydrodynamique// J. Math. Pures et Appl. 1933. - V.12. - P. 1-82.

77. Leray J. Essai sur les mouvements plans d'un liquide visqueux que limitent des parois// J. Math. Pures et Appl. 1934. - V.13. - P. 331418.

78. Leray J. Essai sur le mouvement d'un liquide visqueux emplissant l'espace// Acta Math. 1934. - V.63. - P. 193-248.

79. Le Tallec P. Domain decomposition methods in computational mechanics// Computational mechanics advances. 1994. - V.l, N 2. -P. 121-220.

80. Le Tallec P., Sassi T. Domain decomposition with nonmatching grids: augmented Lagrangian approach// Math. Сотр. 1995. - V.64. -P. 1367-1396.

81. Lichtenstein L. Uber einige existenzprobleme der hydrodynamik// Math. Zeit. 1928. - V.28. - P. 387-415.

82. Lions J.L., Magenes E. Non-homogeneous boundary value problems and applications. Part I. N.Y.: Springer-Verlag, 1972. - 357 p.

83. Little L., Saad Y. Block LU preconditioners for symmetric and nonsy-mmetric saddle point problems. Minnesota: Minnesota Supercompu-ting Inst., 1999. - 23 p.

84. Mansfield L.E. On the optimal approximation of linear functional in spaces of bivariate functions// SIAM J. Numer. Anal. -1971. V.8. -P. 115-126.

85. Mansfield L.E. On the variational characterization and convergence of bivariate splines// Numer. Math. 1972. - V.20. - P. 99-114.

86. Odqvist F. Uber die randwertaufgaben der hydrodynamik zaher flussi-gkeiten// Math. Zeit. 1930. - V.32. - P. 329-375.

87. Osborn J. Regularity of solution of the Stokes problem in a polygonal domain// In: Numerical Solution of Partial Differential Equations -III (SYNPADE 1975) (Hubbard B. eds.). N.Y.: Academic Press, 1976. - P. 393-411.

88. Raviart P.A. Methode des elements finis// These, Universite Paris VI. 1972.

89. Rusten Т., Winther R. A preconditioned iterative method for saddle point problems// SIAM J. Matrix Anal. Appl. 1992. - V.13. -P. 887-904.

90. Saad Y. Iterative methods for sparse linear systems. New Jersey: PWS, 1996. - 450 p.

91. Shi Z., Xu X. Multigrid for the Wilson mortar element method// Сотр. Methods in Appl. Math. 2001. - V.l, N 1. - P. 99-112.

92. Strang G. Variational crimes in the finite element method// In: The Mathematical Foundations of the Finite Element Method with Applications to Partial Differential Equations (Aziz A. eds.). N.Y.: Academic Press, 1972. - P. 689-710.

93. Strang G. Approximation in the finite element method// Numer. Math. 1972. - V.19. - P. 81-98.

94. Turner M.J., Clough R.W., Martin H.C., Topp L.J. Stiffness and deflection analysis of complex structures// J. Aero Sci. 1956. - V.23. -P. 805-823.

95. Varga R.S. Hermite interpolation type Ritz methods for two-point boundary value problems// In: Numerical Solution of Partial Differential Equations (Bramble J.H. eds.). - N.Y.: Academic Press, 1966. -P. 365-373.

96. Verfurth R. Error estimates for a mixed finite element approximation of the Stokes equations// RAIRO Anal. Numer. 1984. - V.18. -P. 175-182.

97. Verfurth R. A combined conjugate gradient-multigrid algorithm for the numerical solution of the Stokes problem// IMA J. of Num. Anal. 1984. - V.4. - P. 441-455.

98. Wohlmuth B. Hierarchical a posteriori error estimators for mortar finite element methods with Lagrange multipliers// SIAM J. Numer. Anal. 1999. - V.36, N 5. - P. 1636-1658.

99. Wohlmuth B. A multigrid method for saddle point problems arising from mortar finite element discretizations// ETNA. 2000. - V.ll. -P. 43-54.

100. Zlamal M. On the finite element method// Numer. Math. 1968. -V.12. - P. 394-402.

101. Zlamal M. A finite element procedure of the second order of accuracy// Numer. Math. 1970. - V.14. - P. 394-402.

102. Zlamal M. Curved elements in the finite element method. Part I// SIAM J. Numer. Anal. 1973. - V.10. - P. 229-240.

103. Zlamal M. Curved elements in the finite element method. Part II// SIAM J. Numer. Anal. 1974. - V.ll. - P. 347-362.

104. Zygmund A., Calderon A. On the existence of certain singular integrals// Acta Math. 1952. - V.88. - P. 85-139.