流场二维空间多尺度测量系统及其测量方法技术方案

本发明专利技术涉及流场解析技术领域，尤其是一种流场二维空间多尺度测量系统及其测量方法，该系统包括入口、出口和底面，该系统中设有二维正交小波计算域，并通过建立系数为12的Daubechies小波基矩阵CN，计算出二维正交小波计算域中瞬态流场流向速度U。本发明专利技术提供了基于空间尺度抽取流场不同尺度湍流结构的方法，客服了传统方法只能抽取流场大尺度结构的弱点，虽然前次提出的一维正交小波变换也能实现对流场的多尺度解析，但是这种方法是基于时间多尺度。通过二维正交小波变换，获取流场中包含的空间多尺度物理信息，为深入研究湍流流场结构提供了新的定性和定量的分析方法。

【技术实现步骤摘要】

本专利技术涉及流场解析
，尤其是一种流场二维空间多尺度测量系统及其测量方法。
技术介绍
湍流是指粘性流体在高雷诺数条件下由于流动失稳而引起的极端复杂的流动状态。作为典型的湍流流场，钝体绕流广泛存在于交通运输、航空航天、水利工程等科学研究和工程实际中。在一定雷诺数范围内，钝体绕流产生的边界层分离和漩涡脱落等复杂流动现象导致阻力、振动和噪音的形成，至今仍然是流体力学领域研究的重点和难点，此外，工程技术中的大量问题与湍流密切相关。从流体力学的观点来看，流体流经钝体表面后会产生由不同尺度漩涡结构组成的湍流流场。作为流体运动的“肌腱”，这些漩涡结构串联了整个流场，对分离区的形成和流动噪音的产生有着决定性的影响。鉴于漩涡运动对湍流流场的重要性，我们需要获取漩涡结构中包含的更细节的物理信息，将其分解成不同的尺度进行研究。作为一种能够识别流场中大尺度湍流结构的解析方法，本征正交分解（POD）被广泛的应用于湍流研究领域。这种方法的主要特点是通过确定最优正交基将流场中不同尺度的结构按照包含能量的大小分解成不同模态，通常被用于分析大尺度的高能流场结构。然而在实际的湍流结构是由不同尺度的漩涡组成，这种多尺度特征不仅体现在时间上，也体现在空间上。各种尺度漩涡结构包含不同大小的能量，除高能的大尺度漩涡外，还存在中小尺度的漩涡结构。在一定条件下，尺度较小的漩涡能够抑制大尺度漩涡的形成和发展，进而减小钝体的阻力和噪音。但是这种方法是基于时间多尺度，在抽取空间流场时，不能严格按照空间尺度为基准对流场进行多尺度解析。在空间上，湍流结构可看成是由不同大小的漩涡结构叠加而成，它们的叠加会使得流场呈现出杂乱无序的状态，这就需要我们从空间上把不同大小的漩涡结构抽取出来。为了克服一维正交小波变换无法准确的抽取流场不同空间尺度漩涡结构的缺点，我们提出了基于二维正交小波的流场空间多尺度解析方法。
技术实现思路
为了克服现有的不足，本专利技术提供了一种流场二维空间多尺度测量系统及其测量方法。本专利技术解决其技术问题所采用的技术方案是：一种流场二维空间多尺度测量系统，该系统包括入口、出口和底面，入口设置为速度入口，系统中速度大小为U0=15m/s，对应的雷诺数大小为10000，出口条件设置为压力出口，该系统中设有二维正交小波计算域。根据本专利技术的另一个实施例，进一步包括，该系统采用非结构性网格进行划分，大小为5LX2LX1.5L，整个系统由542945个网格组成，且两侧面、顶面设置为对称边界条件，底面采用固壁无滑移边界条件。根据本专利技术的另一个实施例，进一步包括，二维正交小波计算域的大小为1.5LX1.5L，由64X64个网格点组成。一种流场二维空间多尺度测量系统的测量方法，首先通过k-ε模型计算得到流场的定常解，得到的定常解通过LES模型进行非定常计算，得到瞬态尾流流场，将二维正交小波计算域中各网格点流场参数用表示，其中i=1,…,nx;j=1,…,ny;nx,ny,的值分别为64和64，参数v代表了流向速度U(x方向)或者纵向(y方向)速度V。，二维正交小波计算域中瞬态流场流向速度U的计算步骤如下：①建立Daubechies小波基矩阵CN，系数设为12；②把瞬态流场中包含的nxXny(64X64)个流向速度U放入大小为2NX2N的二维矩阵UN：(1)在本次计算中上式N的值为6。③首先通过系数为12的Daubechies小波基矩阵CN对初始二维数据矩阵UN沿水平(x)方向进行一维正交小波变换:(2)上式中i=1,…,2N-1,m=1,…,2N，上标T代表了矩阵的转置运算，计算产生了由平顺系数和差分系数组成的列向量和，列向量的大小为2NX1，表示的是第一次分解得到的小波系数矩阵，整个乘积过程相当于沿着水平方向对二维矩阵的列向量进行连续的一维正交小波分解。④通过置换矩阵PN把矩阵中的平顺系数置换到矩阵中的前半部分，差分系数置换到矩阵中的后半部分，这样就完成了对二维矩阵沿着水平方向的初次分解，分解过程如下：(3)⑤把前步得到的转置后可进一步沿垂直方向(y)进行分解就完成了对二维矩阵的初次正交小波变换：(4)从式（3）可以得到由平顺系数ss,差分系数sd,ds和dd组成的小波系数矩阵，为了便于对下一步计算的说明，我们把系数矩阵分为四个部分，第一个部分是由平顺系数ss组成的子带S1，第二部分是由差分系数sd组成的子带，第三部分是由差分系数ds组成的子带，第四部分是由差分系数dd组成的子带，各个子带包含了2N-1X2N-1个元素。由差分系数组成的子带，和被定义为尺度为N-2的小波系数。⑥第二次正交小波变换只针对第一次变换中的子带S1进行，计算过程如下：(5)第二次分解过程也产生了四个不同的子带，各个子带的大小为2N-2X2N-2，其中子带，和被定义为尺度为N-3的小波系数。⑦把得到的光顺系数进一步进行卷积和置换计算，直到最大尺度的小波系数被分解出来，这个过程将重复进行N-3次，最终得到大小为23X23的平顺系数子带SN-3，这个子带被定义为尺度为1的小波系数。下式5中矩阵被称作小波系数矩阵，由N-2个不同尺度的小波系数分量组成，小波分析矩阵W可通过小波基矩阵和置换矩阵的级联算法获得，分解过程示意如图2所示。(6)W=P4C4…PN-1CN-1PNCN,(7)⑧由于小波分析矩阵的正交性(WTW=I，其中I为单位矩阵)，二维离散小波的逆变换可通过下式来表达：(8)⑨基于小波系数矩阵中各个系数分量的空间尺度特性，我们将小波系数分解为各分量之和并对各个分量进行重构：(9).(10)其中称为尺度1小波分量(最大空间尺度)，称为尺度N-2的小波分量(最小空间尺度)。本专利技术的有益效果是，本专利技术提供了基于空间尺度抽取流场不同尺度湍流结构的方法，客服了传统方法只能抽取流场大尺度结构的弱点，虽然前次提出的一维正交小波变换也能实现对流场的多尺度解析，但是这种方法是基于时间多尺度。通过二维正交小波变换，获取流场中包含的空间多尺度物理信息，为深入研究湍流流场结构提供了新的定性和定量的分析方法。附图说明下面结合附图和实施例对本专利技术进一步说明。图1是本专利技术数值计算模型示意图；图2是本专利技术二维正交小波分解示意图；图3是本专利技术不同尺度小波分量空间两点自相关函数分布图；图4是本专利技术瞬态流场流线和对应的涡量云图；图4中a是LES模型计算图；图4本文档来自技高网...

【技术保护点】
一种流场二维空间多尺度测量系统，其特征是，该系统包括入口（1）、出口（2）和底面（3），入口（1）设置为速度入口，系统中速度大小为U0=15m/s，对应的雷诺数大小为10000，出口（2）条件设置为压力出口，该系统中设有二维正交小波计算域（4）。

【技术特征摘要】
1.一种流场二维空间多尺度测量系统，其特征是，该系统包括入口（1）、出口（2）和底面
（3），入口（1）设置为速度入口，系统中速度大小为U0=15m/s，对应的雷诺数大小为10000，出
口（2）条件设置为压力出口，该系统中设有二维正交小波计算域（4）。
2.根据权利要求1所述的流场二维空间多尺度测量系统，其特征是，该系统采用非结构
性网格进行划分，大小为5LX2LX1.5L，整个系统由542945个网格（5）组成，且两侧面、顶
面设置为对称边界条件，底面（3）采用固壁无滑移边界条件。
3.根据权利要求1所述的流场二维空间多尺度测量系统，其特征是，二维正交小波计算
域（4）的大小为1.5LX1.5L，由64X64个网格（5）组成。
4.根据权利要求1所述的流场二维空间多尺度测量系统的测量方法，其特征是，首先通
过k-ε模型计算得到流场的定常解，得到的定常解通过LES模型进行非定常计算，得到瞬态
尾流流场，将二维正交小波计算域（4）中各网格点流场参数用表示，其中i=1,
…,nx;j=1,…,ny;nx,ny,的值分别为64和64，参数v代表了流向速度U(x方向)或
者纵向(y方向)速度V。
5.二维正交小波计算域（4）中瞬态流场流向速度U的计算步骤如下：
①建立Daubechies小波基矩阵CN，系数设为12；
②把瞬态流场中包含的nxXny(64X64)个流向速度U放入大小为2NX2N的二维矩阵
UN：
(1)
在本次计算中上式N的值为6；
③首先通过系数为12的Daubechies小波基矩阵CN对初始二维数据矩阵UN沿水平(x)方
向进行一维正交小波变换:
(2)
上式中i=1,…,2N-1,m=1,…,2N，上标T代表了矩阵的转置运算，计算产生了由平顺系
数和差分系数组成的列向量和，列向量的大小为2NX1，表示的是第一次分解
得到的小波系数矩阵，整个乘积过程相当于沿着水平方向对二维矩阵的列向量进行连续的
一维正...

【专利技术属性】
技术研发人员：郑焱，张丹，贝绍轶，张兰春，李波，韩冰源，王汝佳，
申请(专利权)人：江苏理工学院，
类型：发明
国别省市：江苏;32