Представлены модификация для трехмерного случая и программная реализация алгоритма расчета неустойчивости Рэлея - Тейлора в высоковязкой несжимаемой ньютоновской жидкости. Алгоритм основан на использовании аналитического выражения функции Грина соответствующей краевой задачи, что значительно ускоряет расчет, особенно при использовании параллельных вычислений. Разработанное программное обеспечение может работать как на ресурсах центральных процессоров (CPU) с использованием инструкций SSE и протокола MPI, так и на графических ускорителях (GPU) NVIDIA в случае доступности узлов гибридно-го кластера. Взаимодействие с GPU реализовано с помощью технологий CUDA и MPI. Работа ориентирована на моделирование геологического процесса - соляного диапиризма как частного случая неустойчивости Рэлея - Тейлора, выражающегося во всплывании твердой каменной соли через более плотные перекрывающие осадки в геологическом масштабе времени
