Irrotational Flow with Irregular Boundary Shapes

An algorithm for the numerical solution of two-dimensional irrotational flows based on an inverse transformation is presented as an efficient means of solving flows with irregular boundary shapes. Using the intersections of a set of streamlines and equipotential lines as computational nodal points, the algorithm determines the coordinates of these intersections. The streamlines and equipotential lines can be arbitrarily selected to concentrate the nodal points in areas of high velocity flow so as to increase the accuracy in these areas. The algorithm is particularly well suited to modern high speed computers, and FORTRAN statements are given. Since the program operates on many different boundary shapes with no internal alterations, it is suitable for use by persons with limited background in numerical procedures.