CC BY
98
7
МАТЕМАТИЧЕСКОЕ И КОМПЬЮТЕРНОЕ
МОДЕЛИРОВАНИЕ _MODELING AND SIMULATION
УДК 532.529
МНОГОСЕТОЧНЫЙ МЕТОД УСКОРЕНИЯ СХОДИМОСТИ ПРИ РЕШЕНИИ ЗАДАЧ ГАЗОВОЙ ДИНАМИКИ НА НЕСТРУКТУРИРОВАННЫХ СЕТКАХ
К.Н. Волков^ b
a Университет ИТМО, Санкт-Петербург, Россия
b Университет Кингстона, Лондон, Великобритания, k.volkov@kingston.ac.uk
Аннотация. Разработан подход к эффективной реализации численных методов и вычислительных алгоритмов, предназначенных для моделирования течений невязкого и вязкого сжимаемого газа в областях сложной геометрической конфигурации. Развивается метод ускорения сходимости итерационного процесса, основанный на использовании многосеточных технологий решения систем разностных уравнений большой размерности. Уравнения Эйлера и Навье-Стокса дискретизируются на структурированных и неструктурированных сетках при помощи схем высокой разрешающей способности по времени и по пространству. Многосеточный метод реализуется на основе стандартного C/F разбиения переменных и метода стандартной интерполяции. Предложен подход, позволяющий преодолеть проблемы, связанные с хранением коэффициентов матрицы разных знаков. Возможности разработанных средств численного моделирования демонстрируются на примере задачи обтекания профиля. Приведены результаты расчетов течений невязкого и вязкого газа на структурированных и неструктурированных сетках при использовании различных составных компонентов многосеточной технологии. Сравнение скорости вычислений и факторов сходимости на структурированных и неструктурированных сетках показало экономичность разработанного подхода, а также слабую зависимость показателей качества от числа узлов сетки. Выбор оптимального числа внутренних итераций позволил сократить время счета на 10%. Результаты работы могут быть использованы при разработке вычислительных пакетов, предназначенных для решения задач газовой динамики в областях сложной формы. Ключевые слова: вычислительная газовая динамика, неструктурированная сетка, многосеточный метод, сходимость, ускорение.
MULTI-GRID METHOD OF CONVERGENCE SPEEDING-UP FOR THE SOLUTION OF GAS DYNAMICS PROBLEMS ON UNSTRUCTURED MESHES
K.N. Volkova b
a ITMO University, Saint Petersburg, Russia
b Kingston University, London, United Kingdom, k.volkov@kingston.ac.uk
Abstract. An approach for effective implementation of numerical methods and computational algorithms is developed to simulate flows of non-viscous and viscous compressible gas in the complex domains. The convergence speeding-up method of iterative process is discussed based on the usage of multi-grid technologies to the solution of large systems of finite difference equations. Euler and Navier-Stokes equations are quantized on structured and unstructured meshes with highresolution schemes in time and space. Multi-grid method is implemented on the basis of standard C/F splitting of variables and standard interpolation method. A specific approach is proposed to overcome the problems related to storing the matrix coefficients of different signs. A problem of flow around an airfoil is used to demonstrate the possibilities of computational tools developed. The results computed for non-viscous and viscous gas on structured and unstructured meshes are presented with the usage of various components of multi-grid technology. Comparison of computation speed and convergence factors on structured and unstructured meshes showed the economy of the developed approach and weak dependence of quality characteristics on the number of mesh points. The choice of optimal inner iterations number made it possible to decrease computation time for 10%. Obtained results can be applied while software creation for fluid dynamics problems in the complex configuration domains.
Keywords: computational fluid dynamics, unstructured mesh, multi-grid method, convergence, speeding-up.
Введение
Для расчетов течений жидкости и газа в областях сложной формы, характерных для современных технических устройств и технологических приложений, широкое применение находит метод конечных объемов на неструктурированных сетках. Расчеты течений жидкости и газа на неструктурированных сетках сталкиваются с рядом проблем, связанных с медленной скоростью сходимости итерационного процесса и существенными затратами процессорного времени, необходимого для получения сеточно-независимого и сходящегося решения.
Дискретизация уравнений Эйлера или Навье-Стокса порождает систему разностных уравнений с разреженной матрицей большого размера, которая решается тем или иным методом. При дискретизации уравнений на неструктурированных сетках матрица коэффициентов является несимметричной и не имеет диагонального преобладания. В этом случае классические итерационные методы либо перестают работать, либо приводят к медленной скорости сходимости [1, 2]. Проблемы, возникающие при численном
решении уравнений Эйлера и Навье-Стокса и приводящие к медленной сходимости, показаны на рис. 1. (символы Е и N8 соответствуют уравнениям Эйлера и уравнениям Навье-Стокса). При решении уравнений Эйлера дискретная жесткость является следствием использования скалярного шага по времени. Основной источник трудностей решения уравнений Навье-Стокса состоит в использовании сеток с растянутыми и скошенными ячейками в пограничном слое [3]. Дополнительные трудности возникают при моделировании турбулентных течений. Для улучшения сходимости используется матричный шаг по времени и метод направленного огрубления сетки [4, 5].
E, NS
NS
NS
E, NS
NS
Дискретная жесткость
Направленные связи между переменными
E, NS
Численные проблемы
E, NS NS
Медленная сходимость
Рис. 1. Причины медленной сходимости итерационного процесса
Имеются два подхода к ускорению газодинамических расчетов, один из которых связан с ускорением сходимости [6, 7] (применение многосеточных технологий, методов предобусловливания и стабилизации численного решения), а другой - с использованием средств, предоставляемых современными параллельными системами обработки данных, включая графические процессоры общего назначения [8].
В настоящей работе обсуждается реализация и применение алгебраического многосеточного метода для решения задач газовой динамики, описываемых уравнениями Эйлера и Навье-Стокса. В отличие от геометрического подхода [5], реализация алгебраического многосеточного метода не требует информации о геометрии решаемой задачи, используя лишь структуру матрицы коэффициентов системы разностных уравнений.
В многосеточном методе обычный итерационный процесс комбинируется с коррекцией решения на последовательности грубых сеток. Одна итерация многосеточного метода (многосеточный цикл) включает в себя предварительное сглаживание (рге-БшооШ!^); расчет невязки на текущем уровне сетки; ограничение и коррекцию невязки на грубой сетке; продолжение и интерполяцию ошибки на подробную сетку; коррекцию решения на подробной сетке с использованием поправки, интерполируемой с грубой сетки; заключительное сглаживание для погашения высокочастотных компонент погрешности, появляющихся после интерполяции на подробную сетку.
Многосеточный метод устанавливает лишь структуру вычислительного алгоритма, эффективность которого зависит от адаптации его отдельных компонентов к решаемой задаче [9]. Многообразие имеющихся подходов связывается с поиском компромисса между общей эффективностью алгоритма и объемом требуемой вычислительной работы. Имеющиеся реализации многосеточных методов различаются способом построения операторов ограничения и продолжения, выбором сглаживающей процедуры, реализацией многосеточного цикла, способом построения последовательности сеток. В отличие от имеющихся реализаций, в разработанном подходе используется метод, позволяющий преодолеть проблемы, связанные с наличием как положительных, так и отрицательных связей между переменными. Возможности подхода демонстрируются на примере расчета обтекания профиля на структурированных и неструктурированных сетках.
Реализация многосеточного метода
Конечно-объемная дискретизация уравнений Эйлера или Навье-Стокса на фиксированном сеточном уровне приводит к решению системы разностных уравнений
аУ = f". (1)
Матрица А размером п*п состоит из коэффициентов ау (/,у=1, ..., п), связанных с дискретизацией. Вектор и длиной п состоит из неизвестных узловых значений искомой функции. Вектор f длиной п составлен из коэффициентов, обусловленных дискретизацией, а также из известных значений вектора решения, задаваемых граничными условиями.
Построение сеточных уровней. Вместо матрично-векторной терминологии используются термины, связанные с последовательностью фиктивных сеток, под узлами которых понимаются узлы направ-
ленного графа, соответствующего данной матрице. Узел сетки i ей' (имеет смысл переменной и') связан с переменной j ей', если а' Ф 0. Множество переменных, соседних с переменной i, имеет вид
N' ={j ей' : j Ф i, ah Ф 0} (i ей').
Используя сеточную терминологию, уравнение A^u' = f' можно рассматривать как фиктивное сеточное уравнение на подробной сетке Q', а уравнение AHuH = f - как фиктивное сеточное уравнение на грубой сетке йH с й'. Под Q' = {1, 2, ..., n} понимается множество индексов неизвестных. Множество Qlh имеет меньшее число элементов и соответствует грубой сетке.
Переменные на грубом сеточном уровне строятся как некоторое подмножество переменных, принадлежащих подробной сетке [10-12]. Множество переменных подробной сетки разделяется на два несвязанных подмножества й' = C' U F', где C' - множество переменных, которые принадлежат грубой сетке (C-переменные), - множество переменных, которые принадлежат подробной сетке (F-переменные).
Сеточному уровню Qk (k = 1, 2, ..., M-1) соответствуют подмножества Ck и Каждому сеточному уровню соответствуют сеточные операторы A1, A2, ., AM (при этом A1=A), а также операторы интерполяции Pk = Ik+1 и ограничения Rk = Ik+1, где к = 1, 2, ..., M-1. Оператор ограничения представляет собой транспонированный оператор интерполяции Rk = (Pk)' (под ' понимается операция транспонирования). Для построения матрицы системы на грубом сеточном уровне вычисляется произведение Галеркина RkAkPk. Сглаживание погрешности на каждом сеточном уровне производится при помощи сглаживающего оператора Sk (k = 1, 2, ..., M-1). Процесс огрубления повторяется до тех пор, пока размер системы не уменьшается до уровня, когда становится эффективным прямой метод ее решения (например, метод исключения Гаусса).
Для выбора матричных коэффициентов, используемых для построения грубого сеточного уровня, вводятся понятия сильного влияния переменных и сильной зависимости между переменными [11, 12]. Переменная i сильно зависит от переменной j, а переменная j имеет сильное влияние на переменную i, если матричный коэффициент aj по величине превосходит все внедиагональные коэффициенты матрицы (standard coarsening)
-atj >6max(~ая). (2)
k Ф1
Параметр 0 < 9 < 1 в соотношении (2) контролирует число сильных связей между переменными (обычно 9 = 0,25). Множество Si представляет собой множество всех переменных j, имеющих сильные связи с переменной i (множество переменных, от которых переменная i сильно зависит):
S , = {у: J' * i, - a.j > 9max (~aik )J .
Имеющиеся теоретические подходы относятся к решению скалярных эллиптических дифференциальных уравнений в частных производных [13], дискретизация которых приводит к системе разностных уравнений с М-матрицей (в этом случае скорость сходимости метода не зависит от размера задачи). При наличии как отрицательных, так и положительных коэффициентов, стоящих вне главной диагонали, используются следующие определения:
(а.., {{ а.. < 0, Г 0, 15 а.. < 0,
а" =\ . . а+ =\ у
" I 0, 15 а. > 0, . | а., 15 а. > 0. . . .
При этом имеется два несвязанных подмножества переменных
е N : аI < 0}, м;=у е N. : а* > 0} .
Наличие связей между переменными различного знака в уравнении (1) усложняет реализацию метода и приводит к замедлению сходимости.
Интерполяция на подробную сетку. Интерполяция погрешности решения eh = \кнeH определяется в виде
^ Н^ ^, 15 , е Ек (3)
где Рк с Ск представляет собой множество переменных, участвующих в интерполяции. При интерполяции с грубой на подробную сетку погрешность решения е* полагается равной ен , если переменная , является С-переменной, и взвешенной сумме переменных из множества Р* , если переменная принадлежит подмножеству Е. Множество Рк является малым подмножеством множества С-переменных, находящихся в некоторой близости от переменной , , что гарантирует разреженность матрицы Aн. С другой
nh
стороны, множество р содержит достаточное число переменных, с которыми переменная i является сильно связанной.
Для простоты верхний индекс h опускается, а вместо (3) используется соотношение
=Ё wkek (l еF). (4)
кер
Локальное приближение для алгебраически гладкой погрешности решения имеет вид
e =- — S aijej (i eQ). (5)
aii jeN,
Реализация метода заключается в построении подходящего C/F разбиения, а также в определении множества интерполяционных переменных р с C (i е F) и весовых множителей wlk таких, что соотношение (4) дает приемлемое приближение для любой алгебраически гладкой погрешности.
Точная процедура интерполяции дается соотношением (5), полагая P, = Ni и wk = -aik/aii. Подход требует построения такого C/F разбиения, что для каждой переменной i е F все переменные, соседние с переменной i, содержались бы в подмножестве C, и приводит к построению прямого метода решения системы разностных уравнений [11, 12]. Для построения эффективной вычислительной процедуры требуется построение малых подмножеств P, чтобы обеспечить разреженность оператора Галеркина и достичь высокой скорости сходимости.
Для прямой интерполяции (direct interpolation) р с Nt, а оператор интерполяции строится из соотношения (5) при известном подмножестве переменных, не участвующих в интерполяции (часть суммы, соответствующая алгебраически гладкой ошибке, для которой j е Nj \ р).
Для матриц общего вида используются коэффициенты at и a- и подмножества N* и N- , а положительные коэффициенты добавляются к коэффициентам, стоящим на главной диагонали. Процедура интерполяции определяется соотношением
+ a
il l
keP,
£ alk^k = 0, (6)
где
7 a..
•чт^ + jeN, 4
aii = aii + 7 ak , a, = V-- •
jeN 7 k£pa*
Такой подход приводит к положительным весовым коэффициентам
wik =-a, aaL (i e F, k e p ).
a..
В общем случае при построении процедуры интерполяции учет положительных связей в соотношении (6) оказывается более важным, чем учет отрицательных связей [11, 12] (для положительных связей погрешность решения изменяется наиболее медленно). Для некоторой переменной i e F, имеющей
как отрицательные, так и положительные связи, предполагается, что N- Ф0 и N+ Ф0 , а C/F разбиение является таким, что как минимум одна связь любого знака принадлежит подмножеству C. Выбираются два множества интерполяционных переменных таких, что 0 = P- с C П N- и 0 = P+ с C П N .
Для улучшения процедуры интерполяции предполагается, что имеется такое C/F разбиение, при котором для каждой переменной i e F множество интерполяционных переменных p с Nt П N+ является данным. Вместо непосредственного приближения неинтерполяционной части уравнения для переменной i в соотношении (5), производится исключение погрешностей решения е. для всех переменных j g P,, имеющих сильные связи с переменной i в уравнении j. Стандартная интерполяция (standard interpolation) приводит к расширению множества интерполяционных переменных (уменьшение разреженности оператора Галеркина), но оказывается достаточно эффективной на практике.
Сглаживание. Применение сглаживающей процедуры общего вида дает
ии+, = SuB + (I - S)A-1f . (7)
В качестве сглаживающего оператора выбирается матрица S = I - Q1 A, где Q - нижняя треугольная часть матрицы A. В качестве сглаживающей процедуры в соотношении (7) широкое применение находит метод Гаусса-Зейделя.
Показатели эффективности. Сходимость многосеточного метода зависит от изменения нормы погрешности решения
IIе t =(A 4eh, eh Г
Скорость сходимости характеризуется уменьшением нормы невязки на двух соседних циклах р = ||гЛ+1||/||гЛ||. Фактор сходимости получается при помощи осреднения скорости сходимости по всем многосеточным циклам (при достижении заданного уровня невязки) или при помощи осреднения по заданному числу циклов, если заданный уровень невязки в расчетах не достигается. Фактор асимптотической сходимости рассчитывается по отношению к норме невязки на первом многосеточном цикле.
Фактор сложности выражается в виде числа операций, которые необходимо выполнить на одной итерации, и требуемой для этого памяти. Операторная сложность (operator complexity) служит индикатором памяти, требуемой для реализации многосеточного метода
Г V
° ~ V | A ®Г
где |А(к)| - число ненулевых коэффициентов на уровне к.
Фактор сложности сетки (grid complexity) показывает, насколько быстро производится огрубление
сетки:
rG =V ^,
к П
где nk - число переменных (узлов сетки) на уровне к.
Сложность многосеточного цикла (cycle complexity) характеризуется параметром
| A(k) | п V | A L „к
ГГ = V jA^ У ,
где дк = ц1 + д2. Параметр у соответствует типу многосеточного цикла (у = 1 для V-цикла и у = 2 для W-цикла). Объединение фактора сходимости и сложности многосеточного цикла дает объем вычислительной работы, приходящейся на единицу точности, rW = -rr/log р.
Результаты расчетов
Рассмотрим невязкое и вязкое обтекание профиля NACA0012 при числе Маха невозмущенного потока Мм = 0,8 и при различных углах атаки (а = 0° и а = 1,25°). Расчетная область приводится на рис. 2. На входной границе задаются условия невозмущенного течения, а на выходной границе - условия свободного вытекания. На верхней и нижней границах расчетной области используются условия скольжения. При этом R = 10Г, H = 10Г и L = 30Г, где Г - хорда профиля.
Рис. 2. Расчетная область
Многосеточная процедура реализуется в рамках У-цикла (^ = 2 и ц2 = 1)- Для интегрирования по времени используется трех- (КК3) и пятишаговый (КК5) метод Рунге-Кутты, а для дискретизации невязких и вязких потоков - схема Ми8СЬ и центрированная разностная схема. Детали конечно-объемной дискретизации обсуждаются в работе [14]. Для построения С/Е разбиения применяется метод стандартного огрубления. В качестве интерполяции используется метод стандартной интерполяции, а для сглаживания - метод Гаусса-Зейделя. Для вычисления произведения Галеркина, хранения разреженных матриц и реализации операций матричной алгебры используется формат БЬЬРЛСК.
Расчеты проводятся как на структурированной сетке (вариант 1, сетка типа О), так и на неструктурированной сетке с треугольными ячейками (вариант 2) различной разрешающей способности. Расчетные сетки приводятся на рис. 3 (показана часть расчетной области). При решении вязкой задачи узлы сетки сгущаются около поверхности профиля для надлежащего разрешения пограничного слоя. Максимальное отношение сторон ячеек сетки составляет 500.
а б
Рис. 3. Структурированная (а, вариант 1) и неструктурированная (б, вариант 2) сетки
Для невязкой задачи распределения коэффициента давления по поверхности профиля показаны на рис. 4 и имеют практически одинаковый вид как на структурированной, так и на неструктурированной сетках. Решение дает сильный скачок на верхней поверхности профиля (при xlL = 0,68) и слабый скачок на его нижней поверхности (при xlL = 0,32). Для условий задачи интегрирование давления по поверхности профиля дает коэффициенты трения и подъемной силы (^ = 0,0227 и Cy = 0,3527), хорошо согласующиеся с данными [15].
1,2
0,4
CJ I
-0,4
-1,2
0 0,5 1
xlL
Рис. 4. Распределения коэффициента давления по поверхности профиля, полученные на структурированной (сплошная линия) и неструктурированной (значки □) сетках
История сходимости на структурированной и неструктурированной сетках показана на рис. 5. В отличие от осциллирующей сходимости на структурированной сетке, решение задачи на неструктурированной сетке дает гладкое изменение невязки. Уровень сходимости (минимальный уровень невязки) примерно одинаковый в обоих случаях.
Информация о сходимости многосеточной процедуры приводится в табл 1. Качество огрубления сетки характеризуется параметрами, представляющими собой общее число ячеек и узлов сеток различной разрешающей способности (параметры CGv и CGC). Зависимость скорости сходимости от разрешения сетки относительно слаба.
Моделирование вязкого обтекания профиля проводится при числе Рейнольдса Яе=500 (число Рей-нольдса рассчитывается по хорде профиля). Решение вязкой задачи дает распределение коэффициента давления по поверхности профиля, слабо отличающееся от решения, показанного на рис. 4. Коэффициенты
сопротивления и подъемной силы принимают значения Сх = 0,0225 и Су = 0,3536. Относительное различие результатов решения невязкой и вязкой задачи составляет 0,88% для коэффициента сопротивления и 0,26% для коэффициента подъемной силы. Информация о сходимости многосеточной процедуры приводится в табл. 2. Скорость сходимости при нулевом угле атаки оказывается лучше, чем при угле атаки а = 1,25°.
0 -2
-4
-6
-8
-10
-12
0 60 120 180 240
п
Рис. 5. Изменение невязки в зависимости от числа многосеточных циклов на структурированной (линия 1)
и на неструктурированной (линия 2) сетках
2
1
Размер Число уровней CGv CGc p при a=0° p при a=1,25°
RK3 RK5 RK3 RK5
2607 2 1,58 1,34 0,24 0,18 0,28 0,18
5258 3 1,64 1,36 0,34 0,24 0,25 0,17
10273 4 1,67 1,36 0,25 0,24 0,24 0,17
20621 5 1,64 1,36 0,34 0,21 0,23 0,21
Таблица 1. Фактор сходимости при решении уравнений Эйлера
Размер Число уровней CGv CGc p при a = 0° p при a = 1,25°
RK3 RK5 RK3 RK5
8872 2 1,64 1,37 0,42 0,37 0,26 0,21
18416 3 1,88 1,44 0,43 0,37 0,61 0,57
36388 4 1,67 1,37 0,41 0,35 0,52 0,45
Таблица 2. Фактор сходимости при решении уравнений Навье-Стокса
На подробной сетке применение алгебраического многосеточного метода дает фактор сходимости, равный 0,55 для уравнений Эйлера и 0,81 для уравнений Навье-Стокса. Для геометрического многосеточного метода, реализованного в работе [5], асимптотический фактор сходимости составляет 0,80 для уравнений Эйлера и 0,91 для уравнений Навье-Стокса.
В большинстве реализаций полагается, что 9 = 0,25, а числу уровней сетки присваивается максимальное возможное значение [10, 16]. На практике коэффициент 9 изменяется в достаточно широких пределах (обычно от 0,2 до 0,4). Значение 9 = 0,25 соответствует интерполяции по четырем точкам на подробной сетке (в геометрическом методе этому значению соответствует увеличение шага сетки в два раза при переходе от подробной сетки к грубой). Изменение параметра 9 от 0,2 до 0,4 дает максимальное увеличение времени счета на 7,5% на структурированной сетке. Оптимальным значением на неструктурированной сетке является 0,0625. Выбор оптимального числа внутренних итераций позволяет сократить время счета на 10% в зависимости от граничных условий задачи.
Заключение
Предложен подход к эффективной реализации многосеточной процедуры на неструктурированных сетках. К особенностям реализации относятся подход, связанный с хранением положительных и отрицательных связей между переменными, а также реализация операций с разреженными матрицами в формате ELLPACK. Полученные результаты показывают, что применительно к решению задач механики жидкости алгебраический многосеточный метод позволяет ускорить сходимость итерационного процесса, требуя сравнительно небольших модификаций программного кода. Алгебраический многосеточный метод 162 Научно-технический вестник информационных технологий, механики и оптики
162 Scientific and Technical Journal of Information Technologies, Mechanics and Optics
2014, № 4 (92)
оказывается чувствительным к изменению сетки и ее разрешению на самом грубом сеточном уровне. Оценки показывают, что выбор параметров в алгебраическом многосеточном методе позволяет сократить общее время счета на 10%. Подходящий выбор операторов продолжения и ограничения позволяет сгладить высокочастотные моды погрешности решения при переходе к более грубому сеточному уровню. Учет сильных связей между переменными в процедуре интерполяции приводит к ускорению сходимости, в то время как слабые связи оказывают негативное влияние на сходимость, увеличивая объем вычислительной работы. Зависимость подхода от свойств сглаживающей процедуры сравнительно слаба.
В отличие от геометрического подхода, алгебраический метод допускает сравнительно простую программную реализацию в рамках существующего программного кода. Дополнительной реализации и модификации требуют структуры данных, предназначенные для хранения разреженных матриц, а также реализация операций матричной алгебры с разреженными матрицами. Вместе с тем, метод построения последовательности сеток и его адаптация к решаемой задаче (невязкая/вязкая, ламинарная/турбулентная) во многом предопределяют успех применения многосеточной технологии ускорения сходимости.
Литература
1. Brandt A. Guide to multigrid development // Lecture Notes in Mathematics. 1982. V. 960. P. 220-312.
2. Wesseling P. An introduction to mutigrid methods. Chichester: John Wiley & Sons, 1992. 152 p.
3. Pierce N.A., Giles M.B., Jameson A., Martinelli L. Accelerating three-dimensional Navier-Stokes calculations // Proc. 13th Computational Fluid Dynamics Conference. 1997. P. 676-698.
4. Волков К.Н. Блочное предобусловливание уравнений Эйлера и Навье-Стокса при моделировании низкоскоростных течений на неструктурированных сетках // Журнал вычислительной математики и математической физики. 2009. Т. 49. № 10. С. 1868-1884.
5. Волков К.Н. Многосеточные технологии для решения задач газовой динамики на неструктурированных сетках // Журнал вычислительной математики и математической физики. 2010. Т. 50. № 11. С. 1938-1952.
6. Cagnone J.S., Sermeus K., Nadarajah S.K., Laurendeau E. Implicit multigrid schemes for challenging aerodynamic simulations on block-structured grids // Computers and Fluids. 2011. V. 44. N 1. P. 314-327.
7. Colin Y., Deniau H., Boussuge J.-F. A robust low speed preconditioning formulation for viscous flow computations // Computers and Fluids. 2011. V. 47. N 1. P. 1-15.
8. Brodtkorb A.R., Hagen T.R., Satra M.L. Graphics processing unit (GPU) programming strategies and trends in GPU computing // Journal of Parallel and Distributed Computing. 2013. V. 73. N 1. P. 4-13.
9. Cleary A.J., Falgout R.D., Henson V.E., Jones J.E., Manteuffel T.A., McCormick S.F., Miranda G.N., Ruge J.W. Robustness and scalability of algebraic multigrid // SIAM Journal on Scientific Computing. 2000. V. 21. N 5. P. 1886-1908.
10. Ruge J., Stuben K. Algebraic multigrid (AMG) / In:Multigrid Methods (Frontiers in Applied Mathematics) / Ed. S.F. McCormick. Philadelphia: SIAM, 1987. V. 3. P. 73-130.
11. Stuben K. A review of algebraic multigrid // Journal of Computational and Applied Mathematics. 2001. V. 128. N 1-2. P. 281-309.
12. Stuben K. An introduction to algebraic multigrid. In: Multigrid / Ed. U. Trottenberg, C. Oosterlee, A. Schuller. London: Academic Press, 2001. P. 413-532.
13. Yang U.M. Parallel algebraic multigrid methods - high performance preconditioned // Lecture Notes in Computational Science and Engineering. 2006. V. 51. P. 209-236.
14. Волков К.Н. Дискретизация уравнений Навье-Стокса на неструктурированной сетке при помощи метода контрольного объема и разностных схем повышенной разрешающей способности // Журнал вычислительной математики и математической физики. 2008. Т. 48. № 7. С. 1250-1273.
15. AGARD-AR-211: Test cases for inviscid flow field methods // Advisory Reports. 1986.
16. Suero R., Pinto M.A.V., Marchi C.H., Araki L.K., Alves A.C. Analysis of algebraic multigrid parameters for two-dimensional steady-state heat diffusion equations // Applied Mathematical Modelling. 2012. V. 36. N 7. P. 2996-3006.
Волков Константин Николаевич - доктор физико-математических наук, научный сотрудник, Университет ИТМО, Санкт-Петербург, Россия; старший лектор, Университет Кингстона, Лондон, Великобритания, k.volkov@kingston.ac.uk
Konstantin N. Volkov - D.Sc, scientific researcher, ITMO University, Saint Petersburg, Russia;
Senior Lecturer, Kingston University, London, United Kingdom, k.volkov@kingston.ac.uk
Принято к печати 15.05.14 Accepted 15.05.14
