• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    金屬顆粒燃燒過(guò)程表面曳力變化的數(shù)值研究

    2015-04-22 05:52:08劉叢林
    固體火箭技術(shù) 2015年4期
    關(guān)鍵詞:曳力當(dāng)量流場(chǎng)

    賀 征, 劉叢林,李 卓,顧 璇,郜 冶

    (哈爾濱工程大學(xué) 航天與建筑工程學(xué)院,哈爾濱 150001)

    ?

    金屬顆粒燃燒過(guò)程表面曳力變化的數(shù)值研究

    賀 征, 劉叢林,李 卓,顧 璇,郜 冶

    (哈爾濱工程大學(xué) 航天與建筑工程學(xué)院,哈爾濱 150001)

    為了考察固體火箭發(fā)動(dòng)機(jī)內(nèi)金屬顆粒在生長(zhǎng)過(guò)程中所產(chǎn)生的非球顆粒受力問(wèn)題,采用數(shù)值模擬方法,在通過(guò)和文獻(xiàn)及實(shí)驗(yàn)數(shù)據(jù)對(duì)比驗(yàn)證確定最佳計(jì)算模型的基礎(chǔ)上,對(duì)初始半徑為100 μm的鋁顆粒處于不同相變?nèi)紵A段時(shí)所受到的曳力進(jìn)行對(duì)比分析。結(jié)果表明,在多相流場(chǎng)中,非球顆粒表面可能存在兩處壓力為零的點(diǎn),曳力系數(shù)普遍大于其當(dāng)量直徑顆粒的計(jì)算結(jié)果,當(dāng)顆粒外形嚴(yán)重偏離球體時(shí),經(jīng)驗(yàn)公式計(jì)算所得的曳力系數(shù)有失精準(zhǔn),需進(jìn)行相應(yīng)修正。

    金屬顆粒;非球顆粒;燃燒;當(dāng)量直徑;曳力系數(shù)

    0 引言

    含金屬顆粒的氣固多相流是一種十分典型的多相流動(dòng)。其中,分散相顆粒的動(dòng)力學(xué)特性是工程應(yīng)用中不可忽視的重要問(wèn)題。當(dāng)金屬顆粒在流場(chǎng)中發(fā)生相變?nèi)紵龝r(shí),部分凝相將回落到顆粒表面,令其轉(zhuǎn)變?yōu)椴灰?guī)則的非球體[1-2],導(dǎo)致動(dòng)力學(xué)特性發(fā)生一系列改變,進(jìn)而影響整個(gè)多相流場(chǎng)的流動(dòng)特性。而經(jīng)典的顆粒受力理論是以規(guī)則球形為基礎(chǔ)進(jìn)行推導(dǎo)的,在實(shí)際應(yīng)用中,不可避免地會(huì)產(chǎn)生一定偏差[3]。采用數(shù)值模擬方法,針對(duì)固體火箭發(fā)動(dòng)機(jī)內(nèi),金屬顆粒在生長(zhǎng)過(guò)程中,因形狀變化而引起受力改變的現(xiàn)象進(jìn)行研究,對(duì)工程應(yīng)用有一定參考作用。

    早期的研究集中于球形液滴,并已經(jīng)形成了成熟的經(jīng)驗(yàn)關(guān)系式[4-5]。相對(duì)而言,對(duì)于非球形顆粒的研究開(kāi)展的較少,Komar[6]、Loth[7]開(kāi)展了圓柱體的阻力系數(shù)計(jì)算,并推導(dǎo)其阻力系數(shù)的變化。Sabine[8]等研究了若干顆粒結(jié)合時(shí)的阻力系數(shù)。Youngho and Changhoon[9]模擬計(jì)算了變形顆粒的的阻力和升力,發(fā)現(xiàn)其數(shù)值與顆粒的形變有直接關(guān)系。

    目前,對(duì)于非球形不規(guī)則顆粒的數(shù)值計(jì)算,大多采用近似球形的處理方法[10]。這種近似過(guò)程的實(shí)現(xiàn)主要通過(guò)2個(gè)途徑進(jìn)行——引入等效直徑或引入球形度。其中,前者應(yīng)用的較為廣泛,最常用的等效直徑是等體積直徑dv[11]。引入dv后,認(rèn)為非球顆粒與其等效直徑相同的球形顆粒受力相同,文中對(duì)這2種計(jì)算結(jié)果進(jìn)行了對(duì)比分析。

    1 物理模型

    鋁顆粒初始進(jìn)入發(fā)動(dòng)機(jī)時(shí),其表面將迅速被氧化物所包裹,顆粒形狀已經(jīng)發(fā)生變化[12],顆粒表面積以及同來(lái)流方向相垂直的橫截面面積與球形顆粒不再相同,受力狀態(tài)發(fā)生一定改變。

    Merrill[13]研究了鋁顆粒在固體火箭發(fā)動(dòng)機(jī)環(huán)境中的燃燒過(guò)程,對(duì)Al2O3的生成過(guò)程做了仔細(xì)的推導(dǎo)和計(jì)算。相比于傳統(tǒng)的R2定律,Merrill所得結(jié)果更加符合實(shí)際情況。以Merrill的研究結(jié)果為基礎(chǔ),對(duì)初始半徑為100 μm的鋁顆粒在熱燃?xì)饬髦刑幱诓煌嘧內(nèi)紵A段時(shí)的受力問(wèn)題進(jìn)行了分析,計(jì)算工況列于表1中。

    表1 計(jì)算工況匯總Table1 Summary of computation cases

    其中,t表示顆粒進(jìn)入發(fā)動(dòng)機(jī)燃燒室中的時(shí)間,RAl、Rox分別表示顆粒中鋁的半徑和與其相連的球形氧化物半徑,Req為其相應(yīng)的等體積半徑,Ap、Apeq分別代表真實(shí)顆粒與等體積顆粒垂直于來(lái)流方向的橫截面面積。統(tǒng)計(jì)表明,隨顆粒燃燒狀態(tài)的改變,等體積顆粒的橫截面面積偏離真實(shí)顆粒的水平(Δ%)較大,最大相差49.9%,即便在顆粒初始進(jìn)入流場(chǎng)時(shí),也有1.1%的差別。

    利用Fluent軟件,對(duì)流場(chǎng)中顆粒的受力狀態(tài)進(jìn)行分析。顆粒半徑r取為100 μm,為準(zhǔn)確計(jì)算顆粒附近流場(chǎng)的流動(dòng)狀態(tài),同時(shí)兼顧計(jì)算網(wǎng)格的數(shù)量,計(jì)算區(qū)域中,取平行于來(lái)流方向的長(zhǎng)l為30r,即3 mm,垂直于來(lái)流方向的寬d為20r,即2 mm。顆粒置于流場(chǎng)中,靠近來(lái)流方向,中心與入口邊界相距l(xiāng)1為10r,即1 mm。采用非正交網(wǎng)格進(jìn)行劃分,顆粒的近壁面處做加密處理,最小網(wǎng)格尺寸為5 μm,最大為100 μm,總網(wǎng)格數(shù)約為22萬(wàn),如圖1所示。

    2 控制方程

    假設(shè)來(lái)流是穩(wěn)態(tài)的,其物理性質(zhì)在運(yùn)動(dòng)過(guò)程中不發(fā)生變化,計(jì)算的控制方程為

    ▽·U=0

    (1)

    ▽·(ρUU)=-▽p+μ▽2U

    (2)

    為簡(jiǎn)化模擬條件,不考慮氣固兩相間的傳熱,取空氣為介質(zhì),密度為1.225kg/m3,粘性系數(shù)為1.789×10-5kg/(m·s),流場(chǎng)入口速度取發(fā)動(dòng)機(jī)燃燒室環(huán)境下的均值,為20m/s,流場(chǎng)出口為自由邊界。

    顆粒曳力系數(shù)Cd由顆粒所受到的阻力FD來(lái)定義[14]:

    (3)

    顆粒所受的阻力FD由2部分組成:顆粒表面壓強(qiáng)梯度所產(chǎn)生的壓力阻力,以及流體粘性所產(chǎn)生的粘性阻力。因此,顆粒的曳力系數(shù)也分為2部分——壓力曳力系數(shù)和粘性曳力系數(shù)。

    隨顆粒相變?nèi)紵倪M(jìn)行,尺寸不斷改變,顆粒雷諾數(shù)Rep隨之發(fā)生變化,以當(dāng)量直徑計(jì),對(duì)應(yīng)于t=0、30、60、90 ms時(shí),Rep分別為274.0、116.2、10.4、9.6。

    (a)三維計(jì)算模型

    (b)模型網(wǎng)格分布

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

    3.1 模型驗(yàn)證

    為驗(yàn)證各種湍流模型對(duì)計(jì)算結(jié)果的影響,首先對(duì)顆粒在1≤Rep≤1 000范圍內(nèi)的受力問(wèn)題進(jìn)行計(jì)算,以尋求最佳計(jì)算模型。采用的湍流模型包括realizablek-ε模型、RNGk-ε模型、standardk-ε模型和Spalart-Allmadas模型。

    驗(yàn)證工作共包括32種工況,顆粒雷諾數(shù)Rep分別取為1、10、30、50、100、300、500及1 000,流場(chǎng)工質(zhì)均為空氣。文獻(xiàn)[15]與[16]分別對(duì)Rep在18.7~87.3以及10~809.7范圍內(nèi)的顆粒曳力系數(shù)做了實(shí)驗(yàn),文獻(xiàn)[17]采用直接模型方法對(duì)Rep在10~300范圍內(nèi)的顆粒曳力系數(shù)進(jìn)行了計(jì)算,以上結(jié)果均可作為參考。對(duì)于Rep=1和1 000的工況,尚無(wú)實(shí)驗(yàn)可詢(xún),取標(biāo)準(zhǔn)阻力曲線(xiàn)值相對(duì)比。

    圖2直觀(guān)地反映了各模型的計(jì)算結(jié)果與標(biāo)準(zhǔn)阻力曲線(xiàn)和文獻(xiàn)實(shí)驗(yàn)值的對(duì)比。當(dāng)顆粒雷諾數(shù)Rep<100時(shí),除Spalart-Allmadas模型外,其余3種模型均能很好地與標(biāo)準(zhǔn)阻力曲線(xiàn)相符。其中,realizablek-ε模型下各種情況的計(jì)算值與標(biāo)準(zhǔn)阻力曲線(xiàn)值的均方差為3.63,RNGk-ε模型的均方差為3.77、標(biāo)準(zhǔn)k-ε模型為4.05,而Spalart-Allmadas模型的均方差則高達(dá)14.28。當(dāng)顆粒雷諾數(shù)Rep較高(>100)時(shí),4種模型的計(jì)算值與標(biāo)準(zhǔn)曲線(xiàn)相比,差別都很大,但可較好符合文獻(xiàn)[2-3]的實(shí)驗(yàn)結(jié)果。為兼顧較大范圍內(nèi)的計(jì)算,選擇Realizablek-ε模型進(jìn)行顆粒受力分析的數(shù)值模擬更為適宜。

    圖2 不同湍流模型的計(jì)算結(jié)果與文獻(xiàn)結(jié)果的對(duì)比Fig.2 Computed results of different models vsreference results

    3.2 顆粒周?chē)鲌?chǎng)壓力與速度分布對(duì)比分析

    在多相流場(chǎng)中,顆粒會(huì)隨著流動(dòng)發(fā)生位置偏轉(zhuǎn),為分析真實(shí)顆粒與球形當(dāng)量徑顆粒的差別,取典型狀態(tài)——顆粒在垂直來(lái)流方向投影面最大時(shí)進(jìn)行計(jì)算。不同時(shí)刻下,真實(shí)顆粒與當(dāng)量直徑顆粒流場(chǎng)附近的壓力與速度分布如圖3、圖4所示。在迎風(fēng)方向上,顆粒表面附近的氣流速度由流場(chǎng)入口沿顆粒中心軸線(xiàn)方向迅速降低,在顆粒表面處出現(xiàn)駐點(diǎn),同時(shí)當(dāng)?shù)貕毫_(dá)到最高值。顆粒背風(fēng)面處也存在著一個(gè)速度為零的區(qū)域,同時(shí)當(dāng)?shù)貕毫档阶畹椭?。流?chǎng)中顆粒背風(fēng)域形成長(zhǎng)長(zhǎng)的流動(dòng)尾跡表明,迎風(fēng)面的氣流速度衰減梯度明顯大于背風(fēng)面的速度增長(zhǎng)梯度。

    觀(guān)察顆粒表面的氣流變化可發(fā)現(xiàn),氣流由駐點(diǎn)的零速度向滯止點(diǎn)的零速度轉(zhuǎn)變過(guò)程是氣流速度沿顆粒表面先增加,至顆粒與來(lái)流相垂直的軸線(xiàn)中心處,速度達(dá)到最大值,且此值大于流場(chǎng)的平均速度;之后,因顆粒背風(fēng)面存在渦流區(qū),氣流速度逐漸降低,至顆粒平行于來(lái)流的軸線(xiàn)中心處,再次降為0。在所計(jì)算的4種雷諾數(shù)下,顆粒迎風(fēng)面與背風(fēng)面的壓差達(dá)550 Pa左右,顆粒表面處的最大速度值比流場(chǎng)平均速度高10%。隨著顆粒在流場(chǎng)中相變?nèi)紵倪M(jìn)行,當(dāng)顆粒尺寸減小時(shí),顆粒背風(fēng)域所形成的流動(dòng)尾跡隨之變短,同時(shí)顆粒附近的壓力變化區(qū)(以-50~50 Pa計(jì))半徑減小。

    (a)當(dāng)量直徑顆粒

    (b)真實(shí)顆粒

    對(duì)比圖3、圖4中當(dāng)量直徑顆粒流場(chǎng)與真實(shí)顆粒流場(chǎng)流動(dòng)情況可知,當(dāng)t=0 ms,顆粒剛剛開(kāi)始燃燒時(shí),真實(shí)顆粒軸向長(zhǎng)度大于當(dāng)量直徑下的顆粒,顆粒背風(fēng)域所形成的流動(dòng)尾跡較長(zhǎng)。但此后,因當(dāng)量直徑下顆粒背風(fēng)域所引起的渦流較強(qiáng),氣流需要較長(zhǎng)的路徑來(lái)恢復(fù)原來(lái)的速度。所以,其余3種工況下,真實(shí)顆粒背風(fēng)域所形成的流動(dòng)尾跡相對(duì)較短。因?yàn)轭w粒在相變過(guò)程中呈現(xiàn)出不規(guī)則的幾何形體,真實(shí)顆粒表面附近的壓力分布出現(xiàn)較大波動(dòng),尤其在t=60 ms和t=90 ms時(shí),顆粒外形明顯分成兩個(gè)區(qū)域。因此,在顆粒的表面處存在兩處壓力為0的點(diǎn),這與當(dāng)量直徑下顆粒表面處的壓力分布有很大區(qū)別。

    (a)當(dāng)量直徑顆粒

    (b)真實(shí)顆粒

    3.3 顆粒曳力對(duì)比分析

    圖5反映了鋁顆粒在燃燒過(guò)程中其曳力系數(shù)的變化規(guī)律,以Cd表示;同時(shí),也計(jì)算了當(dāng)量直徑下顆粒曳力系數(shù)變化,以Cdeq表示。由理論分析知,隨著燃燒的進(jìn)行,顆粒直徑越來(lái)越小,將導(dǎo)致顆粒雷諾數(shù)不斷降低,由此也會(huì)使得曳力系數(shù)不斷升高。應(yīng)用當(dāng)量直徑計(jì)算所得的結(jié)果也符合這一規(guī)律,但以真實(shí)顆粒計(jì)算時(shí),二者卻存在一定差別,整體而言,真實(shí)顆粒的曳力系數(shù)普遍大于當(dāng)量直徑下的顆粒。圖5(a)表明,隨顆粒相變?nèi)紵倪M(jìn)行,真實(shí)顆粒表面的壓力曳力系數(shù)發(fā)生很大變化,且無(wú)規(guī)律可循。因?yàn)轭w粒的形狀在燃燒過(guò)程中不斷發(fā)生不規(guī)則的改變,所以其表面壓力分布狀況也表現(xiàn)出很大的波動(dòng)。圖5(b)表明,真實(shí)顆粒表面的粘性曳力系數(shù)一直大于當(dāng)量直徑下的顆粒,尤其當(dāng)燃燒進(jìn)行了60 ms后,以當(dāng)量直徑計(jì)算所得到的粘性曳力系數(shù)僅為真實(shí)顆粒的68.9%,當(dāng)顆粒終止燃燒時(shí),前者仍比后者高18%。圖5(c)為顆粒的總曳力系數(shù)對(duì)比,其變化趨勢(shì)與圖5(b)相似。

    從所計(jì)算的4種工況而言,以當(dāng)量直徑計(jì)算所得的壓力曳力系數(shù)略高于真實(shí)顆粒,而粘性曳力系數(shù)則低于真實(shí)顆粒。所以,總曳力系數(shù)的偏差并不明顯。但在顆粒燃燒進(jìn)行到60 ms時(shí),因外形已嚴(yán)重偏離球體,顆粒表面壓力曳力系數(shù)與粘性曳力系數(shù)均大于當(dāng)量直徑下的計(jì)算值,其總曳力系數(shù)與經(jīng)驗(yàn)公式計(jì)算所得的結(jié)果大不相同。此時(shí),標(biāo)準(zhǔn)阻力曲線(xiàn)已不再適用,需對(duì)其進(jìn)行相應(yīng)的修正。

    (a)壓力曳力系數(shù)對(duì)比

    (b)粘性曳力系數(shù)對(duì)比

    (c)總曳力系數(shù)對(duì)比

    4 結(jié)論

    因?yàn)榻饘兕w粒在相變過(guò)程中呈現(xiàn)出不規(guī)則的幾何形體,真實(shí)顆粒表面附近的壓力分布與周?chē)鲌?chǎng)速度分布同球形顆粒有較大區(qū)別,其曳力系數(shù)亦有較大的變化。

    (1)隨著顆粒在流場(chǎng)中相變?nèi)紵倪M(jìn)行,當(dāng)顆粒尺寸減小時(shí),顆粒背風(fēng)域所形成的流動(dòng)尾跡隨之變短;同時(shí),顆粒附近的壓力變化區(qū)半徑減小。

    (2)在真實(shí)非球體顆粒的表面處,可能存在兩處壓力為0的點(diǎn)。顆粒直徑的減小將導(dǎo)致其雷諾數(shù)降低,曳力系數(shù)不斷升高。整體而言,真實(shí)顆粒的曳力系數(shù)普遍大于其當(dāng)量直徑顆粒。其中,壓力曳力系數(shù)變化很大,但無(wú)規(guī)律可循。

    (3)粘性曳力系數(shù)普遍大于當(dāng)量直徑下的顆粒,當(dāng)燃燒進(jìn)行60 ms后,其值為以當(dāng)量直徑計(jì)算所得結(jié)果的1.46倍。當(dāng)顆粒終止燃燒時(shí),前者仍比后者高出18%。

    (4)當(dāng)顆粒外形嚴(yán)重偏離球體時(shí),經(jīng)驗(yàn)公式計(jì)算所得的阻力系數(shù)有失精準(zhǔn)。此時(shí),標(biāo)準(zhǔn)阻力曲線(xiàn)已不再適用,可考慮對(duì)當(dāng)量直徑的計(jì)算結(jié)果進(jìn)行系數(shù)修正。

    [1] 方丁酉.兩相流體力學(xué)[D].長(zhǎng)沙:國(guó)防科技大學(xué)出版社,1988:55-58.

    [2] Olsen S E,Beckstead M W.Burn time measurements of single aluminum particles in steam and CO2mixtures[J].Journal of Propulsion and Power,1996,12(4):662-671.

    [3] 林建中.超常顆粒多相流體動(dòng)力學(xué)——圓柱狀顆粒兩相流[M].北京:科學(xué)出版社,2008:85-86.

    [4] Khan A R,Richardson J F.The resistance to motion of a solid sphere in a fluid[J].Chem.Eng.Commun,1987,62:135-150.

    [5] Hartman M,Yates J G.Free-fall of solid particles through fluids[J].Collet.Czechoslov.Chem.Commun,1993,58(5):961-982.

    [6] Komar P D.Settling velocities of circular cylinders at low Reynolds number[J].J.Geol.,1980,48:327-328.

    [7] Loth E.Drag of non-spherical solid particles of regular and irregular shape[J].Powder Technology,2008,182:342-353.

    [8] Sabine Tran-cng,Michael Gay,efastathios E Michaelides.Drag coefficients of irregular shaped particles[J].Powder Technology,2004,139:21-32.

    [9] Youngho Suh,Changhoon Lee.A numerical method for the calculation of drag and lift of a deformable droplet in shear flow[J].Journal of Computational Physics,2013,241:35-57.

    [10] Lu Hui-lin,Liu Wen-tie,Zhao Guang-bo.Computational modeling of dense gas particle flow in a pipe:kinetic theory approach of granular flow[J].Journal of Chemical Industry and Engineering (China),2005,5(1):31-38.

    [11] Gan Lin,Xu Mao-sheng,zhu Bing-chen.Heat transfer parameters of packed bed with ring pellet[J].Journal of Chemical Industry and Engineering (China),2000,5(6):778-783.

    [12] Robert Geisler.A global view of the use of aluminum fuel in solid rocket motors[R].AIAA 2002-3748.

    [13] Merrill K King.Aluminum combustion in a solid rocket motor environment[J].Proceedings of The Combustion Institute,2009,32(2):2107-2114.

    [14] Ounis H,Ahmadi G,McLaughlin J B.Brownian diffusion of submicrometer particles in the viscous sublayer[J].Journal of Colloid and Interface Science,1991,143(1):266-277.

    [15] Unnikrisham A,Chhabra R P.An experimental study of motion of cylinders in Newtonian fluids:wall effects and drag coefficient[J].Can.J.Chem.Eng.,1991,9(9):729-735.

    [16] Warnica W D,Renksizbulut M,Strong A B.Drag coefficients of spherical liquid droplets,Part 2:Turbulent gaseous fields[J].Experiments in Fluids.,1995,18:265-276.

    [17] By Prosenjit Bagchi,Balachandar S.Inertial and viscous forces on a rigid sphere in straining flows at moderate Reynolds numbers[J].J.Fluid Mech.,2003,481:105-148.

    (編輯:崔賢彬)

    Numerical study on the change of surface drag force during metal particle combustion

    HE Zheng,LIU Cong-lin,LI Zhuo,GU Xuan,GAO Ye
    (College of Aerospace and Civil Engineering, Heilongjiang, Harbin 150001, China)

    To study the force acted on the non-spherical particle which is produced by burning metal particle in solid ramjet motor,numerical simulation was employed to compare and analyze the drag of aluminum particle whose initial radius is 100 μm in variable phases of the burning process. Before that,the best computational model was determined and validated by comparing with the reference and experimental data.The results show that,in multiphase flow there are two points where the pressure is zero at the non-spherical particle's surface.Compared to the equivalent diameter particle, the real particle's drag coefficient is generally bigger.If the particle is not close to the spherical any more, the drag coefficient computed by the empirical formula is of great roughness,and the formula should be consequentially amended.

    metal particle;non-spherical particle;combustion;equivalent diameter;drag coefficient

    2014-07-04;

    :2014-08-24。

    國(guó)家自然科學(xué)基金(11372079);中央高校基本科研業(yè)務(wù)費(fèi)專(zhuān)項(xiàng)基金(HEUCF130203)。

    賀征(1978—),男,博士/副教授,研究方向?yàn)榘l(fā)動(dòng)機(jī)內(nèi)燃燒與多相流動(dòng)。E-mail:hezheng1978@163.com

    劉叢林(1981—),女,博士。E-mail:383523445@qq.com

    V435

    A

    1006-2793(2015)04-0492-05

    10.7673/j.issn.1006-2793.2015.04.008

    猜你喜歡
    曳力當(dāng)量流場(chǎng)
    預(yù)測(cè)天然氣斜井臨界攜液流量新方法
    循環(huán)流化床鍋爐爐膛流動(dòng)特性數(shù)值模擬進(jìn)展
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場(chǎng)計(jì)算
    轉(zhuǎn)杯紡排雜區(qū)流場(chǎng)與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場(chǎng)統(tǒng)計(jì)分析
    基于EMMS模型的攪拌釜內(nèi)氣液兩相流數(shù)值模擬
    黃河之聲(2016年24期)2016-02-03 09:01:52
    基于瞬態(tài)流場(chǎng)計(jì)算的滑動(dòng)軸承靜平衡位置求解
    超壓測(cè)試方法對(duì)炸藥TNT當(dāng)量計(jì)算結(jié)果的影響
    環(huán)空附加當(dāng)量循環(huán)密度的計(jì)算方法
    斷塊油氣田(2014年5期)2014-03-11 15:33:50
    黄色女人牲交| h日本视频在线播放| 91麻豆av在线| 宅男免费午夜| 亚洲专区字幕在线| 18美女黄网站色大片免费观看| 美女黄网站色视频| 国产一级毛片七仙女欲春2| 窝窝影院91人妻| av片东京热男人的天堂| 51午夜福利影视在线观看| 男女做爰动态图高潮gif福利片| 国产精品1区2区在线观看.| 黑人欧美特级aaaaaa片| 欧美色视频一区免费| 国产成人精品久久二区二区91| 嫩草影院入口| 久久久久久久精品吃奶| 日韩 欧美 亚洲 中文字幕| 超碰成人久久| 亚洲欧洲精品一区二区精品久久久| 久久九九热精品免费| 国产av不卡久久| 99国产极品粉嫩在线观看| 亚洲无线观看免费| 成年女人毛片免费观看观看9| 91久久精品国产一区二区成人 | 丰满人妻熟妇乱又伦精品不卡| 午夜精品久久久久久毛片777| 制服丝袜大香蕉在线| 香蕉av资源在线| 波多野结衣高清无吗| 亚洲av五月六月丁香网| 91老司机精品| 两人在一起打扑克的视频| 欧美黑人巨大hd| 亚洲中文av在线| 啦啦啦免费观看视频1| 国产一级毛片七仙女欲春2| 国产极品精品免费视频能看的| 欧美国产日韩亚洲一区| 国产真实乱freesex| 亚洲在线观看片| 老司机福利观看| 观看美女的网站| 热99在线观看视频| 99久久99久久久精品蜜桃| 嫩草影院入口| 97超视频在线观看视频| av在线蜜桃| 国产激情偷乱视频一区二区| xxxwww97欧美| 狠狠狠狠99中文字幕| 久久精品国产亚洲av香蕉五月| 动漫黄色视频在线观看| 麻豆一二三区av精品| 成人鲁丝片一二三区免费| 最近最新免费中文字幕在线| 天堂动漫精品| 亚洲精品在线观看二区| 757午夜福利合集在线观看| 欧洲精品卡2卡3卡4卡5卡区| 俺也久久电影网| 久久久久久久午夜电影| 美女大奶头视频| 亚洲成人免费电影在线观看| 真人做人爱边吃奶动态| 国产v大片淫在线免费观看| 国产精品 国内视频| 亚洲激情在线av| 999久久久国产精品视频| 91九色精品人成在线观看| 99久久综合精品五月天人人| 最近最新免费中文字幕在线| 国产成年人精品一区二区| 国产极品精品免费视频能看的| 国产成人精品久久二区二区免费| 国产精品亚洲av一区麻豆| 免费搜索国产男女视频| 欧美性猛交╳xxx乱大交人| 18禁黄网站禁片免费观看直播| 九九在线视频观看精品| 美女被艹到高潮喷水动态| 成年女人毛片免费观看观看9| 18禁黄网站禁片午夜丰满| 亚洲精品中文字幕一二三四区| 99热只有精品国产| 99国产极品粉嫩在线观看| 老司机在亚洲福利影院| 国产探花在线观看一区二区| 老司机福利观看| 精品一区二区三区视频在线 | 日本五十路高清| 女同久久另类99精品国产91| 亚洲中文日韩欧美视频| 色尼玛亚洲综合影院| 精品一区二区三区四区五区乱码| 嫩草影院入口| 99久久国产精品久久久| 十八禁网站免费在线| 午夜福利在线在线| 欧美日韩综合久久久久久 | 精品免费久久久久久久清纯| 国产黄a三级三级三级人| 免费av毛片视频| 丁香欧美五月| av中文乱码字幕在线| 久久欧美精品欧美久久欧美| 亚洲熟女毛片儿| 免费看光身美女| 美女被艹到高潮喷水动态| 蜜桃久久精品国产亚洲av| 亚洲精华国产精华精| 最近视频中文字幕2019在线8| 国产精品,欧美在线| 美女大奶头视频| 亚洲专区中文字幕在线| 99久久久亚洲精品蜜臀av| 国产99白浆流出| 成年版毛片免费区| 久久久久免费精品人妻一区二区| 天堂影院成人在线观看| 免费看日本二区| 欧美日韩一级在线毛片| 色精品久久人妻99蜜桃| 亚洲无线在线观看| 国产伦精品一区二区三区视频9 | 日本三级黄在线观看| av片东京热男人的天堂| 在线免费观看不下载黄p国产 | 夜夜爽天天搞| 一夜夜www| 国产精品av久久久久免费| 国产三级黄色录像| 国产av不卡久久| 十八禁网站免费在线| 精品乱码久久久久久99久播| 久久久成人免费电影| 欧美在线黄色| 少妇人妻一区二区三区视频| 婷婷精品国产亚洲av在线| 成人av一区二区三区在线看| 国产亚洲精品久久久久久毛片| 老熟妇乱子伦视频在线观看| 日韩国内少妇激情av| 欧美日韩国产亚洲二区| 国产精品一及| 美女黄网站色视频| 岛国视频午夜一区免费看| 欧美午夜高清在线| 国产精品久久久av美女十八| 99国产极品粉嫩在线观看| 日韩人妻高清精品专区| 国内精品美女久久久久久| 中文亚洲av片在线观看爽| 久久国产精品影院| 一个人看的www免费观看视频| 欧美日韩瑟瑟在线播放| 久久久久国产精品人妻aⅴ院| 无限看片的www在线观看| 成人精品一区二区免费| 国产伦精品一区二区三区四那| 99热这里只有是精品50| xxx96com| 18禁观看日本| 国产精品久久久久久久电影 | 天堂网av新在线| 日本一本二区三区精品| 中文在线观看免费www的网站| 成熟少妇高潮喷水视频| 亚洲成a人片在线一区二区| 99在线视频只有这里精品首页| 十八禁人妻一区二区| 18禁美女被吸乳视频| 叶爱在线成人免费视频播放| 久久精品综合一区二区三区| 变态另类丝袜制服| 亚洲18禁久久av| 美女高潮喷水抽搐中文字幕| 国产精品一及| 亚洲激情在线av| 日韩欧美一区二区三区在线观看| 精品99又大又爽又粗少妇毛片 | 两个人视频免费观看高清| 少妇裸体淫交视频免费看高清| 制服丝袜大香蕉在线| 国产黄色小视频在线观看| 搡老妇女老女人老熟妇| 最近最新中文字幕大全免费视频| а√天堂www在线а√下载| 狠狠狠狠99中文字幕| 亚洲欧美日韩卡通动漫| 欧美zozozo另类| cao死你这个sao货| 久久精品91无色码中文字幕| 又大又爽又粗| 久9热在线精品视频| 国产97色在线日韩免费| 夜夜夜夜夜久久久久| 成年免费大片在线观看| 欧美丝袜亚洲另类 | 免费看十八禁软件| 午夜免费激情av| 99在线视频只有这里精品首页| 热99re8久久精品国产| 色视频www国产| 欧美国产日韩亚洲一区| 午夜精品一区二区三区免费看| 欧美日本视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩精品中文字幕看吧| 婷婷精品国产亚洲av| 女同久久另类99精品国产91| 成人特级黄色片久久久久久久| 身体一侧抽搐| 91字幕亚洲| 日本五十路高清| 97人妻精品一区二区三区麻豆| 精品国产三级普通话版| 久久天堂一区二区三区四区| 在线免费观看的www视频| 国产精品一区二区三区四区久久| 琪琪午夜伦伦电影理论片6080| 久久中文字幕一级| 男女视频在线观看网站免费| 午夜免费激情av| 国产麻豆成人av免费视频| 性欧美人与动物交配| 欧洲精品卡2卡3卡4卡5卡区| 一级毛片精品| 日本熟妇午夜| 亚洲成人免费电影在线观看| 亚洲专区中文字幕在线| 一级a爱片免费观看的视频| 亚洲激情在线av| 亚洲国产欧美人成| 精品免费久久久久久久清纯| 国模一区二区三区四区视频 | 人人妻人人看人人澡| 欧美乱色亚洲激情| 国语自产精品视频在线第100页| 全区人妻精品视频| 一进一出抽搐gif免费好疼| 国产精品影院久久| 亚洲自拍偷在线| 日韩 欧美 亚洲 中文字幕| 最新中文字幕久久久久 | 久久久久久久精品吃奶| 亚洲欧美日韩高清专用| 首页视频小说图片口味搜索| 色综合欧美亚洲国产小说| 国产高清videossex| 色尼玛亚洲综合影院| 99热这里只有是精品50| 国产探花在线观看一区二区| 免费一级毛片在线播放高清视频| 亚洲国产欧美网| 黑人巨大精品欧美一区二区mp4| 嫩草影视91久久| 午夜福利欧美成人| 国产av不卡久久| 男女之事视频高清在线观看| 午夜免费激情av| 国产成人av教育| 午夜福利成人在线免费观看| 精品国产乱子伦一区二区三区| 欧美日韩国产亚洲二区| 国产亚洲精品久久久久久毛片| 人人妻,人人澡人人爽秒播| 久久精品国产亚洲av香蕉五月| 亚洲国产欧美人成| 日本 欧美在线| 久久久久性生活片| 伦理电影免费视频| 欧美乱码精品一区二区三区| 亚洲人成网站在线播放欧美日韩| 日韩欧美国产一区二区入口| 精品国产美女av久久久久小说| 久久久精品欧美日韩精品| 十八禁人妻一区二区| 成年免费大片在线观看| 日本一本二区三区精品| 日本熟妇午夜| 午夜a级毛片| 天堂√8在线中文| 色精品久久人妻99蜜桃| 中文字幕精品亚洲无线码一区| 黄色丝袜av网址大全| 久久精品国产综合久久久| 国产97色在线日韩免费| av黄色大香蕉| 99视频精品全部免费 在线 | 一边摸一边抽搐一进一小说| 久久这里只有精品19| 中文在线观看免费www的网站| 欧美日韩国产亚洲二区| 麻豆久久精品国产亚洲av| 一a级毛片在线观看| 欧美性猛交黑人性爽| www.www免费av| 国产成人欧美在线观看| 免费在线观看日本一区| 精品福利观看| 午夜福利18| 亚洲国产欧美人成| 又大又爽又粗| 99久久综合精品五月天人人| 美女被艹到高潮喷水动态| 久久国产精品影院| 亚洲成人中文字幕在线播放| 久久精品综合一区二区三区| a在线观看视频网站| 亚洲第一欧美日韩一区二区三区| 国产激情偷乱视频一区二区| 久久香蕉精品热| 啦啦啦免费观看视频1| 亚洲av电影在线进入| 香蕉国产在线看| 熟女少妇亚洲综合色aaa.| 日日摸夜夜添夜夜添小说| 在线a可以看的网站| 日韩欧美精品v在线| 亚洲熟妇中文字幕五十中出| 欧美激情在线99| 黄色片一级片一级黄色片| 99久久成人亚洲精品观看| 国产精品亚洲美女久久久| 国内精品美女久久久久久| 男人舔女人下体高潮全视频| 国产三级黄色录像| 亚洲精品久久国产高清桃花| 免费看日本二区| 国产伦精品一区二区三区四那| 黑人巨大精品欧美一区二区mp4| 亚洲美女黄片视频| 国产av不卡久久| 成人无遮挡网站| 亚洲第一欧美日韩一区二区三区| 精品国内亚洲2022精品成人| 成年版毛片免费区| 国产男靠女视频免费网站| 日韩欧美三级三区| 日本a在线网址| 成人性生交大片免费视频hd| 久久伊人香网站| 国产亚洲精品av在线| 亚洲中文字幕一区二区三区有码在线看 | 国产精品亚洲一级av第二区| 色噜噜av男人的天堂激情| 制服丝袜大香蕉在线| 免费在线观看成人毛片| 中文字幕av在线有码专区| 久久久久久人人人人人| 亚洲熟妇中文字幕五十中出| 99riav亚洲国产免费| 美女扒开内裤让男人捅视频| 亚洲专区字幕在线| 国产精品久久久久久亚洲av鲁大| 免费高清视频大片| 这个男人来自地球电影免费观看| 成人三级做爰电影| 国产成+人综合+亚洲专区| 男人舔奶头视频| 国产成人福利小说| 毛片女人毛片| 国产又黄又爽又无遮挡在线| 欧美日韩精品网址| 美女免费视频网站| 免费大片18禁| 97超级碰碰碰精品色视频在线观看| 日韩大尺度精品在线看网址| 国产亚洲精品久久久久久毛片| 亚洲片人在线观看| 69av精品久久久久久| 免费看日本二区| 香蕉久久夜色| 日本五十路高清| 嫩草影院精品99| 国内精品久久久久久久电影| 真人做人爱边吃奶动态| 悠悠久久av| 亚洲熟妇熟女久久| 亚洲色图 男人天堂 中文字幕| 久久热在线av| 午夜两性在线视频| 精品国产美女av久久久久小说| 天天躁日日操中文字幕| 国产精品久久电影中文字幕| 免费无遮挡裸体视频| 最新中文字幕久久久久 | 成年人黄色毛片网站| 制服人妻中文乱码| 国产精品98久久久久久宅男小说| 一a级毛片在线观看| 国产高清有码在线观看视频| 免费看日本二区| 国产伦在线观看视频一区| 精品人妻1区二区| 嫩草影院精品99| 97超级碰碰碰精品色视频在线观看| 99精品在免费线老司机午夜| 99热6这里只有精品| 日本一本二区三区精品| 午夜福利在线观看免费完整高清在 | 国产精品一区二区免费欧美| 国产一级毛片七仙女欲春2| 久久久久久久久久黄片| 1024手机看黄色片| 成人精品一区二区免费| www.自偷自拍.com| 日韩欧美免费精品| 校园春色视频在线观看| 欧美日韩综合久久久久久 | 韩国av一区二区三区四区| 国产男靠女视频免费网站| av在线天堂中文字幕| 最新在线观看一区二区三区| 国产亚洲精品一区二区www| 91av网一区二区| 亚洲国产欧美一区二区综合| 最近最新中文字幕大全电影3| 欧美三级亚洲精品| 国产成人福利小说| 免费在线观看成人毛片| 亚洲专区字幕在线| 男人舔女人下体高潮全视频| 母亲3免费完整高清在线观看| 国产高清videossex| 午夜福利在线在线| 亚洲黑人精品在线| 日本一本二区三区精品| 日本黄色视频三级网站网址| а√天堂www在线а√下载| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品青青久久久久久| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区免费欧美| 午夜免费成人在线视频| a级毛片在线看网站| 99久久精品一区二区三区| 国产真人三级小视频在线观看| 国产亚洲精品久久久com| 曰老女人黄片| 香蕉丝袜av| 一区二区三区国产精品乱码| 国产主播在线观看一区二区| 欧美成狂野欧美在线观看| 久久久久久久久中文| 巨乳人妻的诱惑在线观看| 在线观看午夜福利视频| 亚洲 欧美一区二区三区| 日本五十路高清| 色播亚洲综合网| 午夜精品一区二区三区免费看| 久久久久亚洲av毛片大全| 国产精品一及| 桃红色精品国产亚洲av| 99久久精品一区二区三区| 亚洲成人精品中文字幕电影| 亚洲av第一区精品v没综合| 精品久久久久久成人av| 国产视频一区二区在线看| 五月玫瑰六月丁香| 中文字幕精品亚洲无线码一区| 又黄又粗又硬又大视频| 少妇丰满av| 国产欧美日韩一区二区精品| netflix在线观看网站| 一区福利在线观看| 欧美乱妇无乱码| 亚洲午夜理论影院| 男人舔女人下体高潮全视频| 欧美成狂野欧美在线观看| 欧美乱码精品一区二区三区| 热99在线观看视频| 99在线视频只有这里精品首页| 中文字幕高清在线视频| 五月伊人婷婷丁香| 男插女下体视频免费在线播放| 99国产精品99久久久久| 97人妻精品一区二区三区麻豆| 一边摸一边抽搐一进一小说| 国产精品野战在线观看| 国产一区二区在线观看日韩 | 国产乱人视频| 99国产精品一区二区蜜桃av| 淫妇啪啪啪对白视频| av天堂在线播放| 色播亚洲综合网| 小蜜桃在线观看免费完整版高清| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 欧美绝顶高潮抽搐喷水| 亚洲激情在线av| 狂野欧美白嫩少妇大欣赏| 在线观看66精品国产| 综合色av麻豆| 色尼玛亚洲综合影院| 好男人电影高清在线观看| cao死你这个sao货| 国产精品电影一区二区三区| 18禁黄网站禁片免费观看直播| 禁无遮挡网站| 99久久99久久久精品蜜桃| 热99在线观看视频| 女警被强在线播放| 精品福利观看| 成人特级av手机在线观看| 蜜桃久久精品国产亚洲av| ponron亚洲| 日韩av在线大香蕉| 天堂av国产一区二区熟女人妻| 国产精品久久久久久精品电影| 夜夜夜夜夜久久久久| 午夜两性在线视频| 叶爱在线成人免费视频播放| www.熟女人妻精品国产| 日韩欧美国产一区二区入口| 国产精品 国内视频| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| x7x7x7水蜜桃| 日韩中文字幕欧美一区二区| 免费观看人在逋| 99热这里只有是精品50| av国产免费在线观看| 国产真实乱freesex| h日本视频在线播放| 听说在线观看完整版免费高清| 久久午夜综合久久蜜桃| 美女大奶头视频| 欧美日韩福利视频一区二区| 免费av不卡在线播放| 黄色丝袜av网址大全| 国产精品影院久久| 男插女下体视频免费在线播放| 亚洲va日本ⅴa欧美va伊人久久| 99在线人妻在线中文字幕| 国产久久久一区二区三区| 精品人妻1区二区| 欧美乱妇无乱码| 手机成人av网站| 成在线人永久免费视频| 亚洲人成网站高清观看| 村上凉子中文字幕在线| 男女那种视频在线观看| 亚洲精品在线观看二区| 真人一进一出gif抽搐免费| 婷婷精品国产亚洲av| 一个人免费在线观看电影 | 国产一区在线观看成人免费| 精品一区二区三区四区五区乱码| 夜夜看夜夜爽夜夜摸| 亚洲黑人精品在线| 香蕉久久夜色| 欧美性猛交黑人性爽| 精品久久久久久久久久免费视频| 亚洲成人中文字幕在线播放| 99热6这里只有精品| 久久精品91蜜桃| 色播亚洲综合网| 欧美日本视频| 麻豆国产97在线/欧美| 亚洲人与动物交配视频| 在线a可以看的网站| 国产精品久久久久久人妻精品电影| 国产成+人综合+亚洲专区| 视频区欧美日本亚洲| 高潮久久久久久久久久久不卡| 婷婷六月久久综合丁香| 91麻豆av在线| 午夜福利18| 97超级碰碰碰精品色视频在线观看| 午夜视频精品福利| 久久午夜亚洲精品久久| 99视频精品全部免费 在线 | 99精品在免费线老司机午夜| 婷婷六月久久综合丁香| 国产单亲对白刺激| 免费看a级黄色片| bbb黄色大片| 午夜视频精品福利| 国产亚洲精品久久久com| 亚洲精品一卡2卡三卡4卡5卡| 国产高清有码在线观看视频| 国产成+人综合+亚洲专区| or卡值多少钱| 一级毛片高清免费大全| 哪里可以看免费的av片| 91老司机精品| 欧美黑人欧美精品刺激| 国产高清三级在线| 亚洲av美国av| 久久午夜综合久久蜜桃| 日韩欧美一区二区三区在线观看| 国产激情久久老熟女| 99re在线观看精品视频| 久久草成人影院| 少妇人妻一区二区三区视频| 欧美丝袜亚洲另类 | 成人高潮视频无遮挡免费网站| 国产精品自产拍在线观看55亚洲| 小蜜桃在线观看免费完整版高清| 亚洲精品粉嫩美女一区| 亚洲 欧美一区二区三区| 天堂√8在线中文| 国产精品电影一区二区三区| 99在线人妻在线中文字幕| 免费观看精品视频网站| 免费高清视频大片| 人妻丰满熟妇av一区二区三区| 19禁男女啪啪无遮挡网站| 99国产精品一区二区三区| 禁无遮挡网站| 五月伊人婷婷丁香|