Вариационная постановка

Постановка задачи

МКЭ для двумерной краевой задачи для эллиптического уравнения в полярной вариационная постановка - student2.ru системе координат. Базисные функции билинейные на прямоугольниках. Краевые условия всех типов. Коэффициент диффузии вариационная постановка - student2.ru разложить по биквадратичным базисным функциям. Матрицу СЛАУ генерировать в разреженном строчном формате. Для решения СЛАУ использовать МСГ или ЛОС с неполной факторизацией.

Решаемое уравнение

вариационная постановка - student2.ru -div( вариационная постановка - student2.ru grad вариационная постановка - student2.ru )+ вариационная постановка - student2.ru вариационная постановка - student2.ru

задано в некоторой области вариационная постановка - student2.ru с границей вариационная постановка - student2.ru .

В полярной системе координат:

вариационная постановка - student2.ru .

Краевые условия

вариационная постановка - student2.ru вариационная постановка - student2.ru , вариационная постановка - student2.ru ,

Расчетная область

Расчётная область разбивается на прямоугольники, содержащие по четыре узла.

Локальная нумерация выглядит следующим образом:

Теоретическая часть

Вариационная постановка

Запишем для исходной краевой задачи эквивалентную вариационную постановку в форме уравнения Галёркина, для этого правую и левую часть исходного уравнения домножим на функцию вариационная постановка - student2.ru из пространства пробных функций вариационная постановка - student2.ru и проинтегрируем по вариационная постановка - student2.ru :

вариационная постановка - student2.ru вариационная постановка - student2.ru вариационная постановка - student2.ru (-div( вариационная постановка - student2.ru grad вариационная постановка - student2.ru )+ вариационная постановка - student2.ru вариационная постановка - student2.ru .

Преобразуем полученное уравнение с использованием формулы Грина:

вариационная постановка - student2.ru - вариационная постановка - student2.ru вариационная постановка - student2.ru .

Воспользовавшись краевыми условиями преобразуем интегралы:

вариационная постановка - student2.ru - вариационная постановка - student2.ru вариационная постановка - student2.ru .

В качестве вариационная постановка - student2.ru выберем вариационная постановка - student2.ru пространство пробных функций вариационная постановка - student2.ru , которые на границе S­1 удовлетворяют нулевым первым краевым условиям. При этом будем считать, что вариационная постановка - student2.ru , где вариационная постановка - student2.ru множество функций, имеющих суммируемые с квадратом первые производные и удовлетворяющих только первым краевым условиям на границе S1:

вариационная постановка - student2.ru + вариационная постановка - student2.ru + вариационная постановка - student2.ru = вариационная постановка - student2.ru + вариационная постановка - student2.ru + вариационная постановка - student2.ru . (1.1)

Наши рекомендации