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

In the work, we consider the problem of accelerating the iteration process of the numerical solution of boundary-value problems for partial differential equations (PDE) by the method of collocations and least residuals (CLR). To solve this problem, it is proposed to combine simultaneously three techniques of the iteration process acceleration: the preconditioner, the multigrid algorithm, and the correction of the PDE solution at the intermediate iterations in the Krylov subspace. The influence of all three techniques of the iteration acceleration was investigated both individually for each technique and at their combination. Each of the above techniques is shown to make its contribution to the quantitative figure of iteration process speed-up. The algorithm which employs the Krylov subspaces makes the most significant contribution. The joint simultaneous application of all three techniques for accelerating the iterative solution of specific boundary-value problems enabled a reduction of the CPU time of their solution on computer by a factor of up to 230 in comparison with the case when no acceleration techniques were applied. A two-parameter preconditioner was investigated. It is proposed to find the optimal values of its parameters by the numerical solution of a computationally inexpensive problem of minimizing the condition number of the system of linear algebraic equations. The problem is solved by the CLR method and it is modified by the preconditioner. It is shown that it is sufficient to restrict oneself in the multigrid version of the CLR method only to a simple solution prolongation operation on the multigrid complex to reduce substantially the CPU time of the boundary-value problem solution.Numerous computational examples are presented, which demonstrate the efficiency of the approaches proposed for accelerating the iterative processes of the numerical solution of the boundary-value problems for the two-dimensional Navier–Stokes equations. It is pointed out that the proposed combination of the techniques for accelerating the iteration processes may be also implemented within the framework of other numerical techniques for the solution of PDEs.Рассматривается проблема ускорения итерационного процесса численного решения методом коллокаций и наименьших невязок (КНН) краевых задач для уравнений с частными производными (PDE). Для решения этой проблемы впервые предложено комбинированно приме- нять одновременно три способа ускорения итерационного процесса: предобуславливатель, многосеточный алгоритм и коррекцию решения PDE на промежуточных итерациях в подпространстве Крылова. Исследовано влияние на итерационный процесс всех трех способов его ускорения как по отдельности, так и при их комбинировании. Показано, что каждый из указанных способов вносит свой вклад в количественный показатель ускорения итерационного процесса. При этом наибольший вклад дает применение алгоритма, использующего подпространства Крылова. Комбинированное применение одновременно всех трех способов ускорения итерационного процесса решения конкретных краевых задач позволило уменьшить время их решения на компьютере до 230 раз по сравнению со случаем, когда никакие способы ускорения не применялись. Исследован двухпараметрический предобуславливатель. Предложено находить оптимальные значения его параметров путем численного решения относительно нетрудоемкой задачи минимизации числа обусловленности модифицированной предобуславливателем системы линейных алгебраических уравнений, решаемой в методе КНН. Показано, что в многосеточном варианте метода КНН для существенного уменьшения времени решения краевой задачи достаточно ограничиться только простой операцией продолжения решения на многосеточном комплексе. Приводятся многочисленные примеры расчетов, демонстрирующие эффективность предлагаемых подходов к ускорению итерационных процессов решения краевых задач для двумерных уравнений Навье–Стокса. Указывается, что предложенная комбинация способов ускорения итерационных процессов может быть реализована также в рамках применения других численных методов решения PDE