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

    高頻熱聲發(fā)動(dòng)機(jī)的聲耦合特性

    2017-10-12 08:44:05劉益才莫雙林雷斌義夏淯博馬穎
    關(guān)鍵詞:熱聲諧振腔聲壓

    劉益才,莫雙林,雷斌義,夏淯博,馬穎

    ?

    高頻熱聲發(fā)動(dòng)機(jī)的聲耦合特性

    劉益才,莫雙林,雷斌義,夏淯博,馬穎

    (中南大學(xué)能源科學(xué)與工程學(xué)院,湖南長沙,410083)

    將協(xié)同理論引入高頻熱聲發(fā)動(dòng)機(jī)的耦合研究中,基于協(xié)同理論得出對(duì)高頻熱聲發(fā)動(dòng)機(jī)耦合的重要影響參數(shù);建立2列熱聲振蕩器在同一個(gè)諧振腔中聲耦合的數(shù)值研究模型,分析高頻熱聲發(fā)動(dòng)機(jī)初始相位差及間距、諧振腔長度以及諧振管長度對(duì)聲耦合作用的影響。研究結(jié)果表明:對(duì)高頻熱聲發(fā)動(dòng)機(jī)耦合有重要影響的參數(shù)有熱聲發(fā)動(dòng)機(jī)間的初始相位差、諧振頻率以及耦合強(qiáng)度,高頻熱聲發(fā)動(dòng)機(jī)間距對(duì)協(xié)同影響較??;在保持模型直徑、管壁粗糙度等尺寸結(jié)構(gòu)不變的條件下,當(dāng)諧振腔長度為6 cm時(shí)出現(xiàn)協(xié)同現(xiàn)象;當(dāng)諧振腔長度為6 cm時(shí),改變諧振管的長度,諧振頻率差對(duì)耦合的影響很大,協(xié)同狀態(tài)消失并出現(xiàn)拍頻現(xiàn)象,拍頻與2列熱聲振蕩器的諧振頻率差近似呈線性關(guān)系。

    高頻熱聲發(fā)動(dòng)機(jī);聲耦合;數(shù)值仿真

    近30年里,熱聲熱機(jī)的理論研究與實(shí)際應(yīng)用都取得了很大的發(fā)展[1?2]。常規(guī)的熱聲發(fā)動(dòng)機(jī)因其工作頻率在300 Hz以內(nèi)[3?5],使其具有較大的尺寸結(jié)構(gòu),且聲功密度較低,起振溫差較大,限制了熱聲發(fā)動(dòng)機(jī)在低品位熱源領(lǐng)域的應(yīng)用。近年來,國內(nèi)外眾多學(xué)者采用多級(jí)串聯(lián)環(huán)形熱聲發(fā)動(dòng)機(jī)的方式來降低熱聲發(fā)動(dòng)機(jī)的起振條件及提高聲功密度和輸出功率[6?8]。但隨著熱聲發(fā)動(dòng)機(jī)級(jí)數(shù)的增加,整機(jī)尺寸增大,在一定程度上限制了環(huán)形熱聲發(fā)動(dòng)機(jī)的發(fā)展。隨著人們對(duì)低品位能源的重視,高頻微型化熱聲發(fā)動(dòng)機(jī)成為了重要的研究方向[9?11]。高頻熱聲發(fā)動(dòng)機(jī)的諧振頻率范圍為2~25 kHz。高頻熱聲發(fā)動(dòng)機(jī)的結(jié)構(gòu)大小與諧振頻率有關(guān),諧振頻率越大,熱機(jī)尺寸越小。盡管總聲功隨著熱聲發(fā)動(dòng)機(jī)尺寸的減小而降低,但是單位聲能密度隨著熱聲發(fā)動(dòng)機(jī)尺寸的減小而增大,且微型熱聲發(fā)動(dòng)機(jī)需要的起振溫差更小,故如何將眾多高頻微型熱聲發(fā)動(dòng)機(jī)整合為一個(gè)整體,實(shí)現(xiàn)低起振溫度以及高單位聲功輸出是高頻微型化熱聲發(fā)動(dòng)機(jī)的重點(diǎn)研究目標(biāo)。微型熱聲發(fā)動(dòng)機(jī)陣列式排布是解決該問題的一個(gè)有效途徑。當(dāng)多個(gè)熱聲發(fā)動(dòng)機(jī)處于同一腔體內(nèi)陣列式排布時(shí),它們產(chǎn)生的聲振蕩會(huì)發(fā)生相互作用,各個(gè)熱聲發(fā)動(dòng)機(jī)的聲振蕩頻率和相位角會(huì)發(fā)生改變,當(dāng)這種相互作用效果足夠強(qiáng)時(shí),會(huì)使各個(gè)熱聲發(fā)動(dòng)機(jī)的聲振蕩特性趨于穩(wěn)定,達(dá)到一種同相或反相協(xié)同的狀態(tài)[12]。而對(duì)協(xié)同特性進(jìn)行研究有利于理解高頻熱聲發(fā)動(dòng)機(jī)聲耦合特性。

    1 振蕩耦合理論

    熱聲發(fā)動(dòng)機(jī)是一種聲振蕩器。振蕩器的耦合就是振蕩器之間通過第3種介質(zhì)進(jìn)行相互作用,并且對(duì)彼此產(chǎn)生影響。對(duì)于線性振蕩器來說,耦合是線性的。對(duì)于熱聲發(fā)動(dòng)機(jī)這種非線性的振蕩器來說,其耦合機(jī)理很復(fù)雜,很難用準(zhǔn)確的數(shù)學(xué)模型進(jìn)行描述。耦合是熱聲振蕩器之間相互影響的方式,然而,不是所有的耦合作用都會(huì)使振蕩器達(dá)到一種穩(wěn)定的協(xié)同狀態(tài)。只有振蕩器之間的作用強(qiáng)度即耦合強(qiáng)度足夠大時(shí),才會(huì)使振蕩器產(chǎn)生一種穩(wěn)定的協(xié)同狀態(tài)。

    在全局耦合條件下,2列熱聲振蕩器在相空間上的數(shù)學(xué)描述可以表示為:

    (2)

    圖1 二維相空間的極限環(huán)

    極限環(huán)上的相位角可以表示為:

    假設(shè)耦合強(qiáng)度很小,自變量只有相位角,則式(2)可以改寫為:

    (4)

    假設(shè)2個(gè)振蕩器滿足諧振條件:

    (6)

    其中:和為沒有公約數(shù)的整數(shù)。將傅里葉級(jí)數(shù)項(xiàng)和代入式(4)可得:

    其中:

    (8)

    從式(9)可知:當(dāng)2個(gè)振蕩器通過第3種介質(zhì)相互作用時(shí),會(huì)產(chǎn)生耦合效應(yīng),從相角變化率的角度看,會(huì)改變每個(gè)振蕩器的相位角。然而,這種改變與諧振器各自的固有頻率、初始相位角以及耦合強(qiáng)度有關(guān)。若耦合強(qiáng)度較弱,則耦合作用不穩(wěn)定,相位差隨著初始相位角和固有頻率的改變而改變;若耦合強(qiáng)度較強(qiáng),則會(huì)使振蕩器之間的相位差到達(dá)一種穩(wěn)定狀態(tài),稱為協(xié)同;當(dāng)協(xié)同相位差接近0°時(shí),稱為正向協(xié)同;當(dāng)協(xié)同相位差接近180°時(shí),稱為反向協(xié)同。

    2 熱聲振蕩器耦合的數(shù)值模擬

    由振蕩耦合理論可知,2列熱聲振蕩器的耦合作用主要受到熱聲的振蕩特性和諧振腔結(jié)構(gòu)特性的影響。熱聲自激振蕩其實(shí)質(zhì)是一種聲波振蕩,而聲波的主要特性是壓力的波動(dòng)。振蕩頻率、相位角和振蕩幅值為壓力波動(dòng)的重要參數(shù),這3個(gè)參數(shù)也是熱聲振蕩特性的具體體現(xiàn)。對(duì)于圓柱型諧振腔來說,諧振腔的結(jié)構(gòu)特性主要體現(xiàn)為諧振腔的長度、半徑以及2列熱聲振蕩器的擺放位置。為了研究熱聲振蕩特性和諧振腔結(jié)構(gòu)特性對(duì)2列熱聲振蕩器耦合作用的影響,利用Fluent軟件對(duì)2列高頻熱聲振蕩器的耦合特性進(jìn)行數(shù)值研究。其模型簡圖如圖2所示。圖2中,點(diǎn)1,2和3為所取的3個(gè)監(jiān)測點(diǎn),用于監(jiān)測耦合作用對(duì)2列聲波及其合作用的影響。監(jiān)測點(diǎn)1和2分別距離諧振管開口端距離0.5 cm左右,監(jiān)測點(diǎn)3位于諧振腔端面附近。分別用A和B代表2列熱聲振蕩器。

    圖2 熱聲振蕩器耦合模型

    2.1 初始相位角對(duì)熱聲振蕩器耦合的影響

    近年來,人們對(duì)高頻熱聲發(fā)動(dòng)機(jī)的耦合研究較 多[13?16]。BRENNA等[17]以2個(gè)高頻微型熱聲熱機(jī)為例,對(duì)兩者的聲耦合進(jìn)行了實(shí)驗(yàn)研究。然而,對(duì)于熱聲發(fā)動(dòng)機(jī)來說,其初始相位角是隨機(jī)的,實(shí)驗(yàn)中無法準(zhǔn)確測量2列高頻熱聲發(fā)動(dòng)機(jī)的初始相位差,進(jìn)而無法對(duì)2列高頻熱聲發(fā)動(dòng)機(jī)在諧振腔中的耦合特性進(jìn)行準(zhǔn)確描述。為了研究熱聲振蕩初始相位差對(duì)熱聲耦合作用的影響,在熱聲板疊處施加壓力波的自定義函數(shù)(UDF),以替代熱聲自激振蕩過程。對(duì)于這2列壓力波動(dòng),設(shè)定頻率和幅值相等,初始相位差為別為45°,90°,145°和180°。對(duì)于固定結(jié)構(gòu)的諧振腔來說,其耦合強(qiáng)度是一定的。首先研究一定耦合強(qiáng)度下初始相位差對(duì)2列熱聲振蕩器耦合作用的影響。該模型中諧振管的長度為3 cm,半徑為0.4 cm;諧振腔的長度為 4 cm,半徑為4 cm,2列熱聲振蕩器A和B的間距為2 cm。初始相位差對(duì)耦合后2列熱聲振蕩器相位差及合成聲壓幅值的影響如圖3所示。

    從圖3可以看出:初始相位差不同,耦合后的相位差也不同,耦合后的相位差隨著初始相位差的增大而增大,并且逐漸趨于180°,而2列振蕩器合作用的聲壓幅值由于這種反向作用的增大而逐漸減小。在諧振腔結(jié)構(gòu)尺寸不變也就是耦合強(qiáng)度不變條件下,改變熱聲振蕩的初始相位差,對(duì)熱聲振蕩器的耦合產(chǎn)生不同的影響。耦合后的相位差隨著初始相位角的改變而改變,這說明在諧振腔的長度為4 cm,半徑為4 cm,2列諧振腔的間距為2 cm的條件下,耦合強(qiáng)度較弱,2列熱聲振蕩器雖然發(fā)生了耦合作用,但沒有達(dá)到一種穩(wěn)定的耦合相位差,也就是沒有發(fā)現(xiàn)“鎖相”現(xiàn)象,沒有產(chǎn)生“協(xié)同”作用。

    1—耦合后相位差;2—合成聲壓幅值。

    2.2 熱聲振蕩器距離對(duì)聲耦合的影響

    仍不改變諧振腔的尺寸結(jié)構(gòu),在諧振腔長度為 4 cm、半徑為4 cm的條件下,改變2列熱聲振蕩器在諧振腔中的相對(duì)位置,在耦合強(qiáng)度不變時(shí)觀察熱聲振蕩器的布置位置對(duì)耦合作用的影響。所采用的模型如圖2所示。2列熱聲振蕩器之間的間距分別為2,3和4 cm。設(shè)定這2列熱聲振蕩器的壓力波動(dòng)頻率和幅值相等,初始相位差為別為45°,90°,145°和180°。耦合后聲壓幅值和耦合后相位差隨熱聲振蕩器間距的關(guān)系分別如圖4和圖5所示。

    從圖4可知:當(dāng)2列熱聲振蕩器的間距分別為2,3和4 cm,初始相位差為45°時(shí),耦合后的聲壓幅值分別為8.071 5,9.974 0和10.654 5 kPa;在初始相位差為180°時(shí),耦合后的聲壓幅值分別為1 280,856和207 Pa??梢?,隨著2列熱聲振蕩器間距的增大,初始相位差的變化對(duì)合成聲壓幅值的影響越來越大。從圖5可以看出:在熱聲振蕩器間距不同時(shí),耦合后的相位差隨著初始相位差的增大而增大;當(dāng)2列熱聲振蕩器間距為4 cm時(shí),初始相位差對(duì)耦合后相位差的影響最大??梢姡谥C振腔尺寸結(jié)構(gòu)不變時(shí),改變熱聲振蕩器間距對(duì)協(xié)同作用的影響很小。

    熱聲振蕩器間距/cm:1—2;2—3;3—4。

    熱聲振蕩器間距/cm:1—2;2—3;3—4。

    2.3 諧振腔長度對(duì)聲耦合的影響

    當(dāng)諧振腔結(jié)構(gòu)尺寸固定時(shí),耦合強(qiáng)度也固定。當(dāng)耦合強(qiáng)度較弱時(shí),2個(gè)熱聲振蕩器的耦合作用程度較弱,對(duì)熱聲振蕩器的初始狀態(tài)干擾較弱;當(dāng)2個(gè)熱聲振蕩器的初始相位差不同時(shí),所產(chǎn)生的耦合效果不同。下面研究諧振腔的尺寸結(jié)構(gòu)對(duì)2個(gè)熱聲振蕩器耦合程度的影響。采用圖2所示模型,諧振腔半徑為4 cm,2列諧振腔的間距為2 cm,諧振腔的長度變化范圍為2~6 cm。2列熱聲振蕩器的初始相位差分別為45°,90°,145°和180°。

    隨著諧振腔長度變化,耦合后振蕩幅值的變化趨勢不同,如圖6所示。從圖6可見:當(dāng)初始相位差為45°時(shí),合成聲壓幅值隨諧振腔長度的增加而減?。划?dāng)諧振腔長度為2 cm時(shí),合成聲壓幅值為16.974 5 kPa;而當(dāng)諧振腔長度為6 cm時(shí),合成聲壓幅值僅為3.587 kPa;當(dāng)初始相位差為180°時(shí),合成聲壓幅值變化不大,諧振腔長度為2 cm時(shí),合成聲壓幅值為2.217 kPa;當(dāng)諧振腔長度為6 cm時(shí),合成聲壓幅值為2.333 5 kPa;當(dāng)諧振腔尺寸為6 cm時(shí),合成聲壓幅值隨初始相位差的變化變化平緩,達(dá)到了一種穩(wěn)定的狀態(tài)。這種穩(wěn)定狀態(tài)也可以從耦合后的相位差隨諧振腔長度的變化關(guān)系中得到進(jìn)一步驗(yàn)證,如圖7所示。

    諧振腔長度/cm:1—2;2—3;3—4;4—5;5—6。

    諧振腔長度/cm:1—2;2—3;3—4;4—5;5—6。

    從圖7可以看出:當(dāng)初始相位差為45°時(shí),耦合后的相位差隨著諧振腔長度的不同差別較大;當(dāng)諧振腔的長度為2 cm時(shí),耦合后的相位差為36°;當(dāng)諧振腔的長度為6 cm時(shí),耦合后的相位差為171°;當(dāng)初始相位差為180°時(shí),在不同諧振腔長度下,耦合相位差都接近180°。值得注意的是:當(dāng)諧振腔的長度為 6 cm時(shí),耦合后的相位差隨初始相位差的變化變化不大,穩(wěn)定在180°左右,可以認(rèn)為發(fā)生了“鎖相”現(xiàn)象。由此可知:當(dāng)諧振腔的長度為6 cm時(shí),耦合強(qiáng)度較大,耦合效果好,2列熱聲振蕩器達(dá)到反向協(xié)同狀態(tài)。

    通過上述分析可知:諧振腔的長度對(duì)2列熱聲振蕩器的耦合效果有較大影響;在諧振腔長度為6 cm,半徑為4 cm,2列諧振腔間距為2 cm的條件下,耦合強(qiáng)度最大,合成聲壓幅值和耦合相位差幾乎不隨初始相位差的變化而變化,發(fā)現(xiàn)“鎖相”現(xiàn)象,熱聲振蕩器達(dá)到反向協(xié)同狀態(tài)。這說明當(dāng)2個(gè)或2個(gè)以上的熱聲振蕩器耦合時(shí),諧振腔的尺寸結(jié)構(gòu)對(duì)熱聲振蕩器的耦合效果有至關(guān)重要的影響;當(dāng)耦合強(qiáng)度較弱時(shí),初始相位差對(duì)熱聲振蕩器的耦合有很大影響,由于熱聲振蕩器的初始相位差是隨機(jī)的,所以,宏觀上體現(xiàn)為熱聲振蕩器耦合現(xiàn)象不穩(wěn)定;當(dāng)耦合強(qiáng)度較強(qiáng)時(shí),初始相位差對(duì)熱聲振蕩器的耦合影響較小,耦合后的熱聲振蕩器維持在一種穩(wěn)定狀態(tài),體現(xiàn)為耦合后的2列熱聲振蕩器的相位差維持不變,即出現(xiàn)“鎖相”現(xiàn)象,達(dá)到協(xié)同狀態(tài)。熱聲振蕩器的協(xié)同分為正向協(xié)同和反向協(xié)同。正向協(xié)同耦合的熱聲振蕩器的相位差一般在0°左右,反向協(xié)同耦合的熱聲振蕩器的相位差一般在180°左右。正向協(xié)同使耦合后的聲壓幅值大于單獨(dú)的熱聲振蕩器的聲壓幅值,在理想條件下,耦合后的聲壓幅值為各個(gè)單獨(dú)熱聲振蕩器聲壓幅值之和;反向協(xié)同則會(huì)削弱熱聲振蕩器的聲壓幅值。

    2.4 諧振頻率對(duì)聲耦合的影響

    在諧振腔長度為6 cm時(shí),耦合強(qiáng)度最大,耦合后2個(gè)熱聲振蕩器的相位差不變,幅值穩(wěn)定,2個(gè)熱聲振蕩器達(dá)到了反向協(xié)同狀態(tài)。然而,之前的結(jié)果都是在2個(gè)熱聲振蕩器具有相同的諧振頻率下得出的,根據(jù)耦合理論可知,振蕩頻率對(duì)振蕩耦合也具有重要影響。由熱聲熱機(jī)的振蕩特性可知,諧振頻率主要由諧振管的長度決定,對(duì)于1波長和4波長的駐波型熱聲熱機(jī),諧振頻率隨著諧振管的增大而減小,反之,則增大。在諧振腔長度為6 cm、諧振管間距為2 cm時(shí),通過改變諧振管的長度研究諧振頻率對(duì)耦合作用的影響。熱聲振蕩器A的長度不變,熱聲振蕩器B的長度分別增加0.4,0.6和1.0 mm。經(jīng)過耦合作用后諧振腔端面中心處的振蕩特性如圖8~10所示。

    耦合后的2列熱聲振蕩器的相位差呈現(xiàn)周期性變化,相位差從0°變化到180°。從圖8~10可以看出:在諧振腔端面中心,耦合后的壓力幅值呈現(xiàn)周期性變化,時(shí)強(qiáng)時(shí)弱,出現(xiàn)波包,在聲學(xué)上稱為拍頻現(xiàn)象。這表明當(dāng)諧振管長度差從0.4 mm變化到1.0 mm時(shí),協(xié)同狀態(tài)消失,2個(gè)熱聲振蕩器的相互作用呈現(xiàn)周期性變化。當(dāng)諧振管長度差為0.4,0.6和1.0 mm時(shí),諧振腔端面中心的拍頻與改變諧振管長度產(chǎn)生的頻率差如圖11所示。從圖11可以看出:拍頻與頻率差近似呈線性關(guān)系,拍頻隨著諧振頻率差的增大而增大。由此可見:在諧振腔尺寸結(jié)構(gòu)一定即耦合強(qiáng)度一定時(shí),2個(gè)熱聲振蕩器的諧振頻率對(duì)耦合作用有重要影響。所以,在高頻微型熱聲熱機(jī)耦合設(shè)計(jì)中,應(yīng)盡量保持諧振管尺寸的一致性。

    圖8 諧振管長度差為0.4 mm時(shí)諧振腔端面中心壓力耦合特性

    圖9 諧振管長度差為0.6 mm時(shí)諧振腔端面中心壓力耦合特性

    圖10 諧振管長度差為1.0 mm時(shí)諧振腔端面中心壓力耦合特性

    圖11 拍頻與諧振頻率差的關(guān)系

    3 結(jié)論

    1) 耦合強(qiáng)度較弱時(shí)初始相位差對(duì)熱聲振蕩器耦合的影響較大,耦合結(jié)果隨著初始相位差的改變而改變。

    2) 在諧振腔尺寸結(jié)構(gòu)不變時(shí),熱聲振蕩器間距對(duì)協(xié)同作用的影響很小。

    3) 諧振腔長度對(duì)熱聲振蕩器耦合作用有重要影響,在改變諧振腔長度時(shí)出現(xiàn)鎖相現(xiàn)象。

    4) 諧振頻率差對(duì)耦合作用影響很大,隨著諧振頻率的改變協(xié)同狀態(tài)消失,2個(gè)熱聲振蕩器耦合之后出現(xiàn)拍頻現(xiàn)象。

    [1] BACKHAUSS S, SWIFT G W. A thermoacoustic-Stirling heat engine[J]. Nature, 1999, 399(27): 335?338.

    [2] BAKHAUSS S, SWIFT G W. A thermoacoustic-Stirling heat engine: detailed study[J]. The Journal of the Acoustical Society of America, 2000, 107(6): 3148?3166.

    [3] LI Linyu, WU Zhanghua, HU Jianying. A novel heat-driven thermoacustic natural gas liquefaction system[J]. Energy, 2016, 117: 523?529.

    [4] SKARIA M, RASHEED KKA, SHAFI K A. Simulation studies on the performance of thermoacoustic prime movers and refrigerator[J]. Computers & Fluids, 2015, 111: 127?136.

    [5] 劉益才, 張明研, 黃謙, 等. 熱聲熱機(jī)板疊式回?zé)崞鹘Y(jié)構(gòu)數(shù)值計(jì)算[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010, 41(3): 1186?1189. LIU Yicai, ZHANG Mingyan, HUANG Qian, et al. The numerical calculation of thermoacoustic engine in stack regenerator structure[J]. Journal of Central South University (Science and Technology), 2010, 41(3): 1186?1189.

    [6] HASEGAWA S, YAMAGUCHI T, OSHINOYA Y. A thermoacoustic refrigerator driven by a low temperature differential, high-efficiency multistage thermoacoustic engine[J]. Applied Thermal Engineering, 2013, 58(1/2): 394?399.

    [7] ORINO Y, SAKAMOTO S I, INUI Y, et al. Numerical analysis of the effect of local diameter reduction on the critical temperature of thermoacoustic oscillations in a looped tube[J]. Japanese Journal of Applied Physics, 2011, 264/265(53): 1982?1987.

    [8] UEDA Y, KATO C. Stability analysis of thermally induced spontaneous gas oscillations in straight and looped tubes[J]. The Journal of the Acoustical Society of America, 2008, 124(2): 851?858.

    [9] GARRETT S L. Thermoacoustic engines and refrigerators[J]. American Institute of Physics, 2012, 1440: 9?22.

    [10] GARRETT S L. Thermoacoustic engines and refrigerators[J]. American Journal of Physics, 2004, 72(1): 11?17.

    [11] RODRIGUEZ I A, SYMKO O G. Miniature travelling wave thermoacoustic engine[J]. Journal of Acoustical Society of America, 2009, 125(4): 2562.

    [12] 夏淯博. 高頻熱聲熱機(jī)及其耦合特性研究[D]. 長沙: 中南大學(xué)能源科學(xué)與工程學(xué)院, 2015: 23?32. XIA Yubo. Research on acoustic coupling characteristics of high frequency thermoacoustic engines[D]. Changsha: Central South University. School of Energy Science and Engineering, 2015: 23?32.

    [13] SPOOR P S, SWIFT G W. Mode-locking of acoustic resonators and its applications to vibration cancelation in acoustic heat engines[J]. Journal of Acoustical Society of America, 1999, 106(3): 1353?1362.

    [14] SPOOR P S, SWIFT G W. The huygens entrainment phenomenon and thermocaocustic engines[J]. Journal of Acoustical Society of America, 2000, 108(2): 588?600.

    [15] SURATHU R. A study of coupling two thermoacoustic laser[D]. Utah, USA: University of Utah. Department of Physics, 2013: 80?113.

    [16] 劉曉蕾. 高頻熱聲振蕩器的協(xié)同行為研究[M]. 武漢: 華中科技大學(xué)能源與動(dòng)力工程學(xué)院, 2013: 78?123. LIU Xiaolei. The study on synergistic behaviors of high frequency thermoacoustic oscillators[M]. Wuhan: Huazhong University of Science and Technology. School of Energy and Power Engineering, 2013: 78?123.

    [17] GILLMAN B J. A Study of coupled thermoacoustic engines[D]. Utah, USA: University of Utah. Department of Physics, 2009: 56?123.

    (編輯 陳燦華)

    Acoustic coupling characteristics of high frequency thermoacoustic engines

    LIU Yicai, MO Shuanglin, LEI Binyi, XIA Yubo, MA Ying

    (School of Energy Science and Engineering, Central South University, Changsha 410083, China)

    The synergy theory was introduced into the acoustic coupling research of high-frequency thermoacoustic engines. Based on the synergy theory, the important parameters that had impact on the coupling of high-frequency thermoacoustic engines were obtained. Then the numerical model for two thermoacoustic engines in the same cavity was established. The influence of initial phase difference between the thermoacoustic engines, the distance between the thermoacoustic engines, the length of the cavity and the length of the resonance tube on the acoustic coupling characteristic were analyzed. The results show that the initial phase difference between the thermoacoustic engines, the resonant frequency and the coupling strength have important impact on the coupling of high-frequency thermoacoustic engines. The change of the distance between the thermoacoustic engines has small effect on synergistic. When the length of the cavity changes without changing the other structure parameters, and the length of cavity is 6 cm, the synergistic phenomenon appears and the coupling phase difference and coupling sound pressure amplitude reach a steady state. Changing the length of the resonance tube when the cavity length is 6 cm, the resonant frequency difference between the thermoacoustic oscillators has great influence on the coupling, and the synergistic phenomenon disappears and the phenomenon of beat frequency appears. The beat frequency and frequency difference between the thermoacoustic oscillators is approximately linear.

    high-frequency thermoacoustic engine; acoustic coupling; numerical simulation

    10.11817/j.issn.1672-7207.2017.09.035

    TB65

    A

    1672?7207(2017)09?2529?07

    2016?10?15;

    2016?12?21

    國家自然科學(xué)基金資助項(xiàng)目(51276201) (Project(51276201) supported by the National Natural Science Foundation of China)

    劉益才,博士,教授,從事熱聲熱機(jī)、斯特林制冷機(jī)、高效蓄冷蓄熱材料、微型低溫制冷機(jī)以及微型制冷系統(tǒng)振動(dòng)和噪聲抑制等研究;E-mail: lyccsu@csu.edu.cn

    猜你喜歡
    熱聲諧振腔聲壓
    基于嘴唇處的聲壓數(shù)據(jù)確定人體聲道半徑
    熱驅(qū)動(dòng)熱聲制冷技術(shù)發(fā)展現(xiàn)狀與展望
    非對(duì)稱雙級(jí)環(huán)路行波熱聲熱機(jī)的實(shí)驗(yàn)研究*
    車輛結(jié)構(gòu)噪聲傳遞特性及其峰值噪聲成因的分析
    汽車工程(2018年12期)2019-01-29 06:46:36
    熱聲效應(yīng)及其應(yīng)用研究進(jìn)展
    用于小型銣如原子鐘中介質(zhì)諧振腔激勵(lì)分析
    電子測試(2018年11期)2018-06-26 05:56:12
    蘊(yùn)含在熱聲制冷技術(shù)中的物理知識(shí)
    基于GIS內(nèi)部放電聲壓特性進(jìn)行閃絡(luò)定位的研究
    電測與儀表(2016年9期)2016-04-12 00:30:02
    波導(dǎo)諧振腔Fano共振特性研究
    微波諧振腔模式數(shù)的程序求解法
    国产成人欧美| 午夜免费观看网址| 中文字幕久久专区| 久久香蕉激情| 中文字幕另类日韩欧美亚洲嫩草| 亚洲伊人色综图| 欧美性长视频在线观看| 熟女少妇亚洲综合色aaa.| 国产精品 国内视频| 长腿黑丝高跟| 久久天堂一区二区三区四区| 日韩中文字幕欧美一区二区| 窝窝影院91人妻| 婷婷丁香在线五月| 一本久久中文字幕| 黄色视频,在线免费观看| 99国产极品粉嫩在线观看| 曰老女人黄片| 在线观看舔阴道视频| 久久影院123| 免费在线观看黄色视频的| 日本 欧美在线| 久久久精品欧美日韩精品| 免费av毛片视频| 日韩 欧美 亚洲 中文字幕| 精品国产超薄肉色丝袜足j| 欧美最黄视频在线播放免费| 在线观看日韩欧美| a在线观看视频网站| 亚洲狠狠婷婷综合久久图片| 国产蜜桃级精品一区二区三区| 亚洲一区二区三区色噜噜| 色播在线永久视频| 最新美女视频免费是黄的| 久9热在线精品视频| 免费高清在线观看日韩| 在线视频色国产色| 午夜视频精品福利| 一二三四在线观看免费中文在| 好看av亚洲va欧美ⅴa在| 国产激情久久老熟女| 欧美人与性动交α欧美精品济南到| 99国产精品免费福利视频| 国产精品国产高清国产av| 国产aⅴ精品一区二区三区波| 午夜成年电影在线免费观看| 操美女的视频在线观看| 人成视频在线观看免费观看| 欧美黄色片欧美黄色片| av在线播放免费不卡| 韩国精品一区二区三区| 亚洲成人免费电影在线观看| 可以在线观看毛片的网站| 老司机午夜福利在线观看视频| av超薄肉色丝袜交足视频| 51午夜福利影视在线观看| 欧美精品亚洲一区二区| 亚洲中文日韩欧美视频| 久久伊人香网站| 色av中文字幕| 免费久久久久久久精品成人欧美视频| 国产单亲对白刺激| 可以免费在线观看a视频的电影网站| 日韩高清综合在线| 日本 欧美在线| 欧美成狂野欧美在线观看| 啪啪无遮挡十八禁网站| 麻豆av在线久日| 视频在线观看一区二区三区| 亚洲av第一区精品v没综合| 一夜夜www| 最近最新中文字幕大全免费视频| 丝袜美腿诱惑在线| 久久中文字幕人妻熟女| 黄色视频,在线免费观看| 欧美黄色片欧美黄色片| 久热这里只有精品99| 99国产极品粉嫩在线观看| 老熟妇仑乱视频hdxx| 伊人久久大香线蕉亚洲五| 中文亚洲av片在线观看爽| 变态另类成人亚洲欧美熟女 | av电影中文网址| 一级毛片女人18水好多| 可以在线观看毛片的网站| 午夜日韩欧美国产| 韩国av一区二区三区四区| 亚洲专区字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 曰老女人黄片| www日本在线高清视频| 免费一级毛片在线播放高清视频 | 在线永久观看黄色视频| 久久国产精品男人的天堂亚洲| 女人被躁到高潮嗷嗷叫费观| 精品少妇一区二区三区视频日本电影| 激情在线观看视频在线高清| 可以在线观看毛片的网站| 99久久99久久久精品蜜桃| 精品久久久久久久人妻蜜臀av | 中文亚洲av片在线观看爽| 亚洲人成网站在线播放欧美日韩| 日韩精品中文字幕看吧| 无限看片的www在线观看| 欧美日韩瑟瑟在线播放| 黄色成人免费大全| 欧美精品啪啪一区二区三区| 亚洲av片天天在线观看| 欧美乱码精品一区二区三区| 精品高清国产在线一区| 精品国产乱子伦一区二区三区| 免费久久久久久久精品成人欧美视频| 亚洲天堂国产精品一区在线| 国产av一区在线观看免费| 亚洲avbb在线观看| 亚洲男人的天堂狠狠| 熟女少妇亚洲综合色aaa.| 国产精品 欧美亚洲| 欧美日韩亚洲国产一区二区在线观看| 欧美中文综合在线视频| 老司机午夜福利在线观看视频| 亚洲成国产人片在线观看| 欧美成人免费av一区二区三区| aaaaa片日本免费| 欧美黄色淫秽网站| 欧美中文日本在线观看视频| 国产亚洲欧美98| 一卡2卡三卡四卡精品乱码亚洲| 精品日产1卡2卡| 亚洲国产中文字幕在线视频| 国产高清videossex| 精品久久久久久成人av| 午夜久久久在线观看| 国产精品久久久av美女十八| 少妇熟女aⅴ在线视频| 少妇熟女aⅴ在线视频| 手机成人av网站| 999久久久精品免费观看国产| 欧美日韩黄片免| 亚洲精品一区av在线观看| 国产真人三级小视频在线观看| 久久亚洲精品不卡| 亚洲人成77777在线视频| 色播亚洲综合网| 三级毛片av免费| 国产精品久久电影中文字幕| 高清在线国产一区| 男男h啪啪无遮挡| 757午夜福利合集在线观看| 亚洲欧美激情在线| av欧美777| 黄色视频,在线免费观看| www.www免费av| 波多野结衣一区麻豆| 亚洲第一欧美日韩一区二区三区| 日韩欧美国产在线观看| 人妻丰满熟妇av一区二区三区| 精品久久久精品久久久| 国产精品99久久99久久久不卡| 天天躁夜夜躁狠狠躁躁| 天天躁狠狠躁夜夜躁狠狠躁| 久久亚洲精品不卡| 国产精品日韩av在线免费观看 | 午夜日韩欧美国产| 精品国产亚洲在线| 一区福利在线观看| 亚洲人成77777在线视频| 丝袜人妻中文字幕| 久久久国产成人精品二区| 久久人妻av系列| 一本大道久久a久久精品| 精品国产超薄肉色丝袜足j| 日本五十路高清| 一进一出好大好爽视频| 19禁男女啪啪无遮挡网站| av网站免费在线观看视频| 日韩高清综合在线| 可以免费在线观看a视频的电影网站| 亚洲精品国产色婷婷电影| 看片在线看免费视频| 岛国在线观看网站| 纯流量卡能插随身wifi吗| 午夜福利欧美成人| 亚洲无线在线观看| 日韩大码丰满熟妇| www国产在线视频色| 亚洲欧洲精品一区二区精品久久久| 中文字幕高清在线视频| 欧美激情久久久久久爽电影 | 久久亚洲真实| 久久久久久亚洲精品国产蜜桃av| 国内久久婷婷六月综合欲色啪| 国产激情欧美一区二区| 亚洲最大成人中文| 亚洲专区中文字幕在线| 婷婷六月久久综合丁香| 老汉色av国产亚洲站长工具| 亚洲国产精品成人综合色| 欧美色视频一区免费| 91麻豆av在线| 国产激情欧美一区二区| 欧美 亚洲 国产 日韩一| av电影中文网址| 亚洲国产欧美日韩在线播放| 在线视频色国产色| 91精品三级在线观看| 免费搜索国产男女视频| www.www免费av| 欧美另类亚洲清纯唯美| 亚洲情色 制服丝袜| 久久久久久久久中文| 成年版毛片免费区| 国产精品亚洲一级av第二区| 久久久久久人人人人人| 无遮挡黄片免费观看| 18禁国产床啪视频网站| 中文字幕人妻丝袜一区二区| av超薄肉色丝袜交足视频| 大码成人一级视频| 国产成人免费无遮挡视频| 曰老女人黄片| 久久人人精品亚洲av| 久久亚洲精品不卡| 欧美成人性av电影在线观看| 桃红色精品国产亚洲av| 色综合欧美亚洲国产小说| 高潮久久久久久久久久久不卡| 精品久久久精品久久久| 999精品在线视频| 久久婷婷人人爽人人干人人爱 | 欧美色欧美亚洲另类二区 | 黄色女人牲交| 午夜亚洲福利在线播放| 18禁裸乳无遮挡免费网站照片 | 淫秽高清视频在线观看| 中亚洲国语对白在线视频| 国内精品久久久久久久电影| 亚洲欧洲精品一区二区精品久久久| 免费一级毛片在线播放高清视频 | 成人亚洲精品一区在线观看| 久久中文字幕一级| 91国产中文字幕| 男女之事视频高清在线观看| 久久精品国产亚洲av高清一级| 久久性视频一级片| 中文字幕av电影在线播放| 亚洲精品在线美女| 久久久久久大精品| 国产麻豆成人av免费视频| 制服人妻中文乱码| 欧美乱色亚洲激情| 女人被躁到高潮嗷嗷叫费观| 日本五十路高清| 午夜福利免费观看在线| 十八禁人妻一区二区| 校园春色视频在线观看| 看片在线看免费视频| 国语自产精品视频在线第100页| 成人av一区二区三区在线看| 国产一区在线观看成人免费| 露出奶头的视频| 欧美人与性动交α欧美精品济南到| 嫩草影院精品99| 亚洲 国产 在线| 久久伊人香网站| 久久久久久久精品吃奶| 一边摸一边抽搐一进一出视频| x7x7x7水蜜桃| 中文字幕久久专区| 欧美精品亚洲一区二区| 国产精品久久久久久亚洲av鲁大| 欧美黄色片欧美黄色片| 老鸭窝网址在线观看| 精品久久蜜臀av无| 色老头精品视频在线观看| 亚洲精品一区av在线观看| av有码第一页| 美女 人体艺术 gogo| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产高清在线一区二区三 | 亚洲av片天天在线观看| 亚洲三区欧美一区| 午夜福利,免费看| 真人一进一出gif抽搐免费| 国产日韩一区二区三区精品不卡| 国产精品一区二区精品视频观看| 日韩成人在线观看一区二区三区| 亚洲人成77777在线视频| 国产成人av教育| 色综合婷婷激情| 男女午夜视频在线观看| 久久久久久大精品| 十分钟在线观看高清视频www| 午夜日韩欧美国产| 人人妻人人澡欧美一区二区 | 久久影院123| 999久久久精品免费观看国产| 免费在线观看完整版高清| 久久久久国内视频| 最近最新中文字幕大全电影3 | 亚洲男人的天堂狠狠| 国产成人精品久久二区二区免费| 欧美亚洲日本最大视频资源| 黄片小视频在线播放| 日本 欧美在线| 亚洲天堂国产精品一区在线| 成年版毛片免费区| 在线永久观看黄色视频| 久久久久久久久久久久大奶| 午夜久久久在线观看| 777久久人妻少妇嫩草av网站| 一a级毛片在线观看| а√天堂www在线а√下载| 成年女人毛片免费观看观看9| 亚洲一区二区三区不卡视频| 日韩欧美国产一区二区入口| 欧美激情极品国产一区二区三区| 给我免费播放毛片高清在线观看| 国产精品永久免费网站| 97人妻天天添夜夜摸| 日韩 欧美 亚洲 中文字幕| 制服丝袜大香蕉在线| 很黄的视频免费| 久久精品亚洲精品国产色婷小说| 大型黄色视频在线免费观看| 色综合站精品国产| 非洲黑人性xxxx精品又粗又长| 啦啦啦韩国在线观看视频| 欧美亚洲日本最大视频资源| 他把我摸到了高潮在线观看| 麻豆av在线久日| 欧美大码av| 12—13女人毛片做爰片一| 女生性感内裤真人,穿戴方法视频| 国产亚洲av嫩草精品影院| www.精华液| 国产精品国产高清国产av| 午夜a级毛片| 久久天躁狠狠躁夜夜2o2o| 国产人伦9x9x在线观看| 亚洲色图av天堂| 欧美黑人欧美精品刺激| 久久香蕉激情| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 国产亚洲精品久久久久5区| 免费久久久久久久精品成人欧美视频| 国产精品1区2区在线观看.| 日本 av在线| 十八禁人妻一区二区| 最近最新中文字幕大全电影3 | 在线播放国产精品三级| 一级黄色大片毛片| 久久中文看片网| 久99久视频精品免费| 久久 成人 亚洲| av天堂在线播放| 国产精品久久视频播放| 欧美大码av| 91字幕亚洲| 亚洲欧美一区二区三区黑人| 亚洲一区高清亚洲精品| 国产精品野战在线观看| 欧美精品亚洲一区二区| 免费在线观看完整版高清| 18禁黄网站禁片午夜丰满| 91在线观看av| 少妇被粗大的猛进出69影院| 老鸭窝网址在线观看| 法律面前人人平等表现在哪些方面| av免费在线观看网站| 一二三四社区在线视频社区8| 热99re8久久精品国产| 又黄又粗又硬又大视频| 国产真人三级小视频在线观看| 电影成人av| 久久国产亚洲av麻豆专区| 久久久久亚洲av毛片大全| 90打野战视频偷拍视频| 欧美色欧美亚洲另类二区 | 国产精品爽爽va在线观看网站 | 成在线人永久免费视频| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 久久亚洲精品不卡| 免费在线观看黄色视频的| 9191精品国产免费久久| 在线视频色国产色| 国产成人影院久久av| 中文字幕久久专区| 村上凉子中文字幕在线| 久久亚洲真实| 日日爽夜夜爽网站| 久久精品aⅴ一区二区三区四区| 久久午夜亚洲精品久久| 麻豆一二三区av精品| 午夜影院日韩av| 成人手机av| 欧美 亚洲 国产 日韩一| 日韩欧美一区视频在线观看| 亚洲五月婷婷丁香| 精品久久久精品久久久| 99re在线观看精品视频| 久久精品成人免费网站| 91老司机精品| 欧美在线黄色| 欧美日韩一级在线毛片| 91九色精品人成在线观看| 免费av毛片视频| 亚洲第一欧美日韩一区二区三区| 国产精品影院久久| 国产一区二区在线av高清观看| 99香蕉大伊视频| 少妇熟女aⅴ在线视频| 久久精品国产清高在天天线| 亚洲色图 男人天堂 中文字幕| 色播在线永久视频| 激情在线观看视频在线高清| 久久中文字幕一级| 久9热在线精品视频| 国产单亲对白刺激| 欧美乱色亚洲激情| 男女午夜视频在线观看| 国产三级黄色录像| 久久这里只有精品19| 国产97色在线日韩免费| 大型av网站在线播放| 国内精品久久久久久久电影| 亚洲精品av麻豆狂野| ponron亚洲| 亚洲黑人精品在线| 一级毛片精品| 国产三级在线视频| av超薄肉色丝袜交足视频| 老司机福利观看| 亚洲熟女毛片儿| av有码第一页| 久久草成人影院| 国产午夜福利久久久久久| 男女做爰动态图高潮gif福利片 | 欧美黑人欧美精品刺激| 久9热在线精品视频| 国产亚洲精品av在线| 国产男靠女视频免费网站| 欧美激情高清一区二区三区| 精品高清国产在线一区| 老司机在亚洲福利影院| 精品人妻1区二区| 久久中文字幕一级| 国产激情久久老熟女| 国产精品一区二区精品视频观看| 亚洲成人免费电影在线观看| 日韩欧美一区视频在线观看| 两个人视频免费观看高清| 三级毛片av免费| 成人永久免费在线观看视频| 国产精品久久久久久亚洲av鲁大| 欧美激情极品国产一区二区三区| 18禁裸乳无遮挡免费网站照片 | 美女高潮喷水抽搐中文字幕| 桃红色精品国产亚洲av| 最近最新中文字幕大全电影3 | bbb黄色大片| 久久久久久久久中文| 国产一区二区三区视频了| 亚洲av日韩精品久久久久久密| 久99久视频精品免费| 欧美另类亚洲清纯唯美| 国产成+人综合+亚洲专区| 美女大奶头视频| 免费看十八禁软件| 国内精品久久久久精免费| 一区二区三区高清视频在线| 国产精品,欧美在线| 国产成人精品久久二区二区91| 性欧美人与动物交配| 成人国产综合亚洲| 美女高潮喷水抽搐中文字幕| 两性夫妻黄色片| 久久国产亚洲av麻豆专区| 伊人久久大香线蕉亚洲五| videosex国产| 亚洲精品国产一区二区精华液| 久久国产精品影院| 成人亚洲精品一区在线观看| 在线av久久热| 少妇 在线观看| 天堂√8在线中文| 欧美黄色淫秽网站| 老熟妇乱子伦视频在线观看| 日韩欧美一区二区三区在线观看| 大香蕉久久成人网| 少妇熟女aⅴ在线视频| 一区二区日韩欧美中文字幕| 18美女黄网站色大片免费观看| 国产在线精品亚洲第一网站| 精品国产超薄肉色丝袜足j| АⅤ资源中文在线天堂| 国产精品秋霞免费鲁丝片| 日日夜夜操网爽| 男女下面进入的视频免费午夜 | 亚洲专区中文字幕在线| 国产av精品麻豆| 级片在线观看| 天天躁夜夜躁狠狠躁躁| 久久久久久久午夜电影| www国产在线视频色| 97人妻天天添夜夜摸| 午夜福利一区二区在线看| 男人操女人黄网站| 日韩国内少妇激情av| 1024香蕉在线观看| 午夜老司机福利片| 曰老女人黄片| 久久精品影院6| 狂野欧美激情性xxxx| 国产精品综合久久久久久久免费 | 国产精品永久免费网站| 国产男靠女视频免费网站| 久久香蕉激情| 国产精品综合久久久久久久免费 | 久久久久久久久免费视频了| 午夜a级毛片| 校园春色视频在线观看| 亚洲专区中文字幕在线| 亚洲人成电影观看| 又紧又爽又黄一区二区| 国产精品二区激情视频| 一卡2卡三卡四卡精品乱码亚洲| 日本免费一区二区三区高清不卡 | 在线观看免费视频日本深夜| 91av网站免费观看| 咕卡用的链子| 免费高清在线观看日韩| 国产在线精品亚洲第一网站| 国产精品久久久人人做人人爽| 老司机午夜福利在线观看视频| 国产男靠女视频免费网站| 国产国语露脸激情在线看| 精品一区二区三区视频在线观看免费| 麻豆av在线久日| 精品不卡国产一区二区三区| 久久久久久久精品吃奶| 99国产精品一区二区三区| 日韩欧美免费精品| 成人永久免费在线观看视频| 大陆偷拍与自拍| 国产1区2区3区精品| 亚洲成av片中文字幕在线观看| 国产精品 欧美亚洲| 在线视频色国产色| 色av中文字幕| 一区福利在线观看| 国产精品永久免费网站| 国产一区在线观看成人免费| 可以免费在线观看a视频的电影网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品久久国产高清桃花| 亚洲午夜精品一区,二区,三区| 欧美日韩亚洲综合一区二区三区_| 黄片播放在线免费| 国产成年人精品一区二区| 精品久久久久久久久久免费视频| 国产91精品成人一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产精品sss在线观看| 桃色一区二区三区在线观看| 黄色 视频免费看| 69精品国产乱码久久久| 激情视频va一区二区三区| 欧美在线一区亚洲| 麻豆成人av在线观看| av视频在线观看入口| 不卡av一区二区三区| 免费女性裸体啪啪无遮挡网站| 在线观看舔阴道视频| videosex国产| 搞女人的毛片| 亚洲精华国产精华精| 国产精品98久久久久久宅男小说| 两个人视频免费观看高清| 欧美性长视频在线观看| 国产精品av久久久久免费| 一级毛片精品| 国内毛片毛片毛片毛片毛片| 淫妇啪啪啪对白视频| 国产精品久久久人人做人人爽| 91麻豆av在线| 在线观看免费视频网站a站| 亚洲精品国产一区二区精华液| 精品一品国产午夜福利视频| 亚洲avbb在线观看| 亚洲成人国产一区在线观看| 亚洲av电影不卡..在线观看| 亚洲国产看品久久| svipshipincom国产片| 国产主播在线观看一区二区| 欧美成人免费av一区二区三区| 欧美在线黄色| 黄色丝袜av网址大全| 午夜免费观看网址| 人妻丰满熟妇av一区二区三区| 久久久久久久精品吃奶| 中文字幕精品免费在线观看视频| 天天添夜夜摸| 国产伦一二天堂av在线观看| 国产精品av久久久久免费| 色播在线永久视频| 两个人视频免费观看高清| 手机成人av网站| 国产成人影院久久av| 亚洲成人免费电影在线观看| www.www免费av| 亚洲成人久久性| 日本撒尿小便嘘嘘汇集6|