李昂 沈偉 張宇坤

摘 要：針對各類網(wǎng)格劃分軟件生成的網(wǎng)格特征差異明顯、格式不一、無法直接進(jìn)行拼接處理的問題，本文以CFX-Pre提供的非相容計算域交界面網(wǎng)格處理技術(shù)為基礎(chǔ)，通過對網(wǎng)格進(jìn)行物理變換及相關(guān)技術(shù)處理，在保證計算精度并精簡計算量的前提下生成格式統(tǒng)一的一體化整機(jī)流場計算網(wǎng)格。結(jié)果表明通過合理運用CFX-Pre網(wǎng)格拼接技術(shù)，可以實現(xiàn)不同格式異構(gòu)網(wǎng)格的拼接處理，極大地簡化了復(fù)雜物理模型的網(wǎng)格建立過程，對流場仿真具有重要意義。

關(guān)鍵詞：CFD 異構(gòu)流場網(wǎng)格 網(wǎng)格拼接

中圖分類號：TP391.9 文獻(xiàn)標(biāo)識碼：A 文章編號：1672-3791（2018）08（a）-0082-05

網(wǎng)格劃分是CFD數(shù)值計算的一項重要環(huán)節(jié)。進(jìn)行全流場網(wǎng)格劃分，首先多基于流場物理結(jié)構(gòu)特點，將目標(biāo)流場劃分為若干計算區(qū)域。利用不同格式網(wǎng)格結(jié)構(gòu)特點，通過相應(yīng)軟件分別進(jìn)行各計算域網(wǎng)格劃分，并通過技術(shù)手段進(jìn)行非一致或異構(gòu)網(wǎng)格拼接。高效的拼接方法能夠進(jìn)一步提高CFD在工程領(lǐng)域的應(yīng)用效率。

國外建立了多種針對CFD應(yīng)用的網(wǎng)格劃分軟件，如Gridgen、ICEM CFD、Gambit、CFD-Geom等，基本解決了CFD應(yīng)用中的流場網(wǎng)格劃分問題。但各類網(wǎng)格劃分軟件生成的網(wǎng)格格式不同，存儲形式不一致，且各軟件默認(rèn)基礎(chǔ)坐標(biāo)系及尺寸單位等均有不同，必須通過適當(dāng)?shù)木W(wǎng)格拼接技術(shù)進(jìn)行統(tǒng)一處理。目前關(guān)于不同格式或異構(gòu)網(wǎng)格拼接方法的研究依然有限，如何通過技術(shù)處理實現(xiàn)不同格式網(wǎng)格拼接，是CFD數(shù)值計算亟需解決的問題。

CFX-Pre支持輸入各種類型、各種來源的網(wǎng)格，具有良好的網(wǎng)格兼容性。在此種背景下，利用CFX-Pre前處理功能實現(xiàn)網(wǎng)格拼接是一種可行手段。

本文擬針對某型發(fā)動機(jī)流場計算各類網(wǎng)格格式及結(jié)構(gòu)特征不一致的問題，運用CFX-Pre進(jìn)行技術(shù)處理，生成高質(zhì)量一體化全流場網(wǎng)格，為前處理軟件流場網(wǎng)格拼接提供借鑒。

1 導(dǎo)入網(wǎng)格

發(fā)動機(jī)各核心氣路部件特征結(jié)構(gòu)存在較大差異，通常運用不同網(wǎng)格劃分軟件對相關(guān)部件進(jìn)行網(wǎng)格劃分，進(jìn)而生成不同格式、不同類別的各氣路部件網(wǎng)格。這種劃分方法往往導(dǎo)致各部件網(wǎng)格類別差異較大，無法直接進(jìn)行統(tǒng)一操作。在這種情況下，需通過技術(shù)手段將各個網(wǎng)格文件導(dǎo)入至同一前處理軟件并進(jìn)行相應(yīng)技術(shù)處理，以滿足后續(xù)計算需要。

運用各類網(wǎng)格劃分軟件創(chuàng)建的三角形/四面體網(wǎng)格（非結(jié)構(gòu)化網(wǎng)格）或四邊形/六面體網(wǎng)格（結(jié)構(gòu)化網(wǎng)格），均可通過網(wǎng)格輸入向?qū)?dǎo)入至CFX-Pre前處理軟件中，并進(jìn)行相應(yīng)技術(shù)處理。

CFX-Pre支持以下格式網(wǎng)格文件：

Anysis Meshing（*.cmdb/*.dsdb）；CFX Mesh（*.gtm/*.cfx）；CFX-Solver Input（*.def/*.res/*.cfg

/*.trn/*.bak）；ICEM CFD（*.cfx/*.cfx5/*.msh）；Anysis（*.cdb/*.inp）；Fluent（*.cas/*.msh）；OGNS（*.cgns/*.cgn）；CFX-TASCflow（*.grd）；CFX-4（*.geo）；CFX-BladeGenPlus（*.bgd）等。

針對不同格式網(wǎng)格，CFX-Pre共提供以下兩種輸入向?qū)А?/p>

（1）對于常規(guī)葉輪機(jī)械網(wǎng)格，CFX-Pre在導(dǎo)入過程中可以自動識別，進(jìn)而自動生成網(wǎng)格邊界條件，如圖1所示。

（2）對于Gambit及ICEM等軟件生成的各類異構(gòu)網(wǎng)格，CFX-Pre可以實現(xiàn)直接讀取，便于后續(xù)網(wǎng)格處理，如圖2所示。

通過以上兩種網(wǎng)格輸入向?qū)?，將通過各類網(wǎng)格劃分軟件得到的網(wǎng)格文件導(dǎo)入至CFX-Pre前處理軟件，以便進(jìn)行后續(xù)技術(shù)處理，最終生成高質(zhì)量一體化網(wǎng)格。

2 網(wǎng)格變換

由于各軟件在網(wǎng)格劃分過程中默認(rèn)的基礎(chǔ)坐標(biāo)系及尺寸單位、劃分方法等均有不同，統(tǒng)一導(dǎo)入的各部件網(wǎng)格往往具有錯誤的物理關(guān)系，需對網(wǎng)格進(jìn)行相應(yīng)變換處理，使坐標(biāo)系及尺寸鏈保持一致。

CFX-Pre允許用戶通過不同手段修改導(dǎo)入的網(wǎng)格文件，包括網(wǎng)格平移變換、旋轉(zhuǎn)變換、映射變換及放縮變換等。

2.1 平移變換

通過指定節(jié)點的笛卡爾坐標(biāo)偏移量，或定義平移前后三維坐標(biāo)，均可以進(jìn)行網(wǎng)格平移操作。同時，CFX-Pre允許用戶進(jìn)行多倍復(fù)制，并保留原始網(wǎng)格模型，如圖3所示。

運用網(wǎng)格平移功能，并定義計算域交界面相應(yīng)計算方法，可以將不同計算域進(jìn)行拼接，進(jìn)而生成完整計算區(qū)域，同時需要注意的是，坐標(biāo)偏移量可以為負(fù)值，如圖4所示。

為縮減計算量，進(jìn)而達(dá)到簡化求解過程的目的，在CFX-Pre中可以將多重計算域合并為單一計算域，便于后續(xù)邊界條件的設(shè)定。在邊界條件設(shè)定過程中，針對每個計算域所包含的表面，必須進(jìn)行相關(guān)定義。當(dāng)計算域特征相似或相同時，計算域表面（如壁面、壓力入口、壓力出口等）的設(shè)定將花費較多時間及精力。這種情況下，進(jìn)行計算域合并便可體現(xiàn)其優(yōu)越性。

2.2 旋轉(zhuǎn)變換

由圖5可知，通過定義旋轉(zhuǎn)主軸、旋轉(zhuǎn)方向及相應(yīng)旋轉(zhuǎn)角度可以對目標(biāo)網(wǎng)格進(jìn)行旋轉(zhuǎn)變換。對于旋轉(zhuǎn)主軸，CFX-Pre提供直觀定義及矢量定義兩種方法。

圖6為對某型渦槳發(fā)動機(jī)進(jìn)氣道網(wǎng)格進(jìn)行旋轉(zhuǎn)變換，實現(xiàn)與進(jìn)口導(dǎo)流葉片的網(wǎng)格拼接，生成一體化計算區(qū)域。

2.3 映射變換

CFX-Pre同樣提供針對目標(biāo)網(wǎng)格的映射變換。這里提供多種映射平面定義方法：可以直接選定坐標(biāo)平面也可定義3個坐標(biāo)點確定映射平面，進(jìn)而依據(jù)需求對目標(biāo)網(wǎng)格進(jìn)行映射，如圖7所示。

需注意的是，映射變換通常以坐標(biāo)平面為映射平面，故對目標(biāo)進(jìn)行映射變換后，往往還需對目標(biāo)進(jìn)行平移變換等技術(shù)處理，同時針對計算域網(wǎng)格交界面進(jìn)行合理拼接，從而生成流場計算所需的一體化網(wǎng)格。

2.4 放縮變換

定義放縮因子，及放縮變換參考坐標(biāo)，可對網(wǎng)格進(jìn)行放縮變換，進(jìn)而修正不同軟件所劃分網(wǎng)格存在的尺寸差異。需注意的是，對網(wǎng)格進(jìn)行放縮后，所劃分的網(wǎng)格單元會進(jìn)行同樣的放縮變換，如圖9、圖10所示。

某型渦槳發(fā)動機(jī)各氣路部件網(wǎng)格劃分過程中，為滿足計算精度需要，劃分高質(zhì)量計算域網(wǎng)格，針對燃燒室、螺旋槳計算域及柱體包裹區(qū)域均采用非結(jié)構(gòu)化網(wǎng)格劃分，非結(jié)構(gòu)化網(wǎng)格由于不具備規(guī)則拓?fù)浣Y(jié)構(gòu)、網(wǎng)格節(jié)點分布隨意等特點，可以較好地適應(yīng)具有復(fù)雜結(jié)構(gòu)外形的計算區(qū)域。

進(jìn)氣道、壓氣機(jī)及渦輪等典型葉輪式部件采用結(jié)構(gòu)化網(wǎng)格劃分。結(jié)構(gòu)化網(wǎng)格往往具有更好的網(wǎng)格質(zhì)量，具備網(wǎng)格生成速度快、數(shù)據(jù)結(jié)構(gòu)簡單等特點，可以在保證計算精度的前提下縮減計算量，進(jìn)而提高計算效率，如圖11所示。

各氣路部件結(jié)構(gòu)化及非結(jié)構(gòu)化網(wǎng)格統(tǒng)一導(dǎo)入至CFX-Pre中，并依據(jù)制定的整機(jī)尺寸鏈對各部件網(wǎng)格進(jìn)行相應(yīng)旋轉(zhuǎn)、平移、映射及放縮變換，通過網(wǎng)格拼接生成發(fā)動機(jī)整機(jī)計算域一體化網(wǎng)格，滿足后續(xù)整機(jī)流場數(shù)值計算分析需要，如圖12所示。

3 網(wǎng)格拼接原理

僅對網(wǎng)格進(jìn)行物理變換無法實現(xiàn)交界面數(shù)據(jù)傳遞，還需進(jìn)一步對各計算域網(wǎng)格交界面進(jìn)行技術(shù)處理。

當(dāng)拼接計算域具有非一致網(wǎng)格邊界，即兩個相鄰計算區(qū)域網(wǎng)格節(jié)點位置并不相同時，CFX-Pre首先計算組成邊界的界面區(qū)域交叉點，并依據(jù)交叉情況完成流過網(wǎng)格交界面的數(shù)值計算。由圖13可知，界面區(qū)域由A-B、B-C、D-E、E-F組成，并在交界面處產(chǎn)生面a-d、d-b、b-e、e-c。當(dāng)計算由交界面流入單元Ⅳ時，面D-E即被忽略，而面d-b和b-e被使用，并分別將信息由單元Ⅰ和單元Ⅲ帶入至單元Ⅳ中。

當(dāng)計算域交界面相容性較差，且計算域具有周向運動特性時，CFX-Pre提供另一種拼接方法Mixing-Plane，依據(jù)流場流動方向定義上游交界面及下游交界面，進(jìn)而實現(xiàn)計算域交界面網(wǎng)格拼接，如圖14所示。

在Mixing-Plane拼接方法中，將流體區(qū)域近似為穩(wěn)態(tài)問題進(jìn)行求解，將相鄰區(qū)域流場數(shù)據(jù)在動、靜區(qū)域交界面沿周向進(jìn)行加權(quán)平均，并作為邊界條件進(jìn)行傳遞。同時Mixing-Plane允許設(shè)置周期性邊界條件，可以極大地簡化旋轉(zhuǎn)機(jī)械流場數(shù)值模擬的計算量。

通過上述技術(shù)處理，對發(fā)動機(jī)各氣路部件網(wǎng)格文件進(jìn)行相應(yīng)處理，進(jìn)而生成可以用于流場計算的高質(zhì)量統(tǒng)一前處理網(wǎng)格文件為后續(xù)流場仿真分析做出鋪墊。

4 結(jié)語

網(wǎng)格劃分是航空發(fā)動機(jī)流場計算前處理的重要一環(huán)。發(fā)動機(jī)各核心氣路部件結(jié)構(gòu)特征區(qū)別較大，通過相應(yīng)網(wǎng)格劃分軟件生成的網(wǎng)格往往存在明顯差異。本文通過運用CFX-Pre提供的不相容計算區(qū)域網(wǎng)格交界面處理技術(shù)，針對不同格式及不同類型網(wǎng)格不易進(jìn)行拼接計算的問題提供了一種可行的通用處理方法，通過對網(wǎng)格文件進(jìn)行物理變換，并對計算域交界面進(jìn)行技術(shù)處理，生成數(shù)值計算所需的高質(zhì)量一體化網(wǎng)格文件，進(jìn)而滿足后續(xù)計算需求，同時為Fluent、Anysis等商用軟件網(wǎng)格拼接提供一種可借鑒的技術(shù)手段。

參考文獻(xiàn)

[1] 洪俊武，梁孝平，王光學(xué)，等.多重拼接網(wǎng)格技術(shù)應(yīng)用研究[J].空氣動力學(xué)報，2007，25（1）：45-49.

[2] 侯樹強(qiáng)，王燦星，林建忠.葉輪機(jī)械內(nèi)部流場數(shù)值模擬研究綜述[J].流體機(jī)械，2005，33（5）：30-34.

[3] 詹鴻飛，黃帥，楊欣毅.多級渦輪工作流場的數(shù)值仿真研究[J].系統(tǒng)仿真技術(shù)，2016，12（2）：118-122.