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

    基于四方程模型的LBE環(huán)管湍流換熱數(shù)值研究

    2022-04-25 01:01:32蘇興康
    原子能科學(xué)技術(shù) 2022年4期
    關(guān)鍵詞:圓管無量雷諾數(shù)

    蘇興康,顧 龍,3,*

    (1.中國科學(xué)院 近代物理研究所,甘肅 蘭州 730000;2.中國科學(xué)院大學(xué) 核科學(xué)與技術(shù)學(xué)院,北京 100049;3.蘭州大學(xué) 核科學(xué)與技術(shù)學(xué)院,甘肅 蘭州 730000)

    由于鉛及鉛鉍共晶合金(LBE)熱工水力性能良好、低壓下沸點(diǎn)高、化學(xué)惰性強(qiáng),鉛冷快堆(LFR)被選為第4代核能系統(tǒng)的主力堆型之一[1]。中國加速器驅(qū)動嬗變研究裝置(CiADS)將LBE冷卻快堆作為次臨界堆堆型[2-3]。三角形排列棒束是LFR燃料組件常見的形式,研究LBE冷卻棒束的熱工水力特性對LFR安全性和經(jīng)濟(jì)性具有重要意義。在不考慮燃料棒相互作用的情況下,三角形棒束可簡化為用六邊形傳熱單元等效的環(huán)管[4-6]。因此,環(huán)管內(nèi)的湍流換熱研究有助于分析三角形棒束熱工水力。LBE的不透明、腐蝕性特點(diǎn)增加了LBE環(huán)境中傳熱實(shí)驗(yàn)的難度[7-8],需借助計算流體力學(xué)(CFD)數(shù)值模擬LBE在環(huán)形管道內(nèi)的湍流換熱傳熱過程。然而,LBE的分子普朗特數(shù)(Pr≈0.01~0.03)較低,導(dǎo)致雷諾時均方法(RANS)中的雷諾比擬假設(shè)會影響LBE的傳熱計算[9]。即用于模擬傳統(tǒng)流體(Pr≈1)的恒定湍流普朗特數(shù)(Prt),可能會影響LBE傳熱的評估[10]。在過去的40年里,為提高LBE的計算精度,國內(nèi)外研究了用于封閉能量方程的兩方程kθ-εθ模型。用于封閉動量方程的兩方程k-ε湍流模型和用于封閉能量方程的兩方程kθ-εθ傳熱模型統(tǒng)稱為四方程k-ε-kθ-εθ模型。

    文獻(xiàn)[11-21]對四方程模型的系數(shù)函數(shù)和近壁熱湍流行為進(jìn)行了大量的研究,發(fā)展了多種適用于不同Pr流體的四方程模型。近年來,Manservisi等[22-24]在前人研究的基礎(chǔ)上,發(fā)展了新的四方程LBE湍流換熱模型,并利用直接數(shù)值模擬(DNS)數(shù)據(jù)和LBE傳熱實(shí)驗(yàn)數(shù)據(jù)對模型進(jìn)行了驗(yàn)證,在平板、管道、三角形和四邊形棒束通道內(nèi)的計算結(jié)果較好。文獻(xiàn)[25-30]在Manservisi四方程模型的基礎(chǔ)上,提出了比耗散率k-ω-kθ-ωθ和對數(shù)比耗散率k-Ω-kθ-Ωθ形式的四方程模型,并進(jìn)一步計算和驗(yàn)證了LBE在平面、管道、后臺階和六邊形棒束中的湍流換熱,預(yù)測精度良好。

    盡管人們對低Pr流體的可靠計算工具愈加感興趣,但商業(yè)代碼仍較缺乏,基于四方程模型的LBE環(huán)管內(nèi)的熱工水力研究較少。因此,本文基于開源CFD軟件OpenFOAM,開發(fā)四方程求解器4eqnFoam?;?eqnFoam,計算LBE在平板和圓管內(nèi)的湍流換熱,并與DNS數(shù)據(jù)和實(shí)驗(yàn)關(guān)系式進(jìn)行對比驗(yàn)證。最后,基于該求解器,計算LBE在環(huán)管內(nèi)的流動傳熱,獲得不同加熱條件下LBE的關(guān)鍵參量,為計算鉛鉍流體湍流換熱提供參考。

    1 湍流換熱模型

    環(huán)管幾何模型示于圖1。內(nèi)、外壁面分別受恒熱流qi、qo加熱。邊界條件劃分為入口Γi、出口Γo、外壁面Γw,o和內(nèi)壁面Γw,i。Ri和Ro分別為內(nèi)外半徑。

    圖1 同心環(huán)管示意圖Fig.1 Sketch of concentric annulus

    對于不可壓縮、常物性、無重力的充分發(fā)展雷諾時均方程組,有:

    (1)

    (2)

    (3)

    (4)

    其中,ux為沿流動方向的速度分量。

    (5)

    其中,k和νt分別為湍動能和湍流運(yùn)動黏度。采用引入了Kolmogorov速度尺度uε=(vε)1/4的Abe兩方程k-ε湍流模型計算νt[17]。ε為湍流耗散率。

    (6)

    其中,αt為湍流熱擴(kuò)散系數(shù)。選用Manservisi兩方程kθ-εθ傳熱模型計算αt[22-24]。該模型引入了湍流溫度動能kθ和湍流溫度耗散率εθ,能較好地計算LBE湍流傳熱。

    2 數(shù)值求解器

    OpenFOAM是基于C++編寫的開源CFD類庫[31-32]?;谏鲜鐾牧鱾鳠崮P秃蚈penFOAM,開發(fā)了4eqnFoam求解器。如圖2所示,4eqnFoam主要由動量方程UEqn.H、能量方程TEqn.H、SIMPLE壓力泊松方程pEqn.H、k-ε湍流模型VEqn.H、kθ-εθ換熱模型HEqn.H等子程序組成。前處理時,采用GAMBIT軟件劃分網(wǎng)格得到mesh文件,通過fluentMeshToFoam工具導(dǎo)入OpenFOAM中得到網(wǎng)格文件polyMesh。在前處理時,初始及邊界條件數(shù)據(jù)、網(wǎng)格及流體物性數(shù)據(jù)、離散及求解格式數(shù)據(jù)分別存放于0、constant、system文件夾中。求解器執(zhí)行計算時,前處理數(shù)據(jù)被求解器相關(guān)子程序調(diào)用?;谟邢摅w積法離散并求解穩(wěn)態(tài)控制方程。對流項(xiàng)采用高斯迎風(fēng)格式、擴(kuò)散項(xiàng)采用高斯線性格式。矩陣選用高斯賽德爾進(jìn)行迭代求解,殘差收斂條件為:Max|Qn+1/Qn-1|<10-10,其中Q代表ui、T、p、k、ε、kθ、εθ,n為迭代步數(shù)。

    圖2 4eqnFoam求解器在OpenFOAM中的框架Fig.2 Framework of 4eqnFoam solver in OpenFOAM

    3 計算結(jié)果

    3.1 數(shù)值驗(yàn)證(三維平板)

    圖3 平板無量綱速度、無量綱雷諾應(yīng)力、無量綱溫度和無量綱雷諾熱流Fig.3 Non-dimensional velocity,non-dimensional Reynolds stress,non-dimensional temperature and non-dimensional Reynolds heat flux of plane

    3.2 數(shù)值驗(yàn)證(圓管,Ri=0)

    無量綱速度u+和無量綱雷諾應(yīng)力τ12沿圓管壁面無量綱法向距離y+的分布如圖4所示。由圖4a可見,工況B、D與管道壁面速度分布規(guī)律吻合良好,線性區(qū)滿足u+=y+,且對數(shù)區(qū)滿足u+=2.5lny++5,相對誤差均在4%以內(nèi)。由圖4b可見,工況A無量綱雷諾應(yīng)力分布與工況DNS-A吻合良好,相對誤差在2%以內(nèi)。

    圖4 圓管無量綱速度和無量綱雷諾應(yīng)力Fig.4 Non-dimensional velocity and non-dimensional Reynolds stress of pipe

    圖5 圓管無量綱溫度和無量綱溫度脈動Fig.5 Non-dimensional temperature and non-dimensional temperature fluctuation of pipe

    圓管無量綱雷諾熱流〈uj′θ′〉+的分布如圖6a所示,工況A的雷諾熱流分布預(yù)測結(jié)果與DNS-A分布趨勢相近,整體預(yù)測結(jié)果高于DNS結(jié)果,最大相對誤差約為20%。努塞爾數(shù)Nu隨派克利數(shù)Pe的分布示于圖6b,將努塞爾數(shù)計算結(jié)果與Kirillov推薦的圓管傳熱關(guān)系式(Nu=4.5+0.018Pe0.8,25

    3.3 環(huán)管內(nèi)流動傳熱

    1)環(huán)管Ⅰ:外壁絕熱內(nèi)壁恒熱流加熱(qi=qw,qo/qi=0)

    表1 環(huán)管Ⅰ計算工況Table 1 Calculation condition of annulus Ⅰ

    圖7為環(huán)管Ⅰ無量綱速度u/ub和無量綱雷諾應(yīng)力τ12分別沿外壁離開到內(nèi)壁的無量綱距離y/d的分布。工況A的速度場、雷諾應(yīng)力場均與工況DNS-A吻合良好。最大速度出現(xiàn)在y/d=1.10處,零雷諾應(yīng)力出現(xiàn)在y/d=1.14附近。

    圖7 環(huán)管Ⅰ無量綱速度和無量綱雷諾應(yīng)力Fig.7 Non-dimensional velocity and non-dimensional Reynolds stress in annulus Ⅰ

    圖8 不同雷諾數(shù)下的環(huán)管Ⅰ無量綱溫度和無量綱溫度脈動Fig.8 Non-dimensional temperature and non-dimensional temperature fluctuation for different Re in annulus Ⅰ

    圖9a、b分別為不同雷諾數(shù)下環(huán)管Ⅰ的無量綱雷諾熱流〈uj′θ′〉+和無量綱熱擴(kuò)散系數(shù)α+=αt/α(α+用于表征湍流熱流量與分子熱流量的比值)。由圖9a可見,雷諾數(shù)越大,無量綱雷諾熱流峰值越大。由圖9b可見,在近壁區(qū),α+?1,分子導(dǎo)熱作用占主導(dǎo)地位。當(dāng)LBE的流動雷諾數(shù)較小時,全流場α+均小于1。只有當(dāng)雷諾數(shù)增加到一定程度時,如工況F,開始出現(xiàn)轉(zhuǎn)折點(diǎn)α+≥1,湍流熱擴(kuò)散作用在主流區(qū)開始強(qiáng)于LBE的分子導(dǎo)熱作用。

    環(huán)管Ⅰ中的湍流普朗特數(shù)Prt=vt/αt分布如圖9c所示。主流區(qū)的Prt分布相對近壁區(qū)平滑。隨著雷諾數(shù)的增加,平均湍流普朗特數(shù)〈Prt〉逐漸減小,當(dāng)雷諾數(shù)增加到一定程度時,如工況D、E和F,〈Prt〉變化不明顯。此外,選取了多個外壁絕熱內(nèi)壁加熱的環(huán)形管道內(nèi)液態(tài)金屬傳熱關(guān)系式對比Nu計算結(jié)果,如文獻(xiàn)[38-42]的傳熱關(guān)系式。在低Pe下4eqnFoam預(yù)測的Nu結(jié)果與仇子鋮[39]和Jaeger[40]關(guān)系式符合良好,中等Pe下與朱鋒杰[38]和石雙凱[41]關(guān)系式符合良好,而高Pe下與張貴勤[42]傳熱關(guān)系式符合良好(圖9d)。

    圖9 不同Re下環(huán)管Ⅰ的無量綱雷諾熱流、無量綱熱擴(kuò)散系數(shù)、湍流普朗特數(shù)和努塞爾數(shù)Fig.9 Non-dimensional Reynolds heat flux,non-dimensional thermal diffusivity,turbulent Prandtl number and Nusselt number for different Re in annulus Ⅰ

    2)環(huán)管Ⅱ:內(nèi)外壁面熱流密度相等(qi=qo=qw)

    令qi=qo=qw,環(huán)管Ⅱ其他計算工況與表1一致。內(nèi)外壁面熱流密度相等時,θ+的最大值和零梯度值出現(xiàn)在最大速度值對應(yīng)位置(y/d=1.10)附近(圖10a)。相比于外壁絕熱時,無量綱溫度脈動在外壁附近多出1個峰值,整體波動變化情況更加劇烈,且兩峰之間的谷值位置均偏向內(nèi)壁(圖10b)。無量綱雷諾熱流的分布曲線和雷諾應(yīng)力分布曲線均不對稱(圖10c)。零雷諾熱流的位置出現(xiàn)在y/d=1.07附近,且不隨雷諾數(shù)變化。隨著雷諾數(shù)的增加,平均湍流普朗特數(shù)〈Prt〉逐漸減小,當(dāng)雷諾數(shù)增加到一定程度時,如工況D、E和F,〈Prt〉變化不明顯(圖10d)。

    圖10 不同Re下環(huán)管Ⅱ的無量綱溫度、無量綱溫度脈動、無量綱雷諾熱流和湍流普朗特數(shù)Fig.10 Non-dimensional temperature,non-dimensional temperature fluctuation,non-dimensional Reynolds heat flux and turbulent Prandtl number for different Re in annulus Ⅱ

    3.4 四方程模型求解器性能分析及適應(yīng)性評價

    當(dāng)前,4eqnFoam求解器是基于OpenFOAM V6版本特性編寫的,編譯集成在Ubuntu 18.04系統(tǒng)環(huán)境中。串行運(yùn)行下,求解不同幾何和網(wǎng)格算例時,求解器所消耗的計算時間和內(nèi)存列于表2。測試處理器為Intel(R)Core(TM)i7-9700 CPU @ 3.00 GHz,系統(tǒng)運(yùn)行內(nèi)存為16 GB @ 2667 MHz,測試計算平臺是8核心8線程的HP ProDesk 680 G4 MT。

    表2 4eqnFoam求解器性能分析Table 2 Performance analysis of 4eqnFoam solver

    從平板、圓管和環(huán)管的計算驗(yàn)證結(jié)果來看,本研究開發(fā)的4eqnFoam求解器具備了初步計算不可壓縮、常物性和無重力條件下LBE湍流換熱過程的能力,并能獲得雷諾熱流、溫度脈動和湍流普朗特數(shù)等結(jié)果。由于驗(yàn)證過程中只考慮了Pr=0.025常物性下LBE在簡單幾何中的對流傳熱,對于更廣范圍內(nèi)的低Pr流體、復(fù)雜幾何結(jié)構(gòu)(后臺階、棒束通道等)以及變物性條件下4eqnFoam的適應(yīng)性和計算精度,還需進(jìn)一步研究。

    4 結(jié)論

    為獲得LBE在環(huán)管內(nèi)不同加熱條件下的傳熱現(xiàn)象,基于開源CFD程序OpenFOAM,開發(fā)了四方程模型求解器4eqnFoam?;?eqnFoam求解器,將LBE在平板和圓管內(nèi)的傳熱計算結(jié)果與DNS和實(shí)驗(yàn)關(guān)系式進(jìn)行了對比,驗(yàn)證了4eqnFoam的有效性。將模擬計算擴(kuò)展到LBE在內(nèi)壁加熱外壁絕熱、內(nèi)外壁面等熱流加熱的環(huán)管內(nèi),獲得了不同雷諾數(shù)下LBE充分發(fā)展的無量綱溫度、溫度脈動、雷諾熱流、熱擴(kuò)散系數(shù)和湍流普朗特數(shù)等分布。結(jié)果表明:隨著雷諾數(shù)的增加,兩種加熱條件下LBE在環(huán)管內(nèi)的無量綱溫度脈動、雷諾熱通量、湍流熱擴(kuò)散作用變大,而平均湍流普朗特數(shù)減?。豢拷诿娴耐牧髌绽侍財?shù)較湍流核心處的變化劇烈;當(dāng)雷諾數(shù)增加到一定程度時,平均湍流普朗特數(shù)變化不明顯?;谒姆匠棠P颓蠼馄鳎蔀橛嬎鉒BE熱工水力現(xiàn)象提供更多參考。

    猜你喜歡
    圓管無量雷諾數(shù)
    烏雷:無量之物
    一種方便連接的涂塑鋼管
    鋼管(2021年2期)2021-11-30 02:11:01
    劉少白
    藝術(shù)品(2020年8期)2020-10-29 02:50:02
    一種圓管內(nèi)孔自動打磨機(jī)的設(shè)計
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    論書絕句·評謝無量(1884—1964)
    炳靈寺第70 窟無量壽經(jīng)變辨識
    西藏研究(2017年3期)2017-09-05 09:45:07
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設(shè)計研究
    柔性圓管在渦激振動下的模態(tài)響應(yīng)分析
    国产成人欧美在线观看| 欧美性猛交╳xxx乱大交人| 女同久久另类99精品国产91| av中文乱码字幕在线| 最近最新免费中文字幕在线| 国产一区二区激情短视频| 国产成人av激情在线播放| 亚洲黑人精品在线| 国产一区二区在线观看日韩 | 久久久久九九精品影院| av有码第一页| 18禁黄网站禁片免费观看直播| 变态另类成人亚洲欧美熟女| 看黄色毛片网站| 欧美另类亚洲清纯唯美| 日韩成人在线观看一区二区三区| 国产蜜桃级精品一区二区三区| 在线永久观看黄色视频| 日韩三级视频一区二区三区| 国产亚洲精品久久久久5区| 国产黄色小视频在线观看| а√天堂www在线а√下载| 国产精品久久久久久人妻精品电影| 国产在线观看jvid| 欧美黄色片欧美黄色片| 亚洲人成77777在线视频| 叶爱在线成人免费视频播放| 一边摸一边抽搐一进一小说| 亚洲一卡2卡3卡4卡5卡精品中文| 日本撒尿小便嘘嘘汇集6| av福利片在线观看| 一级毛片女人18水好多| 深夜精品福利| 国产高清视频在线播放一区| 日日爽夜夜爽网站| 蜜桃久久精品国产亚洲av| 午夜福利视频1000在线观看| 国产精品九九99| 欧美成人午夜精品| 美女扒开内裤让男人捅视频| a在线观看视频网站| 国产高清激情床上av| 黄色 视频免费看| 欧美另类亚洲清纯唯美| 精品少妇一区二区三区视频日本电影| 香蕉国产在线看| 久久久久久久久久黄片| 91字幕亚洲| 国产亚洲av高清不卡| 天堂av国产一区二区熟女人妻 | 波多野结衣巨乳人妻| 日本a在线网址| 亚洲一区中文字幕在线| 一级毛片女人18水好多| 久久人妻av系列| 久久久久国产一级毛片高清牌| 1024手机看黄色片| 黄频高清免费视频| 国产精品久久久久久精品电影| 日韩精品青青久久久久久| 亚洲欧美激情综合另类| 国产av又大| 午夜福利在线观看吧| 国产精品 国内视频| 免费观看人在逋| 日韩av在线大香蕉| 非洲黑人性xxxx精品又粗又长| 国产v大片淫在线免费观看| 亚洲天堂国产精品一区在线| 九色国产91popny在线| 伦理电影免费视频| 亚洲午夜理论影院| 精品少妇一区二区三区视频日本电影| 日韩欧美 国产精品| 国产男靠女视频免费网站| 久久香蕉精品热| 久久精品人妻少妇| 国产成人系列免费观看| 琪琪午夜伦伦电影理论片6080| 精品免费久久久久久久清纯| 一级作爱视频免费观看| 欧美色欧美亚洲另类二区| 色尼玛亚洲综合影院| 50天的宝宝边吃奶边哭怎么回事| 国产精品一区二区三区四区免费观看 | 欧美三级亚洲精品| 97人妻精品一区二区三区麻豆| 国产蜜桃级精品一区二区三区| 国产三级黄色录像| 亚洲一码二码三码区别大吗| 1024手机看黄色片| 美女扒开内裤让男人捅视频| 欧美性长视频在线观看| 国产一级毛片七仙女欲春2| 日本五十路高清| av在线播放免费不卡| 草草在线视频免费看| 亚洲熟妇熟女久久| 一夜夜www| 午夜福利在线观看吧| 亚洲五月婷婷丁香| 亚洲国产中文字幕在线视频| 熟女少妇亚洲综合色aaa.| 亚洲欧美日韩高清在线视频| 在线观看免费视频日本深夜| 国产一区在线观看成人免费| 亚洲精品美女久久av网站| 亚洲aⅴ乱码一区二区在线播放 | 中文字幕高清在线视频| 国产亚洲精品久久久久久毛片| 亚洲精品在线美女| 97人妻精品一区二区三区麻豆| 长腿黑丝高跟| 18美女黄网站色大片免费观看| 99热6这里只有精品| 黄色丝袜av网址大全| 麻豆一二三区av精品| 白带黄色成豆腐渣| 国产精品影院久久| 国产一区二区激情短视频| 久久久久久久精品吃奶| 极品教师在线免费播放| 国内精品一区二区在线观看| 精品久久久久久久人妻蜜臀av| 一夜夜www| 久久久精品欧美日韩精品| 国产在线观看jvid| 亚洲国产欧美一区二区综合| 波多野结衣高清作品| 嫩草影视91久久| 国产精品一区二区免费欧美| 毛片女人毛片| 欧美高清成人免费视频www| 国产成人av教育| 精品熟女少妇八av免费久了| 蜜桃久久精品国产亚洲av| 丁香六月欧美| 亚洲专区国产一区二区| 久久国产乱子伦精品免费另类| 亚洲国产欧洲综合997久久,| 久久精品国产亚洲av高清一级| 女警被强在线播放| 日本三级黄在线观看| 欧美日本视频| 在线观看www视频免费| 国产一区二区在线观看日韩 | 夜夜爽天天搞| 少妇的丰满在线观看| 两个人看的免费小视频| 99国产精品一区二区三区| 黑人欧美特级aaaaaa片| 亚洲av日韩精品久久久久久密| 欧美极品一区二区三区四区| 日本一区二区免费在线视频| 69av精品久久久久久| 制服丝袜大香蕉在线| 欧美性猛交╳xxx乱大交人| 天堂√8在线中文| 欧美在线黄色| 婷婷亚洲欧美| 九色国产91popny在线| 国产精品 欧美亚洲| 国产成人啪精品午夜网站| 狂野欧美激情性xxxx| 99精品久久久久人妻精品| 香蕉国产在线看| 色综合站精品国产| 亚洲成人免费电影在线观看| 日韩欧美国产在线观看| 舔av片在线| 日韩欧美精品v在线| 99国产精品一区二区三区| 亚洲第一电影网av| 午夜福利视频1000在线观看| 久久婷婷人人爽人人干人人爱| 1024视频免费在线观看| 露出奶头的视频| 亚洲美女视频黄频| 成人国产一区最新在线观看| 国产久久久一区二区三区| 非洲黑人性xxxx精品又粗又长| 大型av网站在线播放| 免费搜索国产男女视频| 制服诱惑二区| 欧美黄色淫秽网站| 亚洲av电影不卡..在线观看| 精品第一国产精品| www.www免费av| 国产精品一区二区三区四区久久| 亚洲九九香蕉| 黄色视频不卡| 性色av乱码一区二区三区2| 最近最新免费中文字幕在线| 桃红色精品国产亚洲av| 精品熟女少妇八av免费久了| 精品久久久久久成人av| 窝窝影院91人妻| 最近最新中文字幕大全电影3| 啦啦啦韩国在线观看视频| 黑人操中国人逼视频| www.精华液| 一本一本综合久久| av超薄肉色丝袜交足视频| 亚洲精品国产精品久久久不卡| 黄色 视频免费看| 免费在线观看影片大全网站| 男人舔奶头视频| 国产又色又爽无遮挡免费看| 国产私拍福利视频在线观看| 99国产综合亚洲精品| a级毛片a级免费在线| 在线观看一区二区三区| 免费在线观看完整版高清| 一区二区三区高清视频在线| 久久国产精品影院| 亚洲一区二区三区色噜噜| 天天躁狠狠躁夜夜躁狠狠躁| 美女免费视频网站| 欧美日韩黄片免| 日韩欧美国产在线观看| 亚洲av成人av| 日韩精品免费视频一区二区三区| 欧美成人免费av一区二区三区| 91九色精品人成在线观看| 男人舔女人的私密视频| 69av精品久久久久久| 岛国视频午夜一区免费看| 久久这里只有精品中国| 国产一级毛片七仙女欲春2| 一区二区三区国产精品乱码| 丰满人妻一区二区三区视频av | 一本大道久久a久久精品| 巨乳人妻的诱惑在线观看| 99精品久久久久人妻精品| 国产精品久久久av美女十八| 日日爽夜夜爽网站| 这个男人来自地球电影免费观看| 精华霜和精华液先用哪个| 日日夜夜操网爽| 亚洲七黄色美女视频| 国产成人精品久久二区二区免费| a级毛片在线看网站| www.自偷自拍.com| 51午夜福利影视在线观看| 国产探花在线观看一区二区| 亚洲成a人片在线一区二区| 久久久精品欧美日韩精品| avwww免费| 欧美日韩黄片免| 无限看片的www在线观看| 亚洲天堂国产精品一区在线| 1024手机看黄色片| 日本a在线网址| 亚洲熟妇中文字幕五十中出| 午夜免费成人在线视频| 国产精品一区二区精品视频观看| 欧美最黄视频在线播放免费| www国产在线视频色| 久久性视频一级片| 听说在线观看完整版免费高清| 久久精品国产亚洲av香蕉五月| 午夜福利高清视频| 老司机深夜福利视频在线观看| 亚洲人成网站高清观看| 桃色一区二区三区在线观看| 黄片大片在线免费观看| a级毛片a级免费在线| 香蕉av资源在线| 99国产极品粉嫩在线观看| 国产精品免费视频内射| 天堂av国产一区二区熟女人妻 | 久久久久国产一级毛片高清牌| 18禁裸乳无遮挡免费网站照片| 90打野战视频偷拍视频| 亚洲成人国产一区在线观看| 亚洲精品一区av在线观看| 国产视频一区二区在线看| 香蕉久久夜色| 可以在线观看毛片的网站| 搡老熟女国产l中国老女人| 一卡2卡三卡四卡精品乱码亚洲| 午夜精品在线福利| 欧美在线一区亚洲| 欧美日韩中文字幕国产精品一区二区三区| 在线永久观看黄色视频| av中文乱码字幕在线| avwww免费| 色播亚洲综合网| 十八禁网站免费在线| 午夜福利在线在线| www.自偷自拍.com| 久99久视频精品免费| 香蕉丝袜av| 99在线人妻在线中文字幕| 日韩精品免费视频一区二区三区| 国产精品影院久久| 麻豆久久精品国产亚洲av| 久久国产精品人妻蜜桃| 欧美一级毛片孕妇| 亚洲欧美日韩高清专用| 一边摸一边抽搐一进一小说| а√天堂www在线а√下载| 九色成人免费人妻av| 美女黄网站色视频| 色噜噜av男人的天堂激情| 亚洲精品一区av在线观看| 免费在线观看影片大全网站| 精品少妇一区二区三区视频日本电影| 欧美 亚洲 国产 日韩一| 窝窝影院91人妻| 久久久精品大字幕| www.熟女人妻精品国产| 婷婷丁香在线五月| 午夜老司机福利片| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品中文字幕在线视频| 免费看美女性在线毛片视频| www.熟女人妻精品国产| 婷婷亚洲欧美| 国产精品久久久久久亚洲av鲁大| 久久精品国产清高在天天线| 亚洲九九香蕉| 午夜精品久久久久久毛片777| 亚洲第一电影网av| 老熟妇仑乱视频hdxx| 日韩三级视频一区二区三区| 色综合站精品国产| 国产精品自产拍在线观看55亚洲| 亚洲精品粉嫩美女一区| 国产欧美日韩精品亚洲av| 黄色 视频免费看| 中文字幕久久专区| 欧美日韩福利视频一区二区| 长腿黑丝高跟| 999久久久国产精品视频| 一区二区三区国产精品乱码| 国产成人av教育| 90打野战视频偷拍视频| 久久人人精品亚洲av| 老司机深夜福利视频在线观看| 天天躁夜夜躁狠狠躁躁| 12—13女人毛片做爰片一| 日韩欧美国产一区二区入口| 国产av在哪里看| 国产精品影院久久| 亚洲 国产 在线| 国产高清激情床上av| 久久久精品欧美日韩精品| 精品欧美国产一区二区三| 久久热在线av| 又粗又爽又猛毛片免费看| 91av网站免费观看| 1024视频免费在线观看| 亚洲欧美日韩高清专用| 最近在线观看免费完整版| 1024手机看黄色片| 三级毛片av免费| bbb黄色大片| 亚洲一区高清亚洲精品| 又爽又黄无遮挡网站| 成人18禁在线播放| 脱女人内裤的视频| 国产熟女xx| 午夜激情av网站| 在线观看www视频免费| 亚洲 国产 在线| 国产一区二区三区在线臀色熟女| 色综合欧美亚洲国产小说| 欧美日韩亚洲国产一区二区在线观看| 国产高清videossex| 50天的宝宝边吃奶边哭怎么回事| 中文字幕av在线有码专区| 一本综合久久免费| 久热爱精品视频在线9| 国产午夜福利久久久久久| 国产精品影院久久| av欧美777| 长腿黑丝高跟| 少妇裸体淫交视频免费看高清 | 久久精品人妻少妇| 黑人操中国人逼视频| 欧美zozozo另类| 真人一进一出gif抽搐免费| 久久久国产精品麻豆| 成人欧美大片| 中文字幕久久专区| 亚洲av电影不卡..在线观看| 99热这里只有是精品50| 99久久99久久久精品蜜桃| 亚洲最大成人中文| 亚洲中文字幕一区二区三区有码在线看 | 亚洲人成网站在线播放欧美日韩| 日韩免费av在线播放| 美女免费视频网站| 丰满人妻一区二区三区视频av | 久久久水蜜桃国产精品网| 黄色视频,在线免费观看| 午夜福利在线观看吧| 在线观看日韩欧美| 成年女人毛片免费观看观看9| 午夜成年电影在线免费观看| 午夜久久久久精精品| 亚洲午夜理论影院| 麻豆成人av在线观看| 国产精品影院久久| 亚洲激情在线av| a级毛片a级免费在线| 麻豆av在线久日| 99久久国产精品久久久| 欧美另类亚洲清纯唯美| 国产精品野战在线观看| 国产乱人伦免费视频| 欧美+亚洲+日韩+国产| 在线永久观看黄色视频| 亚洲性夜色夜夜综合| 国产精品久久视频播放| 亚洲成a人片在线一区二区| 久久人妻av系列| 老司机深夜福利视频在线观看| 亚洲一区二区三区色噜噜| www日本黄色视频网| 国产黄a三级三级三级人| 母亲3免费完整高清在线观看| 99精品久久久久人妻精品| 日韩欧美在线二视频| 国产亚洲精品av在线| 日韩中文字幕欧美一区二区| 精华霜和精华液先用哪个| 亚洲精品粉嫩美女一区| 最近最新免费中文字幕在线| 丁香六月欧美| 男女下面进入的视频免费午夜| 在线看三级毛片| 国产高清视频在线播放一区| 亚洲成av人片免费观看| 在线十欧美十亚洲十日本专区| 午夜福利免费观看在线| 欧美黄色淫秽网站| 国产av麻豆久久久久久久| 小说图片视频综合网站| 美女免费视频网站| 两个人免费观看高清视频| 听说在线观看完整版免费高清| 午夜福利视频1000在线观看| 成年女人毛片免费观看观看9| 国产精品免费视频内射| 黄色 视频免费看| 久久婷婷人人爽人人干人人爱| 97人妻精品一区二区三区麻豆| 久久人人精品亚洲av| 2021天堂中文幕一二区在线观| 国产私拍福利视频在线观看| 一个人观看的视频www高清免费观看 | 狂野欧美激情性xxxx| 国产三级黄色录像| 一级片免费观看大全| 亚洲第一欧美日韩一区二区三区| 日日摸夜夜添夜夜添小说| 又粗又爽又猛毛片免费看| 啦啦啦免费观看视频1| 99在线人妻在线中文字幕| 黄片大片在线免费观看| 成人18禁高潮啪啪吃奶动态图| 小说图片视频综合网站| 国产蜜桃级精品一区二区三区| 啪啪无遮挡十八禁网站| 国产激情久久老熟女| 亚洲一区中文字幕在线| 一a级毛片在线观看| 亚洲国产精品成人综合色| 欧美极品一区二区三区四区| 又黄又爽又免费观看的视频| 99re在线观看精品视频| 岛国在线免费视频观看| 亚洲精品国产一区二区精华液| 日韩精品中文字幕看吧| 一级a爱片免费观看的视频| 午夜精品久久久久久毛片777| 日本成人三级电影网站| 岛国在线免费视频观看| 中文字幕av在线有码专区| 国产一区在线观看成人免费| 麻豆av在线久日| 久久 成人 亚洲| 国产男靠女视频免费网站| 啪啪无遮挡十八禁网站| 免费观看精品视频网站| 午夜免费成人在线视频| 男人的好看免费观看在线视频 | 一区二区三区高清视频在线| 男女床上黄色一级片免费看| 男人的好看免费观看在线视频 | 国产午夜福利久久久久久| 丰满的人妻完整版| 18美女黄网站色大片免费观看| 色综合欧美亚洲国产小说| 亚洲欧美日韩东京热| 美女免费视频网站| 黄频高清免费视频| x7x7x7水蜜桃| 9191精品国产免费久久| 色哟哟哟哟哟哟| 亚洲精品色激情综合| 18禁裸乳无遮挡免费网站照片| 精品久久久久久久久久久久久| 国产精品一区二区精品视频观看| 国产av一区二区精品久久| 精品久久久久久久末码| 国产aⅴ精品一区二区三区波| 黑人巨大精品欧美一区二区mp4| 久久精品夜夜夜夜夜久久蜜豆 | 精品久久久久久久毛片微露脸| 午夜激情福利司机影院| 亚洲国产欧美网| 男人舔女人的私密视频| 国产亚洲欧美在线一区二区| 欧美日韩一级在线毛片| 天堂av国产一区二区熟女人妻 | 日韩免费av在线播放| 欧美性猛交黑人性爽| 99久久精品国产亚洲精品| 亚洲午夜理论影院| 亚洲成人精品中文字幕电影| 香蕉丝袜av| 亚洲精品国产一区二区精华液| 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美一区二区三区在线观看| 国产精品 国内视频| 制服人妻中文乱码| 国产亚洲精品久久久久久毛片| 欧美黑人欧美精品刺激| 亚洲人成电影免费在线| 啦啦啦观看免费观看视频高清| 在线观看免费视频日本深夜| 日本精品一区二区三区蜜桃| 真人一进一出gif抽搐免费| 精品国内亚洲2022精品成人| 男人舔奶头视频| 精品国产乱码久久久久久男人| 国产单亲对白刺激| 婷婷精品国产亚洲av| 亚洲精品久久成人aⅴ小说| 精品午夜福利视频在线观看一区| 老熟妇仑乱视频hdxx| 麻豆国产av国片精品| 无限看片的www在线观看| 日本成人三级电影网站| 国产av不卡久久| 人妻丰满熟妇av一区二区三区| 91麻豆精品激情在线观看国产| 特大巨黑吊av在线直播| 在线观看一区二区三区| 精品熟女少妇八av免费久了| 欧美三级亚洲精品| 啦啦啦免费观看视频1| 免费观看人在逋| 国产探花在线观看一区二区| 亚洲真实伦在线观看| 国产精品久久久久久人妻精品电影| 91大片在线观看| 嫩草影视91久久| cao死你这个sao货| 人妻丰满熟妇av一区二区三区| www.自偷自拍.com| 69av精品久久久久久| 九色国产91popny在线| ponron亚洲| 大型av网站在线播放| 亚洲 欧美一区二区三区| 搡老岳熟女国产| 国产精品美女特级片免费视频播放器 | 国产一区在线观看成人免费| 欧美日韩一级在线毛片| 国产精品野战在线观看| 人人妻,人人澡人人爽秒播| av超薄肉色丝袜交足视频| 色尼玛亚洲综合影院| 少妇粗大呻吟视频| 麻豆av在线久日| 欧美 亚洲 国产 日韩一| 国产一区二区激情短视频| 日本三级黄在线观看| 悠悠久久av| 99久久无色码亚洲精品果冻| 一区二区三区国产精品乱码| 在线免费观看的www视频| tocl精华| 国产精品一区二区精品视频观看| av在线天堂中文字幕| 一级毛片女人18水好多| 久久99热这里只有精品18| 国内精品久久久久久久电影| 精品少妇一区二区三区视频日本电影| 日韩欧美免费精品| 人妻丰满熟妇av一区二区三区| 日韩成人在线观看一区二区三区| 麻豆国产av国片精品| 久久精品亚洲精品国产色婷小说| 黄色视频不卡| 欧美久久黑人一区二区| 午夜福利欧美成人| 国产精品亚洲av一区麻豆| 久久精品人妻少妇| 欧美人与性动交α欧美精品济南到| 黑人巨大精品欧美一区二区mp4| 日日干狠狠操夜夜爽| 国产成人精品久久二区二区免费| 亚洲国产欧洲综合997久久,| 久久伊人香网站| 一级a爱片免费观看的视频| 国产日本99.免费观看| 大型黄色视频在线免费观看| 亚洲精品中文字幕一二三四区| 国产亚洲精品av在线| 午夜福利免费观看在线|