Исследование точности и применимости разностной схемы для решения задачи диффузии-конвекции при больших сеточных числах Пекле

Работа посвящена исследованию точности разностной схемы, используемой при решении задачи диффузии-конвекции в случае больших сеточных чисел Пекле. Численное решение задачи транспорта взвешенных частиц предлагается проводить на основе усовершенствованной схемы «кабаре». Ее разностный оператор представляет собой линейную комбинацию операторов разностных схем «крест» и «кабаре», при этом модифицированная схема получается из схем с оптимальными весовыми коэффициентами. При определенных значениях весовых коэффициентов объединение приводит к взаимной компенсации ошибок аппроксимации, а результирующая схема приобретает лучшие, чем исходные схемы, свойства. Кроме того, она включает функцию заполненности ячеек, что позволяет естественным образом моделировать задачи в областях со сложной геометрией. Вычислительные эксперименты осуществлены на решении задачи транспорта взвеси, возникающей, например, при распространении шлейфов взвеси в водной среде и изменении рельефа дна в связи с выпадением взвешенных частиц грунта в осадок при выгрузке грунта в водоем (дампинге). Приведены результаты моделирования транспорта взвешенных частиц при различных значениях сеточного числа Пекле. Реализация алгоритма осуществлялось при помощи программно-аппаратной архитектуры параллельных вычислений: на центральном процессоре (Central Processing Unit - CPU) и на графическом ускорителе (Graphics Processing Unit - GPU). Решение прикладной задачи показало свою эффективность на CPU при расчетных сетках небольшого размера, а при необходимости измельчения шагов по пространству предпочтительным оказалось решение на GPU. Выяснено, что при использовании модифицированной схемы «кабаре» увеличение скорости движения водного потока не приводит к потере точности решения за счет диссипативных источников и сопровождается незначительном ростом вычислительных трудозатрат.