• <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
    最近最新中文字幕大全电影3 | 天天操日日干夜夜撸| 在线观看日韩欧美| 正在播放国产对白刺激| a级毛片黄视频| 午夜免费成人在线视频| 中出人妻视频一区二区| 成人永久免费在线观看视频| 亚洲少妇的诱惑av| 国产一区二区三区在线臀色熟女 | 中文字幕最新亚洲高清| 亚洲欧美激情在线| 丰满的人妻完整版| 交换朋友夫妻互换小说| 法律面前人人平等表现在哪些方面| 成人永久免费在线观看视频| 一级a爱片免费观看的视频| 国产精品一区二区在线观看99| 国产av又大| 国产精品免费大片| 欧美不卡视频在线免费观看 | 叶爱在线成人免费视频播放| 久久久国产成人精品二区 | 黑人巨大精品欧美一区二区蜜桃| 人人妻人人澡人人爽人人夜夜| 久9热在线精品视频| 成人18禁高潮啪啪吃奶动态图| 激情在线观看视频在线高清 | 久久国产精品人妻蜜桃| 亚洲成人免费电影在线观看| 一本一本久久a久久精品综合妖精| 国产精品 国内视频| 自拍欧美九色日韩亚洲蝌蚪91| 好看av亚洲va欧美ⅴa在| 久久人人97超碰香蕉20202| 国产精品二区激情视频| 免费观看a级毛片全部| 黄色a级毛片大全视频| 九色亚洲精品在线播放| 91麻豆精品激情在线观看国产 | 成人影院久久| av免费在线观看网站| 少妇粗大呻吟视频| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av高清一级| 久久精品91无色码中文字幕| 久久国产精品人妻蜜桃| 国产欧美日韩一区二区三区在线| 久久久国产精品麻豆| av福利片在线| 黄色怎么调成土黄色| 婷婷成人精品国产| 国产午夜精品久久久久久| 后天国语完整版免费观看| 精品第一国产精品| 黄网站色视频无遮挡免费观看| 丰满的人妻完整版| 免费观看a级毛片全部| 国产99久久九九免费精品| 天堂俺去俺来也www色官网| 99热只有精品国产| 免费在线观看黄色视频的| 免费在线观看日本一区| 99久久人妻综合| 91成年电影在线观看| 国产国语露脸激情在线看| 国精品久久久久久国模美| ponron亚洲| 日韩三级视频一区二区三区| 日本wwww免费看| 中文字幕高清在线视频| 国产在线一区二区三区精| 国产一区二区激情短视频| 亚洲人成电影免费在线| 色播在线永久视频| 村上凉子中文字幕在线| 中文字幕人妻丝袜制服| av天堂久久9| 国产亚洲精品久久久久久毛片 | 精品卡一卡二卡四卡免费| 嫩草影视91久久| 人人妻人人澡人人爽人人夜夜| 色婷婷av一区二区三区视频| 国产精品国产高清国产av | 国精品久久久久久国模美| 精品电影一区二区在线| 亚洲中文av在线| 精品一区二区三区四区五区乱码| 黄色毛片三级朝国网站| 极品少妇高潮喷水抽搐| 男女下面插进去视频免费观看| 成熟少妇高潮喷水视频| 亚洲欧美激情在线| 人人妻人人澡人人爽人人夜夜| 午夜老司机福利片| 丰满迷人的少妇在线观看| 久久久久国产一级毛片高清牌| 国产精品影院久久| 欧美性长视频在线观看| 女人精品久久久久毛片| 岛国毛片在线播放| 一本大道久久a久久精品| 精品国产乱码久久久久久男人| 欧美激情久久久久久爽电影 | 夜夜夜夜夜久久久久| 免费人成视频x8x8入口观看| 国产成人精品无人区| 91国产中文字幕| 首页视频小说图片口味搜索| 99精品在免费线老司机午夜| 老司机在亚洲福利影院| 黄色毛片三级朝国网站| 亚洲专区字幕在线| 精品少妇久久久久久888优播| 69精品国产乱码久久久| 80岁老熟妇乱子伦牲交| 变态另类成人亚洲欧美熟女 | 国产精品免费视频内射| 久久精品国产综合久久久| x7x7x7水蜜桃| 99久久精品国产亚洲精品| 亚洲在线自拍视频| 日本a在线网址| 91av网站免费观看| 18禁裸乳无遮挡免费网站照片 | 18禁国产床啪视频网站| 美女高潮到喷水免费观看| 极品教师在线免费播放| 欧美精品一区二区免费开放| 777米奇影视久久| 99精品在免费线老司机午夜| 国产精品99久久99久久久不卡| 亚洲精品在线观看二区| 亚洲av美国av| 男女下面插进去视频免费观看| √禁漫天堂资源中文www| 国产亚洲精品久久久久久毛片 | 欧美性长视频在线观看| 亚洲av成人av| 精品欧美一区二区三区在线| 亚洲五月色婷婷综合| 人人妻人人爽人人添夜夜欢视频| 国产区一区二久久| 一边摸一边做爽爽视频免费| 免费高清在线观看日韩| 黄色丝袜av网址大全| 女人被躁到高潮嗷嗷叫费观| 亚洲av日韩在线播放| 久久午夜综合久久蜜桃| 涩涩av久久男人的天堂| 亚洲五月婷婷丁香| 99久久99久久久精品蜜桃| 国产淫语在线视频| 午夜精品在线福利| 成人永久免费在线观看视频| 亚洲av成人av| 欧美日韩亚洲国产一区二区在线观看 | 丝袜美足系列| 久久婷婷成人综合色麻豆| 久久精品国产清高在天天线| 大陆偷拍与自拍| 亚洲国产欧美一区二区综合| 国产97色在线日韩免费| 欧美久久黑人一区二区| 超碰成人久久| 老熟女久久久| 在线观看免费视频网站a站| 女性生殖器流出的白浆| a级毛片在线看网站| 首页视频小说图片口味搜索| av国产精品久久久久影院| 亚洲成人免费电影在线观看| 51午夜福利影视在线观看| 777米奇影视久久| 午夜两性在线视频| 91av网站免费观看| 国精品久久久久久国模美| 亚洲欧美色中文字幕在线| 一边摸一边抽搐一进一小说 | 欧美激情高清一区二区三区| 国产精品久久电影中文字幕 | 十八禁人妻一区二区| 久久久久久久久免费视频了| 波多野结衣av一区二区av| 嫁个100分男人电影在线观看| 色综合婷婷激情| 精品国产一区二区三区四区第35| 别揉我奶头~嗯~啊~动态视频| 99国产精品一区二区三区| 99精品欧美一区二区三区四区| e午夜精品久久久久久久| 曰老女人黄片| 成人亚洲精品一区在线观看| 国产人伦9x9x在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产熟女午夜一区二区三区| 少妇 在线观看| 欧美激情极品国产一区二区三区| 正在播放国产对白刺激| 国产一区有黄有色的免费视频| 国产精华一区二区三区| 欧美国产精品va在线观看不卡| 99re6热这里在线精品视频| e午夜精品久久久久久久| 国产视频一区二区在线看| 国产三级黄色录像| 国产高清videossex| 乱人伦中国视频| 亚洲国产欧美一区二区综合| 亚洲国产欧美网| 另类亚洲欧美激情| 女人被狂操c到高潮| a级毛片黄视频| 在线永久观看黄色视频| 亚洲 欧美一区二区三区| 后天国语完整版免费观看| 国产成人欧美在线观看 | 欧美日本中文国产一区发布| 成年人黄色毛片网站| 丁香六月欧美| 老熟女久久久| 欧美日本中文国产一区发布| 久久久久国内视频| 久久精品国产99精品国产亚洲性色 | 欧美成人免费av一区二区三区 | 国产又色又爽无遮挡免费看| 老鸭窝网址在线观看| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕精品免费在线观看视频| 亚洲全国av大片| 一个人免费在线观看的高清视频| 亚洲人成电影免费在线| 久久国产亚洲av麻豆专区| av线在线观看网站| 国产精品一区二区在线观看99| 麻豆av在线久日| 亚洲成国产人片在线观看| 王馨瑶露胸无遮挡在线观看| 久久精品国产a三级三级三级| 大香蕉久久网| 国产精华一区二区三区| 欧美黄色淫秽网站| 老司机亚洲免费影院| 国产成人系列免费观看| 免费女性裸体啪啪无遮挡网站| 99精品久久久久人妻精品| 高潮久久久久久久久久久不卡| 精品国产超薄肉色丝袜足j| 自线自在国产av| 精品国产一区二区三区久久久樱花| 午夜福利在线观看吧| 久久精品亚洲av国产电影网| 亚洲精品久久午夜乱码| 一二三四在线观看免费中文在| 国产午夜精品久久久久久| 国产色视频综合| 国产激情欧美一区二区| 色老头精品视频在线观看| 老汉色av国产亚洲站长工具| 欧美色视频一区免费| 黄片播放在线免费| 日韩中文字幕欧美一区二区| 国产精品.久久久| 亚洲第一欧美日韩一区二区三区| 亚洲午夜精品一区,二区,三区| 一区二区日韩欧美中文字幕| 狠狠狠狠99中文字幕| 精品一区二区三区av网在线观看| 国产精品亚洲一级av第二区| 午夜福利在线免费观看网站| 十分钟在线观看高清视频www| 女人爽到高潮嗷嗷叫在线视频| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 国产精品亚洲av一区麻豆| 水蜜桃什么品种好| 亚洲成av片中文字幕在线观看| 国产精华一区二区三区| 一级毛片精品| 涩涩av久久男人的天堂| 18在线观看网站| 热re99久久国产66热| 精品视频人人做人人爽| 电影成人av| 99热网站在线观看| 手机成人av网站| 精品国产美女av久久久久小说| 波多野结衣一区麻豆| 啦啦啦在线免费观看视频4| 女人高潮潮喷娇喘18禁视频| 久久久水蜜桃国产精品网| av福利片在线| 美女福利国产在线| 一区二区日韩欧美中文字幕| 亚洲黑人精品在线| 黑人猛操日本美女一级片| 久久久久久免费高清国产稀缺| 十八禁网站免费在线| 国产精品久久视频播放| 91九色精品人成在线观看| 精品视频人人做人人爽| 麻豆国产av国片精品| 国产精品国产高清国产av | 老司机午夜福利在线观看视频| 窝窝影院91人妻| 精品国产超薄肉色丝袜足j| 捣出白浆h1v1| 久久久国产成人免费| 一级毛片高清免费大全| 乱人伦中国视频| 中国美女看黄片| 色播在线永久视频| 国产一区有黄有色的免费视频| 国产精品久久久久久人妻精品电影| 国产片内射在线| 精品国产超薄肉色丝袜足j| 亚洲av美国av| 日本撒尿小便嘘嘘汇集6| 女性被躁到高潮视频| 久久久精品免费免费高清| 亚洲aⅴ乱码一区二区在线播放 | 亚洲欧美一区二区三区黑人| 校园春色视频在线观看| 国产精品国产高清国产av | 精品久久久久久久久久免费视频 | 免费在线观看视频国产中文字幕亚洲| 搡老岳熟女国产| 男女之事视频高清在线观看| 国产亚洲欧美精品永久| 国产无遮挡羞羞视频在线观看| 亚洲一码二码三码区别大吗| 欧美久久黑人一区二区| 韩国精品一区二区三区| 亚洲av欧美aⅴ国产| 亚洲国产中文字幕在线视频| 国产成人欧美| 亚洲精品中文字幕一二三四区| 亚洲在线自拍视频| videos熟女内射| 午夜福利在线观看吧| 久久热在线av| tube8黄色片| 操出白浆在线播放| 曰老女人黄片| 岛国毛片在线播放| 欧美成人免费av一区二区三区 | 欧美老熟妇乱子伦牲交| 久久天躁狠狠躁夜夜2o2o| 久久青草综合色| 无遮挡黄片免费观看| 18在线观看网站| 丁香六月欧美| 美女午夜性视频免费| 亚洲熟女精品中文字幕| 超碰成人久久| 国产91精品成人一区二区三区| 狂野欧美激情性xxxx| 老司机在亚洲福利影院| av有码第一页| 制服人妻中文乱码| 侵犯人妻中文字幕一二三四区| 亚洲av电影在线进入| 看免费av毛片| 两个人免费观看高清视频| 国产欧美日韩综合在线一区二区| a级毛片黄视频| 99精品欧美一区二区三区四区| 日本vs欧美在线观看视频| 一级黄色大片毛片| 亚洲久久久国产精品| 一级a爱片免费观看的视频| 制服诱惑二区| 国产精品电影一区二区三区 | 国产精品秋霞免费鲁丝片| 女同久久另类99精品国产91| 亚洲熟妇熟女久久| av一本久久久久| 极品少妇高潮喷水抽搐| 一级,二级,三级黄色视频| 黑人操中国人逼视频| 亚洲男人天堂网一区| www.精华液| 激情视频va一区二区三区| 国产一区二区三区在线臀色熟女 | av一本久久久久| 天天躁夜夜躁狠狠躁躁| 欧美日韩成人在线一区二区| 亚洲中文日韩欧美视频| 国产主播在线观看一区二区| 桃红色精品国产亚洲av| 国产亚洲精品一区二区www | 久久 成人 亚洲| 国产精品久久久av美女十八| 久久精品91无色码中文字幕| 超碰成人久久| 91成年电影在线观看| 免费av中文字幕在线| 日本五十路高清| 夜夜爽天天搞| 两个人看的免费小视频| 午夜影院日韩av| 久久久久久久精品吃奶| 午夜91福利影院| av有码第一页| 中文字幕最新亚洲高清| 好男人电影高清在线观看| 亚洲五月天丁香| 一本综合久久免费| 热99久久久久精品小说推荐| 亚洲成a人片在线一区二区| 亚洲国产欧美日韩在线播放| 精品国产乱码久久久久久男人| av国产精品久久久久影院| 另类亚洲欧美激情| 亚洲成国产人片在线观看| 飞空精品影院首页| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩福利视频一区二区| 精品国产国语对白av| 久久精品国产亚洲av香蕉五月 | 久久精品国产a三级三级三级| 国产精华一区二区三区| 又黄又爽又免费观看的视频| 露出奶头的视频| 亚洲色图综合在线观看| 一本一本久久a久久精品综合妖精| 国产亚洲欧美精品永久| 国产深夜福利视频在线观看| 香蕉国产在线看| 亚洲一区二区三区不卡视频| 精品人妻熟女毛片av久久网站| 成人18禁在线播放| av线在线观看网站| 欧美色视频一区免费| 成人特级黄色片久久久久久久| 51午夜福利影视在线观看| 亚洲欧美日韩高清在线视频| 在线永久观看黄色视频| 多毛熟女@视频| 天天影视国产精品| 热99re8久久精品国产| 国产无遮挡羞羞视频在线观看| 电影成人av| 国产在视频线精品| 黄片小视频在线播放| 欧美日韩亚洲国产一区二区在线观看 | 国产精品 国内视频| 一区二区三区国产精品乱码| 精品国产美女av久久久久小说| av网站在线播放免费| 村上凉子中文字幕在线| 精品熟女少妇八av免费久了| 黄色成人免费大全| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲自偷自拍图片 自拍| 国产精品二区激情视频| 在线播放国产精品三级| 精品亚洲成国产av| 欧美日韩福利视频一区二区| 亚洲三区欧美一区| 69av精品久久久久久| 美国免费a级毛片| 国产免费现黄频在线看| 国产1区2区3区精品| 1024香蕉在线观看| 亚洲综合色网址| av中文乱码字幕在线| x7x7x7水蜜桃| 亚洲全国av大片| 超色免费av| 手机成人av网站| 亚洲av第一区精品v没综合| 啦啦啦在线免费观看视频4| 91麻豆精品激情在线观看国产 | 日韩免费高清中文字幕av| 亚洲国产欧美一区二区综合| 女人被躁到高潮嗷嗷叫费观| 精品久久蜜臀av无| 亚洲男人天堂网一区| 在线国产一区二区在线| 亚洲一码二码三码区别大吗| 超碰97精品在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲精品久久成人aⅴ小说| 夜夜躁狠狠躁天天躁| 极品教师在线免费播放| 亚洲精品成人av观看孕妇| 一区在线观看完整版| 色在线成人网| 久久影院123| 亚洲中文av在线| 欧美激情久久久久久爽电影 | 色精品久久人妻99蜜桃| 纯流量卡能插随身wifi吗| 成人国语在线视频| 国产精品国产av在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲av熟女| 不卡一级毛片| 亚洲成a人片在线一区二区| 国产av一区二区精品久久| 亚洲精品av麻豆狂野| 狠狠婷婷综合久久久久久88av| 久久香蕉精品热| x7x7x7水蜜桃| 少妇被粗大的猛进出69影院| 看黄色毛片网站| 老鸭窝网址在线观看| 国产高清视频在线播放一区| 可以免费在线观看a视频的电影网站| 亚洲成人国产一区在线观看| 韩国av一区二区三区四区| 1024视频免费在线观看| 亚洲精品在线观看二区| 国产片内射在线| 成人亚洲精品一区在线观看| 国产精品免费视频内射| 国产精品秋霞免费鲁丝片| 久热这里只有精品99| 久久国产亚洲av麻豆专区| 天天操日日干夜夜撸| 成人av一区二区三区在线看| 精品久久蜜臀av无| 亚洲精品美女久久久久99蜜臀| 好看av亚洲va欧美ⅴa在| 最新的欧美精品一区二区| 伊人久久大香线蕉亚洲五| 日韩欧美在线二视频 | 欧美性长视频在线观看| 国产av又大| 天天影视国产精品| 欧美激情久久久久久爽电影 | 中国美女看黄片| 亚洲av成人不卡在线观看播放网| 99久久综合精品五月天人人| 深夜精品福利| 欧美性长视频在线观看| 国产熟女午夜一区二区三区| 精品视频人人做人人爽| 国产1区2区3区精品| 一级片免费观看大全| 成人国产一区最新在线观看| 男女高潮啪啪啪动态图| 大陆偷拍与自拍| 成人特级黄色片久久久久久久| 国产97色在线日韩免费| 国产精品免费一区二区三区在线 | 免费在线观看完整版高清| 成人永久免费在线观看视频| 欧美成狂野欧美在线观看| 久久久国产欧美日韩av| 人人妻人人添人人爽欧美一区卜| 午夜亚洲福利在线播放| 久久婷婷成人综合色麻豆| 亚洲综合色网址| 国产一区二区激情短视频| 欧美精品亚洲一区二区| 免费在线观看亚洲国产| 欧美激情 高清一区二区三区| 成人永久免费在线观看视频| 亚洲色图 男人天堂 中文字幕| 欧美在线黄色| 一级毛片高清免费大全| 18在线观看网站| 精品久久久久久,| 一级毛片女人18水好多| 夜夜爽天天搞| 日本五十路高清| 9色porny在线观看| а√天堂www在线а√下载 | 狂野欧美激情性xxxx| 国产免费现黄频在线看| 亚洲专区字幕在线| 成人黄色视频免费在线看| 久久精品熟女亚洲av麻豆精品| 男男h啪啪无遮挡| 欧美亚洲 丝袜 人妻 在线| 亚洲国产精品sss在线观看 | 热99国产精品久久久久久7| 别揉我奶头~嗯~啊~动态视频| 久久精品aⅴ一区二区三区四区| 王馨瑶露胸无遮挡在线观看| 国产色视频综合| 国产又色又爽无遮挡免费看| 亚洲成人手机| 久久性视频一级片| 国产不卡一卡二| 色播在线永久视频| 久久性视频一级片| 精品视频人人做人人爽| 精品国产国语对白av| 欧美成狂野欧美在线观看| 成年动漫av网址| 一级毛片高清免费大全| 亚洲五月婷婷丁香| 精品国产一区二区三区四区第35| 亚洲成人手机| 欧美成狂野欧美在线观看| 在线观看免费视频网站a站| 黑人操中国人逼视频| tocl精华| 国产蜜桃级精品一区二区三区 | 性色av乱码一区二区三区2| 国产欧美日韩一区二区三区在线| 正在播放国产对白刺激| 51午夜福利影视在线观看| 成人国产一区最新在线观看| 欧美久久黑人一区二区| 亚洲aⅴ乱码一区二区在线播放 | 免费看a级黄色片| 国产日韩一区二区三区精品不卡| 亚洲国产精品一区二区三区在线| 老汉色∧v一级毛片| 国产成人欧美在线观看 | 看片在线看免费视频|