Итерационный метод расчета течений вязкопластической среды Бингама

Рассматривается разностная схема для расчета плоских течений вязкопластической среды Бингама. В качестве математической модели среды используется вариационное неравенство Дюво–Лионса. Обе компоненты скорости аппроксимируются в узлах основной сетки, давление и все компоненты тензоров скоростей деформации и напряжений — в узлах полусмещенной сетки. Показано, что итерационный метод типа Узавы, используемый для нахождения решения вариационного неравенства, требует специальной адаптации в случае сеточной задачи. В качестве модельного примера рассматривается численное решение задачи о течении вязкопластической среды в каверне. Полученные результаты сравниваются с известными. Статья подготовлена к печати во время пребывания Е.А. Муравлëвой в Max Planck Institute for Mathematics in the Sciences, 04103, Leipzig, Germany.