MATLAB Answers

Plotting a given velocity field (vortex flow)

조회 수: 123(최근 30일)
Mammadbaghir Baghirzade
Mammadbaghir Baghirzade 2020년 4월 25일
편집: Mammadbaghir Baghirzade 2020년 11월 10일
Hi all,
I am trying to plot following u_x and u_y functions, which are the x and y components of a velocity field.
The results should look like a vortex, as can be seen in the attached picture which I took from "wolframalpha.com" when I tried to compare my results.
However, instead such a simle 2D vortex I get complex velocity field which is also attached here.
I suspect it is due to wrong initialiazing the domain or wrong usage domain values in the "for loop", which is shown in the end of this message.
Could you please suggest me any possible solution.
I appreciate your time.
  댓글 수: 2
Mammadbaghir Baghirzade
Mammadbaghir Baghirzade 2020년 4월 25일
Thank you for your reply.
I did plot these functions, however, I got complex result as in the attached.
But instead I need to get something like following. And I can not get such simple vortex flow.

댓글을 달려면 로그인하십시오.

채택된 답변

Ameer Hamza
Ameer Hamza 2020년 4월 25일
편집: Ameer Hamza 2020년 4월 25일
try this
x = 0:0.1:3;
y = 0:0.1:3;
[X,Y] = meshgrid(x,y);
Psi = (1/pi)*sin(pi*X).^2.*sin(pi*Y).^2;
u_x = 2*(sin(pi*X).^2).*(sin(pi*Y)).*(cos(pi*Y));
u_y = -2*(sin(pi*Y).^2).*(sin(pi*X)).*(cos(pi*X));
figure;
contourf(X,Y,Psi)
title('velocity field intensity')
figure;
quiver(X,Y,u_x,u_y)
title('velocity field direction')
  댓글 수: 4
Ameer Hamza
Ameer Hamza 2020년 4월 26일
I am glad to be of help.

댓글을 달려면 로그인하십시오.

추가 답변(0개)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by