Пакет прикладных программ эра-dd для решения двумерных краевых задач на квазиструктурированных сетках Текст научной статьи по специальности «Математика»

УДК 519.688

Д. О. Беляев, А. Н. Козырев, В. М. Свешников

Институт вычислительной математики и математической геофизики СО РАН пр. Акад. Лаврентьева, 6, Новосибирск, 630090, Россия

E-mail: victor@lapasrv.sscc.ru

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

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

Ключевые слова: пакет прикладных программ, квазиструктурированные сетки, структуры данных, графический интерфейс, метод декомпозиции.

Введение

Необходимость в разработке пакета прикладных программ для решения уравнения Пуассона в плоских или осесимметричных областях практически произвольной формы с граничными условиями Дирихле и Неймана на квазиструктурированных несогласованных сетках возникла при расчете электрических полей в задачах сильноточной электроники, характеризующихся сильной разномасштабностью и неоднородностью потока заряженных частиц. Решение ищется итерационным методом декомпозиции расчетной области на подобласти, сопрягаемые без наложения [1. С. 99]. Граница расчетной области аппроксимируется отрезками прямых или дугами окружностей. Разработан специальный графический интерфейс [2. P. 11], дающий возможность графического и числового ввода информации о геометрии расчетной области, краевых условиях и квазиструктурированной сетке. Вблизи границы исходная прямоугольная сетка подвергается локальной модификации, состоящей в сдвиге приграничных узлов в точки пересечения координатных линий с границей, а в подобластях с сильной неоднородностью решения имеется возможность локального увеличения числа узлов. Тем самым обеспечивается адаптивность сетки. Исходное дифференциальное уравнение в подобластях аппроксимируется методом конечных объемов [3], а уравнение Пуанкаре -Стеклова на границе сопряжения подобластей - системой линейных алгебраических уравнений [4. С. 106], которая решается методом сопряженных градиентов. В статье приводится описание структур данных для работы с квазиструктурированной сеткой пакета прикладных программ и рассматриваются примеры численных расчетов.

Постановка задачи и алгоритм ее решения

Пусть в замкнутой области G = G U Г с границей Г = Г1 U Г2 требуется решить краевую задачу

ди

А и = gi, и\г= g 2,

11 дп

= gз- (1)

3 Г2

Здесь и = и(Т) - искомая функция; g1 = g1(T), g2 = g2(Т), g3 = g3(Т) - заданные функции (Т = (х, у) - текущая точка, где х, у - декартовы или цилиндрические координаты); А -оператор Лапласа; п - нормаль к Г. Будем предполагать, что Г2 состоит из отрезков, параллельных координатным осям.

ISSN 1818-7900. Вестник НГУ. Серия: Информационные технологии. 2010. Том 8, выпуск 1 © Д. О. Беляев, В. М. Козырев, В. М. Свешников, 2010

Проведем декомпозицию расчетной области G. Для этого опишем вокруг нее прямоугольник R = ¡0 < х < Dx 0 < y < Dy }, где Dx, Dy заданы (G с R ). Построим в R прямоугольную равномерную макросетку

qh = |х7 = ihx,yj = JHy, I = J = 0Ny, Hx = ndl,Hy = ,

где Nx, Ny - заданные целые числа, с шагами Ях, Я >> Л (Л - максимальный шаг сетки, на которой аппроксимируется задача (1)). При этом область G разбивается на подобласти Gk . Среди них будем различать внутренние G^, граничные G® и внешние подобласти. Замыкание G^ внутренних подобластей состоит только из точек G, граничных G(2) - из точек G и Г, а внешние подобласти исключаются из расчетов. Граница сопряжения подобластей, которую мы обозначим через у, состоит из отрезков координатных линий макросетки Qя . Точки пересечения ее координатных линий назовем макроузлами. В подобластях Gk построим равномерные прямоугольные подсетки

Qh,k = К = Х/ + ikKk , yjt = YJ + Jkhy,k , ik = 0 nx,k , Jk = 0 ny,k }

X - X Y - Y с шагами hx k = —111--, hyk = -J±1-— , причем будем предполагать без существенного

, nxk ' nyk

ограничения общности, что nxk = 2mx,k, nyk = 2my,k, где mxk , myk > 0 - целые числа. В целях

адаптации в граничных подобластях проведем локальную модификацию подсеток, состоящую в сдвиге приграничных узлов, отстоящих от границы на расстояние, меньшее половины шага сетки, в ближайшие точки пересечения координатных линий сетки с границей. Объединение подсеток образует результирующую квазиструктурированную сетку Qh. Краевая задача (1) методом конечных объемов с точностью o(h2) (h - максимальный шаг результирующей сетки) заменяется на Qh приближенной задачей [3]:

Ahuh = g^ uh|r = g 2 , lhU h |Г = ^ uhl= vh , iii 112 iy

где uh - приближенные значения функции u, а A h, lh - аппроксимации оператора Лапласа в расчетной области и оператора нормальной производной на Г; vh - след uh на у . Ее решение сводится к решению следующих дискретных подзадач:

A? uhk) = gi(k), uhk)| = vhk); (2)

'У k

Ak uhk) = gi(k), uhk)|rik = g2k), lhk) uhk) Г = g3k), uhk)| = vhk) (3)

r2,k ' k

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

- ^ и-, - Р{зик)м^+РЗ(их=/и, (4)

где (/, -) - номер узла подсетки, а индекс к мы опускаем. Для построения сеточных уравнений методом конечных объемов рассмотренная сетка подвергается триангуляции, состоящей в том, что четырехугольные элементы разбиваются кратчайшей диагональю на треугольники. В каждом треугольнике проводятся «серединные» перпендикуляры к сторонам. Вокруг каждого внутреннего сеточного узла строится замкнутая элементарная ячейка Дирихле - Вороного (множество точек, лежащих к данному узлу, ближе, чем к остальным), граница которой состоит из отрезков «серединных» перпендикуляров. Исходное дифференциальное уравнение интегрируется по ячейке, и в полученных соотношениях интегралы заменяются приближенными выражениями, что и дает уравнения (4). Более подробно алгоритм вычисления коэффициентов уравнений (4) приведен в работе [5. С. 51].

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

На границе сопряжения введем сетку , состоящую из узлов подсеток О/) и не содержащую макроузлов. Пусть N - число узлов . В работе [1. С. 99] показано, что функцию Ук можно найти из решения системы линейных алгебраических уравнений вида

ЛУ, + Ь = 0, (5)

_ _ А' V = gh. _(6)

Здесь Л = , г = 1, N, j = 1, N + р} - прямоугольная матрица; Ь = Ь, г = 1, N}, gh = ., г = N +1, N + р}, у/ , * = 1, N + р} - векторы (Р - число макроузлов);

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

Аh Vh = gl, Vh= g2, К Vh|Г= g3, Vh|т= 0

относительно функции V h.

Элементы матрицы Л явно не вычисляются. Решение системы (5), (6) проводится методом сопряженных градиентов по формулам

г (0)= Ь - Лу^ , р (0)= г(0); (7)

>+!)=,,(") +а Л") а(") = (- (»)-(») (п) „(»)'

= ук") + а "Р(п), а(") = (г(п), г(п))/(Лр(п), р(п)); (8)

г ("+1) = г (п)-а пЛр(п), р(п+1) = г(п+1) +р пр(п); (9)

р" =(г("+1),г("+1))/(г(п),г(п)), " = 0,1,..., (10)

где Л - квадратная матрица, состоящая из прямоугольных блоков Л, А(р); Ь - вектор, состоящий из подвекторов - Ь, gh. Вычисление Л^, где под wh понимаются векторы у/0,

/г

р(п), фигурирующие в приведенных ранее формулах, сводится к расчету произведений Лм> и А(р)wh. Рассмотрим первое из них. Решим краевую задачу

Аh фh = gl, Фh= g2, К Ф/|г= gз, Ф/|т= wh.. (11)

Разумеется, данная краевая задача сводится к решению подзадач вида (5), (6). Тогда искомое произведение равно

= sh - Ь,

где sh = .., г = 1, N} - вектор, компоненты которого определяются как

= (й(+Ч+)) + (й/-Ч~)), * = . (12)

Здесь - операторы, аппроксимирующие внутренние нормальные производные на границе сопряжения по формулам

(, ) - 3(ф h )0 + 4(ф h )1 -(Ф h )2

(й/Ф h )0 =-2й-'

где 0, 1, 2 - соседние узлы, разделенные шагом /. Если узел 2 модифицированный, то используется простейшая двухточечная схема, связывающая узлы 0, 1. Вычисление А(р)у/") не вызывает затруднений, так как элементы матрицы А(р известны.

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

АнУн = ^ Т;) Ун 1г= g2(тр), 1нУн|Г= g3(тр),

где A H , lH - аппроксимации оператора Лапласа и оператора граничных условий на данной макросетке; vH - искомая сеточная функция; T' - узлы макросетки Q'H . Функция vH затем интерполируется в узлы сетки ю h на границе сопряжения и макроузлы T'p, что дает функцию

vh0).

Вычислив значения искомой функции на границе сопряжения, функции u f) в подобластях найдем из решения задач (2), (3).

Структура данных для квазиструктурированной сетки

Для каждой k -й внутренней или граничной подобласти хранятся все узлы структурированной подсетки Qhk, включая внешние узлы. Среди узлов, принадлежащих расчетной

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

Массивы для сетки юл на границе сопряжения специально не заводятся, так как данная сетка по сути дела состоит из узлов, являющихся проекциями на границу сопряжения подсе-ток Q h,k.

Приведем описание хранимых скаляров и массивов, придерживаясь фортрановских правил записи массивов и вводя локальную нумерацию.

1. Внешние данные, которые формируются графическим препроцессором.

1.1. Скаляры.

1.1.1. nx, ny - число интервалов макросетки Qя по х, y , т. е. фактически число подобластей в соответствующих направлениях.

1.1.2. xb, yb, xe, ye - координаты левого нижнего и правого верхнего углов прямоугольника, описанного вокруг расчетной области G .

1.2. Массив.

1.2.1. NNSG(nx, ny, 2) - число интервалов подсеток Qhk по x, y.

2. Внутренние данные, которые формируются при обработке квазиструктурированной

сетки.

2.1. Скаляры.

2.1.1. nngi - суммарное число узлов в подобластях, которые участвуют в расчетах. Входят все узлы подобластей. Каждая подобласть представляется как отдельный прямоугольник, для которого хранятся все узлы, включая стороны прямоугольника и макроузлы. На границах сопряжения двух подобластей узлы хранятся для каждой из них, т. е. совпадающие узлы повторяются.

2.1.2. nlmn - число модифицированных узлов.

2.1.3. ncoe1 - длина массива COE1 - коэффициентов сеточных уравнений в подобластях.

2.1.4. ncoe3 - длина массива COE3 - коэффициентов сеточных уравнений в макроузлах.

2.2. Массивы для подсеток в расчетной области.

2.2.1. ITSG(nx, ny) - тип подобласти (-1 - внешняя, 0 - внутренняя, 1 - граничная).

2.2.2. IANG(nx, ny) - адреса начал подмассивов узлов в подобластях: (для массивов ITNG, RUNG).

2.2.3. ITNG(nng7') - тип узлов в расчетной области (-1 - внешний, 0 - внутренний, > 0 - модифицированный, причем данное число - адрес координат в массиве COLM (см. далее)).

2.2.4. RUNG(nngi) - значения искомой функции в узлах в расчетной области (те же узлы, что и в массиве ITNG).

2.2.5. IACO(nngi) - вспомогательный массив для работы итерационных модулей, вычисляющих значения искомой функции в подобластях (-1 - внешний узел, 0 -узел со стандартными коэффициентами, >0 - адрес коэффициентов в массиве COE1). Если узел лежит на границе с условием Дирихле, то он в массиве IACO имеет признак -1, т. е. формально относится к внешнему узлу. Значение функции в нем не вычисляется в итерационном процессе, а просто заносится при расчете коэффициентов. Если же узел лежит на границе с условием Неймана, то его признак в массиве IACO - число > 0 и функция в нем рассчитывается.

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

2.2.7. COE2(nx, ny, 3) стандартные коэффициенты сеточных уравнений в подобластях (для каждого узла - 3 коэффициента: pf], p(1j, pi2] в формуле (4)).

2.2.8. COLM(3*n/mn) - координаты модифицированных узлов. Для каждого узла хранится х, у, ig. Последнее число - номер границы, на которой лежит рассматриваемый узел.

2.3. Массивы для макроузлов.

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

2.3.1. ITMN(nx + 1, ny + 1) - тип макроузлов в расчетной области (-1 - внешний, 0 -внутренний, > 0 - модифицированный, причем данное число - адрес координат в массиве COLM).

2.3.2. RUMN(nx + 1, ny + 1) - значения искомой функции в макроузлах.

2.3.3. IAMN(nx + 1, ny + 1) - вспомогательный массив для вычисления функции в макроузлах (-1 - внешний, 0 - узел со стандартными коэффициентами, > 0 - адрес коэффициентов в массиве COE3 для данного узла). Аналог массива IACO, поэтому к нему относятся те же замечания, что и к IACO.

2.3.4. COE3(ncoe3) - коэффициенты сеточных уравнений для макроузлов, лежащих вблизи границы расчетной области (аналог COE1).

Структура пакета

Пакет состоит из макромодулей Graf, Geom, CaCoe, SolverU, Service, которые работают в приведенной последовательности, записывая результат работы в файлы. Макромодуль Graf является препроцессором пакета, макромодули Geom, CaCoe, SolverU образуют процессор, а макромодуль Service служит постпроцессором. Приведем краткое описание макромодулей.

Макромодуль Graf служит для задания пользователем входной информации средствами графического интерфейса. Выходная информация записывается в следующие файлы:

• era.dat, point.dat - геометрия области и граничные условия;

• fqgrid.dat - квазиструктурированная сетка.

Макромодуль Geom служит для обработки геометрии расчетной области. В его функции входит выполнение следующих процедур: построение сеток, локальная модификация сеток, разметка узлов в подобластях. Исходные данные - информация о геометрии расчетной области, граничных условиях и квазиструктурированной сетке, хранящаяся в файлах era.dat, point.dat, fqgrid.dat. Выходные данные:

• скаляры, которые записываются в файл fscalar.dat;

• массивы для узлов в подобластях: ITSG, IANG, ITNG;

• массив координат локально модифицированных узлов COLM.

Все массивы записываются в файл fsgeom.dat.

Макромодуль CaCoe проводит расчет коэффициентов сеточных уравнений в подобластях. Коэффициенты хранятся для каждого околограничного узла, а также для каждого узла, лежащего на границе с условием Неймана. Подчеркнем, что для граничных узлов с условием Неймана строятся сеточные уравнения и в них рассчитываются значения искомой функции, а для граничных узлов с условием Дирихле в функции данного макромодуля входит засылка заданных на внешней границе значений функции в массив RUNG. Исходные данные содержатся в файлах era.dat, point.dat, fqgrid.dat, fscalar.dat, fgeom.dat. Выходные данные:

• скаляры в файле fscalar.dat, дополненные новыми значениями;

• массивы COE1, COE2;

• массив RUNG;

• вспомогательный массив IACO.

Массив COE1 записывается в файл fcoe1.dat, массив COE2 - в файл fcoe2.dat, массив RUNG - в файл fu.dat, а массивы IACO - в файл fassist.dat.

Макромодуль SolverU реализует итерационный процесс по подобластям, а макромодуль Service служит для изображения квазиструктурированной сетки и вывода изолиний решения.

Графический интерфейс

для автоматизации работы с исходными данными

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

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

Графическое окно визуализации предназначено для отображения геометрии расчетной области.

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

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

Диалоговые окна предназначены для ввода / вывода той или иной информации о задаче. Они состоят из заголовка, кнопок управления окном (свернуть, восстановить, закрыть окно), поля ввода / вывода, кнопок подтверждения выполненных в окне операций.

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

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

Для ввода геометрического объекта имеется три способа: 1) графический; 2) числовой; 3) графико-числовой.

Рассмотрим вначале ввод отрезка прямой. Самый простой способ ввода - графический. Он состоит в том, что на панели инструментов активизируется соответствующая кнопка,

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

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

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

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

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

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

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

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

Ввод квазиструктурированной сетки проводится в специальном диалоговом окне и включает два этапа: 1) построение прямоугольной макросетки Он , определяющей декомпозицию расчетной области; 2) построение подсеток Окк в подобластях, образовавшихся на первом этапе.

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

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

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

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

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

Примеры численных расчетов

Рассматривалась следующая тестовая задача:

А и = 0, и|г= g (Т), ^ 111 дп

= 0,

Г2

о расчете электрического поля в четверти цилиндрического конденсатора. Здесь Т - текущая точка; Г1 = Г11 и Г12, где Г11, Г12 - четверти концентрических окружностей с радиусами

Я1 = 0,1 и Я2 = 1; Г2 - прямолинейные отрезки, параллельные координатным осям,

Г 0при Т еГ11,

g [ 1приТ еГ12

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

Рис. 1. Квазиструктурированная сетка Рис. 2. Картина электрического поля

в цилиндрическом конденсаторе в цилиндрическом конденсаторе

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

Рис. 3. Квазиструктурированная сетка Рис. 4. Картина электрического поля

в электронной пушке в электронной пушке

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

1. Свешников В. М. Построение прямых и итерационных методов декомпозиции // Сиб-ЖИМ. 2009. Т. 12, № 3 (39). С. 99-109.

2. Kozyrev A. N., Sveshnikov V. M. Automation of initial data input in program package ERA-DD for the solving of two-dimensional boundary value problems on quasistructured grids // Bull. of the Novosibirsk Computing Center. Ser. Numerical Analysis. Novosibirsk: ICMMG. I. 2009. Vol. 14. P. 11-23.

3. Ильин В. П. Методы конечных разностей и конечных объемов для эллиптических уравнений. Новосибирск: ИВМиМГ (ВЦ) СО РАН, 2001. 318 с.

4. Свешников В. М. Прямой метод декомпозиции без наложения подобластей для решения краевых задач на прямоугольных квазиструктурированных сетках // Вычисл. технологии. 2008. Т. 13, № 2. С. 106-118.

5. Ильин В. П., Павлов М. В, Свешников В. М. Решение двумерных краевых задач на квазиструктурированных сетках // Тр. Междунар. конф. RDAMM-2001. Новосибирск, 2001. Спецвыпуск. Т. 6, ч. 2. С. 51-59.

Материал поступил в редколлегию 10.12.2009

D. O. Belyajev, A. N. Kozyrev, V. M. Sveshnikov

PROGRAM PACKAGE ERA-DD FOR THE SOLVIN TWO-DIMENSIONAL BOUNDARY VALUE PROBLEMS ON QUASISTRUCTURED GRIDS

The program package for the solving boundary value problems on quasistructured not co-ordinated grids by the iterative domain decomposition method on the subdomains interfaced without overlapping is developed. The description of computational algorithm, the data structures for work with quasistructured grids is resulted. The graphic interface for automation of initial data input and examples of numerical calculations are given.

Keywords: program package, quasistructured grids, data structures, graphic interface, decomposition method.