Аннотация:
Рассматривается проблема ускорения итерационного процесса численного решения методом коллокаций и наименьших невязок (КНН) краевых задач для уравнений с частными производными (PDE). Для решения этой проблемы впервые предложено комбинированно применять одновременно три способа ускорения итерационного процесса: предобуславливатель, многосеточный алгоритм и коррекцию решения PDE на промежуточных итерациях в подпространстве Крылова. Исследовано влияние на итерационный процесс всех трех способов его ускорения как по отдельности, так и при их комбинировании. Показано, что каждый из указанных способов вносит свой вклад в количественный показатель ускорения итерационного процесса. При этом наибольший вклад дает применение алгоритма, использующего подпространства Крылова. Комбинированное применение одновременно всех трех способов ускорения итерационного процесса решения конкретных краевых задач позволило уменьшить время их решения на компьютере до 230 раз по сравнению со случаем, когда никакие способы ускорения не применялись. Исследован двухпараметрический предобуславливатель. Предложено находить оптимальные значения его параметров путем численного решения относительно нетрудоемкой задачи минимизации числа обусловленности модифицированной предобуславливателем системы линейных алгебраических уравнений, решаемой в методе КНН. Показано, что в многосеточном варианте метода КНН для существенного уменьшения времени решения краевой задачи достаточно ограничиться только простой операцией продолжения решения на многосеточном комплексе. Приводятся многочисленные примеры расчетов, демонстрирующие эффективность предлагаемых подходов к ускорению итерационных процессов решения краевых задач для двумерных уравнений Навье–Стокса. Указывается, что предложенная комбинация способов ускорения итерационных процессов может быть реализована также в рамках применения других численных методов решения PDE.
Ключевые слова:предобуславливание, подпространства Крылова, многосеточные алгоритмы, итерации Гаусса–Зейделя, уравнения Навье–Стокса, метод коллокаций и наименьших невязок.