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

    汽液兩相二維壓力波傳播程序開發(fā)及驗證

    2023-06-19 09:31:10龔政宇辜峙钘潘麒文張牧昊戴嘉寧
    核技術 2023年6期
    關鍵詞:液態(tài)水傳質(zhì)蒸汽

    龔政宇 辜峙钘,2 潘麒文 張牧昊 戴嘉寧

    1(成都理工大學 核技術與自動化工程學院 成都 610059)

    2(中國核動力研究設計院 核反應堆系統(tǒng)設計技術重點實驗室 成都 610213)

    壓力波傳播特性是反應堆事故安全分析的關注點之一,反應堆系統(tǒng)內(nèi)壓力波傳播將引發(fā)水力學載荷問題,獲得壓力波的傳播數(shù)據(jù)是開展結構載荷分析、評價的必要條件[1]。核能系統(tǒng)反應堆內(nèi)的壓力波傳播現(xiàn)象一般由瞬態(tài)事故引發(fā),例如,蒸汽輪機/主蒸汽管閥突開、突閉,或諸如冷卻劑喪失(Loss Of Coolant Accident,LOCA),蒸汽發(fā)生器傳熱管破裂(Steam Generator Tube Rupture accident,SGTR)[2]等事故。大部分壓力波傳播數(shù)值模擬采用可壓縮歐拉流體動力學模型[3]。而對于兩相流動、傳熱、傳質(zhì)條件下的壓力波傳播問題,通常選用兩流體六方程模型,被系統(tǒng)分析程序廣泛使用,如RELAP5,TRACE,WAHA[4-5]。Sokolowski等[6]針對典型系統(tǒng)分析程序RELAP5 模擬壓力波傳播的能力開展了分析與評價,包括單相氣體介質(zhì)、單相蒸汽介質(zhì)及汽水兩相介質(zhì)條件下壓力波傳播問題。Barna和Serre使用蒸汽冷凝引發(fā)的水錘效應及壓力波傳播案例,對WAHA3 和CATHARE 程序也進行了分析與評價[7-8]。另外,兩流體六方程模型一致采用平衡態(tài)壓力模型,即氣液兩相的動量守恒方程共用同一壓力,針對該問題,相關研究機構已開始研發(fā)新一代反應堆系統(tǒng)分析程序,如RELAP7,其核心便是雙壓力七方程模型。Delchini等[9]基于單相、兩相激波管基準例題,分析、評價了RELAP7程序捕捉壓力波傳播行為的能力,結果表明:在壓力波傳播模擬方面,該程序的確優(yōu)于RELAP5程序。

    壓力波傳播行為通常存在多維特性,幾何結構越復雜,該特征越顯著,例如SGTR事故下壓力波在傳熱管束范圍內(nèi)的復雜傳播、反射及疊加行為,而前述系統(tǒng)分析程序通常屬于一維程序,無法處理壓力波的多維傳播問題,因此開展壓力波多維傳播模擬技術、方法研究具有重要意義。針對該問題,本文就汽、水兩相流中存在的壓力波二維軸對稱傳播行為,開展了數(shù)值模型與算法研究,控制方程采用二維軸對稱圓柱坐標系,考慮了汽、水兩相非平衡態(tài)相變傳熱、傳質(zhì)行為,控制方程求解使用了時間步分離四步算法,開發(fā)了程序代碼TPFPWPC-2D (Two Dimensional Two Phase Flow Pressure Wave Propagation Code)。采用典型汽、水兩相激波管基準例題開展了程序驗證,同時為考查該程序模擬壓力波二維傳播特性的能力,對圓柱體空間內(nèi)壓力波的軸對稱二維傳播特性進行了數(shù)值模擬與分析研究。

    1 模型和算法

    1.1 控制方程

    本文涉及的控制方程主要包括質(zhì)量守恒方程、動量守恒方程、能量守恒方程,且全部控制方程均基于二維軸對稱圓柱體坐標系導出。

    1) 質(zhì)量守恒方程

    就質(zhì)量組分而言,考慮三種不同質(zhì)量組分,即液態(tài)水、水蒸氣、不可凝氣體,其質(zhì)量守恒方程如式(1)~(3)所示。其中:ΓC、ΓV為相變質(zhì)量源項;αL為液態(tài)水體積份額;αG為蒸汽混合物(水蒸氣和不可凝氣體)體積份額;ρˉL、ρˉGh2o、ρˉGgas分別是液態(tài)水、水蒸氣、不可凝氣體宏觀密度;ρL、ρGh2o、ρGgas為相應理論密度;uL、uG分別是液態(tài)水、蒸汽混合物徑向速度;νL、νG為相應軸向速度。可以看出水蒸氣和不可凝結性氣體共享同一速度場。

    2) 動量守恒方程

    本文考慮了兩類速度場組分,即液態(tài)水速度場和蒸汽混合物速度場,其徑向動量守恒方程由式(4)、(5)描述,軸向動量守恒方程如式(6)、(7)所示。其中:p為網(wǎng)格壓強;FDR,G,R、FDR,L,R表示氣、液兩相徑向曳力;FDR,G,Z、FDR,L,Z表示氣、液兩相軸向拽曳力;FVM,G,R、FVM,L,R為 氣、液 兩 相 徑 向 虛 擬 質(zhì) 量 力;FVM,G,Z、FVM,L,Z表示氣、液兩相軸向虛擬質(zhì)量力;ρˉG為蒸汽混合物宏觀密度,表示為ρˉG=ρˉGh2o+ρˉGgas。

    3) 能量守恒方程

    考慮兩種能量組分,即液態(tài)水、蒸汽混合物,其能量守恒方程如式(8)、(9)所示。其中:QΓ,G、QΓ,L為非平衡態(tài)相變傳熱;QG、QL為平衡態(tài)溫差傳熱;a為接觸面積;h為換熱系數(shù);ΔT為溫差。由于本文模擬對象涉及強傳熱、傳質(zhì),而且對流效應也非常顯著,因此忽略了網(wǎng)格間熱傳導的貢獻,也忽略了相間相互作用力做功的貢獻。

    1.2 數(shù)值算法

    本文模擬對象具有強傳熱、傳質(zhì),源項變化劇烈,壓力波傳播急劇等典型特征,同時其“速度-壓強”耦合尤為緊密,基于完全統(tǒng)一、全隱式求解算法對本模擬對象而言不太現(xiàn)實,主要是收斂性難以保障。針對此問題,本文采用了時間步分割四步算法技術,以提高其數(shù)值計算收斂性與穩(wěn)定性。

    1) 第一步:網(wǎng)格內(nèi)源相計算

    首先,聯(lián)立迭代求解水蒸氣質(zhì)量守恒方程,液態(tài)水、蒸汽混合物能量守恒方程,僅考慮傳熱、傳質(zhì)源項,其核心是基于汽、液兩相狀態(tài)方程(Equation Of State,EOS)的非平衡態(tài)傳熱、傳質(zhì)迭代計算,迭代收斂后,進而更新各組分宏觀密度、比內(nèi)能、溫度及傳質(zhì)速率。然后,根據(jù)傳質(zhì)速率單獨求解各組分動量守恒方程,僅考慮其相變源項,由于當前“速度-壓強”并不匹配,因此相變源項中速度采用顯式格式,最終更新速度場。上述求解過程整體邏輯框架體系如圖1所示,其中上標1 表示第一步值,上標n表示前時間步值。

    圖1 第一步:網(wǎng)格內(nèi)源項計算數(shù)值算法框架體系Fig.1 Step 1: numerical algorithm framework for cell-inside source term calculation

    2) 第二步:網(wǎng)格間對流計算

    首先,單獨求解全部質(zhì)量守恒方程,只考慮對流貢獻,且采用顯式格式;然后,單獨求解全部能量守恒方程,只考慮對流貢獻,且也采用顯式格式;最后,聯(lián)立迭代求解動量守恒方程,僅忽略相變源項,其中壓力項、對流項均采用顯式格式,拽曳力、虛擬質(zhì)量力均采用隱式格式。由于在第2步中,“速度-壓強”依然不匹配,因此其對流項、壓強均采用顯式格式。上述求解過程整體邏輯框架體系如圖2 所示,其中上標2 表示第二步值,上標0 表示前時間步值。

    圖2 第二步:網(wǎng)格間對流計算數(shù)值算法框架體系Fig.2 Step 2: numerical algorithm framework for intercell convection calculation

    3) 第三步:壓力迭代

    首先,確定殘差方程,本文共考慮4個,第1個是壓力殘差(EOS 壓力和網(wǎng)格壓力之差),第2、3 分別是液態(tài)水、蒸汽混合物質(zhì)量守恒方程殘差,其對流項采用半隱格式,第4 個是蒸汽混合物能量守恒方程殘差,其比內(nèi)能對流項采用顯式格式(與膨脹功相比,忽略對流效應影響),空泡份額對流項采用半隱格式。然后,將上述守恒方程殘差對p、pe、pw、pn、ps、ρˉL、ρˉG、TG進行泰勒展開,即可獲得迭代方程組。最后,充分利用狀態(tài)方程EOS 模型,迭代求解即可更新網(wǎng)格壓強,進而更新兩相速度、各組分宏觀密度、各組分比內(nèi)能等其他狀態(tài)參數(shù)。上述求解過程整體邏輯框架體系如圖3 所示,其中上標3 表示第三步值。

    圖3 第三步:壓力迭代算法框架體系Fig.3 Step three: pressure iterative algorithm framework

    4) 第四步:時間步末更新

    重復第二步具體過程,但其網(wǎng)格壓強、對流項速度將采用第三步最新結果,最終更新宏觀密度、比內(nèi)能、速度,即。其具體算法框架與第二步極為相似,限于篇幅不再贅述,其中上標n+ 1 表示當前時間步值。

    1.3 輔助模型

    1) 兩相流流型模型

    兩相流流型決定了各相組分間的接觸方式,繼而影響交界面面積、曳力系數(shù)、傳熱系數(shù)及傳熱、傳質(zhì)計算。文中考慮三種不同流型,包括泡狀流、彌散流和過渡流,具體流型由蒸汽相體積份額判斷。當αG≤0.3時,視為純泡狀流,當αG≥0.7時,視為純彌散流,當0.3 <αG<0.7 時,視為過渡流。對于控制體網(wǎng)格而言,可統(tǒng)一視為由泡狀流區(qū)域、彌散流區(qū)域疊加構成,并采用泡狀流區(qū)域、彌散流區(qū)域體積份額fB、fD進行具體表征,其中fB、fD按式(10)計算,其中αB、αD分別是泡狀流、彌散流特征值,分別默認取0.3、0.7。xB、xD為組分在泡狀流、彌散流區(qū)域內(nèi)的分配系數(shù),具體按式(11)計算。則泡狀流區(qū)域內(nèi)液態(tài)水(連續(xù)相)體積份額可表示為αL,B=αLxB,彌散流區(qū)域內(nèi)液態(tài)水(液滴)體積份額可表示為αL,D=αLxD。蒸汽混合物在泡狀流區(qū)域、彌散流區(qū)域內(nèi)的體積份額按式(12)、(13)計算。

    泡狀流區(qū)域內(nèi)水蒸氣(氣泡)體積份額為:

    彌散流區(qū)域內(nèi)水蒸氣(連續(xù)相)體積份額為:

    2) 相界面輸運模型

    為捕捉在流動,傳熱、傳質(zhì)等過程中相界面演化行為,考慮了相界面輸運方程,具體包括泡狀流區(qū)域內(nèi)蒸汽泡和彌散流區(qū)域內(nèi)液滴的相界面輸運行為,其輸運方程分別如下,其中:aG,B、aD,D分別是泡狀流區(qū)域內(nèi)氣泡相界面積、彌散流區(qū)域內(nèi)液滴相界面積。

    就本文模擬而言,泡狀流區(qū)域內(nèi)蒸汽泡相界面輸運考慮了對流、膨脹/壓縮、成核及破裂行為;彌散流區(qū)域內(nèi)液滴僅考慮了對流效應。其中,氣泡成核即式(14)中源項SG,B,N,具體模型參考了先進流體動力學程序AFDM[10],具體采用平衡態(tài)模型,其平衡態(tài)參數(shù)aeG,B,N、時間常數(shù)τN及源項SG,B,N按式(16)計算。蒸汽泡破裂模型采用了韋伯破裂模型,也使用平衡態(tài)方法計算,其平衡態(tài)參數(shù)aeG,B,W、時間常數(shù)τW及源項SG,B,W按式(17)計算。

    3) 相間曳力模型

    在泡狀流區(qū)域,“液態(tài)水-水蒸氣”接觸方式為連續(xù)液體包裹氣泡,在彌散流區(qū)域,其接觸方式為連續(xù)蒸汽包裹液滴。因此,總體上“液態(tài)水-水蒸氣”接觸方式可按“連續(xù)相-彌散相”處理。下文以連續(xù)相對彌散項施加的徑向曳力進行闡述,正如式(18)所示,其中:KDR,D,R為“連續(xù)相-彌散相”曳力系數(shù);aCD為交界面積;uC、uD分別表示連續(xù)相和彌散相徑向速度;μC為連續(xù)相動力黏度;ρC為連續(xù)相密度;rD為彌散相半徑。系數(shù)CD參考Ishii發(fā)展的相關模型[11]。

    4) 虛擬質(zhì)量力模型

    在兩相流中,當彌散相(如氣泡、液滴)相對連續(xù)相加速時,將產(chǎn)生“虛擬質(zhì)量效應”,由于本文需要模擬“汽-水”兩相流壓力波傳播,因此考慮了虛擬質(zhì)量力模型。具體參考了Fullmer和Bertodano報道的相關模型[12],系統(tǒng)分析程序RELAP5/MOD3.3 也使用了該模型。

    液態(tài)水所受虛擬質(zhì)量力為:

    式中:αL為液態(tài)水體積份額;ρL為液態(tài)水理論密度;αG為蒸汽混合物體積份額;注意到上述虛擬質(zhì)量力采用的是矢量形式,CVM為常數(shù),通常取0.5。

    蒸汽混合物虛擬質(zhì)量力為:

    5) 傳熱傳質(zhì)模型

    就本模擬而言,考慮了“汽-水”兩相傳熱、傳質(zhì)模型,尤其是交界面相變傳熱模型。

    a) 傳質(zhì)模型

    傳質(zhì)考慮的是交界面相變傳質(zhì),“液態(tài)水-水蒸氣”交界面凈傳熱速率可表示為:

    式中:aL,G表示“液態(tài)水-蒸汽混合物”交界面積;RL,Gh2o是“液態(tài)水-水蒸氣”接觸比例;hL為交界面之液態(tài)水側(cè)換熱系數(shù);hG為交界面之蒸汽混合物側(cè)換熱系數(shù);TLi,Gh2o為“液態(tài)水-水蒸氣”交界面溫度,取飽和溫度與平衡態(tài)傳熱交界面溫度之較大者。qiL,Gh2o>0表示交界面存在凈熱量傳出,則交界面發(fā)生了水蒸氣凝結現(xiàn)象,凝結速率按式(22)計算;qiL,Gh2o<0表示交界面存在凈熱量傳入,則交界面發(fā)生了液態(tài)水汽化現(xiàn)象,汽化速率按式(23)計算;qiL,Gh2o= 0 表示交界面處于傳熱平衡態(tài),無相變發(fā)生。

    b)傳熱模型

    就文中模擬對象而言,包括非平衡態(tài)相變傳熱和平衡態(tài)溫差傳熱。液態(tài)水、水蒸氣的非平衡態(tài)相變傳熱分別按式(24)、(25)計算。

    液態(tài)水、水蒸氣的平衡態(tài)溫差傳熱分別按式(26)、(27)計算。

    式中:RL,Ggas表示“液態(tài)水-不可凝氣體”接觸比例;TiL,Ggas表示“液態(tài)水-不可凝氣體”交界面溫度,取其平衡態(tài)溫度。

    6)狀態(tài)方程EOS模型

    狀態(tài)方程EOS 是封閉可壓縮兩相流守恒方程的必備條件,同時也是本文開展傳熱、傳質(zhì)數(shù)值迭代,壓力迭代的重要保障。本文模擬涉及液態(tài)水、水蒸氣、不可凝氣體三種物質(zhì)組分的狀態(tài)方程。

    a)水蒸氣狀態(tài)方程

    水蒸氣狀態(tài)方程采用了MRK(Modified Redlich-Kwong)模型,下式給出了水蒸氣分壓與其溫度和比體積之間的關系。

    其中:x(TG)按式(29)計算;Th2o,crt為水臨界溫度。

    b)不可凝氣體狀態(tài)方程

    不可凝氣體采用理想氣體狀態(tài)方程模型。

    c)液態(tài)水狀態(tài)方程

    本文考慮液態(tài)水為可壓縮流體,式(30)給出了液態(tài)水溫度與其比內(nèi)能和壓強間的關系,其中:TS,L為液態(tài)水飽和溫度,采用液態(tài)水比內(nèi)能計算;?TL?p為液態(tài)水溫度對壓強偏導數(shù),也使用液態(tài)水比內(nèi)能計算。

    類似地,式(31)給出了液態(tài)水比體積與其比內(nèi)能和壓強之間的關系,形式與式(30)類似。

    2 程序驗證

    為驗證本文開發(fā)的汽-水兩相二維軸對稱壓力波傳播模擬程序TPFPWPC-2D,選用了典型的汽-水兩相激波管基準例題,如圖4所示,雖然該激波管屬于一維問題,但可以驗證該程序的關鍵模型與算法。該激波管長100 m,中心插入一隔板,將激波管劃分為兩個不同區(qū)域,左側(cè)為15 MPa的汽水混合物高壓區(qū),右側(cè)為10 MPa 的汽水混合物低壓區(qū),具體初始參數(shù)如圖4 所示。由于該程序具有二維軸對稱特性,恰可建立上述激波管,具體在徑向劃分3 個網(wǎng)格,軸向劃分200 個網(wǎng)格,激波管幾何模型豎直放置,則需忽略重力與基準例題保持一致。另外,建模時需要考慮到徑向均勻性,即壓力波僅沿軸向傳播。激波管上、下兩端采用絕對反射邊界條件。很明顯,汽-水兩相流激波管基準例題并無解析解,因此采用了程序?qū)Ρ确绞竭M行驗證,具體選用了典型系統(tǒng)分析程序RELAP5 公開發(fā)布的模擬結果作為對比基準[6]。

    圖4 兩相激波管示意圖Fig.4 Schematic diagram of two-phase shock tube

    TPFPWPC-2D 程序在方程求解上采用了時間步分割四步算法,即人為將控制方程中傳熱、傳質(zhì)、相變所致動量交換分離處置,這極可能造成時間步長敏感性問題,針對此問題,首先開展了時間步長敏感性分析,共選用15 個時間步長,即5 000 μs、2 500 μs、1 000 μs、750 μs、500 μs、250 μs、100 μs、75 μs、50 μs、25 μs、10 μs、7.5 μs、5 μs、2.5 μs、1 μs,對此基波管進行模擬。圖5給出了不同時間步長下此激波管中心位置0.081 s時刻的壓強值,很明顯時間步長在1~100 μs范圍內(nèi)具有良好的數(shù)值穩(wěn)定性。同時,測試表明若時間步長大于5 ms,計算則發(fā)散(庫朗限制)。

    圖5 汽水兩相激波管TPFPWPC-2D時間步長敏感性分析Fig.5 Time-step sensitivity analysis of TPFPWPC-2D based on water-vapor two-phase shock tube

    根據(jù)上述時間步長敏感性分析結果,最終選用了10-5s 作為該驗證用時間步長,且上述文獻中RELAP5 也是采用了此步長。圖6 分別給出了0.081 s 時刻壓力、汽水兩相溫度、速度、空泡體積份額的分布情況。從圖6(a)可知,隨著隔膜突然撤離,汽-水兩相高壓區(qū)內(nèi)形成沖擊波并傳入汽-水兩相低壓區(qū),致使低壓區(qū)內(nèi)壓力上升;相反汽-水兩相低壓區(qū)內(nèi)形成膨脹波并傳入汽-水兩相高壓區(qū),致使高壓區(qū)內(nèi)壓力急劇下降。從圖6(b)可知,汽-水兩相溫度分布與壓力分布具有一定相似性。從圖6(c)可知,激波管中心兩側(cè)區(qū)域出現(xiàn)了汽-水兩相高速運動,最大速度達到約48 m·s-1。從圖6(d)可知,由于汽-水兩相膨脹、汽化效應,高壓區(qū)偏中心區(qū)域空泡份額出現(xiàn)了急劇上升,可見明顯極大值點;相反,由于汽-水兩相壓縮、凝結效應,低壓區(qū)偏中心位置空泡份額出現(xiàn)了急劇下降,可見明顯極小值點。從壓力、溫度、速度及空泡份額分布特性比較可以看出,TPFPWPC-2D程序與系統(tǒng)分析程序RELAP5模擬結果整體符合較好,但TPFPWPC-2D 與RELAP5 依然存在局部差異,很可能是由兩程序所考慮的模型精細度或個別模型差異造成的,也可能與數(shù)值算法有關系。

    圖6 兩相基波管壓力分布(a)、溫度分布(b)、速度分布(c)和空泡份額分布(d)Fig.6 Pressure profile (a), temperature profile (b), velocity profile (c) and void fraction profile (d) of the two-phase shock tube

    3 壓力波二維軸對稱傳播特性研究

    為研究圓柱體區(qū)域內(nèi)汽-水兩相流動、傳熱、傳質(zhì)條件下壓力波二維軸對稱傳播特性,使用文中開發(fā)的TPFPWPC-2D 程序開展了數(shù)值模擬研究。數(shù)值模擬幾何模型如圖7 所示,建??紤]了徑向和軸向(即對稱軸)兩方向,實際模擬的幾何區(qū)域為一個高1.0 m,半徑為1.0 m 的圓柱體。徑向0~0.2 m、軸向0.1~ 0.3 m 為高壓區(qū)域,其余為低壓區(qū)域,且初始時刻兩區(qū)域處于相互隔離狀態(tài),數(shù)值模擬開始后,兩區(qū)域流體瞬間接觸。

    圖7 二維壓力波傳播實例示意圖Fig.7 Schematic of 2D pressure wave case

    首先,針對該問題開展了網(wǎng)格敏感性分析,選用了20×20、40×40、50×50、60×60、70×70、80×80、90×90、100×100、110×110、120×120、130×130、140×140、150×150,共計13套不同網(wǎng)格方案,均考慮等間距劃分,并以8.0×10-4s時刻計算區(qū)域內(nèi)最大壓強值為觀測點,對網(wǎng)格敏感性進行評價。圖8即為評價結果,很明顯,當網(wǎng)格劃分方案大于100×100時,數(shù)值模擬結果已具較好穩(wěn)定性,據(jù)此下文模擬選用了120×120 網(wǎng)格劃分方案,開展壓力波二維傳播模擬與特性分析。

    圖8 二維壓力波傳播網(wǎng)格敏感性分析Fig.8 Mesh sensitivity analysis of 2D pressure wave propagation

    圖9 給出了壓力波的二維傳播過程,注意該圖已進行對稱性拓撲。從數(shù)值模擬結果可以看出,隨著兩個不同區(qū)域內(nèi)流體被解除隔離,來自高壓區(qū)域的流體向低壓區(qū)域射流,進而產(chǎn)生了向外快速傳播的沖擊波(shock wave),致使外部低壓區(qū)域內(nèi)壓力劇增。同時,來自低壓區(qū)域內(nèi)的膨脹波(rarefaction wave)也隨之產(chǎn)生,并傳播進入高壓區(qū)域,致使其內(nèi)部壓力降低。隨著時間推移,壓力波不斷向邊界傳播,約0.4 ms 后,壓力波已抵達底部邊界,被邊界壁面完全反射,進而產(chǎn)生反向沖擊波,反向沖擊波與迎面沖擊波隨即相遇,繼而出現(xiàn)了疊加效應,致使下邊界附近出現(xiàn)了局部高壓區(qū)。在下一個階段,隨著壓力波不斷向外擴散,沖擊波逐漸減弱。

    圖9 二維壓力波傳播實例數(shù)值模擬結果Fig.9 The simulation results of 2D pressure wave propagation case

    4 結語

    本文基于汽-水兩相歐拉流體動力學理論,開發(fā)了一種基于二維軸對稱圓柱坐標系的壓力波傳播模擬二維程序TPFPWPC-2D,考慮了汽-水兩相非平衡態(tài)相變傳熱、傳質(zhì),偏微分方程求解使用了時間步分離四步算法。并采用了典型的汽水兩相激波管基準例題開展了程序的驗證工作及時間步長敏感性分析,數(shù)值模擬結果與典型的系統(tǒng)分析程序RELAP5進行了對比。最后為考查該程序模擬壓力波二維傳播的能力,對圓柱體空間區(qū)域內(nèi)的壓力波傳播行為進行了二維軸對稱數(shù)值模擬,同時還開展了網(wǎng)格敏感性分析。本文研究得到的結論如下:

    1)針對汽-水兩相激波管基準例題,TPFPWPC-2D 程序計算獲得的壓力、溫度、速度及空泡份額分布數(shù)據(jù)與典型的系統(tǒng)分析程序RELAP5 符合較好,程序總體上可以很好地模擬強傳熱傳質(zhì)兩相條件下壓力波傳播行為,進而從整體上驗證了該程序的結構模型及數(shù)值算法的正確性。同時,從時間步長敏感性分析結果得知,本文采用的時間分割四步算法具有較為顯著的時間步長敏感性。

    2)針對圓柱體空間內(nèi)壓力波二維軸對稱傳播算例,采用TPFPWPC-2D 程序開展了數(shù)值模擬,獲得壓力波在圓柱體內(nèi)的二維軸對稱傳播特性,尤其是捕捉到了壓力波在二維傳播過程中反射、疊加行為,進而從一定層面上確認了該程序模擬壓力波二維傳播問題的能力。

    作者貢獻聲明龔政宇實施研究,分析/解釋數(shù)據(jù),起草撰寫文章;辜峙钘提供研究思路,指導論文寫作與修改,提供研究經(jīng)費支持;潘麒文負責材料支持,提供相關文獻;張牧昊負責文章內(nèi)容的審閱指導;戴嘉寧參與文獻收集整理工作,數(shù)據(jù)整理。

    猜你喜歡
    液態(tài)水傳質(zhì)蒸汽
    核電廠蒸汽發(fā)生器一次側(cè)管嘴堵板研發(fā)和應用
    基于微波輻射計的張掖地區(qū)水汽、液態(tài)水變化特征分析
    Ka/Ku雙波段毫米波雷達功率譜數(shù)據(jù)反演液態(tài)水含量方法研究
    零下溫度的液態(tài)水
    PEMFC氣體擴散層中液態(tài)水傳輸實驗研究綜述
    電源技術(2017年10期)2017-11-09 05:22:21
    一種新型蒸汽發(fā)生器結構設計
    中國核電(2017年1期)2017-05-17 06:10:04
    第一艘蒸汽輪船
    蒸汽閃爆
    氨基酸鹽吸收二氧化碳過程的傳質(zhì)特性
    PTFE膜吸收CO2的工藝參數(shù)對傳質(zhì)性能的影響
    АⅤ资源中文在线天堂| av黄色大香蕉| 国产午夜精品论理片| 中文亚洲av片在线观看爽| 亚洲精品日韩av片在线观看| 最近中文字幕高清免费大全6 | 亚洲国产色片| 国内精品久久久久精免费| 三级男女做爰猛烈吃奶摸视频| 国产精品影院久久| 黄片小视频在线播放| 日本免费a在线| 亚洲最大成人手机在线| 尤物成人国产欧美一区二区三区| 国产毛片a区久久久久| 99国产极品粉嫩在线观看| 色哟哟·www| 观看美女的网站| 熟女人妻精品中文字幕| 9191精品国产免费久久| 日本与韩国留学比较| 亚洲国产欧美人成| 99热6这里只有精品| 国产探花在线观看一区二区| 一进一出抽搐gif免费好疼| 99热只有精品国产| 日韩亚洲欧美综合| 国产色爽女视频免费观看| 变态另类成人亚洲欧美熟女| 色吧在线观看| 亚洲欧美日韩高清专用| 亚洲av中文字字幕乱码综合| 欧美又色又爽又黄视频| 免费大片18禁| 欧美另类亚洲清纯唯美| 免费av不卡在线播放| 少妇被粗大猛烈的视频| 女人十人毛片免费观看3o分钟| xxxwww97欧美| 搡老熟女国产l中国老女人| 日韩欧美国产一区二区入口| 舔av片在线| 精品一区二区三区视频在线| a级一级毛片免费在线观看| 永久网站在线| 亚洲无线在线观看| 欧美丝袜亚洲另类 | 12—13女人毛片做爰片一| 日韩欧美三级三区| 欧美日韩福利视频一区二区| 免费av观看视频| 亚洲 国产 在线| 久久久久久久午夜电影| 亚洲成a人片在线一区二区| 看黄色毛片网站| 在线观看av片永久免费下载| a级毛片a级免费在线| 欧美xxxx性猛交bbbb| 真人做人爱边吃奶动态| 亚洲真实伦在线观看| 热99在线观看视频| 精品免费久久久久久久清纯| 欧美黄色淫秽网站| 天堂网av新在线| 国产一级毛片七仙女欲春2| 久久精品国产清高在天天线| 1024手机看黄色片| 免费av观看视频| 婷婷精品国产亚洲av在线| 亚洲最大成人av| 老司机午夜十八禁免费视频| 成人高潮视频无遮挡免费网站| 精品人妻视频免费看| 一级毛片久久久久久久久女| 国产成人啪精品午夜网站| 国语自产精品视频在线第100页| 一区福利在线观看| 欧美乱妇无乱码| 人妻丰满熟妇av一区二区三区| 永久网站在线| 女生性感内裤真人,穿戴方法视频| 日韩有码中文字幕| 国产精品亚洲美女久久久| 欧美日本亚洲视频在线播放| 狠狠狠狠99中文字幕| 成熟少妇高潮喷水视频| 精品国产亚洲在线| 午夜福利免费观看在线| 色av中文字幕| 日韩中文字幕欧美一区二区| 欧美性猛交黑人性爽| 老司机午夜十八禁免费视频| 网址你懂的国产日韩在线| 九九久久精品国产亚洲av麻豆| 国模一区二区三区四区视频| 欧美最黄视频在线播放免费| 国产午夜精品久久久久久一区二区三区 | netflix在线观看网站| 国产精品影院久久| 丝袜美腿在线中文| 怎么达到女性高潮| 亚洲美女黄片视频| 欧美成人一区二区免费高清观看| 欧美潮喷喷水| 一级a爱片免费观看的视频| av视频在线观看入口| 国产精品永久免费网站| 中文字幕av成人在线电影| 亚洲av电影在线进入| 美女高潮的动态| 亚洲国产精品合色在线| 国产高清视频在线观看网站| 两人在一起打扑克的视频| 特级一级黄色大片| 黄片小视频在线播放| 亚洲欧美日韩高清专用| 日韩国内少妇激情av| 草草在线视频免费看| 日韩高清综合在线| 中文字幕av在线有码专区| av中文乱码字幕在线| 男女之事视频高清在线观看| 成人美女网站在线观看视频| 免费在线观看影片大全网站| 欧美成狂野欧美在线观看| 国产aⅴ精品一区二区三区波| 夜夜看夜夜爽夜夜摸| 亚洲内射少妇av| 一级黄片播放器| 熟妇人妻久久中文字幕3abv| 丰满人妻一区二区三区视频av| 欧美乱妇无乱码| 午夜视频国产福利| 久久精品久久久久久噜噜老黄 | 亚洲精品456在线播放app | 草草在线视频免费看| 桃红色精品国产亚洲av| 亚洲三级黄色毛片| 午夜影院日韩av| 超碰av人人做人人爽久久| 国产精品久久久久久久久免 | 一个人免费在线观看电影| 亚洲精品456在线播放app | 色吧在线观看| 亚洲欧美激情综合另类| 小说图片视频综合网站| 国产一区二区亚洲精品在线观看| 成熟少妇高潮喷水视频| 日本免费一区二区三区高清不卡| 亚洲七黄色美女视频| 精品人妻一区二区三区麻豆 | 日日摸夜夜添夜夜添小说| av在线观看视频网站免费| 免费搜索国产男女视频| 午夜免费男女啪啪视频观看 | 日本五十路高清| 女人十人毛片免费观看3o分钟| 国产又黄又爽又无遮挡在线| 最好的美女福利视频网| 国产精华一区二区三区| 国产探花极品一区二区| 国产av不卡久久| 最新在线观看一区二区三区| 欧美日韩综合久久久久久 | 国产伦人伦偷精品视频| 一本精品99久久精品77| 最好的美女福利视频网| 午夜福利视频1000在线观看| 内射极品少妇av片p| 久久久成人免费电影| 最近在线观看免费完整版| 在线播放无遮挡| 男人舔女人下体高潮全视频| 三级毛片av免费| 激情在线观看视频在线高清| 淫秽高清视频在线观看| 黄色丝袜av网址大全| 成人毛片a级毛片在线播放| 国内毛片毛片毛片毛片毛片| 国产在线精品亚洲第一网站| 久久人人精品亚洲av| 欧美最黄视频在线播放免费| 欧美日韩综合久久久久久 | 制服丝袜大香蕉在线| 日本 欧美在线| 99热6这里只有精品| 国产亚洲精品综合一区在线观看| 久久久久久九九精品二区国产| 九色成人免费人妻av| 美女cb高潮喷水在线观看| 国产大屁股一区二区在线视频| 91在线观看av| 日韩欧美精品v在线| 国产伦人伦偷精品视频| 欧美丝袜亚洲另类 | 亚洲va日本ⅴa欧美va伊人久久| 丁香欧美五月| 国产精品女同一区二区软件 | 一卡2卡三卡四卡精品乱码亚洲| 国产熟女xx| 99国产极品粉嫩在线观看| 麻豆国产97在线/欧美| 青草久久国产| 尤物成人国产欧美一区二区三区| 国产精品1区2区在线观看.| 男女下面进入的视频免费午夜| 欧美zozozo另类| 蜜桃亚洲精品一区二区三区| 色5月婷婷丁香| 国产av在哪里看| 露出奶头的视频| xxxwww97欧美| 91狼人影院| 国产男靠女视频免费网站| 听说在线观看完整版免费高清| 能在线免费观看的黄片| 九色国产91popny在线| 日本成人三级电影网站| 亚洲美女视频黄频| 天天躁日日操中文字幕| 内地一区二区视频在线| 真人做人爱边吃奶动态| 国产乱人视频| 成年女人永久免费观看视频| 简卡轻食公司| 欧美一级a爱片免费观看看| 九九在线视频观看精品| 亚洲综合色惰| 午夜老司机福利剧场| 国内久久婷婷六月综合欲色啪| 激情在线观看视频在线高清| 热99在线观看视频| 亚洲18禁久久av| 亚洲欧美日韩卡通动漫| 国产成人aa在线观看| 又爽又黄a免费视频| 欧美日韩国产亚洲二区| 99久久久亚洲精品蜜臀av| 婷婷六月久久综合丁香| 欧美性猛交╳xxx乱大交人| 禁无遮挡网站| 日本一本二区三区精品| 亚洲成人免费电影在线观看| 草草在线视频免费看| 啪啪无遮挡十八禁网站| 久久精品国产99精品国产亚洲性色| 丁香六月欧美| 内射极品少妇av片p| 亚洲 国产 在线| 色吧在线观看| 色精品久久人妻99蜜桃| 免费看日本二区| 内射极品少妇av片p| 丁香欧美五月| 51午夜福利影视在线观看| 直男gayav资源| 亚洲第一区二区三区不卡| 久久久久久久久大av| 色精品久久人妻99蜜桃| 久久国产精品人妻蜜桃| 国产在线男女| 99国产极品粉嫩在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人与动物交配视频| 欧美最黄视频在线播放免费| 久久久久久九九精品二区国产| 久久久久久久午夜电影| 最近在线观看免费完整版| 免费人成在线观看视频色| av在线观看视频网站免费| 欧美黑人巨大hd| 精品人妻一区二区三区麻豆 | 日本免费a在线| 免费看日本二区| 国产精品98久久久久久宅男小说| 香蕉av资源在线| 99久久无色码亚洲精品果冻| 亚洲精品乱码久久久v下载方式| 99久久九九国产精品国产免费| 搡女人真爽免费视频火全软件 | 老女人水多毛片| 欧美日韩黄片免| 亚洲熟妇熟女久久| bbb黄色大片| 一本综合久久免费| 18禁黄网站禁片免费观看直播| av视频在线观看入口| 中文在线观看免费www的网站| 午夜福利在线观看免费完整高清在 | 一级a爱片免费观看的视频| 最后的刺客免费高清国语| 午夜精品一区二区三区免费看| 又粗又爽又猛毛片免费看| 听说在线观看完整版免费高清| 久久久久久国产a免费观看| 最后的刺客免费高清国语| 九色国产91popny在线| 午夜福利高清视频| 黄色视频,在线免费观看| 国产久久久一区二区三区| 亚洲色图av天堂| 高清日韩中文字幕在线| 国产亚洲精品久久久com| 日韩欧美在线乱码| 亚洲av成人不卡在线观看播放网| 美女xxoo啪啪120秒动态图 | 免费看a级黄色片| 久久亚洲真实| 免费观看精品视频网站| 婷婷精品国产亚洲av在线| 日韩av在线大香蕉| 丰满人妻熟妇乱又伦精品不卡| 国产免费一级a男人的天堂| 免费看美女性在线毛片视频| 蜜桃久久精品国产亚洲av| 麻豆久久精品国产亚洲av| 色哟哟·www| 亚洲一区高清亚洲精品| 国产精品一区二区三区四区免费观看 | 国产成人aa在线观看| 一个人看视频在线观看www免费| 久久精品久久久久久噜噜老黄 | 亚洲综合色惰| 精品熟女少妇八av免费久了| 免费在线观看日本一区| 国产男靠女视频免费网站| 一区二区三区四区激情视频 | 国产精品伦人一区二区| 看免费av毛片| 欧美+亚洲+日韩+国产| 欧美成人a在线观看| 特大巨黑吊av在线直播| bbb黄色大片| 色5月婷婷丁香| 色噜噜av男人的天堂激情| av女优亚洲男人天堂| 伦理电影大哥的女人| 日本精品一区二区三区蜜桃| 成年免费大片在线观看| 又黄又爽又免费观看的视频| 日本精品一区二区三区蜜桃| 一个人观看的视频www高清免费观看| 欧美zozozo另类| 又黄又爽又刺激的免费视频.| 日韩人妻高清精品专区| 亚洲欧美日韩卡通动漫| av在线老鸭窝| 午夜福利在线在线| 亚洲av第一区精品v没综合| 国产大屁股一区二区在线视频| 人人妻,人人澡人人爽秒播| 男女视频在线观看网站免费| 国产精品野战在线观看| 国产真实伦视频高清在线观看 | 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 一区二区三区高清视频在线| 一个人免费在线观看的高清视频| 欧美国产日韩亚洲一区| 国产一区二区三区视频了| 国产亚洲欧美在线一区二区| 极品教师在线视频| 免费在线观看影片大全网站| 日日干狠狠操夜夜爽| 天堂网av新在线| 久久精品国产自在天天线| 久久久国产成人免费| 成年版毛片免费区| 国产欧美日韩一区二区精品| 老女人水多毛片| av天堂在线播放| 亚洲精品久久国产高清桃花| 91在线观看av| 91麻豆精品激情在线观看国产| 中出人妻视频一区二区| 久久精品91蜜桃| 国产伦人伦偷精品视频| 怎么达到女性高潮| 男女下面进入的视频免费午夜| 欧美色视频一区免费| 欧美精品啪啪一区二区三区| 亚洲熟妇中文字幕五十中出| 欧美黄色淫秽网站| 亚洲精品在线美女| 一个人免费在线观看电影| 亚洲在线观看片| 国产午夜精品久久久久久一区二区三区 | 男人和女人高潮做爰伦理| 欧美激情国产日韩精品一区| 国产精品,欧美在线| 免费看日本二区| 1000部很黄的大片| 啪啪无遮挡十八禁网站| 9191精品国产免费久久| 亚洲激情在线av| 在线免费观看不下载黄p国产 | 日韩精品青青久久久久久| 男插女下体视频免费在线播放| 一个人看视频在线观看www免费| 床上黄色一级片| 黄色丝袜av网址大全| 婷婷精品国产亚洲av在线| 一区二区三区免费毛片| 在线十欧美十亚洲十日本专区| 少妇的逼好多水| 久久伊人香网站| 亚洲乱码一区二区免费版| 一级作爱视频免费观看| 欧美精品啪啪一区二区三区| 色综合站精品国产| 九色国产91popny在线| 国产一区二区在线观看日韩| 国产精华一区二区三区| 婷婷亚洲欧美| 极品教师在线免费播放| 在线播放国产精品三级| 色在线成人网| 51国产日韩欧美| 午夜两性在线视频| 一区二区三区四区激情视频 | 91午夜精品亚洲一区二区三区 | 日韩亚洲欧美综合| 少妇熟女aⅴ在线视频| 美女黄网站色视频| 校园春色视频在线观看| 久久精品人妻少妇| 精品免费久久久久久久清纯| 国产精品伦人一区二区| 999久久久精品免费观看国产| 久久6这里有精品| 波多野结衣巨乳人妻| 婷婷六月久久综合丁香| 国产精品亚洲av一区麻豆| 97碰自拍视频| 成人无遮挡网站| 国内精品一区二区在线观看| 成年女人永久免费观看视频| 一本精品99久久精品77| av国产免费在线观看| 无人区码免费观看不卡| 国产一区二区三区视频了| 一区二区三区免费毛片| 国产伦在线观看视频一区| 九九热线精品视视频播放| 国产又黄又爽又无遮挡在线| 最新中文字幕久久久久| 国产麻豆成人av免费视频| 亚洲专区国产一区二区| 男女做爰动态图高潮gif福利片| 成人无遮挡网站| 中文资源天堂在线| 日本精品一区二区三区蜜桃| 久久伊人香网站| 国产真实伦视频高清在线观看 | 国产精品美女特级片免费视频播放器| 麻豆成人av在线观看| 岛国在线免费视频观看| 国产亚洲精品综合一区在线观看| 亚洲内射少妇av| 久久伊人香网站| 亚洲自拍偷在线| 在线免费观看不下载黄p国产 | 波多野结衣高清无吗| 男女视频在线观看网站免费| 在线观看舔阴道视频| 深夜精品福利| 久久伊人香网站| eeuss影院久久| 国产v大片淫在线免费观看| 久久人人精品亚洲av| 午夜激情欧美在线| 精品人妻偷拍中文字幕| 最近最新免费中文字幕在线| 熟妇人妻久久中文字幕3abv| 亚洲av日韩精品久久久久久密| 啦啦啦观看免费观看视频高清| 国产精品免费一区二区三区在线| 亚洲精品乱码久久久v下载方式| 亚洲国产精品999在线| 别揉我奶头~嗯~啊~动态视频| 麻豆久久精品国产亚洲av| 午夜福利视频1000在线观看| 久久国产精品影院| 天天躁日日操中文字幕| 日韩精品中文字幕看吧| а√天堂www在线а√下载| 色5月婷婷丁香| 中文字幕熟女人妻在线| 脱女人内裤的视频| 国产成年人精品一区二区| 观看免费一级毛片| 久久精品国产亚洲av香蕉五月| 亚洲精品亚洲一区二区| 国产三级黄色录像| 脱女人内裤的视频| 天堂√8在线中文| 中文字幕久久专区| 久久久久久久久久黄片| 国产又黄又爽又无遮挡在线| 91麻豆精品激情在线观看国产| 国产av一区在线观看免费| 午夜福利18| 午夜激情欧美在线| 久99久视频精品免费| 久久精品91蜜桃| 最后的刺客免费高清国语| 国产精品乱码一区二三区的特点| 亚洲精品色激情综合| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本与韩国留学比较| 国产视频内射| 97碰自拍视频| 欧美xxxx性猛交bbbb| 午夜福利高清视频| 18禁黄网站禁片午夜丰满| 亚洲精品在线美女| 一本久久中文字幕| .国产精品久久| 天天躁日日操中文字幕| 在线观看66精品国产| 中亚洲国语对白在线视频| 国产毛片a区久久久久| 亚洲黑人精品在线| 伊人久久精品亚洲午夜| 欧美性猛交╳xxx乱大交人| 久久久久久久精品吃奶| 看免费av毛片| АⅤ资源中文在线天堂| 国产在线精品亚洲第一网站| 亚洲av.av天堂| 亚洲国产精品久久男人天堂| 精品国产三级普通话版| 特大巨黑吊av在线直播| 美女大奶头视频| 首页视频小说图片口味搜索| 国产野战对白在线观看| 久久精品国产99精品国产亚洲性色| 亚洲乱码一区二区免费版| 首页视频小说图片口味搜索| 国产野战对白在线观看| 国产亚洲欧美98| 欧美在线黄色| 校园春色视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 色综合欧美亚洲国产小说| 国产主播在线观看一区二区| www.熟女人妻精品国产| 欧美午夜高清在线| 99热这里只有是精品50| 亚洲熟妇熟女久久| 天堂av国产一区二区熟女人妻| 天堂√8在线中文| 在线国产一区二区在线| 国产精品自产拍在线观看55亚洲| 丁香六月欧美| or卡值多少钱| 久久久久久久久久成人| 少妇裸体淫交视频免费看高清| 搞女人的毛片| 欧美在线一区亚洲| 亚洲精品一卡2卡三卡4卡5卡| 女人十人毛片免费观看3o分钟| 搡老熟女国产l中国老女人| 51午夜福利影视在线观看| 免费av不卡在线播放| 青草久久国产| 国产高潮美女av| 在线免费观看不下载黄p国产 | 日韩免费av在线播放| 免费在线观看亚洲国产| 亚洲乱码一区二区免费版| 男人狂女人下面高潮的视频| avwww免费| 在线播放国产精品三级| 欧美黄色片欧美黄色片| 国产精品电影一区二区三区| 亚洲av电影在线进入| 免费av观看视频| 90打野战视频偷拍视频| 国产午夜精品论理片| 两个人视频免费观看高清| 看十八女毛片水多多多| 99riav亚洲国产免费| 精品国产亚洲在线| 无人区码免费观看不卡| 97碰自拍视频| 又爽又黄无遮挡网站| eeuss影院久久| 亚洲国产精品久久男人天堂| 男女床上黄色一级片免费看| 精品一区二区三区av网在线观看| 国产国拍精品亚洲av在线观看| 精品国内亚洲2022精品成人| 亚洲无线在线观看| 久久久久久久久久成人| 女人十人毛片免费观看3o分钟| 全区人妻精品视频| 亚洲av免费高清在线观看| 91午夜精品亚洲一区二区三区 | 麻豆av噜噜一区二区三区| 精品久久久久久久久亚洲 | 别揉我奶头 嗯啊视频| 国产蜜桃级精品一区二区三区| 中文在线观看免费www的网站| 久久人人爽人人爽人人片va | 国内少妇人妻偷人精品xxx网站| 亚洲欧美激情综合另类| 精华霜和精华液先用哪个| 校园春色视频在线观看| 欧美黄色片欧美黄色片| 日韩大尺度精品在线看网址| or卡值多少钱| 激情在线观看视频在线高清| 一本久久中文字幕| 成人国产一区最新在线观看|