Finding streamlines in blood flow

[neochan]

Homework Statement

Hi!

I'm creating my master thesis in computer science. The title is 'simulator of magnetic resonance angiography'. My problem is that I have to create some kind of overlay on my image with streamlines of flowing blood. I have to arrays with velocity vectors in two directions up/down and left/right - each pixel has its two values in these two directions. How can I simply create streamlines using this data?

Homework Equations

I'm using Eulers approach - connecting vectors point by point

The Attempt at a Solution

I've tried to use Eulers approach, but at some point it didn't work (see picture)

as you can see, this line isn't going right way. it should go into right part of the vessel, but it's not. Colours are representation of value of velocity (green means fast, red - slow). Is there any simple solution for my problem?

I'm creating it in c++ and qt if this information can help you :)

 

[marcusl]

If using Matlab is an option to you, then you can use its streamline plot function. You can easily overlay the lines onto your image. You can also overlay arrows representing the vector values.

 

[Chestermiller]

You can calculate the stream function, and plot the lines of constant stream function. This will give the stream lines. The main direction of the velocity in your problem sees to be in the upward (y) direction. The partial derivative of the stream function with respect to x is equal to the velocity in the y-direction. You can integrate with respect to x to get the stream function values along each level of y. The stream function should be constant with respect to x along each boundary of the channel. Take the value at the left boundary equal to zero. Hopefully, the values of the stream function you calculate at the right boundary will all be approximately equal. This will be the case if the flow is incompressible. All this hinges, of course, if the flow is two dimensional (channel flow). If those are tubes, you have to proceed differently.