Matlab实现Q-criterion涡结构可视化

这篇博文主要是因为有个小伙伴表示用Paraview读取CSV格式数据可视化不好用，边界会有阴影，希望换一种方式，我便顺便写一下Q-criterion的实现。

首先，得到速度的数据矩阵以后，可用gradient分别对速度的三个分量计算得到9个分量，这里0.01是因为我的dx是0.01。

[dudx, dudy, dudz] = gradient(ux,0.01);
[dvdx, dvdy, dvdz] = gradient(uv,0.01);
[dwdx, dwdy, dwdz] = gradient(uw,0.01);

得到▽v之后，网上随便查一下Q值的计算，如下：

q=-0.5*(dudx.^2+dvdy.^2+dwdz.^2)-dudy.*dvdx-dudz.*dwdx-dvdz.*dwdy;

如果有需要的话，可以看一下最大值最小值范围。

max(max(max(q)))
min(min(min(q)))

比如取iso_q=0.1，即0.1为等值面。

figure()
p=patch(isosurface(xxx,yyy,zzz,q,iso_q, Vel_magnitude));
set(p,'FaceColor','interp','EdgeColor','none');
axis off
colormap jet

这里的上色为速度大小，值则存储在Vel_magnitude。

后续如何调整色调，光亮不再多言。