基于人工神經(jīng)網(wǎng)絡的非結構網(wǎng)格尺度控制方法1)

王年華 , 魯 鵬 * 常興華 張來平 鄧小剛

* (中國空氣動力研究與發(fā)展中心空氣動力學國家重點實驗室,四川綿陽 621000)

? (西南科技大學信息工程學院,四川綿陽 621010)

** (重慶文理學院智能制造工程學院,重慶 402160)

?? (國防科技創(chuàng)新研究院無人系統(tǒng)技術研究中心,北京 100071)

*** (軍事科學院,北京 100091)

引言

網(wǎng)格生成是計算流體力學(computational fluid dynamics,CFD)數(shù)值計算的第一步,也是未來CFD六大重要研究領域之一[1-2].在現(xiàn)代CFD 應用過程中,自動生成復雜構型的高質量網(wǎng)格(包括網(wǎng)格自適應)依然是一個重大挑戰(zhàn)性問題.自動化程度和網(wǎng)格質量是網(wǎng)格生成過程中最重要的兩個問題[3-4].據(jù)統(tǒng)計,網(wǎng)格生成通常占據(jù)整個計算周期大約60%的人力時間,高度自動化的網(wǎng)格生成方法無疑可以很大程度節(jié)約CFD 計算周期內(nèi)的人工成本.

網(wǎng)格空間尺度分布控制在網(wǎng)格生成中至關重要,對于網(wǎng)格質量和流動求解精度影響較大.通常希望在幾何曲率大、流動梯度大等重點區(qū)域加密網(wǎng)格,而在非重點區(qū)域則希望網(wǎng)格盡可能稀疏并均勻過渡.傳統(tǒng)的非結構網(wǎng)格尺度控制方法主要有:函數(shù)指定法、插值類方法、背景網(wǎng)格法和根據(jù)流場特征進行自適應等.

函數(shù)指定法適用于簡單問題的全場網(wǎng)格尺度控制,比如可以指定線性函數(shù)控制翼型的網(wǎng)格尺度分布;或用于局部網(wǎng)格尺度控制,比如可以指定幾何級數(shù)或指數(shù)函數(shù)對局部網(wǎng)格進行加密[5].

采用背景網(wǎng)格法[6-9]進行網(wǎng)格尺度控制,需要先用規(guī)則的矩形結構網(wǎng)格或采用較稀疏的非結構網(wǎng)格覆蓋全計算域,在計算域內(nèi)分布一定的“點源”“線源”或者“面源”等局部網(wǎng)格分布控制參數(shù),將這些源視作離散的“熱源”,求解熱傳導方程(泊松方程),得到的穩(wěn)態(tài)解即為全計算域的尺度控制參數(shù)分布[6].或者根據(jù)局部幾何特征(曲率、狹縫、窄邊等信息)確定局部網(wǎng)格尺度,再在背景網(wǎng)格上求解網(wǎng)格尺度滿足的梯度限制方程,將局部網(wǎng)格尺度光滑到全場,得到背景網(wǎng)格上的網(wǎng)格尺度分布[7-9].在生成網(wǎng)格時,根據(jù)控制空間中某點在矩形背景網(wǎng)格中所處的單元,通過背景網(wǎng)格的尺度分布插值得到該點處的網(wǎng)格尺度.

除此之外,還有一些學者提出一些其他形式的背景網(wǎng)格法.如Deister 等[10]提出由最大最小尺度及最大曲率角對幾何進行柵格化,計算得到局部網(wǎng)格尺度,并存儲在自適應的背景笛卡爾網(wǎng)格上.Quadros等[11-12]提出采用幾何體離散骨架的幾何臨近信息、特征尺寸、邊界曲率來測量幾何復雜度,并根據(jù)幾何復雜度生成點源,根據(jù)點源確定網(wǎng)格尺度分布,最終將網(wǎng)格尺度存儲在叉樹結構的笛卡爾背景網(wǎng)格上.Ruiz-Girones 等[13]提出通過在背景網(wǎng)格上求解一種新的非線性方程來控制四邊形網(wǎng)格尺度,等等.

徑向基函數(shù)(radial basis function,RBF)[14-15]可用于數(shù)據(jù)插值,在已知邊界上的網(wǎng)格尺度后,可用RBF 方法將邊界網(wǎng)格尺度插值到內(nèi)場,采用貪婪算法還可以一定程度提高RBF 插值的效率,是一種簡單高效的插值類尺度控制方法.

根據(jù)流動特征物理量的梯度量等判據(jù)進行網(wǎng)格自適應[16-17]也是控制網(wǎng)格尺度分布的一種有效方法,能夠根據(jù)流場變化控制網(wǎng)格疏密,在梯度大的區(qū)域生成更密的網(wǎng)格,能夠更精細地捕捉流動特征,具有更高的計算精度和適應性,是來重要發(fā)展方向之一[18].

近年來,基于人工神經(jīng)網(wǎng)絡的深度學習方法在工業(yè)社會甚至流體力學領域得到廣泛研究和應用[19-21].大數(shù)據(jù)驅動的人工智能方法成為除理論分析、數(shù)值計算和實驗技術以外一種新的研究范式,為各個領域帶來了新的研究思路和方法.

在網(wǎng)格生成領域,經(jīng)過多年的工程實踐,已經(jīng)積累了大量各種類型的網(wǎng)格數(shù)據(jù),這些數(shù)據(jù)包含了網(wǎng)格生成規(guī)則及技術人員在網(wǎng)格生成方面的知識和經(jīng)驗,是天然的機器學習訓練樣本.通過機器學習對網(wǎng)格生成規(guī)則進行學習,可以簡化傳統(tǒng)算法,提高網(wǎng)格生成效率[22].而網(wǎng)格尺度分布的控制也需要技術人員根據(jù)對流動問題的分析和經(jīng)驗合理確定,采用機器學習方法對網(wǎng)格尺度分布進行控制,有望減少人工工作量和對人工經(jīng)驗的依賴.

本文從網(wǎng)格質量、效率、靈活性和自動化程度4 個方面綜合分析各類網(wǎng)格尺度控制方法的優(yōu)缺點.為了克服傳統(tǒng)背景網(wǎng)格插值法效率低、自動化程度不高等方面的不足,本文從效率和自動化程度角度提出兩種網(wǎng)格尺度控制方法,首先將RBF 插值方法應用于網(wǎng)格尺度控制,采用貪婪算法對RBF 插值參考點序列進行精簡,實現(xiàn)高效的RBF 網(wǎng)格尺度分布控制方法.進一步將提出一種采用人工神經(jīng)網(wǎng)絡進行非結構網(wǎng)格尺度控制的方法,通過引入相對壁面距離和相對網(wǎng)格尺度,初步確定合理的神經(jīng)網(wǎng)絡輸入輸出參數(shù),建立人工神經(jīng)網(wǎng)絡訓練模型,采用商業(yè)軟件生成二維圓柱和二維翼型非結構三角形網(wǎng)格作為訓練樣本,通過訓練和學習建立起相對壁面距離和相對網(wǎng)格尺度之間的映射關系,進而實現(xiàn)不同密度的二維圓柱、不同二維翼型在不同遠場大小情況下的網(wǎng)格尺度分布控制.

1 基于背景網(wǎng)格法的尺度控制方法回顧

背景網(wǎng)格法可采用規(guī)則的笛卡爾直角結構網(wǎng)格、非規(guī)則的非結構網(wǎng)格或者自適應笛卡爾網(wǎng)格[9],各種類型的背景網(wǎng)格各有優(yōu)缺點,由于規(guī)則結構網(wǎng)格求解和插值效率較高而得到廣泛使用.

在背景網(wǎng)格上布置“點源”、“線源”或者“面源”等局部網(wǎng)格分布控制參數(shù),通過求解熱傳導方程,得到穩(wěn)態(tài)解即為全流場的尺度控制參數(shù)分布[5].該方法能夠生成分布均勻的網(wǎng)格,且能更靈活地考慮流動局部特征對網(wǎng)格分布的影響.

1.1 背景網(wǎng)格法

以二維問題為例,網(wǎng)格尺度S滿足如式(1)所示的穩(wěn)態(tài)熱傳導方程(泊松方程),以確保尺度源項的作用在物理空間內(nèi)光滑分布.

式中G為源項,定義如式(2)所示

式中下標“i”,“j”代表背景網(wǎng)格節(jié)點,N為熱源總數(shù),Ψn為第n個源的強度因子,函數(shù)In及Jn分別定義如下

式中S和f為在各源處的尺度,r為背景網(wǎng)格節(jié)點到源的距離,|ln|為線源長度.

文獻[5]還提出了考慮方向性的網(wǎng)格密度控制方法,可以通過改變源的強度函數(shù)來實現(xiàn),具體可參考文獻,本文不再贅述.

因此,采用背景網(wǎng)格法進行網(wǎng)格尺度控制的步驟為:

(1)在矩形背景網(wǎng)格上采用中心差分離散泊松方程;

(2)采用Gauss-Seidel 超松弛迭代求解離散的泊松方程,得到背景網(wǎng)格上的網(wǎng)格尺度分布;

(3)在網(wǎng)格生成過程中,根據(jù)當?shù)匚恢迷诒尘熬W(wǎng)格中插值得到當?shù)鼐W(wǎng)格尺度,用于控制網(wǎng)格生成過程.

1.2 背景網(wǎng)格法網(wǎng)格生成實例

本節(jié)采用二維圓柱、NACA0012 翼型和30P30N多段翼型作為考核算例,對背景網(wǎng)格法進行實例測試.本文算例中采用的三角形網(wǎng)格生成算法為作者發(fā)展的基于ANN (artificial neural network)的陣面推進法[22],該方法在傳統(tǒng)陣面推進法的基礎上,通過引入ANN 進行生成模式判斷和新點預測,減少了相交性判斷次數(shù),網(wǎng)格生成效率提高30%,其網(wǎng)格生成步驟可簡要概括為:

(1)將幾何邊界離散成初始陣元;

(2)從最小陣面出發(fā),自動選擇網(wǎng)格模板點,人工神經(jīng)網(wǎng)絡根據(jù)網(wǎng)格模板判斷生成模式并預測新點坐標;

(3)根據(jù)生成模式、新點坐標和局部網(wǎng)格尺度生成新網(wǎng)格單元;

(4)判斷新單元是否合適,合適則更新數(shù)據(jù)結構;

(5)回到步驟2,直至所有面變成非活躍面,整個計算域被網(wǎng)格填滿.

圖1～圖3 給出了3 個算例點源設置示意圖和生成的網(wǎng)格,圖中藍色圓點即為點源所在位置.3 個算例人工設置的點源數(shù)量分別為12,24 和44.圖中背景網(wǎng)格僅作為示意,實際背景網(wǎng)格節(jié)點數(shù)量需要根據(jù)邊界網(wǎng)格尺度進行調(diào)整,以確保對網(wǎng)格尺度分布場的分辨率和插值精度.比如對于在前后緣網(wǎng)格尺度較小的翼型,背景網(wǎng)格必須足夠密才可以有效反映出空間網(wǎng)格尺度的變化.

圖1 圓柱算例點源設置及網(wǎng)格生成情況Fig.1 Nodal source settings and corresponding triangular mesh over a 2D cylinder

圖2 NACA0012 算例點源設置及網(wǎng)格生成情況Fig.2 Nodal source settings and corresponding triangular mesh over NACA0012 airfoil

圖3 30P30N 算例點源設置及網(wǎng)格生成情況Fig.3 Nodal source settings and corresponding triangular mesh over 30P30N airfoil

在本文三個算例中,背景網(wǎng)格規(guī)模分別為51 × 51,301 × 301 和401 × 401.由于是在背景網(wǎng)格上迭代求解泊松方程,因此背景網(wǎng)格的網(wǎng)格數(shù)量直接決定了迭代求解的效率.

圖中結果顯示,在人工設置合適的點源參數(shù)后,可以在背景網(wǎng)格上得到恰當?shù)木W(wǎng)格尺度分布.

要說明的是本文暫未采用線源,實際上線源可以看作按線段排列的具有一定強度分布的點源集合,因此在處理本文的簡單問題時,只采用點源進行網(wǎng)格尺度的控制,用以說明背景網(wǎng)格法的優(yōu)缺點.

2 基于RBF 方法的尺度控制方法

徑向基函數(shù)是一種常用的插值函數(shù),常用于網(wǎng)格變形的插值[14-15,23-24],如圖4 所示即為用NACA-0012 翼型的S 型啟動來模擬魚的游動過程,采用RBF 插值方法將物面的變形插值到空間來生成動網(wǎng)格.RBF 插值方法也可以用于網(wǎng)格尺度的插值控制,在已知邊界網(wǎng)格尺度的情況下,求得插值系數(shù)矩陣,可將邊界尺度插值到整個計算域.

圖4 RBF 網(wǎng)格變形方法Fig.4 Mesh deformation controlled by RBF interpolation

2.1 基于貪婪算法的RBF 插值方法

RBF 插值方法是比較成熟的插值方法,對于一個變量場,如網(wǎng)格變形情況下的位移場,或網(wǎng)格尺度控制情況下的網(wǎng)格密度場,可以用RBF 插值公式表示為

其中,N為參考點數(shù)目,f(r) 為某待求點的函數(shù)值(網(wǎng)格尺度),r為待求點的位置矢量,ri為參考點的位置矢量,φ 為RBF 基函數(shù),‖r-ri‖ 為參考點與待求點之間的歐氏距離.wi為第i個參考點的權重系數(shù).

以網(wǎng)格尺度分布控制為例,權重系數(shù)滿足

式中Sp為各個參考點的網(wǎng)格尺度,參考點通常為邊界點,也即邊界網(wǎng)格尺度.基函數(shù)的類型包括全域型基函數(shù)和緊支型基函數(shù),具體可以參考文獻[14-15,23-24],本文不再贅述.

權重系數(shù)的求解通常需要求解系數(shù)矩陣的逆,在參考點數(shù)目較大時,會導致求解效率較低,因此可以引入貪婪算法對參考點進行精簡.

基于貪婪算法的RBF 方法的步驟可簡要歸納為:

(1)初始參考點集為空集,隨機取一邊界點作為初始參考點;

(2)根據(jù)參考點,采用RBF 插值計算其他邊界點處的網(wǎng)格尺度;

(3)由于參考點數(shù)目太少,RBF 插值得到的邊界網(wǎng)格尺度與給定值存在一定誤差,找出誤差最大的點;

(4)若誤差最大的點不為已有參考點,則將該點做為新參考點加入?yún)⒖键c序列,否則重新任選一點加入?yún)⒖键c序列;

(5)重復步驟(2)～ (4),直到最大誤差或者參考點數(shù)目滿足要求,確定最終參考點序列;

(6)根據(jù)最終參考點序列進行RBF 插值,得到空間所有位置的網(wǎng)格尺度分布.

2.2 RBF 方法網(wǎng)格生成實例

本節(jié)仍然采用二維圓柱、NACA0012 翼型和30P30N 多段翼型作為考核算例,對RBF 方法進行實例測試.算例中所采用的三角形網(wǎng)格生成算法仍為基于ANN 的陣面推進法,具體可參考文獻[22].

RBF 基函數(shù)取為Wendland’s C0,緊支半徑取為計算域范圍大小的1/4.網(wǎng)格尺度控制對尺度插值誤差的要求不高,本文以尺度的最大相對誤差不超過5%為標準,進行參考點序列精簡,表1 給出了采用貪婪算法對參考點進行精簡的結果,結果顯示對于翼型算例,精簡序列后參考點數(shù)量至少減少一半,RBF 插值的耗時也減少了近一半,在保證了插值效果的同時提高了插值效率.

表1 RBF 方法參考點的數(shù)目及插值耗時Table 1 Number of reference nodes and time consumption on interpolation

圖5～ 圖7 給出了精簡后的參考點的位置及根據(jù)精簡參考點進行尺度控制生成的非結構網(wǎng)格.結果顯示RBF 方法能夠根據(jù)邊界網(wǎng)格尺度插值得到空間的網(wǎng)格尺度,插值得到的網(wǎng)格尺度分布過渡均勻.由于初始參考點的選擇是隨機的,參考點的分布具有一定的隨機性,但是均能夠保證插值誤差滿足要求.

圖5 圓柱算例精簡后的參考點及網(wǎng)格生成情況Fig.5 Reference nodes corresponding triangular mesh over a 2D cylinder

圖6 NACA0012 算例精簡后的參考點及網(wǎng)格生成情況Fig.6 Reference nodes and corresponding triangular mesh over NACA0012 airfoil

圖7 30P30N 算例精簡后的參考點及網(wǎng)格生成情況Fig.7 Reference nodes and corresponding triangular mesh over 30P30N airfoil

3 基于人工神經(jīng)網(wǎng)絡的尺度控制方法

基于人工神經(jīng)網(wǎng)絡的深度學習具有較強的非線性擬合能力,能夠通過現(xiàn)有樣本數(shù)據(jù)的訓練識別出數(shù)據(jù)中隱含的非線性映射關系.網(wǎng)格尺度分布實際上是一個關于幾何特征、流場特征的非線性映射,如圖8 所示,幾何特征(物面曲率、狹縫、細小結構等)可以直接影響網(wǎng)格尺度分布,幾何特征也可以決定流動特征(梯度量等),從而間接決定網(wǎng)格分布.不同幾何外形的網(wǎng)格尺度分布不同,同一幾何外形在不同來流條件(流場特征)情況下,網(wǎng)格尺度分布也不同.

圖8 網(wǎng)格分布控制與幾何特征和流場特征的關系Fig.8 Relationship between mesh size control,geometry,and flow features

流場特征與幾何外形、來流條件和邊界條件相關,目前采用機器學習方法進行流場預測是一個熱點研究問題[25-26].而根據(jù)流場特征確定網(wǎng)格分布則可以根據(jù)網(wǎng)格自適應的相關準則進行,也可以采用機器學習的方法進行[27].本文初步考慮幾何外形對網(wǎng)格分布的影響,采用ANN 建立幾何特征與網(wǎng)格尺度分布之間的關系.

3.1 ANN 輸入輸出參數(shù)

Chedid 和Najiar 等[28]曾嘗試用ANN 建立起網(wǎng)格密度與幾何特征之間的關系,如圖9 所示,神經(jīng)網(wǎng)絡輸入考慮了空間點gi到邊界點的最小距離 δ1,次小距離 δ2及其對應的夾角a*和a**及到夾角對應兩邊的投影距離j1a,j1b,j2a和j2b共8 個輸入?yún)?shù),輸出該空間點處的網(wǎng)格密度,定義為一定大小區(qū)域內(nèi)節(jié)點的數(shù)量.該方法能夠較為全面地反映空間點與邊界之間的幾何關系(距離及投影距離),以及空間點對應的邊界處的幾何特征(夾角、曲率),但計算量較大不利于提高尺度控制效率,同時由于無法網(wǎng)格局部加密控制,失去了對網(wǎng)格尺度控制的靈活性,因此優(yōu)勢并未得到發(fā)揮.

圖9 文獻[28]中神經(jīng)網(wǎng)絡的輸入?yún)?shù)Fig.9 Input parameters for the artificial neural network in Ref.[28]

為提高效率,本文初步選擇最小壁面距離wdist作為輸入?yún)?shù),網(wǎng)格尺度Sp作為輸出參數(shù).為提高ANN 的泛化性,其輸入輸出參數(shù)通常需要進行歸一化操作,本文采用計算域的大小Lr_d,遠場邊界網(wǎng)格尺度Lr_f和物面邊界網(wǎng)格尺度Lr_w對輸入輸出進行分別歸一化.

同時,由于從物面到遠場,網(wǎng)格尺度變化較大,網(wǎng)格尺度相對值可能在103量級以上,為盡量縮小輸入輸出的值域范圍,提高ANN 訓練效果,本文對輸入輸出參數(shù)進行開根號.此外,在物面附近網(wǎng)格尺度變化快,需要同時考慮物面參考值和遠場參考值進行歸一化,具體輸入輸出參數(shù)形式如表2 所示.

表2 ANN 輸入輸出模型Table 2 Parameter model for artificial neural network

3.2 訓練方法及參數(shù)

本文基于Matlab 神經(jīng)網(wǎng)絡訓練工具設計全連接的人工神經(jīng)網(wǎng)絡,網(wǎng)絡含有1 個輸入層,1 個隱藏層,1 個輸出層.輸入層含有1 個神經(jīng)元,輸出層含有1 個神經(jīng)元,隱藏層的神經(jīng)元數(shù)量為10 個,激活函數(shù)采用Sigmoid 函數(shù),損失函數(shù)為均方誤差函數(shù),訓練方法采用Levenberg-Marquardt 反向傳播方法,該方法比常規(guī)的梯度下降反向傳播算法訓練效率更高,具體可以參考文獻[29].神經(jīng)網(wǎng)絡的結構如圖10所示.

圖10 基于Matlab 的人工神經(jīng)網(wǎng)絡訓練工具Fig.10 Artificial neural network training tool based on Matlab

選擇二維圓柱網(wǎng)格和NACA0012 翼型網(wǎng)格作為訓練樣本,如圖11 所示,網(wǎng)格面?zhèn)€數(shù)分別為4295,3950,每個網(wǎng)格面對應一組樣本數(shù)據(jù)點.網(wǎng)格訓練樣本按70%,15%和15%的比例隨機劃分為訓練集、測試集和驗證集,圖12 給出了在翼型三角形網(wǎng)格訓練集上的Loss 值及在驗證集上的預測精度收斂歷程.結果顯示:經(jīng)過120 次迭代后,Loss 值及預測誤差均下降到了0.001 58 左右.

圖11 網(wǎng)格分布訓練樣本網(wǎng)格Fig.11 Sample grids for ANN training

圖12 訓練Loss 值和精度收斂歷程Fig.12 Convergence of loss and accuracy on sample grids

3.3 ANN 預測結果

基于前述的ANN 訓練結果,分別預測生成了不同密度的二維圓柱網(wǎng)格、NACA0012 密網(wǎng)格、RAE2822 翼型以及三段翼型,在物面附近和遠場均取得了較好的效果,同時還能適應不同遠場大小的情況,如圖13 所示.算例中所采用的三角形網(wǎng)格生成算法仍為基于ANN 的陣面推進法,具體可參考文獻[22].

圖13 ANN 模型各向同性網(wǎng)格預測結果Fig.13 Mesh size controlled by ANN model for isotropic triangular grids

圖13 ANN 模型各向同性網(wǎng)格預測結果(續(xù))Fig.13 Mesh size controlled by ANN model for isotropic triangular grids (continued)

同時,本文還將ANN 應用于各向異性混合網(wǎng)格尺度控制,各向異性四邊形采用層推進逐層推進生成[30-32],而各向同性三角形采用基于ANN 的陣面推進生成[22].層推進的推進方向、多方向推進數(shù)量均采用ANN 預測,同時在凹角處考慮局部推進步長,避免網(wǎng)格相交,具體方法細節(jié)可以參考文獻[30].為與各向同性網(wǎng)格尺度控制相一致,本文將層推進生成的最后一層網(wǎng)格作為虛擬物面,用于計算最小壁面距離,作為ANN 控制網(wǎng)格尺度的輸入?yún)?shù).另外,層推進的網(wǎng)格尺度控制仍然采用指定物面第一層網(wǎng)格高度和增長率的方式給定.圖14 給出了NA0012翼型和30P30N 三段翼型的生成結果.由圖中結果可見,網(wǎng)格分布均勻合理,網(wǎng)格質量滿足要求,說明本文發(fā)展的ANN 方法可應用于各向異性混合網(wǎng)格的尺度控制.

圖14 ANN 模型各向異性混合網(wǎng)格預測結果Fig.14 Mesh size controlled by ANN model for anisotropic hybrid grids

4 背景網(wǎng)格法、RBF 方法與ANN 方法的比較

本文分別采用3 種方法生成了幾個典型幾何外形的非結構網(wǎng)格,其網(wǎng)格生成質量、生成效率、尺度控制靈活性、自動化程度等情況存在一定差別,本節(jié)對3 種方法進行對比分析.

從網(wǎng)格質量的角度,在網(wǎng)格尺度較小的部位,如翼型前后緣、狹縫等,背景網(wǎng)格法要求加密背景網(wǎng)格,以分辨最小網(wǎng)格尺度,而RBF 方法和ANN 方法不依賴于背景網(wǎng)格,因此在尺度較小的部位得到的尺度分布控制效果比采用矩形背景網(wǎng)格要好.

從自動化程度和靈活性角度,采用在背景網(wǎng)格上設置點源控制網(wǎng)格分布,對分布的控制最為靈活,能夠根據(jù)幾何特征和流場特征預先設置點源來改變網(wǎng)格分布.而背景網(wǎng)格法需要人工設置點源,本文算例中,在遠場和翼型附近人工設置數(shù)十個點源,每個點源人為給定強度、位置等參數(shù),需要一定的經(jīng)驗,自動化程度較低.而RBF 方法和ANN 方法只需要給定離散的邊界節(jié)點,就可以得到空間網(wǎng)格尺度分布,自動化程度相對較高,但是靈活性較低.

從網(wǎng)格生成效率的角度,背景網(wǎng)格法需要在笛卡爾網(wǎng)格上迭代求解泊松方程,在背景網(wǎng)格很密時,求解效率較低.而RBF 方法雖然需要求解矩陣的逆,但在引入貪婪算法之后,參考點數(shù)量減少,RBF 方法的效率得到提高.ANN 方法雖然需要求解最小壁面距離,但壁面距離的值并不需要十分精確,因此可以采用一些近似求法,也可以達到更高的效率.

表3 給出了3 種方法在生成網(wǎng)格過程中尺度控制所耗費的時間,由于每種方法生成的網(wǎng)格單元數(shù)存在一定差異,因此表中也給出了網(wǎng)格單元數(shù)量.

表3 3 種方法控制網(wǎng)格尺度耗時對比Table 3 Efficiency comparison of the three methods

背景網(wǎng)格法耗時基本由背景網(wǎng)格的規(guī)模及點源的數(shù)量決定.在本文三個算例中,背景網(wǎng)格規(guī)模為51 × 51,301 × 301 和401 × 401,點源的數(shù)量大致為20 個,因此在采用背景網(wǎng)格法時,翼型算例耗時明顯增加.

由表3 中數(shù)據(jù)可見,在背景網(wǎng)格數(shù)量較少時,背景網(wǎng)格法的效率較高,但是隨著背景網(wǎng)格數(shù)量和點源數(shù)量的增加,背景網(wǎng)格法效率明顯下降.相較于傳統(tǒng)背景網(wǎng)格法,RBF 方法和ANN 方法的耗時明顯減少,耗時僅為背景網(wǎng)格法耗時的1/10～ 1/5,網(wǎng)格尺度控制效率相應提高了5～ 10 倍.而且,在外形相對復雜的情況下,ANN 方法展現(xiàn)了更好的控制效率.可以預見,在三維復雜外形情況下,ANN 的控制效果和效率會更高.

5 總結與展望

本文提出了一種采用ANN 進行網(wǎng)格尺度分布控制的方法,初步確定了合理的神經(jīng)網(wǎng)絡輸入輸出參數(shù),基于Matlab 建立人工神經(jīng)網(wǎng)絡模型,采用商業(yè)軟件生成二維圓柱和NACA0012 翼型非結構三角形網(wǎng)格作為網(wǎng)格密度訓練樣本,通過訓練和學習建立起壁面距離和網(wǎng)格尺度的映射關系,實現(xiàn)了不同密度的二維圓柱和不同二維翼型的網(wǎng)格尺度預測.同時,發(fā)展了基于RBF 方法的網(wǎng)格尺度控制方法,采用貪婪算法對插值參考點序列進行精簡,將RBF 插值效率提高了一倍.

與傳統(tǒng)背景網(wǎng)格法相比,基于ANN 方法和RBF方法的網(wǎng)格尺度預測效率提高5～10 倍,有助于進一步提高網(wǎng)格生成效率.最后將基于ANN 的網(wǎng)格尺度控制方法拓展應用于各向異性混合網(wǎng)格的尺度控制,得到的網(wǎng)格質量滿足要求,證明了方法的實用性.

展望ANN 方法在尺度分布控制領域的應用前景,在以往的CFD 研究和工程實踐中,已經(jīng)積累了大量各種類型的網(wǎng)格,其中包含了網(wǎng)格尺度分布信息.對已有的網(wǎng)格數(shù)據(jù)按照幾何外形進行分類,訓練得到可以處理不同類別幾何外形的神經(jīng)網(wǎng)絡,在需要對新的幾何外形進行網(wǎng)格尺度控制時,只需選擇已經(jīng)訓練好的該類別的神經(jīng)網(wǎng)絡即可高效快速完成尺度控制.進一步還可以考慮流場特征(來流條件),對神經(jīng)網(wǎng)絡的適用范圍進行細分,使得其能夠同時考慮幾何外形和流場特征進行尺度控制.

可以進一步考慮物面幾何曲率、局部特征尺寸(窄邊、狹縫)等參數(shù),提高神經(jīng)網(wǎng)絡的泛化性,同時研究通過神經(jīng)網(wǎng)絡進行流場預測,并根據(jù)流場特征進行網(wǎng)格自適應.