Numerical solution of transonic stream function equation

The stream function equation, in conservation form, looks similar to the full potential equation and existing methods (e.g. artificial compressibility) can be readily applied. Rotational flows can be calculated once the vorticity (due to shocks or nonuniformity) is evaluated. There are, however, two main difficulties: First, the density is not uniquely determined in terms of the flux (there are two solutions; the subsonic and the supersonic branch with a square root singularity at the sonic point). Methods to overcome this difficulty are studied and results are presented with some remarks on inviscid separation and closed stream lines. Second, the need of two stream functions for three dimensional calculations is briefly discussed.