基于內(nèi)外域結(jié)構(gòu)的開式轉(zhuǎn)子流場數(shù)值模擬

吳俁

摘 要：文章采用數(shù)值模擬方法對某對轉(zhuǎn)開式轉(zhuǎn)子的流場進行了研究，以驗證該方法應用于對轉(zhuǎn)開式轉(zhuǎn)子流場模擬的準確性。建立了對轉(zhuǎn)開式轉(zhuǎn)子的三維模型，并通過全三維定常數(shù)值模擬詳細分析了流場參數(shù)。通過拉力、總壓比的徑向分布、出口相對氣流角以及折轉(zhuǎn)角等參數(shù)與實驗數(shù)據(jù)的比較得知，其結(jié)果吻合較好。由此說明，采用內(nèi)外域網(wǎng)格結(jié)構(gòu)的數(shù)值模擬方法對對轉(zhuǎn)開式轉(zhuǎn)子流場的模擬結(jié)果是可靠的。

關鍵詞：對轉(zhuǎn)開式轉(zhuǎn)子；數(shù)值模擬；內(nèi)外域

1 符號表

表1

2 概述

隨著能源價格的不斷上升以及環(huán)境問題的日益突出，經(jīng)濟性和環(huán)保性成為決定新一代民用航空發(fā)動機技術發(fā)展方向的主導因素。開式轉(zhuǎn)子發(fā)動機通過采用兩轉(zhuǎn)子對轉(zhuǎn)的結(jié)構(gòu)，下游轉(zhuǎn)子能夠消除上游轉(zhuǎn)子的滑流損失，可比單轉(zhuǎn)子的推進效率提高6-8%[1-6]。上世紀80年代的對轉(zhuǎn)開式轉(zhuǎn)子發(fā)動機的飛行試驗結(jié)果表明，其燃油經(jīng)濟性比當時先進水平的渦扇發(fā)動機高30%以上[7]，相應地其污染物和溫室氣體排放也比渦扇發(fā)動機低很多。

對轉(zhuǎn)開式轉(zhuǎn)子的研究通常包括數(shù)值模擬和實驗兩種方法，文章采用內(nèi)外域網(wǎng)格結(jié)構(gòu)的數(shù)值模擬方法對其流場進行研究，通過將模擬結(jié)果和實驗結(jié)果進行對比，驗證該方法的可靠性。

3 研究對象

3.1 模型來源及介紹

研究選用文獻[8]中的轉(zhuǎn)子模型，其葉尖半徑為200mm，輪轂半徑為50mm，根部為NACA-65系列翼型，中間為雙圓弧翼型，尖部為NACA-16翼型。葉片為后掠式，后掠角為27°[8]。單轉(zhuǎn)子和對轉(zhuǎn)轉(zhuǎn)子葉型相同，單轉(zhuǎn)子為6葉片，對轉(zhuǎn)轉(zhuǎn)子前后排均為3葉片，其實驗模型如圖1所示。

圖1 單轉(zhuǎn)子和對轉(zhuǎn)轉(zhuǎn)子

3.2 模型建立

根據(jù)文獻[8]中公開的幾何數(shù)據(jù)建模，模型分為單轉(zhuǎn)子和對轉(zhuǎn)轉(zhuǎn)子兩種，如圖2所示。對轉(zhuǎn)轉(zhuǎn)子分為大小間距，小間距前后排轉(zhuǎn)子間距為80mm，大間距為160mm。

圖2 單轉(zhuǎn)子和對轉(zhuǎn)轉(zhuǎn)子幾何模型

4 數(shù)值方法

文章流場計算采用商業(yè)軟件fluent，求解雷諾平均N-S方程，選用k-？著湍流模型以及四面體非結(jié)構(gòu)網(wǎng)格，采用有限體積法進行數(shù)值求解，時間推進采用一階高精度格式。

4.1 單轉(zhuǎn)子計算域、網(wǎng)格及邊界條件

數(shù)值模擬采用單通道計算。將計算域劃分為內(nèi)外域，內(nèi)域為旋轉(zhuǎn)域，外域為靜止域。圖3給出了主要的計算域尺寸和邊界條件，進口為速度進口，出口和圓周表面為壓力出口，兩側(cè)為周期性邊界條件。內(nèi)外域交界面采用凍結(jié)轉(zhuǎn)子法。網(wǎng)格為四面體非結(jié)構(gòu)網(wǎng)格。

圖3 單轉(zhuǎn)子計算域及邊界條件

4.2 對轉(zhuǎn)轉(zhuǎn)子計算域、網(wǎng)格及邊界條件

對轉(zhuǎn)轉(zhuǎn)子的計算方法和單轉(zhuǎn)子基本相同。整體計算域分為三部分，兩個內(nèi)域和一個外域，前后排葉片在兩個不同的內(nèi)域中，如圖4所示。進口為速度進口，出口和圓周表面為壓力出口，兩側(cè)為周期性邊界條件，交界面處采用凍結(jié)轉(zhuǎn)子法。采用四面體非結(jié)構(gòu)網(wǎng)格，網(wǎng)格數(shù)分別為110w和100w。

圖4 對轉(zhuǎn)轉(zhuǎn)子計算域

5 計算結(jié)果及分析

實驗測量了單轉(zhuǎn)子的拉力、出口總壓比、相對氣流角和轉(zhuǎn)折角，對轉(zhuǎn)轉(zhuǎn)子的出口總壓比、前后排轉(zhuǎn)子間切向速度。將數(shù)值模擬得到的上述流場參數(shù)和實驗結(jié)果進行了比較。

5.1 單轉(zhuǎn)子計算結(jié)果及分析

表2列出了單轉(zhuǎn)子拉力、出口總壓比、氣流角和轉(zhuǎn)折角測量的實驗條件，分別對應工況1、2、3。

表2 單轉(zhuǎn)子各實驗工況

實驗測量的拉力指的是轉(zhuǎn)子在前進方向上提供的力，單個葉片的拉力為48.83N，數(shù)值模擬得到的單個葉片的拉力為52N，誤差為6.5%，在合理的范圍內(nèi)。

如圖5展示了出口總壓比的徑向分布。在50%葉高到90%葉高之間數(shù)值模擬和實驗結(jié)果基本一致。圖6展示了氣流轉(zhuǎn)折角和出口相對氣流角的徑向分布。無論是氣流轉(zhuǎn)折角還是出口相對氣流角，數(shù)值模擬結(jié)果和實驗結(jié)果都非常接近。

5.2 對轉(zhuǎn)轉(zhuǎn)子計算結(jié)果及分析

對轉(zhuǎn)轉(zhuǎn)子sec1、sec2、sec3的幾何示意圖如圖7。

對轉(zhuǎn)轉(zhuǎn)子的總壓比、sec1，sec2切向速度、sec3軸向速度的測量中，β0.75=36度，U0=68m/s，N1=9000rpm，后排轉(zhuǎn)子轉(zhuǎn)速則有7000rpm，8000rpm，9000rpm三種。

圖8展示了大小間距下sec3處總壓比的徑向分布。隨著后排轉(zhuǎn)子轉(zhuǎn)速的提高，出口總壓比也相應提高。小間距的數(shù)值模擬結(jié)果更接近實驗值，尤其是在葉中部分，基本相同。大間距中，數(shù)值模擬得到的總壓比則偏大，經(jīng)過分析，認為這是由于數(shù)值模擬得到的前后排轉(zhuǎn)子間的切向速度和實驗有一定偏差所致。大間距和小間距的整體趨勢都和實驗結(jié)果相吻合。

圖9給出了N2=7000rpm時，section1和section2處的切向速度徑向分布。在r/R=0.5處，切向速度均為15m/s左右。差別偏大的地方在葉尖部分，分析其原因是葉尖附近有網(wǎng)格的轉(zhuǎn)靜交界面，給模擬帶來誤差。

圖9 sec.1和sec.2切向速度徑向分布

6 結(jié)束語

文章建立了單轉(zhuǎn)子和對轉(zhuǎn)轉(zhuǎn)子的計算模型，對其流場進行了數(shù)值模擬和分析，結(jié)論如下：采用基于內(nèi)外域網(wǎng)格結(jié)構(gòu)的雷諾平均N-S方程定常數(shù)值模擬能夠比較準確地模擬單轉(zhuǎn)子和對轉(zhuǎn)轉(zhuǎn)子的流場。

參考文獻

[1]孫曉峰，胡宗安，周盛， 反旋槳扇非定常負荷噪聲的研究[J].航空動力學報，1988，3（3）：227-230.

[2]F. Bruce Metzger. A Review of Propeller Noise Prediction Methodology 1919-1994[S]. NASA CR-198156， 1995.

[3]Donald B. Hanson. Compressible Helicoidal Surf ace Theory for Propeller Aerodynamics and Noise[S]. AIAA JOURNAL， 1983，21（6）：881-889.

[4]Donald B. Hanson. Near-Field Frequency-Domain Theory for Propeller Noise[J]. AIAA JOURNAL， 1985，23（4）：499-504.

[5]Donald B. Hanson. Direct Frequency Domain Calculation of Open Rotor Noise[J]. AIAA JOURNAL， 1992， 30（9）：2334-2337.

[6]Donald B. Hanson. Noise of Counter Rotation Propellers[J].AIAA-84-2305， 1984.

[7]What Ever Happened to Propfans[J].Flight Global， 2007.

[8]Fujii， S. Nishiwaki， H.， and Takeda， K.， Aeroacoustics of an Advanced Propeller Design Under Take-Off and Landing Conditions[J].Journal of Aircraft，1986，23（2）：136-141.

[9]Fujii， S. Noise and Performance of a Counter-Rotation Propeller[J].Journal of Aircraft，1986，23（11）：719-724.