Сравнение сеток для моделирования волновых полей в цилиндрических конструкциях

Любое численное моделирование физических процессов начинается с разбиения расчетной области на конечное множество ячеек — дискретизации. Полученная система ячеек называется вычислительной сеткой. Среди сеточных методов выделяют структурированные и неструктурированные сетки. Структурированные сетки проще в реализации и требуют меньше ресурсов, но плохо приспособлены для областей сложной формы. Цилиндрические и круговые области, такие как подземные трубопроводы, скважины, опоры мостов, стволы горных выработок и тоннели, часто встречаются в прикладных задачах, что делает актуальным поиск оптимального способа наложения структурированной сетки на круглое сечение.

Два подхода: блочно-структурированный и химерный

В блочно-структурированной сетке область разбивается на несколько простых блоков, в каждом из которых строится своя регулярная сетка. Для круга это квадрат в центре и четыре криволинейных «лепестка» между квадратом и окружностью. В трехмерном случае двумерная сетка «вытягивается» вдоль оси. Стыковка блоков требует интерполяции, так как узлы на общих границах не совпадают.

Химерная сетка предполагает перекрывающиеся компоненты: центральную квадратную сетку и кольцевую сетку на периферии. Перекрытие обеспечивает зону для обмена данными, и на каждом временном шаге происходит синхронизация. В трехмерном случае квадратная сетка превращается в параллелепипед, а кольцо — в цилиндрический слой.

Сравнение по точности и скорости

Исследователи провели модельный эксперимент: в центре круга прикладывался точечный импульс, а 360 виртуальных приемников регистрировали волну. В изотропной среде цилиндрическая волна должна иметь идеально круговой фронт. Отклонения — численная анизотропия — возникают из-за структуры сетки. Химерная сетка показала значительно меньший разброс амплитуды по углу, что говорит о ее более высокой точности. Однако время расчета на химерной сетке оказалось на 15–20% больше, чем на блочно-структурированной, при одинаковом шаге.

Николай Хохлов, ведущий научный сотрудник лаборатории прикладной вычислительной геофизики МФТИ, объяснил: «Мы хотели дать практикующему вычислителю простую инструкцию: если нужна точность — берите химерную сетку, если скорость — блочно-структурированную». Авторы подчеркивают, что использовали один и тот же алгоритм и одинаковые физические условия для получения надежных количественных данных.

Работа поддержана Российским научным фондом (грант №25-71-10027). В будущем ученые планируют исследовать другие типы сеток, в том числе неструктурированные и адаптивные, а также проверить оба подхода в сочетании с другими численными методами. Кроме того, аналогичный анализ необходим для конических и сферических геометрий, где проблема полюсной особенности стоит еще острее.