How to plot streamlines in matlab ???

Question

ANTONIO LEUZZI 2019-2-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/446892-how-to-plot-streamlines-in-matlab

回答： ANTONIO LEUZZI 2019-3-28

Cartel1.xlsx

Hi everybody, I m new. I ve this question like title. I have to plot streamlines with circular cylinder, from PIV data of velocity. For each data I ve one value of X Y U V (2D model), a matrix with n rows x 4 columns. From this matrix, that represents the average values of velocity, I ve just elaborated a quiver plot with the cylinder but I need of your help to understand how can I use streamline function in matlab or if is necessary to write the strem function like fluid mechanics theory. In annex the matrix

Thanks

6 个评论
显示 4更早的评论隐藏 4更早的评论

Bob Thompson 2019-2-25

How are you identifying which points connect to form a streamline?

ANTONIO LEUZZI 2019-2-25

this is my problem Bob; I tried in different ways to buil the grid first and after to plot the U and V values on meshgrid X and Y, but the result is always a square grid. I Know how to not consider the cylinder (NaN values for U and V in the coordinates X and Y of its position respect the grid) but at this point I m not sure of my grid. Have you confidence with stream function?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Bob Thompson 2019-2-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/446892-how-to-plot-streamlines-in-matlab#answer_362696

Try using the quiver function. It's not perfect streamlines, but it should give you an idea of what the velocity field looks like, which should essentially achieve the same affect. This should still plot in a 'square' because I suspect that is what is available for your (x,y) data, so that will be the bounds of your 2D plot area. Based on the image you attached you should be getting a grid that is roughly rectangular anyway, as it is a cross section of your flow perpendicular to the cylinder you are flowing around.

data = xlsread('Cartel1.xlsx');
quiver(data(:,1),data(:,2),data(:,3),data(:,4));

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

ANTONIO LEUZZI 2019-2-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/446892-how-to-plot-streamlines-in-matlab#answer_362752

Sorry Bob but I think that you don t Know or you aren t understanding my problem!!! If you read with attention my previously messages I ve write that I ve already plotted quiver data!!!!!!!!! My problem is HOW TO PLOT STREAMLINE WITH MY VELOCITY DATA??? HOW TO CREATE A GRID OF POINT ??? Are you understandong now?

12 个评论
显示 10更早的评论隐藏 10更早的评论

Bob Thompson 2019-2-26

I'm going to go through this step by step to make sure I am understanding, if I am not, feel free to correct me.

I know that streamline dericìved by stream function and that stream function has got an equation for two both components of velocity (U and V), but my trouble is how to write it in a code.

You have read the documentation and you understand that the streamline function calculates the streamline values for U and V directions using the stream function, but you are not sure how to write a) the stream function in code, b) the streamline command in code?

The streamline function in Matlab actually just generates a line based on the data field that you feed it. I have not looked in greater detail, but I believe this is more a process of generating a series of points with interpolation than it is spceifically using a fluid dynamics stream function.

To respond to problem (a), I do not know the stream function equations off the top of my head, but I also don't know that they are necessary in this situation. If you are looking to write the code to generate streamlines without using the streamline() command then the basic pseudo code I came up with is this:

data = xlsread('Cartel1.xlsx');
ns = 15; % Number of streamlines
streams = zeros(max(data(:,1))*10,2,ns); % Preallocate streamline array size
                                         % (xlocation,ylocation,streamline number)
streams(1,1,:) = 0;
streams(1,2,:) = [min(data(:,2)):(max(data(:,2))-min(data(:,2)))/(ns-1):max(data(:,2))];
c = 1;
for i = 0.1:0.1:max(data(:,1)); % values from x = 0 to max x value
    c = c + 1; % Array index counter
    streams(c,1,:) = i;  % Set new x values
    for j = 1:ns;        % Loop through each stream y value
        streams(c,2,j) = streams(c-1,2,j)+interp2(data(:,1),data(:,2),data(:,4),i,steams(c-1,2,j));
    end
end

To respond to problem (b), this is actually quite easy, because you already have the x,y,u,v data necessary to create the streamlines, you just need to decide what the starting points are.

streamline(data(:,1),data(:,2),data(:,3),data(:,4),[zeros(1,100)],[min(data(:,2)):(max(data(:,2))-min(data(:,2)))/99:max(data(:,2))]);
% Should plot 100 streamlines starting on left edge

How can I eventually write the stream function by myself , without streamline of Mat?

See above for beginnings of a code. I do wonder why you feel a need to write this yourself though. Is there a specific problem with the built in streamline command?

Or if I follow the stream2D of Mat hoe is necessary meshgrid of data??

I'm assuming you're getting this idea from the examples found in the streamline documentation. It is not necessary for you to generate a mesh grid because you already have the x,y data in your excel document. The examples shown use meshgrid() because they are generating a generic sample set of data so other poeple can recreate their example.

Because if I meshgrid my data onbtain a square grid

You should be getting a square gird because you are looking at a 2D plot, and those are almost always presented as square grids.

Jan 2019-3-7

@ANTONIO LEUZZI: Uppercase means shouting in an internet forum. With bold type also and a bunch of exclamation marks your message looks excited. Calm down.

Post the current code and a copy of the complete error message.

ANTONIO LEUZZI 2019-3-7

There wasn t absolutely my intention to shout but I ve highlighted ineffectiveness for resolving thsi problem despite Bob's help!!!!! I m sorry for this misunderstood and thanks for your suggest

请先登录，再进行评论。

Answer 3

ANTONIO LEUZZI 2019-3-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/446892-how-to-plot-streamlines-in-matlab#answer_367882

Hei guys. Unfortunately I ve not complete this part of code yet. I would still need of your help. I will try to be clearest possible. Than, I ve PIV measurements that consist of X,Y,U,V data, or rather each value has got an X displacement and a Y dispalcement at wich corrispond U velocity component (of dispalcement) and V velocity component (of displacement) respectively. So I ve a matrix of this value. My ask is hoe to plot streamlines with the presence of a cylinder? I ve attached the code with the data. In this code you can find first part in wich the data were arranged and in wich were erased spikes values. After that there s a second part in wich I ve computed main velocity, for U and V, plotted velocity profiles, for U and V and computed Reynolds Shear Stress. I attached even an Excel file with the matrix of data. I know the stream function in matlab bu I think that in my case would need another way to plot the stremlines.

This image is the velocity displacement with a quiver function. Streamlines follow easily the row wise. As you can see from the Excel file, Y length is 166.0957 that is the legth of the regitred PIV image; also for X length. But the process of the data in PIV were made for each Y value constant corrispond all X values and so U values and V values. Than streamlines follow X data posted on the same ordinate.

I don t try to extract them. Please help me. ANd thanks for your help

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

How to plot streamlines in matlab ???

6 个评论
显示 4更早的评论隐藏 4更早的评论

回答（3 个）

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

12 个评论
显示 10更早的评论隐藏 10更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

How to plot streamlines in matlab ???

6 个评论 显示 4更早的评论隐藏 4更早的评论

回答（3 个）

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

12 个评论 显示 10更早的评论隐藏 10更早的评论

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

Community Treasure Hunt

6 个评论
显示 4更早的评论隐藏 4更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

12 个评论
显示 10更早的评论隐藏 10更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论