Научный семинар «Экономичные методы численного интегрирования задачи Коши для систем ОДУ с контрастными структурами и сингулярностями»
28 апреля с 15:00 до 16:20 по московскому времени
Многие проблемы физики и техники приводят к жестким начальным и начально-краевым задачам (сводящиеся к начальным задачам огромной размерности). Примерами являются кинетика реакций, кумуляция, взрыв, пробой в плазме и полупроводниках, выбор оптимальных режимов зажигания термоядерных мишеней и т.д. Другой характерной особенностью таких задач является плохая обусловленность: стремительное нарастание ошибок округления сильно ограничивает точность. В некоторых из перечисленных задач решение имеет сингулярность, то есть обращается в бесконечность за конечное время. Наличие сингулярностей представляет серьезную трудность для численного расчета.
С 1970-х годов для задач Коши предлагались методы автоматического выбора шага, адаптированного к решению. Однако они не гарантировали получения заданной пользователем точности и нередко приводили к грубо ошибочным результатам (выдаваемая оценка погрешности могла отличаться от фактической точности на 9-10 порядков).
В данной работе предложен принципиально новый метод автоматического выбора шага по кривизне интегральной кривой в аргументе длина дуги, названный геометрически-адаптивным. Он прост, имеет наглядный геометрический смысл и значительно снижает трудоемкость вычислений. Доказана оптимальность этого алгоритма в смысле метрики Хаусдорфа. Разработана процедура расчета, позволяющая одновременно с решением он вычисляет асимптотически точное значение погрешности. Предложенный алгоритм выбора шага обеспечивает хорошую надежность расчетов даже при использовании явных схем невысоких порядков точности. Построены экономичные формулы для кривизны интегральной кривой, позволяющие применять этот алгоритм с явными схемами Рунге-Кутты.
Для задачи кинетики реакций традиционно использовались трудоемкие неявные схемы. В данной работе предложена специализированная явная схема, имеющая достаточную надежность и отличающаяся очень малой трудоемкостью. Показано, что эта схема превосходит известные схемы по точности и надежности.
Для задач Коши с сингулярностями решения предложен новый способ численного обнаружения и исследования ближайшей особенности, основанный на использовании длины дуги интегральной кривой в качестве аргумента. Это кардинально повышает надежность и точность исследования. Метод позволяет не только определить тип особенности, но и вычислить ее порядок и момент достижения с гарантированной точностью. Конкретные формулы метода получены для полюсов степенного и логарифмического типов и для произведения степени на логарифм.
Для задач Коши со множественными полюсами целого порядка предложен метод инверсной функции. Он позволяет продолжать решение за полюс, вычисляя само решение и положение полюса с высокой точностью вплоть до ошибок компьютерного округления.
Докладчик
Белов Александр Александрович, к.ф.-м.н., доцент каф. Прикладной информатики и теории вероятностей.