Программная реализация. Программа DifAV (Diffusion on Asymmetric Vortex) реализована в системе MATLAB с помощью инструментария PDEtools

Программа DifAV (Diffusion on Asymmetric Vortex) реализована в системе MATLAB с помощью инструментария PDEtools. Расчет диффузии состоит из двух этапов. На первом этапе рассчитывается уравнение Лапласа (4.18) в заданной геометрии. Геометрические параметры задачи (рис. 4.11), количество точек сетки N и граничное условие (4.21) для уравнения Лапласа для верхней deltaphiup (φ ₁) и нижней deltaphiup (φ ₂) областей движения можно задавать на панели управления под именем vortex (см. рис. 4.12). Решение представляется двумерной функцией φ (x, y), линии уровня которой также показаны на панели vortex (рис. 4.12). Для улучшения cходимости решения уравнения Лапласа (4.18) острые углы внутреннего прямоугольника сглажены окружностью радиуса 0.05 в единицах рисунка (bx=by=1). Частота линий уровня характеризует скорость частицы, согласно системе (4.19). Регулярное движение, таким образом, в верхней широкой части вихря медленное, а в нижней части и по бокам быстрое.

На втором этапе рассчитывается движение частиц в полях (4.19) и (4.22). Параметры расчета можно менять на панели particles (рис. 4.13). На этой панели выдается также количество частиц (particles number), расчетное время (Time) и положение частиц на рисунке линий уровня вихря. Регулярное и тепловое движения можно независимо включать и выключать.

Рис. 4.12

Рис. 4.13

Список литературы к теме 4

Использованная литература: