On the far-field stream function condition for two-dimensional incompressible flows

The present demonstration of the usefulness of the integral series expansion of the stream function as a far-field computational boundary condition shows the method to require only a 10-percent/time-step increase in computational effort over alternative boundary conditions, in the case of implementation of unsteady problems using a direct elliptic solver. So long as the vorticity was encompassed within the computational domain, the method proved sufficiently accurate to yield virtually identical results for two widely different domains. While the integral-series condition yielded the best results for periodic flow, the Neumann condition gave comparable accuracy with less computation time for the steady-flow case despite its inability to treat periodic flow with vortex shedding.