基于CFX-Pre的網(wǎng)格拼接技術(shù)研究

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

2、1（2018）08（a）-0082-05網(wǎng)格劃分是cfd數(shù)值計算的一項重要環(huán)節(jié)。進行全流場網(wǎng)格劃分，首先多基于流場物理結(jié)構(gòu)特點，將目標流場劃分為若干計算區(qū)域。利用不同格式網(wǎng)格結(jié)構(gòu)特點，通過相應(yīng)軟件分別進行各計算域網(wǎng)格劃分，并通過技術(shù)手段進行非一致或異構(gòu)網(wǎng)格拼接。高效的拼接方法能夠進一步提高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)格格式不同，存儲形式不一致，且各軟件默認基礎(chǔ)坐標系及尺寸單位等均有不同，必須通過適當?shù)木W(wǎng)格拼接技術(shù)進

3、行統(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ā)動機流場計算各類網(wǎng)格格式及結(jié)構(gòu)特征不一致的問題，運用cfx-pre進行技術(shù)處理，生成高質(zhì)量一體化全流場網(wǎng)格，為前處理軟件流場網(wǎng)格拼接提供借鑒。1 導入網(wǎng)格發(fā)動機各核心氣路部件特征結(jié)構(gòu)存在較大差異，通常運用不同網(wǎng)格劃分軟件對相關(guān)部件進行網(wǎng)格劃分，進而生成不同格式、不同類別的各氣路部件網(wǎng)格。這種劃分方法往往導致各部

4、件網(wǎng)格類別差異較大，無法直接進行統(tǒng)一操作。在這種情況下，需通過技術(shù)手段將各個網(wǎng)格文件導入至同一前處理軟件并進行相應(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ū)胫羉fx-pre前處理軟件中，并進行相應(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/*.m

5、sh）；anysis（*.cdb/*.inp）；fluent（*.cas/*.msh）；ogns（*.cgns/*.cgn）；cfx-tascflow（*.grd）；cfx-4（*.geo）；cfx-bladegenplus（*.bgd）等。針對不同格式網(wǎng)格，cfx-pre共提供以下兩種輸入向?qū)?。?）對于常規(guī)葉輪機械網(wǎng)格，cfx-pre在導入過程中可以自動識別，進而自動生成網(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)格文件導入至cfx-

6、pre前處理軟件，以便進行后續(xù)技術(shù)處理，最終生成高質(zhì)量一體化網(wǎng)格。2 網(wǎng)格變換由于各軟件在網(wǎng)格劃分過程中默認的基礎(chǔ)坐標系及尺寸單位、劃分方法等均有不同，統(tǒng)一導入的各部件網(wǎng)格往往具有錯誤的物理關(guān)系，需對網(wǎng)格進行相應(yīng)變換處理，使坐標系及尺寸鏈保持一致。cfx-pre允許用戶通過不同手段修改導入的網(wǎng)格文件，包括網(wǎng)格平移變換、旋轉(zhuǎn)變換、映射變換及放縮變換等。2.1 平移變換通過指定節(jié)點的笛卡爾坐標偏移量，或定義平移前后三維坐標，均可以進行網(wǎng)格平移操作。同時，cfx-pre允許用戶進行多倍復制，并保留原始網(wǎng)格模型，如圖3所示。運用網(wǎng)格平移功能，并定義計算域交界面相應(yīng)計算方法，可以將不同計算域進行拼接，進

7、而生成完整計算區(qū)域，同時需要注意的是，坐標偏移量可以為負值，如圖4所示。為縮減計算量，進而達到簡化求解過程的目的，在cfx-pre中可以將多重計算域合并為單一計算域，便于后續(xù)邊界條件的設(shè)定。在邊界條件設(shè)定過程中，針對每個計算域所包含的表面，必須進行相關(guān)定義。當計算域特征相似或相同時，計算域表面（如壁面、壓力入口、壓力出口等）的設(shè)定將花費較多時間及精力。這種情況下，進行計算域合并便可體現(xiàn)其優(yōu)越性。2.2 旋轉(zhuǎn)變換由圖5可知，通過定義旋轉(zhuǎn)主軸、旋轉(zhuǎn)方向及相應(yīng)旋轉(zhuǎn)角度可以對目標網(wǎng)格進行旋轉(zhuǎn)變換。對于旋轉(zhuǎn)主軸，cfx-pre提供直觀定義及矢量定義兩種方法。圖6為對某型渦槳發(fā)動機進氣道網(wǎng)格進行旋轉(zhuǎn)變換，

8、實現(xiàn)與進口導流葉片的網(wǎng)格拼接，生成一體化計算區(qū)域。2.3 映射變換cfx-pre同樣提供針對目標網(wǎng)格的映射變換。這里提供多種映射平面定義方法：可以直接選定坐標平面也可定義3個坐標點確定映射平面，進而依據(jù)需求對目標網(wǎng)格進行映射，如圖7所示。需注意的是，映射變換通常以坐標平面為映射平面，故對目標進行映射變換后，往往還需對目標進行平移變換等技術(shù)處理，同時針對計算域網(wǎng)格交界面進行合理拼接，從而生成流場計算所需的一體化網(wǎng)格。2.4 放縮變換定義放縮因子，及放縮變換參考坐標，可對網(wǎng)格進行放縮變換，進而修正不同軟件所劃分網(wǎng)格存在的尺寸差異。需注意的是，對網(wǎng)格進行放縮后，所劃分的網(wǎng)格單元會進行同樣的放縮變換，

9、如圖9、圖10所示。某型渦槳發(fā)動機各氣路部件網(wǎng)格劃分過程中，為滿足計算精度需要，劃分高質(zhì)量計算域網(wǎng)格，針對燃燒室、螺旋槳計算域及柱體包裹區(qū)域均采用非結(jié)構(gòu)化網(wǎng)格劃分，非結(jié)構(gòu)化網(wǎng)格由于不具備規(guī)則拓撲結(jié)構(gòu)、網(wǎng)格節(jié)點分布隨意等特點，可以較好地適應(yīng)具有復雜結(jié)構(gòu)外形的計算區(qū)域。進氣道、壓氣機及渦輪等典型葉輪式部件采用結(jié)構(gòu)化網(wǎng)格劃分。結(jié)構(gòu)化網(wǎng)格往往具有更好的網(wǎng)格質(zhì)量，具備網(wǎng)格生成速度快、數(shù)據(jù)結(jié)構(gòu)簡單等特點，可以在保證計算精度的前提下縮減計算量，進而提高計算效率，如圖11所示。各氣路部件結(jié)構(gòu)化及非結(jié)構(gòu)化網(wǎng)格統(tǒng)一導入至cfx-pre中，并依據(jù)制定的整機尺寸鏈對各部件網(wǎng)格進行相應(yīng)旋轉(zhuǎn)、平移、映射及放縮變換，通過網(wǎng)

10、格拼接生成發(fā)動機整機計算域一體化網(wǎng)格，滿足后續(xù)整機流場數(shù)值計算分析需要，如圖12所示。3 網(wǎng)格拼接原理僅對網(wǎng)格進行物理變換無法實現(xiàn)交界面數(shù)據(jù)傳遞，還需進一步對各計算域網(wǎng)格交界面進行技術(shù)處理。當拼接計算域具有非一致網(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-e即被忽略，而面d-b和b-e被使用，并分別將信息由單元和單元帶入至單元中。當計算域交界面相容性較差，

11、且計算域具有周向運動特性時，cfx-pre提供另一種拼接方法mixing-plane，依據(jù)流場流動方向定義上游交界面及下游交界面，進而實現(xiàn)計算域交界面網(wǎng)格拼接，如圖14所示。在mixing-plane拼接方法中，將流體區(qū)域近似為穩(wěn)態(tài)問題進行求解，將相鄰區(qū)域流場數(shù)據(jù)在動、靜區(qū)域交界面沿周向進行加權(quán)平均，并作為邊界條件進行傳遞。同時mixing-plane允許設(shè)置周期性邊界條件，可以極大地簡化旋轉(zhuǎn)機械流場數(shù)值模擬的計算量。通過上述技術(shù)處理，對發(fā)動機各氣路部件網(wǎng)格文件進行相應(yīng)處理，進而生成可以用于流場計算的高質(zhì)量統(tǒng)一前處理網(wǎng)格文件為后續(xù)流場仿真分析做出鋪墊。4 結(jié)語網(wǎng)格劃分是航空發(fā)動機流場計算前處理

12、的重要一環(huán)。發(fā)動機各核心氣路部件結(jié)構(gòu)特征區(qū)別較大，通過相應(yīng)網(wǎng)格劃分軟件生成的網(wǎng)格往往存在明顯差異。本文通過運用cfx-pre提供的不相容計算區(qū)域網(wǎng)格交界面處理技術(shù)，針對不同格式及不同類型網(wǎng)格不易進行拼接計算的問題提供了一種可行的通用處理方法，通過對網(wǎng)格文件進行物理變換，并對計算域交界面進行技術(shù)處理，生成數(shù)值計算所需的高質(zhì)量一體化網(wǎng)格文件，進而滿足后續(xù)計算需求，同時為fluent、anysis等商用軟件網(wǎng)格拼接提供一種可借鑒的技術(shù)手段。參考文獻1 洪俊武，梁孝平，王光學，等.多重拼接網(wǎng)格技術(shù)應(yīng)用研究j.空氣動力學報，2007，25（1）：45-49.2 侯樹強，王燦星，林建忠.葉輪機械內(nèi)部流場數(shù)值模擬研究綜述j.流體機械，2005，33（5）：30-34.3 詹鴻飛，黃