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

    考慮熱力學(xué)效應(yīng)的空化模型修正及其適用性研究

    2019-03-06 09:16:54張德勝萬(wàn)福來(lái)王超超施衛(wèi)東
    關(guān)鍵詞:實(shí)驗(yàn)模型

    張德勝 萬(wàn)福來(lái) 許 彬 王超超 施衛(wèi)東

    (1.江蘇大學(xué)流體機(jī)械工程技術(shù)研究中心, 鎮(zhèn)江 212013; 2.南通大學(xué)機(jī)械工程學(xué)院, 南通 226019)

    0 引言

    空化是一種包含湍流、多相流、相變、噪聲、可壓縮和非定常特性的復(fù)雜流動(dòng)現(xiàn)象[1-3]。近年來(lái),低溫介質(zhì)的運(yùn)用促進(jìn)了低溫空化的研究。液化天然氣等低溫介質(zhì)已經(jīng)成為我國(guó)能源領(lǐng)域迅猛發(fā)展的新型產(chǎn)業(yè)[4-5],而其溫度為-160℃左右,極容易導(dǎo)致泵內(nèi)發(fā)生空化現(xiàn)象,這對(duì)低溫介質(zhì)的空化性能提出了較高的研究要求。由于大部分空化都發(fā)生在常溫水中,所以通常假設(shè)空化過程為等溫過程,而低溫流體介質(zhì)物理屬性對(duì)溫度變化敏感,故熱力學(xué)效應(yīng)決定低溫流體空化流數(shù)值模擬的準(zhǔn)確性顯著[6-7]。

    在低溫流體空化的數(shù)值計(jì)算和實(shí)驗(yàn)研究方面,STEPANOFF等[8]最早提出了一種B因子理論,用于預(yù)測(cè)泵中熱力學(xué)效應(yīng)對(duì)揚(yáng)程的影響;HORD[9]以液氮和液氫為介質(zhì),做了較為全面系統(tǒng)的低溫流體的空化實(shí)驗(yàn),進(jìn)一步加深了對(duì)低溫空化特性的認(rèn)識(shí)。CERVONE等[10]對(duì)不同溫度水為介質(zhì)的NACA0015翼形水翼空化開展了一系列實(shí)驗(yàn)研究,分析了熱力學(xué)效應(yīng)對(duì)其空化特性的影響。王巍等[11]對(duì)不同溫度水為介質(zhì)的NACA0066(MOD)翼型開展了翼型表面空化場(chǎng)的流動(dòng)分析。由于低溫流體工作環(huán)境的限制,極大增加了實(shí)驗(yàn)的難度。近年來(lái),數(shù)值計(jì)算成為研究低溫流體空化的重要手段,HOSANGADI等[12]采用Merkle空化模型對(duì)液氫和液氮繞翼型空化流動(dòng)進(jìn)行了模擬計(jì)算,但模擬結(jié)果與實(shí)驗(yàn)對(duì)比誤差較大。為了提高數(shù)值模擬預(yù)測(cè)精度,TANI等[13]在B因子理論基礎(chǔ)上提出了一種新的空化模型,通過與實(shí)驗(yàn)數(shù)據(jù)相比較確定了空化模型的蒸發(fā)凝結(jié)系數(shù),使得計(jì)算結(jié)果較為吻合。國(guó)內(nèi)學(xué)者也對(duì)低溫流體空化做了一些研究[14-16],評(píng)價(jià)了不同空化模型對(duì)低溫空化的數(shù)值計(jì)算情況,并改進(jìn)了個(gè)別空化模型,使得數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)更為接近,但這些研究都是二維翼型的數(shù)值模擬計(jì)算。因此,本文重點(diǎn)開展三維翼型考慮熱力學(xué)效應(yīng)的空化模型改進(jìn)和驗(yàn)證研究。

    本文通過CFX前處理編輯,將3種典型空化模型嵌入CFX軟件中,并將水、液氮熱傳導(dǎo)系數(shù)以及飽和蒸汽壓強(qiáng)等隨溫度變化的物性參數(shù)引入到求解代碼程序。 同時(shí)考慮空化過程汽化潛熱的影響,并在B因子理論的基礎(chǔ)上,提出一種空化模型的修正方法。采用原始的和修正的空化模型,計(jì)算水及液氮繞翼型的空化流動(dòng),并與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比,評(píng)估不同空化模型對(duì)空化流動(dòng)的熱力學(xué)效應(yīng)的敏感性和修正的空化模型的適用性,為開展低溫空化研究提供理論支撐。

    1 控制方程與數(shù)值方法

    1.1 基本方程

    對(duì)于考慮熱力學(xué)效應(yīng)的空化流動(dòng)問題數(shù)值計(jì)算的控制方程,除了連續(xù)性方程和動(dòng)量方程外,還包括包含能量源項(xiàng)的能量方程,依次為

    (1)

    (2)

    (3)

    式中ρm——混合相密度p——壓強(qiáng)

    μ——混合介質(zhì)的動(dòng)力黏性系數(shù)

    keff——熱傳導(dǎo)系數(shù)

    L——汽化潛熱

    Cp——定壓比熱容

    fv——?dú)庀嗟馁|(zhì)量分?jǐn)?shù)

    T——溫度t——時(shí)間

    xi——x在i方向分量

    xj——x在j方向分量

    δij——克羅內(nèi)克函數(shù)

    ui、uj——速度在i、j方向分量

    μt——湍流粘度

    1.2 湍流模型

    計(jì)算采用了由MENTER[17]提出的SSTk-ω兩方程湍流模型,公式為

    (4)

    (5)

    渦黏系數(shù)計(jì)算公式為

    (6)

    式中ρ——密度k——湍動(dòng)能

    ω——湍流頻率

    σk3、σω2、σω3、α3、β3、β′——湍流模型常數(shù)

    F1、F2——混合函數(shù)

    S——剪切力張量的常數(shù)項(xiàng)

    Pkb——浮力湍動(dòng)能生成項(xiàng)

    Pk——粘性力湍動(dòng)能生成項(xiàng)

    Pωb——湍流模型自定義項(xiàng)

    a1——湍流模型常數(shù)

    湍流模型常數(shù)取值參照文獻(xiàn)[18]。

    1.3 空化模型

    (7)

    式中u——速度矢量

    αv——?dú)庀囿w積分?jǐn)?shù)

    ρv——?dú)庀嗝芏?/p>

    1.3.1Zwart空化模型

    Zwart空化模型[19]是基于單個(gè)空泡的生成和發(fā)展時(shí)空泡體積的變化,基于Rayleigh-Plesse空泡生長(zhǎng)方程推導(dǎo)出的蒸發(fā)源項(xiàng)和凝結(jié)源項(xiàng)表達(dá)式

    (8)

    (9)

    式中RB——平均空泡半徑,取1×10-6m

    ρl——液相密度pv——?dú)庀鄩毫?/p>

    αnuc——?dú)夂梭w積分?jǐn)?shù),取5×10-4

    Fvap、Fcond——蒸發(fā)、凝結(jié)系數(shù),取50、0.01

    1.3.2Merkle空化模型

    Merkle空化模型[20]是基于汽液兩相流模型,由混合密度導(dǎo)出相間質(zhì)量傳輸率,公式為

    (10)

    (11)

    式中U∞——參考速度

    t∞——參考時(shí)間

    經(jīng)驗(yàn)系數(shù)Cvap、Ccond分別取1、80。

    1.3.3Singhal空化模型

    Singhal空化模型[21]是基于汽液兩相流模型,綜合考慮了空泡在相變過程中表面張力和非凝結(jié)性氣體的影響,公式為

    (12)

    (13)

    式中fg——非凝結(jié)氣體的質(zhì)量分?jǐn)?shù),取10-8

    σ——液體表面張力系數(shù)

    經(jīng)驗(yàn)系數(shù)Cvap、Ccond分別取0.02、0.01。

    1.4 空化模型的修正方法

    在上述空化模型中,通常認(rèn)為質(zhì)量傳輸只由液相和氣相的壓力差來(lái)驅(qū)動(dòng),但在高溫水和低溫液體中,發(fā)生空化時(shí),液體汽化吸收汽化潛熱,導(dǎo)致空泡附近液體溫度降低,從而形成空泡外的薄液體邊界層,假設(shè)泡內(nèi)的氣體溫度是均勻的,其壓力等于飽和壓力,溫度等于飽和壓力對(duì)應(yīng)下的飽和溫度[13],由于熱力學(xué)邊界層的存在,泡內(nèi)和泡外形成了一溫度差ΔT,則相變的熱平衡公式為

    ρvυvL=ρlυlCplΔT

    (14)

    式中υv、υl——?dú)庀?、液相的體積質(zhì)量流量

    Cpl——液相的定壓比熱容

    B因子方法[22]及理想氣液混合時(shí)B因子的表達(dá)形式為[23]

    (15)

    式中B——無(wú)量綱因子

    國(guó)內(nèi)外化工行業(yè)常用的純液體飽和蒸汽壓的三參數(shù)方程Antoine方程[24]為

    lgpsat=a-b/(T+c)

    (16)

    式中a、b、c——Antoine常數(shù),取值參照文獻(xiàn)[23]

    psat——液相飽和蒸汽壓力

    假定空泡周圍熱力學(xué)狀態(tài)保持為飽和狀態(tài),則可得由于潛熱傳遞而形成的局部飽和蒸汽壓的變化ΔpL,其公式為

    (17)

    除此之外,研究表明湍動(dòng)能對(duì)空化產(chǎn)生重要的影響[24],采用文獻(xiàn)[24]中所提出的方法來(lái)計(jì)算湍動(dòng)能k對(duì)當(dāng)?shù)仄瘔簭?qiáng)的影響,其公式為

    Δpturb=0.195ρmk

    (18)

    式中 Δpturb——湍動(dòng)能引起的壓力變化量

    基于以上分析,可對(duì)1.3節(jié)中3個(gè)空化模型中蒸發(fā)源項(xiàng)和凝結(jié)源項(xiàng)就考慮熱力學(xué)效應(yīng)的影響進(jìn)行修正,從而可得

    pv=psat+ΔpL+Δpturb

    (19)

    2 不同溫度水的空化流數(shù)值模擬

    2.1 計(jì)算區(qū)域及邊界條件

    為了評(píng)估各空化模型的熱力學(xué)修正效果,首先對(duì)以不同溫度水為介質(zhì)的NACA0015翼形的空化流場(chǎng)進(jìn)行了模擬計(jì)算,并與CERVONE等[10]的實(shí)驗(yàn)值進(jìn)行比較。翼型尺寸與計(jì)算域選取與實(shí)驗(yàn)條件保持一致,如圖1所示。設(shè)置邊界條件:進(jìn)口給定速度與溫度,出口為平均靜壓;壁面為無(wú)滑移絕緣壁面條件。為更準(zhǔn)確地計(jì)算空化流場(chǎng),采用結(jié)構(gòu)網(wǎng)格,并在翼型近壁面區(qū)域進(jìn)行網(wǎng)格加密,使分子到最近壁面的無(wú)量綱距離y+為10~50之間,大多數(shù)學(xué)者認(rèn)為y+值不超過60[25],則計(jì)算y+值滿足壁面函數(shù)要求,如圖2所示。為了對(duì)計(jì)算所用網(wǎng)格進(jìn)行驗(yàn)證,本文選取了3套網(wǎng)格,分別模擬計(jì)算了無(wú)空化流場(chǎng),如圖3a所示。所得壓力系數(shù)pc與實(shí)驗(yàn)值[10]進(jìn)行對(duì)比驗(yàn)證,其中Lc為翼型長(zhǎng)度百分?jǐn)?shù)??芍?jì)算所得壓力系數(shù)與實(shí)驗(yàn)較為接近,如圖3b所示,并考慮計(jì)算經(jīng)濟(jì)性,選取方案2進(jìn)行后續(xù)的計(jì)算。不同方案的網(wǎng)格信息如表1所示,計(jì)算工況見表2。

    圖1 三維水體模型及邊界條件Fig.1 Schematic of 3D model with water and boundary conditions of cases

    圖2 NACA0015翼型網(wǎng)格分布Fig.2 Grids distribution near NACA0015 hydrofoil

    圖3 NACA0015網(wǎng)格無(wú)關(guān)性驗(yàn)證Fig.3 Verification of NACA0015 mesh independence

    計(jì)算中采用的無(wú)量綱參數(shù)為無(wú)窮遠(yuǎn)處空化數(shù)σ∞、當(dāng)?shù)乜栈瘮?shù)σc、壓力系數(shù)pc、當(dāng)?shù)乜栈瘮?shù)與無(wú)窮遠(yuǎn)處空化數(shù)之差σ,分別定義為

    (20)

    (21)

    (22)

    σ=σc-σ∞

    (23)

    式中pout——無(wú)窮遠(yuǎn)處壓力

    T∞——無(wú)窮遠(yuǎn)處溫度

    uin——進(jìn)口速度Tc——當(dāng)?shù)販囟?/p>

    pv(Tc) ——溫度為Tc時(shí)的飽和蒸氣壓力

    pv(T∞) ——溫度為T∞時(shí)的飽和蒸氣壓力

    表1 網(wǎng)格信息Tab.1 Information of different grids

    表2 不同溫度水空化數(shù)值模擬邊界條件Tab.2 Boundary conditions of numerical simulation about cavitation in water at different temperatures

    2.2 計(jì)算結(jié)果與討論

    為了分析修正空化模型在不同溫度水中空化數(shù)值模擬中的應(yīng)用情況,分別運(yùn)用1.4節(jié)中考慮熱力學(xué)效應(yīng)修正的Zwart、Merkle及Singhal 3種空化模型,并對(duì)298、323、343 K溫度水繞NACA0015翼型空化進(jìn)行了數(shù)值模擬研究,并與實(shí)驗(yàn)[8]所得翼型吸力面的壓力系數(shù)分布進(jìn)行對(duì)比分析,如圖4所示。結(jié)果表明,在3種溫度下,修正前后Singhal模型及Zwart模型模擬所得翼型吸力面壓力系數(shù)分布無(wú)明顯差別。溫度為298 K和323 K時(shí),水的熱力學(xué)效應(yīng)比較穩(wěn)定,壓力隨溫度變化梯度較小,3種空化模型修正前后的模擬結(jié)果差異較小。而修正的Merkle模型在水溫為343 K所模擬的結(jié)果表現(xiàn)出較強(qiáng)的修正效果,修正前后有較明顯差別。修正Merkle模型數(shù)值模擬得到的翼型吸力面低壓區(qū)范圍較之未修正前的模擬結(jié)果均有所增加,對(duì)應(yīng)壓力有所降低,并隨溫度的增加,降低得越多,體現(xiàn)出不同溫度水中空化吸收潛熱,引起飽和蒸汽壓強(qiáng)不同程度降低的熱力學(xué)效應(yīng)作用,且越靠近臨界溫度,其敏感度越高。而同時(shí)對(duì)應(yīng)空穴閉合區(qū)壓力變化梯度均明顯增大,且隨溫度的增加,其壓力變化梯度也隨之增大,從而更加接近于實(shí)驗(yàn)值。

    圖4 不同溫度下不同修正空化模型模擬得到的NACA0015翼型吸力面壓力系數(shù)分布Fig.4 Pressure coefficient distribution on suction surface of NACA0015 hydrofoil

    3 不同溫度液氮的空化流數(shù)值模擬

    3.1 計(jì)算區(qū)域及邊界條件

    HORD[9]針對(duì)不同溫度液氮繞水翼空化進(jìn)行了較為全面系統(tǒng)的低溫流體空化實(shí)驗(yàn)研究,其實(shí)驗(yàn)段及水翼結(jié)構(gòu)如圖5和圖6所示。

    通過三維造型生成與實(shí)驗(yàn)所對(duì)應(yīng)三維水翼,并根據(jù)實(shí)驗(yàn)中實(shí)驗(yàn)段的尺寸進(jìn)行了實(shí)驗(yàn)流道的建模。

    圖5 實(shí)驗(yàn)段結(jié)構(gòu)Fig.5 Test section structure

    圖6 HORD水翼結(jié)構(gòu)圖Fig.6 Diagram of hydrofoil structure

    同時(shí)由于其對(duì)稱性結(jié)構(gòu),在數(shù)值模擬過程中可進(jìn)行對(duì)稱邊界條件的設(shè)置,故只構(gòu)建了一半的對(duì)稱水體模型。三維對(duì)稱水體模型結(jié)構(gòu)如圖7所示??紤]與實(shí)驗(yàn)的一致性,進(jìn)口設(shè)置為壓力進(jìn)口,其靜壓與實(shí)驗(yàn)對(duì)應(yīng)進(jìn)口壓力相一致,來(lái)流溫度根據(jù)實(shí)驗(yàn)進(jìn)口所測(cè)的溫度來(lái)設(shè)置;出口設(shè)置為速度出口,其速度與實(shí)驗(yàn)保持一致,邊界條件如表3所示。為精準(zhǔn)地模擬計(jì)算空化流場(chǎng),在ICEM 軟件中進(jìn)行六面體結(jié)構(gòu)網(wǎng)格的劃分,考慮到水翼前緣的圓頭形狀的影響,采用C型結(jié)構(gòu)化網(wǎng)格,在水翼近壁面區(qū)域進(jìn)行網(wǎng)格加密如圖8所示。選取3套網(wǎng)格計(jì)算液氮溫度83.06 K時(shí)翼型表面壓力系數(shù)分布情況,如圖9a所示,各方案網(wǎng)格信息如表4所示。模擬所得壓力系數(shù)與實(shí)驗(yàn)值[9]進(jìn)行比較,如圖9b所示??紤]數(shù)值計(jì)算的準(zhǔn)確性與經(jīng)濟(jì)性,選擇方案2作為后續(xù)模擬計(jì)算網(wǎng)格。

    圖7 三維對(duì)稱水體模型及邊界條件Fig.7 Schematic of symmetrical 3D model with water and boundary conditions of cases

    3.2 計(jì)算結(jié)果與討論

    為了分析不同空化模型在低溫流體空化數(shù)值模擬中的應(yīng)用情況,分別運(yùn)用Zwart、Merkle及Singhal空化模型,采用等溫isothermal的流場(chǎng)熱傳輸模式,分別對(duì)83.06、77.64 K兩種不同溫度液氮繞水翼的空化進(jìn)行數(shù)值計(jì)算模擬研究,并與HORD的實(shí)驗(yàn)數(shù)據(jù)[9]進(jìn)行對(duì)比分析,所得水翼表面的壓力分布如圖10所示。在兩種液氮中,3種模型模擬所得低壓區(qū)壓強(qiáng)保持為對(duì)應(yīng)遠(yuǎn)場(chǎng)溫度下的飽和蒸汽壓強(qiáng)不變,且均大于實(shí)驗(yàn)值,說明液氮中熱力學(xué)效應(yīng)明顯,不容忽視,其中空化與常溫水中不同,不能假設(shè)為等溫過程。而模擬所得低壓區(qū)范圍差異較大,Merkle模型模擬結(jié)果較為接近實(shí)驗(yàn)值,而Zwart模型與Singhal模型兩者模擬結(jié)果相近均偏離實(shí)驗(yàn)值較多,過多地預(yù)測(cè)了低壓區(qū)即空化區(qū)域的范圍。

    表3 不同液氮中空化數(shù)值模擬邊界條件Tab.3 Boundary conditions of numerical simulation about cavitation in nitrogen at different temperatures

    圖8 HORD翼型網(wǎng)格分布Fig.8 Grids distribution near HORD hydrofoil

    圖9 HORD翼型網(wǎng)格變化Fig.9 Changes of HORD hydrofoil mesh

    方案網(wǎng)格節(jié)點(diǎn)數(shù)最小角度最大寬高比最小質(zhì)量1325409653.462920.8042388265653.462940.8043486105653.462920.804

    圖10 空化模型模擬所得HORD水翼表面的壓力分布Fig.10 Pressure distribution on surface of HORD hydrofoil with different cavitation models

    綜合上述分析可知,Merkle模型在低溫流體空化數(shù)值模擬中有較好的可行性,其與實(shí)驗(yàn)吻合度較為接近。為了分析修正空化模型在不同溫度液氮空化數(shù)值模擬中的應(yīng)用情況,運(yùn)用考慮熱力學(xué)效應(yīng)修正的Merkle模型,分別對(duì)83.06、77.64 K兩種不同溫度液氮繞水翼的空化進(jìn)行數(shù)值計(jì)算模擬,并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比分析,所得水翼表面的壓力分布如圖11所示。結(jié)果表明,修正前后Merkle模型所模擬的結(jié)果有較明顯差別,修正的Merkle模型模擬所得水翼表面低壓區(qū)范圍較之未修正前的模擬結(jié)果均有所減小,低壓區(qū)壓力有所降低,且83.06 K中降低量較77.64 K下更多,使得兩者模擬結(jié)果更加吻合實(shí)驗(yàn)低壓區(qū)壓力,這體現(xiàn)出低溫流體中空化吸收潛熱,引起飽和蒸汽壓強(qiáng)不同程度降低的熱力學(xué)效應(yīng)作用,也進(jìn)一步說明了越靠近其臨界溫度,其敏感度越高。

    圖11 修正空化模型數(shù)值模擬所得HORD水翼表面的壓力分布Fig.11 Pressure distribution on surface of HORD hydrofoil with modified cavitation model

    采用考慮熱力學(xué)效應(yīng)修正及等溫模式的Merkle模型模擬所得蒸汽體積分?jǐn)?shù)分布情況,如圖12所示,圖中y軸表示翼型絕對(duì)尺寸,x軸表示翼型相對(duì)尺寸。結(jié)果表明,等溫模式下,無(wú)熱力學(xué)效應(yīng)作用,空穴區(qū)域內(nèi)的壓強(qiáng)恒定為對(duì)應(yīng)遠(yuǎn)場(chǎng)溫度下的飽和蒸汽壓強(qiáng),并且83.06 K中對(duì)應(yīng)遠(yuǎn)場(chǎng)空化數(shù)為1.706,而77.64 K中為1.765,故在83.06 K中的空化區(qū)域大于77.64 K的,其空穴長(zhǎng)度明顯較長(zhǎng)。而考慮熱力學(xué)效應(yīng)時(shí),局部的溫降,導(dǎo)致飽和蒸汽壓強(qiáng)的降低,從而降低了空化強(qiáng)度,使得空化區(qū)域的蒸汽體積分?jǐn)?shù)減小,空穴長(zhǎng)度明顯縮短。

    根據(jù)當(dāng)?shù)乜栈瘮?shù)σc的定義,為了更好地表征熱力學(xué)效應(yīng)在低溫流體空化中的影響,計(jì)算了兩種液氮中,修正Merkle模型模擬所得空化流場(chǎng)的當(dāng)?shù)乜栈瘮?shù)與無(wú)窮遠(yuǎn)處空化數(shù)之差σ分布情況,如圖13所示。考慮熱力學(xué)效應(yīng)時(shí),當(dāng)?shù)乜栈瘮?shù)最大的區(qū)域與最大壓降區(qū)、最大溫降區(qū)相對(duì)應(yīng),都在空化核心區(qū)域,其當(dāng)?shù)乜栈瘮?shù)均大于相應(yīng)的遠(yuǎn)場(chǎng)空化數(shù),由此與等溫模式下相比,考慮熱力學(xué)效應(yīng)時(shí)的空化強(qiáng)度降低。

    圖12 不同液氮中模擬所得HORD水翼表面蒸汽體積分?jǐn)?shù)云圖Fig.12 Vapour volume fraction on surface of HORD hydrofoil simulated in different liquid nitrogen

    圖13 不同液氮中模擬所得HOAD水翼表面σ分布云圖Fig.13 Distribution nephogram of σ on surface of HORD hydrofoil simulated in different liquid nitrogen

    4 結(jié)論

    (1)采用Zwart、Merkle及Singhal共3種空化模型及對(duì)應(yīng)修正后的空化模型,分別對(duì)3種溫度水進(jìn)行數(shù)值模擬,并與實(shí)驗(yàn)進(jìn)行對(duì)比,結(jié)果表明修正的Merkle模型能更好地體現(xiàn)熱力學(xué)效應(yīng)對(duì)空化的影響。

    (2)采用3種空化模型分別數(shù)值模擬83.06、77.64 K兩種液氮繞水翼空化流場(chǎng),并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比,發(fā)現(xiàn)Merkle模型數(shù)值模擬結(jié)果與實(shí)驗(yàn)值較為接近,適用性相對(duì)較強(qiáng)??紤]熱力學(xué)效應(yīng)修正的Merkle空化模型,對(duì)兩種不同溫度液氮中空化進(jìn)行了數(shù)值模擬研究,并與實(shí)驗(yàn)所得水翼表面壓力分布數(shù)據(jù)進(jìn)行對(duì)比分析,再次驗(yàn)證了修正空化模型的適用性及空化模型修正方法在低溫流體空化數(shù)值模擬中的可行性。

    (3)在修正Merkle模型數(shù)值模擬結(jié)果的基礎(chǔ)上,與等溫模式的數(shù)值模擬結(jié)果相對(duì)比,其空化區(qū)域由于溫降的產(chǎn)生,使得對(duì)應(yīng)飽和蒸汽壓強(qiáng)有所降低,致使空化強(qiáng)度減弱,空穴長(zhǎng)度縮短,當(dāng)?shù)乜栈瘮?shù)均大于遠(yuǎn)場(chǎng)空化數(shù)。且相比77.64 K,83.06 K的溫度較高,對(duì)于熱力學(xué)效應(yīng)的敏感性更高,其空化區(qū)域較多的壓降使其熱力學(xué)效應(yīng)更顯著,空化受到抑制作用更明顯,空穴縮短更多,空化區(qū)域蒸汽體積分?jǐn)?shù)降低更多。

    猜你喜歡
    實(shí)驗(yàn)模型
    一半模型
    記一次有趣的實(shí)驗(yàn)
    微型實(shí)驗(yàn)里看“燃燒”
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    做個(gè)怪怪長(zhǎng)實(shí)驗(yàn)
    3D打印中的模型分割與打包
    NO與NO2相互轉(zhuǎn)化實(shí)驗(yàn)的改進(jìn)
    實(shí)踐十號(hào)上的19項(xiàng)實(shí)驗(yàn)
    太空探索(2016年5期)2016-07-12 15:17:55
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    超碰av人人做人人爽久久| 久热久热在线精品观看| 久久久精品大字幕| 精品国产露脸久久av麻豆 | 成人国产麻豆网| 免费看美女性在线毛片视频| 亚洲一区高清亚洲精品| 欧美变态另类bdsm刘玥| 边亲边吃奶的免费视频| 久久久久久久久大av| 国产高清有码在线观看视频| 成人二区视频| 欧美激情国产日韩精品一区| 成人三级黄色视频| 嫩草影院精品99| 亚洲国产精品成人久久小说| 免费人成在线观看视频色| av.在线天堂| 亚洲av成人精品一二三区| 校园人妻丝袜中文字幕| 亚洲欧美精品综合久久99| 在线播放无遮挡| 在线观看av片永久免费下载| 淫秽高清视频在线观看| 建设人人有责人人尽责人人享有的 | 99热这里只有精品一区| 亚洲经典国产精华液单| 在线观看一区二区三区| 精品久久久久久久末码| 日本三级黄在线观看| 岛国毛片在线播放| 久久精品夜夜夜夜夜久久蜜豆| 免费大片18禁| 亚洲国产成人一精品久久久| 国产精品1区2区在线观看.| eeuss影院久久| 视频中文字幕在线观看| 高清毛片免费看| 日韩,欧美,国产一区二区三区 | 国产精品永久免费网站| 全区人妻精品视频| 亚洲欧洲日产国产| 国产乱来视频区| 波多野结衣巨乳人妻| 国产女主播在线喷水免费视频网站 | 久久精品影院6| 男人的好看免费观看在线视频| 亚洲国产精品国产精品| 成年女人看的毛片在线观看| 亚洲内射少妇av| 国产精品熟女久久久久浪| 国产综合懂色| 日本免费在线观看一区| 97超视频在线观看视频| 国产片特级美女逼逼视频| 国产精品久久久久久久久免| 老师上课跳d突然被开到最大视频| 免费播放大片免费观看视频在线观看 | 国产精品精品国产色婷婷| 精品人妻视频免费看| 国产成人免费观看mmmm| 久久久久久大精品| 久久综合国产亚洲精品| 一夜夜www| 乱人视频在线观看| 床上黄色一级片| 亚洲精品色激情综合| 国产精品三级大全| 少妇的逼好多水| 免费搜索国产男女视频| 非洲黑人性xxxx精品又粗又长| 久久精品国产亚洲av天美| 国产精品99久久久久久久久| 成人毛片60女人毛片免费| 神马国产精品三级电影在线观看| 日韩中字成人| 久久亚洲精品不卡| 亚洲精品色激情综合| 日日干狠狠操夜夜爽| 韩国av在线不卡| 午夜爱爱视频在线播放| 日本黄色片子视频| 美女被艹到高潮喷水动态| 91精品一卡2卡3卡4卡| 亚洲三级黄色毛片| 女人久久www免费人成看片 | 亚洲精品乱久久久久久| 亚洲av免费在线观看| 中文精品一卡2卡3卡4更新| 69av精品久久久久久| 国产中年淑女户外野战色| 丝袜喷水一区| 波多野结衣巨乳人妻| 国产av码专区亚洲av| 91aial.com中文字幕在线观看| 午夜a级毛片| 国产精品精品国产色婷婷| 老司机影院成人| 久久草成人影院| 国产视频内射| 国产久久久一区二区三区| 男女那种视频在线观看| 国产午夜精品一二区理论片| 99视频精品全部免费 在线| 精品国产三级普通话版| 一夜夜www| 日韩 亚洲 欧美在线| 男女那种视频在线观看| 亚洲国产高清在线一区二区三| 黄片wwwwww| 我要看日韩黄色一级片| 精品无人区乱码1区二区| 久热久热在线精品观看| 黄色欧美视频在线观看| 亚洲电影在线观看av| 国产精品一及| 国产一区二区亚洲精品在线观看| 搡老妇女老女人老熟妇| 精品人妻熟女av久视频| 成年av动漫网址| 久久久久久久久中文| 身体一侧抽搐| 久久欧美精品欧美久久欧美| 免费播放大片免费观看视频在线观看 | 村上凉子中文字幕在线| 国产成年人精品一区二区| 精品免费久久久久久久清纯| 岛国毛片在线播放| 欧美最新免费一区二区三区| 纵有疾风起免费观看全集完整版 | 欧美激情国产日韩精品一区| 国产精品无大码| 97在线视频观看| av视频在线观看入口| 建设人人有责人人尽责人人享有的 | av卡一久久| 美女cb高潮喷水在线观看| 亚洲欧美精品综合久久99| 亚洲国产精品合色在线| 1000部很黄的大片| 村上凉子中文字幕在线| 国产精品女同一区二区软件| 亚洲最大成人手机在线| 国产精品蜜桃在线观看| 日韩在线高清观看一区二区三区| 天天一区二区日本电影三级| 亚洲中文字幕一区二区三区有码在线看| 日本五十路高清| 丰满人妻一区二区三区视频av| 久久99蜜桃精品久久| 亚洲性久久影院| 免费观看在线日韩| av免费观看日本| 九九热线精品视视频播放| 两个人视频免费观看高清| 久久午夜福利片| 我要看日韩黄色一级片| 简卡轻食公司| 欧美一级a爱片免费观看看| h日本视频在线播放| 亚洲精品,欧美精品| 九色成人免费人妻av| 亚洲va在线va天堂va国产| 亚洲av不卡在线观看| 久久亚洲国产成人精品v| 纵有疾风起免费观看全集完整版 | 日本熟妇午夜| 日本免费一区二区三区高清不卡| 少妇猛男粗大的猛烈进出视频 | 亚洲av不卡在线观看| 久久精品91蜜桃| 少妇的逼好多水| 国产单亲对白刺激| 纵有疾风起免费观看全集完整版 | 一本久久精品| 伊人久久精品亚洲午夜| 日韩三级伦理在线观看| 青春草国产在线视频| 少妇高潮的动态图| 久久久久久国产a免费观看| 欧美极品一区二区三区四区| ponron亚洲| 美女xxoo啪啪120秒动态图| 一卡2卡三卡四卡精品乱码亚洲| 在线天堂最新版资源| 中国国产av一级| 国产在线一区二区三区精 | 丰满少妇做爰视频| 熟妇人妻久久中文字幕3abv| 国产高清视频在线观看网站| 亚洲av不卡在线观看| 人妻夜夜爽99麻豆av| 18禁在线播放成人免费| 久久精品熟女亚洲av麻豆精品 | 欧美区成人在线视频| 国产精品无大码| 亚洲高清免费不卡视频| 中文资源天堂在线| 高清av免费在线| 欧美成人免费av一区二区三区| 国产成人精品久久久久久| av在线天堂中文字幕| 亚洲天堂国产精品一区在线| 精品人妻视频免费看| 一级毛片我不卡| 青春草视频在线免费观看| 日本色播在线视频| 亚洲成人久久爱视频| 看十八女毛片水多多多| or卡值多少钱| 青春草视频在线免费观看| av免费观看日本| 精品一区二区三区视频在线| 国产老妇女一区| 亚洲精品国产成人久久av| 国产av一区在线观看免费| 亚洲av中文av极速乱| 欧美成人一区二区免费高清观看| 欧美日韩国产亚洲二区| 变态另类丝袜制服| 又爽又黄无遮挡网站| 亚洲成人av在线免费| 舔av片在线| 高清在线视频一区二区三区 | 国产精品嫩草影院av在线观看| 麻豆成人av视频| 亚洲一级一片aⅴ在线观看| 国产乱来视频区| 国产精品久久久久久久电影| 国产精品精品国产色婷婷| 又黄又爽又刺激的免费视频.| 亚洲精品,欧美精品| 国产亚洲av嫩草精品影院| 麻豆成人av视频| 欧美精品国产亚洲| 亚洲图色成人| 波多野结衣高清无吗| 国产高清有码在线观看视频| a级一级毛片免费在线观看| 九九久久精品国产亚洲av麻豆| 日本免费a在线| 亚洲成人久久爱视频| 午夜a级毛片| 五月玫瑰六月丁香| 国产精华一区二区三区| 欧美97在线视频| 在线观看av片永久免费下载| 男女国产视频网站| 超碰97精品在线观看| 搞女人的毛片| 3wmmmm亚洲av在线观看| 日韩av在线大香蕉| 直男gayav资源| 在线天堂最新版资源| 又粗又硬又长又爽又黄的视频| 欧美色视频一区免费| av专区在线播放| 日韩精品青青久久久久久| 免费观看人在逋| 久久久国产成人精品二区| 久久久久久久久中文| 精品熟女少妇av免费看| 黄片wwwwww| or卡值多少钱| 成人性生交大片免费视频hd| 最近最新中文字幕大全电影3| 黄色配什么色好看| 天天躁夜夜躁狠狠久久av| 国产亚洲午夜精品一区二区久久 | 午夜福利网站1000一区二区三区| 久久人人爽人人爽人人片va| 麻豆成人午夜福利视频| 日韩在线高清观看一区二区三区| 91av网一区二区| 亚洲国产精品成人综合色| av在线老鸭窝| 永久免费av网站大全| 97在线视频观看| 又爽又黄无遮挡网站| 亚洲av日韩在线播放| 一二三四中文在线观看免费高清| 18+在线观看网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 男女国产视频网站| 99久久精品一区二区三区| 最近最新中文字幕免费大全7| 欧美极品一区二区三区四区| 亚洲精品成人久久久久久| 亚洲成人久久爱视频| 夜夜看夜夜爽夜夜摸| 久久韩国三级中文字幕| 国产午夜精品论理片| 少妇的逼好多水| 国产 一区 欧美 日韩| 婷婷色综合大香蕉| 国产精品一二三区在线看| 日本黄大片高清| 亚洲高清免费不卡视频| 亚洲最大成人中文| 成人特级av手机在线观看| 真实男女啪啪啪动态图| 久久精品久久久久久噜噜老黄 | 久久久久久久久大av| 中文字幕制服av| 中文欧美无线码| 久99久视频精品免费| 精品国产露脸久久av麻豆 | 亚洲av免费高清在线观看| 夫妻性生交免费视频一级片| 色播亚洲综合网| 色5月婷婷丁香| 床上黄色一级片| 国产精品麻豆人妻色哟哟久久 | 精品久久久久久久末码| 99久久九九国产精品国产免费| 最近视频中文字幕2019在线8| 欧美一级a爱片免费观看看| 免费av毛片视频| 欧美成人免费av一区二区三区| 国产三级在线视频| 干丝袜人妻中文字幕| 欧美又色又爽又黄视频| 三级国产精品片| 男女啪啪激烈高潮av片| 熟女人妻精品中文字幕| 国产精品综合久久久久久久免费| 成人av在线播放网站| 亚洲中文字幕一区二区三区有码在线看| 午夜福利在线观看吧| 日韩,欧美,国产一区二区三区 | 九九在线视频观看精品| 久久久久九九精品影院| 精品熟女少妇av免费看| 亚洲av免费高清在线观看| 插逼视频在线观看| 亚洲欧美清纯卡通| 亚洲中文字幕日韩| 毛片女人毛片| 国产免费一级a男人的天堂| videossex国产| 国产毛片a区久久久久| 精品免费久久久久久久清纯| 亚洲国产精品合色在线| 免费搜索国产男女视频| av福利片在线观看| 欧美xxxx黑人xx丫x性爽| 久久久久性生活片| 身体一侧抽搐| 精品久久久久久久末码| 美女脱内裤让男人舔精品视频| 日日摸夜夜添夜夜爱| 国产精品99久久久久久久久| 69人妻影院| 不卡视频在线观看欧美| 校园人妻丝袜中文字幕| 久久久久国产网址| 欧美激情在线99| www.av在线官网国产| 国产视频首页在线观看| 亚洲av中文字字幕乱码综合| 色综合站精品国产| 日韩成人伦理影院| 日本爱情动作片www.在线观看| 色5月婷婷丁香| 偷拍熟女少妇极品色| av在线亚洲专区| 深爱激情五月婷婷| 人妻制服诱惑在线中文字幕| 嫩草影院入口| 三级男女做爰猛烈吃奶摸视频| 亚洲精品影视一区二区三区av| 岛国毛片在线播放| 岛国在线免费视频观看| 国产精品一区二区在线观看99 | 国内精品宾馆在线| 欧美一区二区精品小视频在线| 午夜爱爱视频在线播放| 国产极品天堂在线| 色视频www国产| 欧美日韩在线观看h| 婷婷色av中文字幕| 日韩精品有码人妻一区| 中文字幕亚洲精品专区| 亚洲av成人精品一二三区| 国产91av在线免费观看| 久久久久久伊人网av| 国产又黄又爽又无遮挡在线| 长腿黑丝高跟| 18+在线观看网站| 亚洲av电影在线观看一区二区三区 | 国产精品嫩草影院av在线观看| 日本黄大片高清| 国产片特级美女逼逼视频| 日本免费在线观看一区| 成人性生交大片免费视频hd| 国产精品国产三级专区第一集| 秋霞在线观看毛片| 黄片无遮挡物在线观看| 超碰97精品在线观看| 五月伊人婷婷丁香| 亚洲久久久久久中文字幕| 国产不卡一卡二| 亚洲国产精品久久男人天堂| 99久国产av精品| 婷婷色麻豆天堂久久 | 99久久精品热视频| 国产成人精品一,二区| 亚洲自偷自拍三级| 亚洲丝袜综合中文字幕| 国产精品精品国产色婷婷| 国产免费男女视频| 欧美激情在线99| 免费av不卡在线播放| 亚洲精品国产av成人精品| 中文亚洲av片在线观看爽| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕| 国产精品久久久久久久久免| 美女国产视频在线观看| 一级毛片aaaaaa免费看小| 国产视频首页在线观看| 少妇熟女aⅴ在线视频| 久久精品久久精品一区二区三区| 国产黄片美女视频| 2021少妇久久久久久久久久久| 中文天堂在线官网| 欧美丝袜亚洲另类| 日韩av在线免费看完整版不卡| 狂野欧美白嫩少妇大欣赏| 2021天堂中文幕一二区在线观| 如何舔出高潮| 色综合亚洲欧美另类图片| 日本三级黄在线观看| 日韩大片免费观看网站 | 免费av毛片视频| 久久人妻av系列| 亚洲国产欧美在线一区| 国产精品不卡视频一区二区| 天美传媒精品一区二区| 国语自产精品视频在线第100页| 18禁裸乳无遮挡免费网站照片| 三级国产精品欧美在线观看| 午夜激情福利司机影院| 日本色播在线视频| 中文在线观看免费www的网站| 高清日韩中文字幕在线| 成人一区二区视频在线观看| 乱人视频在线观看| 网址你懂的国产日韩在线| 一边摸一边抽搐一进一小说| ponron亚洲| 国产69精品久久久久777片| 免费观看a级毛片全部| 蜜桃亚洲精品一区二区三区| 久久99热这里只有精品18| 丝袜喷水一区| 狂野欧美激情性xxxx在线观看| 欧美日韩国产亚洲二区| 国产精品久久视频播放| 成年女人看的毛片在线观看| 好男人在线观看高清免费视频| 又粗又爽又猛毛片免费看| 色视频www国产| 国内揄拍国产精品人妻在线| 国产一区有黄有色的免费视频 | 少妇裸体淫交视频免费看高清| 亚洲av成人精品一二三区| 亚洲av.av天堂| 久久99热6这里只有精品| 两个人的视频大全免费| 2021少妇久久久久久久久久久| 色网站视频免费| 天美传媒精品一区二区| 一级毛片aaaaaa免费看小| 日韩制服骚丝袜av| 18禁动态无遮挡网站| 日本免费一区二区三区高清不卡| 一区二区三区免费毛片| 一个人看的www免费观看视频| 欧美zozozo另类| 国产精品人妻久久久久久| 亚洲欧美成人综合另类久久久 | 国产探花在线观看一区二区| 国产一区亚洲一区在线观看| 国产乱人偷精品视频| 亚洲av成人精品一二三区| 亚洲人成网站在线播| av免费观看日本| 嘟嘟电影网在线观看| 精华霜和精华液先用哪个| 免费看日本二区| av播播在线观看一区| 啦啦啦啦在线视频资源| 纵有疾风起免费观看全集完整版 | 激情 狠狠 欧美| 麻豆一二三区av精品| 精品久久久久久久久av| 欧美97在线视频| 久久久久久国产a免费观看| 国产成人午夜福利电影在线观看| 有码 亚洲区| 波多野结衣巨乳人妻| 久久6这里有精品| 身体一侧抽搐| 久热久热在线精品观看| 丰满人妻一区二区三区视频av| 精品久久久久久久人妻蜜臀av| 一区二区三区四区激情视频| 成人毛片60女人毛片免费| 国产视频内射| 菩萨蛮人人尽说江南好唐韦庄 | 午夜a级毛片| 女的被弄到高潮叫床怎么办| 久久热精品热| 精品久久久久久久久av| 免费观看的影片在线观看| 日本黄大片高清| 国产美女午夜福利| 又粗又爽又猛毛片免费看| 国产黄色视频一区二区在线观看 | 日日干狠狠操夜夜爽| 网址你懂的国产日韩在线| 一个人观看的视频www高清免费观看| 熟妇人妻久久中文字幕3abv| 久久久久久大精品| 精品久久久噜噜| 最近最新中文字幕免费大全7| 十八禁国产超污无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 热99在线观看视频| 日本三级黄在线观看| 国产精品久久久久久精品电影| videos熟女内射| 亚洲自偷自拍三级| 一级毛片aaaaaa免费看小| 国产精品三级大全| 国产激情偷乱视频一区二区| 成人特级av手机在线观看| 欧美日本视频| 国产精品日韩av在线免费观看| 日韩一区二区视频免费看| 欧美一区二区国产精品久久精品| 中文字幕久久专区| 亚洲国产精品专区欧美| 尤物成人国产欧美一区二区三区| 亚洲精品日韩av片在线观看| 蜜桃久久精品国产亚洲av| 亚洲欧美成人综合另类久久久 | 色噜噜av男人的天堂激情| 欧美成人精品欧美一级黄| 久久久久精品久久久久真实原创| 国产av一区在线观看免费| 美女国产视频在线观看| 久久精品国产亚洲av天美| 最新中文字幕久久久久| 亚洲欧美精品综合久久99| 色噜噜av男人的天堂激情| 在线播放国产精品三级| 七月丁香在线播放| 男女那种视频在线观看| 亚洲内射少妇av| 免费av观看视频| 欧美性猛交╳xxx乱大交人| 亚洲欧美日韩东京热| 99国产精品一区二区蜜桃av| 在线观看66精品国产| 国模一区二区三区四区视频| 国产高潮美女av| 国产精品国产三级国产专区5o | 麻豆av噜噜一区二区三区| 深夜a级毛片| 久久精品国产亚洲av天美| 国产成年人精品一区二区| 亚洲av中文字字幕乱码综合| 免费看美女性在线毛片视频| 女人久久www免费人成看片 | 在线天堂最新版资源| 色视频www国产| 亚洲av成人av| 免费无遮挡裸体视频| 久久久久网色| 在线观看一区二区三区| 国产成人午夜福利电影在线观看| 午夜免费激情av| 亚洲成人精品中文字幕电影| 青春草视频在线免费观看| 国产成人a∨麻豆精品| 亚洲精品自拍成人| 国产在视频线在精品| 美女xxoo啪啪120秒动态图| 日本黄大片高清| 九九爱精品视频在线观看| 91久久精品国产一区二区三区| 国产老妇女一区| 91精品国产九色| 97超视频在线观看视频| 不卡视频在线观看欧美| 国产黄色视频一区二区在线观看 | 国产精品人妻久久久久久| 日本一二三区视频观看| 久久精品国产自在天天线| 欧美高清成人免费视频www| 亚洲四区av| 亚洲欧美成人综合另类久久久 | 亚洲国产最新在线播放| 久久国产乱子免费精品| 可以在线观看毛片的网站| 亚洲aⅴ乱码一区二区在线播放| 神马国产精品三级电影在线观看| videos熟女内射| 成人午夜高清在线视频| 亚洲精品久久久久久婷婷小说 | 最新中文字幕久久久久| 亚洲国产日韩欧美精品在线观看| 欧美色视频一区免费|