Hello, i’m a new user of fenics.

I’d like to get the streamfunction psi of a given 2d thermal heat flux F=k.grad(T) so that the isolines of psi are parallel to F (or orthogonal to isolines of T).

- First i solve the thermal equation to get the vector flux using a P2 interpolation
- Then, streamfunction is a well-known problem in CFD and psi satisfies a vorticity equation defined here

Unfortunatly the streamlines are not parallel to F.

Did i miss something ?

thanks for helping.