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

    艦船典型艙室氣流組織數(shù)值模擬

    2015-02-21 02:17:15權(quán)崇仁王洋于立慶陳乾謝軍龍
    中國(guó)艦船研究 2015年6期
    關(guān)鍵詞:艙室艦船云圖

    權(quán)崇仁,王洋,于立慶,陳乾,謝軍龍

    1海軍裝備部駐沈陽(yáng)地區(qū)軍事代表局,遼寧沈陽(yáng)110031

    2中國(guó)艦船研究設(shè)計(jì)中心,湖北武漢430064

    3大連船舶重工集團(tuán)有限公司,遼寧大連116011

    4華中科技大學(xué)建筑環(huán)境與能源工程系,湖北武漢430074

    艦船典型艙室氣流組織數(shù)值模擬

    權(quán)崇仁1,王洋2,于立慶3,陳乾4,謝軍龍4

    1海軍裝備部駐沈陽(yáng)地區(qū)軍事代表局,遼寧沈陽(yáng)110031

    2中國(guó)艦船研究設(shè)計(jì)中心,湖北武漢430064

    3大連船舶重工集團(tuán)有限公司,遼寧大連116011

    4華中科技大學(xué)建筑環(huán)境與能源工程系,湖北武漢430074

    針對(duì)艦船艙室空調(diào)系統(tǒng)的末端形式不同于民用建筑空調(diào)系統(tǒng)末端形式的特點(diǎn),提出對(duì)采用布風(fēng)器方式的艦船艙室的氣流組織形式進(jìn)行研究。采用計(jì)算流體動(dòng)力學(xué)(CFD)技術(shù),建立實(shí)船典型兩人艙室的數(shù)值試驗(yàn)?zāi)P停鶕?jù)艦船的實(shí)際邊界條件對(duì)該艙室夏季設(shè)計(jì)工況和冬季設(shè)計(jì)工況下的氣流組織進(jìn)行數(shù)值模擬計(jì)算,并對(duì)典型截面的速度場(chǎng)及溫度場(chǎng)進(jìn)行分析。模擬結(jié)果表明:在夏季設(shè)計(jì)工況下,雖然布風(fēng)器周圍風(fēng)速較高、溫度較低,但艙內(nèi)人員活動(dòng)區(qū)域速度場(chǎng)分布較均勻,艙內(nèi)風(fēng)速小于0.3 m/s,溫度場(chǎng)分布也較為均勻,溫度約為26~27℃;在冬季設(shè)計(jì)工況下,除布風(fēng)器周圍風(fēng)速較高、溫度較高外,艙室大部分區(qū)域風(fēng)速較低,小于0.2 m/s,溫度約為20℃。無(wú)論是夏季還是冬季設(shè)計(jì)工況,采用布風(fēng)器末端形式的典型兩人艙室人員活動(dòng)區(qū)域內(nèi)的氣流速度及溫度均滿足舒適性標(biāo)準(zhǔn)要求,結(jié)果驗(yàn)證了該典型艙室空調(diào)系統(tǒng)布置的合理性。

    艦船艙室;布風(fēng)器;數(shù)值模擬;氣流組織;計(jì)算流體動(dòng)力學(xué)

    0 引 言

    與民用船舶相比,軍用艦船艙室相對(duì)封閉、空間狹小,且人員眾多、設(shè)備復(fù)雜,極易造成空氣污染,危害艦員身心健康,影響艦船戰(zhàn)斗力的發(fā)揮[1]。艙室內(nèi)的空氣環(huán)境包括熱濕環(huán)境和空氣品質(zhì)(IAQ),合理的通風(fēng)組織是良好室內(nèi)熱濕環(huán)境和空氣品質(zhì)的根本保障[2]。因此,針對(duì)艦船艙室特點(diǎn),改進(jìn)艦船空調(diào)系統(tǒng)設(shè)計(jì)和優(yōu)化系統(tǒng)運(yùn)行策略,提高艦船艙室空氣品質(zhì),對(duì)艦船空調(diào)設(shè)計(jì)有著重要的實(shí)際意義[3]。

    目前,世界造船業(yè)已經(jīng)逐步進(jìn)入更高層次的“數(shù)字化”階段,以數(shù)字化建模仿真與優(yōu)化為特征,將信息技術(shù)全面應(yīng)用于船舶的產(chǎn)品開(kāi)發(fā)、設(shè)計(jì)、制造等全過(guò)程[4]。計(jì)算流體動(dòng)力學(xué)(CFD)技術(shù)已被廣泛應(yīng)用于氣流組織模擬及艙內(nèi)環(huán)境的研究。Liu[5]采用Airpak軟件模擬船舶艙室在不同送、回風(fēng)形式下的艙內(nèi)環(huán)境,得出下送風(fēng)的空調(diào)送風(fēng)形式有較好的效果。但在該文獻(xiàn)中并未明確送、回風(fēng)口的型式。梁彥超[6]對(duì)機(jī)艙熱環(huán)境進(jìn)行了數(shù)值模擬分析,并針對(duì)艙內(nèi)的高溫區(qū)域提出了通風(fēng)系統(tǒng)的優(yōu)化措施。周俊男[7]采用CFD軟件對(duì)士兵居住艙室進(jìn)行了3種不同氣流組織形式的數(shù)值模擬,并得出適合該艙室的氣流組織形式。但在該研究中采用的送風(fēng)口為散流器方式,而在實(shí)際的艙室中一般采用布風(fēng)器方式。

    布風(fēng)器和散流器的對(duì)比如圖1所示。散流器被廣泛應(yīng)用于民用建筑的空調(diào)系統(tǒng),而布風(fēng)器則被廣泛應(yīng)用于船舶空調(diào)系統(tǒng)。兩者雖然都是送風(fēng)末端裝置,但在送風(fēng)形式上有所區(qū)別,布風(fēng)器的送風(fēng)氣流在下底板的阻礙作用下,呈水平方向沿四周散射出去;散流器的送風(fēng)氣流在導(dǎo)流板的作用下,呈斜向下方向沿四周送風(fēng)。船舶空調(diào)系統(tǒng)常用布風(fēng)器主要是由于船舶艙室空間一般比較低矮。

    圖1 送風(fēng)末端對(duì)比Fig.1 Comparison of air terminals

    本文將對(duì)采用布風(fēng)器的艦船典型兩人艙室氣流組織進(jìn)行仿真模擬,并分析該艙室的速度場(chǎng)分布和溫度場(chǎng)分布,用來(lái)證明采用布風(fēng)器方式的典型艙室大部分區(qū)域速度場(chǎng)及溫度場(chǎng)分布較為均勻,能夠滿足我國(guó)舒適性空調(diào)調(diào)節(jié)室內(nèi)設(shè)計(jì)標(biāo)準(zhǔn)。

    1 艦船艙室描述

    艦船艙室平面圖如圖2所示,三維立體圖如圖3所示。艙室長(zhǎng)4.8 m,寬4.8 m,高2.5 m,總體積為57.6 m3。艙室設(shè)置2個(gè)布風(fēng)器送風(fēng),每個(gè)布風(fēng)器的送風(fēng)量為250 m3/h,總送風(fēng)量為500 m3/h,換氣次數(shù)為8.7次。艙室門上設(shè)有1個(gè)矩形回風(fēng)口,尺寸為500 mm×400 mm,下邊緣距下甲板高度為0.15 m。

    圖2 艙室俯視平面示意圖Fig.2 Schematic planform of the cabin

    圖3 艙室三維圖Fig.3 Three-dimensional model of the cabin

    布風(fēng)器如圖4所示(尺寸單位:mm):上端為靜壓箱,靜壓箱尺寸為0.45 m×0.45 m×0.16 m;靜壓箱下方為圓形風(fēng)管,風(fēng)管尺寸為直徑0.16 m,高0.06 m;出風(fēng)口處設(shè)置2塊擋板,上層擋板為方形,尺寸為0.26 m×0.26 m,下層擋板為圓形,直徑為0.24 m,上、下?lián)醢彘g的豎向距離為0.02 m。送風(fēng)經(jīng)過(guò)靜壓箱,再經(jīng)垂直風(fēng)口,吹到水平擋板后形成水平圓周出流。布風(fēng)器送風(fēng)口距離地板(即艙室底部)的高度為1.92 m。

    圖4 布風(fēng)器尺寸圖Fig.4 Size of the air distributor

    2 艙室數(shù)值模擬模型

    在艙室數(shù)值模擬過(guò)程中,需要保證流體流動(dòng)都符合質(zhì)量守恒、動(dòng)量守恒和能量守恒定律,故一般將連續(xù)性方程、動(dòng)量方程、能量方程稱為控制方程。

    連續(xù)性方程、動(dòng)量方程(僅給出x方向)及能量方程表達(dá)式為:

    式中:u,v,w分別為流體在x,y,z坐標(biāo)方向上的速度分量;ρ為流體密度;p為流體微元上的壓力;Su為x方向上動(dòng)量守恒方程的廣義源項(xiàng);?(ρE)/?t為能量非穩(wěn)態(tài)項(xiàng);div·(ρEV)為對(duì)流項(xiàng);ρf·V為體積力做功;-div·(pV)為壓力做功;div·(τ·V)為粘性力做功;div·(kdivT)為導(dǎo)熱項(xiàng);SE為源項(xiàng)。

    本研究根據(jù)模擬艙室的布置、送排風(fēng)口具體設(shè)置和熱濕負(fù)荷,采用GAMBIT進(jìn)行典型艙室的流體域數(shù)值模型幾何建模,采用FLUENT對(duì)該模型進(jìn)行數(shù)值求解。采用有限體積法,將待解的微分方程對(duì)每個(gè)控制體積積分并得到關(guān)于各節(jié)點(diǎn)上溫度(或者速度、熱流等其他因變量)的離散方程,求解離散方程組就可以得到相關(guān)的量[8]。采用Tecplot進(jìn)行后處理。

    3 計(jì)算模型假設(shè)及邊界條件

    建立典型艙室的流體域數(shù)值試驗(yàn)?zāi)P?,以布風(fēng)器為中心向外圍逐漸劃分不同區(qū)域并進(jìn)行網(wǎng)格加密。本艙室共劃分非結(jié)構(gòu)化網(wǎng)格60萬(wàn)個(gè),對(duì)布風(fēng)器周圍劃分6個(gè)不同的加密區(qū)域,由布風(fēng)器出口向周圍網(wǎng)格逐漸稀疏。布風(fēng)器出口處劃分10層網(wǎng)格,網(wǎng)格尺寸最小為0.2 mm;模型中不需要進(jìn)行加密的區(qū)域,即艙內(nèi)大部分區(qū)域,最大網(wǎng)格尺寸為200 mm。圖5為網(wǎng)格劃分示意圖(60萬(wàn)網(wǎng)格),圖6為艙室布風(fēng)器附近區(qū)域網(wǎng)格劃分圖。

    圖5 艙室網(wǎng)格劃分圖Fig.5 The schematic mesh plan of cabin

    圖6 艙室布風(fēng)器附近網(wǎng)格劃分圖Fig.6 The zoom of meshing near the air distributor

    艙室空間的空氣流動(dòng)是自然對(duì)流和強(qiáng)迫對(duì)流共同作用形成的混合湍流流動(dòng)。由于實(shí)際艙室內(nèi)的設(shè)備布置、空氣流動(dòng)和傳熱非常復(fù)雜,因此需要對(duì)艙室內(nèi)的空氣流動(dòng)及模型做相關(guān)合理的假設(shè)以便進(jìn)行求解:

    1)艙室內(nèi)的流動(dòng)為低速流動(dòng),可視為不可壓縮流體;

    2)艙室內(nèi)的流體屬于牛頓流體,表面應(yīng)力滿足廣義牛頓粘性應(yīng)力公式;

    3)流場(chǎng)具有高湍流雷諾數(shù),流體的湍流粘性系數(shù)具有各向同性;

    4)由于氣體流動(dòng)為低速不可壓縮流動(dòng),可忽略由流體粘性力做功所引起的耗散熱;

    5)考慮艙室流場(chǎng)為穩(wěn)態(tài),可忽略時(shí)間項(xiàng)的影響;

    6)假定艙內(nèi)設(shè)備工況穩(wěn)定,向周圍空氣放出的熱量通過(guò)定義的邊界均勻或局部地進(jìn)入流體區(qū)域,將流體與固體接觸面上的換熱轉(zhuǎn)化為純流體的加熱或?qū)α鲹Q熱,不考慮輻射傳熱;

    7)艙室為密閉空間,也即在送風(fēng)口處流入計(jì)算區(qū)域,在回風(fēng)口處流出計(jì)算區(qū)域,確保空氣不會(huì)從別處流進(jìn)或流出計(jì)算區(qū)域。

    艙室模擬分夏季工況與冬季工況。室內(nèi)設(shè)計(jì)參數(shù)及送風(fēng)參數(shù)等如表1所示。室內(nèi)熱擾包括人員和照明,將室內(nèi)熱擾折合成單位面積熱流指標(biāo)添加到地板上,取典型艙室人員為2人,人員負(fù)荷即單人全熱負(fù)荷為130 W[9],艙室內(nèi)照明面積指標(biāo)為11 W/m2,室內(nèi)總熱擾513 W,地板單位熱流指標(biāo)為22.3 W/m2。其余側(cè)壁與相鄰空調(diào)區(qū)域沒(méi)有熱量交換,設(shè)置為絕熱壁面。單個(gè)布風(fēng)器的送風(fēng)量250 m3/h。布風(fēng)器送風(fēng)管截面設(shè)置為速度進(jìn)口邊界,速度為3.46 m/s。送風(fēng)溫差為9℃。艙室通風(fēng)口設(shè)置為壓力出口邊界。模擬計(jì)算模型選擇標(biāo)準(zhǔn)k-ε模型,考慮重力的影響。模擬收斂條件設(shè)置為能量殘差10-6,各向速度及k-ε殘差為10-3。

    表1 模擬計(jì)算參數(shù)Tab.1 Parameters of simulated calculation

    本文進(jìn)一步分析了網(wǎng)格劃分的疏密對(duì)模擬結(jié)果的影響。表2為冬季工況下不同網(wǎng)格數(shù)量(即40萬(wàn)、60萬(wàn)和150萬(wàn))對(duì)計(jì)算結(jié)果的影響。網(wǎng)格總量為60萬(wàn)和150萬(wàn)的2次模擬結(jié)果較為接近,網(wǎng)格總量為40萬(wàn)的模擬結(jié)果與另外2個(gè)模擬(即60萬(wàn)與150萬(wàn)網(wǎng)格)結(jié)果有所不同,其原因是網(wǎng)格總量為40萬(wàn)時(shí),布風(fēng)器周圍網(wǎng)格尺寸較大,不如60萬(wàn)和150萬(wàn)網(wǎng)格時(shí)網(wǎng)格緊密??紤]到模擬精確性及計(jì)算機(jī)性能與效率,本文采用了60萬(wàn)網(wǎng)格進(jìn)行夏季工況及冬季工況的詳細(xì)模擬。

    表2 網(wǎng)格疏密對(duì)模擬結(jié)果的影響Tab.2 Effects of meshing on simulation results

    4 模擬計(jì)算及結(jié)果分析

    分夏季和冬季2個(gè)設(shè)計(jì)工況對(duì)該模型進(jìn)行數(shù)值模擬,并將數(shù)據(jù)進(jìn)行后處理。對(duì)典型截面的速度場(chǎng)和溫度場(chǎng)分布進(jìn)行分析??紤]的典型截面為:截面1為通過(guò)典型艙室內(nèi)2個(gè)布風(fēng)器中心的垂直截面;截面2為0.6 m高處平面截面(坐姿時(shí)膝蓋高度約為0.6 m);截面3為1.4 m高處平面截面(坐姿時(shí)頭部高度約為1.4 m)。

    4.1 夏季設(shè)計(jì)工況

    夏季工況下,各截面艙室速度場(chǎng)分布云圖分別如圖7~圖9所示。由截面1艙室速度場(chǎng)分布云圖(圖7)可以看出,送風(fēng)經(jīng)垂直風(fēng)口吹到水平擋板后形成圓周出流,出流方向?yàn)樾鄙戏?。垂直方向上,布風(fēng)器周圍風(fēng)速較高,艙室下部氣流速度較低且分布均勻,風(fēng)速在0.1~0.2 m/s之間。截面2艙室速度場(chǎng)分布云圖(圖8)顯示,在人員坐姿時(shí)膝蓋高度所在的水平面上,艙室大部分區(qū)域速度場(chǎng)分布很均勻,風(fēng)速在0.1 m/s左右,只有靠近出風(fēng)口處的風(fēng)速較高,大于0.3 m/s。截面3艙室速度場(chǎng)分布云圖(圖9)表明,在人員坐姿時(shí)頭部高度所在的水平面上,除壁面夾角附近有小部分區(qū)域速度較大,約為0.4 m/s外,艙室絕大部分區(qū)域風(fēng)速在0.1 m/s左右。

    圖7 截面1速度場(chǎng)分布(夏季工況)Fig.7 The distribution of velocity field for section 1(summer)

    圖8 截面2速度場(chǎng)分布(夏季工況)Fig.8 The distribution of velocity field for section 2(summer)

    圖9 截面3速度場(chǎng)分布(夏季工況)Fig.9 The distribution of velocity field for section 3(summer)

    不同截面艙室溫度場(chǎng)分布云圖如圖10~圖12所示。由截面1艙室溫度場(chǎng)分布云圖(圖10)可以看出,由于布風(fēng)器送風(fēng)溫度為18℃,導(dǎo)致布風(fēng)器周圍溫度較低,垂直方向上,艙室上部空氣溫度低于下部,上部空氣溫度在22~24℃之間,下部空氣溫度在26~27℃左右。截面2艙室溫度場(chǎng)分布云圖(圖11)顯示,在人員坐姿時(shí)膝蓋高度所在的水平面上,艙室大部分區(qū)域溫度場(chǎng)分布比較均勻,空氣溫度在26~27℃之間。截面3艙室溫度場(chǎng)分布云圖(圖12)表明,在人員坐姿時(shí)頭部高度所在的水平面上,艙室溫度場(chǎng)分布也較均勻,空氣溫度在25~27℃之間。

    圖10 截面1溫度場(chǎng)分布(夏季工況)Fig.10 The distribution of temperature field for section 1(summer)

    圖11 截面2溫度場(chǎng)分布(夏季工況)Fig.11 The distribution of temperature field for section 2(summer)

    圖12 截面3溫度場(chǎng)分布(夏季工況)Fig.12 The distribution of temperature field for section 3(summer)

    通過(guò)不同典型截面艙室速度場(chǎng)分布云圖及溫度場(chǎng)分布云圖可以看出,夏季設(shè)計(jì)工況下,艙室大部分區(qū)域的速度場(chǎng)分布較為均勻,速度較低,滿足我國(guó)舒適性空調(diào)調(diào)節(jié)室內(nèi)設(shè)計(jì)標(biāo)準(zhǔn)[10]。夏季室內(nèi)風(fēng)速不大于0.3 m/s的要求;艙室大部分區(qū)域溫度場(chǎng)的分布較均勻,溫度約為26~27℃,滿足我國(guó)舒適性空調(diào)調(diào)節(jié)夏季室內(nèi)設(shè)計(jì)標(biāo)準(zhǔn)范圍22~28℃[10]。布風(fēng)器周圍及艙室回風(fēng)口處風(fēng)速較高,約為0.7 m/s;布風(fēng)器周圍溫度較低,在19~21℃之間。

    4.2 冬季設(shè)計(jì)工況

    冬季工況下,各截面艙室速度場(chǎng)分布云圖如圖13~圖15所示。該工況下不同典型截面艙室速度場(chǎng)分布與夏季工況下的艙室速度場(chǎng)分布相似。艙室大部分區(qū)域速度場(chǎng)分布較為均勻,速度較低,滿足我國(guó)舒適性空調(diào)調(diào)節(jié)室內(nèi)設(shè)計(jì)標(biāo)準(zhǔn)冬季室內(nèi)風(fēng)速不大于0.2 m/s的要求。布風(fēng)器周圍及艙室回風(fēng)口處風(fēng)速較高,約為0.7 m/s。另外,局部區(qū)域(如墻角處)風(fēng)速偏高。

    圖13 截面1速度場(chǎng)分布(冬季工況)Fig.13 The distribution of velocity field for section 1(winter)

    圖14 截面2速度場(chǎng)分布(冬季工況)Fig.14 The distribution of velocity field for section 2(winter)

    圖15 截面3速度場(chǎng)分布(冬季工況)Fig.15 The distribution of velocity field for section 3(winter)

    各截面艙室溫度場(chǎng)分布云圖如圖16~圖18所示。由截面1艙室溫度場(chǎng)分布云圖(圖16)可以看出,由于布風(fēng)器送風(fēng)溫度為27℃,導(dǎo)致布風(fēng)器周圍溫度較高,再加上熱空氣上升的緣故,垂直方向上,艙室上部空氣溫度高于下部,上部空氣溫度在21~22℃之間,下部空氣溫度在在19~20℃之間。截面2艙室溫度場(chǎng)分布云圖(圖17)顯示,在人員坐姿時(shí)膝蓋高度所在的水平面上,艙室大部分區(qū)域的溫度場(chǎng)分布比較均勻,空氣溫度在19~20℃之間,而截面圖下方壁面為暴露側(cè)壁,外側(cè)溫度僅為-18℃,導(dǎo)致暴露側(cè)壁附近空氣溫度低于其他區(qū)域空氣溫度,低于18℃。截面3艙室水平溫度場(chǎng)分布云圖(圖18)表明,在人員坐姿時(shí)頭部高度所在的水平面上,艙室大部分區(qū)域溫度場(chǎng)分布也比較均勻,在19~20℃之間。

    通過(guò)不同典型截面艙室溫度場(chǎng)分布云圖可以看出,冬季設(shè)計(jì)工況下,艙室大部分區(qū)域溫度場(chǎng)分布非常均勻,溫度約為20℃,滿足我國(guó)舒適性空調(diào)調(diào)節(jié)冬季室內(nèi)設(shè)計(jì)標(biāo)準(zhǔn)范圍18~24℃[10]。布風(fēng)器周圍溫度較高,在24~26℃之間。艙室下部暴露側(cè)壁附近的空氣溫度低于18℃。

    圖16 截面1溫度場(chǎng)分布(冬季工況)Fig.16 The distribution of temperature field of section 1(winter)

    圖17 截面2溫度場(chǎng)分布(冬季工況)Fig.17 The distribution of temperature field of section 2(winter)

    圖18 截面3溫度場(chǎng)分布(冬季工況)Fig.18 The distribution of temperature field for section 3(winter)

    5 結(jié) 語(yǔ)

    本文建立了采用布風(fēng)器進(jìn)行空調(diào)送風(fēng)的艦船典型兩人艙室的流體域數(shù)值試驗(yàn)?zāi)P停M(jìn)行了數(shù)值模擬計(jì)算。通過(guò)對(duì)典型截面的速度場(chǎng)及溫度場(chǎng)的分析,結(jié)果表明,無(wú)論是在夏季設(shè)計(jì)工況還是冬季設(shè)計(jì)工況,艙內(nèi)人員工作區(qū)的艙室大部分區(qū)域速度場(chǎng)分布較為均勻,速度較低,小于0.3 m/s,滿足我國(guó)舒適性空調(diào)標(biāo)準(zhǔn)給出的速度要求。雖然布風(fēng)器出口附近速度較高,達(dá)到了0.7 m/s,但位于艙內(nèi)人員活動(dòng)區(qū)以上。在夏季設(shè)計(jì)工況下,艙室大部分區(qū)域溫度場(chǎng)分布較為均勻,溫度約為26~27℃;在冬季設(shè)計(jì)工況下,艙室大部分區(qū)域溫度場(chǎng)分布也很均勻,溫度約為20℃,均滿足舒適性空調(diào)標(biāo)準(zhǔn)。通過(guò)數(shù)值模擬計(jì)算,驗(yàn)證了典型艙室空調(diào)系統(tǒng)布置的合理性。

    [1] 艾偉,謝田華,高占勝,等.艦船艙室空氣質(zhì)量評(píng)價(jià)的污染危害普適指數(shù)法研究[C]//航海技術(shù)理論研究論文集.北京:中國(guó)航海學(xué)會(huì),2009:163-167.

    [2] 楊麗,WANG Bing.不同通風(fēng)方式與室內(nèi)空氣環(huán)境質(zhì)量的數(shù)值模擬分析[J].建筑科學(xué),2014,30(4):78-83. YANG Li,WANG Bing.Numerical simulation analy?sis of different modes of ventilation and indoor air qual?ity[J].Building Science,2014,30(4):78-83.

    [3] 周愛(ài)民,余濤,沈旭東.船舶艙室污染物傳播研究進(jìn)展[J].艦船科學(xué)技術(shù),2014,36(1):10-15. ZHOU Aiming,YU Tao,SHEN Xudong.Research ad?vances on air contaminant transport in ship cabin[J]. Ship Science and Technology,2014,36(1):10-15.

    [4] 柳存根,高鵬.數(shù)字化造船與現(xiàn)代造船模式[J].微型機(jī)與應(yīng)用,2007(9):64,69.

    [5] LIU H M.Simulation and optimization of indoor ther?mal environment in a ship air-conditioning system[J]. Procedia Environmental Sciences, 2011, 11:1055-1063.

    [6] 梁彥超.某船機(jī)艙通風(fēng)系統(tǒng)模擬分析及優(yōu)化設(shè)計(jì)[D].上海:上海交通大學(xué),2011.

    [7] 周俊男.艦船艙室氣流組織的數(shù)值與實(shí)驗(yàn)研究[D].哈爾濱:哈爾濱工程大學(xué),2011.

    [8] 王瑞金,張凱,王剛.Fluent技術(shù)基礎(chǔ)與應(yīng)用實(shí)例[M].北京:清華大學(xué)出版社,2007:1-47.

    [9] 陸耀慶.實(shí)用供熱空調(diào)設(shè)計(jì)手冊(cè)[M].北京:中國(guó)建筑工業(yè)出版社,1993:728-730.

    [10] 建筑環(huán)境與節(jié)能研究院.民用建筑供暖通風(fēng)與空氣調(diào)節(jié)設(shè)計(jì)規(guī)范:GB 50736-2012[S].北京:中國(guó)建筑工業(yè)出版社,2012.

    [責(zé)任編輯:喻菁]

    Numerical simulation of the air distribution in a typical ship cabin

    QUAN Chongren1,WANG Yang2,YU Liqing3,CHEN Qian4,XIE Junlong4
    1 Shenyang Military Representative Department,Naval Armament Department of PLAN,Shenyang 110031,China

    2 China Ship Development and Design Center,Wuhan 430064,China

    3 Dalian Shipbuilding Industry Co.Ltd.,Dalian 116011,China

    4 Department of Building Environment and Energy Engineering,Huazhong University of Science and Technology,Wuhan 430074,China

    This paper presents a study of the air distribution of a typical two-person ship cabin with air dis?tributors being the air terminal due to the fact that the air terminal forms used in ship cabins are very differ?ent from those used in air-conditioning systems in civil buildings.The fluid-domain model of the typical two-person ship cabin is developed by using CFD technology,and the air distribution under the summer de?sign condition and the winter design condition are both simulated by imposing the coincident boundary con?ditions of the real ship,and the velocity and temperature fields of typical sections of the ship cabin are ana?lyzed.The simulation results show that under the summer design condition,the velocities and temperatures close to the air distributor are relatively high and low,respectively,while in the occupation area,the veloci?ty distribution is basically uniform,and the values are less than 0.3 m/s.The temperature distribution is al?so basically uniform,and the temperatures are about 26~27℃.Under the winter design condition,the air velocities in most areas are relatively low(less than 0.2 m/s),and the temperatures are approximately 20℃, except for that the velocity and temperature are relatively high around the air distributor.In brief,the pre?sented conditions successfully meet the comfort requirements under both the summer design condition and the winter design condition.The result further demonstrates the reasonability of the air-conditioning sys?tem designed for typical ship cabins.

    ship cabin;air distributor;numerical simulation;air distribution;CFD

    U664.86

    A

    10.3969/j.issn.1673-3185.2015.06.016

    http://www.cnki.net/kcms/detail/42.1755.TJ.20151110.1026.032.html期刊網(wǎng)址:www.ship-research.com

    權(quán)崇仁,王洋,于立慶,等.艦船典型艙室氣流組織數(shù)值模擬分析[J].中國(guó)艦船研究,2015,10(6):107-113. QUAN Chongren,WANG Yang,YU Liqing,et al.Numerical simulation of the air distribution in a typical ship cabin[J].Chinese Journal of Ship Research,2015,10(6):107-113.

    2015-03-23 < class="emphasis_bold"> 網(wǎng)絡(luò)出版時(shí)間:

    時(shí)間:2015-11-10 10:26

    國(guó)家部委基金資助項(xiàng)目

    權(quán)崇仁,男,1982年生,碩士,工程師。研究方向:裝備質(zhì)量管理王洋(通信作者),男,1981年生,碩士,工程師。研究方向:船舶空調(diào)通風(fēng)。E-mail:philips_king@163.com

    猜你喜歡
    艙室艦船云圖
    艦船通信中的噪聲消除研究
    船舶艙室通風(fēng)控制策略研究
    艦船測(cè)風(fēng)傳感器安裝位置數(shù)值仿真
    成都云圖控股股份有限公司
    黃強(qiáng)先生作品《雨后松云圖》
    名家名作(2017年3期)2017-09-15 11:13:37
    基于TV-L1分解的紅外云圖超分辨率算法
    艦船腐蝕預(yù)防與控制系統(tǒng)工程
    云圖青石板
    艦載導(dǎo)彈意外點(diǎn)火時(shí)艙室內(nèi)噴灑抑制的影響研究
    艦艇受損艙室的破損模式識(shí)別算法
    船海工程(2013年6期)2013-03-11 18:57:37
    一级爰片在线观看| 国产老妇女一区| 日日啪夜夜撸| 亚洲国产欧美人成| 国内精品宾馆在线| 精品一区二区三卡| 国产高潮美女av| 亚洲av男天堂| 卡戴珊不雅视频在线播放| 成人二区视频| 91精品一卡2卡3卡4卡| 边亲边吃奶的免费视频| 中文精品一卡2卡3卡4更新| 一二三四中文在线观看免费高清| 中文字幕免费在线视频6| 免费少妇av软件| 久久女婷五月综合色啪小说 | 激情五月婷婷亚洲| 丝袜脚勾引网站| 国产色爽女视频免费观看| 亚洲av二区三区四区| 亚洲成人精品中文字幕电影| 91久久精品国产一区二区三区| av卡一久久| 在线免费十八禁| 日韩免费高清中文字幕av| 日韩免费高清中文字幕av| 99热6这里只有精品| 国产欧美日韩精品一区二区| 神马国产精品三级电影在线观看| av在线天堂中文字幕| 国产人妻一区二区三区在| 丝瓜视频免费看黄片| 久久久久久久久久久免费av| 亚洲av成人精品一二三区| 中国国产av一级| 久久久久国产精品人妻一区二区| 色5月婷婷丁香| 日日啪夜夜撸| av在线天堂中文字幕| 99热6这里只有精品| 亚洲人成网站在线观看播放| 黄片无遮挡物在线观看| 久久97久久精品| 国产精品伦人一区二区| 91狼人影院| 男女边摸边吃奶| 午夜激情久久久久久久| 亚洲欧美中文字幕日韩二区| 男人添女人高潮全过程视频| 日韩一本色道免费dvd| 日本猛色少妇xxxxx猛交久久| 五月天丁香电影| 日本一本二区三区精品| 一级毛片黄色毛片免费观看视频| 九色成人免费人妻av| 国产老妇伦熟女老妇高清| 成人综合一区亚洲| 亚洲一区二区三区欧美精品 | 国产一区二区三区综合在线观看 | 欧美日韩视频高清一区二区三区二| 国产乱来视频区| 国产精品不卡视频一区二区| av国产免费在线观看| 91在线精品国自产拍蜜月| 美女视频免费永久观看网站| 国产男女超爽视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 美女cb高潮喷水在线观看| 国产高清三级在线| 男女啪啪激烈高潮av片| 久久综合国产亚洲精品| 亚洲久久久久久中文字幕| 亚洲成色77777| 九九久久精品国产亚洲av麻豆| 午夜福利视频精品| 久久久久久久亚洲中文字幕| 免费观看av网站的网址| 大码成人一级视频| 免费在线观看成人毛片| 少妇高潮的动态图| 黑人高潮一二区| 大香蕉97超碰在线| 人人妻人人看人人澡| 神马国产精品三级电影在线观看| 一个人观看的视频www高清免费观看| 精品99又大又爽又粗少妇毛片| 女人久久www免费人成看片| 一区二区av电影网| 热99国产精品久久久久久7| 夜夜爽夜夜爽视频| 国产永久视频网站| 舔av片在线| 中文字幕免费在线视频6| 久久精品国产亚洲av涩爱| 欧美日韩在线观看h| 久久久久久久精品精品| 中文在线观看免费www的网站| 日韩欧美一区视频在线观看 | 在线观看人妻少妇| 中文天堂在线官网| 亚洲欧美日韩卡通动漫| 国产 一区 欧美 日韩| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲最大av| 亚洲精品第二区| 麻豆成人av视频| 少妇被粗大猛烈的视频| 国产在视频线精品| 欧美一级a爱片免费观看看| 欧美精品人与动牲交sv欧美| 色视频在线一区二区三区| 国产91av在线免费观看| 大话2 男鬼变身卡| 日韩强制内射视频| 日本三级黄在线观看| 最近中文字幕2019免费版| 成人国产av品久久久| 一级毛片黄色毛片免费观看视频| 国产男人的电影天堂91| 亚洲成人一二三区av| 欧美激情久久久久久爽电影| 亚洲av福利一区| 少妇人妻精品综合一区二区| 十八禁网站网址无遮挡 | 欧美日韩在线观看h| 久久6这里有精品| 久久精品国产鲁丝片午夜精品| 欧美最新免费一区二区三区| 一级毛片电影观看| 亚洲久久久久久中文字幕| av.在线天堂| 日韩av免费高清视频| 久久精品国产亚洲网站| av天堂中文字幕网| 综合色av麻豆| 久久精品国产a三级三级三级| 永久网站在线| 校园人妻丝袜中文字幕| 亚洲在线观看片| 日韩av在线免费看完整版不卡| 免费大片18禁| 国产色爽女视频免费观看| 午夜福利高清视频| 天天躁日日操中文字幕| 小蜜桃在线观看免费完整版高清| 黄片wwwwww| 久久国内精品自在自线图片| 搞女人的毛片| 国产成人精品婷婷| 少妇熟女欧美另类| 成年女人在线观看亚洲视频 | 国产成人福利小说| 一级毛片我不卡| 韩国高清视频一区二区三区| 久久久久久久国产电影| 亚洲,欧美,日韩| 欧美日韩综合久久久久久| 欧美精品人与动牲交sv欧美| 久久精品久久久久久噜噜老黄| 大片电影免费在线观看免费| 69av精品久久久久久| 国产精品99久久99久久久不卡 | 亚洲av不卡在线观看| 免费黄频网站在线观看国产| 日产精品乱码卡一卡2卡三| 啦啦啦中文免费视频观看日本| 日日啪夜夜爽| 一区二区三区精品91| 91久久精品国产一区二区成人| 久久久a久久爽久久v久久| 久久久久久九九精品二区国产| 99九九线精品视频在线观看视频| 久久久久久久久久久丰满| 人妻夜夜爽99麻豆av| av在线天堂中文字幕| 久久精品国产a三级三级三级| 日韩欧美 国产精品| av在线天堂中文字幕| 国产在线一区二区三区精| 两个人的视频大全免费| 青青草视频在线视频观看| videos熟女内射| 日韩av在线免费看完整版不卡| 69av精品久久久久久| 成人综合一区亚洲| 看非洲黑人一级黄片| 毛片女人毛片| 伦理电影大哥的女人| 午夜激情久久久久久久| 亚洲久久久久久中文字幕| 边亲边吃奶的免费视频| 久久97久久精品| 精品久久久噜噜| 久久国产乱子免费精品| 狂野欧美白嫩少妇大欣赏| 网址你懂的国产日韩在线| 国产av不卡久久| 特大巨黑吊av在线直播| 日日啪夜夜撸| 好男人在线观看高清免费视频| 国产成人一区二区在线| 精品久久国产蜜桃| 直男gayav资源| 男女边摸边吃奶| av卡一久久| 波多野结衣巨乳人妻| 亚洲欧洲国产日韩| 建设人人有责人人尽责人人享有的 | 国产午夜精品一二区理论片| 最近最新中文字幕免费大全7| 秋霞在线观看毛片| 亚洲国产高清在线一区二区三| 中国国产av一级| 国产综合懂色| 亚洲熟女精品中文字幕| 亚洲欧美中文字幕日韩二区| 高清午夜精品一区二区三区| av在线播放精品| av专区在线播放| 日本wwww免费看| 少妇 在线观看| 欧美最新免费一区二区三区| 69人妻影院| 免费少妇av软件| 另类亚洲欧美激情| 精品久久久久久电影网| 亚洲成人一二三区av| 人体艺术视频欧美日本| 日韩成人av中文字幕在线观看| 国产一区有黄有色的免费视频| 黄片无遮挡物在线观看| 国产黄色视频一区二区在线观看| 高清午夜精品一区二区三区| 嫩草影院入口| 特级一级黄色大片| av一本久久久久| 久久久久久久久久久丰满| 亚洲精品,欧美精品| 日本一二三区视频观看| 亚洲欧美日韩东京热| 国产亚洲91精品色在线| 欧美极品一区二区三区四区| 亚洲欧美精品专区久久| 久久久久国产精品人妻一区二区| 在线看a的网站| 91在线精品国自产拍蜜月| 亚洲丝袜综合中文字幕| 精品一区二区免费观看| 中文字幕制服av| 在线观看av片永久免费下载| 热99国产精品久久久久久7| 岛国毛片在线播放| 久久人人爽av亚洲精品天堂 | 国产乱人视频| 99久国产av精品国产电影| 一个人看视频在线观看www免费| 九色成人免费人妻av| 伊人久久国产一区二区| 成人国产麻豆网| 韩国av在线不卡| 身体一侧抽搐| 久久久久九九精品影院| 三级经典国产精品| 久久久久久久精品精品| 人人妻人人看人人澡| 免费不卡的大黄色大毛片视频在线观看| 国产v大片淫在线免费观看| 国产永久视频网站| 亚洲国产精品成人久久小说| 成年女人看的毛片在线观看| 亚洲va在线va天堂va国产| 99re6热这里在线精品视频| 九色成人免费人妻av| 极品少妇高潮喷水抽搐| eeuss影院久久| 一级黄片播放器| 亚洲欧美精品自产自拍| videossex国产| 欧美激情国产日韩精品一区| 欧美性感艳星| 中文资源天堂在线| 国产精品爽爽va在线观看网站| 国产综合懂色| 国产精品一及| 亚洲av在线观看美女高潮| 国产乱人视频| 亚洲欧美日韩卡通动漫| 成人一区二区视频在线观看| 亚洲精品成人久久久久久| 好男人视频免费观看在线| 免费少妇av软件| 国产永久视频网站| 亚洲欧美清纯卡通| 国产精品一区二区在线观看99| 久久综合国产亚洲精品| 中国美白少妇内射xxxbb| 中文精品一卡2卡3卡4更新| 久久精品国产鲁丝片午夜精品| 哪个播放器可以免费观看大片| 国产黄色视频一区二区在线观看| 在线看a的网站| 最后的刺客免费高清国语| 日本熟妇午夜| av国产精品久久久久影院| 午夜免费男女啪啪视频观看| 男女边吃奶边做爰视频| 欧美潮喷喷水| 亚洲自偷自拍三级| 国产综合懂色| 国产精品国产三级国产av玫瑰| 狠狠精品人妻久久久久久综合| 亚洲真实伦在线观看| 国产精品一区www在线观看| 全区人妻精品视频| 一级毛片我不卡| 久久久久久久久久久免费av| av国产精品久久久久影院| 男女无遮挡免费网站观看| 国产黄频视频在线观看| 亚洲va在线va天堂va国产| 国产精品国产av在线观看| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 国产 一区 欧美 日韩| 精品国产三级普通话版| 日本三级黄在线观看| 看非洲黑人一级黄片| 成人二区视频| 成人黄色视频免费在线看| 99热这里只有是精品在线观看| 国产爱豆传媒在线观看| 乱系列少妇在线播放| 我的老师免费观看完整版| 久久久久九九精品影院| 亚洲第一区二区三区不卡| 日韩欧美精品免费久久| 在线天堂最新版资源| 97热精品久久久久久| 成人免费观看视频高清| 18禁在线无遮挡免费观看视频| 岛国毛片在线播放| 一个人观看的视频www高清免费观看| 日韩,欧美,国产一区二区三区| 在线观看美女被高潮喷水网站| 中国美白少妇内射xxxbb| 亚洲最大成人av| 另类亚洲欧美激情| 国产亚洲最大av| 赤兔流量卡办理| 男女边吃奶边做爰视频| 免费看不卡的av| 国产亚洲午夜精品一区二区久久 | 欧美日韩视频高清一区二区三区二| 免费看不卡的av| 亚洲精品国产成人久久av| 久久久精品免费免费高清| 精品久久久久久久久av| 成人欧美大片| 日韩强制内射视频| 麻豆成人av视频| 免费大片18禁| 黄色配什么色好看| 亚洲婷婷狠狠爱综合网| 成人亚洲欧美一区二区av| 成人无遮挡网站| av天堂中文字幕网| 免费观看的影片在线观看| 全区人妻精品视频| 国产亚洲av嫩草精品影院| 日本免费在线观看一区| 国内精品宾馆在线| 久久久久久九九精品二区国产| 午夜爱爱视频在线播放| 亚洲精品视频女| 少妇人妻精品综合一区二区| 国产精品偷伦视频观看了| 久久99热6这里只有精品| 久久99蜜桃精品久久| 国产美女午夜福利| 亚洲精品自拍成人| 色视频www国产| 黄片无遮挡物在线观看| 午夜精品国产一区二区电影 | 日本黄色片子视频| 亚洲成人中文字幕在线播放| 精品99又大又爽又粗少妇毛片| 男女下面进入的视频免费午夜| 亚洲久久久久久中文字幕| 草草在线视频免费看| 99久久人妻综合| 国产成人91sexporn| 尤物成人国产欧美一区二区三区| 亚洲自偷自拍三级| 赤兔流量卡办理| 我要看日韩黄色一级片| 国产成人freesex在线| 在线观看美女被高潮喷水网站| 不卡视频在线观看欧美| 亚洲国产精品成人综合色| 不卡视频在线观看欧美| 王馨瑶露胸无遮挡在线观看| 亚洲国产高清在线一区二区三| 久久久色成人| 日韩av在线免费看完整版不卡| 久久精品国产亚洲网站| 婷婷色综合大香蕉| 新久久久久国产一级毛片| 亚州av有码| 免费人成在线观看视频色| 日韩一本色道免费dvd| 一级黄片播放器| 久久精品人妻少妇| 秋霞在线观看毛片| eeuss影院久久| 国产在线一区二区三区精| 最近最新中文字幕免费大全7| 国产日韩欧美在线精品| 麻豆国产97在线/欧美| 国产成人午夜福利电影在线观看| 久久99热这里只频精品6学生| 国产精品.久久久| videossex国产| 日本-黄色视频高清免费观看| av女优亚洲男人天堂| 亚洲精品乱码久久久久久按摩| 麻豆乱淫一区二区| 26uuu在线亚洲综合色| 日本三级黄在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲三级黄色毛片| 青春草亚洲视频在线观看| 一个人看视频在线观看www免费| 色播亚洲综合网| 美女脱内裤让男人舔精品视频| 好男人在线观看高清免费视频| 日韩av在线免费看完整版不卡| 成人鲁丝片一二三区免费| 少妇被粗大猛烈的视频| 国产色爽女视频免费观看| 麻豆国产97在线/欧美| 日日摸夜夜添夜夜添av毛片| 日本av手机在线免费观看| 波野结衣二区三区在线| 老司机影院成人| 国产乱人视频| 久久久久久久大尺度免费视频| 一区二区av电影网| 国产探花在线观看一区二区| 一级毛片久久久久久久久女| 夫妻性生交免费视频一级片| 中文乱码字字幕精品一区二区三区| 亚洲精品亚洲一区二区| 五月伊人婷婷丁香| 久久精品国产a三级三级三级| 亚洲自偷自拍三级| 91精品一卡2卡3卡4卡| 69人妻影院| 美女脱内裤让男人舔精品视频| 亚洲国产精品成人久久小说| 国产真实伦视频高清在线观看| 九九爱精品视频在线观看| 丰满人妻一区二区三区视频av| 久热这里只有精品99| 色视频在线一区二区三区| 免费av不卡在线播放| 久久久色成人| 91aial.com中文字幕在线观看| 日本三级黄在线观看| 久久人人爽人人爽人人片va| 伦精品一区二区三区| 午夜福利视频精品| 日韩欧美 国产精品| 美女高潮的动态| 九九久久精品国产亚洲av麻豆| 夜夜爽夜夜爽视频| 少妇裸体淫交视频免费看高清| 网址你懂的国产日韩在线| 国产精品嫩草影院av在线观看| 国产黄色免费在线视频| 蜜桃久久精品国产亚洲av| 777米奇影视久久| 熟女av电影| 大片免费播放器 马上看| 日本黄大片高清| 久久精品久久久久久久性| 九九久久精品国产亚洲av麻豆| 男女无遮挡免费网站观看| 综合色丁香网| 肉色欧美久久久久久久蜜桃 | 国产精品秋霞免费鲁丝片| 精品一区二区三卡| 嫩草影院新地址| 插逼视频在线观看| 99热6这里只有精品| 日韩精品有码人妻一区| 亚洲丝袜综合中文字幕| 日本爱情动作片www.在线观看| 国产有黄有色有爽视频| 三级国产精品片| 男女下面进入的视频免费午夜| 精品久久久噜噜| 青青草视频在线视频观看| 人体艺术视频欧美日本| 麻豆乱淫一区二区| 久久99热这里只有精品18| 成人欧美大片| 亚洲国产精品成人综合色| 成人毛片60女人毛片免费| 热re99久久精品国产66热6| 免费av毛片视频| 国产精品国产av在线观看| 91精品伊人久久大香线蕉| 欧美精品一区二区大全| 国产精品蜜桃在线观看| 久久人人爽人人爽人人片va| 日本熟妇午夜| 男人爽女人下面视频在线观看| 街头女战士在线观看网站| 乱码一卡2卡4卡精品| 亚洲精品国产色婷婷电影| 亚洲成色77777| 亚洲欧美一区二区三区国产| 国内揄拍国产精品人妻在线| 午夜免费观看性视频| 91久久精品国产一区二区成人| 亚洲精品aⅴ在线观看| 美女视频免费永久观看网站| 在线观看三级黄色| 十八禁网站网址无遮挡 | 男男h啪啪无遮挡| 国产永久视频网站| 日本三级黄在线观看| 欧美另类一区| 国产男女超爽视频在线观看| 久久人人爽人人片av| 亚洲怡红院男人天堂| 国产精品国产av在线观看| 国产免费一区二区三区四区乱码| 神马国产精品三级电影在线观看| av黄色大香蕉| 亚洲av电影在线观看一区二区三区 | 超碰av人人做人人爽久久| 最近中文字幕2019免费版| a级毛色黄片| 99久久中文字幕三级久久日本| 亚洲无线观看免费| 亚洲精品国产av成人精品| 国产精品偷伦视频观看了| 国产色爽女视频免费观看| 午夜亚洲福利在线播放| 午夜福利网站1000一区二区三区| 中文字幕制服av| 男女边吃奶边做爰视频| 老司机影院毛片| 国产视频内射| 国产精品人妻久久久影院| 免费看不卡的av| 免费av毛片视频| 国产av不卡久久| 国产高清三级在线| 午夜激情久久久久久久| 可以在线观看毛片的网站| 黑人高潮一二区| 亚洲人成网站高清观看| 一边亲一边摸免费视频| 精品久久久噜噜| 久久久精品欧美日韩精品| 你懂的网址亚洲精品在线观看| 美女高潮的动态| 国产精品一及| 精品视频人人做人人爽| 中国三级夫妇交换| 18禁在线播放成人免费| 久久久久久久久久成人| 午夜亚洲福利在线播放| 精品酒店卫生间| 王馨瑶露胸无遮挡在线观看| 日本免费在线观看一区| 亚洲国产色片| 91aial.com中文字幕在线观看| 蜜桃久久精品国产亚洲av| 国产欧美另类精品又又久久亚洲欧美| 国产精品国产av在线观看| 日韩免费高清中文字幕av| 蜜臀久久99精品久久宅男| 欧美日韩亚洲高清精品| 一级毛片 在线播放| 成年女人在线观看亚洲视频 | 亚洲精品第二区| 亚洲欧美精品专区久久| 国产成人a区在线观看| 女人十人毛片免费观看3o分钟| 男女啪啪激烈高潮av片| 国产乱人偷精品视频| 午夜免费鲁丝| 久久热精品热| 国产综合懂色| 女人十人毛片免费观看3o分钟| 午夜日本视频在线| 一级爰片在线观看| 国产精品久久久久久av不卡| 一区二区三区精品91| 男女无遮挡免费网站观看| av国产久精品久网站免费入址| 国产午夜精品久久久久久一区二区三区| 联通29元200g的流量卡| 男男h啪啪无遮挡| 国产精品无大码| 80岁老熟妇乱子伦牲交| 18禁裸乳无遮挡免费网站照片| 婷婷色综合大香蕉| 日韩在线高清观看一区二区三区| 好男人在线观看高清免费视频| 久久久成人免费电影|