If anyone can do it please help me in this regard

Stream function is given by:

With the help of MATLAB, plot the streamlines for the following parameters:

1. V = 1 m/s
2. a= 0.1 m
3. Ψ = -0.16 to 0.16 (with a step of 0.02)
4. 0o < θ < 180 o(with a step of 1o)

Set the limits of plot to -0.4 < x < 0.4 and 0 < y < 0.4

1. Plot velocity vectors at a few points and comment on their orientation with respect to streamlines.

[madhan ravi]

ok where is the code for the equation?

[Walter Roberson]

What are you asking of us? it sounds like you are asking us to post the answer to a homework question .

[Muhammad Adeel Ahsan Awan]

i dont know even how to start code for this, actually i am a beginner

[Muhammad Adeel Ahsan Awan]

its just for the sake of understanding how to plot streamlines , actually there is anotherway by manually , not an homework question

[Walter Roberson]

read the Getting Started section of the documentation . Take the free online course .

[Muhammad Adeel Ahsan Awan]

Whats wrong in this code it is not plotting and showing the white blank graph..plz help in this regard

[Mark Sherstan]

Please format your code instead of posting screenshots. Information regarding how to ask good questions and properly posting them can be found here.

[Muhammad Adeel Ahsan Awan]

okay sir

[Muhammad Adeel Ahsan Awan]

a=0.1
v=1

s=-0.16
t=0
   while s<0.16 && t<180
   r=(-s)+((s*s)+(4*a*a*v*v*sin(t)*sin(t))^(1/2))/(2*v*sin(t))
   x=r*cos(t);
   y=r*sin(t);
   s=s+0.02;
   t=t+1;
   plot(x,y,'r-')
   xlim([-0.4 0.4])
   ylim([0 0.4])
   end

Answer 1

[Mark Sherstan]

You were plotting a single point instead of an array of information. Adding the ii term stores the result for each calculation and the plots it. The following code should now work for you:

a=0.1
v=1
s=-0.16
t=0
ii = 1;

while s<0.16 && t<180
   r=(-s)+((s*s)+(4*a*a*v*v*sin(t)*sin(t))^(1/2))/(2*v*sin(t))
   x(ii)=r*cos(t);
   y(ii)=r*sin(t);
   s=s+0.02;
   t=t+1;
   
   ii = ii + 1;
end


plot(x,y,'r-')
xlim([-0.4 0.4])
ylim([0 0.4])

This is also how you should post code so that people can help you (not just copying and pasting the text). If this solves your original question, please accept the answer so that we can close it, otherwise add further comments below.

[madhan ravi]

Muhammad Adeel Ahsan Awan's answer moved here for consistency:

actually we need to convert degree to radians

[Walter Roberson]

sind and cosd

[madhan ravi]

deg2rad()

[Muhammad Adeel Ahsan Awan]

yes i have done it thanks to all of you