How can I get velocity values for stream3 function
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hello,
I have a 3D flow domain in which I have created matrix format of x, y , z coordinates as well as u , v , w for my velocity vector. I would like to do a particle tracing for the particle release from a cetrain coodrnates. Then I would like to have x , y, z coordinates and u, v, w values of the particle along it's trajectory (my final aim is to calculate the particle travel time within the domain released from a specfic point).
I used the following codes to get the coordinates of the particle along it's path (x,y,z,u,v,w are known):
[sx,sy,sz] = meshgrid(2.6,2.6,-5);
XYZ_Data=stream3(x,y,z,u,v,w,sx,sy,sz);
However, I don't know how can I get the velocity vector of the corresponding coordinates of XYZ_Data! It would be really appreciated if anybody could help me in this issue.
Bests,
0 commentaires
Réponse acceptée
darova
le 6 Avr 2020
A short example
clc,clear
cla
load wind
[sx sy sz] = meshgrid(80,20:10:50,0:5:15);
h = streamline(stream3(x,y,z,u,v,w,sx,sy,sz));
xx = get(h,'xdata');
yy = get(h,'ydata');
zz = get(h,'zdata');
ii = 10;
x1 = xx{ii};
y1 = yy{ii};
z1 = zz{ii};
u1 = interp3(x,y,z,u,x1,y1,z1);
v1 = interp3(x,y,z,v,x1,y1,z1);
w1 = interp3(x,y,z,w,x1,y1,z1);
hold on
jj = 1:30:length(x1);
quiver3(x1(jj),y1(jj),z1(jj),u1(jj),v1(jj),w1(jj),'r')
hold off
2 commentaires
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Vector Fields dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!