亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        蔬菜氣霧栽培箱不同氣流循環(huán)方式的流場(chǎng)和溫度場(chǎng)CFD模擬

        2019-10-10 02:26:56楊希文羅亞輝龍莉霞石毅新胡文武
        關(guān)鍵詞:氣霧生菜氣流

        蔣 蘋,楊希文,羅亞輝,龍莉霞,石毅新,胡文武

        蔬菜氣霧栽培箱不同氣流循環(huán)方式的流場(chǎng)和溫度場(chǎng)CFD模擬

        蔣 蘋1,2,楊希文1,羅亞輝1,2,龍莉霞1,石毅新1,2,胡文武1,2

        (1. 湖南農(nóng)業(yè)大學(xué)工學(xué)院,長(zhǎng)沙 410128;2. 南方糧油作物協(xié)同創(chuàng)新中心,長(zhǎng)沙 410128)

        蔬菜氣霧栽培箱內(nèi)的空氣流動(dòng)和溫度對(duì)箱體內(nèi)的環(huán)境調(diào)節(jié)及農(nóng)作物生長(zhǎng)具有重要作用,農(nóng)作物周圍空氣的均勻性流動(dòng)能促進(jìn)農(nóng)作物的生長(zhǎng)速率。為探究蔬菜氣霧栽培箱內(nèi)流場(chǎng)及溫度場(chǎng)分布規(guī)律,基于計(jì)算機(jī)流體力學(xué)(CFD,computational fluid dynamics)方法,利用FLUENT軟件,結(jié)合標(biāo)準(zhǔn)湍流模型、有孔介質(zhì)模型、作物冠層質(zhì)熱交換模型等,建立了蔬菜氣霧栽培箱不同氣流循環(huán)方式下的CFD模型。并對(duì)氣霧栽培箱內(nèi)的環(huán)境進(jìn)行優(yōu)化設(shè)計(jì),設(shè)計(jì)了3種氣流循環(huán)方案:頂面進(jìn)側(cè)面出,側(cè)面進(jìn)頂面出,側(cè)面進(jìn)側(cè)面出。對(duì)送回風(fēng)口的不同位置布局進(jìn)行了研究,并對(duì)3種氣流循環(huán)方案進(jìn)行了數(shù)值模擬。模擬結(jié)果可知:氣流為頂面進(jìn)側(cè)面出方案中,風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值區(qū)域占58.1%,適宜溫度值區(qū)域占93.6%,通風(fēng)死角區(qū)域占比0.844%;氣流為側(cè)面進(jìn)頂面出方案中,生菜生長(zhǎng)適宜風(fēng)速值區(qū)域占59.6%,適宜溫度值區(qū)域占99.98% ,通風(fēng)死角區(qū)域占比0.069%;氣流為側(cè)面進(jìn)側(cè)面出方案中,風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值區(qū)域占54.3%,適宜溫度值區(qū)域占92.4%,通風(fēng)死角區(qū)域占比16.7%。分析對(duì)比后得到側(cè)面進(jìn)頂面出為最佳氣流循環(huán)方案。并對(duì)此進(jìn)行了試驗(yàn)測(cè)試,結(jié)果表明:氣霧栽培箱內(nèi)溫度、風(fēng)速模擬值和實(shí)測(cè)值進(jìn)行對(duì)比,溫度平均相對(duì)誤差為3.9%,均方根誤差為0.86℃。風(fēng)速平均相對(duì)誤差為3.5%,均方根誤差為0.26 m/s,模擬值和實(shí)測(cè)值誤差較小,模擬效果良好,驗(yàn)證了CFD模型的準(zhǔn)確性。該研究為蔬菜氣霧栽培箱內(nèi)的流場(chǎng)及溫度變化規(guī)律,內(nèi)部環(huán)境調(diào)節(jié),裝置優(yōu)化設(shè)計(jì)提供了參考依據(jù)。

        風(fēng)速;溫度;計(jì)算機(jī)流體力學(xué);優(yōu)化設(shè)計(jì);蔬菜氣霧栽培箱

        0 引 言

        農(nóng)作物氣霧栽培是無土栽培領(lǐng)域一種新型栽培方式,它不僅能提高作物水肥利用率,還能改善作物根系缺氧,提高土地利用率,減輕環(huán)境污染等。氣霧栽培模式主要有立柱式、金字塔式、梯形式[1]。目前部分學(xué)者對(duì)氣霧栽培設(shè)施也有相關(guān)研究,如高建民等[2]、劉義飛[3]、劉巖[4]通過設(shè)計(jì)新型的氣霧栽培設(shè)施進(jìn)行農(nóng)作物的氣霧栽培生產(chǎn),以此提高農(nóng)作物產(chǎn)量和品質(zhì)。本課題組設(shè)計(jì)了一種蔬菜氣霧栽培箱裝置,蔬菜氣霧栽培箱是一種密閉性農(nóng)作物栽培生產(chǎn)系統(tǒng),適宜在溫室大棚、植物工廠、家庭陽臺(tái)等場(chǎng)所對(duì)農(nóng)作物進(jìn)行栽培生產(chǎn)。影響蔬菜氣霧栽培箱內(nèi)農(nóng)作物生長(zhǎng)的主要因素有氣流速度、溫度、濕度、光照等。而氣流速度和溫度是影響蔬菜氣霧栽培箱內(nèi)農(nóng)作物冠層生長(zhǎng)的重要因素,氣流速度和溫度對(duì)農(nóng)作物的蒸騰作用及光合作用有著重要作用,影響農(nóng)作物的品質(zhì)和產(chǎn)量。研究氣流速度和溫度的變化規(guī)律對(duì)推廣氣霧栽培技術(shù)有重要意義。目前CFD(CFD,computational fluid dynamics)主要在溫室、塑料大棚、植物工廠等相關(guān)領(lǐng)域有研究。在氣霧栽培領(lǐng)域采用CFD方法進(jìn)行流場(chǎng)和溫度場(chǎng)的數(shù)值模擬缺乏理論研究。采用CFD數(shù)值模擬方法對(duì)提高氣霧栽培箱內(nèi)部農(nóng)作物栽培的生長(zhǎng)效率有著重要作用。目前相關(guān)研究有:王新忠等[5]利用CFD模型對(duì)全開型玻璃溫室自然通風(fēng)進(jìn)行了流場(chǎng)分析,結(jié)果表明側(cè)窗為45°,天窗調(diào)整至60°時(shí),溫室整體溫度從38.4 ℃降至36.9 ℃,能達(dá)到較好的降溫效果。翁衛(wèi)兵等[6]采用CFD方法對(duì)冷藏車廂內(nèi)流場(chǎng)的流動(dòng)特性進(jìn)行了模擬,結(jié)果表明根據(jù)特定夾角開啟隔熱氣簾風(fēng)機(jī)同時(shí)調(diào)節(jié)風(fēng)向,能有效改善溫度場(chǎng)均勻性分布程度。林加勇等[7]對(duì)公豬舍內(nèi)風(fēng)速場(chǎng)和溫度場(chǎng)環(huán)境分布規(guī)律進(jìn)行模擬研究,結(jié)果表明風(fēng)速場(chǎng)模擬值與實(shí)測(cè)值誤差較小,模擬裝豬時(shí)的豬舍,溫度分布結(jié)構(gòu)與裝豬前相似。高云等[8]對(duì)多環(huán)境參數(shù)控制的豬養(yǎng)殖箱內(nèi)氣流場(chǎng)進(jìn)行了數(shù)值模擬,結(jié)果表明養(yǎng)殖箱內(nèi)部的氣流場(chǎng)形成循環(huán),適宜通風(fēng)且無死角,養(yǎng)殖箱能夠在較長(zhǎng)時(shí)間穩(wěn)定運(yùn)行。大量文獻(xiàn)研究都證明了流場(chǎng)和溫度場(chǎng)對(duì)農(nóng)作物的生長(zhǎng)有著極大影響。目前研究中尚未使用CFD方法用于氣霧栽培種植技術(shù)中,本文在前人研究基礎(chǔ)上使用CFD方法確定農(nóng)作物在氣霧栽培中的適宜風(fēng)速大小及溫度范圍。建立的CFD模型能為氣霧栽培技術(shù)的環(huán)境建模和控制提供理論參考,彌補(bǔ)物理試驗(yàn)的不足。

        為了調(diào)節(jié)蔬菜氣霧栽培箱體氣流及溫度分布的均勻性,合理布置氣霧栽培箱進(jìn)風(fēng)口和出風(fēng)口,避免氣流在氣霧栽培箱內(nèi)部及農(nóng)作物葉片周圍產(chǎn)生繞流,影響農(nóng)作物正常的光合作用和蒸騰作用,該文采用計(jì)算機(jī)流體力學(xué)方法對(duì)氣霧栽培箱內(nèi)部流場(chǎng)和溫度場(chǎng)進(jìn)行數(shù)值模擬,并通過試驗(yàn)驗(yàn)證氣霧栽培箱內(nèi)部流場(chǎng)及溫度場(chǎng)分布規(guī)律。同時(shí)通過合理設(shè)置進(jìn)風(fēng)口速度大小及進(jìn)出風(fēng)口位置,使氣霧栽培箱內(nèi)的氣流速度與溫度兩者達(dá)到適宜及均勻的狀態(tài),讓生菜在氣霧栽培箱內(nèi)的生長(zhǎng)環(huán)境處于最優(yōu)狀態(tài),達(dá)到高質(zhì)高產(chǎn)的目的,以期為氣霧栽培箱的設(shè)計(jì)及優(yōu)化提供理論參考。

        1 材料與方法

        1.1 試驗(yàn)裝置

        本課題組設(shè)計(jì)的蔬菜氣霧栽培箱裝置結(jié)構(gòu)如圖1所示。外部形狀為一個(gè)正六棱,尺寸邊長(zhǎng)0.35 m,高度1.6 m。氣霧栽培箱框架由不銹鋼方鋼焊接而成,頂面、底面及四周用透明亞克力板(PMMA,polymeric methyl methacrylate)進(jìn)行封裝。內(nèi)部栽培架則由6塊擠塑板組成,每一塊擠塑板的尺寸長(zhǎng)1.6 m,寬0.35 m。農(nóng)作物種植在擠塑板(XPS,extruded polystyrene)的定植孔內(nèi)。氣霧栽培箱底部營養(yǎng)液槽里面的營養(yǎng)液通過增壓泵增壓,把營養(yǎng)液通過頂部的高壓霧化噴頭噴射到農(nóng)作物根系,以此提供給農(nóng)作物生長(zhǎng)所需的養(yǎng)分和水分。氣霧栽培箱內(nèi)的送風(fēng)系統(tǒng)采用頂面送回風(fēng)或側(cè)面送回風(fēng)的方式,頂面送回風(fēng)口呈120°角均勻布置3個(gè),尺寸長(zhǎng)為6 cm,寬為6 cm,對(duì)應(yīng)為圖1中頂板上3個(gè)正方形開口;側(cè)面送回風(fēng)口呈120°角均勻布置在3塊亞克力板上,每塊亞克力板上均勻布置3個(gè)送回風(fēng)口,尺寸長(zhǎng)6 cm,寬6 cm,對(duì)應(yīng)為圖1中亞克力板上9個(gè)正方形開口。

        1. 亞克力底板 2.過濾網(wǎng) 3.XPS定植板 4.定植孔 5.防水遮陽膜 6.不銹鋼支撐架 7.亞克力板 8.頂板 9.排風(fēng)扇 10.不銹鋼環(huán)形噴霧圓環(huán) 11.霧化噴頭 12.進(jìn)回風(fēng)口 13.PVC管 14.增壓泵 15.營養(yǎng)液槽

        1.2 試驗(yàn)條件

        試驗(yàn)在湖南農(nóng)業(yè)大學(xué)工學(xué)院智能氣候室內(nèi)進(jìn)行,氣霧栽培箱位于智能氣候室內(nèi)。以智能氣候室為環(huán)境背景倉,對(duì)農(nóng)作物生長(zhǎng)環(huán)境參數(shù)進(jìn)行設(shè)置,模擬出外界環(huán)境真實(shí)值。本文以蔬菜氣霧栽培箱為研究對(duì)象,以生菜為試材。對(duì)于生菜生長(zhǎng)所需溫度,根據(jù)張成波等[9]研究表明生菜發(fā)芽期適宜空氣溫度為18~20 ℃,幼苗期光期空氣溫度為15~20 ℃,采收期光期空氣溫度18~20 ℃。根據(jù)Kitaya,Thongbai等[10-12]研究可知,氣流速度在0.3~1.3 m/s范圍內(nèi)適宜植物生長(zhǎng),植物蒸騰作用及光合作用隨著氣流速度的增大而增強(qiáng),植物周圍空氣流動(dòng)能促進(jìn)植物生長(zhǎng)。以智能氣候室為環(huán)境背景倉,設(shè)置適合生菜生長(zhǎng)所需的溫濕度,氣霧栽培箱外部環(huán)境溫度設(shè)置為18 ℃,相對(duì)濕度控制在85%,光照時(shí)間間隔設(shè)置為12 h。

        1.3 試驗(yàn)方法

        本文以蔬菜氣霧栽培箱為研究對(duì)象,以生菜為試材。對(duì)蔬菜氣霧栽培箱送回風(fēng)口的不同位置布局設(shè)計(jì)出了3種氣流循環(huán)方案,利用FLUENT軟件,對(duì)蔬菜氣霧栽培箱內(nèi)的流場(chǎng)和溫度場(chǎng)進(jìn)行數(shù)值模擬,得到最優(yōu)方案。然后通過試驗(yàn)對(duì)最優(yōu)的方案進(jìn)行驗(yàn)證。試材生菜種植在在氣霧栽培箱內(nèi)部擠塑板上的定植孔內(nèi),在擠塑板上的生菜冠層布置風(fēng)速及溫度監(jiān)測(cè)點(diǎn),每個(gè)監(jiān)測(cè)點(diǎn)距離擠塑板水平距離約5 cm,每塊擠塑板上均勻布置3個(gè)監(jiān)測(cè)點(diǎn),從上往下間隔依次為40 cm,左右等距,6個(gè)面總共布置18個(gè)測(cè)量風(fēng)速值及溫度值監(jiān)測(cè)點(diǎn),如圖2所示。使用長(zhǎng)為8 cm,寬為8 cm的軸流風(fēng)機(jī)調(diào)控風(fēng)速,可布置在12個(gè)進(jìn)回風(fēng)口處,如圖2所示,風(fēng)速可控范圍在0~3.8 m/s。溫度采用溫度傳感器測(cè)量(深圳燦輝電子有限公司生產(chǎn)),布置在溫度監(jiān)測(cè)點(diǎn)處,溫度測(cè)量范圍為-50~110℃,測(cè)試精度:±0.5 ℃,采樣周期2 s。風(fēng)速值測(cè)量采用手持?jǐn)?shù)字風(fēng)速儀進(jìn)行測(cè)量,在風(fēng)速監(jiān)測(cè)點(diǎn)處測(cè)量,分辨率為0.001 m/s,風(fēng)速測(cè)量誤差±3%。

        注:l表示18個(gè)風(fēng)速和溫度值監(jiān)測(cè)點(diǎn)?!瓯硎据S流風(fēng)機(jī)安裝位置。

        2 CFD模型

        2.1 數(shù)值模型

        2.2 生菜作物冠層質(zhì)熱交換模型

        氣霧栽培箱內(nèi)部生長(zhǎng)的生菜把截獲的光照輻射以顯熱和潛熱的方式與周圍空氣進(jìn)行質(zhì)熱交換,根據(jù)作物冠層的能量平衡原理建立數(shù)學(xué)模型[14]:

        2.3 生菜作物多孔介質(zhì)模型

        2.4 定義材料和邊界條件

        對(duì)氣霧栽培箱模型進(jìn)行簡(jiǎn)化,該型氣霧栽培箱為封閉式,因在智能氣候室內(nèi)對(duì)真實(shí)環(huán)境進(jìn)行模擬設(shè)置,所以不存在太陽光直射,光照來源于智能氣候室內(nèi)的植物生長(zhǎng)燈,所以忽略太陽輻射,將氣霧栽培箱6個(gè)亞克力板壁面簡(jiǎn)化為溫度面(wall-temperature);將氣霧栽培箱XPS栽培板設(shè)置為絕熱壁面(wall);進(jìn)風(fēng)口設(shè)置為速度入口(velocity-inlet);出風(fēng)口設(shè)置為壓力出口,壓力設(shè)置為0;進(jìn)出風(fēng)溫度設(shè)置為實(shí)測(cè)值18 ℃[18-20]。模擬的主要參數(shù)和取值如表1所示。

        表1 模擬參數(shù)設(shè)置和取值

        2.5 計(jì)算域和網(wǎng)格劃分

        對(duì)氣霧栽培箱模型進(jìn)行簡(jiǎn)化,建模比例為1∶1,選取氣霧栽培箱葉片生長(zhǎng)區(qū)作為計(jì)算域。在Ansys Workben-ch的ICEM模塊中進(jìn)行網(wǎng)格劃分,由于六面體網(wǎng)格節(jié)點(diǎn)易于計(jì)算,結(jié)構(gòu)化網(wǎng)格生成的網(wǎng)格質(zhì)量高,采用結(jié)構(gòu)化六面體方式對(duì)氣霧栽培箱進(jìn)行網(wǎng)格劃分,利用Fluent軟件對(duì)網(wǎng)格文件進(jìn)行后處理。最后計(jì)算出網(wǎng)格節(jié)點(diǎn)數(shù)為531 628,網(wǎng)格數(shù)量為3 070 692,網(wǎng)格劃分質(zhì)量良好。圖3為蔬菜氣霧栽培箱網(wǎng)格劃分圖。

        圖3 氣霧栽培箱網(wǎng)格劃分

        3 CFD模擬和氣流循環(huán)優(yōu)化設(shè)計(jì)

        3.1 模擬優(yōu)化方案設(shè)計(jì)

        為探究氣霧栽培箱內(nèi)部生菜葉片生長(zhǎng)區(qū)域的風(fēng)速及溫度分布,選取氣霧栽培箱內(nèi)部生菜冠層平面進(jìn)行分析,通過相關(guān)文獻(xiàn)[10-12]可知生菜作物冠層平面適宜風(fēng)速處于0.3~1.30 m/s。

        為了使氣霧栽培箱氣流和溫度滿足生菜作物生長(zhǎng)需要,根據(jù)鄧書輝等[21]、姚家君等[22]、胥芳等[23]、呂曉蘭等[24]學(xué)者對(duì)不同通風(fēng)循環(huán)模式下流場(chǎng)和溫度場(chǎng)的研究,本文設(shè)計(jì)了3種氣流循環(huán)方式,即方案1,頂面進(jìn)側(cè)面出;方案2,側(cè)面進(jìn)頂面出;方案3,側(cè)面進(jìn)側(cè)面出。如圖4所示。利用CFD模擬并比較3種方案中生菜生長(zhǎng)冠層平面風(fēng)速及溫度分布規(guī)律,得出最優(yōu)氣流循環(huán)方案。

        注:紅色平面表示截面A、B和C。

        方案1、方案2、方案3三維模型建立都相同,網(wǎng)格劃分也都一致。將回風(fēng)口設(shè)置為壓力出口邊界條件,將進(jìn)風(fēng)口設(shè)置為速度入口邊界條件。為保證生菜生長(zhǎng)區(qū)域在所需的適宜風(fēng)速范圍內(nèi),對(duì)不同方案的通風(fēng)量,入口風(fēng)速依據(jù)下列公式進(jìn)行計(jì)算[25]:

        3.2 方案模擬仿真分析

        3.2.1 頂面進(jìn)側(cè)面出

        方案1是頂面3個(gè)送風(fēng)口進(jìn)風(fēng),3個(gè)側(cè)面9個(gè)回風(fēng)口出風(fēng)。對(duì)溫度分布圖5a和5b分析可知,從送風(fēng)口處向下呈兩邊回流的趨勢(shì)分布?;亓鳒囟鹊姆植寂c氣流的分布存在一定耦合關(guān)系,風(fēng)速值較大區(qū)域溫度值較小,風(fēng)速值較小區(qū)域溫度值較大。高溫區(qū)域主要集中在進(jìn)口截面的兩側(cè),這是因?yàn)榭諝鈴倪M(jìn)口進(jìn)去后會(huì)在兩側(cè)形成兩個(gè)渦,與主流區(qū)域的能量交換較少,所以溫度相對(duì)較高。作物冠層平面最大溫度為24.55 ℃;溫度平均值為19.05 ℃;適宜溫度所占百分比為93.6%。

        對(duì)進(jìn)口截面氣流速度分布圖5c和5d分析可知,氣流經(jīng)由氣霧栽培箱頂面3個(gè)送風(fēng)口進(jìn)入箱體后,主流區(qū)域直接沖向底部,沖擊底部后再向上流動(dòng)。送風(fēng)口下方氣流變化梯度較大,氣流方向垂直向下,并且主要集中在送風(fēng)口的下部,然后沿著豎直方向向下移動(dòng),在氣霧栽培箱內(nèi)形成漩渦流。部分氣流從氣霧栽培箱側(cè)面回風(fēng)口流出,導(dǎo)致氣霧栽培箱內(nèi)部風(fēng)速減小。對(duì)速度分布流線圖5e分析可知,速度從進(jìn)口進(jìn)入后,主流區(qū)域會(huì)在觸底后向上流動(dòng),并從兩側(cè)的出口流出,整體呈現(xiàn)出對(duì)稱流動(dòng),3個(gè)出口中靠近底部區(qū)域的出口流量比較大,而靠近上側(cè)的出口流量比較小。作物冠層平面最大風(fēng)速為1.43 m/s;內(nèi)部風(fēng)速平均值為0.19 m/s;適宜風(fēng)速所占百分比為58.1%。

        3.2.2 側(cè)面進(jìn)頂面出

        方案2是3個(gè)側(cè)面9個(gè)送風(fēng)口進(jìn)風(fēng),頂面3個(gè)回風(fēng)口出風(fēng)。對(duì)溫度分布圖6a和6b分析可知,送風(fēng)口處溫度較低,回風(fēng)口處溫度稍高,底部區(qū)域的溫度略高于中間區(qū)域,這是因?yàn)閺?個(gè)側(cè)面進(jìn)風(fēng)以后,空氣主要向上流動(dòng),造成了底部區(qū)域的溫度略高于中間區(qū)域,風(fēng)速小的地方溫度較高,風(fēng)速大的地方溫度較低,但整體內(nèi)部溫度分布較為均勻,溫度與氣流的分布存在一定的耦合關(guān)系。作物冠層平面最大溫度為21.68 ℃;最小溫度為17.85 ℃;溫度平均值為18.38 ℃;適宜溫度所占百分比為99.98%。

        圖5 方案1溫度和速度分布云圖及速度分布流線圖

        對(duì)氣流速度分布圖6c和6d分析可知,氣流經(jīng)由氣霧栽培箱3個(gè)側(cè)面9個(gè)送風(fēng)口進(jìn)入箱體后,進(jìn)口速度較高,沖擊內(nèi)壁面后向兩側(cè)流動(dòng),速度分布相對(duì)比較均勻。從進(jìn)口截面上的溫度分布云圖可以看出,溫度分布比較均勻,沒有明顯的高溫與低溫區(qū)域。兩個(gè)進(jìn)口面中間的區(qū)域溫度略高于進(jìn)口面進(jìn)風(fēng)口上下方向。氣流方向朝遠(yuǎn)離送風(fēng)口的方向速度逐漸變小。氣流從氣霧栽培箱頂面回風(fēng)口流出時(shí)氣流速度有變大的趨勢(shì)。對(duì)速度分布流線圖6e分析可知,流動(dòng)整體分布比較均勻,沒有明顯的高速與低速。底部區(qū)域的流動(dòng)比較小,這是因?yàn)槌隹谠谏蟼?cè),流體進(jìn)去后,主要向上流動(dòng)。作物冠層平面最大風(fēng)速為1.62 m/s;內(nèi)部風(fēng)速平均值為0.19 m/s;適宜風(fēng)速所占百分比為59.6%。

        圖6 方案2溫度和速度分布云圖及速度分布流線圖

        3.2.3 側(cè)面進(jìn)側(cè)面出

        方案3是1個(gè)側(cè)面3個(gè)送風(fēng)口進(jìn)風(fēng),2個(gè)側(cè)面6個(gè)回風(fēng)口出風(fēng)。對(duì)溫度分布圖7a和7b分析可知,送風(fēng)口處溫度向前逐漸變高,兩邊逐漸降低?;仫L(fēng)口處溫度稍高,靠近進(jìn)口端面的溫度比較低,出口端面的溫度比較高??拷隹诮孛嫔系臏囟缺容^高。溫度整體呈現(xiàn)出明顯的左右不均勻性。作物冠層平面最大溫度為25.95 ℃;最小溫度為17.65 ℃;溫度平均值為19.2 ℃;適宜溫度所占百分比為92.4%。

        對(duì)氣流速度分布圖7c和7d分析可知,進(jìn)風(fēng)口上下方向氣流變化梯度較大,進(jìn)口面上的速度比較高,其余5個(gè)面域的速度比較低,速度分布的均勻性比較差。送風(fēng)口出氣流速度較大,氣流方向朝遠(yuǎn)離送風(fēng)口的方向速度向四周逐漸變小,3個(gè)送風(fēng)口的氣流分布上下不均勻,存在回風(fēng)口的2個(gè)側(cè)面氣流分布梯度較小,氣流從回風(fēng)口流出時(shí)氣流速度有變大的趨勢(shì),靠近出口的截面上整體速度比較低。對(duì)速度分布流線圖7e分析可知,進(jìn)口段流線分布比較密集,該段的流速較高,進(jìn)去后從兩側(cè)出口流出,進(jìn)口端正對(duì)面的區(qū)域中,流動(dòng)情況比較差。作物冠層平面最大風(fēng)速為1.41 m/s;內(nèi)部風(fēng)速平均值為0.08 m/s;適宜風(fēng)速所占百分比為54.3%。

        圖7 方案3溫度和速度分布云圖及速度分布流線圖

        3.2.4 方案優(yōu)化討論

        對(duì)3個(gè)方案溫度和風(fēng)速分布云圖分析及表2、3可知。風(fēng)速送風(fēng)口處風(fēng)速一般比較大,溫度較低,所以溫度的分布與氣流的分布存在一定耦合關(guān)系,風(fēng)速較大區(qū)域溫度較小。3個(gè)方案風(fēng)速最小值都為0,方案1風(fēng)速最大值為1.43 m/s,風(fēng)速平均值為0.19 m/s。風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值0.3~1.3 m/s的區(qū)域占58.1%,溫度平均值為19.05 ℃,溫度位于生菜生長(zhǎng)適宜溫度值15~20 ℃的區(qū)域占93.6%。方案2風(fēng)速最大值為1.62 m/s,風(fēng)速平均值為0.19 m/s。風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值0.3~1.3 m/s的區(qū)域占59.6%,溫度平均值為18.38 ℃,溫度位于生菜生長(zhǎng)適宜溫度值15~20 ℃的區(qū)域占99.98%。方案3風(fēng)速最大值為1.41 m/s,風(fēng)速平均值為0.08 m/s。風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值0.3~1.3 m/s的區(qū)域占54.3%,溫度平均值為19.2℃,溫度位于生菜生長(zhǎng)適宜溫度值15~20℃的區(qū)域占92.4%。綜合比較3種方案,方案2為最優(yōu)方案。

        表2 不同氣流循環(huán)方案作物冠層平面風(fēng)速

        表3 不同氣流循環(huán)方案作物冠層平面溫度

        3.2.5 方案通風(fēng)死角討論

        氣霧栽培箱內(nèi)部進(jìn)行農(nóng)作物栽培時(shí),由于部分區(qū)域通風(fēng)不均勻會(huì)產(chǎn)生死角,導(dǎo)致沒有氣流速度分布,影響溫度的分布均勻性,會(huì)影響農(nóng)作物的正常生長(zhǎng)。所以對(duì)進(jìn)出風(fēng)口的設(shè)計(jì)布局顯得尤為重要,包括進(jìn)出風(fēng)口大小、數(shù)量、布局等都會(huì)影響氣流的均勻流動(dòng)。設(shè)計(jì)合理的進(jìn)出風(fēng)口可避免農(nóng)作物生長(zhǎng)過程中引起的一系列生產(chǎn)問題。對(duì)文中3種方案的死角問題進(jìn)行討論,通過數(shù)值模型截面分析可知(如表4和圖8所示)。

        注:紅框表示通風(fēng)死角區(qū)域

        方案1中,通風(fēng)死角區(qū)域占比0.844%,主要分布在通風(fēng)進(jìn)口面和通風(fēng)出口面之間的交界處靠近上側(cè)的部分,這是因?yàn)轱L(fēng)從上側(cè)進(jìn)入后,到達(dá)底部區(qū)域后方向旋轉(zhuǎn),向上流動(dòng),從出口面的3個(gè)出口流出,越向上風(fēng)流量越小,所以在上側(cè)區(qū)域會(huì)出現(xiàn)部分死角。方案2中,通風(fēng)死角區(qū)域占比僅0.069%,沒有明顯的通風(fēng)死角,主要原因是進(jìn)出風(fēng)口布局較為合理。方案3中,通風(fēng)死角區(qū)域較大,占16.7%,主要分布在2個(gè)出口所在面的中間區(qū)域,這是因?yàn)轱L(fēng)從進(jìn)口進(jìn)去后,直接從兩側(cè)流向出口,只有極少空氣流向兩個(gè)出口的中間區(qū)域。通過對(duì)3中方案進(jìn)行對(duì)比可知,方案2的通風(fēng)性較好,通風(fēng)死角區(qū)域占比最小。

        表4 不同方案通風(fēng)死角所占百分比及溫度平均值

        3.3 試驗(yàn)結(jié)果分析與CFD模型驗(yàn)證

        為驗(yàn)證仿真模型的可靠性和正確性,在蔬菜氣霧栽培箱內(nèi)種植生菜的情況下,對(duì)最優(yōu)方案2各點(diǎn)溫度和風(fēng)速的實(shí)測(cè)值與模擬值進(jìn)行試驗(yàn)驗(yàn)證,如圖9所示。對(duì)圖9a的18個(gè)溫度監(jiān)測(cè)點(diǎn)的模擬值與實(shí)測(cè)值進(jìn)行對(duì)比分析可知,測(cè)點(diǎn)11的誤差值最大,最大誤差為1.35 ℃,相對(duì)誤差最大6.9%,平均相對(duì)誤差為3.9%,均方根誤差為0.86 ℃。對(duì)圖9b的18個(gè)風(fēng)速監(jiān)測(cè)點(diǎn)的模擬值與實(shí)測(cè)值進(jìn)行對(duì)比分析可知,測(cè)點(diǎn)6的絕對(duì)誤差最大,最大誤差為0.54 m/s,相對(duì)誤差最大9.5%,平均相對(duì)誤差為3.5%,均方根誤差為0.26 m/s。

        從圖中對(duì)數(shù)據(jù)比較可以發(fā)現(xiàn),實(shí)測(cè)值與模擬值存在偏差,但整體的溫度和風(fēng)速分布較為均勻,所建立的CFD模型是準(zhǔn)確的,可用此模型對(duì)其他方案進(jìn)行驗(yàn)證分析。

        4 結(jié) 論

        1)本文中設(shè)計(jì)了3種氣流循環(huán)方式:方案1為頂面進(jìn)側(cè)面、方案2為側(cè)面進(jìn)頂面出、方案3為側(cè)面進(jìn)側(cè)面出。對(duì)3個(gè)方案溫度和風(fēng)速分布云圖分析后得到風(fēng)速送風(fēng)口處風(fēng)速一般比較大,溫度較低,風(fēng)速較大區(qū)域溫度較小,溫度的分布和氣流的分布存在一定耦合關(guān)系。方案1風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值區(qū)域占58.1%,適宜溫度值區(qū)域占93.6%,通風(fēng)死角區(qū)域占比0.844%。方案2風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值區(qū)域占59.6%,適宜溫度值區(qū)域占99.98%,通風(fēng)死角區(qū)域占比0.069%。方案3風(fēng)速位于生菜生長(zhǎng)適宜風(fēng)速值區(qū)域占54.3%,,適宜溫度值區(qū)域占92.4%,通風(fēng)死角區(qū)域占比16.7%。綜合比較3種方案并結(jié)合生菜生長(zhǎng)的溫度和風(fēng)速,方案2為最優(yōu)方案。

        2)將18個(gè)點(diǎn)的模擬值與實(shí)測(cè)值進(jìn)行對(duì)比,溫度平均相對(duì)誤差為3.9%,均方根誤差為0.86 ℃。風(fēng)速平均相對(duì)誤差為3.5%,均方根誤差為0.26 m/s。模擬值和實(shí)測(cè)值誤差較小,模擬效果良好。建立的CFD模型可準(zhǔn)確的描述蔬菜氣霧栽培箱內(nèi)部的流場(chǎng)和溫度場(chǎng)變化規(guī)律和空間分布特性。

        圖9 氣霧栽培箱溫度及風(fēng)速監(jiān)測(cè)點(diǎn)模擬值與實(shí)測(cè)值比較

        3)本文中模擬了3種氣流循環(huán)方案,得到側(cè)面進(jìn)頂面為最佳氣流循環(huán)模式,并對(duì)送回風(fēng)口的不同位置布局進(jìn)行了研究。如果繼續(xù)進(jìn)行優(yōu)化研究,可結(jié)合正交試驗(yàn)方法,將送回風(fēng)口大小、風(fēng)速、數(shù)量及位置等設(shè)置為不同因素,對(duì)該類蔬菜氣霧栽培箱的送回風(fēng)口的通風(fēng)進(jìn)行優(yōu)化設(shè)計(jì),避免通風(fēng)死角區(qū)域,模擬并設(shè)計(jì)出最佳通風(fēng)模式。

        4)所建立的CFD模型,模擬預(yù)測(cè)了不同氣流循環(huán)方式下蔬菜氣霧栽培箱內(nèi)部溫度及風(fēng)速的變化規(guī)律。本文對(duì)單個(gè)的蔬菜氣霧栽培箱流場(chǎng)和溫度場(chǎng)進(jìn)行了模擬,可對(duì)大規(guī)模商業(yè)化種植,包括在溫室大棚內(nèi)氣霧栽培設(shè)備布局、大型霧化工廠種植、家庭陽臺(tái)霧化裝置等,氣霧栽培技術(shù)的環(huán)境建模和控制提供更多理論依據(jù)和信息,彌補(bǔ)物理試驗(yàn)的不足。

        [1] 張建國,何春梅,凌敏,等. 氣霧栽培的應(yīng)用與研究綜述[J]. 林業(yè)與環(huán)境科學(xué),2017,33(4):130-134.

        Zhang Jianguo, He Chunmei, Ling Min, et al. Review on the application and research progress of aeroponic cultivation[J]. Forestry and Environmental Science, 2017, 33(4): 130-134. (in Chinese with English abstract)

        [2] 高建民,黃桂珍,尹文楚,等. 桁架式超聲霧化栽培器的霧滴沉降和根際溫濕度變化規(guī)律[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(6): 185-192. Gao Jianmin, Huang Guizhen, Yin Wenchu, et al. Variation of root-zone droplet deposition, temperatureand humidity of truss ultrasonic atomization cultivator[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(6): 185-192. (in Chinese with English abstract)

        [3] 劉義飛. 基于LabVIEW的溫室番茄霧培控制系統(tǒng)研究[D]. 北京:中國農(nóng)業(yè)科學(xué)院,2014. Liu Yifei. Aeroponics Control System for Tomato Cultivation in Greenhouse Based on LabVIEW[D]. Beijing: Chinese Academy of Agricultural Sciences, 2014. (in Chinese with English abstract)

        [4] 劉巖. 植物工廠中霧培成套裝置的設(shè)計(jì)與應(yīng)用[D]. 杭州:浙江大學(xué),2014. Liu Yan. Design and Application of Aeroponics Package Unit in Plant Factory[D]. Hangzhou: Zhejiang University, 2014. (in Chinese with English abstract)

        [5] 王新忠,張偉建,張良,等. 基于CFD的夏季屋頂全開型玻璃溫室自然通風(fēng)流場(chǎng)分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào),2016,47(10): 332-337. Wang Xinzhong, Zhang Weijian, Zhang Liang, et al. Analysis of flow field for full open-roof glass greenhouse with nature ventilation in summer based on CFD[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(10): 332-337. (in Chinese with English abstract)

        [6] 翁衛(wèi)兵,房殿軍,李強(qiáng),等. 冷藏運(yùn)輸廂內(nèi)流場(chǎng)和溫度場(chǎng)協(xié)同控制[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014, 45(6): 260-265. Weng Weibing, Fang Dianjun, Li Qiang, et al. Cooperative control of flow field and temperature field in refrigerated transport carriage[J]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(6): 260-265. (in Chinese with English abstract)

        [7] 林加勇,劉繼軍,孟慶利,等. 公豬舍夏季溫度和流場(chǎng)數(shù)值CFD模擬及驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào),2016,32(23): 207-212. Lin Jiayong, Liu Jijun, Meng Qingli, et al. Numerical CFD simulation and verification of summer indoor temperature and airflow field in boar building[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(23): 207-212. (in Chinese with English abstract)

        [8] 高云,陳震撼,王瑜,等. 多環(huán)境參數(shù)控制的豬養(yǎng)殖箱設(shè)計(jì)及箱內(nèi)氣流場(chǎng)分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(2):203-212. Gao Yun, Chen Zhenhan, Wang Yu, et al. Design for pig breeding chamber under multiple environment variable control and analysis of internal flow field[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(2): 203-212. (in Chinese with English abstract)

        [9] 張成波,楊其長(zhǎng). 植物工廠概論[M].北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社,2005:56-156.

        [10] Kitaya Y, Tsuruyama J, Shibuya T, et al. Effects of air current speed on gas exchange in plant leaves and plant canopies[J]. Advances in Space Research the Official Journal of the Committee on Space Research, 2003, 31(1): 177-182.

        [11] Thongbai P, Kozai T, Ohyama K. CO2and air circulation effects on photosynthesis and transpiration of tomato seedlings[J].Scientia Horticulturae, 2010, 126(3): 338-344.

        [12] Kitaya Y. Effects of air velocity on photosynthesis of plant canopies under elevated CO2levels in a plant culture system[J]. Advances in Space Research the Official Journal of the Committee on Space Research, 2004, 34(7): 1466-1469.

        [13] Versteeg H K, Malalasekera W. An Introduction to Computational Fluid Dynamics[M]. London: Longman Group Ltd. , 1995.

        [14] Boulard T, Wang S. Experimental and numerical studies on the heterogeneity of crop transpiration in a plastic tunnel[J]. Computers and Electronics in Agriculture, 2002, 34(1/2/3): 173-190.

        [15] Bruse M. Development of a Numerical Model for the Simulation of Exchange Processes Between Small Scale Environmental Design and Microclimate in Urban Areas[D]. Bochum: University of Bochum, 1998.

        [16] Bruse M.Development of a Micro-scale Model for the Calculation of Surface Temperature in Structured Terrain[D]. Bochum: University of Bochum,1995.

        [17] Kichah A, Bournet P E, Migeon C, et al. Measurement and CFD simulation of microclimate characteristics and transpiration of an impatiens pot plant crop in a greenhouse[J]. Biosystems Engineering, 2012, 112(3): 252 -252.

        [18] Wilson J D. Numerical studies of flow through a windbreak[J]. Journal of Wind Engineering & Industrial Aerodynamics, 1985, 21(2): 119-154.

        [19] 王福軍. 計(jì)算流體動(dòng)力學(xué)分析[M].北京:清華大學(xué)出版社,2004: 1-57.

        [20] 張師帥. 計(jì)算流體動(dòng)力學(xué)及其應(yīng)用[M].武漢:華中科技大學(xué)出版社,2011:1-110.

        [21] 鄧書輝,施正香,李保明,等. 擋風(fēng)板對(duì)低屋面橫向通風(fēng)牛舍內(nèi)空氣流場(chǎng)影響的PIV測(cè)試[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(1):188-194. Deng Shuhui, Shi Zhengxiang, Li Baoming, et al. Influence of baffle on air flow field in low profile cross ventilated cattle barn analysis based on PIV test[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(1): 188-194. (in Chinese with English abstract)

        [22] 姚家君,郭彬彬,丁為民,等. 基于鵝舍氣流場(chǎng)CFD模擬的通風(fēng)系統(tǒng)結(jié)構(gòu)優(yōu)化與驗(yàn)證[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(3):214-220. Yao Jiajun, Guo Binbin, Ding Weimin, et al. Structure optimization and validation of goose house ventilation system based on airflow field simulation by CFD[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(3): 214- 220. (in Chinese with English abstract)

        [23] 胥芳,蔡彥文,陳教料,等. 濕簾-風(fēng)機(jī)降溫下的溫室熱/流場(chǎng)模擬及降溫系統(tǒng)參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(9):201-208. Xu Fang, Cai Yanwen, Chen Jiaoliao, et al. Temperature/ flow field simulation and parameter optimal design for greenhouses with fan-pad evaporative cooling system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(9): 201-208. (in Chinese with English abstract)

        [24] 呂曉蘭,張美娜,常有宏,等. 果園風(fēng)送噴霧機(jī)導(dǎo)流板角度對(duì)氣流場(chǎng)三維分布的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(15):81-87. Lü Xiaolan, Zhang Meina, Chang Youhong, et al. Influence of deflector angles for orchard air-assisted sprayer on 3D airflow distribution[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(15): 81-87. (in Chinese with English abstract)

        [25] 蔡增基,龍?zhí)煊? 流體力學(xué)泵與風(fēng)機(jī)[M].北京:中國建筑工業(yè)出版社,2009:156-179.

        CFD simulation of flow field and temperature field of different air flow circulation modes in vegetable aeroponics cultivation box

        Jiang Pin1,2, Yang Xiwen1, Luo Yahui1,2, Long Lixia1, Shi Yixin1,2, Hu Wenwu1,2

        (1.410128,; 2.410128,)

        The air flow and temperature in the vegetable aeroponics cultivation box play an important role in the internal environment regulation and crop growth. The uniform air flow around crops can promote the growth rate of the crop. In order to study the distribution of flow field and temperature field in vegetable aeroponics cultivation box, CFD models of vegetable aeroponics cultivation box under different air flow circulation modes were established based on computational fluid dynamic (CFD) method, FLUENT software, standard turbulence model, porous medium model, crop canopy heat exchange model. The environment in the aeroponics cultivation box was optimized, and three air flow circulation schemes were designed(air flows from top to side, side to top and side to side). The layout of the different positions of the return air outlets was studied, and the numerical simulations of the three airflow circulation schemes were carried out. The simulation results showed that in the scheme of air flow entering from the top and exiting from the side, the suitable wind speed value area of lettuce accounts for 58.1%, the suitable temperature value area accounts for 93.6%, and the ventilation dead zone area accounts for 0.844%.In the scheme of air flow side entering and exiting the top surface, the suitable wind speed value area of lettuce accounts for 59.6%, the suitable temperature value area accounts for 99.98%, and the ventilation dead zone area accounts for 0.069%. In the scheme of air flow side entering side exit, the suitable wind speed value area for lettuce growth accounts for 54.3%, the suitable temperature value area accounts for 92.4%, and the ventilation dead zone area accounts for 16.7%. After comparison, the side entering the top surface was obtained as the optimal air circulation scheme. The temperature and wind speed simulation values and measured values in the aeroponics cultivation box were compared. The results showed that the average temperature relative error was 3.9% and the root mean square error was 0.86 ℃. The average relative error of wind speed was 3.5%, and the root mean square error was 0.26 m/s. The error between simulated and measured values was small, and the simulation effect was good. The accuracy of the CFD model was verified. In crop cultivation, aeroponics was a new cultivation method in the field of soilless cultivation, which could be applied to greenhouses, plant factories, family balconies and other places through microclimate control technology. In this paper, by regulating the wind speed in the aeroponics box, the climate growth environment and temperature uniformity of the crops can be ensured, and the problems of the roots and canopies of the crops can be effectively solved, and the crop yield and quality can be improved. Airflow velocity and temperature are important factors for crop growth in aeroponics boxes, which affect the transpiration and photosynthesis of crops. Therefore, it is of great significance to study the changes of air velocity and temperature to promote air aeroponics technology. This study provides a reference for the flow field and temperature change law in vegetable aeroponics cultivation box, internal environment adjustment and device optimization design.

        wind speed; temperature; computational fluid dynamics; optimization design; vegetable aeroponics cultivation box

        2019-02-20

        2019-07-16

        國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFD0700903-2);湖南省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2018NK2063)

        蔣 蘋,教授,博士,博士生導(dǎo)師,主要從事農(nóng)業(yè)機(jī)械智能控制研究。Email:teacher_jp@163.com

        10.11975/j.issn.1002-6819.2019.16.026

        S317

        A

        1002-6819(2019)-16-0233-09

        蔣 蘋,楊希文,羅亞輝,龍莉霞,石毅新,胡文武.蔬菜氣霧栽培箱不同氣流循環(huán)方式的流場(chǎng)和溫度場(chǎng)CFD模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(16):233-241. doi:10.11975/j.issn.1002-6819.2019.16.026 http://www.tcsae.org

        Jiang Pin, Yang Xiwen, Luo Yahui, Long Lixia, Shi Yixin, Hu Wenwu. CFD simulation of flow field and temperature field of different air flow circulation modes in vegetable aeroponics cultivation box[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(16): 233-241. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.16.026 http://www.tcsae.org

        猜你喜歡
        氣霧生菜氣流
        脆嫩爽口的生菜
        中老年保健(2022年6期)2022-08-19 01:44:10
        氣流的威力
        廢棄氣霧罐回收裝置的機(jī)械結(jié)構(gòu)設(shè)計(jì)
        生菜怎么吃更健康
        生菜?你愛吃圓的還是散葉兒的?
        食品與健康(2017年3期)2017-03-15 18:07:48
        生菜有道
        美食堂(2017年2期)2017-02-14 00:16:54
        工廠化蔬菜氣霧栽培技術(shù)
        固體運(yùn)載火箭變軌發(fā)動(dòng)機(jī)噴管氣流分離研究
        飛片下的空氣形成的“超強(qiáng)高速氣流刀”
        基于停留時(shí)間分布的氣流床氣化爐通用網(wǎng)絡(luò)模型
        中文字幕人妻互换av| 荡女精品导航| 日韩欧美在线观看成人| 免费人成年小说在线观看| 少妇内射视频播放舔大片| 国产成人精品免费视频大全| 亚洲黄片av在线免费观看| 久久久精品人妻一区二区三区游戏 | 色诱视频在线观看| 国内精品视频在线播放不卡 | 一个少妇的淫片免费看| 中文字幕中文有码在线| 日本不卡在线视频二区三区| 国产性一交一乱一伦一色一情| 日日噜噜夜夜久久密挑| 人妻少妇69久久中文字幕| 曰韩少妇内射免费播放| 成年美女黄的视频网站| 久久国内精品自在自线图片| 国产专区国产av| 精品免费一区二区三区在| 国产精品二区三区在线观看| 亚洲一区二区三区日本久久九| 又粗又黄又猛又爽大片免费| a级毛片成人网站免费看 | 亚洲中文字幕无码专区| 欧美在线观看www| 精品人妻av中文字幕乱| 精品香蕉一区二区三区| 中文字幕精品一区二区2021年| 久久综合视频网站| 亚洲性爱区免费视频一区| 中文字幕高清视频婷婷| 国产在线一区二区三区乱码| 夹得好湿真拔不出来了动态图| 中文字幕精品一二三四五六七八| 亚洲色大成网站www尤物| 国产91在线精品福利| 亚洲国产精品嫩草影院久久av| 一本色道久在线综合色| 97久久精品人妻人人搡人人玩|