基于經(jīng)驗(yàn)?zāi)B(tài)分解的粒子圖像測(cè)速流場(chǎng)重構(gòu)研究

王曉宇,賀靜,柯漢兵,林梅

(1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.武漢第二船舶設(shè)計(jì)研究所,430205,武漢)

粒子圖像測(cè)速(PIV)技術(shù)是一種誕生于20世紀(jì)80年代的無擾動(dòng)流場(chǎng)測(cè)試技術(shù),它克服了傳統(tǒng)接觸式流場(chǎng)測(cè)量中單點(diǎn)測(cè)量的局限性,可以對(duì)存在變化流速及各類渦結(jié)構(gòu)的復(fù)雜流動(dòng)實(shí)現(xiàn)全流場(chǎng)瞬時(shí)非接觸測(cè)量[1-2]。在PIV流場(chǎng)圖像后處理中常采用光學(xué)楊氏條紋法、自相關(guān)法或互相關(guān)法,逐點(diǎn)處理PIV圖像,從而獲得流場(chǎng)速度數(shù)據(jù)[3-4]。由于在實(shí)驗(yàn)中可能存在示蹤粒子分布不均勻、PIV圖像變形或扭曲、PIV圖像背景噪聲、成像系統(tǒng)變形等因素,獲得的流速數(shù)據(jù)中會(huì)存在干擾正常流場(chǎng)分析的錯(cuò)誤數(shù)據(jù)。當(dāng)錯(cuò)誤數(shù)據(jù)造成的誤差超過2%時(shí),便會(huì)影響正常的流場(chǎng)分析研判[5]。

針對(duì)PIV流場(chǎng)數(shù)據(jù)中錯(cuò)誤數(shù)據(jù)的消除方法,國內(nèi)諸多學(xué)者都進(jìn)行了相關(guān)研究。王平讓總結(jié)了錯(cuò)誤數(shù)據(jù)產(chǎn)生的原因,并從PIV圖像濾波處理的角度提出了消除錯(cuò)誤數(shù)據(jù)的方法[6]。王燦星等建立了基于互相關(guān)技術(shù)的連續(xù)PIV圖像之間的粒子對(duì)應(yīng)方法,提出了識(shí)別錯(cuò)誤數(shù)據(jù)的方法和消除準(zhǔn)則[7]。高琪等基于速度場(chǎng)的本征正交分解后處理技術(shù),通過迭代方法有效地實(shí)現(xiàn)了錯(cuò)誤數(shù)據(jù)的識(shí)別與修復(fù)[8]。吳龍華等在對(duì)人腦判別PIV錯(cuò)誤數(shù)據(jù)方式進(jìn)行模擬的基礎(chǔ)上,建立了可以進(jìn)行準(zhǔn)確錯(cuò)誤數(shù)據(jù)識(shí)別的多證據(jù)推理Hopfield神經(jīng)網(wǎng)絡(luò)模型[9]。王晉軍等引入不可壓縮連續(xù)性方程對(duì)PIV測(cè)得的流場(chǎng)數(shù)據(jù)進(jìn)行后處理,使得修正后的流場(chǎng)嚴(yán)格滿足差分形式的連續(xù)性方程,實(shí)現(xiàn)了流場(chǎng)數(shù)據(jù)的準(zhǔn)確修正[10]。阮曉東等采用Delaunay網(wǎng)格技術(shù)建立了基于流體連續(xù)性假設(shè)的自動(dòng)檢測(cè)PIV數(shù)據(jù)中錯(cuò)誤數(shù)據(jù)的算法,實(shí)現(xiàn)了錯(cuò)誤數(shù)據(jù)的識(shí)別與替換[11]。張工利用Visual C++開發(fā)了通用的PIV圖像后處理軟件,實(shí)現(xiàn)了PIV圖像的增強(qiáng)等功能,可以從PIV圖像中消除部分錯(cuò)誤數(shù)據(jù)[12]。王元等廣泛研究了各類提高流場(chǎng)測(cè)量及數(shù)據(jù)后處理精度的粒子追蹤算法,減少了錯(cuò)誤數(shù)據(jù)的產(chǎn)生[13-16]。上述方法對(duì)PIV流場(chǎng)數(shù)據(jù)中錯(cuò)誤數(shù)據(jù)的剔除和修正都有較好的效果,但也存在處理流程復(fù)雜、效率不高、引入非物理偏差等不足。

本研究提出了一種基于經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)的PIV流場(chǎng)數(shù)據(jù)后處理方法,將數(shù)據(jù)分解得到的波形異常的本征模態(tài)分量(IMF)進(jìn)行濾波處理并與其余波形平滑的正常本征模態(tài)分量進(jìn)行反向疊加,重構(gòu)流場(chǎng)數(shù)據(jù),使錯(cuò)誤數(shù)據(jù)得到有效剔除。當(dāng)分解得到的本征模態(tài)分量較多且錯(cuò)誤數(shù)據(jù)集中在某一個(gè)本征模態(tài)分量時(shí),直接將與原始流場(chǎng)數(shù)據(jù)相關(guān)的本征模態(tài)分量反向疊加重構(gòu)流場(chǎng)數(shù)據(jù),摒棄不相關(guān)的本征模態(tài)分量,實(shí)現(xiàn)錯(cuò)誤數(shù)據(jù)的直接剔除。本方法已經(jīng)在多個(gè)模擬或?qū)嶒?yàn)獲取的變化較為緩和的平穩(wěn)流場(chǎng)數(shù)據(jù)處理中得到了良好的使用,本文利用該方法對(duì)2個(gè)典型的人為添加錯(cuò)誤矢量的模擬標(biāo)準(zhǔn)流場(chǎng)進(jìn)行數(shù)據(jù)處理,同時(shí)對(duì)某實(shí)驗(yàn)的原始流場(chǎng)數(shù)據(jù)進(jìn)行修正,以上述3個(gè)流場(chǎng)為例,說明并驗(yàn)證該方法的可行性。

1 基于經(jīng)驗(yàn)?zāi)B(tài)分解的流場(chǎng)數(shù)據(jù)重構(gòu)原理

1.1 經(jīng)驗(yàn)?zāi)B(tài)分解概述

經(jīng)驗(yàn)?zāi)B(tài)分解是Huang等于1998年提出的一種無需基函數(shù)的自適應(yīng)性信號(hào)數(shù)據(jù)處理新方法,在各工程領(lǐng)域得到了廣泛、高效的應(yīng)用[17-18]。該方法利用數(shù)據(jù)自身的尺度特征進(jìn)行數(shù)據(jù)分解,特別適用于非平穩(wěn)、非線性數(shù)據(jù)的平穩(wěn)化處理。經(jīng)驗(yàn)?zāi)B(tài)分解處理信號(hào)的基本步驟如下。

(1)確定隨機(jī)信號(hào)x(t)在整個(gè)時(shí)間尺度上的所有極大值點(diǎn)和極小值點(diǎn),并用三次樣條差值曲線分別連結(jié)所有極大值點(diǎn)和極小值點(diǎn),形成上包絡(luò)線a(t)和下包絡(luò)線b(t),求取上、下包絡(luò)線的平均值

(1)

(2)求取原始隨機(jī)信號(hào)x(t)與m(t)之間的差值

h(t)=x(t)-m(t)

(2)

(3)如果h(t)的極值點(diǎn)與過零點(diǎn)的個(gè)數(shù)最多相差一個(gè)且h(t)的上、下包絡(luò)線的均值為0,則h(t)為第一個(gè)本征模態(tài)分量c1。求取x(t)與c1的差值

r(t)=x(t)-c1

(3)

將r(t)視為最新的隨機(jī)信號(hào),重新開始上述分解過程(1)～(3)求取剩余的本征模態(tài)分量,直至r(t)變?yōu)閱握{(diào)信號(hào)。如果h(t)不能滿足上述條件,則將h(t)視為最新的隨機(jī)信號(hào),重新開始上述分解過程(1)～(3)。

關(guān)于經(jīng)驗(yàn)?zāi)B(tài)分解更加具體的數(shù)學(xué)化描述及證明過程請(qǐng)參見文獻(xiàn)[19-20]。經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解后,一組數(shù)據(jù)可以被分解為若干個(gè)表征不同數(shù)據(jù)尺度和數(shù)據(jù)特征的本征模態(tài)分量,它們是組成原始信號(hào)的不同尺度的一系列分量,是對(duì)原始信號(hào)特征的一系列表達(dá)。理論上連續(xù)的兩個(gè)本征模態(tài)分量的標(biāo)準(zhǔn)差在0.2～0.3之間時(shí),則分解得到的本征模態(tài)分量既能滿足線性條件,穩(wěn)定性良好,又能夠較好地反映原始隨機(jī)信號(hào)的性質(zhì)。經(jīng)驗(yàn)?zāi)B(tài)分解具有正交性和完備性的特征,因此將分解得到的本征模態(tài)分量反向疊加即可得到原始數(shù)據(jù)。在處理現(xiàn)代流動(dòng)測(cè)試技術(shù)所獲得的流場(chǎng)數(shù)據(jù)時(shí),數(shù)據(jù)量的大小并不會(huì)直接影響所分解出的本征模態(tài)分量數(shù)量,而原始流場(chǎng)的復(fù)雜程度則對(duì)本征模態(tài)分量的數(shù)量影響較大。原始數(shù)據(jù)經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解后一般可以得到5～15個(gè)本征模態(tài)分量,其中本征模態(tài)分量在5～10個(gè)時(shí)可以認(rèn)為較少,超過10個(gè)時(shí)可以認(rèn)為分解得到的本征模態(tài)分量較多。

1.2 相關(guān)系數(shù)概述

在數(shù)據(jù)處理中常用相關(guān)系數(shù)來表征兩組數(shù)據(jù)之間的相關(guān)程度,在本研究中通過求解不同的本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)系數(shù),可以判斷該分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)程度。相關(guān)系數(shù)的劃分并沒有固定的標(biāo)準(zhǔn),而是根據(jù)所處理數(shù)據(jù)的具體特征來決定的。表1所示的是一種相關(guān)性的劃分方式,本研究根據(jù)這種相關(guān)系數(shù)大小的劃分方式將本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的關(guān)系歸類為不相關(guān)、微相關(guān)、實(shí)相關(guān)、高度相關(guān)4類,其中微相關(guān)、實(shí)相關(guān)、高度相關(guān)統(tǒng)稱為相關(guān)。

表1 相關(guān)性劃分

1.3 流場(chǎng)數(shù)據(jù)處理流程

圖1是本研究所提出的基于經(jīng)驗(yàn)?zāi)B(tài)分解的PIV流場(chǎng)數(shù)據(jù)后處理方法流程圖,其中實(shí)線表示的流程是通用方法,虛線表示的流程是在原始流場(chǎng)較為復(fù)雜、原始流場(chǎng)數(shù)據(jù)經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解后得到的本征模態(tài)分量較多且錯(cuò)誤數(shù)據(jù)集中在某一個(gè)本征模態(tài)分量時(shí)對(duì)通用方法的一種簡(jiǎn)化。

圖1 基于EMD的PIV流場(chǎng)數(shù)據(jù)后處理流程圖

錯(cuò)誤數(shù)據(jù)在數(shù)值大小和方向上與正確數(shù)據(jù)有明顯的差異,在原始流場(chǎng)數(shù)據(jù)沒有經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解時(shí),錯(cuò)誤數(shù)據(jù)摻雜在正確數(shù)據(jù)中難以識(shí)別。當(dāng)用經(jīng)驗(yàn)?zāi)B(tài)分解將原始流場(chǎng)數(shù)據(jù)分解為若干個(gè)本征模態(tài)分量時(shí),錯(cuò)誤數(shù)據(jù)就會(huì)造成本征模態(tài)分量曲線出現(xiàn)明顯的局部不平滑。由于錯(cuò)誤數(shù)據(jù)造成本征模態(tài)分量曲線局部的不平滑十分明顯,因此常用的濾波方法就可以將錯(cuò)誤數(shù)據(jù)點(diǎn)在本征模態(tài)分量中濾除而不誤刪正確數(shù)據(jù)。根據(jù)經(jīng)驗(yàn)?zāi)B(tài)分解的正交性和完備性,波形異常的本征模態(tài)分量經(jīng)濾波處理后與波形平滑正常的本征模態(tài)分量反向疊加就可以重構(gòu)剔除了錯(cuò)誤數(shù)據(jù)的流場(chǎng)數(shù)據(jù)。上述的原始流場(chǎng)數(shù)據(jù)處理思路可表示為圖1中實(shí)線所示的流程。

原始流場(chǎng)數(shù)據(jù)經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解后得到的本征模態(tài)分量越多,每個(gè)本征模態(tài)分量攜帶的原始流場(chǎng)數(shù)據(jù)就越少,所占原始流場(chǎng)數(shù)據(jù)的比例就越小。在有些實(shí)驗(yàn)中,粒子圖像測(cè)速技術(shù)所獲得的流場(chǎng)數(shù)據(jù)相對(duì)比較復(fù)雜且錯(cuò)誤數(shù)據(jù)比較明顯,這時(shí)分解得到的本征模態(tài)分量較多且錯(cuò)誤數(shù)據(jù)會(huì)集中在某一個(gè)本征模態(tài)分量。摻雜了錯(cuò)誤數(shù)據(jù)的本征模態(tài)分量波形異常且與原始流場(chǎng)數(shù)據(jù)的相關(guān)性較低(呈不相關(guān)),因此在該流程中可以通過求解各本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)性大小來判斷錯(cuò)誤數(shù)據(jù)所在的本征模態(tài)分量。為了提高處理速度和效率,可以按照?qǐng)D1中虛線所示的處理流程,省略相關(guān)性低(呈不相關(guān))的本征模態(tài)分量的濾波過程,直接將該分量刪除并將其余相關(guān)性較高的本征模態(tài)分量反向疊加直接重構(gòu)流場(chǎng)數(shù)據(jù)。隨著錯(cuò)誤數(shù)據(jù)的刪除,也會(huì)損失該本征模態(tài)分量中的正確數(shù)據(jù),但在分解得到的本征模態(tài)分量較多的情況下,損失的小部分正確數(shù)據(jù)不會(huì)對(duì)整體流場(chǎng)數(shù)據(jù)造成明顯的影響,流場(chǎng)的特性仍然可以清楚地表達(dá)。

2 算例驗(yàn)證與應(yīng)用分析

2.1 模擬流場(chǎng)中的驗(yàn)證與分析

圖2是用計(jì)算機(jī)模擬的流速為1 m/s的單向均勻流場(chǎng),在4 000個(gè)原始流場(chǎng)數(shù)據(jù)中加入了明顯的錯(cuò)誤數(shù)據(jù),人為添加的誤差約為1.7%,因此速度矢量場(chǎng)中出現(xiàn)了數(shù)值大小和方向明顯錯(cuò)誤的速度矢量。圖3是原始流場(chǎng)數(shù)據(jù)按圖1中實(shí)線流程處理后得到的速度矢量場(chǎng)。圖4和圖5是對(duì)該原始流場(chǎng)x、y方向全部的4 000個(gè)流場(chǎng)數(shù)據(jù)進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解后得到的本征模態(tài)分量,分別為6個(gè)和8個(gè),不符合使用簡(jiǎn)化流程的要求,因此采用的是通用方法。由圖可知x、y方向的錯(cuò)誤數(shù)據(jù)都集中在了第1個(gè)本征模態(tài)分量中,因此該本征模態(tài)分量的曲線中出現(xiàn)了不平滑的跳變點(diǎn),也正是錯(cuò)誤數(shù)據(jù)存在的位置。對(duì)x、y方向上錯(cuò)誤數(shù)據(jù)集中的第1個(gè)本征模態(tài)分量進(jìn)行濾波處理,即可把錯(cuò)誤數(shù)據(jù)存在的跳變點(diǎn)濾除而不誤刪正確數(shù)據(jù)。將x、y方向上濾波處理后的第1個(gè)本征模態(tài)分量與其他本征模態(tài)分量反向疊加,即可重構(gòu)流場(chǎng)數(shù)據(jù)。圖2與圖3對(duì)比表明,實(shí)線表示的處理流程能有效實(shí)現(xiàn)該單向均勻流場(chǎng)錯(cuò)誤數(shù)據(jù)的剔除和速度矢量場(chǎng)的修正,修正并重構(gòu)后的流場(chǎng)數(shù)據(jù)與正確流場(chǎng)數(shù)據(jù)的誤差約為0.002%。

圖2 原始單向流場(chǎng)數(shù)據(jù)得到的速度矢量場(chǎng)

圖3 處理后的單向流場(chǎng)數(shù)據(jù)得到的速度矢量場(chǎng)

圖4 原始單向流場(chǎng)x方向數(shù)據(jù)的本征模態(tài)分量

圖5 原始單向流場(chǎng)y方向數(shù)據(jù)的本征模態(tài)分量

為了驗(yàn)證所提出的方法在處理相對(duì)復(fù)雜的流場(chǎng)數(shù)據(jù)時(shí)的有效性,本研究模擬了速度大小不定的單向旋轉(zhuǎn)流場(chǎng),并在原始流場(chǎng)數(shù)據(jù)中加入了一些明顯的錯(cuò)誤數(shù)據(jù),人為添加的誤差約為3.3%。圖6和圖7是對(duì)該原始流場(chǎng)x、y方向的流場(chǎng)數(shù)據(jù)進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解后得到的本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)系數(shù),本征模態(tài)分量的數(shù)量和特點(diǎn)符合上文提出的虛線表示的簡(jiǎn)化處理流程的要求,因此處理該流場(chǎng)數(shù)據(jù)時(shí)采用的是簡(jiǎn)化流程。圖8和圖9分別是由原始流場(chǎng)數(shù)據(jù)和按圖1中虛線表示的簡(jiǎn)化流程處理后的流場(chǎng)數(shù)據(jù)得到的速度矢量場(chǎng)。對(duì)比表明,圖1中虛線表示的簡(jiǎn)化處理流程能夠有效地實(shí)現(xiàn)該速度大小不定的單向旋轉(zhuǎn)流場(chǎng)中錯(cuò)誤數(shù)據(jù)的剔除和速度矢量場(chǎng)的修正,修正并重構(gòu)后的流場(chǎng)數(shù)據(jù)與正確流場(chǎng)數(shù)據(jù)的誤差約為0.18%。

圖6 旋轉(zhuǎn)流場(chǎng)x方向流場(chǎng)數(shù)據(jù)的各本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)系數(shù)

圖7 旋轉(zhuǎn)流場(chǎng)y方向流場(chǎng)數(shù)據(jù)的各本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)系數(shù)

圖8 原始旋轉(zhuǎn)流場(chǎng)數(shù)據(jù)得到的速度矢量場(chǎng)

圖9 處理后的旋轉(zhuǎn)流場(chǎng)數(shù)據(jù)得到的速度矢量場(chǎng)

2.2 實(shí)驗(yàn)流場(chǎng)中的應(yīng)用與分析

為了驗(yàn)證本文所提出的基于經(jīng)驗(yàn)?zāi)B(tài)分解的PIV流場(chǎng)數(shù)據(jù)后處理方法在實(shí)驗(yàn)流場(chǎng)數(shù)據(jù)處理中的準(zhǔn)確性與實(shí)用性,本研究對(duì)低速風(fēng)洞中某一T型通道的PIV流場(chǎng)原始數(shù)據(jù)進(jìn)行了處理與分析。關(guān)于該T型通道PIV流動(dòng)測(cè)試實(shí)驗(yàn)的具體細(xì)節(jié)請(qǐng)參見文獻(xiàn)[21-22]。

正如上文所述,實(shí)驗(yàn)流場(chǎng)的復(fù)雜性一般相對(duì)較高,流場(chǎng)數(shù)據(jù)經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解之后得到的本征模態(tài)分量較多。該實(shí)驗(yàn)流場(chǎng)x、y方向的原始流場(chǎng)數(shù)據(jù)經(jīng)過經(jīng)驗(yàn)?zāi)B(tài)分解之后分別得到了10個(gè)本征模態(tài)分量和12個(gè)本征模態(tài)分量,可以用圖1中的簡(jiǎn)化流程進(jìn)行數(shù)據(jù)處理。圖10和圖11分別表示了x方向流場(chǎng)數(shù)據(jù)、y方向流場(chǎng)數(shù)據(jù)的本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)之間的相關(guān)程度。由圖可知,x方向流場(chǎng)數(shù)據(jù)的第1個(gè)本征模態(tài)分量和y方向流場(chǎng)數(shù)據(jù)的第6個(gè)本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的整體相關(guān)性很小(呈不相關(guān)),是錯(cuò)誤數(shù)據(jù)集中的具體表現(xiàn)。該流場(chǎng)數(shù)據(jù)可以按照?qǐng)D1中虛線表示的流程進(jìn)行處理,直接刪除這兩個(gè)本征模態(tài)分量,將x、y方向上剩余的9個(gè)、11個(gè)本征模態(tài)分量反向疊加重構(gòu)流場(chǎng)數(shù)據(jù)。圖12、圖13分別表示了該流場(chǎng)數(shù)據(jù)處理前后的速度矢量場(chǎng),可知原始流場(chǎng)數(shù)據(jù)中位于速度矢量場(chǎng)上部的錯(cuò)誤數(shù)據(jù)(圖12中區(qū)域1)得到了準(zhǔn)確的剔除,一些明顯的錯(cuò)誤矢量得到了有效的消除,速度矢量場(chǎng)更加平滑,品質(zhì)得到了顯著改善,速度矢量場(chǎng)上部的流動(dòng)特征(圖13中區(qū)域1＇)也更加趨近于文獻(xiàn)[22]報(bào)道的各工況下的真實(shí)流動(dòng)情況,但在個(gè)別不存在錯(cuò)誤矢量的區(qū)域出現(xiàn)了流場(chǎng)特征的輕微變化。這是隨著簡(jiǎn)化流程中呈不相關(guān)的本征模態(tài)分量的直接刪除,造成少部分正確數(shù)據(jù)損失所導(dǎo)致的。如果對(duì)局部流場(chǎng)特性的準(zhǔn)確性要求較高,可以不采用簡(jiǎn)化流程,而采用相對(duì)較為煩瑣的通用方法來進(jìn)行處理。

圖10 實(shí)驗(yàn)流場(chǎng)x方向流場(chǎng)數(shù)據(jù)的各本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)系數(shù)

圖11 實(shí)驗(yàn)流場(chǎng)y方向流場(chǎng)數(shù)據(jù)的各本征模態(tài)分量與原始流場(chǎng)數(shù)據(jù)的相關(guān)系數(shù)

圖12 原始實(shí)驗(yàn)流場(chǎng)數(shù)據(jù)得到的速度矢量場(chǎng)

圖13 處理后的實(shí)驗(yàn)流場(chǎng)數(shù)據(jù)得到的速度矢量場(chǎng)

3 結(jié) 論

本研究提出了一種基于經(jīng)驗(yàn)?zāi)B(tài)分解的PIV流場(chǎng)數(shù)據(jù)后處理方法,對(duì)存在錯(cuò)誤數(shù)據(jù)的本征模態(tài)分量進(jìn)行濾波處理或直接刪除,將其余本征模態(tài)分量反向疊加,實(shí)現(xiàn)了錯(cuò)誤數(shù)據(jù)的剔除與流場(chǎng)數(shù)據(jù)的重構(gòu)。利用該方法處理了計(jì)算機(jī)模擬流場(chǎng)和某實(shí)驗(yàn)流場(chǎng)的原始數(shù)據(jù),結(jié)果表明錯(cuò)誤數(shù)據(jù)得到了有效剔除,重構(gòu)的流場(chǎng)數(shù)據(jù)可更加真實(shí)地反映流場(chǎng)特征。

目前任何一種PIV流場(chǎng)數(shù)據(jù)后處理方法都不能完全將錯(cuò)誤數(shù)據(jù)剔除而不影響正確數(shù)據(jù),本研究所提出的方法也是如此,但本方法具備以下兩個(gè)優(yōu)點(diǎn):首先,該方法最大限度地從原始流場(chǎng)數(shù)據(jù)中剔除了錯(cuò)誤數(shù)據(jù),并非通過各類算法在計(jì)算時(shí)“避開”錯(cuò)誤數(shù)據(jù),具備處理的徹底性;其次,該方法簡(jiǎn)單易行,沒有復(fù)雜的錯(cuò)誤數(shù)據(jù)識(shí)別流程與標(biāo)準(zhǔn),可以快捷簡(jiǎn)便地提高流場(chǎng)分析的準(zhǔn)確性。在PIV技術(shù)中,錯(cuò)誤數(shù)據(jù)超過2%就已經(jīng)會(huì)影響流場(chǎng)分析,證明實(shí)驗(yàn)中存在的誤差太大,此時(shí)解決誤差大的最好辦法是重新進(jìn)行實(shí)驗(yàn),而不是消除已有的錯(cuò)誤數(shù)據(jù)。就這一點(diǎn)而言,本方法對(duì)誤差在2%以內(nèi)的流場(chǎng)數(shù)據(jù)是有效的,但本方法也有自身的缺點(diǎn),當(dāng)所處理的流場(chǎng)變化非常劇烈、雜亂(例如射流撞擊壁面)時(shí),本方法容易把部分真實(shí)反映流場(chǎng)的數(shù)據(jù)視為錯(cuò)誤數(shù)據(jù),因此本方法只適用于處理變化較為緩和的平穩(wěn)流場(chǎng)。

受該方法的啟發(fā),也可以將原始的PIV圖像視為一個(gè)二維數(shù)組,利用經(jīng)驗(yàn)?zāi)B(tài)分解的方法通過上述流程來對(duì)其進(jìn)行處理。這時(shí)處理的對(duì)象是PIV原始粒子圖像,各個(gè)本征模態(tài)分量也是以圖像或數(shù)組的形式來呈現(xiàn)的,最終通過本征模態(tài)分量的反向疊加得到處理后的PIV圖像,并從中提取流場(chǎng)數(shù)據(jù)。此外,本方法在三維流場(chǎng)數(shù)據(jù)中的應(yīng)用及可行性還有待驗(yàn)證,在接下來的研究工作中,上述兩部分內(nèi)容將是研究的重點(diǎn)。