chm_3

Классификация методов взвешенных невязок

Классификацию методов взвешенных невязок рассмотрим на примере уравнения Пуассона

(1.0)

с граничными условиями

, (1.1)

(1.2)

где – граница области . Решение задачи ( 1 .0) – ( 1 .2) будем искать в виде конечной суммы

(1.3)

с использованием пробных функций ; – коэффициенты, подлежащие определению. Пусть функции _i(x) выбраны таким образом, что решение, представленное в форме ( 1 .3), точно удовлетворяет граничным условиям ( 1 .1) и ( 1 .2). В этом случае коэффициенты должны быть найдены из условия удовлетворения решения ( 1 .3) исходному дифференциальному уравнению^¹.

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

распределенную в области .

Выберем систему взвешивающих функций , с помощью которой взвесим невязку в области , и потребуем выполнения условия

откуда следует

. (1.4)

Если система взвешивающих функций обладает свойством полноты, то из условия ортогональности невязки _m ко всем функциям _k следует, что . Если, кроме этого, система пробных функций _i(x) обладает свойством замкнутости, то , где u – решение задачи ( 1 .0).

Пусть пробные функции _i(x) выбраны так, что не удовлетворяют точно граничным условиям ( 1 .1) и ( 1 .2) на границах Г_U и Г_Q. Взвесим невязки

на соответствующих границах области,

Теперь, как и в предыдущем случае, можно потребовать выполнения условий ортогональности невязок _U и _Q взвешивающим функциям:

, (1.5)

(1.6)

и решать совместно систему уравнений ( 1 .4) – ( 1 .6) для поиска коэффициентов a_i разложения ( 1 .3) искомого решения. Однако, целесообразно объединить соотношения ( 1 .4) и ( 1 .5), ( 1 .6).

Для дальнейших выкладок воспользуемся теоремой Грина^¹ [9],

. (1.7)

Первоначально положим , то есть ( 1 .7) принимает вид

С учетом этого соотношение ( 1 .4) преобразуется к выражению

Согласно ( 1 .6),

откуда следует

. (1.8)

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

Вновь воспользуемся выражением ( 1 .7), полагая теперь ,

Подставляя полученное соотношение в ( 1 .8), получаем

Поскольку, согласно ( 1 .5),

после преобразований приходим к выражению

. (1.9)

Следует обратить внимание, что в полученном выражении искомая функция вынесена из под дифференциального оператора^². В то же время под знаком оператора Лапласа^³ находятся взвешивающие функции, что ужесточает предъявляемые к ним требования.

Методы (): Методы ():

Рис. 1.1. Классификация методов взвешенных невязок

Если взвешивающие функции подобраны так, так они удовлетворяют уравнению

из выражения ( 1 .9) следует соотношение, связывающие лишь граничные значения искомой функции u_m и ее производной ,

. (1.10)

Полученные соотношения являются основой граничных методов для решения уравнений в частных производных, например, метода граничных элементов, метода граничных интегральных уравнений, метода Трефца и ряда других. На рис. 1 .1 приведена схема классификации методов взвешенных невязок, предложенная в монографии [2].

Содержание