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

    變截面微管道中高zeta 勢下冪律流體的旋轉(zhuǎn)電滲滑移流動*

    2022-07-22 05:56:36張?zhí)禅?/span>任美蓉崔繼峰陳小剛王怡丹
    物理學(xué)報 2022年13期
    關(guān)鍵詞:磁場

    張?zhí)禅?任美蓉 崔繼峰 陳小剛 王怡丹

    (內(nèi)蒙古工業(yè)大學(xué)理學(xué)院,呼和浩特 010051)

    1 引言

    微生物傳感器、微機(jī)電微生物傳感器系統(tǒng)(MEMS) 等微流控器件廣泛應(yīng)用于生物醫(yī)學(xué)和生化分析領(lǐng)域,是研究基本物理過程的有力工具之一[1].在這些過程中,壓力梯度、電場、磁場或它們的適當(dāng)組合是常用的驅(qū)動機(jī)制.與以往單一的壓力驅(qū)動流動模式相比,電滲和電磁驅(qū)動機(jī)制近年來受到越來越多的關(guān)注.

    電滲流 (EOF) 是由于電場在管道中的作用以及雙電層 (EDL) 在通道壁上的存在而引起的流動.當(dāng)微通道的固體表面與電解質(zhì)溶液接觸時,它們之間的化學(xué)平衡導(dǎo)致界面獲得凈固定電荷,之后,電解質(zhì)溶液中的自由離子聚集在溶液表面附近的區(qū)域,形成雙電層 (EDL).近年來,人們逐漸認(rèn)識到非牛頓流體在各種工業(yè)和工程應(yīng)用中比牛頓流體更為迫切.利用電滲流的微流控器件已經(jīng)被證明在醫(yī)學(xué)應(yīng)用中是有效和實用的.冪律流體模型是最簡單、最常見的一種描述非牛頓流體的方法,受到該領(lǐng)域研究者的特別關(guān)注.例如,Patel 等[2]研究了通過聚電解質(zhì)接枝 (軟) 窄通道旋轉(zhuǎn)的電滲流中冪律流體的流動行為;Srinivas[3]研究了由冪律模型控制的非牛頓流體的二維電滲流,用于圓形和橢圓形微通道;Nekoubin[4]對非牛頓流體在彎曲矩形微通道中的電滲流動進(jìn)行了數(shù)值研究,發(fā)現(xiàn)在剪切稀化流體流動中可以通過提高 zeta 電勢來顯著提高體積流量和混合性能;Ba?os 等[5]研究了由于服從冪律的振蕩電滲流而引起的流體動力學(xué);Ba?os 等[6]研究了冪律流體中 Navier 滑移條件的影響下的振蕩電滲流 (OEOF),并發(fā)現(xiàn)了 (OEOF)在微流控設(shè)備 (例如微混合器) 中很有用;Jiang和Qi[7]研究了具有 Navier 滑移邊界微平行管道內(nèi) Eyring 流體在外加電場力和壓力作用下的電滲流動.對于微流控器件,旋轉(zhuǎn)電滲流 (EOF) 可以提高不同液體的混合效率,從而縮短微流體的混合時間.旋轉(zhuǎn)電滲流中渦流的產(chǎn)生取決于微管道的幾何形狀;Ajdari[8]報道了微通道中正弦表面勢電滲流中渦流的產(chǎn)生;Chang 和 Wang[9]考慮了旋轉(zhuǎn)的EOF,發(fā)現(xiàn)旋轉(zhuǎn)降低了外加電場方向上的 EO 流速;Song 等[10]研究了三級流體在兩平行板通道中的旋轉(zhuǎn)電滲流動;Shit 等[11]研究了滑動速度對非均勻微通道中旋轉(zhuǎn)電滲流的影響;劉全生等[12]研究了平行微管道中Jeffrey 流體的非定常電滲流動;段娟等[13]研究了有限長微擴(kuò)張管道內(nèi)冪律流體在電滲驅(qū)動下非穩(wěn)態(tài)的流動特性.

    除上述電滲流外,磁流體力學(xué) (MHD) 流動也因其在工程和醫(yī)學(xué)領(lǐng)域的應(yīng)用前景而備受關(guān)注[14].為了實現(xiàn)更有效的流動控制,電磁流體力學(xué) (EMHD)流動也受到了廣泛的關(guān)注,即對導(dǎo)電流體施加外加電場和磁場.電場和橫向磁場的相互作用可以產(chǎn)生洛倫茲力.大量的理論和實驗工作在文獻(xiàn)中可用于分析EMHD 流動的行為.Jian 和 Chang[15]獲得了非均勻磁場影響下 EMHD 速度分布的近似解析解;Xie 和 Jian[16]研究了通過窄微通道的冪律流體的旋轉(zhuǎn)電磁流體力學(xué) (EMHD) 流,結(jié)果可以用作開發(fā)精巧而有效的電磁設(shè)備的初步藍(lán)圖,用于涉及流量控制或物質(zhì)分離的應(yīng)用;Habib 等[17]研究了冪律流體的非定常層流自由對流的熵產(chǎn)生分析,流體流存在熱輻射和磁場;Sarkar 和Ganguly[18]研究了在疊加磁場的存在下,電動調(diào)制的軸向壓力驅(qū)動的冪律流體通過微通道的傳輸;Yang 等[19]討論了在平行微通道中的電磁流體力學(xué) (EMHD) 電滲流 (EOF) 和不可壓縮的三級流體的熵產(chǎn)生;Xie和 Jian[20]進(jìn)行了兩層通過微平行通道的磁流體動力電滲流的熵產(chǎn)生分析,并發(fā)現(xiàn)該發(fā)明可以用于設(shè)計有效的熱微設(shè)備.

    對于冪律流體在高 zeta 勢、滑移邊界條件下對微流控器件 EOF 的影響研究已經(jīng)進(jìn)行了很多,但很少有對具有旋轉(zhuǎn)框架的非牛頓流體進(jìn)行同樣的科學(xué)研究.Wang 等[21]研究了高zeta 勢下冪律流體在非均勻微通道中的旋轉(zhuǎn)電滲流動 (EOF) 受滑移速度的影響;Xie 和 Jian[22]討論了冪律流體在微通道中高 zeta 勢下的旋轉(zhuǎn)電滲流動,利用非線性 Poisson-Boltzmann 模型得到了雙電層電勢分布的解析解.

    受上述研究的啟發(fā),本研究旨在研究高 zeta勢下具有滑移邊界條件的冪律流體在壁面周期性變化的變截面微管道中外加磁場的旋轉(zhuǎn)電滲流.利用Navier-Stokes 方程和非線性 Poisson-Boltzmann模型,描述電勢的分布和雙層結(jié)構(gòu),并用有限差分法數(shù)值計算外加磁場的旋轉(zhuǎn)電滲流的速度分布.討論相關(guān)物理參數(shù)對冪律流體無量綱速度的影響.

    2 模型描述

    如圖1 所示,考慮不可壓縮冪律流體在有限長L的變截面對稱微管道中的非定常外加磁場的旋轉(zhuǎn)電滲流.假設(shè)變截面微管道高度z′=±h′(x′) 呈余弦變化,其微管道壁的表達(dá)式為

    圖1 變截面微通道中流體流動示意圖Fig.1.Schematic view of the flow in a variable cross-section microchannel.

    其中a為微管道邊界高度的振幅.

    使用以下坐標(biāo)變換:

    (1)式轉(zhuǎn)換為以下形式:

    其中 (x′,y′,z′) 為笛卡爾坐標(biāo)系;H為管道固定高度;a為微管道壁的振幅;L是微管道的有限長度,其長度遠(yuǎn)大于板間高度邊界的高度.

    假設(shè)微管道系統(tǒng)中充滿了電解質(zhì)溶液,具有介電常數(shù)ε且微管道系統(tǒng)圍繞z軸旋轉(zhuǎn)的角速Ω=(0,0,Ω),假設(shè)流體在微通道中的流動繞x′軸對稱,并沿x′方向流動.流體同時受到電場E=(Ex,Ey,0) 和磁場B=(0,0,Bz) 的作用,使微管道壁上均勻地充滿了 zeta 電勢ψω,根據(jù)靜電學(xué)理論,靜電荷密度ρf與電勢ψ滿足 Poisson-Boltzmann方程,即

    由方程(4)和方程(5) 聯(lián)立可以直接得到:

    相應(yīng)的邊界條件可以表示為

    方程 (6) 在無Debye-Hückel 線性近似下,根據(jù)邊界條件(7)式和(8)式,方程(6)的解析解如下[21]:

    其中ψω為壁面電勢;且κ為雙電層厚度的倒數(shù).n0為電解質(zhì)離子密度,z0為離子的化合價,e為單位電荷,kB為Boltzmann常數(shù),T是絕對溫度.

    整個系統(tǒng)圍繞z軸旋轉(zhuǎn),在非慣性坐標(biāo)系下水平旋轉(zhuǎn)的流體層中,所有的速度都是相對速度,對于不可壓縮流體,給出了連續(xù)性方程:

    其中V是外加磁場的電滲流的速度且V=(u,v,0) .在具有恒定角速度的旋轉(zhuǎn)系統(tǒng)中,柯西動量方程表示為

    其中τ是應(yīng)力張量;ρ為冪律流體的密度;t為時間;P為離心力修正的壓力,P=p-p|Ω×r|2/2,r=(x,y,z) ;p為壓強(qiáng);f為體積力,等于電場 力ρe(x)E和洛倫茲力J×B的和,即

    其中電流J滿足歐姆定律:

    在典型的微通道流動中,磁雷諾數(shù)的量級約為10—5,因此感應(yīng)磁場可以忽略不計,這也意味著感應(yīng)磁場與流速無關(guān).

    對于冪律流體[23],動力黏度μ與剪切變形速率Γ相關(guān)表達(dá)式為

    這里,η是黏性系數(shù),n是冪律流體的行為指數(shù).根據(jù)n的取值,冪律流體可分為三類:1) 膨脹塑性流體 (n>1),其黏度隨剪切速率的增大而增大;2)假塑性流體 (n< 1),其黏度隨剪切速率的增大而減小;3)n=1 時,得到的是牛頓流體.由于黏性冪律流體的動態(tài)黏度和剪切變形速率,黏性應(yīng)力張量和應(yīng)變張量之間的關(guān)系可以表示為

    其中?V是速度梯度的張量;?V ′是它的轉(zhuǎn)置.

    利用文獻(xiàn)中Chang 和Wang[9]以及Xie 和Jian[22]的方法,可以找到u=u(z,t),v=v(z,t) 形式的解.在這種情況下,連續(xù)性方程 (10) 將自動滿足,冪律流體的動力粘度可簡化為

    則剪切應(yīng)力可以寫成:

    因此,柯西動量方程可以簡化為

    特別地,n=1 時,方程 (20)和方程(21) 中的μ(Δ)表示旋轉(zhuǎn)環(huán)境下電動勢驅(qū)動的牛頓流體流動.考慮了Navier 滑移邊界條件,給出了相應(yīng)的初始條件和邊界條件:

    其中b是滑移長度參數(shù).

    為了方便計算,引入無量綱變量:

    這里Ueo=-εExψω/η;Ha稱為哈特曼數(shù),表示磁力和黏滯力的比值;S是一個無量綱參數(shù),用于估計電場強(qiáng)度;K為電動寬度,被定義為微通道半高H與EDL 厚度κ-1的比值.

    將方程 (25) 代入方程 (17) 中,無量綱后冪律流體的動力黏度變?yōu)?/p>

    同樣將方程 (25) 代入方程 (9) 中,得到無量綱后的Poisson-Boltzmann 方程為

    相應(yīng)的邊界條件(7)式和(8)式化為

    因此,無量綱后的柯西動量方程(20)和方程(21) 轉(zhuǎn)化為

    Navier 滑移邊界條件及相應(yīng)的初始條件和邊界條件(22)式、(23)式和(24)式可轉(zhuǎn)化為

    3 數(shù)值算法

    利用有限差分法求解外加磁場的旋轉(zhuǎn)系統(tǒng)中冪律流體的旋轉(zhuǎn)電滲流,考慮域 (tˉ,zˉ) 屬于 [0,T]×[0,h] 速度分量u和v在離散后的網(wǎng)格上的任意一點 (i,m) 上的表達(dá)式為

    考慮方程 (30) 和 方程(31) 是拋物型方程,對時間的偏導(dǎo)數(shù)進(jìn)行一階向前差分,其表達(dá)式為

    柯西動量方程 (30) 和方程 (31) 最后的差分格式為

    離散化后,相應(yīng)的初始條件和邊界條件為

    4 結(jié)果和討論

    本文研究中,對于所有計算,空間步長和時間步長選取分別為 Δz=10-2和 Δt=5×10-6;在目前的計算中,典型微通道流動的參數(shù)值如下:電解質(zhì)溶液的介電常數(shù)ε=709×10-12F/m ;半高H=100 μm ;電解質(zhì)的濃度n0=NAc,其中c=10-5mol/L為摩爾濃度,NA=6.02×1023/mol ;ρ=1.0×103kg/m3;離子價態(tài)z0=1 ;e=1.6×10-19C ;kB=1.38×10-23J/K ;電場強(qiáng) 度Ex=104V/m ;絕對溫度T=293 K ;η=0.9×10-3;σe的變化范圍為 2 .2×10-4—106S/m ;磁場強(qiáng)度B的取值范圍為 0 .02—1 T ;交變電 場強(qiáng)度Ey的取值范圍為0—1 V/m .

    為了驗證本研究數(shù)值方法的正確性,當(dāng)微管道高度波動變化的振幅a=0 時,微管道的形狀為平行管道;當(dāng)n=1 和β=0 時(即對于牛頓流體和無滑移邊界條件),和K=30,Ω=100 rad/s,=1 V的情況下,給出了相應(yīng)的速度分量u,v的數(shù)值結(jié)果與Chang 和Wang[9]的解析結(jié)果如圖2 所示.結(jié)果表明,它們之間的一致性很好.因此,本文所建立的數(shù)值方法是有效的,可用于計算冪律流體在n≠1 時的外加磁場的旋轉(zhuǎn)電滲速度分布.

    圖2 目前數(shù)值解與Chang 和Wang[9]解析解的比較,其中 β=0,K=30,Ω=100 rad/s,=1 V,a=0,Ha=0,S=0Fig.2.Comparison of the current numerical solution with the analytical solution of Chang and Wang [9],β=0,K=30,Ω=100 rad/s,=1 V,a=0,Ha=0,S=0 .

    如圖3 和圖4 所示為不同的冪律流體行為指數(shù)n(0.4,0.5,0.8,1,1.2)的無量綱外加磁場旋轉(zhuǎn)電滲流速度的數(shù)值變化.圖3 和圖4 分別呈現(xiàn)了滑移參數(shù)β=0 和β/=0 (即,無滑移邊界條件和存在滑移邊界條件) 時冪律流體行為指數(shù)n對x方向和y方向的速度刨面的影響.從圖3(a)和圖4(a)可以看出,對于假塑性流體 (n<1),它們的形狀更接近于插頭狀的外加磁場電滲流速度剖面,對于膨脹塑性流體 (n>1),它們更接近于拋物線形的速度剖面.除墻體上的滑移距離外,這與Wang 等[21]與 Xie 和 Jian[22]的研究結(jié)果一致.圖3(a)和圖4(a)中還可以看出,x方向的速度大小隨著冪律行為指數(shù)n的減小而增大,其原因在于冪律行為指數(shù)的增加導(dǎo)致流體的粘度較大,這使得在微通道中流動更加困難.圖3(a)和圖4(a)的對比可以得出,滑移邊界條件存在時對壁面處流體黏性的減小,也導(dǎo)致相比于不存在滑移條件時,到達(dá)了更大的速度最大值和距離壁面更短的距離.同樣的現(xiàn)象也可以在圖3(b)和圖4(b)中找到.

    圖3 當(dāng)無滑移邊界條件時,冪律流體行為指數(shù)n 對外加磁場的旋轉(zhuǎn)電滲流速度的影響,其中 β=0,K=10,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1Fig.3.When there is a no-slip boundary condition,the influence of power-law fluid behavior index n on rotating electroosmotic flow velocity with the external magnetic field,β=0,K=10,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1 .

    圖5 中給出了哈特曼數(shù)Ha對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響.圖5 和圖6 展示了當(dāng)冪律流體行為指數(shù)n<1 及n>1 時,哈特曼數(shù)Ha對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響.由圖5(a)和圖5(b)與圖6(a)和圖6(b)可以看出,當(dāng)哈特曼數(shù)Ha>1 時,速度隨著哈特曼數(shù)Ha的增加而減小,因為當(dāng)哈特曼數(shù)Ha大時,洛倫茲力里的阻力部分σeuB2變大,總的洛倫茲力變成遠(yuǎn)遠(yuǎn)小于電滲力.但當(dāng)哈特曼數(shù)Ha<1 時,x方向速度u的大小隨著Ha的增加而增加,如圖5(c)和圖6(c)所示,這意味著在哈特曼數(shù)Ha小的情況下,洛倫茲力可以加快流體流動.

    圖6 哈特曼數(shù)Ha 對外加磁場的旋轉(zhuǎn)電滲流速度的影響,其中n=1.2,K=10,β=0.1,Ω=100 rad/s,=5 V,a=0.05,S=1Fig.6.The influence of Hartmann number Ha on the velocity of rotating electroosmotic flow with external magnetic field,n=1.2,K=10,β=0.1,Ω=100 rad/s,=5 V,a=0.05,S=1 .

    如圖7 和圖8 所示,圖中分別呈現(xiàn)了行為指數(shù)n=0.8 和n=1.2 時,旋轉(zhuǎn)角速度 Ω 對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響.從圖7(a)、圖7(c)、圖7(d)和圖8(a)、圖8(c)、圖8(d)中可以看出,當(dāng)沒有旋轉(zhuǎn)角速度時,x方向的u速度在內(nèi)部有一個平坦的速度平臺.隨著Ω的增大,外加磁場的旋轉(zhuǎn)電滲流速度u減小,與Xie 和 Jian[22]的結(jié)果相似.中心速度甚至在大的旋轉(zhuǎn)角速度時變?yōu)樨?fù)值(Ω=400,800).由圖7(b)和圖8(b) 表明,隨著旋轉(zhuǎn)角速度Ω的進(jìn)一步增大,y方向的速度v的中心速度表現(xiàn)出減小的特性.其結(jié)果與Chang 和Wang[9]相似.

    圖7 旋轉(zhuǎn)角速度 Ω 對外加磁場的旋轉(zhuǎn)電滲流速度的影響,其中 n=0.8,=5 V,a=0.05,Ha=1,S=1 (a) K=10,β=0.1;(b) K=10,β=0.1;(c) K=10,β=0;(d) K=20,β=0.1.Fig.7.The influence of the rotational angular velocity Ω on the rotational electroosmotic flow velocity of the external magnetic field,n=0.8,=5 V, a=0.05,Ha=1,S=1 (a) K=10,β=0.1;(b) K=10,β=0.1;(c) K=10,β=0;(d) K=20,β=0.1 .

    圖8 旋轉(zhuǎn)角速度 Ω 對外加磁場的旋轉(zhuǎn)電滲流速度的影響,其中 n=1.2,=5 V,a=0.05,Ha=1,S=1 (a) β=0.1,K=10. (b) β=0.1,K=10. (c) β=0,K=10. (d) β=0.1,K=30Fig.8.The influence of the rotational angular velocity Ω on the rotational electroosmotic flow velocity of the external magnetic field,n=1.2,=5 V, a=0.05,Ha=1,S=1 (a) β=0.1,K=10. (b) β=0.1,K=10. (c) β=0,K=10. (d) β=0.1,K=30 .

    圖9、圖10 和圖11 分別呈現(xiàn)了當(dāng)行為指數(shù)n=0.8 和n=1.2 時電動寬度K對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響.從圖9 和圖10 可以看出電動寬度K的增大導(dǎo)致外加磁場的旋轉(zhuǎn)電滲流速度的增加,并隨著n的減小電動寬度K影響的效果更加顯著.從圖10 和圖11 中可以看到,在滑移邊界條件下 (見圖10) 出現(xiàn)與不存在滑移邊界條件下 (Xie 和Jian[22],見圖11) 不同的結(jié)論.此外,從圖9(b)中,當(dāng)n=0.8 時,速度v呈碗狀,其大小隨電動寬度K增加.圖10(b)和圖11(b) 中,當(dāng)n=1.2 時,速度v呈拋物線形狀,在滑移邊界條件下速度隨電動寬度K增大而增大;無滑移邊界條件下隨電動寬度K增大導(dǎo)致速度減小.

    圖9 電動寬度 K 對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響,其中 n=0.8,β=0.1,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1Fig.9.The influence of the electric width K on the velocity distribution of rotating electroosmotic flow with external magnetic field,n=0.8,β=0.1,Ω=100 rad/s,=5 V, a=0.05,Ha=1,S=1 .

    圖10 電動寬度 K 對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響,其中 n=1.2,β=0.1,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1Fig.10.The influence of the electric width K on the velocity distribution of rotating electroosmotic flow with external magnetic field,n=1.2,β=0.1,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1 .

    圖11 電動寬度 K 對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響,其中 n=1.2,β=0,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1Fig.11.The influence of the electric width K on the velocity distribution of rotating electroosmotic flow with external magnetic field,n=1.2,β=0,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1 .

    當(dāng)冪律流體行為指數(shù)分別為n=0.8 和n=1.2 時,不同的滑移參數(shù)β對外加磁場的旋轉(zhuǎn)電滲流速度分布的影響如圖12 和圖13 所示.從圖12和圖13 可以發(fā)現(xiàn),隨著滑移參數(shù)β的增加,x方向和y方向的速度均相應(yīng)的增大.由于壁面的滑移效應(yīng),變截面微管道的中心速度要較為平穩(wěn),接近壁面的速度受到β值的影響出現(xiàn)劇烈變化,并快速減小.從圖中可以發(fā)現(xiàn)由于滑移邊界條件的存在,壁面上的速度分布出現(xiàn)差異,β值越小邊界附近的速度下降的越劇烈.當(dāng)滑移邊界條件不存在時,即β=0,邊界處的速度將下降到0.

    圖12 滑移參數(shù) β 對外加磁 場的旋轉(zhuǎn)電滲 流速度的影響,其 中 n=0.8,K=10,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1Fig.12.The influence of the slip parameter β on the rotating electroosmotic flow velocity with an external magnetic field,n=0.8,K=10,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1 .

    圖13 滑移參數(shù) β 對外加磁 場的旋轉(zhuǎn)電滲 流速度的影響,其 中 n=1.2,K=10,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1Fig.13.The influence of the slip parameter β on the rotating electroosmotic flow velocity with an external magnetic field,n=1.2,K=10,Ω=100 rad/s,=5 V,a=0.05,Ha=1,S=1 .

    5 結(jié)論

    本文建立了高zeta 勢下具有滑移邊界條件的冪律流體在變截面微管道中的外加磁場的旋轉(zhuǎn)電滲流的數(shù)學(xué)模型,利用有限差分法,對冪律流體進(jìn)行了數(shù)值研究,得到了速度分量的數(shù)值解.結(jié)果表明,在相同的條件下,變截面微管道中假塑性流體(n<1) 的速度比膨脹塑性流體 (n>1) 的速度高.當(dāng)哈特曼數(shù)Ha>1 時,速度隨著哈特曼數(shù)Ha的增加而減小;但當(dāng)哈特曼數(shù)Ha<1 時,x方向速度u的大小隨著Ha的增加而增加.旋轉(zhuǎn)角速度 Ω 越大,外加磁場旋轉(zhuǎn)電滲流的速度逐漸減小,影響效果由壁面逐漸向中心區(qū)域增加,且達(dá)到最大值的距離逐漸靠近壁面.旋轉(zhuǎn)系統(tǒng)產(chǎn)生的科里奧利力對外加磁場的電滲流起主導(dǎo)作用.電動寬度K的增大導(dǎo)致外加磁場的旋轉(zhuǎn)電滲流速度的增加,并隨著n的減小電動寬度K影響的效果更加顯著.滑移參數(shù)的增大使變截面微管道內(nèi)的流動增強(qiáng),最大值速度接近微管道壁面.

    附錄A

    兩邊同時乘以 dψ/dz得到:

    對上述表達(dá)式積分并應(yīng)用邊界條件(8)式和公式cosh 2x=2sinh2x+1,可得到:

    其中α=z0e/kBT,κ2=2n0z02e2/εkBT.

    由于 sinhx是一個奇函數(shù),并考慮邊界條件(7)式和(8)式,電勢ψ由以下方程給出:

    利用不定積分表達(dá)式:

    并對 (A4)式的等價形式進(jìn)行變換,d(αψ/2)/sinh(αψ/2)=κdz,在區(qū)間 [z,h] 上積分,可以得到電勢ψ分布的表達(dá)式,即(9)式:

    猜你喜歡
    磁場
    西安的“磁場”
    帶磁場的廣義Zakharov模型的奇異解
    為什么地球有磁場呢
    文脈清江浦 非遺“磁場圈”
    華人時刊(2020年13期)2020-09-25 08:21:42
    《磁場》易錯易混知識剖析
    帶電粒子在交變電、磁場中的運(yùn)動
    磁場的性質(zhì)和描述檢測題
    磁場測量儀的設(shè)計與實現(xiàn)
    電子制作(2016年1期)2016-11-07 08:42:44
    2016年春季性感磁場
    Coco薇(2016年1期)2016-01-11 16:53:24
    磁場對于電子分布的影響
    精品国产超薄肉色丝袜足j| 久99久视频精品免费| bbb黄色大片| 久久中文字幕人妻熟女| 国产91精品成人一区二区三区| 久久久久久免费高清国产稀缺| 欧美大码av| 国产日本99.免费观看| 一级a爱片免费观看的视频| 亚洲av电影不卡..在线观看| 毛片女人毛片| 天天添夜夜摸| 国产一区二区在线观看日韩 | 国产成+人综合+亚洲专区| 欧美色视频一区免费| 免费在线观看完整版高清| 97人妻精品一区二区三区麻豆| 中文亚洲av片在线观看爽| 日本成人三级电影网站| 日本一本二区三区精品| 熟女少妇亚洲综合色aaa.| 久久这里只有精品19| 成年人黄色毛片网站| 久久久久久大精品| 一区福利在线观看| 成人国产一区最新在线观看| 一级黄色大片毛片| 精品国产乱子伦一区二区三区| 午夜老司机福利片| 日本a在线网址| 久久精品综合一区二区三区| 亚洲精品美女久久久久99蜜臀| 国产精品久久久人人做人人爽| 精品电影一区二区在线| 日本黄大片高清| 久热爱精品视频在线9| 国产亚洲精品综合一区在线观看 | 久久久久国产精品人妻aⅴ院| 黄色a级毛片大全视频| 校园春色视频在线观看| 亚洲av电影在线进入| 天天躁夜夜躁狠狠躁躁| 激情在线观看视频在线高清| 亚洲成av人片在线播放无| 国产人伦9x9x在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人视频免费观看高清| 亚洲精品粉嫩美女一区| 又粗又爽又猛毛片免费看| 国产精品亚洲一级av第二区| 日韩欧美在线乱码| 欧洲精品卡2卡3卡4卡5卡区| 亚洲色图 男人天堂 中文字幕| 51午夜福利影视在线观看| 亚洲成人国产一区在线观看| av天堂在线播放| 成年版毛片免费区| 色综合婷婷激情| svipshipincom国产片| 在线观看免费午夜福利视频| 欧美成人性av电影在线观看| 哪里可以看免费的av片| www.自偷自拍.com| 日本熟妇午夜| 人妻丰满熟妇av一区二区三区| 操出白浆在线播放| 大型av网站在线播放| 亚洲av成人av| 男女午夜视频在线观看| 老汉色∧v一级毛片| 国产精品乱码一区二三区的特点| 首页视频小说图片口味搜索| 国产真实乱freesex| 午夜亚洲福利在线播放| 久久中文看片网| 男女下面进入的视频免费午夜| 亚洲精品中文字幕在线视频| 搡老熟女国产l中国老女人| 成人18禁在线播放| 97超级碰碰碰精品色视频在线观看| 淫秽高清视频在线观看| 免费电影在线观看免费观看| 国产午夜福利久久久久久| 一本综合久久免费| 亚洲人成网站高清观看| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲av嫩草精品影院| 日本 av在线| 国内少妇人妻偷人精品xxx网站 | 国模一区二区三区四区视频 | 人妻丰满熟妇av一区二区三区| 午夜a级毛片| 欧美成人性av电影在线观看| 亚洲真实伦在线观看| 99在线视频只有这里精品首页| 少妇被粗大的猛进出69影院| 中文字幕人成人乱码亚洲影| 97人妻精品一区二区三区麻豆| 妹子高潮喷水视频| 老熟妇乱子伦视频在线观看| 免费一级毛片在线播放高清视频| 人妻久久中文字幕网| 久久香蕉激情| 俺也久久电影网| 母亲3免费完整高清在线观看| 久久中文字幕一级| 欧美黑人精品巨大| 全区人妻精品视频| 不卡av一区二区三区| 亚洲人成网站高清观看| 国产亚洲欧美98| 给我免费播放毛片高清在线观看| a级毛片a级免费在线| 九色成人免费人妻av| 亚洲国产看品久久| 人人妻,人人澡人人爽秒播| 欧美成人午夜精品| 日日摸夜夜添夜夜添小说| 丰满人妻熟妇乱又伦精品不卡| 三级毛片av免费| 亚洲色图av天堂| 国产单亲对白刺激| 老鸭窝网址在线观看| 19禁男女啪啪无遮挡网站| 午夜免费观看网址| 女人高潮潮喷娇喘18禁视频| 成人特级黄色片久久久久久久| 香蕉久久夜色| 国产激情欧美一区二区| 午夜久久久久精精品| 黄色片一级片一级黄色片| 一本一本综合久久| 1024视频免费在线观看| 性色av乱码一区二区三区2| 天天添夜夜摸| 淫秽高清视频在线观看| 国产一区二区激情短视频| 美女黄网站色视频| 亚洲精品久久成人aⅴ小说| 美女大奶头视频| 无遮挡黄片免费观看| 色在线成人网| 日本黄色视频三级网站网址| 亚洲专区国产一区二区| 国产精品乱码一区二三区的特点| 国产亚洲精品一区二区www| 国产亚洲精品综合一区在线观看 | 欧美日韩瑟瑟在线播放| 熟妇人妻久久中文字幕3abv| 女人爽到高潮嗷嗷叫在线视频| 制服人妻中文乱码| 午夜精品一区二区三区免费看| 成人永久免费在线观看视频| 成人特级黄色片久久久久久久| 99国产极品粉嫩在线观看| 日韩欧美一区二区三区在线观看| 天堂√8在线中文| 一级片免费观看大全| 亚洲人成电影免费在线| 中出人妻视频一区二区| 免费观看精品视频网站| 看片在线看免费视频| 精品一区二区三区av网在线观看| 国产人伦9x9x在线观看| 黑人欧美特级aaaaaa片| 淫妇啪啪啪对白视频| 色综合欧美亚洲国产小说| 亚洲天堂国产精品一区在线| 看免费av毛片| 欧美精品亚洲一区二区| 欧美乱码精品一区二区三区| 精品久久久久久成人av| 午夜精品在线福利| 大型av网站在线播放| 男女床上黄色一级片免费看| 欧美午夜高清在线| 夜夜看夜夜爽夜夜摸| 精品一区二区三区视频在线观看免费| 久久香蕉精品热| 可以在线观看毛片的网站| 色尼玛亚洲综合影院| 狂野欧美激情性xxxx| 精品久久蜜臀av无| 人成视频在线观看免费观看| 欧美一区二区国产精品久久精品 | 国产乱人伦免费视频| 日日夜夜操网爽| 亚洲精品av麻豆狂野| 又黄又粗又硬又大视频| 91九色精品人成在线观看| 首页视频小说图片口味搜索| 国产熟女午夜一区二区三区| 99在线视频只有这里精品首页| 午夜免费激情av| 久久精品影院6| 国产又色又爽无遮挡免费看| 免费av毛片视频| 亚洲精品美女久久av网站| 亚洲av电影在线进入| 一边摸一边抽搐一进一小说| 香蕉国产在线看| 亚洲成a人片在线一区二区| 老司机福利观看| 久久精品国产99精品国产亚洲性色| 亚洲精品国产一区二区精华液| 日本三级黄在线观看| xxx96com| 色综合婷婷激情| 成人高潮视频无遮挡免费网站| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区视频了| 国产精品香港三级国产av潘金莲| 99精品欧美一区二区三区四区| 免费av毛片视频| 成人av一区二区三区在线看| 老汉色∧v一级毛片| 亚洲av成人一区二区三| 亚洲av电影不卡..在线观看| www.精华液| 成人一区二区视频在线观看| 久久久久免费精品人妻一区二区| 国产精品乱码一区二三区的特点| 欧美不卡视频在线免费观看 | 中文字幕人成人乱码亚洲影| 无遮挡黄片免费观看| 一a级毛片在线观看| 亚洲av成人不卡在线观看播放网| 一级a爱片免费观看的视频| 美女午夜性视频免费| 精品乱码久久久久久99久播| 亚洲美女视频黄频| 久久久久久久久久黄片| 国内毛片毛片毛片毛片毛片| 嫁个100分男人电影在线观看| av片东京热男人的天堂| 亚洲第一电影网av| www国产在线视频色| 国内精品久久久久精免费| 欧美绝顶高潮抽搐喷水| 欧美精品啪啪一区二区三区| 12—13女人毛片做爰片一| 亚洲成人精品中文字幕电影| 波多野结衣高清作品| 黄色视频,在线免费观看| 亚洲精华国产精华精| 99久久国产精品久久久| 亚洲自偷自拍图片 自拍| 国产精品久久久久久亚洲av鲁大| 久久精品国产综合久久久| 男人舔奶头视频| 亚洲专区中文字幕在线| 人妻丰满熟妇av一区二区三区| 又大又爽又粗| 在线观看日韩欧美| 国产成人av教育| 午夜精品久久久久久毛片777| 国产精品美女特级片免费视频播放器 | 制服丝袜大香蕉在线| 国内毛片毛片毛片毛片毛片| 精品免费久久久久久久清纯| 两个人视频免费观看高清| 搡老妇女老女人老熟妇| 大型av网站在线播放| 亚洲国产精品成人综合色| 国产成人精品久久二区二区91| 大型av网站在线播放| tocl精华| 99热这里只有精品一区 | 欧美zozozo另类| 禁无遮挡网站| 非洲黑人性xxxx精品又粗又长| 亚洲熟妇熟女久久| 国产亚洲精品久久久久5区| 夜夜爽天天搞| 国产成人精品无人区| 毛片女人毛片| 成年女人毛片免费观看观看9| 成年女人毛片免费观看观看9| 国产激情久久老熟女| 黄色a级毛片大全视频| 欧美久久黑人一区二区| 国产激情偷乱视频一区二区| 黄色毛片三级朝国网站| 岛国在线观看网站| 国产99白浆流出| 天天躁狠狠躁夜夜躁狠狠躁| 国产不卡一卡二| 国产av一区二区精品久久| 久久久水蜜桃国产精品网| 99国产综合亚洲精品| 高清毛片免费观看视频网站| av福利片在线| 很黄的视频免费| 久久久国产成人精品二区| 国产区一区二久久| 午夜两性在线视频| 日韩免费av在线播放| 好男人电影高清在线观看| 国产成人系列免费观看| 琪琪午夜伦伦电影理论片6080| 亚洲精品中文字幕一二三四区| 1024手机看黄色片| 很黄的视频免费| www.熟女人妻精品国产| 成人精品一区二区免费| 精品福利观看| 午夜视频精品福利| 国产野战对白在线观看| 久久婷婷人人爽人人干人人爱| 亚洲精品色激情综合| 欧美性猛交╳xxx乱大交人| 亚洲国产精品成人综合色| 欧美人与性动交α欧美精品济南到| 国产成人欧美在线观看| 成熟少妇高潮喷水视频| 丰满人妻熟妇乱又伦精品不卡| 可以在线观看毛片的网站| 亚洲国产精品合色在线| 动漫黄色视频在线观看| 少妇人妻一区二区三区视频| 国产成年人精品一区二区| 中文字幕熟女人妻在线| 在线观看午夜福利视频| 免费在线观看成人毛片| 高清在线国产一区| 人人妻人人看人人澡| 又黄又爽又免费观看的视频| 精品国产亚洲在线| 男人舔女人下体高潮全视频| 久久久精品欧美日韩精品| 高清毛片免费观看视频网站| 免费看日本二区| netflix在线观看网站| 在线a可以看的网站| 国产一区在线观看成人免费| 亚洲熟妇熟女久久| 国产精品久久视频播放| 亚洲国产日韩欧美精品在线观看 | 国产蜜桃级精品一区二区三区| 舔av片在线| 成年人黄色毛片网站| 老司机福利观看| 久久精品亚洲精品国产色婷小说| 好男人电影高清在线观看| 久久精品国产亚洲av香蕉五月| 琪琪午夜伦伦电影理论片6080| 亚洲成av人片免费观看| 老司机福利观看| 国产精品 欧美亚洲| 女人爽到高潮嗷嗷叫在线视频| 99热这里只有是精品50| av福利片在线观看| 两个人免费观看高清视频| 日韩av在线大香蕉| 精品日产1卡2卡| 激情在线观看视频在线高清| 18禁观看日本| 国产高清激情床上av| 啪啪无遮挡十八禁网站| svipshipincom国产片| 国产午夜福利久久久久久| 变态另类丝袜制服| 18禁黄网站禁片午夜丰满| 一本综合久久免费| 日韩成人在线观看一区二区三区| 国产成人精品无人区| 久久国产精品人妻蜜桃| ponron亚洲| 日日干狠狠操夜夜爽| xxx96com| 日韩免费av在线播放| 一个人免费在线观看的高清视频| 亚洲欧美激情综合另类| 757午夜福利合集在线观看| 久久天堂一区二区三区四区| 可以在线观看的亚洲视频| 亚洲国产日韩欧美精品在线观看 | √禁漫天堂资源中文www| 午夜两性在线视频| 欧美一区二区国产精品久久精品 | 韩国av一区二区三区四区| 97人妻精品一区二区三区麻豆| 亚洲精华国产精华精| 50天的宝宝边吃奶边哭怎么回事| www.自偷自拍.com| 久久久久久久精品吃奶| 亚洲在线自拍视频| 亚洲国产精品sss在线观看| 亚洲一码二码三码区别大吗| 亚洲天堂国产精品一区在线| 亚洲成人久久性| 免费观看人在逋| 成人18禁在线播放| 欧美日韩黄片免| 日韩欧美免费精品| 国产高清videossex| 亚洲全国av大片| 搡老妇女老女人老熟妇| www.熟女人妻精品国产| 久久久久久人人人人人| 精品久久久久久久毛片微露脸| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成网站在线播放欧美日韩| 99在线人妻在线中文字幕| 天天添夜夜摸| 在线看三级毛片| 婷婷精品国产亚洲av| 亚洲男人天堂网一区| 免费在线观看成人毛片| 又黄又爽又免费观看的视频| 亚洲精品色激情综合| 首页视频小说图片口味搜索| 亚洲熟女毛片儿| 亚洲狠狠婷婷综合久久图片| www.熟女人妻精品国产| 日韩免费av在线播放| 两个人视频免费观看高清| 精品久久久久久久久久免费视频| av天堂在线播放| a级毛片在线看网站| 91九色精品人成在线观看| 国产精品免费一区二区三区在线| 18禁黄网站禁片午夜丰满| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲欧美在线一区二区| www.熟女人妻精品国产| 精品人妻1区二区| 人人妻,人人澡人人爽秒播| 国产精品一及| 亚洲人与动物交配视频| 国产精品免费视频内射| 国产精品九九99| 99久久综合精品五月天人人| 人妻夜夜爽99麻豆av| 嫩草影院精品99| 两个人看的免费小视频| 午夜福利在线在线| 日韩高清综合在线| 91九色精品人成在线观看| 免费电影在线观看免费观看| 欧美色欧美亚洲另类二区| 一级黄色大片毛片| 欧美性猛交黑人性爽| 在线播放国产精品三级| 男女午夜视频在线观看| 18禁裸乳无遮挡免费网站照片| 中文资源天堂在线| 黄色女人牲交| 18禁黄网站禁片午夜丰满| 欧美精品亚洲一区二区| 国产黄a三级三级三级人| 校园春色视频在线观看| 精品一区二区三区视频在线观看免费| 一夜夜www| 国产激情偷乱视频一区二区| 无遮挡黄片免费观看| 黄频高清免费视频| 黄色a级毛片大全视频| 丁香欧美五月| 久9热在线精品视频| 日韩国内少妇激情av| 国产午夜精品久久久久久| 国产片内射在线| 久久国产乱子伦精品免费另类| 亚洲成av人片免费观看| 日韩国内少妇激情av| 一二三四社区在线视频社区8| 老司机福利观看| 欧美又色又爽又黄视频| 亚洲成av人片免费观看| 久久久久免费精品人妻一区二区| 国产精品九九99| 最新美女视频免费是黄的| 美女免费视频网站| 亚洲成人久久爱视频| 国产一级毛片七仙女欲春2| 成人午夜高清在线视频| 亚洲av成人一区二区三| 丰满人妻一区二区三区视频av | 亚洲专区国产一区二区| 国内毛片毛片毛片毛片毛片| 国产探花在线观看一区二区| 婷婷六月久久综合丁香| 国产97色在线日韩免费| 人妻夜夜爽99麻豆av| 99riav亚洲国产免费| 国内精品一区二区在线观看| 成在线人永久免费视频| 看免费av毛片| 在线永久观看黄色视频| 国产亚洲欧美98| av天堂在线播放| 好男人在线观看高清免费视频| 久热爱精品视频在线9| 亚洲狠狠婷婷综合久久图片| 中文在线观看免费www的网站 | 狠狠狠狠99中文字幕| 三级毛片av免费| 久久久久久久久免费视频了| 欧美日韩精品网址| 人人妻人人澡欧美一区二区| 十八禁网站免费在线| 日日干狠狠操夜夜爽| 午夜福利视频1000在线观看| 久久久久久亚洲精品国产蜜桃av| 国产视频一区二区在线看| 免费搜索国产男女视频| 黑人操中国人逼视频| 天堂动漫精品| 高清在线国产一区| 男女下面进入的视频免费午夜| 免费看日本二区| 亚洲美女黄片视频| 国产69精品久久久久777片 | 大型黄色视频在线免费观看| 欧美色视频一区免费| 久久久久久免费高清国产稀缺| 久久精品91无色码中文字幕| 国产伦一二天堂av在线观看| 亚洲人成77777在线视频| 黄色a级毛片大全视频| 他把我摸到了高潮在线观看| √禁漫天堂资源中文www| 欧美黑人精品巨大| 色综合欧美亚洲国产小说| 最好的美女福利视频网| 中文字幕精品亚洲无线码一区| 久久久久国产精品人妻aⅴ院| 极品教师在线免费播放| 天堂av国产一区二区熟女人妻 | 亚洲人成电影免费在线| 黄色 视频免费看| 国产av又大| 国产视频一区二区在线看| 免费人成视频x8x8入口观看| 超碰成人久久| av在线天堂中文字幕| 在线观看舔阴道视频| av在线天堂中文字幕| 女人爽到高潮嗷嗷叫在线视频| 香蕉丝袜av| 亚洲国产日韩欧美精品在线观看 | 亚洲精品中文字幕一二三四区| 中文字幕av在线有码专区| 久久久水蜜桃国产精品网| 日韩大尺度精品在线看网址| 亚洲av美国av| 国产人伦9x9x在线观看| 国产欧美日韩一区二区三| 久久热在线av| 丰满的人妻完整版| 国产人伦9x9x在线观看| 深夜精品福利| 久久久精品大字幕| 亚洲av美国av| 久久热在线av| 长腿黑丝高跟| 看片在线看免费视频| 久久香蕉国产精品| 亚洲av中文字字幕乱码综合| 国产精品日韩av在线免费观看| 国产1区2区3区精品| 变态另类丝袜制服| 国产1区2区3区精品| 亚洲av电影在线进入| 亚洲国产日韩欧美精品在线观看 | 亚洲欧美精品综合一区二区三区| 国产成人精品久久二区二区91| 亚洲午夜精品一区,二区,三区| 色播亚洲综合网| 又紧又爽又黄一区二区| 91国产中文字幕| 妹子高潮喷水视频| 国产精品美女特级片免费视频播放器 | 好男人在线观看高清免费视频| 国产高清视频在线播放一区| 欧美日韩福利视频一区二区| 日韩欧美 国产精品| 免费在线观看黄色视频的| 男人舔奶头视频| av片东京热男人的天堂| 美女大奶头视频| 国产激情久久老熟女| 丰满的人妻完整版| 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 免费在线观看日本一区| 19禁男女啪啪无遮挡网站| 精品久久久久久久人妻蜜臀av| 中文字幕熟女人妻在线| 少妇熟女aⅴ在线视频| cao死你这个sao货| 一二三四在线观看免费中文在| 这个男人来自地球电影免费观看| 三级男女做爰猛烈吃奶摸视频| 久久久久九九精品影院| 999久久久精品免费观看国产| 久99久视频精品免费| 波多野结衣高清无吗| 中出人妻视频一区二区| av在线播放免费不卡| 国产人伦9x9x在线观看| 美女高潮喷水抽搐中文字幕| 91老司机精品| 久久精品综合一区二区三区| 亚洲av片天天在线观看| 欧美日韩乱码在线| 在线十欧美十亚洲十日本专区| 亚洲黑人精品在线| 国产精品免费视频内射| 色精品久久人妻99蜜桃| 亚洲中文字幕一区二区三区有码在线看 | 国产精品久久久久久精品电影| 日韩欧美国产在线观看| 一二三四社区在线视频社区8|