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

    周期性分流微通道的結構設計及散熱性能*

    2021-06-01 08:32:00王晗袁禮王超王如志
    物理學報 2021年10期
    關鍵詞:微結構熱阻冷卻液

    王晗 袁禮 王超 王如志?

    1) (北京工業(yè)大學材料與制造學部,新能源材料與技術研究所,新型功能材料教育部重點實驗室,北京 100124)

    2) (北京經(jīng)濟管理職業(yè)學院信息學院,北京 100102)

    微通道散熱器在集成電路中具有重要應用, 但目前傳統(tǒng)的長直微通道散熱過程導致溫度不均勻, 散熱效率較低.本文設計了一種周期性分流微結構并與傳統(tǒng)微通道進行集成, 實現(xiàn)了一種高效率的周期性分流微通道散熱器.基于以上周期性分流微通道, 系統(tǒng)研究了單根微通道內(nèi)微結構數(shù)目、微結構的排布方式及結構參數(shù)對其散熱性能的影響.結果表明, 引入的分流微結構可增大換熱面積、打破原有層流邊界層、促進冷/熱冷卻液混合、顯著改善微通道散熱性能.在100 W/cm2的熱流密度下, 入口端冷卻液流速為1.18 m/s時, 單根微通道內(nèi)引入9組微結構后, 其最高溫度下降約24 K, 熱阻下降約44%, 努塞爾數(shù)增大約124%, 整體傳熱性能(PEC)達1.465.進一步地, 微結構采用交錯漸變的周期排布方式, 沿流動方向逐漸變寬的擾流元使得冷卻液被充分利用, 減少了高/低溫區(qū)的存在且緩解了散熱面沿流動方向存在的溫度梯度, 壓降損失相較于均勻排布也有一定程度的降低, 有效提升了散熱效率.本文提出的周期性分流微通道將在大功率集成電路及電子冷卻領域中具有廣闊的應用前景.

    1 引 言

    如今, 微電子技術發(fā)展迅猛, 電子芯片朝著小尺寸、高集成度的方向發(fā)展[1].也意味著在有限的體積內(nèi), 越來越多的電子元器件集成其中[2].在電子芯片性能提高的同時, 也產(chǎn)生了更多的熱量, 使芯片溫度升高.

    溫度對電子芯片的性能及可靠度有著直接影響.10 ℃法則表明: 在室溫附近溫度每提高10 ℃,電子元件使用壽命將減少一半; 在70—80 ℃時,溫度升高1 ℃則會導致電子器件可靠性下降5%[3].因此, 材料與器件的傳熱及散熱問題成為了微電子領域不可忽視的問題之一.對于材料的傳熱, 當下更多學者著重于對微尺度傳熱以及固液界面?zhèn)鳠犷I域的研究.對微尺度傳熱領域, Zhang等[4]總結了低維納米材料中與尺寸有關傳熱行為的最新研究, 討論了各種重要尺寸的相關熱傳輸行為, 且提出了幾種重要的熱傳遞行為和獨特的潛在物理機制.Xu等[5]總結了導電聚合物及其復合材料中熱傳遞的基本原理, 以及調(diào)節(jié)其導熱性的最新進展.對于此類問題的研究完善了熱傳遞的理論體系, 加深了人們對微尺度傳熱的理解, 但仍然存在諸多特性及理論等待人們?nèi)ヌ剿骷巴晟?對于固液界面?zhèn)鳠? Ma等[6]研究了界面電荷修飾對石墨烯和水之間的卡皮察熱阻的影響, 發(fā)現(xiàn)與不進行電荷修飾的情況相比, 卡皮察熱阻降低幅度達97%.Dmitry等[7]使用分子動力學模擬研究了多層石墨烯與水之間的卡皮察熱阻, 并通過施加整體橫斷面壓力和改變石墨烯的疏水性, 控制了微流體和納米流體系統(tǒng)中的散熱.對此類問題的研究大大推動了對跨固-液界面?zhèn)鳠峒吧釂栴}的研究及發(fā)展, 同時對于電子冷卻領域具有一定程度的借鑒意義.

    與此同時, 在電子冷卻領域, 更多學者關注于散熱方法的研究以及散熱結構的設計.于是, 各種芯片散熱技術應運而生, 主要包括自然對流、強制風冷以及微通道散熱器等[8], 其中微通道散熱器因為具有優(yōu)異的散熱性能, 自1981年被Tuckerman和Pease[9]首次提出后就備受關注.在過去的幾十年中, 有大量的研究人員對其進行研究, 致力于微通道散熱器的完善和發(fā)展.部分學者對微通道散熱器的材料進行研究, 裘騰威等[10,11]采用多孔銅作為襯底制備出的微通道散熱器散熱能力優(yōu)異.Qi等[12]制備出的金剛石微通道散熱器具有優(yōu)異的傳熱系數(shù).還有學者將不同種類[13]以及不同顆粒形狀[14]的納米流體代替水作為冷卻劑以此增強換熱.但此類方式大大增加了微通道的制造難度及材料成本且材料導熱性能發(fā)展有極限, 更多的研究著重于對微通道散熱器的結構進行改進.有學者受仿生學的啟發(fā), 對微通道的整體結構進行設計以改善其熱性能.Rubin-Jimenez等[15]及俞煒等[16]設計了y型微通道并研究了其流動機理; Ghaedamini等[17]對樹形微通道進行研究; 還有部分學者用蜂巢結構對微通道進行設計[18].除此之外, 還有學者同樣對微通道的整體結構進行設計改進, Abo-Zahhad等[19]設計出的菱形微通道、Wang等[20]設計出的肋間斷微通道以及Bhandari和Prajapati[21]設計出的混合微通道均較大程度改善了換熱性能.但由于此類結構往往較為復雜, 制備難度較大, 大部分還停留在理論及數(shù)值仿真階段.還有學者以傳統(tǒng)矩形長直微通道為基礎, 著重于微通道內(nèi)局部微結構的設計[22], Xia等[23]在微通道側(cè)壁加入周期性三角形的凹穴, 但熱性能提升較小; 有學者在長直微通道中加入周期性矩形[24]、水滴形[25]擾流元, 提升了熱性能但造成了較大壓降損失.Xie等[26]將不同形狀微針鰭引入微通道中; 文獻[27, 28]在傳統(tǒng)矩形微通道中加入垂直以及Y形分叉; Li等[29]在微通道中加入散熱片以及空腔;還有學者將設計的雙分裂圓柱體引入微通道中[30].此類方式相比于材料的更換以及微通道整體結構的設計對微通道的熱性能改善更大, 且節(jié)約了成本、制造難度低, 成為了最具前景的微通道散熱器性能提升方式.但往往由于微結構對冷卻液流動的阻礙導致壓降損失升高, 以致犧牲了較大的能耗,且微結構的均勻排布使得基底上沿流動方向出現(xiàn)明顯的溫度梯度或局部高/低溫區(qū), 溫度分布的不均勻?qū)е孪到y(tǒng)穩(wěn)定性降低.所以如何降低壓降損失以及保證溫度的均勻性成為了提升此類微通道散熱器熱性能的關鍵.

    針對微結構的引入造成壓降損失較大的問題,本研究系統(tǒng)設計與仿真了一種周期性分流微結構,發(fā)現(xiàn)其在提升散熱效率的同時, 避免微結構處通道過窄而產(chǎn)生較大壓降的損失.進一步針對微通道散熱器溫度分布不均的問題, 打破原有微通道內(nèi)微結構的均勻排布方式, 通過考慮分流微結構在微通道內(nèi)的排布方式以及結構參數(shù)的變化, 設計與仿真了一種交錯漸變分流微通道, 發(fā)現(xiàn)此排布方式在進一步提升熱性能的同時避免了換熱面上沿流動方向出現(xiàn)的溫度梯度, 減少了換熱面上的高/低溫區(qū),使換熱面溫度變得更加均勻.此設計將為高集成電路芯片散熱問題提供一種新的高效散熱結構.

    2 物理模型描述

    2.1 長直微通道散熱器

    圖1為矩形長直微通道散熱器示意圖, 將作為對照組出現(xiàn)在之后的研究中.其由40根截面為矩形的微通道組成, 等距排列于20 mm × 20 mm的基底上.在本文的計算中, 熱源的最大功率為400 W(100 W/cm2), 換熱面為其底面(下表面).假定每根微通道內(nèi)部高Hc為400 μm, 寬Wc為300 μm, 微通道上下壁厚Tc均為50 μm, 側(cè)壁厚Tw為100 μm.所以, 單根微通道的整體高度H便為500 μm, 寬度Wa為500 μm, 兩根微通道之間的壁厚為200 μm.

    圖1 長直微通道散熱器結構示意圖 (a)長直微通道散熱器; (b)單根微通道截面Fig.1.Schematic diagram of the long straight microchannel heat sink: (a) Long straight microchannel heat sink;(b) cross section of a single microchannel.

    2.2 分流微通道

    為提高微通道的散熱性能, 基于矩形長直微通道, 本文提出了一種具有分流微結構(在下文中簡稱為微結構)的微通道散熱器, 其示意圖如圖2所示.微結構由兩部分組成, 其一是微通道中間的長方體擾流元, 其二是微通道側(cè)壁的凹穴結構.由于微通道間的壁厚限制, 考慮將凹穴的深度c設置為75 μm, 為盡量減少擾流元對流體流動的阻礙,本文將擾流元的寬度b設置為75 μm.同時, 將凹穴的長度d設置為800 μm, 為保證冷卻液在通道中的流動暢通, 將擾流元的長度a設置成3/5d即480 μm.

    圖2 分流微通道結構示意圖 (a)分流微通道散熱器;(b)分流微結構局部俯視圖Fig.2.Schematic diagram of the split-flow microchannel structure: (a) Split-flow microchannel heat sink; (b) Partial top view of the split-flow microstructure.

    2.3 周期性分流微通道

    基于矩形微通道以及分流微通道, 本文設計了一種具有周期性微結構的周期性分流微通道.其結構是將分流微結構周期性排布于每根矩形微通道上(除入口端).為探究每根分流微通道上分流微結構的個數(shù)對微, 道散熱器散熱性能的影響, 將單根微通道上微結構的數(shù)目分別設置為0, 3, 9, 15個,并將其等距分布于每根微通道中, 單根微通道的結構示意圖如圖3中SM1—SM4所示.

    為改善微通道散熱面溫度分布不均的狀況, 本研究設計了四種不同的排布方式, 如圖3中DM1—DM4所示.DM1為交錯排布, 即在單根微通道內(nèi)微結構仍為等距排布, 但對于相鄰的兩根微通道,微結構交錯分布.DM2為漸密排布, 從微通道的入口端至出口端, 兩組微結構之間的間距逐漸變小,微結構越來越密集.而DM3則在等距排列的基礎上, 沿流動方向逐漸加寬了擾流元的寬度b, 在此稱之為漸變排布.DM4將DM1與DM3的排布方式進行結合, 即交錯漸變排布.將微結構數(shù)量均設置成9組.

    基于集成電路的基材大都為硅, 因此在本文中選用硅作為微通道散熱器的基底材料.考慮水為最常用的散熱流體, 本文中水被選作為冷卻介質(zhì).

    圖3 含有不同數(shù)量及排布方式微結構的單/雙根微通道示意圖: SM1 (0組); SM2 (3組); SM3 (9組); SM4 (15組);DM1 (交錯排布); DM2 (漸密排布); DM3 (漸變排布); DM4(交錯漸變排布)Fig.3.Schematic diagram of single/double microchannels with different numbers and arrangements of microstructures: SM1 (0 group); SM2 (3 groups); SM3 (9 groups);SM4 (15 groups); DM1 (staggered arrangement); DM2(gradually arranged); DM3 (gradient arrangement); DM4(staggered gradient arrangement).

    3 數(shù)學模型及計算方法

    3.1 流體流動模型

    在本研究中, 由于流體流的克努森數(shù)小于10–3, 因此假定所有微通道散熱器中的流體流都是連續(xù)的.對于連續(xù)流體, 在單位時間內(nèi)流經(jīng)單位體積空間輸入與輸出的質(zhì)量差與其內(nèi)部質(zhì)量變化的代數(shù)和等于零.同時, 假設本研究中的冷卻液為不可壓縮流體, 單位時間單位體積空間內(nèi)流入與流出的液體體積之差等于零, 所以模型的質(zhì)量守恒方程如下[31]:

    其中ρ為物質(zhì)密度;u代表流體速度;?u代表速度梯度.

    納維-斯托克斯(Navier-Stokes)方程是普遍用于描述黏性不可壓縮流體動量守恒的運動方程, 適用于此分析.同時, 計算中做出如下假設: (a)流體內(nèi)部的體積力被忽略; (b)流體為層流; (c)流體為不可壓縮且穩(wěn)定; (d)考慮黏性耗散項, 冷卻液的動力學黏度如表1所列.所以流體的動量守恒可做如下表示[31]:

    其中,μ代表動力學黏度;p為流體靜壓.

    表1 不同溫度下水的物理參數(shù)Table 1.Physical parameters of water at different temperatures.

    在本研究中, 模型的邊界條件如下:

    1)對于流場, 認為入口速度是均勻的, 溫度為室溫(293.15 K);

    其中,u0為冷卻液入口端流速;Tf代表流體溫度;T0代表入口端流體溫度.

    2)在出口處施加了一定的壓力條件(1 atm(1 atm=101325 Pa)):

    其中,Pf代表流體靜壓;Pout為流體出口端靜壓.

    在以上的邊界條件下, 通過對方程(1)以及方程(2)進行求解便可以求得三維模型中的速度場以及壓力場.

    3.2 流/固體傳熱模型

    廣義傳熱方程(能量守恒方程)可用于描述熱傳導過程.本研究涉及流體以及固體傳熱, 在流體域中, 熱量主要通過流體對流傳熱以及熱傳導進行傳播, 為簡化計算, 忽略對傳熱影響較小的黏性發(fā)熱項、壓力功項以及累積項.同時, 冷卻液的物性會隨著溫度的變化而變化, 如表1所列[32], 所以流體域中的能量守恒方程如下:

    其中ρ為物質(zhì)密度;cp為材料比熱容;T代表溫度,?T代表溫度梯度;k代表導熱系數(shù).

    在固體域中, 只考慮熱傳導項, 同時, 假設材料物性不隨時間及溫度的變化而變化, 如表2所列[32], 能量守恒方程如下:

    其邊界條件設定如下:

    1)底面施加100 W/cm2的廣義向內(nèi)熱通量邊界條件:

    其中,qm為換熱面上的熱通量.

    2)側(cè)壁使用對稱邊界條件:

    其中,qs為側(cè)壁上的熱通量.

    3)其余外表面均為熱絕緣狀態(tài)

    其中,qf為其余表面的熱通量.

    表2 硅的物理參數(shù)Table 2.Physical parameters of silicon.

    將以上邊界條件以及層流模塊求出的速度場結果代入能量守恒方程進行求解, 便可求得模型固體以及流體域中的溫度場.溫度場同時也將作為自變量對材料的物性造成影響, 在此形成雙向耦合關系.計算中考慮將流/固體傳熱模塊進行耦合, 通過穩(wěn)態(tài)方程求解在穩(wěn)定工作狀態(tài)下的壓力場、流體速度場、流/固體溫度場分布.

    本文所涉及的基材物理參數(shù)及散熱流體參數(shù)從用表1與表2參數(shù).

    3.3 流體參數(shù)

    雷諾數(shù)是一種表征流體流動情況的無量綱數(shù),在管道流中, 雷諾數(shù)小于2300的流動為層流.微通道入口端雷諾數(shù)由以下公式定義:

    其中,Re代表雷諾數(shù),D是水力直徑, 在本文中,假設流體為層流, 計算過程中, 入口端雷諾數(shù)取200—600.

    通過模型壓力場計算結果可以求得冷卻液在微通道內(nèi)的壓降損失, 壓降是流體流動前后的壓力差, 代表了流體流動過程中能量損失的大小, 由以下公式定義:

    其中, ?P代表壓降,P0為冷卻液入口端平均壓力,Pout為冷卻液出口端平均壓力.

    通過壓降損失的計算結果, 以及流體速度場的求解結果可以通過(12)式和(13)式求得系統(tǒng)所需泵送功率的大小以及達西摩擦因子的大小.其中泵送功率決定了系統(tǒng)的能耗, 摩擦因子反映了微通道對流體流動的阻力.

    泵送功率:

    其中,P代表泵送功率,Aa是通道的入口橫截面積,n0是微通道的數(shù)量,u0是流體在通道入口處的平均速度.

    摩擦因子:

    其中,f代表泵送功率,D為通道長度.

    通過流/固體溫度場的求解結果可以通過(14)式和(15)式求得模型的總熱阻以及平均努塞爾數(shù).其中熱阻反映了模型傳熱能力的大小, 努塞爾數(shù)反映了對流換熱的強烈程度.

    散熱器的總熱阻:

    其中,R代表熱阻,Tm是換熱面最高溫度,Ta是上游流體溫度,A是換熱面面積.

    微通道的平均努塞爾數(shù)[33]:

    其中,Nu為平均努塞爾數(shù),k是導熱率,Ab是每個通道的計算表面面積,Tb是通道壁表面的平均溫度,Tm代表冷卻液的平均溫度.

    在此基礎上, 采用整體傳熱性能(performance evaluation criterion, PEC)來比較周期性分流微通道與傳統(tǒng)微通道之間傳熱性能及流動阻力的提高[34].

    其中,Nu0為SM1的平均努塞爾數(shù),f0為SM1的達西摩擦因子.若PEC > 1, 則表示周期性分流微通道的整體傳熱性能優(yōu)于傳統(tǒng)長直微通道.

    3.4 網(wǎng)格設置

    本文的模擬仿真主要采用基于有限元方法進行仿真計算的COMSOL軟件進行.有限元方法的核心是通過網(wǎng)格劃分的方式將復雜的幾何實體離散成有限個規(guī)則的幾何單元.網(wǎng)格的質(zhì)量以及數(shù)量直接決定了計算結果的準確性, 使用常見的單根矩形長直微通道(SM1)進行網(wǎng)格獨立性測試, 對流體域部分采用質(zhì)量較高的六面體結構化網(wǎng)格進行劃分, 對于剩余固體域部分, 采用自由四面體網(wǎng)格進行劃分.并使用6種不同數(shù)量的網(wǎng)格來檢查網(wǎng)格獨立性, 網(wǎng)格數(shù)目分別為893083, 1128905, 1414841,1702482, 1916125, 3475672.入 口 速 度 固 定 為1.18 m/s.表3列出了以3475672個網(wǎng)格為基準的計算結果.

    如表3中所示, 隨著網(wǎng)格逐漸細化, 計算結果中壓降以及最高溫度的誤差逐漸變小, 對于網(wǎng)格5,其最高溫度和壓降的誤差很小, 可以忽略不計.因此, 兼顧計算準確性與效率, 選擇網(wǎng)格5所對應的網(wǎng)格劃分方式及大小進行后續(xù)所有仿真計算.

    4 結果與討論

    基于以上物理及數(shù)學模型與網(wǎng)格設置, 分別計算了 SM1—SM4以及DM1—DM4的壓力場、速度場以及溫度場特性.通過壓力特性、流動特性以及溫度特性三方面對設計的不同微通道散熱器進行分析, 以驗證微結構對微通道熱性能提升的有效性.同時, 探究了微結構數(shù)目對其熱性能的影響,探索其最優(yōu)微結構數(shù)目.此外, 也研究了不同的排布方式及結構參數(shù)對微通道散熱器熱性能的影響.

    表3 網(wǎng)格獨立性研究Table 3.Grid independence research.

    4.1 壓力特性

    圖4(a)給出了在入口流速為1.18 m/s的情況下, SM1~SM4 四種情況下在流體流動方向上的壓力變化曲線.其中SM1在流體進入微通道直至流出微通道的過程中, 壓力下降穩(wěn)定.可以預料到,如果流動是層流且流動充分展開, 則通常的平滑通道內(nèi)的壓降會由于通道壁面對流體的阻力顯示出線性下降, 這與圖4(a)所示的結果相接近.在引入微結構后, 流體壓力在呈現(xiàn)線性下降的基礎上出現(xiàn)了波動.圖5給出了SM2中流體流經(jīng)微結構附近時的壓力分布云圖.可以看到, 流體在途經(jīng)微結構區(qū)域前壓力分布均勻, 流體流經(jīng)微結構時首先進入突擴段, 即微通道突然變寬段, 突擴型通道中流體在流動過程中遇到管道截面突然變寬, 流體具有一定速度, 不能立即隨通道形狀的變化而變化, 從而產(chǎn)生局部紊流, 造成阻力增大.同時, 由于擾流元側(cè)壁對流體流動的阻礙作用, 壓力升高.之后流體被分為兩束流入左右兩側(cè)凹穴后繼續(xù)向前流動, 此時微通道寬度較無微結構段略寬, 由于壁面摩擦阻力使得壓力迅速下降.經(jīng)過擾流元后流體再次進入突擴段, 產(chǎn)生局部紊流, 同時受到微通道左右凹穴邊緣的阻礙導致壓力再次上升, 隨后進入突縮段,即微通道突然變窄段, 此時被打破的層流邊界層正處于恢復階段, 且由于流出微結構后流體流動方向與微通道側(cè)壁不平行, 存在一定夾角, 導致低壓區(qū)出現(xiàn), 使得流體局部壓力迅速下降.隨后, 層流邊界層恢復, 層流壓力再次呈現(xiàn)線性下降, 這與圖4(a)中SM2—SM4所示的流體沿流動方向流經(jīng)任意一組微結構時壓力出現(xiàn)的先升高, 再降低, 再升高,再降低的趨勢一致.

    由圖4(a)可知, 由于微結構的引入使得分流微通道相較于傳統(tǒng)長直微通道的壓降損失更大, 也產(chǎn)生了更多的壓力突變區(qū), 使系統(tǒng)穩(wěn)定性產(chǎn)生一定程度的降低, 這也是在所難免的.隨著微結構數(shù)量增加產(chǎn)生了更大的壓降損失, 但壓降損失的增長速度越來越慢, 這是由于流體在流經(jīng)微結構時原本層流的邊界層被打破, 經(jīng)過微結構后又逐漸恢復, 但隨著微結構數(shù)量增加, 兩組微結構之間距離變短,層流邊界層在未完全恢復時就已經(jīng)到達下一組微結構處, 減少了每組微結構所產(chǎn)生的壓力損失, 但總體來說仍然是更多組的微結構產(chǎn)生更大的壓降損失.圖4(b)給出了5種微結構數(shù)量相同但排布方式不同的微通道在不同入口端流速情況下的壓降損失.對于SM3, DM1及DM2, 雖然3種情況下壓力的波動變化有所差異, 但3種不同的微結構排布方式的總壓降損失十分接近.可以得出在微結構數(shù)目相同的情況下, 不同的微結構排布方式對微通道整體的壓降損失沒有明顯的影響.但對于DM3及DM4, 壓降損失出現(xiàn)了減少, 在入口端流速為1.88 m/s時最為明顯, 壓降減小了1200 Pa左右, 這是由于減薄了部分擾流元寬度, 降低了其對流體流動的阻礙所導致的.

    圖4 (a) SM1—SM4微通道內(nèi)主流線方向流體的壓力變化; (b) SM3及DM1—DM4在不同入口端流速下的壓降損失Fig.4.(a)Pressure change of the fluid in the direction of the main flow line in the SM1–SM4 microchannel; (b) the pressure drop loss of SM3 and DM1–DM4 at different inlet flow rates.

    圖5 SM2中微結構附近局部流體的壓力變化切面云圖Fig.5.Cross-sectional cloud diagram of pressure change of local fluid near the microstructure in SM2.

    壓降的增加將會提高微通道散熱器所需的泵送功率.泵送功率越高, 能耗越大, 微通道系統(tǒng)的穩(wěn)定性就越低.圖6(a)給出了SM1—SM44種情況下的泵送功率與總熱阻之間的關系.結果表明,整體熱阻隨著泵送功率的增加而減小.對于給定的泵送功率, 微結構的存在會導致整體熱阻下降, 且微結構數(shù)量越多, 熱阻越低, 這表明了微結構對于提升微通道熱性能的有效性.但是, 從圖6(a)可以看到, 當微結構到達9組后, 繼續(xù)增加微結構數(shù)目,相同泵送功率下的熱阻變化變得微小, 不能繼續(xù)提升微通道的散熱性能, 初步確定9組為單根微通道內(nèi)的最佳組元數(shù)目.在微結構數(shù)目相同但排布方式及結構參數(shù)不同的情況下, DM1—DM4相同泵送功率下的熱阻幾乎相同, 且均優(yōu)于SM3, 其中DM4的熱阻最低, 如圖6(b)所示.換言之, 在幾種情況熱阻相同時, 交錯漸變式排布微通道所需的泵送功率更低, 能耗損失更小, 系統(tǒng)具有更高的穩(wěn)定性.

    4.2 流動特性

    圖7所示為SM2中流體在流經(jīng)微結構前后(x–y)平面中的局部速度變化.可以看到, 微通道入口與進入微結構之前的位置流體速度相對穩(wěn)定, 流體的流動由于壁面的摩擦, 均呈現(xiàn)沿流動方向外側(cè)流速變慢而內(nèi)部流動變快的趨勢, 如圖7(a)所示.但是如圖7(b)所示, 微結構的引入改變了微通道內(nèi)的冷卻液流動, 流體在流經(jīng)微結構時內(nèi)部擾流元的存在可以有效地促進邊界層的生長, 在流體中間部分產(chǎn)生新的邊界層, 且在局部位置產(chǎn)生了較低流速的二次流及回流區(qū), 促進了冷卻液的混合, 相應地增強熱傳遞.流場在離開微結構后開始逐漸恢復到其原始狀態(tài), 如圖7(c)所示, 速度逐漸趨于穩(wěn)定且中間的邊界層逐漸消失.

    圖6 整體熱阻與泵送功率的關系 (a) SM1—SM4; (b) SM3, DM1—DM4Fig.6.Relationship between overall thermal resistance and pumping power: (a) SM1–SM4; (b) SM3, DM1–DM4.

    圖7 流體在SM2微通道內(nèi)不同位置的流速分布Fig.7.Flow velocity distribution of fluid at different positions in the SM2 microchannel.

    圖8更加直觀給出了流體流經(jīng)微結構前后邊界層生長與恢復的過程, 當流體進入微結構區(qū)域但未到達擾流元時, 其流速切面如圖8(b)所示.此時, 由于受到擾流元的干擾, 中間流體開始向左右兩側(cè)流動, 新的層流邊界層開始在中間位置產(chǎn)生.在流體經(jīng)過擾流元后, 如圖8(c)所示, 左右兩側(cè)流體開始向中間方向流動, 中部邊界層開始逐漸消失.在邊界層生長與消失的過程中, 冷卻液的橫向流動也促進了微通道內(nèi)冷/熱冷卻液的混合, 促進了熱交換.

    圖8 不同位置流速切面圖Fig.8.Cross-sectional view of flow velocity at different locations.

    4.3 溫度特性

    微通道底面的最高溫度與入口端流速之間的關系如圖9(a)所示.無論是否引入微結構, 通過改變?nèi)肟诙死鋮s液的流速均可對其表面最高溫度進行控制.冷卻液流速越大時, 散熱器的表面最高溫度越低.而在相同的流速下, 微結構的引入使微通道底面的最高溫度得以下降, 這直觀證明了微結構改善微通道散熱器熱性能的有效性.此外, 擁有更多組微結構的微通道散熱器表面最高溫度更低.在入口端流速為1.18 m/s時, SM2比SM1的最高溫度下降了13.1 K, SM3同樣比SM2最高溫度下降11.26 K, 而繼續(xù)增加微結構數(shù)量時, 最高溫度的降低變得緩慢.對于5種微結構數(shù)量相同但排布方式不同的情況, DM1~DM4的底面最高溫度均略低于均勻排布(SM3), 其中DM4最高溫度最低, 在不同流速下相比SM3均下降了3~5 K.這也證明了交錯漸變排布方式有利于進一步降低散熱器表面結溫.

    圖9 微通道底面最高溫度與入口端流速的關系 (a) SM1—SM4; (b) SM3及DM1—DM4Fig.9.Relationship between the maximum temperature on the bottom of the microchannel and the flow rate at the inlet: (a) SM1–SM4; (b) SM3 and DM1–DM4.

    圖10 (a)給出了SM1~SM4了在入口速度為1.18 m/s的情況下, 沿著主流動方向的微通道底部中心線上的溫度分布.SM1的底面溫度沿著流動方向呈現(xiàn)持續(xù)上升的趨勢, 這是由于冷卻液在流動過程中與微通道壁面進行熱交換, 外部流體被逐漸降溫, 導致?lián)Q熱效果下降所致.而引入微結構后溫度則呈現(xiàn)出周期性變化的趨勢, SM2—SM4底面溫度沿流動方向在微結構處出現(xiàn)波動, 呈現(xiàn)先下降再上升的波動, 延緩了溫度的持續(xù)上升.這主要是由于微結構的引入打破了流體的邊界層, 增大了換熱面積及促進了冷/熱冷卻液的混合所導致的.但微結構在提高微通道散熱性能的同時, 也使其溫度分布變得不均勻.如圖11所示, SM1—SM4基底出現(xiàn)多處低溫區(qū), 且換熱面均沿流動方向出現(xiàn)較大的溫度梯度, 溫度分布的不均勻?qū)a(chǎn)生熱應力, 使系統(tǒng)穩(wěn)定性降低.

    從圖11可以發(fā)現(xiàn), 不同的微結構排布方式對微通道的溫度分布有較大影響.對于DM2, 由于入口一端微結構排布稀疏, 導致其散熱面沿流動方向上溫度變化波動性較大, 散熱面溫度均勻性較低,這也是不難想象的.此難以避免的缺點也是建立對照組時并未將其與其他排布方式進行結合并進行分析的原因.但對于DM1及DM4由于采用了交錯式排布, 避免了在兩組微結構中間出現(xiàn)長距離的無微結構段, 在一定程度上避免了高/低溫區(qū)的產(chǎn)生.在圖10(b)中可以更直觀地看出DM1及DM4在主流動方向上的溫度波動變得微小.且對于DM4,其沿流動方向逐漸變寬的擾流元使得冷卻液被充分利用, 換熱面上除入口段外沿流動方向溫度始終維持于較小范圍內(nèi), 說明漸變的排布方式使微通道散熱面上的溫度變得較為均勻, 減少了高/低溫區(qū)域、減小了溫度梯度, 使系統(tǒng)變得更加穩(wěn)定.綜上所述, DM4的溫度特性最優(yōu).

    圖10 不同情況下底面上沿主流動方向上溫度變化 (a) SM1—SM4; (b) SM3及DM1—DM4Fig.10.Temperature changes along the main flow direction on the bottom surface under different conditions:(a) SM1–SM4; (b) SM3 and DM1–DM4.

    圖11 不同情況下?lián)Q熱面溫度分布云圖Fig.11.Cloud diagram of temperature distribution of heat exchange surface under different conditions.

    圖12 給出了不同情況下入口雷諾數(shù)與整體熱阻之間的關系.顯然, 不同情況下的熱阻均隨著雷諾數(shù)的增加而減小.這是由于流量的增加導致雷諾數(shù)增加, 改善了微通道的熱性能.在雷諾數(shù)相同的情況下, 引入了微結構的微通道的整體熱阻明顯低于傳統(tǒng)長直微通道.且隨著微結構數(shù)量的增加熱阻逐漸降低.對于入口端雷諾數(shù)為400的微通道, 相比于SM1, SM2熱阻下降約25%, SM3熱阻下降了約44%.這表明微結構大大提高了微通道散熱器的傳熱能力, 但在超過9組時, 微結構數(shù)量的增加對其熱阻的減小變得緩慢.同時, 在數(shù)目相同的情況下改變排布方式其熱阻也會有一定程度的降低,其中DM4熱阻最低, 相比于SM3下降了12%左右, 也證明了交錯漸變排布方式的微通道具有更優(yōu)的傳熱能力.

    圖12 微通道散熱器整體熱阻與入口雷諾數(shù)的關系Fig.12.The relationship between the overall thermal resistance of the microchannel radiator and the entrance Reynolds number.

    圖13 (a)給出了8種情況下平均努塞爾數(shù)與雷諾數(shù)的關系.在所有情況下, 努塞爾數(shù)均隨進口雷諾數(shù)的增加而增加.內(nèi)部包含微結構的微通道的努塞爾數(shù)均高于矩形長直微通道(SM1), 表明微結構的引入將增強微通道內(nèi)的對流換熱.且微結構數(shù)目為9組時努塞爾數(shù)的提升最為明顯.在改變微結構的排布方式時, 從圖13(a)可以發(fā)現(xiàn), 采取不同的排布方式努塞爾數(shù)有著不同程度的提升, 其中DM4努塞爾數(shù)最高, 相比均勻排布的SM3努塞爾數(shù)提高約5%, 故其內(nèi)部對流換熱最為劇烈.

    圖13 (a)不同情況下努塞爾數(shù)與入口端雷諾數(shù)的關系;(b)不同情況PEC與入口雷諾數(shù)的關系Fig.13.(a) The relationship between Nusselt number and inlet Reynolds number under different conditions; (b) the relationship between performance evaluation criterion and inlet Reynolds number under different conditions.

    從圖13(b)可以發(fā)現(xiàn), 在不同的入口端雷諾數(shù)下, 包含微結構的微通道散熱器的整體傳熱性能(PEC)均大于1, 證明了微結構對綜合熱性能提升的有效性.其中微結構采用交錯漸變方式排布的DM4 PEC值最高, 在微結構數(shù)量為9組、雷諾數(shù)為700時達到1.79, 展現(xiàn)出最優(yōu)的綜合熱性能.

    5 結 論

    本文基于周期性分流微通道的結構設計及其散熱仿真, 系統(tǒng)研究了微結構數(shù)量及微結構排布方式對其壓降、流動及傳熱性能的影響, 主要結論如下:

    1)分流微結構的引入可增加微通道內(nèi)部換熱面積, 打破原有層流邊界層且局部產(chǎn)生二次流動, 增強冷/熱冷卻液的混合, 降低散熱面上最高結溫, 有效提升微通道散熱器的散熱性能.

    2)微結構的數(shù)量對微通道的散熱性能有著直接的影響, 綜合考慮微通道散熱器的熱性能、壓降損失以及系統(tǒng)穩(wěn)定性, 最終確定9組為單根微通道內(nèi)最佳微結構數(shù)目.在入口流速為1.18 m/s時, 相較于傳統(tǒng)長直微通道散熱器基底表面最高結溫下降24.5 K左右, 熱阻下降44 %左右, 努塞爾數(shù)增大約124 %, 整體傳熱性能(PEC)達1.465.

    3)優(yōu)化排布方式對微通道散熱器的熱性能有進一步的提高, 其中采用交錯漸變的微結構排布方式對微通道熱性能提升最為明顯, 在微結構數(shù)量為9組, 入口端流速為1.18 m/s時, 相比于微結構均勻排布的微通道散熱器, 其壓降損失降低約500 Pa;表面最高結溫下降約3.5 K; 熱阻下降約12 %; PEC提升15 %左右.

    4)微結構采用交錯漸變的周期排布方式避免了兩組微結構間的長距離無微結構段, 沿流動方向逐漸變寬的擾流元使得冷卻液被充分利用, 從而減少微通道散熱器換熱面上的高/低溫區(qū)域, 降低沿流動方向出現(xiàn)的溫度梯度, 使散熱面上溫度分布均勻, 減少熱應力的產(chǎn)生, 一定程度上提高了冷卻效率以及系統(tǒng)穩(wěn)定性.

    本文提出的周期性分流微通道散熱結構, 為高密集電路的高效散熱問題提供了新的解決思路與方法, 有望用于超大規(guī)模集成電路芯片的高效散熱基底.

    猜你喜歡
    微結構熱阻冷卻液
    豎直單U型地埋管換熱器埋管間負熱阻現(xiàn)象的參數(shù)化研究與分析
    太陽能學報(2024年3期)2024-06-12 00:00:00
    及時換管子捂緊錢袋子
    改變缸體進排氣側(cè)冷卻液分布的方法和應用
    界面熱阻對L型鎂合金鑄件凝固過程溫度場的影響
    金屬微結構電鑄裝置設計
    用于視角偏轉(zhuǎn)的光學膜表面微結構設計
    換熱設備污垢熱阻和腐蝕監(jiān)測技術綜述
    粘結型La0.8Sr0.2MnO3/石墨復合材料的微結構與電輸運性質(zhì)
    寶馬530車冷卻液液位過低報警
    桑塔納3000車冷卻液溫度偏高
    免费看不卡的av| 中文字幕人妻丝袜一区二区| 黄片小视频在线播放| 尾随美女入室| 国产精品九九99| 国产男女超爽视频在线观看| 国产亚洲精品久久久久5区| 老鸭窝网址在线观看| 免费看十八禁软件| 欧美人与性动交α欧美精品济南到| 少妇被粗大的猛进出69影院| 18禁国产床啪视频网站| 黄色视频不卡| 国产精品av久久久久免费| 一级毛片黄色毛片免费观看视频| 欧美激情极品国产一区二区三区| 18禁国产床啪视频网站| 国产精品九九99| 你懂的网址亚洲精品在线观看| 人人妻人人澡人人爽人人夜夜| 国产伦理片在线播放av一区| 女人久久www免费人成看片| 人人妻人人澡人人爽人人夜夜| 国产有黄有色有爽视频| 纵有疾风起免费观看全集完整版| 最近中文字幕2019免费版| 国产免费福利视频在线观看| 亚洲精品日韩在线中文字幕| 观看av在线不卡| 视频在线观看一区二区三区| 国产成人精品无人区| 国产成人一区二区在线| 尾随美女入室| 人体艺术视频欧美日本| 伊人亚洲综合成人网| 男人操女人黄网站| 极品人妻少妇av视频| 欧美另类一区| 国产99久久九九免费精品| 一本久久精品| 亚洲午夜精品一区,二区,三区| 9热在线视频观看99| 美女视频免费永久观看网站| 十八禁高潮呻吟视频| 欧美黑人精品巨大| www.999成人在线观看| 午夜福利免费观看在线| 中文字幕另类日韩欧美亚洲嫩草| 男女免费视频国产| 亚洲色图 男人天堂 中文字幕| 69精品国产乱码久久久| 亚洲精品乱久久久久久| 久久久久精品国产欧美久久久 | 热99久久久久精品小说推荐| videosex国产| 男的添女的下面高潮视频| 色精品久久人妻99蜜桃| 国产黄色免费在线视频| 操美女的视频在线观看| 激情五月婷婷亚洲| 精品少妇内射三级| 中文乱码字字幕精品一区二区三区| 在线 av 中文字幕| 国产精品一二三区在线看| 自线自在国产av| 一区福利在线观看| 亚洲国产精品999| 成人亚洲欧美一区二区av| 老鸭窝网址在线观看| 日本一区二区免费在线视频| 亚洲欧美成人综合另类久久久| 亚洲av男天堂| 黑人猛操日本美女一级片| 日日摸夜夜添夜夜爱| 亚洲欧美中文字幕日韩二区| 欧美亚洲日本最大视频资源| 国产在线观看jvid| 久久久久精品人妻al黑| 天天操日日干夜夜撸| 观看av在线不卡| 最新在线观看一区二区三区 | 夫妻午夜视频| 欧美久久黑人一区二区| 9191精品国产免费久久| 我要看黄色一级片免费的| 一区二区三区激情视频| 乱人伦中国视频| 欧美在线黄色| 日本黄色日本黄色录像| 欧美日韩亚洲综合一区二区三区_| 你懂的网址亚洲精品在线观看| 母亲3免费完整高清在线观看| 十八禁人妻一区二区| 极品人妻少妇av视频| 搡老乐熟女国产| 无限看片的www在线观看| 国产视频一区二区在线看| 欧美+亚洲+日韩+国产| 美女午夜性视频免费| 久久精品国产亚洲av涩爱| 中文字幕人妻丝袜制服| 欧美精品一区二区大全| 国产精品一国产av| 国产97色在线日韩免费| 中国国产av一级| 大片免费播放器 马上看| 国产精品一二三区在线看| 亚洲第一av免费看| 少妇被粗大的猛进出69影院| 色网站视频免费| 汤姆久久久久久久影院中文字幕| 久久99精品国语久久久| 亚洲成国产人片在线观看| 大话2 男鬼变身卡| 久久久亚洲精品成人影院| 免费不卡黄色视频| 2018国产大陆天天弄谢| 久久这里只有精品19| 成人国产av品久久久| 一级黄片播放器| 免费在线观看影片大全网站 | 免费女性裸体啪啪无遮挡网站| 99国产精品一区二区蜜桃av | 国产成人啪精品午夜网站| 黑人欧美特级aaaaaa片| 午夜福利影视在线免费观看| 精品国产超薄肉色丝袜足j| 新久久久久国产一级毛片| 欧美变态另类bdsm刘玥| av有码第一页| 麻豆av在线久日| 日本vs欧美在线观看视频| 久久久久久久久免费视频了| 五月天丁香电影| 欧美成人精品欧美一级黄| 黄片小视频在线播放| 99国产精品一区二区蜜桃av | 日韩中文字幕欧美一区二区 | 国产成人av教育| 日韩一卡2卡3卡4卡2021年| 黑人欧美特级aaaaaa片| 97人妻天天添夜夜摸| 国产精品久久久av美女十八| av有码第一页| 久久九九热精品免费| 国产男人的电影天堂91| 大片免费播放器 马上看| 亚洲精品国产区一区二| 赤兔流量卡办理| 人体艺术视频欧美日本| 最黄视频免费看| av天堂久久9| 九色亚洲精品在线播放| videos熟女内射| 日韩熟女老妇一区二区性免费视频| 香蕉丝袜av| 如日韩欧美国产精品一区二区三区| 亚洲熟女精品中文字幕| 91麻豆av在线| 天天躁日日躁夜夜躁夜夜| 中文字幕制服av| 18禁国产床啪视频网站| 欧美在线一区亚洲| 黄色片一级片一级黄色片| 色综合欧美亚洲国产小说| 国产成人a∨麻豆精品| 狠狠精品人妻久久久久久综合| av欧美777| 精品久久久久久久毛片微露脸 | 亚洲精品美女久久av网站| 一级,二级,三级黄色视频| 两人在一起打扑克的视频| 久久久久久久精品精品| 久久久久国产精品人妻一区二区| 成人国产一区最新在线观看 | 日本猛色少妇xxxxx猛交久久| 真人做人爱边吃奶动态| 亚洲精品成人av观看孕妇| 国产免费又黄又爽又色| 夫妻性生交免费视频一级片| 69精品国产乱码久久久| 中国国产av一级| 亚洲精品国产区一区二| 亚洲精品久久成人aⅴ小说| 捣出白浆h1v1| 1024视频免费在线观看| 十分钟在线观看高清视频www| 丰满迷人的少妇在线观看| 午夜激情av网站| 国产日韩欧美在线精品| 免费高清在线观看日韩| 菩萨蛮人人尽说江南好唐韦庄| 久久久久国产一级毛片高清牌| 亚洲,一卡二卡三卡| 天天躁日日躁夜夜躁夜夜| 欧美av亚洲av综合av国产av| 国产在线免费精品| 国产成人影院久久av| www.精华液| www日本在线高清视频| 女人高潮潮喷娇喘18禁视频| 午夜福利免费观看在线| 日本一区二区免费在线视频| 啦啦啦中文免费视频观看日本| 国产老妇伦熟女老妇高清| 欧美在线黄色| 51午夜福利影视在线观看| 麻豆av在线久日| 欧美大码av| 亚洲成人免费电影在线观看 | 18禁黄网站禁片午夜丰满| 50天的宝宝边吃奶边哭怎么回事| 一级黄片播放器| 另类亚洲欧美激情| 免费女性裸体啪啪无遮挡网站| av福利片在线| 久久99热这里只频精品6学生| 国产精品久久久久久精品电影小说| 国产爽快片一区二区三区| 岛国毛片在线播放| 亚洲精品自拍成人| 国产男人的电影天堂91| 久久久精品94久久精品| av国产精品久久久久影院| 久久久久久免费高清国产稀缺| 免费少妇av软件| 亚洲av电影在线进入| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产中文字幕在线视频| 国产亚洲欧美精品永久| 久久精品国产综合久久久| 黑人巨大精品欧美一区二区蜜桃| 国产91精品成人一区二区三区 | 国产一级毛片在线| 深夜精品福利| 国产精品一国产av| 天堂8中文在线网| 晚上一个人看的免费电影| 日本黄色日本黄色录像| av福利片在线| 少妇猛男粗大的猛烈进出视频| 国产在视频线精品| av线在线观看网站| 久久女婷五月综合色啪小说| 久久这里只有精品19| 欧美日韩国产mv在线观看视频| 亚洲精品第二区| 亚洲精品一二三| 亚洲一区二区三区欧美精品| 亚洲精品国产色婷婷电影| www日本在线高清视频| 国产免费又黄又爽又色| av网站免费在线观看视频| www.精华液| 亚洲熟女毛片儿| 一本一本久久a久久精品综合妖精| 男女之事视频高清在线观看 | 在线 av 中文字幕| 国产一区有黄有色的免费视频| 久久人妻熟女aⅴ| 色综合欧美亚洲国产小说| 欧美 亚洲 国产 日韩一| 免费看av在线观看网站| 免费在线观看完整版高清| 午夜av观看不卡| 天堂8中文在线网| 欧美亚洲日本最大视频资源| 国产一区二区三区综合在线观看| 婷婷色麻豆天堂久久| 伦理电影免费视频| 一级a爱视频在线免费观看| 久久狼人影院| 最近最新中文字幕大全免费视频 | 精品人妻1区二区| 国产在线免费精品| 国产成人av教育| 蜜桃国产av成人99| 精品人妻熟女毛片av久久网站| 久久精品成人免费网站| 男女边吃奶边做爰视频| 亚洲欧美中文字幕日韩二区| 亚洲欧美色中文字幕在线| 久久人妻福利社区极品人妻图片 | 国产爽快片一区二区三区| 欧美激情极品国产一区二区三区| 国产色视频综合| 校园人妻丝袜中文字幕| 久久精品久久久久久久性| 捣出白浆h1v1| 少妇人妻 视频| 咕卡用的链子| 久久久国产欧美日韩av| 51午夜福利影视在线观看| 免费久久久久久久精品成人欧美视频| 亚洲精品国产av成人精品| 丁香六月天网| 高清黄色对白视频在线免费看| 大片电影免费在线观看免费| 一级a爱视频在线免费观看| 国产精品香港三级国产av潘金莲 | 99九九在线精品视频| 午夜影院在线不卡| 欧美精品一区二区大全| 免费av中文字幕在线| 欧美日韩黄片免| 国产在线视频一区二区| 性少妇av在线| 国精品久久久久久国模美| 十八禁高潮呻吟视频| 国产成人精品无人区| 亚洲av电影在线观看一区二区三区| 国产亚洲av高清不卡| 国产精品.久久久| 欧美日韩一级在线毛片| 国产成人a∨麻豆精品| 日本欧美视频一区| 亚洲精品日韩在线中文字幕| 亚洲,一卡二卡三卡| www.自偷自拍.com| 欧美激情高清一区二区三区| 亚洲精品国产一区二区精华液| 色婷婷久久久亚洲欧美| 国产日韩欧美在线精品| 一级,二级,三级黄色视频| 黄色怎么调成土黄色| 嫩草影视91久久| 大话2 男鬼变身卡| √禁漫天堂资源中文www| 黄网站色视频无遮挡免费观看| 另类亚洲欧美激情| 亚洲天堂av无毛| 精品久久久精品久久久| 亚洲人成网站在线观看播放| 亚洲色图 男人天堂 中文字幕| 国产免费视频播放在线视频| 999久久久国产精品视频| 国产精品九九99| 在线av久久热| 亚洲一区中文字幕在线| 国产主播在线观看一区二区 | 一级毛片我不卡| 欧美在线黄色| 91国产中文字幕| 脱女人内裤的视频| 国产在线一区二区三区精| 大型av网站在线播放| 精品一区二区三区av网在线观看 | 国产精品av久久久久免费| 黄色片一级片一级黄色片| 色婷婷久久久亚洲欧美| 蜜桃在线观看..| 只有这里有精品99| 啦啦啦在线免费观看视频4| 搡老岳熟女国产| 一区在线观看完整版| 一边摸一边抽搐一进一出视频| 久久av网站| 日韩av免费高清视频| 女人被躁到高潮嗷嗷叫费观| 国产成人精品无人区| 国产99久久九九免费精品| 亚洲成人免费av在线播放| 最新的欧美精品一区二区| 欧美乱码精品一区二区三区| 欧美97在线视频| 麻豆国产av国片精品| 老司机亚洲免费影院| 国产伦理片在线播放av一区| av片东京热男人的天堂| 日韩一区二区三区影片| 成人亚洲精品一区在线观看| 国产片内射在线| 午夜免费成人在线视频| 婷婷色麻豆天堂久久| 一本综合久久免费| 国产一卡二卡三卡精品| 国产欧美日韩综合在线一区二区| 美国免费a级毛片| 欧美老熟妇乱子伦牲交| 国产麻豆69| 性色av一级| 在线看a的网站| 精品一区二区三区av网在线观看 | 欧美亚洲 丝袜 人妻 在线| 美国免费a级毛片| 国产精品二区激情视频| 一边亲一边摸免费视频| 悠悠久久av| 老鸭窝网址在线观看| 精品久久久精品久久久| 亚洲第一青青草原| 狂野欧美激情性xxxx| 女人被躁到高潮嗷嗷叫费观| 纵有疾风起免费观看全集完整版| 国产精品一区二区免费欧美 | 波多野结衣av一区二区av| 国产一区二区三区av在线| 欧美精品啪啪一区二区三区 | 国产日韩欧美亚洲二区| 亚洲美女黄色视频免费看| 一级毛片我不卡| 99re6热这里在线精品视频| 亚洲精品国产一区二区精华液| 国产精品久久久久久精品古装| 老司机影院毛片| 9色porny在线观看| 搡老岳熟女国产| 女警被强在线播放| 精品福利观看| 大香蕉久久网| 久久ye,这里只有精品| 精品亚洲乱码少妇综合久久| 国产精品成人在线| 不卡av一区二区三区| 岛国毛片在线播放| 国产精品免费大片| 精品亚洲乱码少妇综合久久| 91国产中文字幕| 曰老女人黄片| 国产av精品麻豆| 啦啦啦在线免费观看视频4| 久久毛片免费看一区二区三区| 久久精品久久久久久噜噜老黄| www日本在线高清视频| 亚洲成色77777| 国产高清国产精品国产三级| 国产一区有黄有色的免费视频| 午夜av观看不卡| 免费日韩欧美在线观看| 成人亚洲精品一区在线观看| 亚洲av电影在线观看一区二区三区| 免费黄频网站在线观看国产| 日本午夜av视频| 又紧又爽又黄一区二区| 国产精品av久久久久免费| 免费一级毛片在线播放高清视频 | 日韩大码丰满熟妇| 一本综合久久免费| 国产在线免费精品| 国产高清国产精品国产三级| 国产xxxxx性猛交| 女人精品久久久久毛片| 嫩草影视91久久| 99香蕉大伊视频| 欧美成人午夜精品| 狠狠婷婷综合久久久久久88av| 欧美少妇被猛烈插入视频| 国产一区有黄有色的免费视频| 久久精品亚洲av国产电影网| 欧美日韩精品网址| 免费在线观看视频国产中文字幕亚洲 | 男的添女的下面高潮视频| 午夜91福利影院| 国产在线免费精品| 黑人欧美特级aaaaaa片| 考比视频在线观看| 国产精品人妻久久久影院| 亚洲精品国产av蜜桃| 在线观看免费日韩欧美大片| 欧美精品亚洲一区二区| 免费观看人在逋| 99久久99久久久精品蜜桃| 欧美黄色片欧美黄色片| 熟女少妇亚洲综合色aaa.| 无遮挡黄片免费观看| 亚洲av成人精品一二三区| 国产成人欧美| 一级,二级,三级黄色视频| 一边摸一边抽搐一进一出视频| av有码第一页| 亚洲欧美日韩高清在线视频 | 亚洲精品久久午夜乱码| 久久天堂一区二区三区四区| 成年女人毛片免费观看观看9 | 不卡av一区二区三区| 人人妻人人澡人人爽人人夜夜| 精品第一国产精品| av天堂在线播放| 99精品久久久久人妻精品| 少妇 在线观看| 男女下面插进去视频免费观看| 日本av手机在线免费观看| 在线看a的网站| 我要看黄色一级片免费的| 国产精品av久久久久免费| av在线播放精品| 中文字幕人妻丝袜一区二区| √禁漫天堂资源中文www| 婷婷色av中文字幕| 久久99精品国语久久久| 侵犯人妻中文字幕一二三四区| 黄色视频不卡| 校园人妻丝袜中文字幕| 人人妻人人爽人人添夜夜欢视频| 中文精品一卡2卡3卡4更新| 日本五十路高清| 两人在一起打扑克的视频| av天堂在线播放| 亚洲国产av新网站| 国产在线一区二区三区精| 亚洲欧美成人综合另类久久久| 妹子高潮喷水视频| 欧美日韩成人在线一区二区| 精品少妇一区二区三区视频日本电影| 51午夜福利影视在线观看| 无限看片的www在线观看| 最近中文字幕2019免费版| 欧美乱码精品一区二区三区| 五月开心婷婷网| 中文字幕av电影在线播放| 老汉色av国产亚洲站长工具| 黄色怎么调成土黄色| 捣出白浆h1v1| xxxhd国产人妻xxx| 一二三四社区在线视频社区8| 成人国产av品久久久| 99香蕉大伊视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲激情五月婷婷啪啪| 九色亚洲精品在线播放| 久久天躁狠狠躁夜夜2o2o | 一级毛片黄色毛片免费观看视频| 国产av国产精品国产| 99国产精品免费福利视频| 一本—道久久a久久精品蜜桃钙片| 久久精品国产a三级三级三级| 国产精品一国产av| 另类精品久久| 免费观看a级毛片全部| 丝袜在线中文字幕| 欧美日韩亚洲综合一区二区三区_| 国精品久久久久久国模美| 国产三级黄色录像| 中文字幕av电影在线播放| 大香蕉久久网| 91精品国产国语对白视频| 看十八女毛片水多多多| 少妇猛男粗大的猛烈进出视频| 19禁男女啪啪无遮挡网站| 国产成人免费无遮挡视频| 国产一区二区激情短视频 | 校园人妻丝袜中文字幕| 亚洲国产精品999| 欧美大码av| 亚洲精品美女久久av网站| 日日爽夜夜爽网站| 亚洲中文av在线| 精品国产国语对白av| 一级毛片黄色毛片免费观看视频| 日日爽夜夜爽网站| 岛国毛片在线播放| 中文字幕av电影在线播放| av有码第一页| 大片免费播放器 马上看| 色网站视频免费| 精品久久久精品久久久| 久久99一区二区三区| 国产片特级美女逼逼视频| 亚洲黑人精品在线| 国产精品久久久久久精品电影小说| 成人亚洲欧美一区二区av| 亚洲欧美一区二区三区久久| 国产精品av久久久久免费| 国产精品熟女久久久久浪| 久久免费观看电影| 黄色视频在线播放观看不卡| 男女边摸边吃奶| 精品福利观看| 久久精品成人免费网站| 妹子高潮喷水视频| 视频区图区小说| 亚洲伊人色综图| 精品欧美一区二区三区在线| 国产淫语在线视频| av网站在线播放免费| 考比视频在线观看| 黄片播放在线免费| www.av在线官网国产| 99国产精品99久久久久| 成在线人永久免费视频| 在现免费观看毛片| 久久精品成人免费网站| 啦啦啦啦在线视频资源| 操出白浆在线播放| 国产高清视频在线播放一区 | 国产色视频综合| 久久久久久久国产电影| av片东京热男人的天堂| 国产精品香港三级国产av潘金莲 | 精品一区二区三卡| 国产亚洲午夜精品一区二区久久| 欧美变态另类bdsm刘玥| 丰满迷人的少妇在线观看| 国产在线免费精品| 欧美乱码精品一区二区三区| 亚洲av日韩在线播放| 欧美精品亚洲一区二区| 在线观看免费日韩欧美大片| 中文字幕制服av| 日本wwww免费看| 少妇精品久久久久久久| 一级片'在线观看视频| h视频一区二区三区| 久久久久精品国产欧美久久久 | 中文字幕制服av| 极品人妻少妇av视频| 男人添女人高潮全过程视频| 久久久精品94久久精品| 久久鲁丝午夜福利片| 午夜福利影视在线免费观看| 一级片免费观看大全| 欧美成人午夜精品| 一区在线观看完整版| 99久久人妻综合| 精品久久蜜臀av无|