Q准则的涡结构

军火库

各位大神， 你们好！
     我老师让我做一个Q准则的三维涡量图，由于学的很肤浅，都不知道怎么做，还希望各位大神指点一二。
     希望能具体到步奏，最好是有步奏截图，由于老师催得紧，哪位大神能帮小弟一下的话，那真是感激不尽！
C:\Users\maidou\Desktop

huyanbaopeng

我也是从网上找的，不知道能不能帮助到你
第一步，搜一下tecplot公式编辑规则，网上到处有；
第二步，在tecplot的data->alter->specify equation里面用如下:
{wx}=ddy({w})-ddz({v})
{wy}=ddz({u})-ddx({w})
{wz}=ddx({v})-ddy({u})
这三个是随便一本流体力学书都会提到的涡量计算公式，其中u,v,w分别为三个方向速度，你看了第一步中公式编辑规则就知道,他们的变量符号大概是V5，V6，V7，将V5，V6，V7替代u,v,w带入上三式；
第三步，点击计算，如果网格很多，可能要计算很久。计算完，你就可以向显示速度压力等常规变量一样显示涡量了。


流场计算完成后处理中，有个不变量Q的等值面来描述瞬间涡流结构的方法，就是用三维等势面反映流场结构。Q=0.5×(W×W-S×S)。
其中：W为涡量幅值；S为应变率幅值。
查阅了文献，发现Q=0.5×（ui,j的平方-u i,j×uj,i)=0.5×(W×W-S×S)，对于u，W和S的定义很有一些困惑，
涡量的幅值好理解，就是涡量的大小。但是这个应变率指的是啥，我也不知道。
涡量(wx,wy,wz)是用这个定义
{wx}=ddy({w})-ddz({v})
{wy}=ddz({u})-ddx({w})
{wz}=ddx({v})-ddy({u})
幅值就是这三个的合值
这个涡量是速度梯度的反对称张量对应与fluent中的vorticity，而应变率是速度梯度的对称张量，对应于fluent中的derivative-strain rate。个人认为只是运算时间差别，因此我一般直接导出数据从fluent,包括自定义函数进行Q criterion.比如说，三维fluent导出omega 和Q 只要一分钟左右
假定速度为uvw，用下式显示Q的等值面。
{omega}=0.5*((ddy({w})-ddz({v}))**2+(ddz({u})-ddx({w}))**2+(ddx({v})-ddy({u}))**2)
{S}=ddx({u})**2+ddy({v})**2+ddz({w})**2
{S}={S}+0.5*((ddy({w})+ddz({v}))**2+(ddz({u})+ddx({w}))**2+(ddx({v})+ddy({u}))**2)
{Q}=0.5*({omega}-{S})

周华

UDF中有求导数的方法。

louql

fieldView可以写个Function函数，函数文件我这有