Вычислительные технологии
Том 5, № 1, 2000
МЕТОД НЕОРТОГОНАЛЬНЫХ РЯДОВ ДЛЯ НЕИДЕАЛЬНО СОПРЯЖЕННЫХ
ЭЛЛИПТИЧЕСКИХ ЗАДАЧ С РАЗРЫВНЫМИ КОЭФФИЦИЕНТАМИ
С. В. КостючЕнко Кибернетический центр Томского политехнического университета
Россия e-mail: meko@cc.tpu.edu.ru
A method of nonorthogonal functions is developed to solve a class of nonideally conjugate elliptic problems. The notation of nonideal conjugation conditions is modified. It allowed us to improve the method and develop the universal algorithm using the principle of the calculated domain decomposition into nonoverlapping subdomains. Sewing of solutions in the subdomains is carried out without iterations. Mathematical support is developed to solve problems with concentrated sources. Test problems are solved.
Краевые задачи для эллиптического оператора с разрывными коэффициентами описывают физические поля различной физической природы. К ним относятся поля температур в задачах теплопроводности, напоров и давлений в задачах стационарной фильтрации жидкости в пористой среде и т. д. Известно, что эти и многие другие прикладные задачи являются сопряженными и их решение, особенно для областей сложных конфигураций и различных условий сопряжения, представляет существенные, порой принципиальные, трудности. В постановках таких задач присутствуют не только условия идеального (условия четвертого рода), но и неидеального сопряжения решений вдоль границ разрыва коэффициентов оператора уравнений. О последних идет речь, когда в объектах математического моделирования присутствуют внутренние тонкие слои с аномальными свойствами или когда внутри исследуемых объектов имеются поверхности (границы) — источники или стоки поля.
Такие условия, например, необходимо учитывать при расчетах полей температур в объектах, элементы которых имеют различные теплофизические свойства с неидеальным тепловым контактом — термическим сопротивлением — между ними. Условия неидеального сопряжения необходимы также в случае, когда внутри исследуемого объекта происходят фазовые переходы: кипение, отвердевание, плавление, конденсация и другие процессы, связанные с поглощением или выделением тепла. Аналогичные проблемы имеют место при решении задач прикладной теории фильтрации жидкости в пористых средах. Так, неидеальное сопряжение имеет место в случае природных (например, тектонических) или
© С. В. Костюченко, 2000.
техногенных нарушений строения или фильтрационных свойств водо- и нефтеносных пластов, при наличии межпластовых перетоков по трещинам или "линзам" и т. д. "Неидеально" сопрягать решения краевых задач в подобластях необходимо в случаях формирования геологических карт при имеющихся нарушениях в строении недр. В этом случае строение и изменчивость параметров, характеризующих земную толщу, интерпретируются как некоторые потенциальные поля и описываются краевыми задачами для эллиптического уравнения.
Ниже будет показано, что и при решении "классических", идеально сопряженных, краевых задач методом разделения по подобластям может возникнуть необходимость перехода к неидеально сопряженным задачам.
Проблемам развития методов, создания математического и программного обеспечения для решения неидеально сопряженных эллиптических задач, на наш взгляд, уделено недостаточное внимание [3, 6, 8, 9]. Поэтому разработка и развитие численных и численно-аналитических методов для решения таких задач остаются актуальными.
Одним из эффективных методов решения краевых задач показал себя метод разложения по неортогональным функциям [1, 7] — неортогональных рядов [4]. Метод имеет теоретическое обоснование, но для решения задач в кусочно-однородных областях даже с идеальным сопряжением практически не применялся.
В данной статье этот метод развит для решения класса неидеально сопряженных эллиптических задач. Применение модифицированной записи условий неидеального сопряжения позволило усовершенствовать метод и разработать достаточно универсальный алгоритм, использующий принцип разделения расчетной области на совокупность непересекающихся подобластей со сшивкой решений без каких-либо итераций.
1. Математическая постановка задачи
Пусть в пространстве Rn, n = 2, 3 некоторая многосвязная расчетная область П с границей дП определена совокупностью непересекающихся подобластей П1,, таких, что П = иП с границами дПи, v = 1,... ,n. И пусть в каждой из них поставлена краевая задача с уравнением
Lv [Uv(X)] = -Gvdiv[gradUv(X)] = Fv(x), X g П. (1)
Искомое решение Uv (X) в каждой из подобластей и во всей расчетной области должно находиться с учетом условий, заданных как на внешней границе дnv0 расчетной области П, так и на внутренних границах дП^, разделяющих смежные подобласти nv и Пм с номерами v, ^ соответственно. Коэффициенты Gv принимают в подобластях Пv постоянные значения и претерпевают разрывы на границах Пvß раздела подобластей. Здесь дПм = дП0 u дП1, — покрытие границы П, Fv (X) — заданная функция.
На искомую функцию Uv(X) на внешней границе дП^, в общем случае, накладываются условия I - III родов и смешанные. Они могут быть записаны единой формулой:
ди
lvUv Wo = \ Av(Y) ^ + Bv(Y)Uv(X)
(Y), Y g д^0. (2)
Здесь Ли(V), Би(V) и Фи(V) — известные кусочно-гладкие функции, ии — вектор внешней нормали к границе , 1и — дифференциальный оператор краевых условий.
На границах дП^ = ди дQv смежных подобластей О и (V, ^ = 1,т; V = условия сопряжения могут быть записаны в общем виде:
п
v0
иV (X )|
о,
д!Л
д—
дП
дQv ,
и,(х) + я(у )о, дц-
дП
о,
ди
ди,
+ д(У), У е дП
V,-
(3)
дП
Условия (3) определяют условия неидеального сопряжения решений иу и и, вдоль границ дПу, подобластей Пу и П,, являются обобщением известных условий четвертого рода, весьма распространенных в постановках задач для эллиптических уравнений, и задаются функциями Я(У) и Q(У). В частности:
— при Я = 0, Q = 0 это условия "идеального" сопряжения, они обеспечивают нераз-
ди
рывность решения и и его функции О--— вдоль границы дПу,;
— при Я = 0, Q = 0 для решения и условия (3) определяют разрыв первого рода вдоль границы дПу, и имеют физический смысл тонкого слоя — сопротивления;
о т т
— при Я = 0, Q = 0 для производной д— условия (3) определяют разрыв второго рода вдоль границы Пу, и имеют смысл тонкого слоя — источника.
Известно, что решение и, удовлетворяющее уравнению (1) и краевым условиям (2)-(3), существует и оно единственно.
2. Модифицированные условия неидеального
сопряжения и краевые условия для подобластей
Условия сопряжения в форме (3) не всегда удобны для построения универсальных алгоритмов. Поэтому запишем их в модифицированном "симметричном" виде, определяемом следующей леммой [5].
Лемма. Условия (3) и условия
иу + Оу
иу + Оу
Я 2~
Я 2~
+ С
- С
диу
ди V
диу
д—у
+
Я
2
0_
2
и, + О,
и, + о.
Я 2~
Я 2~
-С
+ С
ди,
ди
, + ^ + 2
ди,
0.
2
(4)
для У е дПу, эквивалентны, если Я(У) = 0 и Q(У) = 0 одновременно, а — константа, Я
= 0, С = ±т
Доказательство этого утверждения производится почленным сложением, вычитанием и последующим упрощением получаемых выражений.
Проведем необходимые для дальнейшего изложения преобразования. Известно, что решение иу задачи (1)-(3), линейной относительно входящих в нее переменных, может быть представлено в виде суперпозиции иу = Уу + Шу, где Уу (X), Шу (X) — непрерывные дважды дифференцируемые в Пу функции. Пусть Шу (X) — одно из известных частных решений уравнения задачи, для которого по определению Ьу(Шу) = Ьу. Отметим, что для рассматриваемого оператора Ьу и правых частей Ьу некоторых специальных видов частные решения Шу могут быть записаны явно [1, ?].
В результате очевидных преобразований переходим к постановке задачи для новой искомой функции Уу (X):
Ьу Уу (X) = 0,
¡V Уи (X )|до = ^ (V),
(5)
где
Фи(V) = ф,(V) - ¡V—(X)|
до,,,
(6)
Для задачи (5) получим такой дифференциальный оператор 1и и запишем такую функции ^(V), которые, задавая краевые условия для решений Уи(X) задач в подобластях по отдельности, априори обеспечивали бы выполнение условий (2) - (3) для решения и(X) и решений и(X) в составной расчетной области 0. Предположим, что некоторая подобласть О такова, что ее граница
дПv = дПv0l и дПv02 и дПvoз и дП^-
На участках дПV0г, ъ = 1, 2, 3 границы дПV в результате преобразования (6) запишем условия первого, второго и третьего родов соответственно:
У' (x )1ап„о1 дК (X)
д^ + К V
Фv0l(V) - — |до^01 = 01 (V), д—
Фv02(V) -
„ + hV у V
ди
дО
дП^оз
дп
= <^02 (V),
дП„
Фv0з(V) -{ дд-и^ + ^-V
^03^ )
(7)
дП^оз
Здесь Ф^^) — заданные функции, ъ =1, 2, 3.
На участке дП^,, границы дПV аналогично запишем условия сопряжения в форме, определяемой Леммой (4), что обеспечит выполнение условий сопряжения (3):
ди
дп„
+
Uv
Gv[ § + С
и, + Сц 2
дП„
Я - С^ + сд
ди,
Я + С
= )
дП
или
где
дК + К
дп
Р
у, + сп ^ - С1 ^ + сд
ди,
Я
'И 2
дП„
Р„
+ Я
, и, + Ч | - - и.
¿V = Р
Р V
дО
^ р = + с
(8)
Очевидно, для задач в подобластях возможен произвольный набор из названных краевых условий.
Таким образом, из (7) - (8) следует, что дифференциальный оператор ^ задачи (5) может быть записан в общем виде в форме условий (2), где кусочно-гладкие функции AV (V) и BV (V) принимают значения, соответствующие заданному типу условий (табл. 1).
Таблица 1
Компоненты краевых условий для задач в подобластях
Тип условий А (Г) ^ + в (у ^ (у)
А» (У) В (У)
I 0 1 <Ри0\(У)
II 1 0 ^02(У )
III 1 ¡V ^оз(У)
IV 1 1 Ч!+с) )
Отметим, что переход от Фи(V) к (V) по формуле (6) приводит, в общем случае, к ду д V
К(V) = У'(У) и = — 00'даже при К = 0 и ^ = 0, т.е. в результате перехода к
однородным уравнениям задач (5) в подобластях получаем также некоторые неидеальные условия сопряжения.
3. Метод неортогональных рядов для задач с разрывными коэффициентами
Эффективность метода решения задачи (5) - (8) в кусочно-однородных расчетных областях сложной конфигурации определяется эффективностью метода решения задач в каждой из подобластей. Исходя из условий (8), метод решения, прежде всего, должен обеспечи-
дУ
вать получение искомых функций У'(V) и в аналитическом виде для точек V € д
а также нахождения функции Уи (X), X € дО,и и с условиями (7) - (8). Переход от Уи (X) к искомому решению (X), X € ди будем считать тривиальным. Поэтому целесообразно, чтобы метод позволял понижать размерность решаемой задачи при переходе от подобластей к их границам, при этом упрощая и обработку геометрической компоненты постановки краевой задачи. К таким методам относятся методы граничных интегральных уравнений, граничных элементов, метод Трефтца [10] и их модификации, одним из которых является метод неортогональных рядов [1].
3.1. Метод неортогональных рядов для задач в подобластях
Основная идея метода неортогональных рядов состоит в построении некоторой вспомогательной поверхности дО,и*, описывающей ди и не имеющей с ней общих точек, введении системы функций Фу = Фу (Zj ^) — фундаментальных решений оператора Ьи и поиске решения Уи (X) задачи (5) в виде разложения:
ри
у^ (X ) = £; (ZjV X). j=l
(9)
Здесь Фу(Zj^) — фундаментальные решения оператора Ьи, Zk € д*, {Zk} € ^ и д, к = 1, то — всюду плотное на д* множество точек, CjV — неизвестные коэффициенты, ри — число аппроксимирующих функций. Для многосвязных областей необходимо построение нескольких поверхностей д*.
Известно [1, 7], что система функций Ф^ полна в классе функций 2(Пг,) и линейно независима, т. е. при достаточно большом разложение (9) может со сколь угодно высокой точностью аппроксимировать решение задачи в подобласти П.
По определению, каждая функция Ф1У удовлетворяет однородному уравнению задачи (5). Разложение (9) также удовлетворяют уравнению ¿(К,) = 0 при любом выборе постоянных С и. Их поиск производится исходя из условия наилучшего удовлетворения граничным условиям и минимизации функционала (10) для системы из т подобластей:
2
I
р V
£ Р V ) [С V1 V (Ф^- V ,У ))]
V=1
, (10)
^=1 > до
Метод неортогональных рядов, использующий системы базисных функций Ф1У — фундаментальных решений, которые по определению удовлетворяют уравнению задачи, можно рассматривать как одну из модификаций метода Трефтца.
Переход от функционала (10) к соответствующей системе линейных алгебраических уравнений [А^] • {С,-} = {В^} может быть произведен, например, по методу наименьших квадратов или методу коллокаций. В результате решения системы уравнений находится вектор коэффициентов {С}. После их подстановки в (9) и перехода к искомой функции и задача считается решенной.
3.2. Апостериорные оценки погрешностей численных решений
Возможность получения достоверных и объективных оценок точности решения краевых задач важна при использовании любого численного или численно-аналитического метода. В условиях, когда получение реальной априорной оценки погрешности численного решения получить чрезвычайно трудно, не меньшее значение приобретают апостериорные оценки, если, конечно, их получение возможно.
В методе неортогональных рядов основным источником погрешности решения задачи является погрешность аппроксимации краевых условий, что позволяет получать эффективную апостериорную оценку погрешности полученного численного решения задачи. В самом деле, для задачи в подобласти П и краевых условий на внешней границе расчетной области дП = идП^ такую оценку дает невязка Е, вычисляемая для узлов У' £ дП, не совпадающих, в общем случае, с узлами У:
2
Е =
pv
р V (у ') -£ (С V Фv(z V ,У'))
¿=1
дПь
В соответствии с известным принципом максимума для гармонических функций [10] погрешность решения задачи в точках X £ П не может быть больше Е : < Е.
Для условий сопряжения (8), а следовательно, и условий (3), аналогичную оценку дает невязка XV,, вычисляемая для узлов У" £ дП также не совпадающих с узлами У:
V,
р V
р V, (у ") - е ¿=1
С
¿V
д Фv(Zj V ,У") + Фv(Zj V ,У ")
дп
дП„
Все используемые здесь обозначения описаны выше.
Таким образом, невязки Е и XV,, вычисленные для всех подобластей с номерами, V, ^ = 1, т; V = ^ в совокупности дают эффективную апостериорную оценку погрешности решения задачи (1)-(3).
2
3.3. Обработка геометрической информации и аналитические компоненты метода
Основными аналитическими компонентами, применяемыми в данной реализации метода неортогональных функций, являются частные решения (X), фундаментальные реше-
.т, (7 ^ яп Жу(X) ф„(£.,• V,х) ния ,Х) и производные от них по нормали к границе дПу:—^^— ,-- в точ-
ках У Е дПу, а также кусочно-гладкие функции А(У), В (У), Д(У), ф(У), к у (У), ^ у0г(У), где г = 1, 2, 3.
Формирование и использование этих компонент при численной реализации метода связано с общей проблемой обработки геометрической информации о форме расчетных областей и их границ при решении краевых задач математической физики. Вероятно, с наиболее общих позиций эту проблему удается решить при применении конструктивных средств теории Д-функций [11].
Эти средства позволяют построить такие аналитические описания (X) подобластей Пу, что Пу = (X) > 0], Е Ст(П) и (X) = 0 — уравнение дПу. Отдельные фрагменты границ дП, имеющие аналитическое описание, используются здесь в качестве элементарных функций. Из этих фрагментов нетрудно сконструировать также и такие аналитические функции — описания элементов границ д П вдоль которых функции А(У), В (У), Д(У), ^(У), ку(У), ^у0г(У) сохраняют значения в классе 1(Пу). Тогда "склейка" значений этих функций, заданных на участках дПграницы дПу = идП, г = 1, / , и продолжение их внутрь расчетной области П могут быть произведены, например, по формуле
I
Ф = ^-, (11)
Е-
Тг
¿=1
где тг = ¿, ^^^) = 0 — уравнения элементов границ, на которых заданы значения функции /¿.
Имея уравнение подобласти Пу, нетрудно вычислить производную от некоторой аналитической функции по внешней нормали к границе дПу:
д/
— = -ЕГвЛ(/)grad(wv).
В качестве функции / здесь могут использоваться Ф^, У.) — фундаментальное решение или их суперпозиция и (X, У) — одно из частных решений оператора Ьу.
Для двухмерной расчетной области и Ь у — оператора Лапласа известное фундаментальное решение Ф^, У.) = Ьп(Д)/(2 * п), где Д = ^((х — ¿.¿) • (х — ¿.¿)), суммирование по г = 1, 2, а
д Ф д Ф
gradФ((X,Уj)) = — г + —= — [(Х1 — У. )/Дг + (х — )/Д?].
Для решения класса прикладных задач нам потребуется частное решение (X) для
— суперпозиции ¿-функций — точечных источников поля и (X), каждый из которых имеет удельную мощность Qi,г = 1,¿:
ху = — £ — X¿).
По определению, фундаментальное решение ^^) дифференциального оператора Ьи — это функция, для которой ЬиФV(X, V) = — V). Поэтому, если правая часть уравнения (1) представлена дельта-функцией, т. е. является точечным источником поля и, то фундаментальное решение само является одним из частных решений оператора Ьи:
№ ^^) = ФУ (X — V).
Очевидно, суперпозиция фундаментальных решений также является частным решением для правой части, представимой в виде суперпозиции дельта-функций:
г г
IV[№(X, V)] = ^ Qj5^ — Vj), где №(X, V) = ^ QjФ(X — V). j=l j=l
Эти частные решения будут использованы при решении сингулярно возмущенной краевой задачи с сосредоточенными источниками.
4. Результаты вычислительных экспериментов
4.1. Результаты решения тестовой задачи
В качестве тестовой рассматривалась задача для уравнения Лапласа в сечении бесконечной двухслойной трубы с заданными условиями Дирихле (первого рода) на внутренней и внешней поверхностях. Радиусы внутренней трубы: Т\ = 1.0 и г2 = 1.5; внешней — г2 = 1.5 и г3 = 2.0. Краевые условия: и(п) = — 1п2.0/2 = —0.3466, и(г3) = 0.0. Параметры 01, 02, К, Q — см. (1) - (3) — считались известными. Для случая К = 0, Q = 0 известно аналитическое решение задачи. Решения задачи для различных значений параметров 01, 02, К, Q приведены в табл. 2.
Результаты получены при следующих параметрах: = 1 , число аппроксимирующих функций для каждой из подобластей щ = 20, V =1, 2, число узлов коллокации на каждой из границ 10. Вспомогательные границы дП^, * — для каждой из подобластей по две окружности радиусов т\* = 0, 3 и г2* = 10, 0. Узлы коллокации и узлы на вспомогательных границах д^* размещались равномерно. Поиск коэффициентов CjV разложения (9) при минимизации функционала (10) производился методом наименьших квадратов. Результаты вычислительных экспериментов подтвердили корректность разработанного математического обеспечения и работоспособность алгоритмов. В частности, относительная ошибка решения тестовой задачи в сравнении с аналитическим решением (см. табл. 2), N = 1 — 8 — менее 0, 0001 %, что подтверждает хорошие аппроксимационные свойства систем фундаментальных решений и эффективность получаемых апостериорных оценок. Эти оценки использованы при анализе погрешностей решения задачи при неидеальном сопряжении решений в подобластях N = 9 — 16.
4.2. Расчет поля давлений в нефтяном пласте с системой скважин
Разработанное и описанное выше математическое и программное обеспечение было апробировано при решении задач расчета поля пластовых давлений в кусочно-однородном нефтяном пласте с системой нагнетательных и нефтяных скважин. Поле пластовых давлений описывается уравнением Пуассона с правой частью в виде суперпозиции дельта-функций — точечных источников процесса фильтрации. Расчетная область составлена из двух подобластей ^ и П2 с гидропроводностями 01 и 02 соответственно. Подобласти
Таблица 2
Результаты решения тестовой задачи при различных параметрах
Численное решение Точное решение
N К я С1 С2 г Ev | т^/и | дUv/д«V и дUv/д«^
1 0.0 0.0 1.0 1.0 1 1.0 10-6 — —0.3466 —0.500 —0.3466 —0.500
2 1.5 — 10-9 —0.1438 0.333 —0.1438 0.333
3 2 1.5 — 10-9 —0.1438 —0,333 —0,1438 —0,333
4 2.0 10-4 — 0.0000 0.250 0.0000 0.250
5 0.0 0.0 1.0 10.0 1 1.0 10-4 — —0.3466 —0.798 —0.3466 —0.798
6 1.5 — 10-5 —0.0230 0.532 —0.0230 0.532
7 2 1.5 — 10-6 —0.0230 —0.0532 —0.0230 —0.0532
8 2.0 10-7 — 0.0000 0.0399 0.0000 0.0399
9 1.0 0.0 1.0 10.0 1 1.0 10-5 — —0.3466 —0.315 — —
10 1.5 — 10-5 —0.2189 0.210 — —
11 2 1.5 — 10-6 —0.0091 —0.0211 — —
12 2.0 10-8 — 0.0000 0.0157 — —
13 0.0 1.0 1.0 10.0 1 1.0 10-4 — —0.3466 —0.898 — —
14 1.5 — 10-4 0.0173 0.598 — —
15 2 1.5 — 10-5 0.0173 0.0402 — —
16 2.0 10-7 — 0.0000 —0.0301 — —
имеют общую границу дП12. Конфигурация расчетной области и схема расположения 17 нагнетательных и 18 нефтяных скважин приведены на рис. 1. Все расстояния и размеры указаны в километрах. На внешних границах дП01 и дП02 заданы условия первого рода и =15 МПа. Среднесуточные дебиты нагнетательных скважин +0.008 м3/с, нефтяных —0.003 м3/с. Условия сопряжения решений и1 и и2 в подобластях и П2 вдоль границы дП12 имеют вид (3). Параметры К и Q определяют "идеальность" и "неидеальность" сопряжения решений.
Вариантные расчеты проводились при различных значениях параметров 01, 02, К и Q. Результаты некоторых вычислительных экспериментов приведены на рис. 2-4. На рис. 2 показано расчетное поле давлений в однородном нефтяном пласте с идеальными условиями сопряжения вдоль границы дП12(К = 0, Q = 0). Здесь гидропроводности подобластей 01 = 02 = 0.001 км3/МПа/с. Такие значения гидропроводностей соответствуют проницаемости пласта 0.5 мкм2, вязкости нефти 5 МПа-с, толщине пласта 10 м. На рис. 3 показано расчетное поле для вариантов К = 10, Q = 0 и К = 0^ = 10, на рис. 4 — расчетное поле для составных условий сшивки на границе дП12 и кусочно-однородной расчетной области с соотношением 01/02 = 2. На всех рисунках идеальные и неидеальные условия сопряжения прослеживаются по характеру прохождения и разрывам линий равных значений пластовых давлений при согласовании с физическим смыслом этих условий.
2.00 -1.50 -1.00 -0.50 0.00 0.50 1.00 1.50 2.00
Рис. 1. Конфигурация расчетной области для задачи расчета поля пластовых давлений в кусочно-однородном нефтяном пласте с системой нефтяных (о) и нагнетательных ( ) скважин.
2.00 -1.50 -1.00 -0.50 0.00 0.50 1.00 1.50 2.00
Рис. 2. Поле пластовых давлений в однородном нефтяном пласте с системой скважин. На внутренней границе заданы условия идеального сопряжения: К = 0, Q = 0.
i i i i i i i i i
2.00 -1.50 -1.00 - 0.50 0.00 0.50 1.00 1.50 2.00
2.00 -1.50 -1.00 - 0.50 0.00 0.50 1.00 1.50 2.00
Рис. 3. Поле пластовых давлений в однородном нефтяном пласте с системой скважин. На внутренней границе заданы условия неидеального сопряжения: К = 10.0, Я = 0 (вверху); Я = 10.0, К = 0 (внизу).
Рис. 4. Поле пластовых давлений в кусочно-однородном нефтяном пласте с 17 нагнетательными
и 18 нефтяными скважинами. Соотношение гидропроводностей подобластей = 2. На внутренней границе ЛБСР заданы сме-
шанные условия сопряжения: на дуге ЛС — идеальное сопряжение подобластей: К = 0, Q = 0; на АВ — неидеальное сопряжение: К = 0, Q = — 20сов(/) вт(/), 270° < / < 360°; на ВС — неидеальное сопряжение:
Q = 0, К = 20соэ(/) яш(/), 0° < / < 90°.
Все результаты получены при разложении решений в подобластях по 60 базисным функциям для каждой из подобластей. В качестве базисных функций использовались фундаментальные решения оператора Лапласа — 1п 1 / Д. Полюса фундаментальных решений размещались на окружностях: для П1 — Д = 0.45 и Д2 = 3.5 с центрами в точках (—0.6, 0.0) и (0.0, 0.0); для П2 — Д = 0.25 и Д2 = 1.3 с центрами (—0.8, 0.0) и (—0.6, 0.0) соответственно. Такой выбор базисных функций обеспечил приемлемые значения невязок Е Т
5. Заключение
В данной работе метод неортогональных функций развит для решения класса неидеально сопряженных эллиптических задач в кусочно-однородных расчетных областях. Существенное развитие метода удалось произвести при использовании модифицированной записи условий неидеального сопряжения решений в подобластях. Разработанный универсальный алгоритм использует принцип декомпозиции расчетной области на непересекающиеся подобласти. Сшивка решений в подобластях производится без каких-либо итераций. Разработано также математическое обеспечение для решения задач с сосредоточенными
источниками в правой части уравнений. Решены тестовые задачи.
За рамками статьи остались некоторые проблемы, свойственные методу неортогональных рядов: трудности обработки плотно заполненных плохо обусловленных матриц, неоднозначность размещения полюсов фундаментальных решений — базисных функций и другие. Однако, как показала практика, многие из этих проблем успешно преодолеваются в результате проведения эффективных апостериорных оценок получаемых численных решений и минимизации соотвествующих невязок.
Список литературы
[1] АлЕксидзЕ М. А. Решение граничных задач методом разложения по неортогональным функциям. Наука, М., 1978.
[2] Бутковский А. Г. Характеристики систем с распределенными параметрами: Справ. пособие. Наука, М., 1979.
[3] Карташов Э. М. Аналитические методы в теории теплопроводности твердых тел. Высш. шк., М., 1985.
[4] Ковылянский Ю.В., ПопковА Т. Л., Слепян Г. Я. О методе неортогональных рядов для численного решения задач математической физики. ЖВМиМФ, 28, №2, 1988, 237-246.
[5] КостюЧЕНко С. В. Математическое и программное обеспечение для решения задач пространственно-распределенной стационарной диффузии и тепломассообмена методом композиции. Дис. ... канд. техн. наук. Томск, 1987.
[6] Кузнецов С. Б. Алгоритм решения квазилинейных эллиптических уравнений в областях с источниками на внутренних границах. Новосибирск, Изд. ВЦ СО АН СССР, Препринт №476, 1983, 21 с.
[7] КУПРАДЗЕ В. Д. О приближенном решении задач математической физики. Усп. ма-тем. наук, 1967, 22, №2 (134), 59-107.
[8] ЛАРДЖЕР Т. И. Приближенные решения задач переноса тепла при наличии фазовых превращений. Ракетная техника и космонавтика, 1967, №11.
[9] Макаров А. М., Романовский В. Р. О неидеально сопряженной задаче теплопроводности. Изв. АН СССР, Энергетика и транспорт, 1981, №5, 130-137.
[10] Михлин С. Г. Вариационные методы в математической физике. Наука, М., 1970.
[11] РвАЧЕв В. Л. Теория Я-функций и некоторые ее приложения. Наук. думка, Киев, 1982.
Поступила в редакцию 20 декабря 1998 г.