При моделировании физических процессов и явлений повсеместно встает задача численного решения дифференциальных задач в частных производных. Одним из способов численного решения уравнений математической физики является подход, основанный на конечноразностной аппроксимации. Однако, существенный недостаток данного метода – необходимость использовать ортогональную гексаэдрическую сетку. При решении задач на практике, порой используется сложная геометрическая конфигурация расчетной области (см. например, рисунок 1), и поэтому упомянутый недостаток является достаточно критичным.
Рисунок 1 – Пример сложной поверхности, как части расчетной области