超聲速降落傘系統(tǒng)的氣動(dòng)干擾數(shù)值模擬研究

薛曉鵬溫志湧

?

超聲速降落傘系統(tǒng)的氣動(dòng)干擾數(shù)值模擬研究

薛曉鵬1溫志湧2

（1 中南大學(xué)航空航天學(xué)院，長(zhǎng)沙410083）（2 香港理工大學(xué)深圳研究院，深圳 518057）

文章基于一種簡(jiǎn)易“浸入邊界技術(shù)”與流固耦合方法對(duì)超聲速來(lái)流條件下的三維降落傘系統(tǒng)進(jìn)行了數(shù)值模擬。文章分析了降落傘系統(tǒng)（太空艙和降落傘衣）流場(chǎng)產(chǎn)生的氣動(dòng)干擾對(duì)傘衣外形變形動(dòng)力學(xué)的影響以及降落傘設(shè)計(jì)參數(shù)（如拖拽距離參數(shù)）對(duì)降落傘性能表現(xiàn)的影響，同時(shí)檢驗(yàn)了分離渦模擬方法在三維超聲速降落傘數(shù)值模擬中的適用性。結(jié)果表明，當(dāng)拖拽距離參數(shù)很小時(shí)，傘衣出現(xiàn)了嚴(yán)重的收縮現(xiàn)象；當(dāng)距離增大時(shí)，觀測(cè)到了傘開(kāi)口的“呼吸現(xiàn)象”，阻力系數(shù)獲得了較大改善。另外，在拖拽距離參數(shù)較大的情況下，由于湍流尾流與傘前激波的氣動(dòng)干擾影響明顯，采用分離渦模擬方法，可以獲得較層流尾流條件下更大的阻力系數(shù)。

氣動(dòng)干擾分離 渦模擬 流固耦合 超聲速降落傘

0 引言

在NASA的“火星科學(xué)實(shí)驗(yàn)室”（MSL）任務(wù)中，“好奇號(hào)”探測(cè)器已經(jīng)成功地著陸在火星表面?！昂闷嫣?hào)”在進(jìn)入火星大氣層以后的從超聲速到亞聲速的減速過(guò)程是由一個(gè)超聲速降落傘系統(tǒng)所完成的[1]。

從20世紀(jì)60年代末開(kāi)始，超聲速降落傘的研究工作吸引了全世界的研究學(xué)者。文獻(xiàn)[2]用風(fēng)洞試驗(yàn)方法檢驗(yàn)了前置體尾流影響下的多種類(lèi)型的超聲速盤(pán)縫帶傘系統(tǒng)的氣動(dòng)特性。文獻(xiàn)[3-4]利用計(jì)算流體力學(xué)方法和流固耦合方法（Fluid Structure Interaction，F(xiàn)SI）首次對(duì)柔性超聲速降落傘進(jìn)行了數(shù)值模擬，并觀測(cè)到了其周?chē)晕擦鳌げㄏ嗷プ饔脼橹饕卣鞯膹?fù)雜流場(chǎng)結(jié)構(gòu)，同時(shí)分析了馬赫數(shù)和拖拽距離參數(shù)對(duì)傘衣阻力系數(shù)的影響。近年來(lái)，美國(guó)NASA及其合作小組對(duì)MSL降落傘模型進(jìn)行了超聲速條件下試驗(yàn)和數(shù)值的充分調(diào)查[5-7]，發(fā)現(xiàn)降落傘不穩(wěn)定主要來(lái)自于降落傘前激波與太空艙尾流的氣動(dòng)干擾，其將導(dǎo)致傘前激波形狀的變化，進(jìn)而減少質(zhì)量流流進(jìn)傘內(nèi)部，并且該非定常氣動(dòng)干擾受馬赫數(shù)、雷諾數(shù)、太空艙形狀以及艙與傘的距離等因素的影響。文獻(xiàn)[8]使用大渦模擬方法和基于有限元模型的結(jié)構(gòu)膜對(duì)大尺度超聲速盤(pán)縫帶傘進(jìn)行了數(shù)值模擬，并成功觀測(cè)到了降落傘體的“呼吸現(xiàn)象”（開(kāi)口面積出現(xiàn)在一定范圍內(nèi)波動(dòng)）。文獻(xiàn)[9]使用浸入邊界技術(shù)以及流固耦合方法對(duì)三維超聲速降落傘進(jìn)行數(shù)值模擬，進(jìn)一步分析了拖拽距離參數(shù)、太空艙以及傘體形狀、馬赫數(shù)對(duì)降落傘性能的影響。近年來(lái)，我國(guó)有關(guān)超聲速降落傘及火星環(huán)境下的減速著陸研究也開(kāi)始啟動(dòng)：文獻(xiàn)[10]利用徑向—軸向動(dòng)量守恒充氣模型研究大氣密度對(duì)降落傘的開(kāi)傘過(guò)程的影響，結(jié)果表明降落傘的充氣時(shí)間和充氣距離隨大氣密度的減小而增大；文獻(xiàn)[11]進(jìn)一步計(jì)算了低密度大氣中盤(pán)縫帶降落傘的開(kāi)傘動(dòng)載；文獻(xiàn)[12]進(jìn)一步簡(jiǎn)化了降落傘的軸向—徑向動(dòng)量方程，分析了不同大氣密度、初始開(kāi)傘速度下的充氣環(huán)境對(duì)降落傘開(kāi)傘過(guò)程的影響；文獻(xiàn)[13]探討了結(jié)構(gòu)參數(shù)對(duì)火星用傘開(kāi)傘性能的影響，發(fā)現(xiàn)當(dāng)盤(pán)縫帶傘的帶寬和縫寬增加，該傘的阻力系數(shù)隨之減少；文獻(xiàn)[14]通過(guò)高速風(fēng)洞試驗(yàn)以及高速飛行投放試驗(yàn)對(duì)超聲速半流傘進(jìn)行了結(jié)構(gòu)設(shè)計(jì)和氣動(dòng)特性分析，得到了半流傘的擺角參數(shù)和阻力系數(shù)變化規(guī)律。

隨著計(jì)算機(jī)技術(shù)的發(fā)展，降落傘問(wèn)題的數(shù)值模擬已經(jīng)取得了顯著的進(jìn)步。一些成功的方法，比如變空間域/穩(wěn)定時(shí)間—空間（Deforming-spatial-domain/stabilized Space-time，DSD/SST）方法[15]、任意拉格朗日—?dú)W拉法（Arbitrary Langrangian Eulerian，ALE）[16-17]、浸入邊界法（Immersed Boundary Method，IBM）[18]等，已經(jīng)成功應(yīng)用到亞聲速降落傘數(shù)值模擬中。但是，至今還沒(méi)有合適的數(shù)值方法求解超聲速降落傘問(wèn)題。Lingard在超聲速降落傘的數(shù)值模擬中采用了ALE方法，然而未能展示復(fù)雜流動(dòng)現(xiàn)象影響下傘衣的外形變化[3-4]。本研究中我們提出了一種簡(jiǎn)易浸入邊界技術(shù)（Immersed Boundary Technique，IBT）應(yīng)用到涉及超聲速降落傘模型的流固耦合問(wèn)題中。該方法已經(jīng)被成功應(yīng)用到亞聲速超聲速降落傘系統(tǒng)的二維和三維的數(shù)值模擬[19]中，模擬所得的開(kāi)傘過(guò)程與試驗(yàn)結(jié)果非常吻合。該方法也被應(yīng)用到了超聲速降落傘的數(shù)值模擬[9]，并得到了與美國(guó)NASA的噴氣推進(jìn)實(shí)驗(yàn)室的超聲速降落傘試驗(yàn)一致的結(jié)果。另外，該技術(shù)也被應(yīng)用在降落傘傘繩在超聲速流場(chǎng)中的影響的數(shù)值模擬中，并觀測(cè)到了傘繩激波，其與美國(guó)NASA噴氣推進(jìn)實(shí)驗(yàn)室對(duì)該研究的試驗(yàn)結(jié)果保持一致[20]。

在高馬赫數(shù)、高雷諾數(shù)情況下，太空艙尾流將是復(fù)雜的湍流流動(dòng)，降落傘周?chē)鲌?chǎng)將出現(xiàn)強(qiáng)烈的太空艙湍流尾流與傘前激波的氣動(dòng)干擾并伴有復(fù)雜的流場(chǎng)結(jié)構(gòu)。本研究將對(duì)三維超聲速降落傘進(jìn)行數(shù)值模擬，分析氣動(dòng)干擾對(duì)降落傘傘衣變形動(dòng)力學(xué)的影響，以及降落傘設(shè)計(jì)參數(shù)如拖拽距離參數(shù)對(duì)降落傘性能表現(xiàn)的影響，同時(shí)檢驗(yàn)分離渦模擬（Detached Eddy Simulation，DES）方法在三維超聲速降落傘數(shù)值模擬中的適用性。

1 三維降落傘模型

本研究中的三維降落傘系統(tǒng)包括太空艙和降落傘傘衣，如圖1所示。傘衣和太空艙之間用傘繩連接，這里傘繩的作用僅是維持傘體和太空艙之間的受力平衡。傘繩以及透氣量對(duì)流場(chǎng)的影響暫且不考慮。

圖1 三維降落傘模型

傘體是一個(gè)直徑為的半球簡(jiǎn)化模型，太空艙是一個(gè)呈錐形的簡(jiǎn)化模型，半錐角為20°，前表面的直徑是=24mm。是從太空艙前表面到傘體入口的直線距離。是降落傘的關(guān)鍵設(shè)計(jì)參數(shù)，即拖拽距離參數(shù)。本研究中模型A和B的拖拽距離參數(shù)如表1所示，其數(shù)值小于NASA的MSL降落傘模型試驗(yàn)中的名義值（大約為10），這是因?yàn)樵谛∮?0的情況下，更加復(fù)雜的氣動(dòng)干擾在以前剛性處理的降落傘模型研究中被觀測(cè)到[21]，因此發(fā)生在更復(fù)雜情況下的柔性降落傘的性能表現(xiàn)是值得研究的。

表1 降落傘系統(tǒng)外形參數(shù)

Tab.1 Specifications for the different parachute systems

2 數(shù)值方法

2.1 計(jì)算條件和計(jì)算網(wǎng)格

本研究計(jì)算中來(lái)流條件如表2所示，來(lái)自于文獻(xiàn)[22]中的試驗(yàn)條件。由于降落傘系統(tǒng)的軸對(duì)稱(chēng)特性，計(jì)算網(wǎng)格由子午線斷面旋轉(zhuǎn)而成，如圖2所示，圖中紅色部分是傘體，藍(lán)色部分是太空艙。網(wǎng)格單元數(shù)約1×106個(gè)。降落傘傘體的控制點(diǎn)數(shù)目為902×441，其在傘體區(qū)域被搜索定義，網(wǎng)格依賴性檢驗(yàn)參考文獻(xiàn)[9]。

表2 本研究計(jì)算所用的來(lái)流條件

Tab.2 Free stream conditions employed in this study

圖2 降落傘的計(jì)算網(wǎng)格

2.2 數(shù)值方法

（1）流場(chǎng)計(jì)算

本研究采用三維可壓縮Navier-Stokes方程數(shù)值求解三維降落傘系統(tǒng)周?chē)某曀倭鲌?chǎng)，控制方程采用有限體積法進(jìn)行離散。無(wú)粘通量選擇SHUS（Simple High-resolution Upwind Scheme）[23]格式，并采用Van Albada限制器通過(guò)3階MUSCL（Monotone Upstream-centered Schemes for Conservation Law）方法進(jìn)行差值提高精度。粘性通量采用2階中心格式離散，時(shí)間推進(jìn)則采用LU-SGS（Lower-upper Symmetric Gauss-seidel Method）格式進(jìn)行。邊界條件的處理上，太空艙面采用無(wú)滑移、等壁溫。入口邊界賦來(lái)流參數(shù)，出口邊界采用外插處理。

本研究采用了“浸入邊界技術(shù)”來(lái)處理三維降落傘的柔性傘體部分。該技術(shù)是由Ochi[24]等人提出，由Miyoshi、Xue等人發(fā)展而成適用于低速、高速流場(chǎng)中降落傘的流固耦合研究中的一種簡(jiǎn)易技術(shù)[9,19]。該浸入邊界技術(shù)有別于傳統(tǒng)的“浸入邊界法”[25]。因?yàn)閭鹘y(tǒng)的浸入邊界法受限于低雷諾數(shù)和不可壓流，且在流場(chǎng)控制方程中增加了力源項(xiàng)，該力源項(xiàng)在嚴(yán)格限制中常常表現(xiàn)不佳。而本浸入邊界技術(shù)未在控制方程中引入一個(gè)力源項(xiàng)，而是通過(guò)流體網(wǎng)格和虛擬網(wǎng)格的關(guān)系來(lái)求得虛擬網(wǎng)格的速度向量，如圖3所示，其求解方程[9,19]為：

V=V-2(V·n)n+V（1）

式中V,V, V分別代表流體網(wǎng)格、虛擬網(wǎng)格和傘控制點(diǎn)的速度；V在壁面上可以分解為垂直方向和平行方向分量V和V；n代表傘表面的單位法向量；下標(biāo)和分別表示流體網(wǎng)格和虛擬網(wǎng)格的變量。式（1）展開(kāi)以后正如圖3中右圖所示的關(guān)系[9,19]。方程的詳細(xì)推導(dǎo)過(guò)程參考文獻(xiàn)[24]。

圖3 虛擬網(wǎng)格和流體網(wǎng)格的速度之間的關(guān)系[9,19]

太空艙的尾流將呈以大尺度渦為特征的高度湍流狀態(tài)，并與來(lái)自傘體的激波發(fā)生激烈氣動(dòng)干擾[5]。本文采用DES方法[26-27]對(duì)其進(jìn)行湍流模擬。DES是一種混合模型，其在計(jì)算時(shí)，近壁面湍流結(jié)構(gòu)的計(jì)算求解雷諾平均N-S方程（Reynolds Averaged Navier-stokes Equation，RANS）得到耗散能量的小渦，同時(shí)遠(yuǎn)離壁面的分離渦明顯的區(qū)域則采用大渦模擬方法（Large Eddy Simulation, LES）求解。本研究中的RANS模型采用的是方程Spalart-Allmaras（S-A）模型。DES處理RANS和LES過(guò)渡區(qū)域的方法是借用S-A模型控制方程中的參數(shù)，這一參數(shù)表示了網(wǎng)格節(jié)點(diǎn)與最近壁面的距離。因?yàn)樾枰鶕?jù)不同距壁面距離來(lái)確定相應(yīng)的求解方法，DES方法對(duì)參數(shù)進(jìn)行了修正，定義了新的長(zhǎng)度尺度,其由壁面距離和當(dāng)?shù)鼐W(wǎng)格尺度共同確定：

（2）結(jié)構(gòu)計(jì)算

在本研究中，降落傘模型的結(jié)構(gòu)計(jì)算采用了質(zhì)量?彈簧?阻尼模型[9,19,28]來(lái)模擬傘體的結(jié)構(gòu)動(dòng)力學(xué)，如圖4所示。在該模型中，傘結(jié)構(gòu)被處理為質(zhì)點(diǎn)與彈簧和阻尼的聯(lián)合體，特別在傘的邊緣，需要考慮來(lái)自傘繩的拉力。其控制方程是基于作用在傘的每個(gè)控制點(diǎn)的牛頓第二定律。傘體結(jié)構(gòu)計(jì)算的時(shí)間推進(jìn)算法采用顯性2階Runge-Kutta 格式。

（a）傘內(nèi)部的MSD模型 （b）傘邊緣的MSD模型

（a）General MSD model （b）Special MSD model on the edge of canopy

圖4 傘的質(zhì)量?彈簧?阻尼結(jié)構(gòu)模型[9,19]

Fig.4 The MSD model

（3）流固耦合格式

流固耦合方法如圖5所示，傘體表面上的壓力分布做為流體力，被應(yīng)用到每一個(gè)傘體控制點(diǎn)的位移和速度計(jì)算中[9,19]，然后這些計(jì)算結(jié)果作為邊界條件再通過(guò)耦合方法“浸入邊界技術(shù)（IBT）”[9,19]被傳送到流場(chǎng)計(jì)算中。

為了同時(shí)計(jì)算流體和結(jié)構(gòu)，弱耦合方法（如圖6所示，圖中為時(shí)間步長(zhǎng)）被應(yīng)用到該計(jì)算中，因?yàn)閭泱w變形和非定常流場(chǎng)的相互影響非常敏感[9]。

圖5 流固耦合方法[9]

圖6 弱耦合格式[9]

3 計(jì)算結(jié)果與分析

3.1 模型A：=57mm，=2.375

圖7 模型A傘衣的外形變化

圖8 模型A降落傘系統(tǒng)的瞬時(shí)馬赫數(shù)等值線

強(qiáng)烈的流場(chǎng)不穩(wěn)定會(huì)導(dǎo)致傘衣的阻力d出現(xiàn)很大的不穩(wěn)定狀態(tài)，該阻力作用在來(lái)流方向，可由傘衣表面的內(nèi)外壓力差計(jì)算。然而，為了更精確的定義和比較，阻力系數(shù)被用來(lái)表征降落傘系統(tǒng)的性能。阻力系數(shù)d[7]定義如下：

圖9比較了模型A的湍流模型（DES）和層流模型（Laminar）的阻力系數(shù)曲線，圖中層流結(jié)果（Laminar）來(lái)自文獻(xiàn)[9]。從圖中可以發(fā)現(xiàn)，大約0.003~0.004s開(kāi)始，其阻力系數(shù)已經(jīng)基本上減小至零，這是因?yàn)樵撃Ｐ蛡阋略诔掷m(xù)地收縮，使得內(nèi)外壓差持續(xù)減小至零。在0.003s之前，由于0s時(shí)候的傘衣處于半球形狀，傘前激波與太空艙尾流發(fā)生激烈氣動(dòng)干擾，導(dǎo)致傘內(nèi)壓力急劇上升，并得到極大的阻力系數(shù)。當(dāng)傘內(nèi)壓力達(dá)到波谷的時(shí)候，傘的阻力變得很小，嚴(yán)重的傘衣收縮使得阻力系數(shù)逐漸減小至零。

圖9 模型A（X=57mm）降落傘系統(tǒng)的阻力系數(shù)

另外值得注意的是，與文獻(xiàn)[9]中層流計(jì)算結(jié)果相比較，傘衣外形變化、降落傘周?chē)鲌?chǎng)結(jié)構(gòu)與氣動(dòng)干擾位置和阻力系數(shù)變化受湍流計(jì)算影響較小。所以，拖拽系數(shù)較小時(shí)，太空艙尾流與傘前激波始終在傘端部發(fā)生氣動(dòng)干擾，導(dǎo)致傘衣嚴(yán)重收縮，進(jìn)而使其阻力系數(shù)急劇減小，降落傘性能失效。

3.2 模型B：=171mm，=7.125

圖10為模型B的傘衣外形變化，與模型A相比，可以觀察到傘衣開(kāi)口呈現(xiàn)出“呼吸現(xiàn)象”。圖11為模型B周?chē)牧鲌?chǎng)分布，傘前激波在傘開(kāi)口的逆流方向與來(lái)自太空艙的湍流尾流發(fā)生劇烈氣動(dòng)干擾，這是降落傘周?chē)鲌?chǎng)不穩(wěn)定的主要來(lái)源[7-8]。因?yàn)橥献Ь嚯x參數(shù)增大，所以太空艙前激波呈現(xiàn)穩(wěn)定狀態(tài)，未與傘前激波發(fā)生干擾。圖12為該模型的阻力系數(shù)曲線（圖中層流結(jié)果（Laminar）來(lái)自文獻(xiàn)[9]），與圖9模型A相比，可以發(fā)現(xiàn)模型B產(chǎn)生了較高的阻力系數(shù)，并出現(xiàn)周期變化。所以，適當(dāng)增加拖拽距離系數(shù)，降落傘系統(tǒng)可以得到較高的阻力系數(shù)，這與Reichenau[2]的實(shí)驗(yàn)發(fā)現(xiàn)是吻合的，與文獻(xiàn)[9]中層流情況下的計(jì)算結(jié)果也是保持一致的。

另外，與文獻(xiàn)[9]層流結(jié)果相比，傘衣外形變化與流場(chǎng)結(jié)構(gòu)特征保持一致，從圖12中層流（Laminar）與湍流（DES）模型結(jié)果相比較，可以發(fā)現(xiàn)穩(wěn)定階段（0.003s以后）湍流計(jì)算所得阻力系數(shù)與層流結(jié)果變化趨勢(shì)基本保持一致，但是幅值要較層流大。這說(shuō)明太空艙后的湍流尾流與激波氣動(dòng)干擾較層流尾流情況下增強(qiáng)，阻力系數(shù)顯著增大。

圖10 模型B傘衣的外形變化

圖11 模型B降落傘系統(tǒng)的瞬時(shí)馬赫數(shù)等值線

圖12 模型B（X=171mm）降落傘系統(tǒng)的阻力系數(shù)

4 結(jié)束語(yǔ)

本文主要分析了降落傘系統(tǒng)周?chē)鷱?fù)雜的氣動(dòng)干擾對(duì)降落傘傘體變形動(dòng)力學(xué)的影響，以及降落傘設(shè)計(jì)參數(shù)，如拖拽距離參數(shù)對(duì)降落傘性能表現(xiàn)的影響。在高馬赫數(shù)、高雷諾數(shù)的條件下，檢驗(yàn)了分離渦模擬方法在三維超聲速降落傘數(shù)值模擬中的必要性與適用性。結(jié)果表明，影響降落傘性能表現(xiàn)的兩大關(guān)鍵因素是降落傘的非定常變形和太空艙尾流與傘前激波的氣動(dòng)干擾。當(dāng)拖拽距離參數(shù)很小時(shí)，傘出現(xiàn)了嚴(yán)重的收縮現(xiàn)象；當(dāng)該距離增大時(shí)，觀測(cè)到了傘開(kāi)口的“呼吸現(xiàn)象”，并且阻力系數(shù)獲得了較大改善。分離渦模擬方法對(duì)大拖拽距離系數(shù)的降落傘系統(tǒng)周?chē)鲌?chǎng)復(fù)雜氣動(dòng)干擾更敏感，進(jìn)而影響實(shí)際阻力系數(shù)的精確度。

保護(hù)區(qū)，指對(duì)水資源保護(hù)、自然生態(tài)系統(tǒng)及珍稀瀕危物種的保護(hù)有重要意義，需劃定進(jìn)行保護(hù)的水域。禁止在飲用水水源一級(jí)保護(hù)區(qū)、自然保護(hù)區(qū)核心區(qū)等范圍內(nèi)新建、改建、擴(kuò)建與保護(hù)無(wú)關(guān)的建設(shè)項(xiàng)目和從事與保護(hù)無(wú)關(guān)的涉水活動(dòng)。

[1] SENGUPTA A, STELTZNER A. Results from the Mars Science Laboratory Parachute Decelerator System Supersonic Qualification Program[C]//Aerospace Conference, Big Sky, MT: IEEE, 2008. DOI: 10.1109/AERO.2008.452684.

[2] REICHENAU D E. Aerodynamic Characteristics of Disk-gap-band Parachutes in the Wake of Viking Entry Forebodies at Mach Numbers from 0.2 to 2.6[R]. AEDC-TR-72-78, Tennessee: Arnold AFB, 1972.

[3] LINGARD J, DARLEY M. Simulation of Parachute Fluid Structure Interaction in Supersonic Flow[C]//18th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, AIAA Paper 2005-1607, Munich, Germany, 2005.

[4] LINGARD J, DARLEY M, UNDERWOOD J C. Simulation of Mars Supersonic Parachute Performance and Dynamics[C]// 19th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, AIAA Paper 2007-2507, Williamsburg, VA, 2007.

[5] SENGUPTA A. Fluid Structure Interaction of Parachutes in Supersonic Planetary Entry[C]//21st AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, AIAA Paper 2011-2541, Dublin, Ireland, 2011.

[6] SENGUPTA A, KELSCH R. Supersonic Performance of Disk-gap-band Parachutes Constrained to a 0-Degree Trim Angle[J]. Journal of Spacecraft and Rockets, 2009, 46(6): 1155-1163.

[7] SENGUPTA A, STELZNER A. Findings from the Supersonic Qualification Program of the Mars Science Laboratory Parachute System[C]//20th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, AIAA Paper 2009-2900, Seattle, Washington, 2009.

[8] KARAGIOZIS K, KAMAKOTI R, CIRAK F, et al. A Computational Study of Supersonic Disk-gap-band Parachutes Using Large-eddy Simulation Coupled to a Structural Membrane[J]. Journal of Fluids and Structures, 2011, 27(2): 175-192.

[9] XUE X, NAKAMURA Y. Numerical Simulation of a Three-dimensional Flexible Parachute System under Supersonic Conditions[J]. Trans. JSASS Aerospace Tech. Japan, 2013, 11: 99-108.

[10] 榮偉, 陳旭, 陳國(guó)良. 大氣密度對(duì)降落傘充氣性能的影響[J]. 航天返回與遙感, 2006, 27(3): 11-16. RONG Wei, CHEN Xu, CHEN Guoliang. The Effect of Atmospheric Density on Parachute Inflation Performances[J]. Spacecraft Recovery & Remote Sensing, 2006, 27(3): 11-16. (in Chinese)

[11] 榮偉, 陳旭, 陳國(guó)良. 低密度大氣中降落傘開(kāi)傘動(dòng)載的研究[J]. 航天返回與遙感, 2006, 27(4): 7-11. RONG Wei, CHEN Xu, CHEN Guoliang. The Study of the Parachute Opening Load in Low Atmospheric Density[J]. Spacecraft Recovery & Remote Sensing, 2006, 27(4): 7-11. (in Chinese)

[12] 郭樹(shù)偉, 董楊彪, 王海濤, 等. 降落傘充氣環(huán)境對(duì)充氣性能的影響[J]. 中國(guó)空間科學(xué)技術(shù), 2008, 28(6): 45-51. GUO Shuwei, DONG Yangbiao, WANG Haitao, et al. Influence of Inflation Conditions on Parachute Inflation[J]. Chinese Space Science and Technology, 2008, 28(6): 45-51. (in Chinese)

[13] 連亮, 張紅英, 李方. 結(jié)構(gòu)參數(shù)對(duì)火星探測(cè)用傘開(kāi)傘性能的影響[J]. 航天返回與遙感, 2012, 33(6): 24-29. LIAN Liang, ZHANG Hongying, LI Fang. Influence of Structure Parameters on Inflation Performance of Parachute Used for Mars Exploration[J]. Spacecraft Recovery & Remote Sensing, 2012, 33(6): 24-29. (in Chinese)

[14] 韓晉陽(yáng), 徐宏, 高峰. 超聲速半流傘設(shè)計(jì)與分析[J]. 航天返回與遙感, 2013, 34(5): 20-28. HAN Jinyang, XU Hong, GAO Feng. Design and Analysis of Supersonic Half-flow Parachute[J]. Spacecraft Recovery & Remote Sensing, 2013, 34(5): 20-28. (in Chinese)

[15] TEZDUYAR T E, BEHR M, LIOU J. A New Strategy for Finite Element Computations Involving Moving Boundaries and Interfacesthe Deforming-spatial-domain/space-time Procedure: I. the Concept and the Preliminary Numerical Tests[J]. Computer Methodsin Applied Mechanics and Engineering, 1992, 94(3): 339-351.

[16] 陳涵, 余莉, 李少騰, 等. 折疊降落傘展開(kāi)過(guò)程研究[J]. 航天返回與遙感, 2012, 33(2): 1-6. CHEN Han, YU Li, LI Shaoteng, et al. A Study on the Opening Processs of Folded Parachute[J]. Spacecraft Recovery & Remote Sensing, 2012, 33(2): 1-6.(in Chinese)

[17] 陳猛, 王璐, 陳涵, 等. 降落傘充氣過(guò)程流固耦合方法的數(shù)值預(yù)測(cè)分析[J]. 南京航空航天大學(xué)學(xué)報(bào), 2013, 45(4): 515-520. CHEN Meng, WANG Lu, CHENG Han, et al. Numerical Prediction Analysis of Parachute Inflation Process Using Fluid-structure Interaction Method[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2013, 45(4): 515-520. (in Chinese)

[18] KIM Y, PESKIN C S. 3-D Parachute Simulation by the Immersed Boundary Method[J]. Computers & Fluid, 2009, 38(6): 1080-1090.

[19] MIYOSHI M, ISHII T, HASHIMOTO A, et al. Computational Analysis of Parachute Motion Using the Immersed Boundary Method[C]//46th AIAA Aerospace Science Meeting and Exhibit, AIAA Paper 2008-311, Reno, Nevada, 2008.

[20] XUE X, KOYAMA H, NAKAMURA Y, et al. Effects of Suspension Line on Flow Field around a Supersonic Parachute[J]. Aerospace Science and Technology, 2015, 43: 63-70.

[21] XUE X, KOYAMA H, NAKAMURA Y. Numerical Simulation of Supersonic Aerodynamic Interaction of Aparachute System[J]. Trans. JSASS Aerospace Tech. Japan, 2013, 11: 33-42.

[22] NISHIYAMA Y. Aerodynamic Characteristics of the Supersonic Parachute with Its Opening Process[D]. Nagoya: Master Thesis of Nagoya University, 2013.

[23] SHIMA E, JOUNOUCHI T. Roe of CFD in Aeronautical Engineering (No.14) -AUSM type Upwind Schemes-[C]// Proceedings of 14th NAL Symposium on Aircraft Computational Aerodynamics, NAL-SP-34, Special publication of National Aerospace Laboratory, 1997: 7-12.

[24] OCHI A, NAKAMURA Y. A Development of Aerodynamics Analysis Tool Using Cartesian System (First Report)[C]//The Japan Society of Fluid Mechanics 19th CFD symposium, 2005.

[25] PESKIN C S. Flow Patterns Around Heart Valves: a Numerical Method [J]. Journal of Computational Physics, 1972, 10(2): 252-271.

[26] BAMHARDT M, DRAYNAT T, et al. Detached Eddy Simulations of the MSL Parachute at Supersonic Conditions[C]//19th AIAA Aerodynamic Decelerator Systems Technology Conference and Seminar, AIAA Paper 2007-2529, Williamsburg, VA, 2007.

[27] El-GENDI M. Effect of Trailing Edge Geometry on the Flow Behavior through Rectilinear Turbine Cascades[D]. Nagoya: Doctor Thesis of Nagoya University, 2010.

[28] BENNEY R J, STEIN K R. Computation Fluid-structure Interaction Model for Parachute Inflation[J]. Journal of Aircraft 1996, 33(4): 730-736.

Numerical Simulation of Aerodynamic Interaction of Supersonic Parachute System

XUE Xiaopeng1WEN Chihyung2

(1 School of Aeronautics and Astronautics, Central South University, Changsha 410084, China)(2 Shenzhen Research Institute, The Hong Kong Polytechnic University, Shenzhen 518057, China)

In the present study, the supersonic flow over a three-dimensional flexible parachute system is numerically simulated by using a simple “immersed boundary technique” together with the fluid-structure coupling scheme. The parachute system employed here consists of a capsule and a canopy. The objective of this study is to investigate the effects of aerodynamic interaction on the dynamics of the canopy behavior, to clarify the performance of the supersonic parachute system in terms of design parameters such as trailing distance between the capsule and canopy, and furthermore to confirm whether “Detached Eddy Simulation”（DES）method can be applied to solve the supersonic flexible three-dimensional parachute problem. As a result, as the trailing distance is rather small, the canopy undergoes severe shrink. As this distance increases, the area oscillation of the canopy shape is observed, leading to a larger drag coefficient. Moreover, when Detached Eddy Simulation method is applied in the case with larger trailing distance, the capsule turbulent wake interacts with the canopy shock more strongly, leading to a larger drag coefficient.

aerodynamic interaction; detached eddy simulation; fluid-structure interaction; supersonic parachute

（編輯：陳艷霞）

V445

A

1009-8518(2016)03-0009-10

10.3969/j.issn.1009-8518.2016.03.002

薛曉鵬，男，1982年生，2009年獲日本名古屋大學(xué)航空航天工程專(zhuān)業(yè)博士學(xué)位，現(xiàn)任中南大學(xué)航空航天學(xué)院講師。研究方向?yàn)闅鈩?dòng)力減速器技術(shù)數(shù)值計(jì)算方法。E-mail：xuexiaopeng@csu.edu.cn。

溫志湧，男，1964年生，1994年獲美國(guó)加州理工學(xué)院航天博士學(xué)位，現(xiàn)任香港理工大學(xué)機(jī)械工程系深圳研究院教授。主要研究方向?yàn)楦叱諝鈩?dòng)力學(xué)，實(shí)驗(yàn)流體力學(xué)，智能流體。

2015-09-23