Added in edit: I just went back and reread your post. In the "specific case of linear systems", after you have found the eigenvalues and eigenvectors, draw straight lines along the direction of the eigenvectors, including their directions as t increases. Draw the other curves close to those lines following their direction and "asymptotic" to those lines.
If your eigenvalues are complex, then, of course, you will have circles or spirals as phase lines.
I assume you have equations of the form dx/dt= f(x,y) and dy/dt= g(x,y). At each point, (x,y), then dy/dx= g(x,y)/f(x,y). Choose a number of points in the plane, and calculate dy/dx for each of them. That gives you the slope of the line through that point that is tangent to the solution curve through that point. Draw a short line through the point with that slope (you might find it easier to use f(x,y) and g(x,y) separately to fine the "run" and "rise"). Oh, and be sure to put a little "arrowhead" on the line to show the direction of "flow" as t increases- use the signs of f(x,y) and g(x,y) to get that.
The more general case, dx/dt= f(x,y,t) and dy/dt= g(x,y,t) is harder- you have a different phase plane for every value of t.