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

    管殼式太陽能相變儲(chǔ)熱器傳熱特性的數(shù)值研究

    2023-04-12 00:00:00韓濤馬彥花方嘉賓閆皓冰魏進(jìn)家董思遠(yuǎn)
    太陽能學(xué)報(bào) 2023年3期
    關(guān)鍵詞:相變材料肋片太陽能

    DOI:10.19912/j.0254-0096.tynxb.2021-1322 文章編號(hào):0254-0096(2023)03-0525-08

    摘 要:通過焓-孔隙法建立數(shù)值模型,研究石蠟RT50在管殼式太陽能儲(chǔ)熱器的傳熱特性,并分析儲(chǔ)熱器安裝形式、不同管徑和管外肋片結(jié)構(gòu)對(duì)相變儲(chǔ)熱過程的影響規(guī)律。結(jié)果表明,立式儲(chǔ)熱器在管徑較大時(shí)(相變材料區(qū)域厚度較薄時(shí))熔化速度比臥式儲(chǔ)熱器快;而在管徑較小時(shí)(相變材料區(qū)域較厚時(shí))熔化速度比臥式儲(chǔ)熱器慢。臥式儲(chǔ)熱器相變材料熔化過程的初始時(shí)間段和末尾時(shí)間段熱傳導(dǎo)占主導(dǎo)地位,而中間時(shí)間段自然對(duì)流換熱占主導(dǎo)地位,且熔化速度主要由自然對(duì)流換熱決定。此外,管外側(cè)安裝肋片能夠顯著提高傳熱系數(shù),縮短熔化時(shí)間。相較于無肋片管道,臥式儲(chǔ)熱器采用環(huán)肋管道可縮短熔化時(shí)長(zhǎng)31.6%;采用直肋管道可縮短熔化時(shí)長(zhǎng)42.1%。直肋比環(huán)肋具有更好的傳熱效果,主要原因?yàn)橹崩吣軌驈?qiáng)化熔化過程中的自然對(duì)流。

    關(guān)鍵詞:太陽能;相變材料;熔化;肋片;自然對(duì)流;熱能儲(chǔ)存

    中圖分類號(hào):TK513.5 " "文獻(xiàn)標(biāo)志碼:A

    0 引 言

    近年來太陽能的開發(fā)利用成為國(guó)內(nèi)外的研究熱點(diǎn)。由于太陽能分布不均、能流密度低以及間歇性日照的特點(diǎn),儲(chǔ)能環(huán)節(jié)成大多數(shù)太陽能應(yīng)用過程中的關(guān)鍵[1]。國(guó)家能源局在2016年啟動(dòng)的中國(guó)第一批太陽能熱發(fā)電示范項(xiàng)目中明確要求太陽能儲(chǔ)熱設(shè)備必須具備最少4 h的儲(chǔ)熱能力?,F(xiàn)有的儲(chǔ)能方式主要分為3種:顯熱儲(chǔ)能、潛熱儲(chǔ)能、熱化學(xué)儲(chǔ)能。其中,潛熱儲(chǔ)能是利用相變材料(phase change material, PCM)發(fā)生相變時(shí)吸熱或放熱的特性儲(chǔ)存或釋放能量,也叫相變儲(chǔ)熱。相比之下,潛熱儲(chǔ)能能量?jī)?chǔ)存密度較大、化學(xué)穩(wěn)定性好、應(yīng)用溫度范圍較廣,被認(rèn)為是最有應(yīng)用前景的儲(chǔ)能方式之一[2]。

    管殼式太陽能相變儲(chǔ)熱系統(tǒng)由于結(jié)構(gòu)簡(jiǎn)單,應(yīng)用十分廣泛。近年來,眾多學(xué)者通過理論、實(shí)驗(yàn)和數(shù)值等方法對(duì)管殼式相變儲(chǔ)熱系統(tǒng)傳熱性能進(jìn)行了大量研究。Lacroix[3]建立了一個(gè)理論模型來預(yù)測(cè)帶有PCM管殼式儲(chǔ)熱系統(tǒng)的瞬態(tài)行為,結(jié)果表明PCM的導(dǎo)熱系數(shù)低是潛熱儲(chǔ)能系統(tǒng)的一個(gè)基本缺點(diǎn),這導(dǎo)致PCM系統(tǒng)的低傳熱速率和延遲響應(yīng)。劉洋等[4]通過實(shí)驗(yàn)研究了不同熱源溫度與環(huán)境溫度下儲(chǔ)能散熱過程中儲(chǔ)能單元的熱性能,結(jié)果發(fā)現(xiàn)添加泡沫銅對(duì)儲(chǔ)能單元的儲(chǔ)散熱具有較好的優(yōu)化效果。

    對(duì)于相變儲(chǔ)熱來說,通過數(shù)值建模的方式,易于優(yōu)化容器的尺寸和形狀以及評(píng)估儲(chǔ)熱系統(tǒng)的性能,且無需大量的研究經(jīng)費(fèi)[5]。張仲彬等[6]建立相變蓄熱膠囊的三維有序及無序堆積數(shù)值模型,發(fā)現(xiàn)當(dāng)采用有序堆積時(shí)系統(tǒng)的各項(xiàng)性能指標(biāo)最優(yōu)。李椿等[7]建立殼管式相變儲(chǔ)能換熱器仿真模型,研究了蓄熱過程中熱源管壁溫、熱管口徑規(guī)格對(duì)傳熱效率的影響。

    與此同時(shí),也有一些學(xué)者研究了PCM流動(dòng)及肋片強(qiáng)化對(duì)相變儲(chǔ)熱過程的影響。Hosseini等[8]采用實(shí)驗(yàn)和數(shù)值相結(jié)合的方法研究了管殼式換熱器內(nèi)石蠟RT50熔化過程中浮力驅(qū)動(dòng)對(duì)流的作用,結(jié)果表明適當(dāng)提高進(jìn)水溫度可有效減少傳熱時(shí)間。Pahamli等[9]通過數(shù)值模擬的方法研究了偏心度、流量大小以及進(jìn)口水溫等參數(shù)對(duì)PCM熔化過程的影響。田楊等[10]數(shù)值模擬研究不同轉(zhuǎn)速下加肋旋轉(zhuǎn)的蓄熱器的蓄熱特性,研究表明加肋并旋轉(zhuǎn)可顯著提升相變蓄熱器的蓄熱效果,縮短蓄熱時(shí)間。Sciacovelli等[11]使用不同的樹形肋片來提高管殼式潛熱儲(chǔ)能裝置的性能,同時(shí)提出一種改進(jìn)瞬態(tài)工況系統(tǒng)的形狀優(yōu)化策略。

    本文通過焓-孔隙法建立數(shù)值模型,研究不同形式太陽能儲(chǔ)熱器(立式和臥式)、不同管徑及管外加裝不同肋片等參數(shù)對(duì)傳熱過程PCM傳熱熔化行為的影響,進(jìn)而闡明相變儲(chǔ)熱過程流動(dòng)傳熱強(qiáng)化機(jī)理,為相變儲(chǔ)熱過程強(qiáng)化設(shè)計(jì)提供參考。

    1 物理及數(shù)值模型

    1.1 物理模型

    1.1.1 立式、臥式換熱器儲(chǔ)熱單元

    目前大多數(shù)國(guó)內(nèi)外對(duì)管殼式儲(chǔ)能換熱器的研究主要集中在設(shè)計(jì)及運(yùn)行參數(shù),對(duì)于儲(chǔ)熱器形式的影響研究缺乏對(duì)比。因此本文針對(duì)不同形式的管殼式儲(chǔ)能換熱器對(duì)相變儲(chǔ)熱過程的影響進(jìn)行探究。管殼式換熱器儲(chǔ)熱單元形狀如圖1所示,其整體長(zhǎng)度為600 mm,根據(jù)擺放形式分為臥式(圖1a)和立式(圖1b)。重力作用為豎直向下。傳熱流體(heat transfer fluid,HTF)為水工質(zhì),走管程,其中立式為下口進(jìn),上口出。截面尺寸如圖1c所示。殼側(cè)填充的PCM材料為石蠟RT50。表1為RT50的各項(xiàng)熱物性參數(shù)。

    RT50黏度[μ]與其溫度[T]的關(guān)系為式(1):

    式中:[A=0.819],[B=-0.01546]。

    1.1.2 翅片形狀和尺寸

    對(duì)于管殼式換熱器,安裝肋片增加換熱面積以達(dá)到強(qiáng)化傳熱的目的是以往研究中常用的方法。尤其對(duì)于無機(jī)鹽類以及有機(jī)相變材料來說,其導(dǎo)熱系數(shù)通常較小,這也是相變儲(chǔ)熱應(yīng)用中的短板之一。以往研究表明,肋片管能有效提高相變過程中的傳熱效率。有學(xué)者對(duì)肋片的高度間距進(jìn)行探究,也有學(xué)者設(shè)計(jì)了各類肋片結(jié)構(gòu)以提高傳熱效率。但大多新型肋片設(shè)計(jì)較為復(fù)雜,難以得到廣泛應(yīng)用;而對(duì)于常見的一些工程應(yīng)用中的肋片類型,其傳熱強(qiáng)化效果缺乏定性定量的認(rèn)識(shí)。本文在維持可力片表面積與體積不變的前提下,分別研究直肋與環(huán)肋對(duì)相變材料熔化過程的強(qiáng)化作用,并進(jìn)行對(duì)比分析。為了便于比較,在圖1a臥式儲(chǔ)熱器管上加裝直肋與環(huán)肋,其形狀尺寸分別如圖2、圖3所示。

    1.2 數(shù)值模型

    1.2.1 焓-孔隙率法模型

    已有研究表明,在相變儲(chǔ)熱過程中,對(duì)流傳熱效果不可忽略。在確保結(jié)果準(zhǔn)確性的基礎(chǔ)上,為了簡(jiǎn)化數(shù)值模型,本文提出以下假設(shè):1)考慮儲(chǔ)熱器保溫效果,儲(chǔ)熱器殼壁面設(shè)置為絕熱,并考慮管壁厚度;2)PCM熔化與凝固過程中考慮自然對(duì)流效應(yīng);3)水工質(zhì)和相變材料的熱物理性質(zhì)被認(rèn)為與溫度無關(guān),均勻且各向同性。

    對(duì)于相變過程的求解,應(yīng)用較為廣泛的是焓-孔隙法模型。其通過對(duì)式(4)~式(7)焓值的分布進(jìn)行求解,從而確定相界面??刂品匠贪ㄟB續(xù)性方程、動(dòng)量方程和能量方程。

    連續(xù)性方程如式(2)所示:

    式中:[▽]——矢量微分算子;[V]——速度矢量。

    動(dòng)量方程如式(3)所示:

    式中:[t]——時(shí)間,s;[ρ]——密度,kg/m3;p——壓強(qiáng),Pa;[μ]——?jiǎng)恿︷ざ?,Pa·s;[β]——熱膨脹系數(shù),1/K;[g]——重力加速度,m/s2;[T]——溫度,K;[Tref]——參考溫度,K;[S]——源項(xiàng)。

    能量方程如式(4)所示:

    式中:[cp]——比熱容,J/(kg·K);[k]——導(dǎo)熱系數(shù),W/(m·K);[hsens]——顯熱焓,J/kg;[hlat]——潛熱焓,J/kg。

    顯熱焓可表示為式(5):

    潛熱焓可表示為式(6):

    式中:[L]——完全熔化潛熱;[λ]——液相分?jǐn)?shù),[λ]可表示為式(7):

    式中:[Ts]——固體溫度,K;[Tliq]——液體溫度,K。動(dòng)量方程中的源項(xiàng)[S]為達(dá)西定律的阻尼項(xiàng),其表達(dá)式為式(8):

    式中:[Amush]——糊狀區(qū)常數(shù),其值通常在104~107之間。本文取106計(jì)算求解。

    1.2.2 邊界及初始條件

    PCM的初始溫度[T0]為[25 ℃];HTF進(jìn)口溫度[Tf,i]設(shè)為70 ℃,入口流速為0.044 m/s;儲(chǔ)能單元?dú)け诿嬖O(shè)置為絕熱。

    2 結(jié)果與討論

    2.1 立、臥式儲(chǔ)熱單元傳熱過程

    對(duì)于立式和臥式儲(chǔ)熱單元,采用相同數(shù)量的結(jié)構(gòu)化網(wǎng)格及時(shí)間步長(zhǎng)計(jì)算,其結(jié)果如圖4所示。圖4中立式儲(chǔ)熱單元內(nèi)PCM熔化時(shí)間比臥式更短,但并不意味立式儲(chǔ)熱單元的絕對(duì)優(yōu)越性。當(dāng)采用管內(nèi)徑[di]為6 mm、外徑[do]為10 mm,而其他條件不變的儲(chǔ)熱單元做實(shí)驗(yàn)對(duì)比,此時(shí)PCM熔化的過程如圖5所示。對(duì)比圖4、圖5可看到,當(dāng)采用更小的管徑時(shí)(相當(dāng)于增加相變材料區(qū)域厚度),臥式儲(chǔ)熱單元內(nèi)RT50熔化速度更快。這說明立式儲(chǔ)熱單元內(nèi)相變材料熔化過程受相變材料區(qū)域厚度影響較大。對(duì)于立式儲(chǔ)熱單元來講,由于自然對(duì)流發(fā)生在豎直方向,取豎直方向截面流動(dòng)及液相云圖如圖6所示。由圖6可看到,整個(gè)相變?nèi)刍^程中流動(dòng)區(qū)域、流動(dòng)趨勢(shì)未發(fā)生明顯變化。但由于對(duì)流傳熱主要發(fā)生在固液兩相區(qū)及其附近,而導(dǎo)熱主要將熱量通過管外壁導(dǎo)向固液兩相區(qū)及固相區(qū),熔化初始階段固液兩相區(qū)及固相區(qū)較大,導(dǎo)熱作用較強(qiáng),而固液兩相區(qū)通過導(dǎo)熱快速產(chǎn)生后,在熔化過程中間階段兩相區(qū)域大小并未明顯增大。因此在熔化前期,導(dǎo)熱效果減弱,對(duì)流有所增強(qiáng),但是熔化速度并未顯著放緩,而在60 min后,固相區(qū)消失,固液兩相區(qū)逐漸減小,但導(dǎo)熱與對(duì)流作用同時(shí)減弱,因而熔化速度迅速減小。而對(duì)于臥式儲(chǔ)熱單元,其PCM熔化速度經(jīng)歷了復(fù)雜的變化過程。結(jié)合圖7中心截面熔化及流動(dòng)情況分析,發(fā)現(xiàn)熔化過程按傳熱情況可分為3個(gè)階段:初始時(shí)間段導(dǎo)熱主導(dǎo)區(qū)、中間時(shí)間段自然對(duì)流主導(dǎo)區(qū)及末尾時(shí)間段導(dǎo)熱主導(dǎo)區(qū)。

    1)初始階段,直接接觸HTF的PCM薄層由于熱傳導(dǎo)作用快速升溫熔化,這個(gè)過程持續(xù)時(shí)間很短。由于PCM的導(dǎo)熱系數(shù)較小,當(dāng)液相層到達(dá)一定厚度時(shí),其熱傳導(dǎo)過程發(fā)生急劇惡化,液相變化趨于平緩。

    2)中間階段,隨著傳熱過程的進(jìn)行,液相區(qū)逐漸擴(kuò)大。由于溫度差導(dǎo)致的熱膨脹效果,液相區(qū)及部分固液相混合區(qū)發(fā)生自然對(duì)流效應(yīng)。在密度差的驅(qū)動(dòng)下,液相沿著溫度較高的管壁面一側(cè)爬升,沿著溫度較低的固相區(qū)域一側(cè)下降,形成環(huán)形對(duì)流,如圖7所示。在固液混合區(qū)邊界處的自然對(duì)流極大地增加了RT50的熔化速度,同時(shí)固液兩相區(qū)持續(xù)擴(kuò)大。

    3)末尾階段(60 min以后),隨著液相區(qū)的擴(kuò)大,液相線的位置持續(xù)降低到下半?yún)^(qū),固液混合區(qū)與液相之間邊界減小,同時(shí)液相區(qū)內(nèi)密度差減小,對(duì)流流動(dòng)速度減小,且此時(shí)腔體上半部液相區(qū)內(nèi)的流動(dòng)對(duì)傳熱失去強(qiáng)化效果,因此熔化速率急劇減小。最終,隨著固相區(qū)減薄,熱傳導(dǎo)效果增強(qiáng),熔化速率略微增大直至完全熔化。

    2.2 直肋、環(huán)肋傳熱強(qiáng)化效果

    對(duì)于管內(nèi)徑[di]為10 mm、外徑[do]為14 mm的臥式儲(chǔ)熱單元,在管上加裝肋片可極大強(qiáng)化相變儲(chǔ)熱過程,提高傳熱效果與PCM熔化速率。從圖8可看到,用直肋管代替無肋管后PCM熔化速率提高了42.1%,而環(huán)肋管取代無肋管后PCM熔化速率提高了31.6%。顯然,直肋管的效果更好。圖9為直肋管儲(chǔ)熱單元中心截面處的熔化以及流動(dòng)情況示意圖,可看到,加裝直肋之后其內(nèi)部PCM流動(dòng)區(qū)域明顯擴(kuò)大,基本覆蓋整個(gè)區(qū)域,對(duì)流傳熱得到顯著增強(qiáng)。因此PCM熔化速度大幅度提高,極大減少了熔化時(shí)間。圖10為直肋儲(chǔ)熱單元中心截面溫度分布云圖。

    對(duì)于加裝環(huán)肋的情況,考慮到三維對(duì)稱性,在離中心位置0、3、6、9 mm處取4個(gè)截面進(jìn)行觀察,結(jié)果如圖11所示。圖12為各個(gè)截面在傳熱進(jìn)行到15 min時(shí)的液相分布及流動(dòng)情況??煽吹?,取肋片所在的A截面,無流動(dòng)跡線。B截面液相流動(dòng)區(qū)域最大。但隨著距離增大,液相流動(dòng)區(qū)域有一定程度縮小,這與環(huán)肋肋片的附近導(dǎo)熱效果更好有關(guān)。取B截面觀察其不同時(shí)刻的流動(dòng)情況,如圖13所示??煽吹?,其在不同時(shí)刻的流動(dòng)狀況與無肋管的類似,環(huán)肋肋片的存在并沒有擴(kuò)大其對(duì)流流動(dòng)域,因此對(duì)于對(duì)流傳熱的強(qiáng)化十分有限,可得出,環(huán)肋對(duì)于傳熱過程傳熱強(qiáng)化主要在熱傳導(dǎo)方面,而直肋可強(qiáng)化熱對(duì)流及熱傳導(dǎo)過程。

    2.3 儲(chǔ)熱器內(nèi)部自然對(duì)流形成機(jī)理

    對(duì)于直肋管儲(chǔ)熱單元來說,其截面流動(dòng)情況可簡(jiǎn)化為圖14所示,開始傳熱后,熱的肋片和管壁與冷的殼壁形成溫度勢(shì)差,PCM熔化后沿?zé)岜诿媾郎?,而遇到冷的壁面下降,在相?duì)封閉的空間形成一個(gè)完整的流動(dòng)循環(huán),從而大幅度強(qiáng)化對(duì)流傳熱。而對(duì)于環(huán)肋管儲(chǔ)熱單元來說,其內(nèi)部PCM流動(dòng)情況可簡(jiǎn)化為圖15所示。由于溫度高的肋片彼此平行,因此

    無法形成溫度勢(shì)差驅(qū)動(dòng)。熱流沿著肋片爬升后,遇到頂部冷的殼壁,只能被迫改變方向,無法在當(dāng)前截面內(nèi)完成循環(huán)流動(dòng),因此與無肋管內(nèi)PCM流動(dòng)狀態(tài)類似。所以在管上加裝環(huán)肋并不能有效強(qiáng)化對(duì)流傳熱過程。

    綜上所述,對(duì)于常見的導(dǎo)熱系數(shù)較小的相變材料,加裝肋片對(duì)熱傳導(dǎo)強(qiáng)化效果有限,應(yīng)當(dāng)考慮其對(duì)自然對(duì)流傳熱過程的強(qiáng)化。設(shè)計(jì)原則應(yīng)當(dāng)包括,在封閉空間內(nèi),考慮熱壁面與冷壁面豎直相對(duì),以便于形成流動(dòng)循環(huán)。同時(shí),當(dāng)冷、熱壁面有一定高度時(shí),對(duì)于受熱膨脹的流體來說,熱壁面與冷壁面的溫度梯度應(yīng)當(dāng)豎直向上(如圖16所示),以便熱流可以沿著壁面爬升。

    3 結(jié) 論

    具有儲(chǔ)熱系統(tǒng)的太陽能熱發(fā)電技術(shù)是未來可再生能源中最具應(yīng)用前景的發(fā)電技術(shù)之一,因此對(duì)于儲(chǔ)熱系統(tǒng)的研究和改進(jìn)具有重要意義。本文通過焓-孔隙法建立管殼式太陽能儲(chǔ)熱器數(shù)值模型,在考慮PCM密度隨溫度變化的基礎(chǔ)上,研究RT50在不同管殼式儲(chǔ)熱器中的熔化行為,并探究相變儲(chǔ)能傳熱過程的結(jié)構(gòu)優(yōu)化方法。得到的結(jié)論如下:

    1) 對(duì)比研究立式和臥式換熱器儲(chǔ)熱單元內(nèi)PCM傳熱熔化過程,發(fā)現(xiàn)臥式儲(chǔ)熱單元內(nèi)PCM熔化過程可分為3個(gè)階段:初始時(shí)間段導(dǎo)熱主導(dǎo),中間時(shí)間段自然對(duì)流主導(dǎo)及末尾時(shí)間段導(dǎo)熱主導(dǎo)。而立式儲(chǔ)熱單元內(nèi)PCM熔化速度對(duì)相變材料區(qū)域的厚度更加敏感。

    2) 對(duì)于翅片強(qiáng)化過程研究表明,在臥式管殼換熱器內(nèi)管上加裝直肋后PCM完全熔化時(shí)間縮短了42.1%,而加裝環(huán)肋完全熔化時(shí)間縮短了31.6%,直肋管的效果比環(huán)肋管更好,這是由于直肋強(qiáng)化了PCM熔化過程中的自然對(duì)流傳熱。

    3) 對(duì)于一些常見的相變儲(chǔ)熱材料,其導(dǎo)熱系數(shù)通常很小,肋片等結(jié)構(gòu)對(duì)傳熱導(dǎo)熱過程的強(qiáng)化效果有限,應(yīng)當(dāng)考慮其自然對(duì)流傳熱強(qiáng)化傳熱過程。

    [參考文獻(xiàn)]

    [1] 韓宗偉, 王一茹, 阿不來提·依米提, 等. 太陽能熱泵相變蓄熱供暖系統(tǒng)參數(shù)影響研究[J]. 太陽能學(xué)報(bào), 2015,36(8): 2028-2035.

    HAN Z W, WANG Y R, YIMIT A, et al. Study on influence of parameters of solar energy heat pump system on" phase" change" heat" storage[J]." Acta" energiae" solaris sinica, 2015, 36(8): 2028-2035.

    [2] AGYENIM F, HEWITT N, EAMES P, et al. A review of materials, heat transfer and phase change problem formulation for latent heat thermal energy storage systems (LHTESS)[J]." " Renewable" "and" "sustainable" "energy reviews, 2010, 14(2): 615-628.

    [3] LACROIX M. Numerical simulation of a shell-and-tube latent heat thermal energy storage unit[J]. Solar energy, 1993, 50(4): 357-367.

    [4] 劉洋, 段建國(guó), 賀秀芬, 等. 低溫相變儲(chǔ)能單元強(qiáng)化傳熱特性實(shí)驗(yàn)[J]. 太陽能學(xué)報(bào), 2020, 41(6): 335-340.

    LIU Y, DUAN J G, HE X F, et al. Experiments on enhanced heat transfer performance of low temperature latent thermal energy storage unit[J]. Acta energiae solaris sinica, 2020, 41(6): 335-340.

    [5] VISKANTA R. Phase-change heat transfer//Solar heat storage: latent heat materials[M]. CRC Press, 2018: 153-222.

    [6] 張仲彬, 孟雨欣, 劉永強(qiáng). 相變膠囊堆積形態(tài)對(duì)蓄熱系統(tǒng)性能的影響分析[J]. 太陽能學(xué)報(bào), 2021, 42(2): 37-42.

    ZHANG Z B,MENG Y X,LIU Y Q. Influence of phase change capsule packing configuration on performance of" "heat" storage" system[J]." Acta" energiae" solaris" sinica, 2021, 42(2): 37-42.

    [7] 李椿, 王志華, 王建春, 等. 殼管式相變儲(chǔ)能換熱器性能研究與場(chǎng)協(xié)同效應(yīng)分析[J]. 太陽能學(xué)報(bào), 2020, 41(3): 226-233.

    LI C, WANG Z H, WANG J C, et al. Performance study and field synergy analysis of shell and tube phase change energy" storage" heat" exchanger[J]. Acta" energiae" solaris sinica, 2020, 41(3): 226-233.

    [8] HOSSEINI M J, RANJBAR A A, SEDIGHI K, et al. A combined experimental and computational study on the melting behavior of a medium temperature phase change storage material inside shell and tube heat exchanger[J]. International communications in heat and mass transfer, 2012, 39(9): 1416-1424.

    [9] PAHAMLI Y, HOSSEINI M J, RANJBAR A A, et al. Analysis of the effect of eccentricity and operational parameters in PCM-filled single-pass shell and tube heat exchangers[J]. Renewable energy, 2016, 97: 344-357.

    [10] 田揚(yáng), 趙明, 胡明禹, 等. 加肋旋轉(zhuǎn)對(duì)相變蓄熱器蓄熱性能的影響及場(chǎng)協(xié)同分析[J]. 太陽能學(xué)報(bào), 2021, 42(3): 395-400.

    TIAN Y, ZHAO M, HU M Y, et al. Effect of ribbed rotation on heat storage performance of phase change thermal storage unit and field synergy analysis[J]. Acta energiae solaris sinica, 2021, 42(3): 395-400.

    [11] SCIACOVELLI A, GAGLIARDI F, VERDA V. Maximization of performance of a PCM latent heat storage system with innovative fins[J]. Applied energy, 2015, 137: 707-715.

    NUMERICAL SIMULATION STUDY OF HEAT TRANSFER CHARACTERISTICS ON SOLAR TUBE-AND-SHELL PHASE

    CHANGE HEAT STORAGE UNIT

    Han Tao1,Ma Yanhua2,F(xiàn)ang Jiabin1,Yan Haobing1,Wei Jinjia1,Dong Siyuan2

    (1. School of Chemical Engineering and Technology, Xi’an Jiaotong University, Xi’an 710049, China;

    2. Faculty of Mechanical and Electrical Engineering, Yunnan Agricultural University, Kunming 650201, China)

    Abstract:A numerical model is established using the enthalpy-porosity approach to study the heat transfer characteristics of a tube-and-shell phase-change heat exchanger filled with paraffin wax RT50. The influence of exchanger placement forms, tube diameters and fin structures on the phase change process of RT50 is analyzed. The results depicts that the vertical heat exchanger has a faster melting rate than the horizontal one as the tube diameter is large. However, the opposite results are obtained in case the tube diameter is small. For the horizontal exchanger, the heat conduction is dominant at the beginning and end of the melting process, while the natural convection plays a more important role at the intermediate stage of melting. Besides, the duration of the melting is mainly determined by the natural convection. In addition, the installation of fins on the outside of the tube can significantly increase the heat transfer coefficient and shorten the melting time. Compared with finless tube, the utilization of ring-fin tube in the horizontal heat storage device can shorten the melting time by 31.6%, and the utilization of straight-fin tube can shorten the melting time by 42.1%. Therefore, the straight fins have a better effect than the ring fins. The main reason is that straight fins can strengthen the natural convection during the melting process.

    Keywords:solar energy; phase change material; melting; fins; natural convection; thermal energy storage

    猜你喜歡
    相變材料肋片太陽能
    應(yīng)用廣泛的太陽能無人機(jī)
    基于Icepak的IGBT模塊水冷散熱特性仿真分析
    肋片形式對(duì)Hampson型換熱器流動(dòng)換熱特性影響
    真空與低溫(2021年2期)2021-03-29 11:25:38
    三維矩形外肋管流動(dòng)傳熱特性數(shù)值模擬
    山東化工(2020年17期)2020-10-23 04:31:40
    太陽能可以這樣玩
    微膠囊相變材料的研究進(jìn)展
    關(guān)于太陽能相變蓄熱系統(tǒng)的研究與分析
    相變材料應(yīng)用于建筑節(jié)能探討
    電化學(xué)法納米線的制備
    太陽能虛擬窗
    亚洲伊人色综图| 视频区欧美日本亚洲| 精品一区二区三区视频在线观看免费 | 最好的美女福利视频网| 色尼玛亚洲综合影院| 亚洲国产精品一区二区三区在线| 神马国产精品三级电影在线观看 | 免费高清在线观看日韩| 久久久久久人人人人人| 一级a爱视频在线免费观看| 欧美一级毛片孕妇| 欧美日韩中文字幕国产精品一区二区三区 | 丁香欧美五月| netflix在线观看网站| 日韩有码中文字幕| 亚洲熟妇熟女久久| ponron亚洲| 亚洲精品国产一区二区精华液| 丝袜美腿诱惑在线| 国产av在哪里看| 欧美乱色亚洲激情| 久久久久国内视频| 精品一区二区三卡| 少妇 在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲一区二区三区不卡视频| 成人手机av| 国产成年人精品一区二区 | 一区福利在线观看| 欧美日韩瑟瑟在线播放| 麻豆久久精品国产亚洲av | 水蜜桃什么品种好| 高潮久久久久久久久久久不卡| 国产精品偷伦视频观看了| 亚洲av日韩精品久久久久久密| 免费看a级黄色片| 在线观看66精品国产| 国产亚洲欧美精品永久| 99国产精品一区二区蜜桃av| 宅男免费午夜| 天堂中文最新版在线下载| 91成人精品电影| 精品福利观看| 亚洲av成人一区二区三| 精品欧美一区二区三区在线| 亚洲五月婷婷丁香| 久久精品人人爽人人爽视色| 美女大奶头视频| 麻豆一二三区av精品| 日本免费a在线| 巨乳人妻的诱惑在线观看| 岛国在线观看网站| 人妻丰满熟妇av一区二区三区| 久久精品国产99精品国产亚洲性色 | 波多野结衣高清无吗| 亚洲午夜理论影院| 国产高清videossex| aaaaa片日本免费| 久久中文字幕一级| av视频免费观看在线观看| 久久久久久人人人人人| 欧美不卡视频在线免费观看 | 久久久久国内视频| av有码第一页| 亚洲国产看品久久| 久久久久久人人人人人| 色老头精品视频在线观看| 亚洲情色 制服丝袜| 国产区一区二久久| 99在线视频只有这里精品首页| 露出奶头的视频| 韩国精品一区二区三区| 侵犯人妻中文字幕一二三四区| 我的亚洲天堂| 在线观看免费午夜福利视频| 最近最新免费中文字幕在线| 久久精品影院6| 久久亚洲精品不卡| 18禁观看日本| 韩国av一区二区三区四区| 久久草成人影院| 两性夫妻黄色片| 极品人妻少妇av视频| 88av欧美| 又紧又爽又黄一区二区| 热99re8久久精品国产| 午夜福利在线观看吧| 黄色视频,在线免费观看| 国产免费男女视频| 亚洲国产中文字幕在线视频| 一级,二级,三级黄色视频| 亚洲欧美一区二区三区久久| 国产一区二区三区综合在线观看| 婷婷精品国产亚洲av在线| 日韩精品中文字幕看吧| 日日夜夜操网爽| 正在播放国产对白刺激| 欧美日韩亚洲国产一区二区在线观看| 这个男人来自地球电影免费观看| 国产亚洲精品久久久久5区| 亚洲成av片中文字幕在线观看| av天堂在线播放| 国产99白浆流出| 久久久久国产精品人妻aⅴ院| 在线观看舔阴道视频| 精品一区二区三区av网在线观看| 久久久久国产一级毛片高清牌| 精品无人区乱码1区二区| 久久天躁狠狠躁夜夜2o2o| 国产伦人伦偷精品视频| 国产欧美日韩精品亚洲av| 一级a爱片免费观看的视频| 亚洲精品在线观看二区| 色哟哟哟哟哟哟| 夜夜看夜夜爽夜夜摸 | 久久精品91蜜桃| 九色亚洲精品在线播放| 亚洲自偷自拍图片 自拍| 国产亚洲精品综合一区在线观看 | 亚洲欧美激情综合另类| www.精华液| 99精品欧美一区二区三区四区| 国产黄a三级三级三级人| 18禁观看日本| 国产亚洲精品第一综合不卡| 久久这里只有精品19| 欧美日韩乱码在线| 精品熟女少妇八av免费久了| 69av精品久久久久久| 亚洲九九香蕉| 欧美日韩亚洲综合一区二区三区_| 人成视频在线观看免费观看| 热99国产精品久久久久久7| 99精品在免费线老司机午夜| 怎么达到女性高潮| 高潮久久久久久久久久久不卡| 免费久久久久久久精品成人欧美视频| 亚洲第一欧美日韩一区二区三区| 99久久精品国产亚洲精品| 黄色片一级片一级黄色片| 久久这里只有精品19| 极品人妻少妇av视频| 一进一出抽搐动态| 一区在线观看完整版| 香蕉国产在线看| 大香蕉久久成人网| 高清毛片免费观看视频网站 | 亚洲国产欧美一区二区综合| 国产精品一区二区三区四区久久 | 国产一卡二卡三卡精品| 免费日韩欧美在线观看| 色播在线永久视频| 日韩大尺度精品在线看网址 | 亚洲欧美一区二区三区黑人| 女生性感内裤真人,穿戴方法视频| 亚洲男人天堂网一区| 亚洲视频免费观看视频| 精品一品国产午夜福利视频| 美女午夜性视频免费| 精品国产一区二区三区四区第35| 日韩精品青青久久久久久| 亚洲成人久久性| 999久久久精品免费观看国产| 亚洲色图av天堂| 日韩国内少妇激情av| 一本大道久久a久久精品| 亚洲欧美一区二区三区久久| 黄色怎么调成土黄色| 中文字幕最新亚洲高清| 校园春色视频在线观看| 亚洲色图 男人天堂 中文字幕| 电影成人av| 国产有黄有色有爽视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩中文字幕国产精品一区二区三区 | 十分钟在线观看高清视频www| 欧美中文日本在线观看视频| 大码成人一级视频| 正在播放国产对白刺激| 夜夜躁狠狠躁天天躁| 丝袜美足系列| 亚洲自偷自拍图片 自拍| 亚洲午夜精品一区,二区,三区| 亚洲av成人一区二区三| 夜夜躁狠狠躁天天躁| 黄色毛片三级朝国网站| 亚洲午夜理论影院| 亚洲午夜理论影院| 亚洲一区中文字幕在线| 亚洲欧美激情在线| 老汉色∧v一级毛片| 99国产极品粉嫩在线观看| 搡老岳熟女国产| 国产精品免费一区二区三区在线| 国产三级黄色录像| 亚洲性夜色夜夜综合| 日本a在线网址| 在线观看午夜福利视频| 在线观看免费午夜福利视频| 久久久久久亚洲精品国产蜜桃av| 日韩欧美在线二视频| 黄片播放在线免费| 91老司机精品| 麻豆一二三区av精品| 成人永久免费在线观看视频| 亚洲av五月六月丁香网| 国产精品98久久久久久宅男小说| 午夜日韩欧美国产| 欧美日韩国产mv在线观看视频| 韩国av一区二区三区四区| 久久中文字幕人妻熟女| 亚洲aⅴ乱码一区二区在线播放 | 成人亚洲精品一区在线观看| 琪琪午夜伦伦电影理论片6080| 极品教师在线免费播放| 香蕉丝袜av| 亚洲精华国产精华精| 国产精品 国内视频| 一区在线观看完整版| bbb黄色大片| 国产区一区二久久| 国产亚洲欧美精品永久| 亚洲精华国产精华精| 99国产精品一区二区蜜桃av| 久久亚洲真实| 中亚洲国语对白在线视频| 成年人黄色毛片网站| 国产麻豆69| 亚洲成国产人片在线观看| 亚洲精品国产精品久久久不卡| 大陆偷拍与自拍| 亚洲一区二区三区欧美精品| 无人区码免费观看不卡| 久久久久国产一级毛片高清牌| 久久热在线av| 无限看片的www在线观看| 精品少妇一区二区三区视频日本电影| ponron亚洲| 国产主播在线观看一区二区| 亚洲精品国产精品久久久不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产成人免费| 久久精品国产99精品国产亚洲性色 | 亚洲中文日韩欧美视频| 免费在线观看日本一区| 亚洲自偷自拍图片 自拍| 天堂俺去俺来也www色官网| 婷婷精品国产亚洲av在线| 午夜福利一区二区在线看| 欧美黑人欧美精品刺激| 99精国产麻豆久久婷婷| 亚洲五月婷婷丁香| 亚洲免费av在线视频| 人人妻人人添人人爽欧美一区卜| 超碰成人久久| 国产乱人伦免费视频| 一进一出抽搐gif免费好疼 | 久久精品成人免费网站| 女警被强在线播放| 成人亚洲精品av一区二区 | 男女下面进入的视频免费午夜 | 美女高潮到喷水免费观看| 91九色精品人成在线观看| 两个人看的免费小视频| 搡老熟女国产l中国老女人| 国产精品免费视频内射| 女人被躁到高潮嗷嗷叫费观| 日韩 欧美 亚洲 中文字幕| 日韩三级视频一区二区三区| 国产有黄有色有爽视频| 嫁个100分男人电影在线观看| 色婷婷av一区二区三区视频| 国产日韩一区二区三区精品不卡| 亚洲午夜理论影院| 亚洲成人免费电影在线观看| 亚洲美女黄片视频| 亚洲男人天堂网一区| 亚洲aⅴ乱码一区二区在线播放 | 免费在线观看完整版高清| 91成年电影在线观看| а√天堂www在线а√下载| ponron亚洲| 亚洲人成电影观看| 日韩有码中文字幕| 亚洲激情在线av| 一本大道久久a久久精品| 亚洲人成网站在线播放欧美日韩| 精品久久久久久,| 宅男免费午夜| 国产精品久久久人人做人人爽| 国产精品免费一区二区三区在线| 欧美午夜高清在线| av在线天堂中文字幕 | 国产亚洲欧美98| 国产99久久九九免费精品| 岛国在线观看网站| 久久久久久人人人人人| 亚洲精品一二三| 欧美在线一区亚洲| 亚洲精品国产区一区二| 黄网站色视频无遮挡免费观看| 美女大奶头视频| 国产精品日韩av在线免费观看 | 午夜福利,免费看| 深夜精品福利| 亚洲男人天堂网一区| 国产精品九九99| 精品国产国语对白av| 亚洲熟女毛片儿| 亚洲人成网站在线播放欧美日韩| 欧美黑人精品巨大| 不卡av一区二区三区| 国产精品久久久久成人av| 看免费av毛片| 精品卡一卡二卡四卡免费| 不卡av一区二区三区| 黄色怎么调成土黄色| 精品久久久久久成人av| 女警被强在线播放| 亚洲一区高清亚洲精品| 黑人欧美特级aaaaaa片| 黄色毛片三级朝国网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人久久性| 久久亚洲精品不卡| a级片在线免费高清观看视频| 久久久久久大精品| 色精品久久人妻99蜜桃| 国产成人影院久久av| 青草久久国产| 午夜精品在线福利| 久久亚洲真实| 老司机福利观看| 我的亚洲天堂| av国产精品久久久久影院| 一本综合久久免费| 成人国产一区最新在线观看| 亚洲国产精品sss在线观看 | 神马国产精品三级电影在线观看 | 日韩欧美免费精品| 一个人观看的视频www高清免费观看 | 欧美一区二区精品小视频在线| 91字幕亚洲| 巨乳人妻的诱惑在线观看| 少妇裸体淫交视频免费看高清 | 精品一区二区三区四区五区乱码| 欧美日韩一级在线毛片| 女警被强在线播放| 国产亚洲精品久久久久5区| 国产精品爽爽va在线观看网站 | 一级毛片精品| 久久青草综合色| 91在线观看av| 日韩中文字幕欧美一区二区| 在线永久观看黄色视频| 国产亚洲精品久久久久久毛片| 12—13女人毛片做爰片一| 久久人人精品亚洲av| av在线播放免费不卡| 黄色女人牲交| 国产成年人精品一区二区 | 69精品国产乱码久久久| 欧美日韩一级在线毛片| 午夜成年电影在线免费观看| 欧美国产精品va在线观看不卡| 亚洲一区二区三区色噜噜 | 少妇 在线观看| 亚洲精品中文字幕在线视频| 亚洲自拍偷在线| 精品人妻1区二区| 久久中文字幕人妻熟女| 老鸭窝网址在线观看| 亚洲色图 男人天堂 中文字幕| 桃色一区二区三区在线观看| 夫妻午夜视频| 日韩精品中文字幕看吧| 纯流量卡能插随身wifi吗| 一区福利在线观看| 满18在线观看网站| 18禁黄网站禁片午夜丰满| 欧美日韩亚洲国产一区二区在线观看| 精品欧美一区二区三区在线| 日韩国内少妇激情av| 色哟哟哟哟哟哟| 久久人人精品亚洲av| 男人舔女人下体高潮全视频| 中文字幕最新亚洲高清| 在线观看免费高清a一片| 欧美精品亚洲一区二区| www.自偷自拍.com| 99精品欧美一区二区三区四区| 两性夫妻黄色片| 成人亚洲精品av一区二区 | 一级毛片高清免费大全| 亚洲专区中文字幕在线| 欧美黑人精品巨大| 久久中文字幕人妻熟女| 69精品国产乱码久久久| 一级毛片精品| 激情在线观看视频在线高清| 亚洲 国产 在线| 国产精品二区激情视频| 欧美日韩福利视频一区二区| www.熟女人妻精品国产| 女人爽到高潮嗷嗷叫在线视频| 777久久人妻少妇嫩草av网站| 午夜福利一区二区在线看| 大陆偷拍与自拍| 国产精华一区二区三区| 亚洲男人的天堂狠狠| 999久久久精品免费观看国产| 久久中文看片网| 久久久久久久午夜电影 | 99riav亚洲国产免费| 国产精品亚洲av一区麻豆| 国产精品亚洲av一区麻豆| 在线播放国产精品三级| 女性被躁到高潮视频| 婷婷六月久久综合丁香| 国产av精品麻豆| 男人的好看免费观看在线视频 | 日韩精品免费视频一区二区三区| 丰满的人妻完整版| 手机成人av网站| 日韩免费av在线播放| 欧美老熟妇乱子伦牲交| 久久久久国内视频| 精品国产乱子伦一区二区三区| 大香蕉久久成人网| 精品久久蜜臀av无| 欧美久久黑人一区二区| 亚洲欧洲精品一区二区精品久久久| 精品久久久精品久久久| 91字幕亚洲| 最好的美女福利视频网| 亚洲少妇的诱惑av| 久久精品aⅴ一区二区三区四区| 成人免费观看视频高清| 中文欧美无线码| av网站在线播放免费| 国产成人免费无遮挡视频| 国产精品成人在线| 老鸭窝网址在线观看| 日本vs欧美在线观看视频| 多毛熟女@视频| 国产欧美日韩一区二区三区在线| 美女高潮喷水抽搐中文字幕| 免费女性裸体啪啪无遮挡网站| 成人三级黄色视频| 亚洲专区中文字幕在线| 黄色 视频免费看| 精品人妻在线不人妻| 少妇 在线观看| 成人av一区二区三区在线看| 免费在线观看完整版高清| 美女 人体艺术 gogo| 久久热在线av| 欧美不卡视频在线免费观看 | 国产精品久久久久久人妻精品电影| 男人舔女人下体高潮全视频| 老鸭窝网址在线观看| 中文字幕最新亚洲高清| 色老头精品视频在线观看| 日韩免费av在线播放| 午夜免费观看网址| 日本免费a在线| 美女 人体艺术 gogo| 色老头精品视频在线观看| 久久久久久久久中文| 色老头精品视频在线观看| 国产精品香港三级国产av潘金莲| 99精品久久久久人妻精品| 成熟少妇高潮喷水视频| 在线观看舔阴道视频| 在线观看午夜福利视频| 精品久久蜜臀av无| 咕卡用的链子| 黄网站色视频无遮挡免费观看| 99国产精品免费福利视频| 在线天堂中文资源库| 成人三级黄色视频| 欧美大码av| 国产精品 国内视频| 亚洲中文字幕日韩| 国产aⅴ精品一区二区三区波| 美女福利国产在线| 久久香蕉精品热| 精品人妻1区二区| 久久精品国产清高在天天线| 亚洲一卡2卡3卡4卡5卡精品中文| 别揉我奶头~嗯~啊~动态视频| 久久国产精品人妻蜜桃| 国产精品影院久久| 国产精品爽爽va在线观看网站 | 老司机亚洲免费影院| 国产精品香港三级国产av潘金莲| 黑丝袜美女国产一区| 欧美精品一区二区免费开放| 精品第一国产精品| 免费日韩欧美在线观看| 在线观看免费视频网站a站| 午夜精品在线福利| 久久久久亚洲av毛片大全| 久久精品aⅴ一区二区三区四区| 国产精品爽爽va在线观看网站 | 美女午夜性视频免费| 日韩三级视频一区二区三区| 男人操女人黄网站| 怎么达到女性高潮| 日韩 欧美 亚洲 中文字幕| 在线看a的网站| 欧美日韩亚洲国产一区二区在线观看| 中文字幕精品免费在线观看视频| 99re在线观看精品视频| av欧美777| 日本撒尿小便嘘嘘汇集6| 一级毛片女人18水好多| 伊人久久大香线蕉亚洲五| 亚洲精品美女久久av网站| 精品日产1卡2卡| 不卡av一区二区三区| 亚洲 国产 在线| 中文字幕人妻丝袜一区二区| 亚洲aⅴ乱码一区二区在线播放 | 桃色一区二区三区在线观看| 国产高清videossex| 久热爱精品视频在线9| 两个人免费观看高清视频| 亚洲人成电影观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产毛片av蜜桃av| 最近最新中文字幕大全免费视频| 国产乱人伦免费视频| 精品国产美女av久久久久小说| 久久九九热精品免费| 成人18禁高潮啪啪吃奶动态图| 国产xxxxx性猛交| 大型av网站在线播放| 成年女人毛片免费观看观看9| 在线国产一区二区在线| 嫩草影视91久久| 国产精品免费一区二区三区在线| 久久九九热精品免费| 日韩一卡2卡3卡4卡2021年| 亚洲中文字幕日韩| 露出奶头的视频| 亚洲精品久久午夜乱码| 999精品在线视频| 成人免费观看视频高清| 人妻丰满熟妇av一区二区三区| 色哟哟哟哟哟哟| 免费人成视频x8x8入口观看| 两个人看的免费小视频| 日本 av在线| 欧美人与性动交α欧美精品济南到| 侵犯人妻中文字幕一二三四区| 男女下面插进去视频免费观看| xxx96com| 国产精品国产av在线观看| 久久性视频一级片| 免费在线观看日本一区| 男人舔女人的私密视频| 制服诱惑二区| 国产成人精品久久二区二区免费| 亚洲av成人不卡在线观看播放网| 性色av乱码一区二区三区2| 色综合欧美亚洲国产小说| 女人高潮潮喷娇喘18禁视频| 亚洲精华国产精华精| 久久 成人 亚洲| 国产真人三级小视频在线观看| 国产成人av教育| 日韩国内少妇激情av| 亚洲第一青青草原| 国产精品永久免费网站| 日韩欧美国产一区二区入口| 真人一进一出gif抽搐免费| 精品卡一卡二卡四卡免费| 在线观看一区二区三区激情| 两人在一起打扑克的视频| 看免费av毛片| 日本欧美视频一区| 久久九九热精品免费| 99国产精品99久久久久| 久久狼人影院| 国产成人精品久久二区二区91| 国产熟女xx| 丁香欧美五月| 18禁裸乳无遮挡免费网站照片 | 18禁黄网站禁片午夜丰满| 午夜激情av网站| 69av精品久久久久久| 高潮久久久久久久久久久不卡| a级片在线免费高清观看视频| 桃色一区二区三区在线观看| 人人妻,人人澡人人爽秒播| 国产av又大| 亚洲精品在线美女| 男女之事视频高清在线观看| 国产伦一二天堂av在线观看| 精品电影一区二区在线| 天堂中文最新版在线下载| 91精品国产国语对白视频| 国产一区二区在线av高清观看| 国产人伦9x9x在线观看| 另类亚洲欧美激情| 成人国语在线视频| 一级a爱视频在线免费观看| 国产精品爽爽va在线观看网站 | 国产黄色免费在线视频| 热re99久久国产66热| 欧美黄色片欧美黄色片| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久av美女十八| 久久午夜亚洲精品久久|