• <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
    国产精品久久视频播放| 欧美黄色淫秽网站| 操出白浆在线播放| 老司机午夜福利在线观看视频| 国产人伦9x9x在线观看| 久久 成人 亚洲| 亚洲成av人片在线播放无| 老司机深夜福利视频在线观看| 一级毛片精品| 亚洲成av人片在线播放无| 亚洲va日本ⅴa欧美va伊人久久| 又黄又粗又硬又大视频| 超碰成人久久| 成人高潮视频无遮挡免费网站| 91九色精品人成在线观看| 亚洲午夜理论影院| 国产亚洲精品一区二区www| 成人永久免费在线观看视频| 精品福利观看| 亚洲中文字幕一区二区三区有码在线看 | 国产成人精品久久二区二区91| 在线观看午夜福利视频| 青草久久国产| 成人手机av| 亚洲一卡2卡3卡4卡5卡精品中文| 99热6这里只有精品| 首页视频小说图片口味搜索| 熟女少妇亚洲综合色aaa.| 成人亚洲精品av一区二区| 制服人妻中文乱码| 国产激情欧美一区二区| 9191精品国产免费久久| 成人国产综合亚洲| 国产一级毛片七仙女欲春2| 亚洲电影在线观看av| 久久久久久九九精品二区国产 | 中国美女看黄片| a在线观看视频网站| 欧美成狂野欧美在线观看| 亚洲成人久久爱视频| 在线播放国产精品三级| 脱女人内裤的视频| 亚洲成av人片在线播放无| 级片在线观看| 欧美一级a爱片免费观看看 | 人人妻,人人澡人人爽秒播| 国产成人系列免费观看| 国产精品久久久人人做人人爽| 亚洲18禁久久av| 久久香蕉国产精品| 国产伦人伦偷精品视频| 久久久久久大精品| 亚洲精品中文字幕一二三四区| 午夜福利在线观看吧| 给我免费播放毛片高清在线观看| 午夜激情福利司机影院| 欧洲精品卡2卡3卡4卡5卡区| 国产av麻豆久久久久久久| 特级一级黄色大片| 又紧又爽又黄一区二区| 国产精品综合久久久久久久免费| 无遮挡黄片免费观看| 国产精品久久久久久精品电影| 久久久久久久久久黄片| 国产一区二区三区在线臀色熟女| 露出奶头的视频| 两个人视频免费观看高清| 日韩av在线大香蕉| 日本一二三区视频观看| 日韩有码中文字幕| 可以免费在线观看a视频的电影网站| 特大巨黑吊av在线直播| 高清在线国产一区| 国产99久久九九免费精品| 久久精品国产综合久久久| 亚洲国产精品合色在线| 日本成人三级电影网站| 久久婷婷成人综合色麻豆| 美女免费视频网站| 在线永久观看黄色视频| 国产精品一区二区三区四区久久| 亚洲人成77777在线视频| 一个人观看的视频www高清免费观看 | 黑人操中国人逼视频| 欧美绝顶高潮抽搐喷水| 熟女电影av网| 国产成人精品久久二区二区91| 人人妻人人看人人澡| 50天的宝宝边吃奶边哭怎么回事| 成人欧美大片| 国产精品亚洲美女久久久| 麻豆成人av在线观看| 身体一侧抽搐| 亚洲乱码一区二区免费版| 亚洲精品久久成人aⅴ小说| 国产精品爽爽va在线观看网站| av超薄肉色丝袜交足视频| bbb黄色大片| 日韩欧美国产在线观看| 真人做人爱边吃奶动态| 国产精品99久久99久久久不卡| 午夜免费成人在线视频| 岛国视频午夜一区免费看| 一个人免费在线观看的高清视频| av有码第一页| 丝袜美腿诱惑在线| 中亚洲国语对白在线视频| 香蕉国产在线看| 性欧美人与动物交配| 制服人妻中文乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 无限看片的www在线观看| 久久久久久免费高清国产稀缺| 亚洲成人中文字幕在线播放| 色综合站精品国产| 婷婷丁香在线五月| 欧美成人午夜精品| 亚洲狠狠婷婷综合久久图片| 最新在线观看一区二区三区| 一级黄色大片毛片| 亚洲狠狠婷婷综合久久图片| 亚洲欧洲精品一区二区精品久久久| 精品国产超薄肉色丝袜足j| 19禁男女啪啪无遮挡网站| 99热6这里只有精品| 在线十欧美十亚洲十日本专区| 婷婷六月久久综合丁香| 日本一区二区免费在线视频| 成人国产综合亚洲| 国产精品一区二区三区四区久久| 亚洲第一欧美日韩一区二区三区| 久久久久久免费高清国产稀缺| 久久久久久国产a免费观看| 国产视频内射| 小说图片视频综合网站| 久久久精品欧美日韩精品| 又爽又黄无遮挡网站| 国产精品98久久久久久宅男小说| 国产三级黄色录像| 欧美日韩亚洲国产一区二区在线观看| 老汉色∧v一级毛片| 国产精品日韩av在线免费观看| e午夜精品久久久久久久| 欧美精品啪啪一区二区三区| 亚洲熟妇熟女久久| 久久久久国产一级毛片高清牌| 久热爱精品视频在线9| 国产成年人精品一区二区| 非洲黑人性xxxx精品又粗又长| 又黄又粗又硬又大视频| 天堂动漫精品| 中出人妻视频一区二区| 麻豆久久精品国产亚洲av| 最近视频中文字幕2019在线8| 亚洲av电影不卡..在线观看| 每晚都被弄得嗷嗷叫到高潮| av视频在线观看入口| 男女那种视频在线观看| 午夜影院日韩av| 一个人观看的视频www高清免费观看 | 国产av一区在线观看免费| 国产午夜福利久久久久久| 亚洲av日韩精品久久久久久密| 性欧美人与动物交配| 亚洲精品一卡2卡三卡4卡5卡| www日本黄色视频网| 麻豆久久精品国产亚洲av| 亚洲专区字幕在线| 动漫黄色视频在线观看| 久久久国产欧美日韩av| 男女视频在线观看网站免费 | 日韩精品免费视频一区二区三区| 久9热在线精品视频| 国产成+人综合+亚洲专区| 日韩欧美一区二区三区在线观看| 免费在线观看亚洲国产| 亚洲中文字幕日韩| 在线观看美女被高潮喷水网站 | 一级a爱片免费观看的视频| 1024视频免费在线观看| 欧美乱码精品一区二区三区| 美女高潮喷水抽搐中文字幕| 色综合欧美亚洲国产小说| 亚洲国产高清在线一区二区三| 久久精品综合一区二区三区| 久久午夜综合久久蜜桃| 日韩成人在线观看一区二区三区| 欧美乱妇无乱码| 在线观看舔阴道视频| 久久这里只有精品中国| 18禁黄网站禁片免费观看直播| 大型av网站在线播放| 免费在线观看日本一区| 久久精品夜夜夜夜夜久久蜜豆 | 婷婷精品国产亚洲av| 女生性感内裤真人,穿戴方法视频| 欧美成人免费av一区二区三区| 不卡一级毛片| 国产精品久久视频播放| 搡老熟女国产l中国老女人| 99久久精品热视频| 成年免费大片在线观看| 99久久久亚洲精品蜜臀av| 国产黄色小视频在线观看| 观看免费一级毛片| 亚洲真实伦在线观看| 99riav亚洲国产免费| 午夜日韩欧美国产| 国产高清视频在线观看网站| 亚洲男人的天堂狠狠| 国内少妇人妻偷人精品xxx网站 | 91大片在线观看| 2021天堂中文幕一二区在线观| 欧美成人一区二区免费高清观看 | 日本五十路高清| 这个男人来自地球电影免费观看| netflix在线观看网站| 他把我摸到了高潮在线观看| 精品久久蜜臀av无| 国产精品影院久久| 后天国语完整版免费观看| 亚洲自偷自拍图片 自拍| 波多野结衣高清作品| 精品欧美一区二区三区在线| 又粗又爽又猛毛片免费看| 1024视频免费在线观看| 欧美激情久久久久久爽电影| 淫秽高清视频在线观看| 久久久精品大字幕| 少妇被粗大的猛进出69影院| 久久草成人影院| 久久久久久人人人人人| 日韩高清综合在线| 婷婷亚洲欧美| 欧美日韩中文字幕国产精品一区二区三区| 高清在线国产一区| 一a级毛片在线观看| 黄色片一级片一级黄色片| 亚洲欧美精品综合一区二区三区| 叶爱在线成人免费视频播放| 波多野结衣巨乳人妻| 欧美午夜高清在线| 免费看日本二区| 999精品在线视频| 99久久无色码亚洲精品果冻| 国产伦在线观看视频一区| 免费高清视频大片| 天堂影院成人在线观看| 一a级毛片在线观看| 国产成人一区二区三区免费视频网站| 我的老师免费观看完整版| 欧美黑人巨大hd| 久久久久免费精品人妻一区二区| 久久香蕉国产精品| 国产成人av激情在线播放| 国产不卡一卡二| 国产精品美女特级片免费视频播放器 | 欧美三级亚洲精品| 午夜福利在线观看吧| 欧美乱色亚洲激情| 国产精品自产拍在线观看55亚洲| 成人一区二区视频在线观看| 国产蜜桃级精品一区二区三区| 国产视频内射| 日韩欧美国产一区二区入口| 欧美人与性动交α欧美精品济南到| 婷婷精品国产亚洲av| 亚洲九九香蕉| 波多野结衣巨乳人妻| 免费在线观看完整版高清| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 黄片大片在线免费观看| 又紧又爽又黄一区二区| 老鸭窝网址在线观看| 国产亚洲精品av在线| 欧美黑人欧美精品刺激| 精品第一国产精品| 国产精品久久久人人做人人爽| 久久久国产成人精品二区| 国产高清视频在线播放一区| tocl精华| 亚洲精品中文字幕在线视频| 国产精华一区二区三区| 黑人巨大精品欧美一区二区mp4| 我要搜黄色片| 欧美日韩国产亚洲二区| 久久中文字幕人妻熟女| 色综合婷婷激情| 日本精品一区二区三区蜜桃| 日本一本二区三区精品| 久久精品国产清高在天天线| 国产高清videossex| 亚洲成人中文字幕在线播放| 久久精品国产99精品国产亚洲性色| 一进一出抽搐gif免费好疼| 国产成人影院久久av| 亚洲av五月六月丁香网| 两个人视频免费观看高清| netflix在线观看网站| 黄色视频不卡| 91老司机精品| 波多野结衣高清作品| 中文字幕最新亚洲高清| 91麻豆精品激情在线观看国产| 久久草成人影院| 婷婷精品国产亚洲av在线| 欧美日韩一级在线毛片| 老司机靠b影院| 日本一二三区视频观看| 成人特级黄色片久久久久久久| 久久久久久久精品吃奶| 曰老女人黄片| 欧美在线一区亚洲| 老司机在亚洲福利影院| 欧美黑人巨大hd| 久久久久国产精品人妻aⅴ院| 久久中文字幕一级| 一本大道久久a久久精品| 精品久久久久久久毛片微露脸| 国产99久久九九免费精品| 国产av一区二区精品久久| 成熟少妇高潮喷水视频| www日本黄色视频网| 国产区一区二久久| 欧美日本亚洲视频在线播放| 欧美丝袜亚洲另类 | or卡值多少钱| 亚洲av中文字字幕乱码综合| 亚洲精品国产精品久久久不卡| xxx96com| 国产精品一区二区三区四区久久| 一级毛片女人18水好多| 床上黄色一级片| 国内精品久久久久精免费| 美女黄网站色视频| 淫妇啪啪啪对白视频| 美女黄网站色视频| 久久人妻av系列| 在线永久观看黄色视频| 国产又黄又爽又无遮挡在线| 午夜福利欧美成人| av福利片在线| 日日夜夜操网爽| 波多野结衣巨乳人妻| 欧美av亚洲av综合av国产av| 亚洲五月婷婷丁香| 午夜日韩欧美国产| 国产精品久久久久久精品电影| 久久久久性生活片| 亚洲精品中文字幕在线视频| 日本免费a在线| 最近最新中文字幕大全电影3| 午夜福利高清视频| 国产精品 国内视频| 一级毛片精品| 久久中文字幕一级| 亚洲第一电影网av| 一进一出抽搐gif免费好疼| 午夜免费激情av| 亚洲va日本ⅴa欧美va伊人久久| 国产精品香港三级国产av潘金莲| 无人区码免费观看不卡| 天天一区二区日本电影三级| 黄片小视频在线播放| 亚洲欧美激情综合另类| 波多野结衣高清无吗| 日本一区二区免费在线视频| 精品人妻1区二区| 午夜老司机福利片| 国产av在哪里看| e午夜精品久久久久久久| 我的老师免费观看完整版| 国产欧美日韩一区二区精品| 伊人久久大香线蕉亚洲五| 无人区码免费观看不卡| 午夜福利18| 国产区一区二久久| 日韩 欧美 亚洲 中文字幕| 欧美极品一区二区三区四区| 亚洲精品久久成人aⅴ小说| 久久伊人香网站| 他把我摸到了高潮在线观看| 国产成人av激情在线播放| 十八禁网站免费在线| 久久久久精品国产欧美久久久| 两个人免费观看高清视频| 日本三级黄在线观看| 久久这里只有精品中国| xxx96com| av有码第一页| 国产伦在线观看视频一区| 性色av乱码一区二区三区2| 超碰成人久久| 国产精品久久视频播放| 欧美日韩黄片免| 麻豆成人午夜福利视频| 久久精品人妻少妇| 精品无人区乱码1区二区| 波多野结衣高清无吗| 国产区一区二久久| av欧美777| 久久久久精品国产欧美久久久| 国产精品免费一区二区三区在线| 精品国产乱子伦一区二区三区| 一二三四社区在线视频社区8| 可以在线观看毛片的网站| 久久久久久久午夜电影| 欧美乱码精品一区二区三区| 亚洲av熟女| 免费无遮挡裸体视频| 欧美黑人巨大hd| 在线观看免费视频日本深夜| 国产av又大| 亚洲 国产 在线| 两人在一起打扑克的视频| 美女免费视频网站| 国产人伦9x9x在线观看| 国产av不卡久久| 国产成人精品久久二区二区免费| 精品免费久久久久久久清纯| 日本成人三级电影网站| 午夜精品久久久久久毛片777| 一进一出抽搐gif免费好疼| x7x7x7水蜜桃| 97人妻精品一区二区三区麻豆| 色综合欧美亚洲国产小说| 两个人免费观看高清视频| xxx96com| 露出奶头的视频| 两性午夜刺激爽爽歪歪视频在线观看 | 精品一区二区三区四区五区乱码| 巨乳人妻的诱惑在线观看| 国产久久久一区二区三区| 国产91精品成人一区二区三区| 黄色a级毛片大全视频| 亚洲中文av在线| 啦啦啦观看免费观看视频高清| 国产精品av视频在线免费观看| 国产aⅴ精品一区二区三区波| 免费高清视频大片| 91九色精品人成在线观看| 亚洲一码二码三码区别大吗| 亚洲国产高清在线一区二区三| 好男人在线观看高清免费视频| 日韩精品青青久久久久久| 国产单亲对白刺激| 亚洲激情在线av| 欧美久久黑人一区二区| 床上黄色一级片| 好看av亚洲va欧美ⅴa在| 69av精品久久久久久| 麻豆成人av在线观看| 最近在线观看免费完整版| 午夜视频精品福利| 国产97色在线日韩免费| 欧美黑人精品巨大| 欧美最黄视频在线播放免费| 欧美日韩亚洲综合一区二区三区_| 手机成人av网站| 一本久久中文字幕| 91字幕亚洲| 国产单亲对白刺激| 国产精品,欧美在线| 国产探花在线观看一区二区| 精品国产乱码久久久久久男人| 亚洲一区中文字幕在线| 制服诱惑二区| 亚洲真实伦在线观看| 亚洲精品国产一区二区精华液| 黄片小视频在线播放| 男男h啪啪无遮挡| 欧美成人性av电影在线观看| 中文亚洲av片在线观看爽| 可以在线观看的亚洲视频| 日韩欧美精品v在线| 人妻久久中文字幕网| 久久天堂一区二区三区四区| 亚洲中文日韩欧美视频| 露出奶头的视频| 精品久久久久久久人妻蜜臀av| 中文字幕熟女人妻在线| 国产精品99久久99久久久不卡| 欧美黑人巨大hd| 老熟妇乱子伦视频在线观看| 久久天躁狠狠躁夜夜2o2o| 免费在线观看影片大全网站| 狂野欧美激情性xxxx| 国产精品,欧美在线| 亚洲中文日韩欧美视频| 母亲3免费完整高清在线观看| 国产99白浆流出| 成年免费大片在线观看| 1024香蕉在线观看| 久久精品国产综合久久久| 一级片免费观看大全| 亚洲精品中文字幕在线视频| 性欧美人与动物交配| 91麻豆精品激情在线观看国产| 成人精品一区二区免费| 人成视频在线观看免费观看| 999精品在线视频| 久久久国产成人精品二区| 宅男免费午夜| 香蕉久久夜色| 午夜福利视频1000在线观看| 村上凉子中文字幕在线| 国产人伦9x9x在线观看| 成在线人永久免费视频| 亚洲人成网站在线播放欧美日韩| 在线看三级毛片| 日韩欧美一区二区三区在线观看| 看黄色毛片网站| 母亲3免费完整高清在线观看| 国产高清视频在线播放一区| 色av中文字幕| 亚洲精品国产精品久久久不卡| 国内精品一区二区在线观看| 好男人在线观看高清免费视频| 成人三级黄色视频| 国产成人av激情在线播放| 免费观看人在逋| 精品国产超薄肉色丝袜足j| 亚洲一区高清亚洲精品| 国产成年人精品一区二区| 特级一级黄色大片| 一本一本综合久久| 在线视频色国产色| 高潮久久久久久久久久久不卡| 少妇裸体淫交视频免费看高清 | 免费搜索国产男女视频| 国产爱豆传媒在线观看 | АⅤ资源中文在线天堂| 久久久久国产精品人妻aⅴ院| 国产欧美日韩一区二区精品| 99re在线观看精品视频| 丁香六月欧美| 国产精品 国内视频| 伊人久久大香线蕉亚洲五| svipshipincom国产片| 露出奶头的视频| 黑人巨大精品欧美一区二区mp4| 国内精品久久久久久久电影| 99久久精品热视频| 亚洲成av人片在线播放无| 99riav亚洲国产免费| 亚洲乱码一区二区免费版| 久久久久久久久免费视频了| 男女之事视频高清在线观看| 亚洲中文字幕一区二区三区有码在线看 | 日韩欧美在线二视频| 女生性感内裤真人,穿戴方法视频| 免费无遮挡裸体视频| 99在线视频只有这里精品首页| 男人舔女人的私密视频| 午夜激情福利司机影院| 2021天堂中文幕一二区在线观| 成年女人毛片免费观看观看9| 搡老岳熟女国产| 亚洲av电影不卡..在线观看| 久久99热这里只有精品18| 一本精品99久久精品77| 免费在线观看完整版高清| 香蕉久久夜色| 欧美日韩乱码在线| 日本黄大片高清| 一卡2卡三卡四卡精品乱码亚洲| 日韩三级视频一区二区三区| 亚洲色图av天堂| 精品无人区乱码1区二区| 美女黄网站色视频| 亚洲最大成人中文| 99久久国产精品久久久| 免费av毛片视频| 人人妻人人看人人澡| 最近最新中文字幕大全免费视频| 欧美精品亚洲一区二区| 香蕉国产在线看| 大型黄色视频在线免费观看| 两人在一起打扑克的视频| 一级a爱片免费观看的视频| svipshipincom国产片| 久久中文看片网| 99精品在免费线老司机午夜| 亚洲无线在线观看| 久久久久国产精品人妻aⅴ院| 黄色片一级片一级黄色片| 国产91精品成人一区二区三区| 一本大道久久a久久精品| 别揉我奶头~嗯~啊~动态视频| 一区福利在线观看| 中文资源天堂在线| 欧美乱妇无乱码| 久久精品夜夜夜夜夜久久蜜豆 | 免费av毛片视频| 欧美 亚洲 国产 日韩一| 国产熟女午夜一区二区三区| 99久久精品热视频| 熟妇人妻久久中文字幕3abv| 欧美日韩福利视频一区二区| 国产免费av片在线观看野外av| 欧美乱色亚洲激情| 亚洲成人久久爱视频| 国产精品免费一区二区三区在线| www日本在线高清视频| 99热6这里只有精品| 欧美黑人欧美精品刺激| 老司机在亚洲福利影院| 国产亚洲欧美98| 国产精品九九99| 小说图片视频综合网站| 99久久综合精品五月天人人| 90打野战视频偷拍视频| 日本黄大片高清| 国产av一区在线观看免费|