Numerical simulation of flow around a system of bodies in stream function–vorticity variables

A numerical algorithm for simulating non-stationary viscous 2D flows based on solving Navier–Stokes equations in stream function–vorticity variables is developed. The algorithm allows to model the forced flows as well as the complex convective flows around a system of bodies in channel. Problem definition is amplified with nonlocal integral equalities based on Pearson conditions in order to obtain boundary values of stream function.