CN111414720A - 一种基于神经网络的流场旋涡检测方法 - Google Patents

Abstract

本发明一种基于神经网络的流场旋涡检测方法，先对流场内所有物理网格点打标签，标记其是否属于旋涡区域，得到标签数据，对涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。根据计算网格下的的涡量场，将每个网格上的涡量值减去涡量场的均值后除以标准差，得到归一化后的涡量场；对得到的归一化后的涡量场和标签数据同时随机采样，采样后对区域内的网格点对应的标签进行判断，构建神经网络，使用随机采样得到的带标签局部涡量场数据训练该神经网络，固定神经网络参数，得到固定参数的神经网络；利用第五步得到的固定参数的神经网络，对待检测涡区的流场进行涡区检测。

Description

一种基于神经网络的流场旋涡检测方法

技术领域

本发明涉及一种基于神经网络的流场旋涡检测方法，属于流场物体检测技术领域。

背景技术

旋涡是流体在流动过程产生的一种螺旋状形态，常见的旋涡包括飞机机翼后的涡、龙卷风等。旋涡在工程中发挥重要的作用，在流体力学界被认为是流体运动中的肌腱，飞机飞行时，旋涡往往产生在机翼后方，旋涡的发生会减少飞机的升力，造成飞机失速，产生消极作用；在流体混合的过程中，旋涡能够加速流体的混合，发挥着积极的作用；旋涡还与气动噪声的产生和消除有着密切的联系。计算流体力学(CFD)通过计算机和数值方法来求解流体力学的控制方程，对流体力学问题进行模拟和分析，CFD是迭代计算的，因此一次计算会产生数以万计的流场。流场是指运动流体所占有的空间区域，包括速度场、压力场、密度场等多个不同的物理变量场。流场往往通过网格点划分，本发明中所说的流场内的点就是指流场内的网格点，CFD计算就是通过对初始场进行计算，不断迭代，产生多个中间流场，模拟出流场的变化。CFD一次计算会产生海量的流场数据，流场内产生旋涡的位置很重要，因此需要对所有流场数据进行判断其是否包含旋涡。

目前尚无现有技术针对流场进行涡结构检测，现有技术只有针对图像领域的物体检测，由于流场与图像完全不同，因此无法直接适用于流场旋涡检测。

关于现有技术中已有的涡检测方法包括全局检测方法、局部检测方法、机器学习方法和深度学习方法。全局方法误检和漏检较少，结果较为可信，但是耗时较长，计算流体力学计算过程往往会产生海量的流场，如果用全局方法一一检测，时间太长，无法接受；局部方法简单快速，易于实现，但存在大量的误检和漏检，结果不可靠；机器学习算法高效精确，但通用性不高，泛化性较差，可扩展性太弱；基于深度学习的方法很少，已有的方法结合全局和局部的优势，针对流场内每个点逐一进行判断其是否属于涡区域，但这些方法得到的结果精度比全局方法低，速度比局部方法慢，只是对两种方法性能的折中。已有的深度学习方法具有以下缺点：(1)针对流场内的点进行逐一判断，这会导致数据的重复计算；(2)检测时间与流场内网格点数目成正比；(3)深度学习模型中使用了全连接层，导致网络参数太多，深度神经网络计算一次的时间较长；(4)预处理时，需要将不规则的物理网格转换为规则的计算网格，使用的输入是计算网格下的速度场，丢失了物理网格的特性。

发明内容

本发明解决的技术问题为：克服现有技术不足，提供一种基于神经网络的流场旋涡检测方法，解决流场内涡结构的自动快速检测问题，实现流场内旋涡的高效精确检测。

本发明解决的技术方案为：一种基于分割网络的流场旋涡检测方法，步骤如下：

第一步：对流场内所有物理网格点打标签，标记其是否属于旋涡区域，若属于则标记为1，否则标记为0，得到标签数据，即标记后的物理网格点；

第二步：根据第一步的流场确定出涡量场，对涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。流场的物理网格和计算网格如图3(a)、(b)所示。

第三步，根据第二步得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场；

第四步，对第三步得到的归一化后的涡量场和第一步得到的标签数据同时随机采样，如图4(a)所示，对归一化后的涡量场按照设定的区域大小进行采样，采样后对区域内的网格点对应的标签进行判断，如果区域内标签全部为0或者全部为1，则丢弃该采样数据，剩余的数据为带标签局部涡量场数据；

第五步，构建卷积神经网络CNN，使用第四步的随机采样得到的带标签局部涡量场数据训练该CNN，固定CNN网络参数，得到固定参数的CNN；

第六步，利用第五步得到的固定参数的CNN，对待检测涡区的流场进行涡区检测，步骤如下：

优选的，第一步：对流场内所有物理网格点打标签，标记其是否属于旋涡区域，若属于则标记为1，否则标记为0，得到标签数据，即标记后的物理网格点，具体如下：

步骤1.1，对流场的物理网格中的每个网格点(i,j,k)计算其涡量场,假设用x,y,z表示笛卡尔坐标系中的三个坐标轴方向,则i表示网格点在x方向上处于第i列，j表示网格点在y方向上处于第j行，k表示网格点在z方向处于第k层,i∈{1,2,…,N_x},j∈{1,2,…,N_y},k∈{1,2,…,N_z},N_x,N_y,N_z分别表示在x，y，z方向上网格点的数量。在网格点(i,j,k)处的速度值为v_i,j,k＝(v_xv_yv_z)，v_x,v_y,v_z分别表示在该网格点处的速度在三个方向上的值，网格点处的速度值是流场中已知的变量，无需计算。使用公式(1)，根据流场内网格点(i,j,k)的速度v_i,j,k，计算网格点(i,j,k)处的涡量值ω_i,j,k；

式中，

表示速度分量v_z在y方向的偏导数，使用公式(1-1)计算，

表示速度分量v_y在z方向的偏导数，使用公式(1-2)计算，

表示速度分量v_x在z方向的偏导数，使用公式(1-3)计算，

表示速度分量v_z在x方向的偏导数，使用公式(1-4)计算，

表示速度分量v_y在x方向的偏导数，使用公式(1-5)计算，

表示速度分量v_x在y方向的偏导数，使用公式(1-6)计算。下面分别说明其计算公式。

公式(1-1)～(1-6中)ξ,η,ζ是曲面坐标系下的三个方向。

表示v_x，v_y，v_z对ξ方向求偏导数，计算公式为(1-7)～(1-9)；

表示v_x，v_y，v_z对η方向求偏导数，计算公式为(1-10)～(1-12)；

表示v_x，v_y，v_z对ζ方向求偏导数，计算公式为(1-13)～(1-15)；ξ_x,ξ_y,ξ_z表示曲面坐标系下ξ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(1-16)～(1-18)。η_x,η_y,η_z表示曲面坐标系下η方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(1-19)～(1-21)。ζ_x,ζ_y,ζ_z表示曲面坐标系下ζ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(1-22)～(1-24)。

公式(1-7)～(1-15)中，(i-1,j,k)，(i+1,j,k)表示当前网格点(i,j,k)在x方向上的前一个和后一个网格点，(i,j-1,k)，(i,j+1,k)表示当前网格点(i,j,k)在y方向上的前一个和后一个网格点，(i,j,k-1)，(i,j,k+1)表示当前网格点(i,j,k)在z方向上的前一个和后一个网格点.

表示网格点(i-1,j,k)处在x,y,z三个方向的速度值，

表示网格点(i+1,j,k)处在x,y,z三个方向的速度值。

表示网格点(i,j-1,k)处在x,y,z三个方向的速度值，

表示网格点(i,j+1,k)处在x,y,z三个方向的速度值。

表示网格点(i,j,k-1)处在x,y,z三个方向的速度值，

表示网格点(i,j,k+1)处在x,y,z三个方向的速度值。上述18个速度值均为流场内已知量，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

公式(1-16)～(1-24)中，J^-1可以采用公式(1-25)表示。

公式(1-16)～(1-25)中，x_ξ,y_ξ,z_ξ,x_η,y_η,z_η,x_ζ,y_ζ,z_ζ使用公式(1-26)～(1-34)计算。

公式(1-26)～(1-34)中，x_i-1,j,k，y_i-1,j,k，z_i-1,j,k表示网格点(i-1,j,k)在笛卡尔坐标系的物理坐标值。x_i+1,j,k，y_i+1,j,k，z_i+1,j,k表示网格点(i+1,j,k)在笛卡尔坐标系的物理坐标值。x_i,j-1,k，y_i,j-1,k，z_i,j-1,k表示网格点(i,j-1,k)在笛卡尔坐标系的物理坐标值。x_i,j+1,k，y_i,j+1,k，z_i,j+1,k表示网格点(i,j+1,k)在笛卡尔坐标系的物理坐标值。x_i,j,k-1，y_i,j,k-1，z_i,j,k-1表示网格点(i,j,k-1)在笛卡尔坐标系的物理坐标值。x_i,j,k+1，y_i,j,k+1，z_i,j,k+1表示网格点(i,j,k+1)在笛卡尔坐标系的物理坐标值。上述物理坐标值均可以直接得到，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

至此，计算得到网格点(i,j,k)的涡量值ω_i,j,k，逐个网格点计算后得到所有网格点的涡量值，进而形成涡量场。

步骤1.2，根据步骤1.1得到的涡量场计算平均涡量场。具体计算方式如公式(2)所示。

公式(2)中，ω_μ表示待计算的平均涡量值，ω_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值。N_x,N_y,N_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

使用公式(3)计算瞬时涡量偏差IVD值。

IVD_i,j,k＝|ω_i,j,k-ω_μ| 公式(3)

公式中ω_μ表示待计算的平均涡量值，ω_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值，IVD_i,j,k表示网格点(i,j,k)处的瞬时涡量偏差值。逐个网格点计算后得到所有网格点的IVD值。

步骤1.3，找到流场中IVD值的所有局部最大值，并得到这些局部最大值的位置，获取包含IVD局部最大值的所有等值线，得到这些等值线的凸度偏差值。假设第i条等值线所包围的区域面积为S_i，该等值线包围的区域中必然存在凸包区域，这是已知的定理，假设该凸包区域面积为S_cvi,则可计算出该等值线的凸度偏差值为

计算出所有等值线的凸度偏差值后，得到最大的凸度偏差值对应的等值线，则认为该等值线包围的区域属于涡区。将凸度偏差值最大的等值线包围的所有网格点标记为1，其余位置标记为0，至此获得流场内所有网格点处的标签信息即标签数据。

优选的，第二步：对步骤1.1中得到的涡量场进行网格转换，将物理网格中的网格点一一映射到规则的计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。流场的物理网格和计算网格如图3(a)、(b)所示。

优选的，第三步，根据第二步得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场，具体如下：

根据第二步得到的计算网格下的的涡量场，使用公式(6)与公式(7)计算涡量场的均值ω_μ和标准差ω_δ，将涡量场减去均值后除以标准差，如公式(8)所示，之后使用公式(9)进行归一化操作；

其中ω_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值，N_x,N_y,N_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

根据公式(6)和公式(7)得到的均值和标准差，使用公式(8)计算其去中心化后网格点(i,j,k)处的正则化涡量值

其中i＝{1,2,…,N_x},j＝{1,2,…,N_y},k＝{1,2,…,N_z}。

之后使用公式(9)计算归一化后的涡量值

遍历所有网格点，计算归一化后的涡量值，形成归一化后的涡量场。

优选的，第四步，对第三步得到的归一化后的涡量场和第一步得到的标签数据同时随机采样，对归一化后的涡量场按照设定的区域大小进行采样，采样后对所有采样块内的网格点对应的标签进行判断，如果采样块内标签全部为0或者全部为1，则丢弃该采样数据，剩余的采样块为带标签局部涡量场数据。随机采样如图4(a)所示。

优选的，第五步，构建卷积神经网络CNN，使用第四步的随机采样得到的带标签局部涡量场数据训练CNN，得到固定参数的CNN，具体如下：

假设CNN共有n层，则CNN的函数表达为：

O＝f_n(f_n-1(f_n-2(…f₁(ω,W₁)…),W_n-1),W_n)，其中f₁…f_n表示第一层到第n层的运算，如果第i层为卷积层，则f_i表示卷积运算；W₁,…,W_n表示第一层到第n层的待训练的参数矩阵,O表示CNN的输出，输出维度与采样块大小相同，ω表示输入，即采样涡量块。在确定了层数及对应层的类型后，通过带标签的采样块可以不断调整各层的参数矩阵。长时间训练后，可以固定网络参数，得到固定参数的CNN。

优选的，第六步，利用第五步得到的固定参数的CNN，对待检测涡区的流场进行涡区检测，步骤如下：6.1、根据待检测涡区的流场，确定出涡量场，对该涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。

6.2根据步骤6.1得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场；

6.3对归一化后的涡量场按照设定的区域大小和位置进行顺序采样，如图4(b)所示，得到采样后的局部涡量场数据；

6.4将步骤6.3采样的局部涡量场数据，送入第五步的固定参数的CNN中，得到CNN输出，CNN输出为数据块；根据神经网络输出中每个网格点的值，判断每个网格点否属于涡，若网格点的值为1，则判定属于涡，打上标签，否则判定为不属于涡，打上标签，形成标签块。

6.5将所有标签块按照位置顺序合并，得到整个流场内所有网格点的标签。

其中，6.1、根据待检测涡区的流场的速度值，确定出待检测流场的涡量场，优选方案如下：

对待检测流场的物理网格中的每个网格点(i,j,k)计算其涡量场,假设用x,y,z表示笛卡尔坐标系中的三个坐标轴方向,则i表示网格点在x方向上处于第i列，j表示网格点在y方向上处于第j行，k表示网格点在z方向处于第k层。在网格点处的速度值为v'_i,j,k＝(v'_xv'_yv'_z)，v'_x,v'_y,v'_z分别表示在该网格点处的速度在三个方向上的值，网格点处的速度值是流场中已知的变量，无需计算。使用公式(10)，根据流场内网格点(i,j,k)的速度v'_i,j,k，计算网格点(i,j,k)处的涡量值ω'_i,j,k；

式中，

表示速度分量v'_z在y方向的偏导数，使用公式(10-1)计算，

表示速度分量v'_y在z方向的偏导数，使用公式(10-2)计算，

表示速度分量v'_x在z方向的偏导数，使用公式(10-3)计算，

表示速度分量v'_z在x方向的偏导数，使用公式(10-4)计算，

表示速度分量v'_y在x方向的偏导数，使用公式(10-5)计算，

表示速度分量v'_x在y方向的偏导数，使用公式(10-6)计算。下面分别说明其计算公式。

公式(10-1)～(10-6中)ξ,η,ζ是曲面坐标系下的三个方向。

表示v'_x，v'_y，v'_z对ξ方向求偏导数，计算公式为(10-7)～(10-9)；

表示v'_x，v'_y，v'_z对η方向求偏导数，计算公式为(10-10)～(10-12)；

表示v'_x，v'_y，v'_z对ζ方向求偏导数，计算公式为(10-13)～(10-15)；ξ_x,ξ_y,ξ_z表示曲面坐标系下ξ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(10-16)～(10-18)。η_x,η_y,η_z表示曲面坐标系下η方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(10-19)～(10-21)。ζ_x,ζ_y,ζ_z表示曲面坐标系下ζ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(10-22)～(10-24)。

公式(10-7)～(10-15)中，(i-1,j,k)，(i+1,j,k)表示当前网格点(i,j,k)在x方向上的前一个和后一个网格点，(i,j-1,k)，(i,j+1,k)表示当前网格点(i,j,k)在y方向上的前一个和后一个网格点，(i,j,k-1)，(i,j,k+1)表示当前网格点(i,j,k)在z方向上的前一个和后一个网格点.

表示网格点(i-1,j,k)处在x,y,z三个方向的速度值，

表示网格点(i+1,j,k)处在x,y,z三个方向的速度值。

表示网格点(i,j-1,k)处在x,y,z三个方向的速度值，

表示网格点(i,j+1,k)处在x,y,z三个方向的速度值。

表示网格点(i,j,k-1)处在x,y,z三个方向的速度值，

表示网格点(i,j,k+1)处在x,y,z三个方向的速度值。上述18个速度值均为流场内已知量，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

公式(10-16)～(10-24)中，J^-1可以采用公式(10-25)表示。

公式(10-16)～(10-25)中，x_ξ,y_ξ,z_ξ,x_η,y_η,z_η,x_ζ,y_ζ,z_ζ使用公式(10-26)～(10-34)计算。

公式(10-26)～(10-34)中，x_i-1,j,k，y_i-1,j,k，z_i-1,j,k表示网格点(i-1,j,k)在笛卡尔坐标系的物理坐标值。x_i+1,j,k，y_i+1,j,k，z_i+1,j,k表示网格点(i+1,j,k)在笛卡尔坐标系的物理坐标值。x_i,j-1,k，y_i,j-1,k，z_i,j-1,k表示网格点(i,j-1,k)在笛卡尔坐标系的物理坐标值。x_i,j+1,k，y_i,j+1,k，z_i,j+1,k表示网格点(i,j+1,k)在笛卡尔坐标系的物理坐标值。x_i,j,k-1，y_i,j,k-1，z_i,j,k-1表示网格点(i,j,k-1)在笛卡尔坐标系的物理坐标值。x_i,j,k+1，y_i,j,k+1，z_i,j,k+1表示网格点(i,j,k+1)在笛卡尔坐标系的物理坐标值。上述物理坐标值均可以直接得到，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

至此，计算得到待检测流场在网格点(i,j,k)处的涡量值ω'_i,j,k，逐个网格点计算后得到所有网格点的涡量值，进而形成待检测流场的涡量场。对该涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。

6.2根据步骤6.1得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场，优选方案如下：

根据步骤6.1得到的计算网格下的的涡量场，使用公式(11)与公式(12)计算涡量场的均值ω'_μ和标准差ω'_δ，将涡量场减去均值后除以标准差，如公式(13)所示，之后使用公式(14)进行归一化操作；

其中ω'_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值，N'_x,N'_y,N'_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

根据公式(11)和公式(12)得到的均值和标准差，使用公式(13)计算其去中心化后网格点(i,j,k)处的正则化涡量值

其中i＝{1,2,…,N_x’},j＝{1,2,…,N_y’},k＝{1,2,…,N_z’}。

之后使用公式(14)计算归一化后的涡量值

遍历所有网格点，计算归一化后的涡量值，形成归一化后的涡量场。

6.3对归一化后的涡量场按照设定的区域大小m*m和位置进行顺序采样，得到采样后的局部涡量场数据(不带标签)，顺序采样方式如图4(b)所示。

6.4将步骤6.3采样的局部涡量场数据，送入第五步的固定参数的神经网络中，得到神经网络输出，大小与输入大小一致，均为m*m；根据神经网络输出的值，判断采样块中每个网格点是否属于涡，若网格点的输出值为1，则判定属于涡，打上标签，否则判定为不属于涡，打上标签，形成局部标签块。

6.5将所有标签块按照位置顺序合并，得到整个流场内所有网格点的标签，优选方案如下：

按照采样局部涡量场数据时的顺序，将局部涡量场对应的局部标签块放置到整个流场中，对于重叠的部分使用覆盖策略判断，即总是使用后面标签块中的标签覆盖前面标签块中的标签。最终得到待检测流场内所有网格点的标签。

本发明与现有技术相比的优点在于：

(1)本发明设计了流场中涡量场的计算方式。相比于速度场，涡量场融合了速度信息和物理网格信息，可以保持物理网格特性；

(2)本发明逐块判断减少了数据重复计算，减少了单个流场的检测时间；

(3)本发明提出了涡量场的归一化方式，通过这种归一化方式，能够保证将涡量场去中心化，得到无偏差归一化后的涡量场，有利于神经网络的训练；

(4)本发明首次提出用采样后的局部涡量场做神经网络输入，这使神经网络的输入与流场本身大小无关，因此对于任意大小的流场，均可以采用本发明的方法进行涡区检测，具有较强的泛化性。

附图说明

图1为本发明整体架构图；

图2为本发明训练和测试流程图，其中(a)为训练过程流程示意图，(b)为测试过程流程示意图；

图3为本发明不规则的物理网格图和转换后的规则计算网格图，其中(a)为不规则的物理网格图，(b)为转换后的规则计算网格图；

图4为本发明随机采样和顺序采样示意图，其中(a)为随机采样示意图；(b)为顺序采样示意图；

图5为本发明的方法流程图。

具体实施方式

下面结合附图和具体实施例对本发明做进一步详细描述。

本发明一种基于神经网络的流场旋涡检测方法，对流场内所有物理网格点打标签，标记其是否属于旋涡区域，得到标签数据；根据第一步的流场确定出涡量场，对涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。根据第二步得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场；对第三步得到的归一化后的涡量场和第一步得到的标签数据同时随机采样，对归一化后的涡量场按照设定的区域大小进行采样，采样后对区域内的网格点对应的标签进行判断，得到带标签局部涡量场数据；构建神经网络，使用第四步的随机采样得到的带标签局部涡量场数据训练该神经网络，固定神经网络参数，得到固定参数的神经网络；利用第五步得到的固定参数的神经网络，对待检测涡区的流场进行涡区检测。

本发明所涉及的计算流体力学(CFD)中的旋涡结构非常重要，及时的发现和跟踪能够有效减少一些负面影响，例如飞机飞行时，产生了旋涡会减少飞机的升力。如果能在CFD计算中分析得到在何种情况下产生旋涡，就可以有效避免飞机失速等消极作用的产生。在CFD计算中，为防止丢掉关键信息，总是间隔一定时间步就保存一次结果，因此一次CFD计算可能需要保存几万到几十万个流场。同时，每个流场数据比较大，单个流场数据量最大可达10GB，在如此大量的数据中找哪些时间步存在旋涡是很困难的。已有的分析工具采用传统算法从海量数据中提取旋涡，算法不够智能，时间开销大，结果不够精确。要想快速准确的掌握海量数据规律，必须结合机器学习技术，由机器代替人去挖掘信息，获取知识。本发明结合机器学习技术，采用卷积神经网络CNN，能够高效精确地从大量流场中检测出旋涡，从而大大解放了劳动力，节省了时间。

如图1所示，为本发明采用卷积申请网络进行流场旋涡检测的原理，预处理是在训练过程中步骤1至步骤4的方案，在测试过程中步骤6.1到步骤6.3的的方案，后处理是指在测试过程中步骤6.5的方案，可视化是指将流场的标签数据形成可视的表达方式。

如图5所示，本发明的一种基于分割网络的流场旋涡检测方法，进一步的优方案如下

第一步：对流场内所有物理网格点打标签，标记其是否属于旋涡区域，若属于则标记为1，否则标记为0，得到标签数据，即标记后的物理网格点；

步骤1.1，对流场的物理网格中的每个网格点(i,j,k)计算其涡量场,假设用x,y,z表示笛卡尔坐标系中的三个坐标轴方向,则i表示网格点在x方向上处于第i列，j表示网格点在y方向上处于第j行，k表示网格点在z方向处于第k层,i∈{1,2,…,N_x},j∈{1,2,…,N_y},k∈{1,2,…,N_z},N_x,N_y,N_z分别表示在x，y，z方向上网格点的数量。在网格点(i,j,k)处的速度值为v_i,j,k＝(v_xv_yv_z)，v_x,v_y,v_z分别表示在该网格点处的速度在三个方向上的值，网格点处的速度值是流场中已知的变量，无需计算。使用公式(15)，根据流场内网格点(i,j,k)的速度v_i,j,k，计算网格点(i,j,k)处的涡量值ω_i,j,k；

式中，

表示速度分量v_z在y方向的偏导数，使用公式(15-1)计算，

表示速度分量v_y在z方向的偏导数，使用公式(15-2)计算，

表示速度分量v_x在z方向的偏导数，使用公式(15-3)计算，

表示速度分量v_z在x方向的偏导数，使用公式(15-4)计算，

表示速度分量v_y在x方向的偏导数，使用公式(15-5)计算，

表示速度分量v_x在y方向的偏导数，使用公式(15-6)计算。下面分别说明其计算公式。

公式(15-1)～(15-6)中ξ,η,ζ是曲面坐标系下的三个方向。

表示v_x，v_y，v_z对ξ方向求偏导数，计算公式为(15-7)～(15-9)；

表示v_x，v_y，v_z对η方向求偏导数，计算公式为(15-10)～(15-12)；

表示v_x，v_y，v_z对ζ方向求偏导数，计算公式为(15-13)-(15-15)；ξ_x,ξ_y,ξ_z表示曲面坐标系下ξ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(15-16)～(15-18)。η_x,η_y,η_z表示曲面坐标系下η方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(15-19)-(15-21)。ζ_x,ζ_y,ζ_z表示曲面坐标系下ζ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(15-22)-(15-24)。

公式(15-7)～(15-15)中，(i-1,j,k)，(i+1,j,k)表示当前网格点(i,j,k)在x方向上的前一个和后一个网格点，(i,j-1,k)，(i,j+1,k)表示当前网格点(i,j,k)在y方向上的前一个和后一个网格点，(i,j,k-1)，(i,j,k+1)表示当前网格点(i,j,k)在z方向上的前一个和后一个网格点.

表示网格点(i-1,j,k)处在x,y,z三个方向的速度值，

表示网格点(i+1,j,k)处在x,y,z三个方向的速度值。

表示网格点(i,j-1,k)处在x,y,z三个方向的速度值，

表示网格点(i,j+1,k)处在x,y,z三个方向的速度值。

表示网格点(i,j,k-1)处在x,y,z三个方向的速度值，

表示网格点(i,j,k+1)处在x,y,z三个方向的速度值。上述18个速度值均为流场内已知量，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

公式(15-16)～(15-24)中，J^-1可以采用公式(15-25)表示。

公式(15-16)～(15-25)中，x_ξ,y_ξ,z_ξ,x_η,y_η,z_η,x_ζ,y_ζ,z_ζ使用公式(15-26)～(15-34)计算。

公式(15-26)～(15-34)中，x_i-1,j,k，y_i-1,j,k，z_i-1,j,k表示网格点(i-1,j,k)在笛卡尔坐标系的物理坐标值。x_i+1,j,k，y_i+1,j,k，z_i+1,j,k表示网格点(i+1,j,k)在笛卡尔坐标系的物理坐标值。x_i,j-1,k，y_i,j-1,k，z_i,j-1,k表示网格点(i,j-1,k)在笛卡尔坐标系的物理坐标值。x_i,j+1,k，y_i,j+1,k，z_i,j+1,k表示网格点(i,j+1,k)在笛卡尔坐标系的物理坐标值。x_i,j,k-1，y_i,j,k-1，z_i,j,k-1表示网格点(i,j,k-1)在笛卡尔坐标系的物理坐标值。x_i,j,k+1，y_i,j,k+1，z_i,j,k+1表示网格点(i,j,k+1)在笛卡尔坐标系的物理坐标值。上述物理坐标值均可以直接得到，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

至此，计算得到网格点(i,j,k)的涡量值ω_i,j,k，逐个网格点计算后得到所有网格点的涡量值，进而形成涡量场。

步骤1.2，根据步骤1.1得到的涡量场计算平均涡量场。具体计算方式如公式(16)所示。

公式(16)中，ω_μ表示待计算的平均涡量值，ω_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值。N_x,N_y,N_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

使用公式(17)计算瞬时涡量偏差IVD值。

IVD_i,j,k＝|ω_i,j,k-ω_μ| 公式(17)

公式中ω_μ表示待计算的平均涡量值，ω_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值，IVD_i,j,k表示网格点(i,j,k)处的瞬时涡量偏差值。逐个网格点计算后得到所有网格点的IVD值。

步骤1.3，找到流场中IVD值的所有局部最大值，并得到这些局部最大值的位置，获取包含IVD局部最大值的所有等值线，得到这些等值线的凸度偏差值。假设第i条等值线所包围的区域面积为S_i，该等值线包围的区域中必然存在凸包区域，这是已知的定理，假设该凸包区域面积为S_cvi,则可计算出该等值线的凸度偏差值为

计算出所有等值线的凸度偏差值后，得到最大的凸度偏差值对应的等值线，则认为该等值线包围的区域属于涡区。在具体实现过程中，可以直接调用Matlab仿真软件中的convhull函数计算凸度偏差值。将凸度偏差值最大的等值线包围的所有网格点标记为1，其余位置标记为0，至此获得流场内所有网格点处的标签信息即标签数据。

第二步：第二步：对步骤1.1中得到的涡量场进行网格转换，将物理网格中的网格点一一映射到规则的计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。

第三步，根据第二步得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场，具体如下：

根据第二步得到的计算网格下的的涡量场，使用公式(18)与公式(19)计算涡量场的均值ω_μ和标准差ω_δ，将涡量场减去均值后除以标准差，如公式(20)所示，之后使用公式(21)进行归一化操作；

其中ω_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值，N_x,N_y,N_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

根据公式(18)和公式(19)得到的均值和标准差，使用公式(20)计算其去中心化后网格点(i,j,k)处的正则化涡量值

其中i＝{1,2,…,N_x},j＝{1,2,…,N_y},k＝{1,2,…,N_z}。

之后使用公式(21)计算归一化后的涡量值

遍历所有网格点，计算归一化后的涡量值，形成归一化后的涡量场。

第四步，第四步，对第三步得到的归一化后的涡量场和第一步得到的标签数据同时随机采样，对归一化后的涡量场按照设定的区域大小进行采样，采样后对所有采样块内的网格点对应的标签进行判断，如果采样块内标签全部为0或者全部为1，则丢弃该采样数据，剩余的采样块为带标签局部涡量场数据。

第五步，构建卷积神经网络CNN，使用第四步的随机采样得到的带标签局部涡量场数据训练CNN，得到固定参数的CNN，如图2(a)所示，优选方案如下：

假设CNN共有n层，则CNN的函数表达为：

O＝f_n(f_n-1(f_n-2(…f₁(ω,W₁)…),W_n-1),W_n)，其中f₁…f_n表示第一层到第n层的运算，如果第i层为卷积层，则f_i表示卷积运算；W₁,…,W_n表示第一层到第n层的待训练的参数矩阵,O表示CNN的输出，输出维度与采样块大小相同，ω表示输入，即采样涡量块。在确定了层数及对应层的类型后，通过带标签的采样块可以不断调整各层的参数矩阵。长时间训练后，可以固定网络参数，得到固定参数的CNN。

第六步，利用第五步得到的固定参数的CNN，对待检测涡区的流场进行涡区检测，如图2(b)所示，步骤如下：

6.1、根据待检测涡区的流场，确定出涡量场，对该涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。

对待检测流场的物理网格中的每个网格点(i,j,k)计算其涡量场,假设用x,y,z表示笛卡尔坐标系中的三个坐标轴方向,则i表示网格点在x方向上处于第i列，j表示网格点在y方向上处于第j行，k表示网格点在z方向处于第k层。在网格点处的速度值为v'_i,j,k＝(v'_xv'_yv'_z)，v'_x,v'_y,v'_z分别表示在该网格点处的速度在三个方向上的值，网格点处的速度值是流场中已知的变量，无需计算。使用公式(10)，根据流场内网格点(i,j,k)的速度v'_i,j,k，计算网格点(i,j,k)处的涡量值ω'_i,j,k；

式中，

表示速度分量v'_z在y方向的偏导数，使用公式(22-1)计算，

表示速度分量v'_y在z方向的偏导数，使用公式(22-2)计算，

表示速度分量v'_x在z方向的偏导数，使用公式(22-3)计算，

表示速度分量v'_z在x方向的偏导数，使用公式(22-4)计算，

表示速度分量v'_y在x方向的偏导数，使用公式(22-5)计算，

表示速度分量v'_x在y方向的偏导数，使用公式(22-6)计算。下面分别说明其计算公式。

公式(22-1)～(22-6中)ξ,η,ζ是曲面坐标系下的三个方向。

表示v'_x，v'_y，v'_z对ξ方向求偏导数，计算公式为(22-7)～(22-9)；

表示v'_x，v'_y，v'_z对η方向求偏导数，计算公式为(22-10)～(22-12)；

表示v'_x，v'_y，v'_z对ζ方向求偏导数，计算公式为(22-13)～(22-15)；ξ_x,ξ_y,ξ_z表示曲面坐标系下ξ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(22-16)-(22-18)。η_x,η_y,η_z表示曲面坐标系下η方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(22-19)-(22-21)。ζ_x,ζ_y,ζ_z表示曲面坐标系下ζ方向对笛卡尔坐标系三个方向x，y，z的偏导数，计算公式为公式(22-22)-(22-24)。

公式(22-7)～(22-15)中，(i-1,j,k)，(i+1,j,k)表示当前网格点(i,j,k)在x方向上的前一个和后一个网格点，(i,j-1,k)，(i,j+1,k)表示当前网格点(i,j,k)在y方向上的前一个和后一个网格点，(i,j,k-1)，(i,j,k+1)表示当前网格点(i,j,k)在z方向上的前一个和后一个网格点。

表示网格点(i-1,j,k)处在x,y,z三个方向的速度值，

表示网格点(i+1,j,k)处在x,y,z三个方向的速度值。

表示网格点(i,j-1,k)处在x,y,z三个方向的速度值，

表示网格点(i,j+1,k)处在x,y,z三个方向的速度值。

表示网格点(i,j,k-1)处在x,y,z三个方向的速度值，

表示网格点(i,j,k+1)处在x,y,z三个方向的速度值。上述18个速度值均为流场内已知量，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

公式(22-16)～(22-24)中，J^-1可以采用公式(22-25)表示。

公式(22-16)～(22-25)中，x_ξ,y_ξ,z_ξ,x_η,y_η,z_η,x_ζ,y_ζ,z_ζ使用公式(22-26)-(22-34)计算。

公式(22-26)～(22-34)中，x_i-1,j,k，y_i-1,j,k，z_i-1,j,k表示网格点(i-1,j,k)在笛卡尔坐标系的物理坐标值。x_i+1,j,k，y_i+1,j,k，z_i+1,j,k表示网格点(i+1,j,k)在笛卡尔坐标系的物理坐标值。x_i,j-1,k，y_i,j-1,k，z_i,j-1,k表示网格点(i,j-1,k)在笛卡尔坐标系的物理坐标值。x_i,j+1,k，y_i,j+1,k，z_i,j+1,k表示网格点(i,j+1,k)在笛卡尔坐标系的物理坐标值。x_i,j,k-1，y_i,j,k-1，z_i,j,k-1表示网格点(i,j,k-1)在笛卡尔坐标系的物理坐标值。x_i,j,k+1，y_i,j,k+1，z_i,j,k+1表示网格点(i,j,k+1)在笛卡尔坐标系的物理坐标值。上述物理坐标值均可以直接得到，无需计算。Δξ表示网格点(i-1,j,k)和网格点(i+1,j,k)在曲面坐标系下ξ方向的距离，Δη表示网格点(i,j-1,k)和网格点(i,j+1,k)在曲面坐标系下η方向的距离，Δζ表示网格点(i,j,k-1)和网格点(i,j,k+1)在曲面坐标系下ζ方向的距离，Δξ，Δη，Δζ都默认设置为1。

至此，计算得到待检测流场在网格点(i,j,k)处的涡量值ω'_i,j,k，逐个网格点计算后得到所有网格点的涡量值，进而形成待检测流场的涡量场。对该涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场。

6.2.根据步骤6.1得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场，优选方案如下：

根据步骤6.1得到的计算网格下的的涡量场，使用公式(23)与公式(24)计算涡量场的均值ω'_μ和标准差ω'_δ，将涡量场减去均值后除以标准差，如公式(25)所示，之后使用公式(26)进行归一化操作；

其中ω'_i,j,k表示步骤1.1中得到的网格点(i,j,k)的涡量值，N'_x,N'_y,N'_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

根据公式(23)和公式(24)得到的均值和标准差，使用公式(25)计算其去中心化后网格点(i,j,k)处的正则化涡量值

其中i＝{1,2,…,N_x’},j＝{1,2,…,N_y’},k＝{1,2,…,N_z’}。

之后使用公式(26)计算归一化后的涡量值

遍历所有网格点，计算归一化后的涡量值，形成归一化后的涡量场。

6.3对归一化后的涡量场按照设定的区域大小和位置进行顺序采样，位置选择与区域大小相关，由人为指定，得到采样后的局部涡量场数据；

6.4将步骤6.3采样的局部涡量场数据，送入第五步的固定参数的神经网络中，得到神经网络输出，大小与输入大小一致，均为m*m；根据神经网络输出的值，判断采样块中每个网格点是否属于涡，若网格点的输出值为1，则判定属于涡，打上标签，否则判定为不属于涡，打上标签，形成局部标签块。

6.5将所有标签块按照位置顺序合并，得到整个流场内所有网格点的标签，优选方案如下：

按照采样局部涡量场数据时的顺序，将局部涡量场对应的局部标签块放置到整个流场中，对于重叠的部分使用覆盖策略判断，即总是使用后面标签块中的标签覆盖前面标签块中的标签。最终得到待检测流场内所有网格点的标签。

本发明进一步优选的方案为：对步骤一中的公式(2)进行改进。流场中一般存在物体，例如机翼等，这些物体的存在对流场的影响巨大。因此在设计物理网格的时候，这些物体的边界处网格密度较大，即相邻物理网格点的距离较小，远离物体边界的空间网格密度较小，也就是相邻网格间距较大。这种改进使整个方法集中检测重要性更大的区域中的旋涡，具备更强的实际应用性。

考虑到网格点的这种不平等性质，可以在计算平均涡量场时加入贡献值。如公式(27)所示。

α_i,j,k表示网格点P_i,j,k的贡献值，也称为权重。

d(P_i,j,k,B)表示网格点P_i,j,k和边界点集合B的距离，该距离定义如公式(29)。

边界点集合B在流场和物体给定的情况下可以直接得到。公式(29)表示，网格点P_i,j,k与物体边界集合的距离为该网格点与所有物体边界点距离的最小值。

本发明进一步优选的方案为：对步骤三中的公式(6)～公式(8)进行改进。在计算平均涡量场时加入网格点的贡献值。

α_i,j,k表示网格点P_i,j,k的贡献值，也称为权重。

d(P_i,j,k,B)表示网格点P_i,j,k和边界点集合B的距离，该距离定义如公式(29)。

边界点集合B在流场和物体给定的情况下可以直接得到。公式(34)表示，网格点P_i,j,k与物体边界集合的距离为该网格点与所有物体边界点距离的最小值。

本发明的一种基于分割网络的流场旋涡检测方法，进行了仿真，充分体现了本发明的技术效果，以双三角翼流场为例。该双三角翼流场大小为301*201*101，流场的雷诺数为1000000，双三角翼的攻角为22.5°。对该案例进行CFD仿真，可以得到上万个中间流场，我们选择最后一个流场进行旋涡检测。使用精度、召回率和执行时间作为判断依据，精度和召回率如公式(35)和(36)所示：

p表示精度，其物理含义为我们的方法判断为涡区的网格点中判断正确的比例是多少，比例越大越好。r表示召回率，物理含义为，本发明找到的涡区网格点数量占所有属于涡区的网格点数量的比例，该比例越大越好。

公式(35)和(36)中，TP表示真实标签为1，且使用本发明判断为1的网格点数量；FP表示真实标签为0，本发明判断为1的网格点数量；FN表示真实标签为0，本发明判断为1的网格点数量。这些数量可以统计出来。

本发明在上述双三角翼流场中涡检测精度达到90％，与背景技术中的工作相比提升了6个百分点，召回率r与相关工作比较提升了5个百分点，达到86％。同时，本发明执行时间较短，仅为150秒，与相关工作比较，加速比达到10倍。

本发明设计了流场中涡量场的计算方式。相比于速度场，涡量场融合了速度信息和物理网格信息，可以保持物理网格特性，而且本发明逐块判断减少了数据重复计算，减少了单个流场的检测时间；进一步的本发明提出了涡量场的归一化方式，通过这种归一化方式，能够保证将涡量场去中心化，得到无偏差归一化后的涡量场，有利于神经网络的训练；

本发明首次提出用采样后的局部涡量场做神经网络输入，这使神经网络的输入与流场本身大小无关，因此对于任意大小的流场，均可以采用本发明的方法进行涡区检测，具有较强的泛化性。

Claims (10)

1.一种基于神经网络的流场旋涡检测方法，其特征在于步骤如下：

第一步：对流场内所有物理网格点打标签，标记其是否属于旋涡区域，若属于则标记为1，否则标记为0，得到标签数据，即标记后的物理网格点；

第二步：根据第一步的流场确定出涡量场，对涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场；

第三步，根据第二步得到的计算网格下的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场；

第四步，对第三步得到的归一化后的涡量场和第一步得到的标签数据同时随机采样，对归一化后的涡量场按照设定的区域大小进行采样，采样后对区域内的网格点对应的标签进行判断，得到带标签局部涡量场数据；

第五步，构建卷积神经网络CNN，使用第四步的随机采样得到的带标签局部涡量场数据训练该CNN，从而固定CNN网络参数，得到固定参数的CNN；

第六步，利用第五步得到的固定参数的CNN，对待检测涡区的流场进行涡区检测。

2.根据权利要求1所述的一种基于神经网络的流场旋涡检测方法，其特征在于：第一步：对流场内所有物理网格点打标签，标记其是否属于旋涡区域，若属于则标记为1，否则标记为0，得到标签数据，即标记后的物理网格点，具体如下：

步骤1.1，对流场的物理网格中的每个网格点(i,j,k)计算其涡量场；

步骤1.2，根据步骤1.1得到的涡量场，计算平均涡量场；根据平均涡量场，计算流场IVD值；

步骤1.3，找到流场中IVD值的所有局部最大值，并得到这些局部最大值的位置，获取包含IVD局部最大值的所有等值线，得到这些等值线的凸度偏差值，将凸度偏差值最大的等值线包围的所有网格点进行标记，至此获得流场内所有网格点处的标签信息即标签数据。

3.根据权利要求2所述的一种基于神经网络的流场旋涡检测方法，其特征在于：步骤1.3中，将凸度偏差值最大的等值线包围的所有网格点进行标记，至此获得流场内所有网格点处的标签信息即标签数据，具体为：

将凸度偏差值最大的等值线包围的所有网格点标记为1，其余位置标记为0，至此获得流场内所有网格点处的标签信息即标签数据。

4.根据权利要求1所述的一种基于神经网络的流场旋涡检测方法，其特征在于：第三步中，根据第二步得到的计算网格下的涡量场，计算该涡量场的均值ω_μ，公式如下：

式中，ω_i,j,k表示网格点(i,j,k)的涡量值，N_x,N_y,N_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

5.根据权利要求1所述的一种基于神经网络的流场旋涡检测方法，其特征在于：第三步中，根据第二步得到的计算网格下的的涡量场，计算该涡量场的标准差，公式如下：

式中，ω_i,j,k表示网格点(i,j,k)的涡量值，N_x,N_y,N_z分别表示流场在x，y，z方向上网格点的维度，是已知的量。

6.根据权利要求1所述的一种基于神经网络的流场旋涡检测方法，其特征在于：第四步中，采样后对区域内的网格点对应的标签进行判断，得到带标签局部涡量场数据，具体为：

采样后对所有采样块内的网格点对应的标签进行判断，如果采样块内标签全部为0或者全部为1，则丢弃该采样数据，剩余的采样块为带标签局部涡量场数据。

7.根据权利要求1所述的一种基于神经网络的流场旋涡检测方法，其特征在于：第五步，构建卷积神经网络CNN，使用第四步的随机采样得到的带标签局部涡量场数据训练CNN，得到固定参数的CNN，具体如下：

首先确定CNN层数、CNN的函数表达式以及确定CNN的输入和输出，对于CNN的输出，输出维度与采样块大小相同，在确定了CNN层数及对应层的类型后，通过带标签的采样块不断调整各层的参数矩阵，训练后，固定网络参数，得到固定参数的CNN。

8.根据权利要求1所述的一种基于神经网络的流场旋涡检测方法，其特征在于：第六步，利用第五步得到的固定参数的CNN，对待检测涡区的流场进行涡区检测，步骤如下：

6.1、根据待检测涡区的流场，确定出涡量场，对该涡量场进行网格转换，将物理网格中的网格点一一映射到计算网格的网格点，抛弃物理坐标信息，得到计算网格下的涡量场；

6.2根据步骤6.1得到的计算网格下的的涡量场，计算该涡量场的均值和标准差；将每个网格上的涡量值减去涡量场的均值后除以标准差，再进行归一化操作，得到归一化后的涡量场；

6.3对归一化后的涡量场按照设定的区域大小和位置进行顺序采样，得到采样后的局部涡量场数据；

6.4将步骤6.3采样的局部涡量场数据，送入第五步的固定参数的CNN中，得到CNN输出，CNN输出为数据块；根据神经网络输出中每个网格点的值，判断每个网格点否属于涡，若网格点的值为1，则判定属于涡，打上标签，否则判定为不属于涡，打上标签，形成标签块；

6.5将所有标签块按照位置顺序合并，得到整个流场内所有网格点的标签。

9.根据权利要求8所述的一种基于神经网络的流场旋涡检测方法，其特征在于：6.4中，判断每个网格点否属于涡，打上标签，形成标签块，具体如下：

判断每个网格点否属于涡，若网格点的值为1，则判定属于涡，打上标签，否则判定为不属于涡，打上标签，形成标签块。

10.根据权利要求8所述的一种基于神经网络的流场旋涡检测方法，其特征在于：6.5中，将所有标签块按照位置顺序合并，得到整个流场内所有网格点的标签，具体如下：

按照采样局部涡量场数据时的顺序，将局部涡量场对应的局部标签块放置到整个流场中，对于重叠的部分使用设定的覆盖策略判断，即总是使用后面标签块中的标签覆盖前面标签块中的标签；最终得到待检测流场内所有网格点的标签。

Images