Контроль точности при численном интегрировании жестких систем

Ранее для численного решения жестких систем обыкновенных дифференциальных уравнений было предложено а) использовать в качестве аргумента длину дуги интегральной кривой и б) выбирать оптимальный шаг интегрирования по кривизне интегральной кривой. В данной работе построена тестовая задача, в которой точное решение представлялось через элементарные функции как аргумента времени t, так и аргумента дуги l. Это позволило провести количественное сравнение различных разностных схем. Показано, что при расчетах с оптимальным шагом удается использовать даже явные схемы Рунге-Кутты. При этом схема первого порядка давала невысокую точность, но очень высокую надежность даже при огромной жесткости. С повышением порядка точности надежность схем ухудшалась. Предложена смешанная стратегия. На первом этапе по надежной схеме первого порядка строится оптимальная сетка, адаптированная к решению. На втором этапе эта сетка сгущается по правилу дробления квазиравномерных сеток, а расчет выполняется по схеме четвертого порядка точности. Смешанная стратегия дает одновременно хорошую надежность и высокую точность расчета.

Previously, for numerical solution of stiff systems of ordinary differential equations, it was proposed to a) use length of the integral curve as the argument and b) choose optimal integration step cording to curvature of the integral curve. In this work, we construct a test problem in which the exact solution is expressed via elementary functions for both time and arc length arguments. This permitted quantitative comparison of various differential schemes. We show that even explicit Runge-Kutta methods are applicable in calculations with optimal step. The first order scheme provides low accuracy but very high reliability even at enormous stiffness. As order of accuracy increases, reliability of the schemes decreases. We propose mixed computation strategy. At the first stage, an optimal mesh adapted to solution is built via the first order scheme. At the second stage, this mesh is thickened according to the rule of quasi-uniform meshes splitting and the calculation is done via the scheme with the fourth order of accuracy. The mixed strategy allowed to achieve both good reliability and high accuracy of calculation.