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

    激波管內(nèi)黏性流場流動數(shù)值模擬

    2021-07-07 09:41:14汪東旭張尊華梁俊杰衡怡君李格升
    關(guān)鍵詞:邊界層激波管內(nèi)

    汪東旭 張尊華 梁俊杰 衡怡君 李格升

    (武漢理工大學(xué)能源與動力工程學(xué)院 武漢 430063)

    0 引 言

    激波管是一種利用激波實現(xiàn)對工質(zhì)進行加熱獲得目標(biāo)范圍內(nèi)溫度和壓力的實驗裝置[1-2],其包括兩端封閉的圓管和膜片分隔開的驅(qū)動段和被驅(qū)動段.激波管有結(jié)構(gòu)簡單,可在較寬范圍調(diào)節(jié)試驗區(qū)的溫度和壓力的優(yōu)點,廣泛應(yīng)用于研究可壓縮流體現(xiàn)象、高溫氣體、氣相燃燒反應(yīng)和超音速燃燒等領(lǐng)域.

    隨著國內(nèi)外流體數(shù)值模擬軟件的發(fā)展,計算流體力學(xué)(CFD)數(shù)值計算技術(shù)在激波與介質(zhì)相互作用機理研究中占據(jù)較為重要的地位,可應(yīng)用計算流體動力學(xué)流場可視化方式研究和解釋激波管內(nèi)的流場結(jié)構(gòu)[3].激波管在設(shè)計優(yōu)化前,需預(yù)先了解激波管中的流動和激波運行過程,由于數(shù)值模擬可以在較廣的參數(shù)范圍內(nèi)較快給出定量效果,所以可在設(shè)計初期利用數(shù)值模擬來優(yōu)化激波管性能參數(shù),選擇最佳結(jié)構(gòu)和驅(qū)動方案;使用CFD數(shù)值模擬方法可設(shè)置多測點監(jiān)測激波的傳播及參數(shù)變化曲線,能夠直觀預(yù)測激波運行過程中內(nèi)部氣流之間、流體與管壁之間和流動隨時間變化過程中氣體間的相互作用而產(chǎn)生的熱力學(xué)特性,獲得空氣動力參數(shù),如溫度、壓力、密度等參數(shù)分布;利用數(shù)值方法對激波管內(nèi)流場流動進行計算對實際設(shè)備的傳感器測量過程也有參考意義.

    近年來,國內(nèi)外學(xué)者針對激波管內(nèi)激波和流體運動特性進行了一系列試驗與模擬研究,無黏模型和層流黏性模型在激波管的流動數(shù)值計算中得到廣泛應(yīng)用.Lamnaouer等[4]針對高壓激波管,模擬了激波在非反應(yīng)流和反應(yīng)流中的傳播和反射,對激波管側(cè)壁傳熱和反射激波/邊界層相互作用等非理想狀態(tài)進行了量化研究,模型能夠準(zhǔn)確地模擬激波和膨脹波的傳播和反射以及激波反射后的流動不均勻性;Amrollah等[5]在兩種不同驅(qū)動氣體壓力和膜片厚度的情況下,基于激波管二維軸對稱模型研究了正激波在層流粘性流體中的運行過程,仿真激波通過的時間間隔數(shù)據(jù)與試驗數(shù)據(jù)一致性良好;Davidson等[6]利用激光診斷技術(shù)對激波管邊界層效應(yīng)進行了試驗研究,結(jié)果表明:邊界層與激波分岔發(fā)生在湍流區(qū)域,它們受激波管結(jié)構(gòu)和激波前后條件影響;何中偉等[7]通過試驗研究發(fā)現(xiàn)實際激波管中管壁壁面邊界層基本上是湍流邊界層占主導(dǎo),激波與邊界層干擾中,大多是激波與湍流邊界層的相互干擾.為更準(zhǔn)確反映激波管中的物理現(xiàn)象,有必要對激波管進行湍流流動模擬研究并驗證其結(jié)果的可信度,其結(jié)果可為下一步開展激波管流體力學(xué)與化學(xué)動力學(xué)耦合的參數(shù)化研究提供重要理論依據(jù).

    利用數(shù)值模擬的方法,模擬理想氣體空氣在總長為2.08 m、直徑為67.00 mm的激波管內(nèi)湍流粘性流場的流動過程,并與相關(guān)試驗數(shù)據(jù)進行對比,驗證模擬方法的可靠性.依據(jù)不同時刻下激波傳播位移和入射激波后溫度、壓力等參數(shù)隨時間的變化,總結(jié)激波運行規(guī)律,獲得激波管最大試驗時間限度.

    1 數(shù)值模擬方法和邊界條件

    1.1 幾何模型描述

    研究對象為軸對稱結(jié)構(gòu)激波管,采用二維模型對計算區(qū)域進行建模,見圖1.該模型參考了阿米爾卡比爾理工大學(xué)的小尺寸激波管全幾何形狀,激波管長度為2.08 m,直徑為67.00 mm,主要由高壓段、中間膜片、低壓段三個部分組成,破膜方式為壓差破膜,當(dāng)高壓段和低壓段達到一定壓差時膜片破裂產(chǎn)生激波.

    圖1 激波管計算區(qū)域

    X=0處定義為膜片位置,高壓驅(qū)動段長度為0.64 m,低壓被驅(qū)動段長度為1.44 m,為記錄不同位置的流體壓力隨時間的變化,在高壓段靠近膜片處至低壓段尾端依次分別設(shè)置了9個壓力傳感器監(jiān)測點Sn,監(jiān)測點位置見表1.

    表1 監(jiān)測點位置

    1.2 控制方程

    為描述激波管內(nèi)的流體流動,采用非定常雷諾平均納維-斯托克斯方程進行數(shù)學(xué)建模.其流體力學(xué)控制方程包括連續(xù)性方程、動量方程和能量方程,為

    連續(xù)方程:

    (1)

    動量方程:

    (2)

    能量方程:

    (3)

    式中:keff為有效熱傳導(dǎo)系數(shù);Jj為組分j的擴散流量;Sh為其他體積熱源項,不考慮化學(xué)反應(yīng)的情況下,Sh=0,以上方程右邊的前三項分別表示熱傳導(dǎo)、組分?jǐn)U散和黏性耗散引起的能量運輸.

    控制方程通過有限體積法轉(zhuǎn)換為可以用數(shù)值方法求解的代數(shù)方程,在每個控制體內(nèi)積分控制方程,進而產(chǎn)生基于控制體的每一個變量都守恒的離散方程.對與任意控制體積V的積分形式的方程有:

    (4)

    式中:ρ為密度;v為速度矢量;A為曲面面積矢量;Γφφ為φ的擴散系數(shù);φ為φ的梯度,對二維問題即為(?φ/?x)i+(?φ/?y)j;Sφ為單位體積φ的源項.對于每個網(wǎng)格區(qū)域:

    (5)

    式中:Nfaces為封閉單元面的數(shù)量;φf為通過表面f的對流量;ρfvf·Af為通過表面的質(zhì)量流量;Af為表面f的面積;φf為表面f上φ的梯度;V為單元體積.

    1.3 模型模擬方法的選擇

    湍流是一種高度復(fù)雜的三維非穩(wěn)態(tài)、帶旋轉(zhuǎn)的不規(guī)則流體運動[8].超聲速流動模擬中常用的有Baldwin-Lomax代數(shù)渦黏性模型,Spalart-Allmaras湍流模型,k-ε雙方程模型,為仿真激波管內(nèi)多尺度復(fù)雜流動的激波及非定常流場流動,選用模擬復(fù)雜流動和高壓力梯度流動均有較好性能的雙方程湍流模型,該模型用于計算有旋的均勻剪切流,平面混合流,平面射流,圓形射流,管道內(nèi)充分發(fā)展流動,都取得了與實驗數(shù)據(jù)比較一致的結(jié)果[9].壓力基求解器主要適用于低速不可壓縮流體,故采用適用于高速可壓縮流體的密度基隱式求解器.求解器的時間類型設(shè)定為瞬態(tài),F(xiàn)lux-Type選擇Roe-FDS差分格式,Roe格式可在流動不連續(xù)性附近產(chǎn)生振蕩解,標(biāo)量梯度采用基于格林高斯網(wǎng)格節(jié)點方法,空間離散化采用二階迎風(fēng)格式,為提高計算精度,選擇雙精度并行模式進行數(shù)值模擬.在求解初始化過程中,利用patch命令對高低壓區(qū)分別設(shè)置相應(yīng)的初始壓力和溫度,仿真迭代時間步長為1×10-7s,最大迭代次數(shù)為20.

    1.4 網(wǎng)格劃分和網(wǎng)格敏感性驗證

    運用前處理軟件ICEM劃分網(wǎng)格,為了加速計算收斂,采用良好的結(jié)構(gòu)化四邊形網(wǎng)格對計算區(qū)域進行離散化,在壁面附近創(chuàng)建精細(xì)的網(wǎng)格層以捕捉邊界層效應(yīng),局部計算網(wǎng)格見圖2.

    圖2 局部計算網(wǎng)格

    FLUENT軟件數(shù)值模擬結(jié)果與網(wǎng)格劃分質(zhì)量有很強的關(guān)聯(lián)性,為減少數(shù)值誤差對模擬結(jié)果的影響,進行了網(wǎng)格敏感性研究.在長為2.08 m、直徑為67 mm的模型中,繪制網(wǎng)格數(shù)目分別為700×15,1 000×20,2 000×25,2 500×40和2 999×80的五種計算網(wǎng)格,X軸方向基礎(chǔ)網(wǎng)格尺寸分別約為2.971,2.080,1.040,0.832和0.694 mm,Y方向全局加密比例因子設(shè)置為1.111,以上網(wǎng)格正交質(zhì)量接近1,網(wǎng)格最大長寬比均低于邊界值80,處于理想范圍內(nèi)[10].選用相同的邊界條件(air-air,p4=500 kPa,p1=50 kPa,T4,1=300 K)和數(shù)值方法進行數(shù)值模擬,記錄t=1.0 ms時刻激波管徑向Y=0.03 m位置所對應(yīng)的3區(qū)和2區(qū)接觸面處溫度邊界層內(nèi)等值線軸向坐標(biāo)值,見圖3.由圖3可知:發(fā)現(xiàn)隨著網(wǎng)格數(shù)增大,徑向Y=0.03 m位置所對應(yīng)的等溫線軸向坐標(biāo)值差異越小,當(dāng)網(wǎng)格數(shù)增大到2 000×25時,激波管溫度邊界層內(nèi)溫度分布幾乎不受網(wǎng)格變化影響,綜合考慮計算成本和求解精確性等方面,選擇網(wǎng)格數(shù)為2 500×40(X軸方向基礎(chǔ)網(wǎng)格尺寸為0.832 mm)的網(wǎng)格劃分方案開展二維激波管內(nèi)流體流動數(shù)值模擬研究.

    圖3 徑向Y=0.03 m時,激波管內(nèi)溫度等值線軸向坐標(biāo)分布

    實際激波管中可觀測到湍流和層流兩種邊界層,層流邊界層向湍流邊界層轉(zhuǎn)變是通過臨界雷諾數(shù)定義的,層流邊界層只有在被驅(qū)動段初始壓力很低、激波較弱等條件下可能存在,邊界層大部分區(qū)域是湍流的.y+值為管壁有界流動的無量綱壁面距離,是評判各種模型壁面邊界層網(wǎng)格分辨率的一個重要參數(shù)[11].采用標(biāo)準(zhǔn)壁面函數(shù)劃分湍流邊界層網(wǎng)格,合理的y+值應(yīng)在(30,300)范圍內(nèi),以保證流場結(jié)構(gòu)解析到粘性子層,壁面邊界層網(wǎng)格y+值計算公式為

    (6)

    式中:Δy為第一層網(wǎng)格與壁面之間的距離;ρ為流體密度;UT為流體速度;μ為黏性系數(shù),隨著網(wǎng)格分辨率的增加,y+值會減小,利用這種規(guī)律使網(wǎng)格分辨率達到理想值.通過調(diào)整管壁第一層網(wǎng)格高度后經(jīng)計算獲得.由于本文流場是瞬態(tài)流動,需利用不同時刻壁面邊界層y+值判斷邊界層網(wǎng)格設(shè)置質(zhì)量,值設(shè)定為0.05 mm時,記錄Case1激波運行全過程壁面邊界層y+值,不同時刻壁面邊界層y+值均在(30,36)范圍內(nèi),故認(rèn)為本文湍流邊界層網(wǎng)格設(shè)置是合理的,第一層網(wǎng)格與壁面之間的距離值設(shè)置為5×10-5.

    1.5 初始邊界條件

    由圖2可知,模型共有壁面(wall)、軸對稱(axis-symmetry)、內(nèi)部面(interior)三種邊界類型,忽略浮力和重力影響:①激波管左右兩端端壁和兩側(cè)管壁均設(shè)定為絕熱、無滑移wall邊界;②管體設(shè)定為axis-symmetry邊界;③膜片部分設(shè)定為interior邊界.

    數(shù)值模擬初始條件見表2.本文研究背景為無反應(yīng)流場流動,為避免非理想介質(zhì)對模擬結(jié)果精確性的影響,高低壓段氣體工質(zhì)均設(shè)定為理想空氣,黏性常數(shù)為1.789 4×10-5kg/(m·s).基于以下邊界條件,研究激波的傳播規(guī)律,分析激波運行和激波后溫度邊界層的分布規(guī)律.

    表2 數(shù)值模擬的邊界條件

    2 數(shù)值模擬方法的驗證

    為驗證湍流粘性流體數(shù)值模擬方法的適用性,需將數(shù)值計算得到的結(jié)果與實驗數(shù)據(jù)進行對照分析,文獻[5]中激波管幾何模型與本文相同,實驗設(shè)定的邊界條件見表3.本文激波管湍流模型采用相同的實驗條件進行模擬計算,湍流模擬算例記為Case2和Case3,激波位移計算結(jié)果情況見圖4.

    表3 文獻實驗條件

    圖4 不同工況下入射激波所抵達的位置

    由圖4可知,Case2條件下,t=0.002 2 s時,正激波抵達位置的橫軸坐標(biāo)為X1=0.966 m,入射激波速度為vs=452.73 m/s;Case3條件下,t=0.002 0 s時,正激波抵達位置的橫軸坐標(biāo)為X2=0.991 m,入射激波速度為vs=495.5 m/s.

    將湍流模擬不同時刻入射激波速度與文獻數(shù)據(jù)進行對比,見表4.由表4可知,本文湍流數(shù)值模擬算例2在0.002 2 s時刻得到的激波速度相對文獻試驗數(shù)據(jù)的誤差為0.291%;算例3在0.002 0 s時刻激波速度相對試驗數(shù)據(jù)的誤差為2.641%.數(shù)值模擬結(jié)果與試驗結(jié)果存在一定的誤差,考慮到試驗過程中激波管管壁為非絕熱壁面,同時實際流體的黏性效應(yīng)和激波管壁面的摩擦效應(yīng)會使激波強度有一定的損失,而模擬時假定激波管管壁為絕熱恒溫?zé)o滑移壁面,因此模擬結(jié)果與試驗結(jié)果必然會產(chǎn)生一定的誤差.湍流模擬結(jié)果相對于層流模擬更接近試驗數(shù)據(jù),且兩個算例模擬結(jié)果產(chǎn)生的誤差數(shù)值小于5%,均在數(shù)值模擬誤差可接受范圍內(nèi),綜上所述,可認(rèn)為該湍流粘性流體數(shù)值模擬方法適合用于預(yù)測激波管內(nèi)激波的運動特性.

    表4 湍流模擬結(jié)果與文獻試驗、模擬值對比

    3 模擬結(jié)果與討論

    3.1 激波的產(chǎn)生與傳播

    在仿真初始狀態(tài),高低壓段充入不同壓力的空氣,此時管內(nèi)包含兩個區(qū)域:初始高壓氣體區(qū)域(4區(qū))、初始低壓氣體區(qū)域(1區(qū)).模擬過程從膜片破裂時刻(t=0 s)開始,假設(shè)X=0 m處的膜片瞬時破裂消失,見圖5.由圖5可知,初始時刻,4區(qū)與1區(qū)之間的壓力間斷面稱為接觸面.

    圖5 t=0 s時刻,膜片破裂時激波管內(nèi)瞬時壓力(Pa)云圖(Case4:p4=500 kPa,p1=20 kPa,T4,1=300 K)

    高低壓氣體之間的壓差會使氣體進行非定常運動,同時形成正激波向低壓段傳播,入射激波和反射激波對工質(zhì)進行壓縮加熱加壓,由于激波對介質(zhì)的加熱時間極短,可認(rèn)為整個壓縮加熱過程是絕熱、非等熵、均勻加熱的.圖6為t=1.50 ms時刻沿軸向位置管內(nèi)入射激波后管體溫度局部分布云圖.由圖6可知:入射激波壓縮區(qū)域溫度高于初始溫度,膨脹波壓縮區(qū)域溫度低于初始溫度,其中具有溫度梯度的薄層稱為溫度邊界層[12].壁面邊界層溫度高于中心主流區(qū)溫度,這是由于粘性影響,氣體在邊界層內(nèi)減速,動能轉(zhuǎn)化成熱能,邊界層內(nèi)的流動為非等熵過程,同時存在動能交換和熱交換過程,因此溫度邊界層處的溫度要高于中心主氣流.

    圖6 t=1.50 ms時,沿軸向位置管內(nèi)溫度(K)局部分布(Case4:p4=500 kPa,p1=20 kPa,T4,1=300 K)

    此時,激波管內(nèi)形成四個區(qū)域,沿軸向位置溫度分布見圖7,膜片破裂后,產(chǎn)生的正激波由膜片位置向低壓段傳播,正激波與高低壓段氣體接觸面之間的區(qū)域稱為2區(qū),產(chǎn)生的膨脹波向高壓段傳播,膨脹波與接觸面之間的區(qū)域稱為3區(qū);正激波到達右側(cè)尾端端壁后會發(fā)生反射產(chǎn)生反射激波,反射激波經(jīng)過的區(qū)域稱為5區(qū),即試驗區(qū).粘性流體流經(jīng)管體壁面時,入射激波后壁面附近的流體區(qū)域會形成一定厚度的邊界層,邊界層內(nèi)激波速度會迅速降低,在管壁處的速度為零,激波管內(nèi)2區(qū)、3區(qū)及5區(qū)壁面因流體粘性效應(yīng)附近均會出現(xiàn)邊界層,1區(qū)和4區(qū)的氣體處于靜止?fàn)顟B(tài)不存在邊界層.

    圖7 t=1.50 ms時,沿軸向位置管內(nèi)溫度(K)分布范圍(Case4:p4=500 kPa,p1=20 kPa,T4,1=300 K)

    采用不同時刻的壓力和溫度云圖來反映激波移動的過程,見圖8~9,給出了激波產(chǎn)生、傳播及反射過程中流體的壓力和溫度變化云圖.t=0.1 ms時,膜片消失后,入射激波向低壓段傳播,經(jīng)過的2區(qū)壓力和溫度迅速升高,膨脹波向高壓段傳播,膨脹波后的區(qū)域壓力和溫度會降低;t=1.0 ms時,由壓力云圖可以清晰觀察到膨脹波和入射激波面的運行位置,根據(jù)溫度云圖可觀察到入射激波與接觸面之間的2區(qū),隨著時間推移,2區(qū)范圍會逐漸增大,這是因為激波運行速度大于接觸面前進速度;當(dāng)t=2.5 ms時,入射激波抵達右尾端端壁并產(chǎn)生反射激波向左傳播,經(jīng)過的5區(qū)壓力和溫度再一次升高達到穩(wěn)定狀態(tài)并持續(xù)一段時間,膨脹波到達左端端壁后發(fā)生反射產(chǎn)生反射膨脹波向右傳播;t=3.0 ms時,反射激波向左傳播與接觸面相遇并產(chǎn)生二次反射激波,二次反射激波向右傳播;t=3.2 ms時刻,二次反射激波會對5區(qū)進行壓縮作用,使5區(qū)壓力和溫度不均衡上升,試驗區(qū)的氣體狀態(tài)發(fā)生改變,試驗區(qū)環(huán)境被破壞,距尾端最遠的測點壓力穩(wěn)定狀態(tài)最先被破壞,當(dāng)二次反射激波行至尾端端壁,5區(qū)有效試驗時間結(jié)束;t=3.5 ms時,二次反射激波已到達尾端端壁,再次反射產(chǎn)生三次反射激波,這個過程會持續(xù)多次,直至激波管內(nèi)整體壓力處于平衡狀態(tài).

    圖8 Case4激波管內(nèi)流體瞬時壓力(Pa)云圖

    圖9 Case4激波管內(nèi)流體溫度(K)云圖

    3.2 激波管5區(qū)試驗時間限度的計算

    以空氣(理想狀態(tài))為例,計算移除高壓段與低壓段之間膜片后可以獲得的激波管5區(qū)有效試驗時間,將反射激波后的流動特性穩(wěn)定的時間定義為試驗時間,試驗時間是評估激波管設(shè)備性能的重要參數(shù),可為進行著火延時測量提供重要理論參考依據(jù).

    激波管中流體運動位置-時間平面圖是一種重要的工具,它可以很好地估計在移除膜片后可獲得的有效試驗測試時間.圖10~11為Case4(air-air,p4=500 kPa,p1=20 kPa,T4,1=300 K)運行時壓力和溫度分布的預(yù)測x-t圖.圖10壓力分布x-t圖中接觸面沒有出現(xiàn),因為它連續(xù)的跟隨入射激波后,圖11為整個激波管設(shè)備的溫度變化,激波、膨脹波和接觸面都呈現(xiàn)在圖中,激波之后是接觸面,直到反射激波與接觸面相互作用,然后被進一步反射.膨脹波到達左端壁面會產(chǎn)生反射膨脹波,由圖11可知,膨脹波1曲線斜率小于和接觸面的曲線斜率,即接觸面?zhèn)鞑ニ俣却笥诜瓷渑蛎洸ǎ虼?,高壓段長度足夠條件下,反射膨脹波到達低壓段末端之前,激波就已與接觸面相交.二次反射激波到達右端壁面時,標(biāo)志激波管測試時間結(jié)束.

    圖10 壓力云圖中,激波的位置隨時間變化曲線

    圖11 溫度云圖中,激波和高低壓段間的接觸面的位置隨時間變化曲線

    在激波管壁面設(shè)置了8個監(jiān)測點,另外在低壓段端面中心設(shè)置了1個監(jiān)測點,圖12為Case4工況下激波管壁面監(jiān)測點壓力隨時間變化情況.開始時,高壓段測點S1保持初始壓力不變,膨脹波經(jīng)過后壓力持續(xù)降低,待反射膨脹波達到時再次降低;低壓段壁面測點S2~S8入射激波到達時壓力均有階躍升高,S6~S8處反射激波經(jīng)過,S6處距尾端較遠,反射激波已遇到接觸面產(chǎn)生二次反射激波,反射強度下降,因此壓力二次升高值低于測點S7、S8.測點S9不同于側(cè)壁測點S8,只有一次壓力階躍過程,入射激波與尾端壁面相遇時刻就是5區(qū)的開始.

    圖12 激波管壁面監(jiān)測點壓力變化

    實際試驗過程中,通常有兩種方式表示測量時間,分別為側(cè)壁傳感器壓力穩(wěn)定時間t側(cè)和端壁傳感器壓力穩(wěn)定時間t端,壓力測點相當(dāng)于實際設(shè)備上的壓力傳感器,以上方式也是5區(qū)試驗時間限度的重要參考參數(shù).表5為激波管試驗時間限度計算結(jié)果比較.

    表5 激波管試驗時間限度計算結(jié)果比較

    由表5可知,該工況下幾種試驗時間限度估計方式得出的結(jié)果差異較小,側(cè)壁測點S8得到的壓力穩(wěn)定時間為1.039 ms,比端壁測點時間S9壓力穩(wěn)定時間1.073 ms要短,這是因為邊界層效應(yīng)的影響,反射激波會發(fā)生激波分岔和激波衰減,離反射端面處的距離越近,受邊界層影響越小,側(cè)壁測點壓力開始上升至穩(wěn)定的時間稍長于尾端端壁測點.然而,激波管在該算例下x-t圖得到的測試時間為1.169 ms,大于端壁測點的壓力穩(wěn)定時間,綜上所述,該算例理論有效試驗時間為1.169 ms,實際測量有效試驗時間最大為1.073 ms.

    4 結(jié) 論

    1) 采用FLUENT中的湍流模型進行激波管流體流動數(shù)值模擬,計算結(jié)果與層流模型模擬結(jié)果以及實驗結(jié)果進行比較,與實驗結(jié)果吻合較好,表明湍流模型對激波管流動具有較好的模擬效果,本文CFD數(shù)值方法是高效可靠的.

    2) 詳細(xì)模擬了激波的產(chǎn)生和傳播過程,給出了激波形成、傳播、反射和激波與接觸面相互作用的整個過程,接觸面與反射激波相遇時刻和二次反射激波行至端壁的時間是有效試驗時間限度的重要影響因素.

    3) 本數(shù)值方法能夠提供x-t關(guān)系圖,從該圖中可以確定激波管實驗測量有效持續(xù)時間,激波管5區(qū)內(nèi)側(cè)壁、端壁壓力傳感器測得的試驗時間長度均小于x-t圖獲得的最大的測試時間,且差值較小,符合預(yù)期結(jié)果.

    猜你喜歡
    邊界層激波管內(nèi)
    一種基于聚類分析的二維激波模式識別算法
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    上傾管內(nèi)油水兩相流流型實驗研究
    斜激波入射V形鈍前緣溢流口激波干擾研究
    適于可壓縮多尺度流動的緊致型激波捕捉格式
    揚礦管內(nèi)高速螺旋流的數(shù)值模擬與仿真
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    腕管內(nèi)原發(fā)性腱鞘結(jié)核誤診1例
    非特征邊界的MHD方程的邊界層
    微細(xì)管內(nèi)CO2流動沸騰換熱特性研究
    99热6这里只有精品| 女同久久另类99精品国产91| 国产在视频线在精品| 天堂影院成人在线观看| 日本三级黄在线观看| 桃红色精品国产亚洲av| 精品欧美国产一区二区三| 亚洲精品久久国产高清桃花| 噜噜噜噜噜久久久久久91| 老汉色av国产亚洲站长工具| av国产免费在线观看| 精品午夜福利视频在线观看一区| 在线观看免费午夜福利视频| 久久久久九九精品影院| 丝袜美腿在线中文| 一卡2卡三卡四卡精品乱码亚洲| 51国产日韩欧美| 午夜福利在线观看吧| 色综合亚洲欧美另类图片| 国产真实伦视频高清在线观看 | 亚洲精品美女久久久久99蜜臀| 久久久久久国产a免费观看| 免费高清视频大片| 国产黄a三级三级三级人| 热99在线观看视频| 久久久成人免费电影| 国产蜜桃级精品一区二区三区| 色噜噜av男人的天堂激情| 日本精品一区二区三区蜜桃| 淫秽高清视频在线观看| 性色av乱码一区二区三区2| 久久婷婷人人爽人人干人人爱| x7x7x7水蜜桃| 免费观看的影片在线观看| 欧美成狂野欧美在线观看| 18禁裸乳无遮挡免费网站照片| 国产成人影院久久av| 久久久久九九精品影院| 久久久久国内视频| 亚洲欧美一区二区三区黑人| 亚洲激情在线av| 人妻夜夜爽99麻豆av| 成人国产一区最新在线观看| 欧美+日韩+精品| 国产精品98久久久久久宅男小说| 在线天堂最新版资源| 1000部很黄的大片| 久久精品夜夜夜夜夜久久蜜豆| 国产精品爽爽va在线观看网站| a在线观看视频网站| 1024手机看黄色片| x7x7x7水蜜桃| 欧美一级a爱片免费观看看| 99久久99久久久精品蜜桃| 九九在线视频观看精品| 一区二区三区高清视频在线| 亚洲av成人不卡在线观看播放网| 亚洲真实伦在线观看| 91av网一区二区| 免费电影在线观看免费观看| 天天添夜夜摸| 亚洲aⅴ乱码一区二区在线播放| 免费观看的影片在线观看| 很黄的视频免费| tocl精华| 精品熟女少妇八av免费久了| 国产精品日韩av在线免费观看| 又黄又粗又硬又大视频| 一个人看的www免费观看视频| 日韩免费av在线播放| 此物有八面人人有两片| 黄色女人牲交| 嫩草影院精品99| 18美女黄网站色大片免费观看| tocl精华| 国产成人av教育| 欧美一级毛片孕妇| 国产av一区在线观看免费| 窝窝影院91人妻| 波野结衣二区三区在线 | 亚洲精品一卡2卡三卡4卡5卡| 欧美黄色淫秽网站| 国产高潮美女av| 免费在线观看亚洲国产| 两个人视频免费观看高清| 18禁美女被吸乳视频| 欧美成人一区二区免费高清观看| 美女 人体艺术 gogo| 国产久久久一区二区三区| 亚洲天堂国产精品一区在线| 97超视频在线观看视频| 国产爱豆传媒在线观看| 国产成人系列免费观看| 国产成年人精品一区二区| 搡老妇女老女人老熟妇| 男人的好看免费观看在线视频| 国产美女午夜福利| 免费看美女性在线毛片视频| 亚洲精品成人久久久久久| 他把我摸到了高潮在线观看| 99在线人妻在线中文字幕| 国内揄拍国产精品人妻在线| 日本一本二区三区精品| 国产激情欧美一区二区| 国产免费男女视频| 国产毛片a区久久久久| 国产毛片a区久久久久| 国产极品精品免费视频能看的| 久久久久亚洲av毛片大全| 久久这里只有精品中国| 亚洲av成人av| 日本a在线网址| 欧美激情久久久久久爽电影| 99久久精品一区二区三区| 成人永久免费在线观看视频| 内射极品少妇av片p| 国产精品久久久久久人妻精品电影| 91久久精品国产一区二区成人 | 麻豆国产av国片精品| 三级毛片av免费| 婷婷丁香在线五月| x7x7x7水蜜桃| 国产精品一及| 嫩草影院精品99| x7x7x7水蜜桃| 国产激情欧美一区二区| 亚洲av五月六月丁香网| av国产免费在线观看| 少妇裸体淫交视频免费看高清| 亚洲五月婷婷丁香| 成人亚洲精品av一区二区| 亚洲中文日韩欧美视频| 波野结衣二区三区在线 | 老熟妇仑乱视频hdxx| 亚洲乱码一区二区免费版| 精品乱码久久久久久99久播| 免费大片18禁| 一二三四社区在线视频社区8| 中文亚洲av片在线观看爽| 午夜福利在线在线| 他把我摸到了高潮在线观看| 噜噜噜噜噜久久久久久91| 久久6这里有精品| 亚洲欧美日韩高清专用| 精品久久久久久成人av| 小蜜桃在线观看免费完整版高清| 亚洲一区二区三区色噜噜| 草草在线视频免费看| 一夜夜www| 18禁黄网站禁片午夜丰满| 人妻久久中文字幕网| 日本三级黄在线观看| 色播亚洲综合网| 99在线视频只有这里精品首页| 午夜激情欧美在线| 国产精品 欧美亚洲| 国产精品99久久久久久久久| 国内精品美女久久久久久| 美女大奶头视频| 免费高清视频大片| 国产欧美日韩一区二区精品| 国产精品一区二区免费欧美| 久久久久国产精品人妻aⅴ院| 99国产极品粉嫩在线观看| 黄色视频,在线免费观看| 黄色片一级片一级黄色片| 国产色婷婷99| 亚洲色图av天堂| 在线观看一区二区三区| 最近视频中文字幕2019在线8| 欧美日韩福利视频一区二区| 日日夜夜操网爽| 在线观看66精品国产| 在线天堂最新版资源| 亚洲精品日韩av片在线观看 | 国内精品久久久久久久电影| www国产在线视频色| 日本一本二区三区精品| 搞女人的毛片| 人妻夜夜爽99麻豆av| 久久99热这里只有精品18| 偷拍熟女少妇极品色| 国产av在哪里看| 久久香蕉精品热| 91av网一区二区| 制服丝袜大香蕉在线| 99精品久久久久人妻精品| 日韩欧美在线二视频| 91九色精品人成在线观看| 一区二区三区高清视频在线| 91麻豆精品激情在线观看国产| 母亲3免费完整高清在线观看| 国产成人欧美在线观看| 久久久久久久亚洲中文字幕 | 啦啦啦韩国在线观看视频| 国产爱豆传媒在线观看| 国产精品精品国产色婷婷| 亚洲欧美激情综合另类| 最近最新中文字幕大全电影3| 天天一区二区日本电影三级| 久久人人精品亚洲av| 舔av片在线| 成年女人毛片免费观看观看9| 一个人看视频在线观看www免费 | 99热6这里只有精品| 日日夜夜操网爽| 国产高潮美女av| 99视频精品全部免费 在线| 免费在线观看日本一区| 亚洲色图av天堂| 久久这里只有精品中国| 亚洲欧美日韩卡通动漫| 夜夜夜夜夜久久久久| 亚洲成av人片在线播放无| 两个人的视频大全免费| 国产欧美日韩一区二区三| 90打野战视频偷拍视频| av天堂在线播放| 91麻豆精品激情在线观看国产| 老鸭窝网址在线观看| 蜜桃亚洲精品一区二区三区| 宅男免费午夜| 日本免费a在线| 美女cb高潮喷水在线观看| 色av中文字幕| 国产精品99久久99久久久不卡| 99久久精品热视频| 俺也久久电影网| 老司机福利观看| 老司机在亚洲福利影院| 男人舔女人下体高潮全视频| 极品教师在线免费播放| 最近在线观看免费完整版| 国产成人av激情在线播放| av在线蜜桃| 久久久久久大精品| 无人区码免费观看不卡| 9191精品国产免费久久| 成人三级黄色视频| 欧美成人性av电影在线观看| 可以在线观看的亚洲视频| 婷婷六月久久综合丁香| 亚洲国产欧美人成| 好男人在线观看高清免费视频| 国产aⅴ精品一区二区三区波| 99久久九九国产精品国产免费| 欧美在线黄色| 久久久久久久精品吃奶| 欧美中文综合在线视频| 特级一级黄色大片| 亚洲成人精品中文字幕电影| 黄色视频,在线免费观看| 一级黄片播放器| 日日夜夜操网爽| 12—13女人毛片做爰片一| 国产免费av片在线观看野外av| 99在线人妻在线中文字幕| 在线a可以看的网站| 国内精品美女久久久久久| 免费电影在线观看免费观看| 国产亚洲av嫩草精品影院| 欧美3d第一页| 天天一区二区日本电影三级| 久久精品夜夜夜夜夜久久蜜豆| 欧美成人一区二区免费高清观看| 老司机福利观看| 亚洲欧美一区二区三区黑人| 狠狠狠狠99中文字幕| 色综合欧美亚洲国产小说| av福利片在线观看| 国产精品一及| 黄片大片在线免费观看| 色尼玛亚洲综合影院| 午夜福利在线观看免费完整高清在 | 国产精品永久免费网站| 精品一区二区三区av网在线观看| 亚洲最大成人中文| 亚洲乱码一区二区免费版| www.色视频.com| 手机成人av网站| 欧美又色又爽又黄视频| 少妇的逼好多水| 国产成人aa在线观看| 成人欧美大片| 久久久久亚洲av毛片大全| 亚洲av第一区精品v没综合| www国产在线视频色| 免费av不卡在线播放| 欧美午夜高清在线| 国产精品 国内视频| 狂野欧美激情性xxxx| 欧美日本视频| 亚洲乱码一区二区免费版| 高清日韩中文字幕在线| 国产极品精品免费视频能看的| 欧美黑人欧美精品刺激| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| 久久欧美精品欧美久久欧美| 51午夜福利影视在线观看| 精品久久久久久成人av| 激情在线观看视频在线高清| 国产精品一及| 欧美一区二区亚洲| 少妇的逼好多水| 欧美在线一区亚洲| 亚洲精品粉嫩美女一区| 在线观看日韩欧美| 欧美日本视频| 国产91精品成人一区二区三区| 一级作爱视频免费观看| 欧美乱码精品一区二区三区| 男女下面进入的视频免费午夜| 欧美日韩一级在线毛片| 成年版毛片免费区| 一二三四社区在线视频社区8| 又爽又黄无遮挡网站| 久久国产精品影院| 十八禁网站免费在线| 精品日产1卡2卡| 国产极品精品免费视频能看的| 香蕉丝袜av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 九九久久精品国产亚洲av麻豆| 亚洲av五月六月丁香网| 最新中文字幕久久久久| 久久九九热精品免费| 中文字幕av在线有码专区| 2021天堂中文幕一二区在线观| 日韩欧美精品免费久久 | 成人鲁丝片一二三区免费| 国产探花极品一区二区| 亚洲av美国av| 日本一二三区视频观看| 久久精品夜夜夜夜夜久久蜜豆| 一级作爱视频免费观看| 亚洲国产精品合色在线| 中文字幕人妻丝袜一区二区| 99久久综合精品五月天人人| 又黄又粗又硬又大视频| 男女下面进入的视频免费午夜| 亚洲成人久久爱视频| 午夜两性在线视频| 乱人视频在线观看| 99热这里只有是精品50| 精品国产三级普通话版| 麻豆成人午夜福利视频| 一个人看的www免费观看视频| 国产一区二区三区在线臀色熟女| 精品日产1卡2卡| 国产三级中文精品| 狠狠狠狠99中文字幕| 午夜精品一区二区三区免费看| 亚洲美女视频黄频| 成年免费大片在线观看| 国产一区二区三区视频了| 国产精品 国内视频| 好看av亚洲va欧美ⅴa在| 男女下面进入的视频免费午夜| 禁无遮挡网站| 国产在线精品亚洲第一网站| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美午夜高清在线| 中文字幕高清在线视频| 亚洲人成网站在线播放欧美日韩| 国产精品爽爽va在线观看网站| 亚洲男人的天堂狠狠| 麻豆成人av在线观看| 老汉色av国产亚洲站长工具| 特大巨黑吊av在线直播| 熟女电影av网| 乱人视频在线观看| 亚洲avbb在线观看| 日韩中文字幕欧美一区二区| 久久久久久大精品| 午夜久久久久精精品| 国产av在哪里看| 免费搜索国产男女视频| 成人欧美大片| 国产一区二区在线观看日韩 | 亚洲人成电影免费在线| 国产精品99久久久久久久久| 国产三级在线视频| 欧美日韩瑟瑟在线播放| 国产主播在线观看一区二区| 日韩亚洲欧美综合| 亚洲无线观看免费| 国产私拍福利视频在线观看| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美 | 99热只有精品国产| 国产一级毛片七仙女欲春2| 国产91精品成人一区二区三区| av女优亚洲男人天堂| 亚洲最大成人中文| 国产蜜桃级精品一区二区三区| 国产麻豆成人av免费视频| 日本成人三级电影网站| 国产三级黄色录像| 老司机午夜福利在线观看视频| 身体一侧抽搐| av在线蜜桃| 国产极品精品免费视频能看的| 丰满乱子伦码专区| 中出人妻视频一区二区| 不卡一级毛片| 18禁美女被吸乳视频| 午夜免费男女啪啪视频观看 | 国产精品亚洲一级av第二区| 国产精品亚洲一级av第二区| 九九热线精品视视频播放| 熟女少妇亚洲综合色aaa.| 一进一出抽搐gif免费好疼| 成人无遮挡网站| 无限看片的www在线观看| 国产成年人精品一区二区| 亚洲一区高清亚洲精品| 久久国产精品影院| 亚洲欧美一区二区三区黑人| 一区二区三区国产精品乱码| 九九久久精品国产亚洲av麻豆| 老司机午夜十八禁免费视频| 两个人的视频大全免费| 色播亚洲综合网| 天堂影院成人在线观看| 欧美大码av| 神马国产精品三级电影在线观看| 麻豆国产97在线/欧美| 成人无遮挡网站| 99久久99久久久精品蜜桃| 午夜福利欧美成人| 女警被强在线播放| 国产高潮美女av| 成年人黄色毛片网站| 欧美区成人在线视频| 国产成人影院久久av| 网址你懂的国产日韩在线| 俄罗斯特黄特色一大片| 欧美一级a爱片免费观看看| 亚洲av美国av| 在线播放国产精品三级| 日本免费一区二区三区高清不卡| 午夜免费激情av| 在线观看免费视频日本深夜| 午夜老司机福利剧场| 全区人妻精品视频| 99riav亚洲国产免费| 99久久九九国产精品国产免费| 欧美激情在线99| www日本在线高清视频| av视频在线观看入口| www日本在线高清视频| 日本a在线网址| 久久久久久久久大av| 成人18禁在线播放| 他把我摸到了高潮在线观看| 久久久久精品国产欧美久久久| 男女午夜视频在线观看| aaaaa片日本免费| 亚洲av成人不卡在线观看播放网| 欧美日韩一级在线毛片| 成人av一区二区三区在线看| 国产亚洲av嫩草精品影院| 精品国产亚洲在线| 别揉我奶头~嗯~啊~动态视频| 午夜两性在线视频| 午夜福利在线观看吧| 在线十欧美十亚洲十日本专区| 色综合站精品国产| 啦啦啦观看免费观看视频高清| 国产私拍福利视频在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费搜索国产男女视频| 波多野结衣巨乳人妻| 岛国视频午夜一区免费看| 法律面前人人平等表现在哪些方面| 国产三级中文精品| 国产淫片久久久久久久久 | 亚洲人成网站在线播放欧美日韩| 最新中文字幕久久久久| xxxwww97欧美| 国产爱豆传媒在线观看| 国产精华一区二区三区| 久久国产精品影院| 日韩高清综合在线| 网址你懂的国产日韩在线| 一个人免费在线观看电影| 国产亚洲精品av在线| 淫秽高清视频在线观看| 亚洲久久久久久中文字幕| 久久精品国产99精品国产亚洲性色| 国产高潮美女av| 免费看十八禁软件| 特级一级黄色大片| 超碰av人人做人人爽久久 | 亚洲专区国产一区二区| 一个人观看的视频www高清免费观看| 欧美日韩国产亚洲二区| 九九热线精品视视频播放| 窝窝影院91人妻| 日本精品一区二区三区蜜桃| 99热这里只有是精品50| 五月玫瑰六月丁香| 真实男女啪啪啪动态图| 搡女人真爽免费视频火全软件 | 国产午夜精品论理片| 国产av不卡久久| 国产欧美日韩精品一区二区| 变态另类丝袜制服| av专区在线播放| 深夜精品福利| 99热这里只有是精品50| 19禁男女啪啪无遮挡网站| 脱女人内裤的视频| 国产精品免费一区二区三区在线| 色尼玛亚洲综合影院| 国产主播在线观看一区二区| 亚洲第一电影网av| 高清日韩中文字幕在线| 观看免费一级毛片| АⅤ资源中文在线天堂| 麻豆一二三区av精品| 久久精品影院6| 国产探花在线观看一区二区| 国产av不卡久久| 国产激情欧美一区二区| 女生性感内裤真人,穿戴方法视频| 脱女人内裤的视频| 国产av不卡久久| 99久久精品一区二区三区| 90打野战视频偷拍视频| 国产成人a区在线观看| 国产成年人精品一区二区| 成人18禁在线播放| 免费电影在线观看免费观看| 国产精品久久电影中文字幕| 人人妻人人澡欧美一区二区| 少妇的丰满在线观看| 欧美在线黄色| 国产色婷婷99| 丝袜美腿在线中文| av在线蜜桃| 欧美日韩亚洲国产一区二区在线观看| 九色成人免费人妻av| 白带黄色成豆腐渣| 国产中年淑女户外野战色| 国产精品av视频在线免费观看| 一级毛片女人18水好多| 亚洲黑人精品在线| 波多野结衣高清无吗| 国内久久婷婷六月综合欲色啪| 色综合站精品国产| 国产精品亚洲一级av第二区| 美女高潮的动态| 欧美黑人欧美精品刺激| 久久久久国产精品人妻aⅴ院| 成人性生交大片免费视频hd| 国产精品av视频在线免费观看| 在线国产一区二区在线| 久久久久九九精品影院| 亚洲精品日韩av片在线观看 | 一a级毛片在线观看| 脱女人内裤的视频| 天天躁日日操中文字幕| 婷婷精品国产亚洲av在线| 久久人妻av系列| 一夜夜www| ponron亚洲| 在线播放国产精品三级| 免费看美女性在线毛片视频| 动漫黄色视频在线观看| 老司机午夜福利在线观看视频| 床上黄色一级片| 少妇高潮的动态图| 亚洲人成伊人成综合网2020| 99热只有精品国产| 岛国在线观看网站| 午夜激情福利司机影院| 亚洲av美国av| 婷婷亚洲欧美| 久久久久性生活片| 国产成年人精品一区二区| 久久久久国产精品人妻aⅴ院| 黄色片一级片一级黄色片| 国产91精品成人一区二区三区| 国产探花极品一区二区| 熟女人妻精品中文字幕| 久久精品国产亚洲av香蕉五月| 亚洲无线观看免费| 国产淫片久久久久久久久 | 国产亚洲欧美在线一区二区| 国产精品亚洲一级av第二区| 蜜桃亚洲精品一区二区三区| 99精品欧美一区二区三区四区| 少妇的逼水好多| 欧美激情在线99| 精品国产超薄肉色丝袜足j| 三级毛片av免费| av天堂在线播放| 国产亚洲欧美98| 亚洲 国产 在线| 日韩精品青青久久久久久| 久久久色成人| 亚洲真实伦在线观看| 99久国产av精品| 18禁国产床啪视频网站| 午夜免费男女啪啪视频观看 | 毛片女人毛片| 国产精品久久久久久亚洲av鲁大| 欧美成狂野欧美在线观看| 精品电影一区二区在线| 夜夜夜夜夜久久久久| 国产欧美日韩一区二区三|