Compute 3D stream surfaces?
 

It is stated that 3D flow can be represented by two scalar stream surfaces. From these, the velocity field can be computed, and stream lines are the intersection of the surfaces. For instance in "2D flow", one surface might be the stream function and the other a set of parallel planes, with intersections projected onto one of the planes showing the stream lines via contour plots.

But for general 3D flow, how might one compute the two stream surfaces on a mesh, either from the flow, or from the governing equations?

Hi,
this is an interesting question for me, just have some ideas to work on.
Starting by the definition you have a potential vector such that

v = curl Psi

that involves 3 scalar functions. Defining a sarface that has zero normal velocity component would lead to satisfy the constraint

n.curl Psi = 0

Now I am thinking to the analogy with a vortical tube, where the surface of the tube is constrained to have zero normal vorticity (see the book of Chorin and Marsden).
Could we construct the stream surface in the same way?