Velocityfield and stream function

  • Thread starter Thread starter MaxManus
  • Start date Start date
  • Tags Tags
    Function Stream
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
1 reply · 4K views
MaxManus
Messages
268
Reaction score
1

Homework Statement


A velocity field is given by v = ui + vj
where:
u = cos(x)sin(y)
v = -sin(x)cos(y)
x = y = [-pi/2,pi/2]


the stream function psi is:

psi = cos(x)cos(y)

Matlab function:
--------------------------------------
function[X,Y,psi] = streamfun(n)

if nargin < 1;
n=20;
end

x=linspace(-0.5*pi,0.5*pi,n);
[X,Y] = meshgrid(x,x);
psi=cos(X).*cos(Y);

------------------------------------------

Homework Equations




calculate the velocity numerically with [tex]\nabla[/tex]psi


The Attempt at a Solution


I am starting with the assumption, the velocity is [tex]\nabla[/tex]psi, is this correct?

--------------------------------------

[x,y,psi] = streamfun(20)
[u,v] = gratient(20)

quiver(x,y,u,v)

and I get numfield.jpg, see attachment, but this is not the same as my velocityfield.jpg

which I got from:

-----------------------------------------------
x = linspace(-0.5*pi,0.5*pi,n);
[x,y] = meshgrid(x,x);

u = cos(x).*sin(y);
v = -sin(x).*cos(y);

quiver(x,y,u,v)
----------------------------------------------------------
Can someone please help me and tell me what I am doing wronge.

Edit:
nabla*psi = -sin(x)cos(y)i - cos(x)*sin(y)j
 

Attachments

  • numfield.jpg
    numfield.jpg
    37.8 KB · Views: 510
  • velocityfield.jpg
    velocityfield.jpg
    37.6 KB · Views: 482
Last edited:
Physics news on Phys.org
One more question:

A velocity field is given by v = ui + vj
where:
u = cos(x)sin(y)
v = -sin(x)cos(y)

Draw stream vectors along the x and y axes.

What is a stream vector?
Is it just like arrows with magnitude and direction v?