基于流线分析涡旋识别.pdf

基于流线分析的涡旋识别1 邹文楠 赵勇 董荭 南昌大学工程力学研究所, 江西南昌(330029) email: zouwn@ 摘 要：涡旋识别是流场分析中的核心课题之一。虽然目前还没有明确的涡旋定义，人们还 是发展了多种基于特征量(feature-based)的涡旋识别技术。本文以流线为研究对象，分析 流线特征并流线之间的缠绕关系进行涡旋的识别，确立速率、曲率、挠率作为流线特征识别 的参数系统，并推导了相邻流线之间缠绕关系的特征判据来定义涡旋(核)区域，建立了新的 涡旋识别的三参数系统，最后以4个典型涡旋流场的涡旋分析验证了识别方法的有效性。 关键词：弯扭张量，流线，涡旋区，涡核区，参数系统 1．引言 涡旋在流场分析和流动理论中一直缺乏明确的定义（formal definition），Lugt (1972)[1] 曾建议将涡旋定义为众多流体微粒围绕一个公共中心的旋转运动，但由于无法预先知晓旋转 中心的定向和运动，使得这样一个似乎与我们的视觉感受一致的定义难以操作，后来 Roth (2000)[2]列出多种域型(region-type)涡旋和线型(line-type)涡核的可操作定义，并针对涡 [3] 轮机内流提出新的涡旋（核）表述，Jiang et al (2004) 选取九种典型的定义，按空间形 态（域型或线型）、伽利略不变性、检测过程的局部性和全局性进行分类，最近 Haller (2004)[4] 还在讨论涡旋的客观性问题。可见学术界关于涡旋的定义远没有定论，但由于它在各种流动 现象中的真实性、在流动（湍流）物理直观中的核心地位、在流场计算数据分析中的必要性 等等，涡旋已经成为是描述流动的最基本概念之一。 涡旋的概念十分复杂。在理论分析中，涡旋常常用速度场的特征量来表征，如涡量、螺 旋度等，它们不仅与速度方向的旋转变化有关、还与速率的分布有关，人们还经常强调涡旋 表征的伽利略不变性。而真实流动显示中的涡旋是曲折流动的流体在当前时刻（或短时间） 留下的某种空间斑图，反映流体相互滑移和滑过的关系，这样的斑图可以是瞬时流线、脉线 和同源物质构型，它们有的与当前速度分布有关、有的与流动的历史有关。从定义来看：（1） 瞬时流线和脉线与参照系有关；（2）同源物质构型在不同惯性参照系下具有同步平移的协变 性；（3）在相对流体系统静止的参照系不同定义具有一致性。定常流动目前是检验各种涡旋 定义的有效性的平台。因为在定常流动中历史的影响完全不反映出来，涡旋的含义在这种流 动中可以各取所需地进行阐释。但是，无论怎样理解，我们真正关心的实际上是对流动演化 动力学起作用的那种涡旋概念（关于这一问题我们将另文进行讨论）。 近期关于流动涡旋研究的一大动力是源于对数值模拟得到的大量数据的后处理，即如何 从这些计算结果提取关于流体运动形态的信息、进行涡旋结构的识别与显示 （Detection and 1 本文工作得到中国自然科学基金（项目编号和中国博士后基金的支持 - 1 - Visualization of Vortices）的问题。通常，需要引入一个或多个特征量（Feature）来建 立一种涡旋识别算法，文献采用过的有涡量ω ∇× u 、螺旋度h V −2 u ⋅ω （V 是速率，也有 用速率与涡量幅值的乘积来标准化的）、速度梯度J ∇u 及其特征值、组合张量Λ S2 =+ Ω 2 （其中S, Ω 分别是速度梯度的对称部分和反对称部分）及其特征值、压力梯度∇p 、加速度 2 2 a du dt 和速度的二阶时间导数b d u dt 等特征量，也有通过追溯流线累积总体参数进行 识别的。类似于此，本文只讨论基于