袁 武，閻 超，楊 威

（1.中國(guó)運(yùn)載火箭技術(shù)研究院 空間物理重點(diǎn)實(shí)驗(yàn)室，北京 100076；2.北京航空航天大學(xué) 國(guó)家計(jì)算流體力學(xué)實(shí)驗(yàn)室，北京 100191）

0 引 言

重疊網(wǎng)格［1］方法在初始挖洞結(jié)束后，若挖洞曲面距離物面很近或是直接將物面默認(rèn)為挖洞曲面，則可能導(dǎo)致洞邊界貼近物面，使網(wǎng)格重疊區(qū)域龐大，同時(shí)插值區(qū)流場(chǎng)梯度很大，影響流場(chǎng)求解的效率，且降低解的精度。

人們希望能將插值區(qū)從物面附近移開(kāi)，以避免插值對(duì)流場(chǎng)計(jì)算的不利影響，因此，出現(xiàn)了一些在挖洞結(jié)束后能對(duì)洞面進(jìn)行優(yōu)化的技術(shù)，如陣面推進(jìn)技術(shù)［2］和割補(bǔ)法［3］等。其中，割補(bǔ)法建立在迭代優(yōu)化思想的基礎(chǔ)上，洞邊界點(diǎn)離散的做推進(jìn)運(yùn)動(dòng)，經(jīng)切割和填補(bǔ)兩個(gè)階段改善洞面的質(zhì)量，由于無(wú)需人工控制洞邊界點(diǎn)的推進(jìn)距離，因此自動(dòng)化程度較高，是目前最有代表性的洞面優(yōu)化技術(shù)之一。

本文對(duì)傳統(tǒng)的割補(bǔ)法技術(shù)進(jìn)行了研究，指出其在處理復(fù)雜網(wǎng)格時(shí)可靠性較差，并針對(duì)割補(bǔ)法的缺陷提出了改進(jìn)的方法，通過(guò)數(shù)值算例對(duì)新方法進(jìn)行了驗(yàn)證。

1 割補(bǔ)法簡(jiǎn)介

割補(bǔ)法將挖洞結(jié)束后形成的網(wǎng)格洞邊界經(jīng)切割和填補(bǔ)兩步，逐漸推離物面，最終在遠(yuǎn)離物面的地方形成網(wǎng)格重疊區(qū)域。在切割階段，由洞邊界點(diǎn)構(gòu)成的陣面逐步遠(yuǎn)離物面向外擴(kuò)張，直至再也沒(méi)有網(wǎng)格單元重疊；填補(bǔ)階段，洞邊界面逐步收縮，直至出現(xiàn)洞外網(wǎng)格單元（即正常點(diǎn)）的重疊。填補(bǔ)結(jié)束后，如果希望得到兩層插值邊界，以確保流場(chǎng)求解時(shí)的二階精度，需要將鄰近洞邊界的洞內(nèi)點(diǎn)轉(zhuǎn)變?yōu)椴逯迭c(diǎn)。

割補(bǔ)法在優(yōu)化過(guò)程中，洞邊界點(diǎn)沿網(wǎng)格線離散地作推進(jìn)運(yùn)動(dòng)，每次推進(jìn)的距離就是該段網(wǎng)格線的長(zhǎng)度。當(dāng)網(wǎng)格尺度在各方向上相差較大，或網(wǎng)格形狀匹配性差時(shí)，推進(jìn)陣面在空間上參差不齊，容易有個(gè)別洞邊界點(diǎn)異常，即產(chǎn)生“孤點(diǎn)”，嚴(yán)重時(shí)能導(dǎo)致整個(gè)推進(jìn)陣面崩潰，使重疊失敗。

2 割補(bǔ)法的研究和改進(jìn)

2.1 縫隙問(wèn)題和兩類(lèi)洞內(nèi)點(diǎn)

在割補(bǔ)法的實(shí)際應(yīng)用中發(fā)現(xiàn)［4］，當(dāng)物面距離很近時(shí)，填補(bǔ)過(guò)程容易因?yàn)樯贁?shù)洞邊界點(diǎn)找點(diǎn)過(guò)程的失敗而導(dǎo)致網(wǎng)格進(jìn)入物面內(nèi)。圖1是前緣襟翼網(wǎng)格重疊失敗的示例，切割結(jié)束后，主翼網(wǎng)格的洞邊界緊挨襟翼后緣，在填補(bǔ)過(guò)程中，洞邊界向相鄰的洞內(nèi)點(diǎn)收縮，在襟翼后緣A點(diǎn)處，洞邊界推進(jìn)物面，并迅速將襟翼內(nèi)的全部洞內(nèi)點(diǎn)轉(zhuǎn)換為正常點(diǎn)，使重疊失敗。

圖1 前緣襟翼網(wǎng)格重疊失敗示例Fig.1 Failure of the grid overset for the multi－element airfoil

有經(jīng)驗(yàn)的設(shè)計(jì)人員在搭建網(wǎng)格時(shí)就會(huì)考慮如何避免狹縫重疊失敗，例如將縫隙適當(dāng)增大，或者在縫隙附近加密網(wǎng)格等。但人們更希望通過(guò)對(duì)重疊網(wǎng)格方法的改進(jìn)來(lái)解決這一問(wèn)題。注意到在網(wǎng)格區(qū)域中，洞內(nèi)點(diǎn)可以劃分為兩類(lèi)（參見(jiàn)圖2）Ⅰ型洞內(nèi)點(diǎn)是由挖洞產(chǎn)生的，落入物面內(nèi)部的點(diǎn)；Ⅱ型洞內(nèi)點(diǎn)是洞面優(yōu)化中，在洞邊界逐漸推離物面的過(guò)程中產(chǎn)生的。圖1示例重疊失敗的直接原因就是，落入襟翼內(nèi)部的Ⅰ型洞內(nèi)點(diǎn)在填補(bǔ)過(guò)程中被錯(cuò)誤修改。

圖2 兩類(lèi)洞內(nèi)點(diǎn)Fig.2 Two kinds of hole point

因此，一種十分自然的思想就是在優(yōu)化過(guò)程中對(duì)兩類(lèi)洞內(nèi)點(diǎn)區(qū)別對(duì)待，Ⅰ型洞內(nèi)點(diǎn)被物面屏蔽，不在流場(chǎng)中，洞內(nèi)點(diǎn)屬性應(yīng)受“保護(hù)”不被修改；Ⅱ型洞內(nèi)點(diǎn)是洞面向遠(yuǎn)離物面方向擴(kuò)展的結(jié)果，在填補(bǔ)過(guò)程中，洞面收縮，Ⅱ型洞內(nèi)點(diǎn)屬性是允許修改的。因此，在程序中對(duì)兩類(lèi)洞內(nèi)點(diǎn)進(jìn)行區(qū)分，填補(bǔ)時(shí)，若相鄰的洞內(nèi)點(diǎn)為Ⅰ型洞內(nèi)點(diǎn)，則不修改相鄰點(diǎn)，洞邊界在該處停止移動(dòng)。采用新方法研究前緣襟翼重疊問(wèn)題，參見(jiàn)圖3，狹縫處挖進(jìn)物面的問(wèn)題得到了解決。

2.2 孤點(diǎn)和雙陣面推進(jìn)

目前，二階空間離散格式在工程應(yīng)用中較為廣泛，對(duì)重疊網(wǎng)格，希望在重疊區(qū)域能建立兩層插值邊界，以滿足二階格式在重疊邊界處求解精度的要求。然而在某些情況下，這一要求并不容易滿足。一種常見(jiàn)的情況仍然是縫隙問(wèn)題，圖4是彈身和尾舵網(wǎng)格重疊示例，切割結(jié)束后，彈身網(wǎng)格洞邊界在尾舵前緣處沿軸向向外推出了一層，就因與尾舵網(wǎng)格沿法向推出的洞邊界相遇而停止移動(dòng)，因此彈身網(wǎng)格在尾舵前緣處最多只能建立一層洞邊界，導(dǎo)致局部插值精度降低。

圖3 前緣襟翼網(wǎng)格重疊示例Fig.3 Illustration of the grid overset for the multi－element airfoil

圖4 彈身和尾舵網(wǎng)格切割結(jié)束時(shí)示意Fig.4 Illustration of cutting result for the missile body and rudder configuration

另一種情況，在填補(bǔ)結(jié)束后建立兩層插值邊界時(shí)，第二層邊界上可能出現(xiàn)有插值點(diǎn)找不到有效貢獻(xiàn)單元的問(wèn)題，特別是當(dāng)重疊區(qū)落在網(wǎng)格扭曲或長(zhǎng)細(xì)比較大的區(qū)域時(shí)。圖5是二維情況的示意，由于網(wǎng)格因素使洞邊界參差不齊，但已建立的一層洞邊界，各點(diǎn)均能找到有效貢獻(xiàn)單元，滿足填補(bǔ)過(guò)程迭代終止的條件，因此是一種可能的情況。第二層插值邊界由洞邊界向洞內(nèi)點(diǎn)一側(cè)建立，mesh1中與洞邊界點(diǎn)A相鄰的洞內(nèi)點(diǎn)B，即默認(rèn)的第二層插值點(diǎn)，而B(niǎo)點(diǎn)貢獻(xiàn)單元頂點(diǎn)中包含對(duì)方的洞內(nèi)點(diǎn)C，是一個(gè)無(wú)效的插值點(diǎn)。文獻(xiàn)［5］稱(chēng)之為孤點(diǎn)，一些求解器中也有針對(duì)孤點(diǎn)的解決方案，但流場(chǎng)求解精度在局部的降低是不可避免的。需要注意的是，切割階段也會(huì)產(chǎn)生孤點(diǎn)，文獻(xiàn)［6］給出了有效的解決措施，本文主要討論填補(bǔ)時(shí)的孤點(diǎn)問(wèn)題。

圖5 二維網(wǎng)格填補(bǔ)階段結(jié)束后重疊區(qū)示意Fig.5 The overset area after the paste procedure for the two－dimensional grid

上述兩種情況，前者是兩物體物面距離太近，縫隙處在交叉方向上只覆蓋了一層網(wǎng)格，理論上已經(jīng)無(wú)法建立兩層插值邊界，只能通過(guò)增大縫隙或加密網(wǎng)格來(lái)改善；后者是由于洞邊界落在網(wǎng)格扭曲或網(wǎng)格差異性較大區(qū)域，割補(bǔ)法不能正確處理，因此是本文希望解決的一種情況。

傳統(tǒng)的割補(bǔ)法在優(yōu)化過(guò)程中，只有一層洞邊界“活動(dòng)”，對(duì)默認(rèn)的第二層邊界即相鄰的洞內(nèi)點(diǎn)沒(méi)有約束，因此優(yōu)化結(jié)束時(shí)無(wú)法確保第二層邊界有效。本文提出了一種在填補(bǔ)階段雙陣面推進(jìn)的思想，即將與洞邊界相鄰的洞內(nèi)點(diǎn)作為第二層洞邊界，兩層邊界同時(shí)進(jìn)行檢查和推進(jìn)，當(dāng)兩層邊界均滿足迭代終止的條件時(shí)，填補(bǔ)階段結(jié)束，此時(shí)保留的兩層邊界即最終的插值邊界。

雙陣面推進(jìn)需解決的關(guān)鍵問(wèn)題是當(dāng)?shù)谝粚舆吔琰c(diǎn)滿足要求，而第二層邊界點(diǎn)錯(cuò)誤時(shí)，如何驅(qū)動(dòng)陣面跨過(guò)第一層邊界點(diǎn)。方法的實(shí)現(xiàn)過(guò)程簡(jiǎn)述如下：

1）用切割階段結(jié)束后產(chǎn)生的洞邊界點(diǎn)和相鄰的洞內(nèi)點(diǎn)構(gòu)造陣面，分別記作frg1和frg2，陣面點(diǎn)處于“活動(dòng)”狀態(tài)；

2）尋找陣面點(diǎn)的網(wǎng)格貢獻(xiàn)單元，并依次進(jìn)行以下判斷：

i.若貢獻(xiàn)單元頂點(diǎn)不含洞內(nèi)點(diǎn)，該陣面點(diǎn)進(jìn)入非活動(dòng)狀態(tài)；

ii.若貢獻(xiàn)單元任一頂點(diǎn)是洞內(nèi)點(diǎn)，則若該陣面點(diǎn)是frg1，標(biāo)記為正常點(diǎn)；若該陣面點(diǎn)是frg2，標(biāo)記為臨時(shí)屬性frg＿tmp。

3）檢查陣面點(diǎn)，進(jìn)行以下判斷：

i.若frg1與frg＿tmp相鄰，修改frg1為正常點(diǎn)；

ii.若frg2與正常點(diǎn)相鄰，修改frg2為frg1。

4）再次檢查陣面點(diǎn)，將frg＿tmp修改為活動(dòng)的frg1；若陣面點(diǎn)是frg1，將相鄰的洞內(nèi)點(diǎn)標(biāo)記為活動(dòng)的frg2。

5）重復(fù)上述步驟2至4，直到再也沒(méi)有活動(dòng)的陣面點(diǎn)。

方法在步驟2和步驟3中，當(dāng)frg1滿足要求，而frg2是無(wú)效插值點(diǎn)時(shí)，將frg1修改為正常點(diǎn)，frg2經(jīng)frg＿tmp再修改為frg1，成功的實(shí)現(xiàn)了洞邊界的調(diào)整。

2.3 改進(jìn)的割補(bǔ)法

本文在傳統(tǒng)割補(bǔ)法的基礎(chǔ)上，引入了兩類(lèi)洞內(nèi)點(diǎn)的概念，和雙陣面推進(jìn)的思想，提出了一種改進(jìn)的割補(bǔ)法，其思路參見(jiàn)圖6。新方法將初始挖洞產(chǎn)生的洞內(nèi)點(diǎn)，標(biāo)記為Ⅰ型洞內(nèi)點(diǎn)，并將相鄰的正常點(diǎn)修改為初始洞邊界。在切割過(guò)程中，保持一層洞邊界向遠(yuǎn)離物面的方向推進(jìn)，產(chǎn)生的洞內(nèi)點(diǎn)標(biāo)記為Ⅱ型洞內(nèi)點(diǎn)。填補(bǔ)過(guò)程中，將與第一層洞邊界相鄰的Ⅱ型洞內(nèi)點(diǎn)修改為第二層洞邊界，兩層邊界同時(shí)保持“活動(dòng)”狀態(tài)，進(jìn)行貢獻(xiàn)單元屬性的檢查和判斷，優(yōu)化結(jié)束的條件是兩層邊界點(diǎn)同時(shí)滿足有效插值點(diǎn)的要求，此時(shí)保留的兩層邊界點(diǎn)即為最終的插值邊界。

圖6 改進(jìn)的割補(bǔ)法一維示意Fig.6 One dimensional procedure for the improved cut－paste method

填補(bǔ)過(guò)程中，若洞邊界點(diǎn)有相鄰的Ⅰ型洞內(nèi)點(diǎn)，認(rèn)為洞邊界在局部區(qū)域已推至物面附近，在此處應(yīng)停止推進(jìn)。有一種可能的情況是圖4所示縫隙問(wèn)題，此時(shí)由于兩類(lèi)洞內(nèi)點(diǎn)被區(qū)分，不會(huì)出現(xiàn)挖進(jìn)物面的問(wèn)題，但局部可能無(wú)法建立兩層插值邊界，程序中會(huì)作一些適當(dāng)?shù)奶幚?，以確保在以局部精度損失為代價(jià)的條件下，使重疊過(guò)程得以繼續(xù)。

3 計(jì)算方法

3.1 重疊網(wǎng)格方法

本文改進(jìn)的割補(bǔ)法技術(shù)已在MI－GRID中得到應(yīng)用。MI－GRID是北航閻超課題組研制的重疊網(wǎng)格軟件，核心模塊采用洞映射和割補(bǔ)法技術(shù)，包含了孤點(diǎn)清除、體積優(yōu)化、物面重疊、動(dòng)態(tài)重疊等方法［7－8］。該重疊網(wǎng)格軟件計(jì)算效率高、可靠性好、使用方便，先后參與了國(guó)內(nèi)多個(gè)航空航天型號(hào)研制工作，包括復(fù)雜外形飛行器、子母彈拋撒、助推級(jí)分離、折疊翼打開(kāi)等項(xiàng)目，得到了有效考核。

3.2 數(shù)值方法

本文算例均求解雷諾平均NS控制方程，空間離散采用Roe的FDS格式，MUSCL插值方法和Van Albada限制器用于獲得二階空間離散精度；湍流模型采用SST剪切應(yīng)力輸運(yùn)模型；時(shí)間離散采用穩(wěn)定性高的LU－SGS隱式計(jì)算方法。關(guān)于數(shù)值方法具體可參考文獻(xiàn)［9］。

4 算例研究

4.1 Titan IV運(yùn)載火箭

本文對(duì)美國(guó)大力神四號(hào)（TitanⅣ）大型捆綁式運(yùn)載火箭的超聲速繞流問(wèn)題［10］進(jìn)行數(shù)值模擬研究。計(jì)算條件為：M∞＝1.6，ReL＝1.1×107，α＝0°，由于流動(dòng)條件對(duì)稱(chēng)，故使用半模計(jì)算。采用重疊網(wǎng)格方法分別生成芯級(jí)和助推級(jí)的計(jì)算網(wǎng)格，網(wǎng)格數(shù)目為190萬(wàn)和40萬(wàn)。

TitanⅣ運(yùn)載火箭采用并聯(lián)式布局，芯級(jí)和助推級(jí)距離較近，是典型的狹縫重疊問(wèn)題。圖7是本文方法建立的對(duì)稱(chēng)面重疊網(wǎng)格和流場(chǎng)結(jié)果，芯級(jí)和助推級(jí)之間的縫隙得到了準(zhǔn)確描述，重疊結(jié)果正確，流場(chǎng)結(jié)構(gòu)清晰合理。圖8是本文方法和傳統(tǒng)割補(bǔ)法重疊結(jié)果的比較，顯然傳統(tǒng)割補(bǔ)法沒(méi)有區(qū)分兩類(lèi)洞內(nèi)點(diǎn)，在狹縫重疊時(shí)因填補(bǔ)中陣面點(diǎn)判斷錯(cuò)誤，使陣面推進(jìn)物面，導(dǎo)致重疊失敗。

圖9是火箭芯級(jí)中心線上的壓力分布與實(shí)驗(yàn)值的比較，計(jì)算結(jié)果與實(shí)驗(yàn)值吻合很好，說(shuō)明本文改進(jìn)的方法重疊結(jié)果較好。圖9中，因?yàn)樵陲L(fēng)洞實(shí)驗(yàn)中，芯級(jí)與助推級(jí)間存在連接機(jī)構(gòu)，導(dǎo)致實(shí)驗(yàn)的峰值點(diǎn)略高。

圖7 對(duì)稱(chēng)面重疊網(wǎng)格和等馬赫線流場(chǎng)圖Fig.7 Illustration of the overset grid and Mach contours for the symmetry plane

圖8 不同方法重疊結(jié)果比較Fig.8 Comparison of the different overset methods

圖9 芯級(jí)中心線壓力分布Fig.9 Pressure distribution of the central stage

4.2 子母彈拋撒算例

子母彈拋撒的空氣動(dòng)力學(xué)問(wèn)題是一類(lèi)典型的超聲速或高超聲速多體干擾、非定常復(fù)雜流動(dòng)問(wèn)題。采用重疊網(wǎng)格方法求解子母彈拋撒問(wèn)題，已成為國(guó)內(nèi)外最常用的方法［11－12］。本文對(duì)典型外形的子母彈拋撒問(wèn)題進(jìn)行了研究，母彈模型包含彈倉(cāng)，子彈模型包含呈十字布局的四片尾翼。計(jì)算條件為：拋撒馬赫數(shù)3.0，拋撒高度10km，母彈迎角0°。母彈、子彈和四片尾翼分別生成計(jì)算網(wǎng)格，其中母彈網(wǎng)格數(shù)目為320萬(wàn)，子彈網(wǎng)格數(shù)目為42萬(wàn)，單個(gè)尾翼網(wǎng)格數(shù)目30萬(wàn)。

圖10、圖11是采用本文方法建立的子母彈模型重疊網(wǎng)格示意。在彈倉(cāng)附近，由于網(wǎng)格過(guò)渡需要，母彈網(wǎng)格沿徑向較密，使重疊區(qū)附近母彈網(wǎng)格各方向尺度，及母彈網(wǎng)格與子彈網(wǎng)格之間的差異較大，這對(duì)割補(bǔ)法一類(lèi)按邏輯坐標(biāo)推進(jìn)的迭代優(yōu)化方法是不利的，但本文方法仍能正確生成重疊網(wǎng)格。圖12對(duì)子彈模型重疊結(jié)果的比較則更為直觀，傳統(tǒng)割補(bǔ)法產(chǎn)生的重疊區(qū)在上述不利區(qū)域容易產(chǎn)生孤點(diǎn)，而改進(jìn)方法在同樣情況下，在第二層邊界建立失敗時(shí)，能由第二層邊界驅(qū)動(dòng)重疊陣面繼續(xù)推進(jìn)，避免了孤點(diǎn)的產(chǎn)生，使重疊結(jié)果較好，重疊邊界整潔。

圖10 對(duì)稱(chēng)面重疊網(wǎng)格示意Fig.10 Illustration of the overset grid for the symmetry plane

圖11 截面重疊網(wǎng)格示意Fig.11 Illustration of the overset grid for the truncated plane

圖12 子彈模型不同方法重疊結(jié)果比較Fig.12 Comparison of the different overset methods for missile configuration

圖13是由數(shù)值計(jì)算獲得的對(duì)稱(chēng)面和截面流場(chǎng)圖，圖中流場(chǎng)等值線在重疊區(qū)銜接較為光滑，流場(chǎng)結(jié)構(gòu)清晰，說(shuō)明本文使用的重疊網(wǎng)格方法和CFD求解方法對(duì)子母彈拋撒流場(chǎng)具有較強(qiáng)的解算能力。

圖13 對(duì)稱(chēng)面等馬赫線流場(chǎng)圖Fig.13 Mach contours of the symmetry plane

5 結(jié) 論

本文對(duì)割補(bǔ)法進(jìn)行了研究，指出傳統(tǒng)的割補(bǔ)法在處理縫隙、網(wǎng)格扭曲、網(wǎng)格差異性大等情況時(shí)可靠性較差，在此基礎(chǔ)上進(jìn)行了改進(jìn)，主要有：

（1）提出了兩類(lèi)洞內(nèi)點(diǎn)的概念，將在物理上被屏蔽的洞內(nèi)點(diǎn)和切割過(guò)程中產(chǎn)生的洞內(nèi)點(diǎn)區(qū)分，前者在填補(bǔ)階段不允許修改，從根本上避免了挖進(jìn)物面的問(wèn)題；

（2）針對(duì)網(wǎng)格扭曲或網(wǎng)格差異大，可能無(wú)法建立兩層有效插值邊界的問(wèn)題，提出了雙陣面聯(lián)動(dòng)推進(jìn)的思想，當(dāng)?shù)诙舆吔琰c(diǎn)不滿足要求時(shí)，能有效驅(qū)動(dòng)陣面繼續(xù)推進(jìn)，避免了第二層插值邊界無(wú)法建立的問(wèn)題；

（3）經(jīng)算例研究，本文針對(duì)割補(bǔ)法提出的改進(jìn)能有效解決狹縫重疊、網(wǎng)格扭曲等復(fù)雜網(wǎng)格情況。

［1］RALPH W N，JEFFREY P S.A summary of the 2004 overset symposium on composite grids and solution technology［R］.AIAA 2005－921.

［2］WEY T C.Development of a mesh interface generator for overlapped structured grids［R］.AIAA 1994－1924.

［3］CHO K W，KWON J H.Development of a fully systemized chimera methodology for steady／unsteady problems［J］.JournalofAircraft，1999，36（6）：973－980.

［4］FAN J J，YAN C.Enhancement and application of overset grid assembly［J］.ChineseJournalofAeronautics，2010，23（6）：631－638.

［5］RALPH W N，DAVID A B，ROBERT F K.Suggar＋＋：an improved general overset grid assembly capability［R］.AIAA 2009－3992.

［6］李亭鶴，閻超，李躍軍.重疊網(wǎng)格技術(shù)中割補(bǔ)法的研究與改進(jìn)［J］.北京航空航天大學(xué)學(xué)報(bào)，2005，31（4）：402－406.

［7］李亭鶴.重疊網(wǎng)格自動(dòng)生成方法研究［D］.北京航空航天大學(xué)，2004.

［8］范晶晶.復(fù)雜重疊網(wǎng)格方法研究及多體運(yùn)動(dòng)的非定常流動(dòng)模擬［D］.北京航空航天大學(xué)，2010.

［9］閻超.計(jì)算流體力學(xué)方法及應(yīng)用［M］.北京航空航天大學(xué)出版社：2006，245－248.

［10］STEPHEN T，JOHNSON C T.Launch－vehicle simulations using a concurrent，implicit Navier－Stokes solver［R］.AIAA 95－0223.

［11］SPINETTI R L，JOLLY B A.Time－accurate numerical simulation of GBU－38sseparating from the B－1Baircraft with various ejector forces，store properties，and loadout configurations［R］.AIAA 2008－187.

［12］張輝.重疊網(wǎng)格技術(shù)在飛行器多體分離問(wèn)題中的應(yīng)用［D］.北京航空航天大學(xué)，2010.