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

    高通量連續(xù)流微反應(yīng)器的流動(dòng)與傳熱數(shù)值模擬

    2021-12-16 01:36:00王麗軍成有為
    關(guān)鍵詞:連續(xù)流徑向反應(yīng)器

    朱 杰,潘 騰,王麗軍,成有為,李 希

    浙江大學(xué)化學(xué)工程與生物工程學(xué)院,浙江 杭州 310027

    隨著連續(xù)化、自動(dòng)化和智能化生產(chǎn)觀念的深入,連續(xù)流反應(yīng)技術(shù)在精細(xì)化工和藥物化工領(lǐng)域快速發(fā)展[1-2]。作為一種典型的連續(xù)流反應(yīng)器,微反應(yīng)器憑借自身高效的傳熱和傳質(zhì)性能在有機(jī)合成反應(yīng)中得到廣泛應(yīng)用[3-4]。然而,微尺度在帶來(lái)優(yōu)越的傳遞性能的同時(shí)也限制了流動(dòng)通量,微反應(yīng)器的大規(guī)模工業(yè)應(yīng)用仍有許多挑戰(zhàn)[5-6]。

    微反應(yīng)器的放大策略主要有三種[7]:一是平行數(shù)目放大,該方法的優(yōu)點(diǎn)是可以保持傳遞性能不變,難點(diǎn)在于如何實(shí)現(xiàn)成千上萬(wàn)臺(tái)設(shè)備的流量均布;二是串聯(lián)數(shù)目放大,該方法無(wú)須流量分配器,可以應(yīng)用于較長(zhǎng)反應(yīng)停留時(shí)間的反應(yīng)以及多步反應(yīng),但是能夠?qū)崿F(xiàn)的通量有限,鮮有超過(guò)100 mL/min處理量的報(bào)道;三是幾何尺寸放大,該方法可以實(shí)現(xiàn)通量的數(shù)量級(jí)增長(zhǎng),降低通道堵塞風(fēng)險(xiǎn),難點(diǎn)在于保留微尺度的傳遞性能??祵幑净趲缀纬叽绶糯蟮牟呗栽O(shè)計(jì)了一種先進(jìn)流反應(yīng)器(Advanced Flow Reactor,AFR)[8-10],該反應(yīng)器的特征幾何尺寸為寬1 mm、高1.1 mm,流動(dòng)通量達(dá)到10~200 mL/min,氣-液相界面積為160~1 300 m-1[11],液-液相界面積為1 000~10 000 m-1[12],僅需七個(gè)AFR組合就可實(shí)現(xiàn)每分鐘公斤級(jí)的原料處理量,展現(xiàn)出了一定的工業(yè)應(yīng)用潛力。Fang等[13]自行設(shè)計(jì)了一種可以連續(xù)合成季銨鹽的反應(yīng)器,該反應(yīng)器由內(nèi)徑為9 mm、總長(zhǎng)為200 m的鈦材蛇形圓管構(gòu)成,管內(nèi)容積為15.7 L,在管內(nèi)插入組合式靜態(tài)混合器以增強(qiáng)傳熱和傳質(zhì)效果,結(jié)果表明該反應(yīng)器可以在幾分鐘內(nèi)實(shí)現(xiàn)相當(dāng)于間歇反應(yīng)器10 h的生產(chǎn)能力,年生產(chǎn)能力超過(guò)1 000 t。

    本研究基于擴(kuò)大通道尺寸并內(nèi)置三維螺旋構(gòu)件的思路設(shè)計(jì)了一種高通量連續(xù)流反應(yīng)器(High-throughput Continuous-flow Reactor,HTCR),該反應(yīng)器在保持較優(yōu)的混合與傳熱性能的前提下,通量比常見(jiàn)微反應(yīng)器增加1~3個(gè)數(shù)量級(jí),為微反應(yīng)器的幾何放大提供了一條新的思路。本工作首先介紹HTCR的幾何構(gòu)型,再采用計(jì)算流體力學(xué)(CFD)的方法對(duì)反應(yīng)器基本單元的對(duì)流傳熱過(guò)程進(jìn)行模擬,分析其流動(dòng)與傳熱規(guī)律,隨后采用單因素分析的方法考察雷諾數(shù)、普朗特?cái)?shù)、螺旋片長(zhǎng)徑比和扭轉(zhuǎn)角對(duì)傳熱的影響,最后通過(guò)數(shù)值擬合得到努塞爾數(shù)與流動(dòng)摩擦系數(shù)的計(jì)算關(guān)聯(lián)式。

    1 反應(yīng)器幾何構(gòu)型設(shè)計(jì)

    流體的旋轉(zhuǎn)、拉伸及剪切,對(duì)強(qiáng)化流體傳熱、混合和傳質(zhì)有重要作用。借鑒SK靜態(tài)混合器的結(jié)構(gòu),提出了一種新的連續(xù)流反應(yīng)器型式,其基本單元為內(nèi)置SK螺旋片的六角直管。如圖1所示,六角管橫截面為正六邊形,SK螺旋片為左旋、右旋交替出現(xiàn)的周期性結(jié)構(gòu),相鄰元件頭尾交叉角90°,將六角管與SK螺旋片組合在一起,即構(gòu)成HTCR的基本單元。

    圖1 內(nèi)置SK螺旋片的六角管示意Fig.1 Schematic diagram of hexagonal tube embedded SK helical plates

    將基本單元沿六個(gè)外壁面法向陣列再切割形成反應(yīng)管束,如圖2所示。將反應(yīng)管束軸向的兩個(gè)端面套上固定板,在板上按一定的規(guī)則開(kāi)槽使每根六角管依次連通,連通后的槽口布置如圖3所示。將槽口表面封閉后再將z方向兩個(gè)端面用矩形板封閉,并在第一根六角管上開(kāi)若干進(jìn)料口,在最后一根六角管上開(kāi)一個(gè)出料口,此時(shí)反應(yīng)器的流體通道構(gòu)建完畢,反應(yīng)器內(nèi)物料流向如圖4所示。

    圖2 反應(yīng)管束示意Fig.2 Schematic diagram of reaction tubes bunch

    反應(yīng)管束中六角管外的區(qū)域?yàn)閾Q熱介質(zhì)區(qū)域,在x方向上兩個(gè)端面分別設(shè)置換熱介質(zhì)進(jìn)、出口,在進(jìn)口與反應(yīng)管束間設(shè)置分布器,各部件的裝配關(guān)系如圖5所示。

    圖5 零部件裝配Fig.5 Assembly diagram of parts

    反應(yīng)器運(yùn)行時(shí),反應(yīng)物分別通過(guò)不同的進(jìn)料口進(jìn)入第一根六角管,在螺旋片的作用下進(jìn)行快速混合與高效傳熱,然后沿立體蛇形路徑流至最后一根六角管再?gòu)某隽峡诹鞒觥Q熱介質(zhì)通過(guò)管線進(jìn)入分布器,從分布器孔流入反應(yīng)管束殼程,通過(guò)間壁傳熱與管內(nèi)反應(yīng)物進(jìn)行熱交換,兩流體錯(cuò)流換熱。

    2 數(shù)學(xué)模型

    2.1 幾何模型

    抽取流體通道的最小重復(fù)單元作為幾何模型,如圖6所示,l為一個(gè)螺片的長(zhǎng)度,橫截面正六邊形邊長(zhǎng)為5 mm,內(nèi)置n組螺旋片。

    圖6 幾何模型Fig.6 Geometric model

    2.2 控制方程

    對(duì)于不可壓縮牛頓流體的三維穩(wěn)態(tài)流動(dòng),忽略質(zhì)量力,以張量形式表示的控制方程如下:

    連續(xù)性方程

    動(dòng)量方程

    能量方程

    由于模擬的雷諾數(shù)(Re)在湍流區(qū)(對(duì)于SK螺旋結(jié)構(gòu),Re>150即為湍流區(qū)[14])。根據(jù)以往對(duì)于三維剪切流的研究,剪切應(yīng)力傳輸模型(SST)k-ω湍流模型可以更準(zhǔn)確地描述流體行為[15],故本研究選用該模型進(jìn)行計(jì)算。

    2.3 網(wǎng)格劃分及無(wú)關(guān)性分析

    采用軟件Star ccm+的網(wǎng)格劃分功能,貼近壁面設(shè)置5層棱柱層網(wǎng)格以解析邊界層,保證算例中壁面y+值為0.1~1.2,網(wǎng)格劃分情況如圖7所示。

    圖7 網(wǎng)格劃分進(jìn)口段(a)和SK段(b)Fig.7 Mesh of inlet section(a) and SK structure section(b)

    網(wǎng)格無(wú)關(guān)性分析的邊界條件設(shè)置:進(jìn)口速度為0.3 m/s,進(jìn)口溫度為60 ℃;出口為壓力邊界,其值為0 Pa;幾何壁面為無(wú)滑移邊界條件,外壁面溫度設(shè)為等溫80 ℃,螺旋片壁面設(shè)置為絕熱,流動(dòng)介質(zhì)設(shè)置為水。由疏到密選取了163 659,262 874,418102,566 706,866 350和14 490 616六個(gè)網(wǎng)格總數(shù)目進(jìn)行計(jì)算,以模擬對(duì)象的努塞爾數(shù)(Nu)和摩擦系數(shù)(f)作為判斷依據(jù),結(jié)果如圖8所示。結(jié)果顯示,f在整個(gè)考察范圍內(nèi)均保持穩(wěn)定,當(dāng)網(wǎng)格數(shù)大于56.6萬(wàn)時(shí),Nu幾乎不隨網(wǎng)格數(shù)的增長(zhǎng)而變化,此時(shí)網(wǎng)格尺寸已經(jīng)可以消除其本身對(duì)計(jì)算結(jié)果的影響,本文案例均采用該尺寸進(jìn)行計(jì)算。

    圖8 網(wǎng)格無(wú)關(guān)性分析Fig.8 Grid independence analysis

    為了排除進(jìn)出口對(duì)計(jì)算結(jié)果的干擾,進(jìn)一步考察螺旋片組數(shù)(n)對(duì)流動(dòng)和傳熱的影響,模擬結(jié)果如圖9所示。當(dāng)n≥6時(shí),Nu和f趨于平緩,說(shuō)明此時(shí)螺旋片組數(shù)對(duì)于傳熱和流動(dòng)阻力已經(jīng)幾乎沒(méi)有影響,故本模擬計(jì)算的n均以6為基準(zhǔn)。

    圖9 螺旋片組數(shù)對(duì)Nu和f的影響Fig.9 Influence of group number of spiral pieces on Nu and f

    3 結(jié)果與討論

    3.1 速度場(chǎng)與二次流

    圖10展示了模擬對(duì)象的流線圖??梢钥闯隽黧w在管內(nèi)呈現(xiàn)出不斷分裂合并的流動(dòng)狀況,與此同時(shí),流體在徑向上還有大量的穿插和交互,這些徑向上的二次流對(duì)于強(qiáng)化傳熱具有重要意義。

    圖10 模擬對(duì)象的流線Fig.10 Streamline of simulation object

    為了細(xì)致分析流體在徑向上的流動(dòng),考察了第七個(gè)螺旋片對(duì)應(yīng)區(qū)域的橫截面速度場(chǎng)的變化情況,結(jié)果如圖11所示。z為模型橫截面距入口的長(zhǎng)度,則z/l可便捷地表示橫截面所處位置,比如z/l為49/8就代表第七個(gè)螺旋片軸向1/8長(zhǎng)度位置的橫截面,黑色箭頭矢量表示速度在橫截面上的投影;顏色等值線表示軸向速度大小。

    圖11 第七個(gè)螺旋片不同橫截面的速度場(chǎng)Fig.11 Velocity field of different cross sections of the 7th helical plate

    可以發(fā)現(xiàn)每個(gè)截面上都存在徑向二次流,這是因?yàn)榱黧w在螺旋前進(jìn)的過(guò)程中,軸心附近的流速高于管壁附近的流速,不同流速的流體微團(tuán)的離心力不同,徑向上產(chǎn)生了逆壓梯度,引起流體從軸心向螺旋外側(cè)管壁的流動(dòng),同時(shí)還會(huì)引起流體沿管壁從螺旋外側(cè)向內(nèi)側(cè)的附加流動(dòng),形成一對(duì)反向的雙旋流渦,這對(duì)旋渦又被稱為Dean渦[16]。圖11中螺旋片兩側(cè)各有一對(duì)Dean渦,圖11(a)和圖11(b)中每一對(duì)都是兩個(gè)同向渦組成,這是因?yàn)橄噜彽膬蓚€(gè)螺旋片頭尾呈90°交叉布置,第六個(gè)螺旋片產(chǎn)生的兩對(duì)Dean渦進(jìn)入第七個(gè)螺旋片時(shí),四個(gè)渦被重組,旋向相同的兩個(gè)渦進(jìn)入了第七個(gè)螺旋片的同一側(cè)。由于同向的兩個(gè)渦的動(dòng)能會(huì)相互抵消,可以觀察到:從圖11(a)到圖11(b),渦的強(qiáng)度在逐漸變?nèi)?;圖11(c)時(shí),其中一個(gè)渦消失;圖11(d)時(shí),渦消失的區(qū)域形成了一個(gè)與消失前反向的渦;再?gòu)膱D11(e)到圖11(h),兩對(duì)反向渦不斷發(fā)展、平衡,如此再進(jìn)入第八個(gè)螺旋片。整個(gè)流體區(qū)域的二次流都表現(xiàn)出這種周期往復(fù)的特性。

    3.2 溫度場(chǎng)及傳熱過(guò)程強(qiáng)化

    在二次流的影響下,反應(yīng)器的傳熱過(guò)程也得到了強(qiáng)化,圖12所示為第七個(gè)螺旋片不同橫截面的溫度場(chǎng)情況。顏色等值線表示溫度場(chǎng),線狀的紋理表示速度矢量在橫截面投影的線積分卷積。

    圖12 第七個(gè)螺旋片不同橫截面的溫度場(chǎng)Fig.12 Temperature field of different cross sections of the 7th helical plate

    由圖12可知,熱邊界層很薄,整個(gè)橫截面的絕大部分區(qū)域溫度是均一的,這是由于在二次流的作用下,流動(dòng)主體區(qū)的冷流體不斷侵?jǐn)_邊界層的熱流體,邊界層的熱流體不斷被卷入流動(dòng)主體區(qū)的冷流體,這個(gè)過(guò)程中不僅流體微團(tuán)的位置發(fā)生了交錯(cuò),流體微團(tuán)的溫度也在快速交換和平衡,形成了大面積溫度均一的區(qū)域。

    為了更加直觀地展現(xiàn)SK螺旋結(jié)構(gòu)強(qiáng)化傳熱的能力,以相同邊界條件下無(wú)構(gòu)件空管的溫度場(chǎng)為對(duì)照,三個(gè)截面的中軸線溫度分布如圖13所示,橫坐標(biāo)χ為橫截面徑向上的無(wú)量綱長(zhǎng)度??梢钥闯鲭S著軸向距離的增大,SK結(jié)構(gòu)對(duì)應(yīng)的主流區(qū)溫度快速上升并逐漸趨于均一,而空管對(duì)應(yīng)的主流區(qū)溫度僅緩慢上升,并且徑向分布很不均勻。在該條件下,SK螺旋結(jié)構(gòu)的傳熱系數(shù)為空管的4.2倍。

    圖13 不同橫截面中軸線的溫度分布Fig.13 Temperature distribution of the central axis with different cross sections

    3.3 場(chǎng)協(xié)同效應(yīng)分析

    將對(duì)流傳熱比作有熱源的熱傳導(dǎo),傳熱速率不僅與流體物性和流速有關(guān),還與速度場(chǎng)與溫度梯度的夾角有關(guān)。當(dāng)流體被加熱時(shí),應(yīng)使夾角盡量接近0°,當(dāng)流體被冷卻時(shí)應(yīng)使夾角接近180°,這一思想被稱為“場(chǎng)協(xié)同原理”[17-18]。

    通過(guò)式(4)計(jì)算不同截面的協(xié)同角,結(jié)果如圖14所示。可以看出有較大面積區(qū)域的協(xié)同角偏離90°,這說(shuō)明SK螺旋結(jié)構(gòu)產(chǎn)生的二次流對(duì)于強(qiáng)化傳熱起著重要的作用。觀察圖14(a)~圖14(c)紅色與藍(lán)色區(qū)域的演變,可以發(fā)現(xiàn)綠色區(qū)域占比逐漸增多,這表明流過(guò)第六個(gè)螺旋片的流體二次流在第七個(gè)螺旋片的反旋向?qū)Я飨轮饾u減弱。在第七個(gè)螺旋片的后半段,即圖14(d)~圖14(h)所示的截面位置,圖中紅藍(lán)色區(qū)域的占比明顯增多,這是由于流體在第七個(gè)螺旋片的主導(dǎo)下逐漸形成新的Dean渦,二次流再次發(fā)展增強(qiáng)。

    圖14 不同截面的協(xié)同角分布Fig.14 Synergy angle distribution of different cross sections

    3.4 雷諾數(shù)與普朗特?cái)?shù)對(duì)傳熱的影響

    雷諾數(shù)(Re)和普朗特?cái)?shù)(Pr)是影響對(duì)流傳熱的重要因素,通過(guò)單因素分析的方式考察了兩者對(duì)傳熱過(guò)程的影響。為了綜合考慮傳熱與流動(dòng)阻力兩方面的結(jié)果,引入性能評(píng)價(jià)指標(biāo)(PEC)來(lái)量化等泵功約束條件下強(qiáng)化傳熱的程度[19],計(jì)算方法見(jiàn)式(5)。

    選取等水力學(xué)直徑的中空?qǐng)A管作為參照結(jié)構(gòu),不同Re時(shí)的Nu和f值如圖15所示,對(duì)PEC的影響如圖16所示??梢钥闯?,在其他因素保持不變時(shí),Nu隨著Re的增大而增大,f隨著Re的增大而減小,PEC呈現(xiàn)先增大再減小的規(guī)律。在考察范圍內(nèi),PEC數(shù)值始終大于1,說(shuō)明SK螺旋結(jié)構(gòu)實(shí)現(xiàn)了強(qiáng)化傳熱。

    圖15 Re對(duì)Nu和f的影響Fig.15 The effect of Re on Nu and f

    圖16 Re對(duì)PEC的影響Fig.16 The effect of Re on PEC

    通過(guò)擬合得到:

    Pr對(duì)Nu和f的影響如圖17所示??梢钥吹剑琋u隨Pr的增大而增大,f數(shù)值上與Pr無(wú)關(guān)。擬合得到Nu與Pr的關(guān)系見(jiàn)式(8),Pr項(xiàng)指數(shù)為0.406,與空管在完全湍流區(qū)的Nu實(shí)驗(yàn)關(guān)聯(lián)式的Pr項(xiàng)指數(shù)0.4基本一致。

    圖17 Pr對(duì)Nu和f的影響Fig.17 The effect of Pr on Nu and f

    3.5 螺旋片長(zhǎng)徑比與扭轉(zhuǎn)角對(duì)傳熱的影響

    幾何結(jié)構(gòu)也是影響傳熱過(guò)程的因素之一,對(duì)于SK螺旋結(jié)構(gòu),特征幾何參數(shù)為螺旋片長(zhǎng)徑比(Ar)和扭轉(zhuǎn)角(φ)。不同Ar對(duì)Nu和f的影響如圖18所示,對(duì)PEC的影響如圖19所示。

    圖18 Ar對(duì)Nu和f的影響Fig.18 The effect of Ar on Nu and f

    圖19 Ar對(duì)PEC的影響Fig.19 The effect of Ar on PEC

    隨著Ar的增大,Nu和f都呈現(xiàn)減小的趨勢(shì),這是由于Ar增大會(huì)使得單位軸向長(zhǎng)度螺旋片旋轉(zhuǎn)角度的下降,二次流的強(qiáng)度減弱。PEC則隨著Ar的增大緩慢減小,當(dāng)Ar>2.4時(shí)減小速度變快,因此在實(shí)際應(yīng)用中Ar設(shè)計(jì)值為1.7~2.4較合適。通過(guò)擬合得到如下的關(guān)聯(lián)式:

    不同φ對(duì)Nu和f的影響如圖20所示,對(duì)PEC的影響如圖21所示,其中φ0為180°??梢钥闯鯪u和f隨著φ/φ0的增大而增大,PEC隨φ/φ0的增大而顯著增大,這說(shuō)明增大扭轉(zhuǎn)角是一種強(qiáng)化SK螺旋結(jié)構(gòu)傳熱的理想方式,在實(shí)際應(yīng)用中選擇180°扭轉(zhuǎn)角的螺旋片是合適的。

    圖20 φ對(duì)Nu和f的影響Fig.20 The effect of φ on Nu and f

    圖21 φ對(duì)PEC的影響Fig.21 The effect of φ on PEC

    擬合得到Nu和f與無(wú)因次扭轉(zhuǎn)角的關(guān)系式分別為:

    3.6 計(jì)算關(guān)聯(lián)式

    取3.4和3.5的所有數(shù)據(jù)進(jìn)行擬合,得到SK螺旋結(jié)構(gòu)的Nu及f的總計(jì)算關(guān)聯(lián)式:

    式中:866<Re<4 330,1<Pr<214,1.7<Ar<2.9,0.5<φ/φ0<1。

    圖22和圖23為總關(guān)聯(lián)式計(jì)算得到的值與模擬值的比較。

    圖22 Nu計(jì)算值與擬合值的對(duì)比Fig.22 Comparison of calculated Nu and fitted Nu

    圖23 f計(jì)算值與擬合值的對(duì)比Fig.23 Comparison of calculated f and fitted f

    由圖22和圖23可知,計(jì)算值和擬合值的最大偏差都在10%以內(nèi)。Nu計(jì)算關(guān)聯(lián)式的相關(guān)系數(shù)R2為0.998 1,顯著性檢驗(yàn)中F值為23 089,F(xiàn)0.95=4.07,F(xiàn)10F0.95;f計(jì)算關(guān)聯(lián)式的相關(guān)系數(shù)為0.999 7,顯著性檢驗(yàn)中F值為111 353,F(xiàn)0.95為4.16,F(xiàn)10F0.95。兩個(gè)關(guān)聯(lián)式都有著較好的符合性和預(yù)測(cè)性。

    4 結(jié) 論

    a)基于擴(kuò)大通道尺寸并內(nèi)置三維螺旋內(nèi)構(gòu)件的思路設(shè)計(jì)了一種高通量連續(xù)流反應(yīng)器,通量可達(dá)1000 mL/min,比傳統(tǒng)連續(xù)流微反應(yīng)器高1~2個(gè)數(shù)量級(jí);

    b)在周期性螺旋片的作用下,流動(dòng)徑向上產(chǎn)生大量二次流,這些漩渦顯著了強(qiáng)化傳熱過(guò)程;

    c)熱邊界層很薄,徑向上溫度近乎均一,SK結(jié)構(gòu)的傳熱系數(shù)為空管的4.2倍;

    d)雷諾數(shù)(Re)對(duì)性能評(píng)價(jià)指標(biāo)PEC影響較大,合適的Re為2 000~3 000;

    e)螺旋片長(zhǎng)徑比(Ar)和扭轉(zhuǎn)角(φ)對(duì)PEC影響較大,合適的Ar為1.7~2.4,φ為180°;

    f)反應(yīng)器內(nèi)努塞爾數(shù)(Nu)和摩擦系數(shù)(f)計(jì)算關(guān)聯(lián)式如下:

    式中:866<Re<4 330,1<Pr<214,1.7<Ar<2.9,0.5<φ/φ0<1。

    符號(hào)說(shuō)明

    Ar—— 螺旋片長(zhǎng)徑比vi—— 速度張量,m/s

    cp—— 流體比熱容,J/(kg·K)z—— 模型橫截面距入口的長(zhǎng)度,mm

    f—— 流動(dòng)摩擦系數(shù)θ—— 協(xié)同角,°

    l—— 螺片的長(zhǎng)度,mmλ—— 流體熱導(dǎo)率,W/(m·K)

    Nu—— 努塞爾數(shù)μ—— 流體黏度,Pa·s

    Pr—— 普朗特?cái)?shù)ρ—— 流體密度,kg/m3

    Re—— 雷諾數(shù)φ—— 扭轉(zhuǎn)角,°

    T—— 流體溫度,Kχ—— 橫截面徑向上的無(wú)量綱長(zhǎng)度

    猜你喜歡
    連續(xù)流徑向反應(yīng)器
    聚丙烯環(huán)管反應(yīng)器升密操作與控制
    云南化工(2021年11期)2022-01-12 06:06:36
    淺探徑向連接體的圓周運(yùn)動(dòng)
    RN上一類Kirchhoff型方程徑向?qū)ΨQ正解的存在性
    基于PID+前饋的3MN徑向鍛造機(jī)控制系統(tǒng)的研究
    一類無(wú)窮下級(jí)整函數(shù)的Julia集的徑向分布
    EGSB反應(yīng)器處理阿維菌素廢水
    上旋流厭氧反應(yīng)器在造紙廢水處理中的應(yīng)用
    硅酸鋅催化臭氧氧化凈水效能連續(xù)流實(shí)驗(yàn)研究
    TMT公司生產(chǎn)管理存在問(wèn)題及精益生產(chǎn)管理改進(jìn)措施探析
    費(fèi)托合成微反應(yīng)器研究進(jìn)展
    王馨瑶露胸无遮挡在线观看| 最近的中文字幕免费完整| 秋霞在线观看毛片| 美女视频免费永久观看网站| 丝袜喷水一区| videosex国产| 美女脱内裤让男人舔精品视频| 日韩大码丰满熟妇| 满18在线观看网站| 亚洲精品日韩在线中文字幕| 久久免费观看电影| 少妇被粗大猛烈的视频| 丰满迷人的少妇在线观看| 大片免费播放器 马上看| 国产日韩欧美视频二区| 嫩草影视91久久| 国产成人午夜福利电影在线观看| 这个男人来自地球电影免费观看 | 精品亚洲成国产av| 黑丝袜美女国产一区| 一区二区三区精品91| 大片电影免费在线观看免费| 水蜜桃什么品种好| 国产精品久久久av美女十八| 亚洲成人av在线免费| 欧美日韩成人在线一区二区| 日韩制服骚丝袜av| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 国产精品久久久av美女十八| 你懂的网址亚洲精品在线观看| 国产精品免费视频内射| 如何舔出高潮| 国产在视频线精品| 成年动漫av网址| 亚洲五月色婷婷综合| 丁香六月欧美| 永久免费av网站大全| 亚洲成色77777| 捣出白浆h1v1| 国产毛片在线视频| 亚洲成人国产一区在线观看 | 9191精品国产免费久久| 自拍欧美九色日韩亚洲蝌蚪91| 一本—道久久a久久精品蜜桃钙片| 女性生殖器流出的白浆| 国产精品国产av在线观看| 午夜激情av网站| 久久久久久免费高清国产稀缺| 国产乱来视频区| 一级爰片在线观看| 在线天堂最新版资源| 成人黄色视频免费在线看| 成人三级做爰电影| 丝袜美足系列| 99久久人妻综合| 一区福利在线观看| 久久人妻熟女aⅴ| 午夜日本视频在线| 欧美日韩亚洲国产一区二区在线观看 | 国产精品久久久av美女十八| 国产成人精品久久久久久| 在线观看免费午夜福利视频| 色视频在线一区二区三区| 777久久人妻少妇嫩草av网站| av在线app专区| 亚洲在久久综合| 伊人亚洲综合成人网| 女人精品久久久久毛片| 最近手机中文字幕大全| 日本91视频免费播放| 九草在线视频观看| 观看av在线不卡| 久久久久人妻精品一区果冻| 少妇人妻精品综合一区二区| 晚上一个人看的免费电影| 午夜日韩欧美国产| 青春草视频在线免费观看| 久热这里只有精品99| 免费少妇av软件| 午夜福利乱码中文字幕| 十八禁高潮呻吟视频| 日韩中文字幕欧美一区二区 | 欧美人与性动交α欧美精品济南到| 午夜av观看不卡| 久久性视频一级片| 哪个播放器可以免费观看大片| 亚洲精品久久久久久婷婷小说| 久久精品亚洲av国产电影网| 美女主播在线视频| 亚洲精品视频女| 亚洲精品久久久久久婷婷小说| 亚洲av中文av极速乱| 免费在线观看黄色视频的| 亚洲成国产人片在线观看| www.熟女人妻精品国产| 色精品久久人妻99蜜桃| 99久久人妻综合| 婷婷色av中文字幕| 国产无遮挡羞羞视频在线观看| 老汉色∧v一级毛片| 99久久人妻综合| 日韩视频在线欧美| 免费少妇av软件| 热re99久久国产66热| 色吧在线观看| 伦理电影免费视频| 午夜福利影视在线免费观看| 黄网站色视频无遮挡免费观看| 男人舔女人的私密视频| 久久久精品免费免费高清| 国产在线免费精品| 亚洲国产欧美一区二区综合| 91aial.com中文字幕在线观看| 成年女人毛片免费观看观看9 | 天天躁狠狠躁夜夜躁狠狠躁| 伊人亚洲综合成人网| 久久午夜综合久久蜜桃| 亚洲精品国产av蜜桃| 国产成人精品在线电影| 熟女少妇亚洲综合色aaa.| 亚洲欧美精品综合一区二区三区| 国产成人一区二区在线| 日本vs欧美在线观看视频| 日本欧美国产在线视频| 日韩欧美一区视频在线观看| 国产深夜福利视频在线观看| 成人免费观看视频高清| 国产一区亚洲一区在线观看| 性色av一级| 国产成人精品福利久久| 99久久精品国产亚洲精品| 国语对白做爰xxxⅹ性视频网站| 国产成人系列免费观看| av在线观看视频网站免费| 黄色怎么调成土黄色| 乱人伦中国视频| 中文字幕人妻丝袜制服| av一本久久久久| 热re99久久精品国产66热6| www.av在线官网国产| 一级毛片 在线播放| av女优亚洲男人天堂| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 少妇被粗大猛烈的视频| av一本久久久久| 老司机影院成人| 老司机靠b影院| 国产免费又黄又爽又色| 一级毛片黄色毛片免费观看视频| 男女高潮啪啪啪动态图| 高清在线视频一区二区三区| 久久天躁狠狠躁夜夜2o2o | 国产乱来视频区| 美女高潮到喷水免费观看| 十分钟在线观看高清视频www| 亚洲欧美中文字幕日韩二区| 岛国毛片在线播放| svipshipincom国产片| 妹子高潮喷水视频| 亚洲七黄色美女视频| www.熟女人妻精品国产| 乱人伦中国视频| 亚洲精品aⅴ在线观看| 黑人猛操日本美女一级片| 飞空精品影院首页| 欧美精品人与动牲交sv欧美| 最近中文字幕高清免费大全6| 久久久久久人妻| 国产一区二区三区综合在线观看| 日韩大码丰满熟妇| 国产麻豆69| 青青草视频在线视频观看| 国产精品熟女久久久久浪| 男人舔女人的私密视频| 亚洲男人天堂网一区| 成人毛片60女人毛片免费| 色婷婷av一区二区三区视频| 观看美女的网站| 久热爱精品视频在线9| 国产成人系列免费观看| 欧美日本中文国产一区发布| a级毛片在线看网站| 久久久精品国产亚洲av高清涩受| 久久久久久久久免费视频了| 国产 精品1| 国产精品久久久av美女十八| 国产av码专区亚洲av| 在线观看免费午夜福利视频| 亚洲av电影在线观看一区二区三区| 亚洲自偷自拍图片 自拍| 黑人巨大精品欧美一区二区蜜桃| 国产毛片在线视频| 欧美精品人与动牲交sv欧美| 最近最新中文字幕免费大全7| 久久青草综合色| 国产亚洲一区二区精品| 赤兔流量卡办理| 亚洲av日韩在线播放| 又大又黄又爽视频免费| 日韩成人av中文字幕在线观看| 侵犯人妻中文字幕一二三四区| 在线观看一区二区三区激情| 国产熟女欧美一区二区| 国产欧美日韩综合在线一区二区| 午夜福利视频精品| 欧美亚洲日本最大视频资源| 国产免费又黄又爽又色| 国产熟女午夜一区二区三区| 少妇 在线观看| 五月开心婷婷网| 日韩av免费高清视频| 国产欧美日韩一区二区三区在线| 亚洲精品自拍成人| 2021少妇久久久久久久久久久| 一级爰片在线观看| 99热国产这里只有精品6| 女性被躁到高潮视频| 亚洲精品美女久久av网站| 又粗又硬又长又爽又黄的视频| 999精品在线视频| 一级毛片黄色毛片免费观看视频| 丁香六月欧美| 成人午夜精彩视频在线观看| 精品久久蜜臀av无| 伊人久久大香线蕉亚洲五| 久久天躁狠狠躁夜夜2o2o | 欧美成人精品欧美一级黄| 男女之事视频高清在线观看 | 亚洲欧美一区二区三区久久| 99国产精品免费福利视频| 国产成人免费观看mmmm| 中文字幕制服av| 亚洲七黄色美女视频| 国产av国产精品国产| 精品国产一区二区三区四区第35| 69精品国产乱码久久久| 韩国av在线不卡| 叶爱在线成人免费视频播放| av国产精品久久久久影院| 久久国产精品男人的天堂亚洲| 日本午夜av视频| a级毛片黄视频| 久久久国产精品麻豆| 久久影院123| 国产亚洲精品第一综合不卡| 欧美在线一区亚洲| 亚洲伊人久久精品综合| 波野结衣二区三区在线| 又大又爽又粗| 大片电影免费在线观看免费| 超碰成人久久| tube8黄色片| 香蕉丝袜av| 少妇的丰满在线观看| 欧美国产精品va在线观看不卡| 免费少妇av软件| 久久精品国产亚洲av高清一级| 中文字幕色久视频| 美女福利国产在线| 18禁观看日本| 欧美xxⅹ黑人| av片东京热男人的天堂| 人妻人人澡人人爽人人| 精品一区二区三卡| 成人毛片60女人毛片免费| www.精华液| 热re99久久精品国产66热6| 成年动漫av网址| 伊人亚洲综合成人网| 欧美少妇被猛烈插入视频| 伊人久久国产一区二区| 国产国语露脸激情在线看| 青春草国产在线视频| 国产成人精品久久二区二区91 | 婷婷色av中文字幕| 色精品久久人妻99蜜桃| 午夜激情av网站| 国产有黄有色有爽视频| netflix在线观看网站| 人妻一区二区av| 免费久久久久久久精品成人欧美视频| 一级爰片在线观看| 国产成人欧美在线观看 | 午夜免费观看性视频| 亚洲国产av新网站| 建设人人有责人人尽责人人享有的| 国产福利在线免费观看视频| 在现免费观看毛片| 99国产综合亚洲精品| 综合色丁香网| 另类亚洲欧美激情| www.精华液| 新久久久久国产一级毛片| 亚洲熟女精品中文字幕| 99热全是精品| 电影成人av| 午夜福利影视在线免费观看| 波野结衣二区三区在线| 嫩草影院入口| 伦理电影大哥的女人| 国产精品女同一区二区软件| 国产成人午夜福利电影在线观看| 午夜久久久在线观看| 日本91视频免费播放| 日韩不卡一区二区三区视频在线| 晚上一个人看的免费电影| 久久天躁狠狠躁夜夜2o2o | 一级爰片在线观看| 狂野欧美激情性xxxx| av线在线观看网站| 波多野结衣一区麻豆| 侵犯人妻中文字幕一二三四区| 在线看a的网站| 国产成人a∨麻豆精品| 欧美人与善性xxx| 高清av免费在线| 国产野战对白在线观看| 操美女的视频在线观看| 亚洲av综合色区一区| 欧美日韩一区二区视频在线观看视频在线| 精品一区在线观看国产| 黄色 视频免费看| 午夜福利免费观看在线| 欧美亚洲 丝袜 人妻 在线| 亚洲欧洲日产国产| 国产野战对白在线观看| 桃花免费在线播放| 欧美 日韩 精品 国产| 亚洲第一区二区三区不卡| 亚洲国产中文字幕在线视频| 国产乱来视频区| 91精品三级在线观看| 午夜av观看不卡| 久久99热这里只频精品6学生| 国产爽快片一区二区三区| 欧美日韩综合久久久久久| 日韩中文字幕视频在线看片| 欧美日韩精品网址| 久久av网站| 精品一区二区三区四区五区乱码 | 可以免费在线观看a视频的电影网站 | 免费看av在线观看网站| 午夜福利免费观看在线| 久久久久国产一级毛片高清牌| 一区二区三区激情视频| 亚洲成色77777| av网站在线播放免费| 亚洲欧美激情在线| 最近中文字幕2019免费版| 国产一区亚洲一区在线观看| 天堂俺去俺来也www色官网| 人人妻人人澡人人爽人人夜夜| 精品午夜福利在线看| 国产精品一区二区精品视频观看| 80岁老熟妇乱子伦牲交| 日本午夜av视频| 国产精品一二三区在线看| 天天躁夜夜躁狠狠久久av| 日日摸夜夜添夜夜爱| 人妻一区二区av| 中文字幕最新亚洲高清| av在线老鸭窝| 777米奇影视久久| 亚洲国产毛片av蜜桃av| 成人国语在线视频| 在线观看一区二区三区激情| 精品久久久久久电影网| 熟妇人妻不卡中文字幕| 日韩制服骚丝袜av| 亚洲国产精品国产精品| 日日撸夜夜添| 丝袜喷水一区| 午夜激情久久久久久久| 男女无遮挡免费网站观看| 国产亚洲欧美精品永久| 国产免费又黄又爽又色| 日韩av免费高清视频| 午夜免费观看性视频| 国产激情久久老熟女| 午夜影院在线不卡| 一边亲一边摸免费视频| avwww免费| 高清av免费在线| 日韩 欧美 亚洲 中文字幕| 精品久久蜜臀av无| 免费少妇av软件| 大码成人一级视频| 国产淫语在线视频| 欧美中文综合在线视频| avwww免费| 午夜久久久在线观看| 日韩精品免费视频一区二区三区| 黑人猛操日本美女一级片| 人妻 亚洲 视频| 日韩av免费高清视频| 日韩精品免费视频一区二区三区| 黑人猛操日本美女一级片| 考比视频在线观看| 成人漫画全彩无遮挡| 宅男免费午夜| 国产成人av激情在线播放| 人人妻人人爽人人添夜夜欢视频| 日韩成人av中文字幕在线观看| 亚洲在久久综合| 日本爱情动作片www.在线观看| 老汉色av国产亚洲站长工具| 国产av国产精品国产| 欧美激情极品国产一区二区三区| 丁香六月欧美| 国产高清不卡午夜福利| 2021少妇久久久久久久久久久| 制服丝袜香蕉在线| 国产片特级美女逼逼视频| 国产成人精品无人区| 亚洲欧美中文字幕日韩二区| 一级片免费观看大全| 水蜜桃什么品种好| 免费人妻精品一区二区三区视频| 国产无遮挡羞羞视频在线观看| 男人操女人黄网站| 狠狠精品人妻久久久久久综合| bbb黄色大片| 99久久人妻综合| 国产97色在线日韩免费| 男女无遮挡免费网站观看| 国产精品 欧美亚洲| 国产成人午夜福利电影在线观看| 国产成人av激情在线播放| 国产一区二区激情短视频 | 满18在线观看网站| 国产亚洲午夜精品一区二区久久| 九色亚洲精品在线播放| 久久久国产精品麻豆| 久久午夜综合久久蜜桃| 99久久人妻综合| 亚洲精品久久久久久婷婷小说| av不卡在线播放| 国产欧美亚洲国产| 黑人猛操日本美女一级片| 无遮挡黄片免费观看| 精品国产乱码久久久久久小说| 97人妻天天添夜夜摸| 丝袜美腿诱惑在线| 啦啦啦在线免费观看视频4| 两个人免费观看高清视频| 亚洲成国产人片在线观看| 999久久久国产精品视频| 欧美日韩成人在线一区二区| 美女国产高潮福利片在线看| 国产精品av久久久久免费| 日韩成人av中文字幕在线观看| 赤兔流量卡办理| 大片电影免费在线观看免费| 国产亚洲精品第一综合不卡| 观看av在线不卡| 国产精品蜜桃在线观看| 日韩制服骚丝袜av| 亚洲熟女毛片儿| av网站免费在线观看视频| 国产精品一二三区在线看| 日本一区二区免费在线视频| 久久国产精品男人的天堂亚洲| 精品一区在线观看国产| 亚洲中文av在线| 99久国产av精品国产电影| 这个男人来自地球电影免费观看 | 99久久精品国产亚洲精品| 美女午夜性视频免费| 男女国产视频网站| 欧美日本中文国产一区发布| 亚洲欧美清纯卡通| 亚洲激情五月婷婷啪啪| 悠悠久久av| 最近最新中文字幕免费大全7| 美女中出高潮动态图| 黑人欧美特级aaaaaa片| 大陆偷拍与自拍| 黄色视频在线播放观看不卡| 观看av在线不卡| 午夜福利视频在线观看免费| 久久久久精品国产欧美久久久 | 亚洲欧美激情在线| 久久久久精品久久久久真实原创| 亚洲国产毛片av蜜桃av| 久久影院123| 人人妻人人澡人人爽人人夜夜| 在线 av 中文字幕| 久久性视频一级片| videos熟女内射| 在线观看三级黄色| 亚洲成人免费av在线播放| 少妇人妻久久综合中文| 伊人久久国产一区二区| 菩萨蛮人人尽说江南好唐韦庄| 午夜精品国产一区二区电影| 国产极品天堂在线| 最近2019中文字幕mv第一页| 久久免费观看电影| 亚洲专区中文字幕在线 | 国产又色又爽无遮挡免| 毛片一级片免费看久久久久| 永久免费av网站大全| 美女视频免费永久观看网站| 日韩中文字幕视频在线看片| 亚洲,欧美,日韩| 精品久久久精品久久久| 国产精品久久久久久人妻精品电影 | 19禁男女啪啪无遮挡网站| 日韩欧美精品免费久久| 1024视频免费在线观看| 女人被躁到高潮嗷嗷叫费观| 国产极品粉嫩免费观看在线| 国产野战对白在线观看| 又粗又硬又长又爽又黄的视频| 99国产精品免费福利视频| 最近手机中文字幕大全| 国产女主播在线喷水免费视频网站| 婷婷色av中文字幕| 久久久久国产一级毛片高清牌| 深夜精品福利| 丰满少妇做爰视频| 亚洲国产欧美一区二区综合| 最近2019中文字幕mv第一页| 午夜老司机福利片| 欧美变态另类bdsm刘玥| 午夜福利视频精品| 亚洲激情五月婷婷啪啪| 亚洲伊人色综图| 夫妻午夜视频| 久久久久久久久久久免费av| 亚洲欧美精品综合一区二区三区| 国产麻豆69| 国产精品国产av在线观看| 王馨瑶露胸无遮挡在线观看| 18在线观看网站| av线在线观看网站| 一级毛片电影观看| www.av在线官网国产| 在现免费观看毛片| 最黄视频免费看| 欧美av亚洲av综合av国产av | 亚洲av日韩精品久久久久久密 | 国产一级毛片在线| 又粗又硬又长又爽又黄的视频| 中文字幕人妻丝袜一区二区 | av.在线天堂| 国产亚洲av高清不卡| 少妇人妻久久综合中文| 黄片小视频在线播放| 亚洲国产欧美网| 亚洲精品久久久久久婷婷小说| videos熟女内射| 亚洲情色 制服丝袜| 国产伦人伦偷精品视频| 黄色 视频免费看| 2018国产大陆天天弄谢| 成人免费观看视频高清| 亚洲精品成人av观看孕妇| 日日撸夜夜添| 亚洲久久久国产精品| 在线 av 中文字幕| bbb黄色大片| 午夜免费男女啪啪视频观看| 日本av手机在线免费观看| av有码第一页| 国产精品久久久久成人av| 亚洲少妇的诱惑av| 如日韩欧美国产精品一区二区三区| 在线亚洲精品国产二区图片欧美| 狂野欧美激情性xxxx| 亚洲视频免费观看视频| 一二三四中文在线观看免费高清| 免费女性裸体啪啪无遮挡网站| 一级a爱视频在线免费观看| 欧美日韩成人在线一区二区| 黄色 视频免费看| 成人国语在线视频| 久久毛片免费看一区二区三区| 天美传媒精品一区二区| 性高湖久久久久久久久免费观看| 99九九在线精品视频| a级毛片在线看网站| 欧美 日韩 精品 国产| 亚洲成国产人片在线观看| 亚洲人成网站在线观看播放| 午夜福利视频精品| 永久免费av网站大全| 欧美亚洲日本最大视频资源| 亚洲一级一片aⅴ在线观看| 国产午夜精品一二区理论片| 亚洲,一卡二卡三卡| 国产福利在线免费观看视频| 亚洲精品中文字幕在线视频| 色婷婷av一区二区三区视频| 人妻 亚洲 视频| 国产1区2区3区精品| 欧美最新免费一区二区三区| 老司机在亚洲福利影院| 国产精品无大码| 精品亚洲成国产av| 制服人妻中文乱码| 久久鲁丝午夜福利片| 久久韩国三级中文字幕| 男的添女的下面高潮视频| 悠悠久久av| 各种免费的搞黄视频| 宅男免费午夜| 大片免费播放器 马上看| 久久免费观看电影| 熟妇人妻不卡中文字幕| 午夜久久久在线观看| 五月开心婷婷网| 丝袜脚勾引网站|