第十六章 核心機 理論計算

　　結(jié)構(gòu)的選擇只是核心機設(shè)計的第一步，陳東風(fēng)要根據(jù)他們開始設(shè)定的航模發(fā)動機的基本參數(shù)來把這些結(jié)構(gòu)的數(shù)學(xué)模型給搭建出來。這樣他才能徹底了解設(shè)計的核心機的理論參數(shù)和實際的要求是否一致。而且理論模型的建立，也是可以大大的減輕日后核心機試車的時候的調(diào)試工作。雖然他們的航模能不能從圖紙上變成實物還不一定，但夢想總要有的，萬一實現(xiàn)了呢？

　　陳東風(fēng)由于大腦開發(fā)過度，對各條公式定理的掌握也比較熟悉，他的數(shù)學(xué)模型雖然不會有很大的開創(chuàng)性，但是他的數(shù)學(xué)模型有著較高的置信度。

　　陳東風(fēng)首先考慮到壓氣機是離心式壓氣機，這種壓氣機一級的增壓比就可以達到5:1，再考慮到航模的速度，就可以得到氣體的流量。如果選擇汽油作為燃料那油氣比理想情況下要達到14：1，就可以粗略的計算出燃燒室的總體積。再由熱力學(xué)的知識把燃燒后氣體的焓轉(zhuǎn)變?yōu)閯幽?，就可以得出噴管處的推力了?；具^程就是這樣，但是其中有很多衰減的因素、效率、損失等需要考慮。

　　陳東風(fēng)主要的計算難點在燃燒室氣動熱力性能數(shù)值分析上。這首先需要建立基本守恒方程組從流體力學(xué)、熱力學(xué)、傳熱傳質(zhì)學(xué)和燃燒學(xué)等基本原理出發(fā)，建立質(zhì)量、動量、能量和組分等守恒方程。

　　再就是要方程確定定解條件，定解條件包括初始條件和邊界條件。邊界條件是指燃燒室進出口，軸線（或?qū)ΨQ面），壁面，自由表面處條件。對于非定常流還要給定初始條件。

　　還要為燃燒室的選擇物理模型。湍流兩相燃燒的基本方程組中有許多項是未知的，因此方程組是不封閉的。為了使方程組封閉必須構(gòu)造或選擇相應(yīng)的物理模型即湍流兩相燃燒模型。

　　最后是為建立離散化方程。湍流兩相燃燒流動的基本方程具有數(shù)目多，相互耦合和非線性的特點，這就決定了在一般情況下只能用數(shù)值方法求解。偏微分方程組的離散化是使用計算機求解的前提。目前有許多離散化方法，如有限差分法等；有中心差分、高階差分等格式。

　　陳東風(fēng)選擇的是有限差分法，而對于湍流兩相燃燒，常用的是有限差分中的控制容積法，一階或高階迎風(fēng)差分格式和隱式方案。求解區(qū)域離散化，又稱計算網(wǎng)格的剖分是微分方程離散的基礎(chǔ)。網(wǎng)格劃分關(guān)系到離散方程的形式、收斂性、經(jīng)濟性及準確性。對于簡單幾何形狀可采用直角坐標系或圓柱坐標系。對于復(fù)雜幾何形狀需要采用曲線坐標系。對于單項流動的離散化方程組，陳東風(fēng)采用的是一種進積分算法，用于求解拋物型問題（如邊界層、射流、管流等）。

　　好在陳東風(fēng)計算過硬，硬是靠一把筆和一疊草稿紙計算了出來。

　　燃燒室的氣動熱力性能數(shù)值分析室基礎(chǔ)，還要在這之上，對葉輪機包括壓氣機和渦輪進行理論計算。

　　葉輪機的理論計算需要用到基本焓方程：Δh*=Lu =ω(r2C2u -r1C1u)。

　　壓氣機設(shè)計轉(zhuǎn)速確定后，可以先選定動葉幾何出口角β2，再根據(jù)加功量（渦輪輸出功率）計算流量，最佳的β2值需經(jīng)過較詳細的計算才能確定，要從流動效率高和易于制造兩個方面考慮來選擇合適的值。

　　確定轉(zhuǎn)速、轉(zhuǎn)子葉片幾何出口角β2和加功量后，可以求出壓氣機的流量和氣流流入擴壓器的速度。

　　由于氣流離開轉(zhuǎn)子不是完全以葉型的幾何出氣角流出，而總是有一點“滑移”，造成實際的C2u值小于理想值（氣流以葉型的幾何出氣角流出轉(zhuǎn)子時的C2u）。通過引入滑移因子σ，可以計算實際的C2u。

　　至于渦輪，由于其性能受渦輪葉片幾何形狀細節(jié)的影響小，這一點同壓氣機不一樣，因為氣流流經(jīng)渦輪時是加速流動，邊界層分離的可能性小，但是，應(yīng)爭取獲得盡可能高的效率。

　　渦輪設(shè)計應(yīng)使渦輪轉(zhuǎn)子出口處氣流的Cu接近與零，而壓氣機轉(zhuǎn)子進口氣流的Cu也為零，由能量守恒方程：mCLK=mTLT（mC為壓氣機流量，mT為渦輪流量）可得：mT（rCu）渦輪轉(zhuǎn)子進口＝mC（rCu）壓氣機轉(zhuǎn)子出口。等式右邊由壓氣機設(shè)計參數(shù)確定。

　　由渦輪轉(zhuǎn)子進口速度三角形得：tanα3，rel =(C3u-u3)/C3r，C3r由連續(xù)方程確定，轉(zhuǎn)子進口設(shè)計攻角為零時，轉(zhuǎn)子進口葉型幾何角β3=α3rel 。轉(zhuǎn)子出口C4u =0，由轉(zhuǎn)子出口速度三角形得：tanα4，rel =- u4/C4r，C4r由連續(xù)方程確定，α4，rel為負值。

　　由于轉(zhuǎn)子出口氣流角和葉型幾何角之間存在著落后角，目前尚無有關(guān)該落后角的準確估算方法，若轉(zhuǎn)子葉片數(shù)目較多，并且氣流在轉(zhuǎn)子通道中的轉(zhuǎn)角也比較小，可以認為落后角為零度，這樣轉(zhuǎn)子出口葉型幾何角β4即為轉(zhuǎn)子出口相對氣流角α4，rel。

　　渦輪轉(zhuǎn)子的喉道面積是另一個影響流量的重要參數(shù)，應(yīng)該使得氣體在喉道處的速度與在渦輪出口處的速度相同，即滿足下列關(guān)系：N葉片數(shù)×喉道寬度=2πr4cos(α4，rel)。

　　這樣就要求轉(zhuǎn)子葉片在出口附近可能曲率比較大，準確的喉道面積比葉型角β4更為重要，但是兩者都應(yīng)盡可能地滿足各自的設(shè)計要求。

　　渦輪轉(zhuǎn)子葉片的數(shù)目可以用與壓氣機相同的方法在中徑處（rmid）估算，由于前緣附近氣流容易分離，建議按比估算出的葉片數(shù)目多25%進行設(shè)計。如果采用大小葉片設(shè)計，則小葉片放置在渦輪轉(zhuǎn)子進口部分。

　　轉(zhuǎn)子葉型可以按與壓氣機相同的方法進行設(shè)計，即沿半徑負荷ΔPb不變?nèi)~型或圓弧葉型，但是，所要求的喉道面積必須保證，喉道寬度可以從計算出。

　　海量的計算后，陳東風(fēng)終于把核心機的理論模型建立起來，包括壓氣機和渦輪葉片的彎曲角度設(shè)計，燃燒室燃燒熱值分析，尾噴口動量分析等。

　　然而這就完成了嗎？不，遠不夠...