Q准则涡识别方法


Q准则(Q Criterion)在涡识别中计算效率高,效果也不错,是一种常用的涡提取/识别方法。

了解Q准则需要从速度梯度张量入手,而速度梯度张量可以分解成两部分

\[\frac{\partial u_i}{\partial x_j} = 0.5\left[ \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right] + 0.5\left[ \frac{\partial u_i}{\partial x_j} - \frac{\partial u_j}{\partial x_i}\right] \]

其中,对称部分记作\(S\),通常被称为应变速率张量;反对称部分记作\(\Omega\),通常被称为旋转速率或涡量张量。

\[S = 0.5\left[ \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right] \\ \Omega = 0.5\left[ \frac{\partial u_i}{\partial x_j} - \frac{\partial u_j}{\partial x_i}\right] \]

反观粘性应力张量的定义

\[\bm{\tau} = \mu\left[ \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right] \]

只是应变速率张量的具体函数。

考虑到这一点,Q值的定义为速度梯度张量的第二不变量

\[Q = 0.5\left[ ||\Omega||^2_F - ||S||^2_F\right] \]

可以看出,Q的正值表示流场中涡量占主导的区域,负值表示应变速率或粘性应力占主导的区域。


参考资料:Q-criterion for Vortex Visualization


免责声明!

本站转载的文章为个人学习借鉴使用,本站对版权不负任何法律责任。如果侵犯了您的隐私权益,请联系本站邮箱yoyou2525@163.com删除。



 
粤ICP备18138465号  © 2018-2025 CODEPRJ.COM