Метод сеток
англ. Grid method
Краевая задача для урчп решается на сетке: производные заменяются разностями, уравнение — системой алгебраических уравнений для узловых значений; главный инструмент вычислительной матфизики.
Непрерывную задачу компьютер переваривает лишь после перевода на язык таблицы чисел. Метод сеток — общий приём решения краевых задач для уравнений в частных производных: область покрывается сеткой узлов, производные заменяются разностными отношениями по соседним узлам, и вместо дифференциального уравнения возникает система алгебраических уравнений для значений в узлах. Набор узлов, через которые записана каждая разность, называют шаблоном, а полученную конструкцию — разностной схемой. Подход универсален: остывающий стержень, бегущая волна, стационарное поле — везде уравнение превращается в сеточную задачу, которую решают методами линейной алгебры. Особенно естественно он ложится на постановки, где данные заданы на границе целиком: типы условий — Дирихле, Нейман и третий род — собраны в статье краевая задача.
Классический пример — пятиточечный шаблон для уравнения Лапласа: значение в узле выражается через четырёх соседей, и на сетке задача превращается в систему из уравнения. Для остывающего стержня длины с шагом получается одиннадцать узлов и девять внутренних неизвестных; явная схема по времени требует , то есть около двухсот шагов на одну секунду процесса. Типичные грабли: нарушенная устойчивость — решение взрывается осцилляциями; неточная передача граничных условий — ломается симметрия, обещанная физикой; слишком грубая сетка — резкие фронты превращаются в лесенку. Стационарные задачи, где время отсутствует, — территория уравнения Лапласа: его решения не имеют внутренних экстремумов, и сеточный максимум обязан сидеть на границе.
Частые вопросы
Чем метод сеток отличается от метода конечных элементов?
Метод сеток делит область на ячейки правильной формы и заменяет производные разностями в узлах — идеально для прямоугольных областей. Метод конечных элементов режет область на треугольники и ищет решение комбинацией базисных функций на элементах. Плата за простоту разностей — трудности на криволинейных границах; зато элементы подстраиваются под геометрию естественно. Поэтому для сложных областей выбирают конечные элементы, для регулярных расчётов — разности: меньше кода и выше скорость.
Как выбирать шаг по времени?
Ограничения два: устойчивость и точность. Параболическим уравнениям явные схемы диктуют квадратично малый шаг: — при и это . Гиперболическим работает условие Куранта : возмущение не должно пробегать за шаг дальше соседней ячейки. Неявные схемы снимают границу устойчивости ценой решения системы уравнений на каждом шаге. По точности шаг выбирают, удваивая разрешение, пока ответ не перестанет меняться.