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

    基于格子Boltzmann 方法的冪律流體二維頂蓋驅(qū)動(dòng)流轉(zhuǎn)捩研究*

    2021-10-08 08:55:54張恒任峰胡海豹
    物理學(xué)報(bào) 2021年18期
    關(guān)鍵詞:牛頓流體雷諾數(shù)二階

    張恒 任峰 胡海豹

    (西北工業(yè)大學(xué)航海學(xué)院,西安 710072)

    研究非牛頓流體轉(zhuǎn)捩問(wèn)題,可為調(diào)控非牛頓流體動(dòng)力特性提供理論基礎(chǔ).相對(duì)于牛頓流體轉(zhuǎn)捩問(wèn)題,非牛頓流體轉(zhuǎn)捩研究較少,缺乏轉(zhuǎn)捩雷諾數(shù)精細(xì)預(yù)報(bào)方法.論文以格子Boltzmann 方法為核心求解器,以典型非牛頓流體冪律模型為例,開(kāi)展了冪律流體二維頂蓋驅(qū)動(dòng)流轉(zhuǎn)捩模擬,給出剪切變稀和剪切增稠流體的第一轉(zhuǎn)捩雷諾數(shù),并分析了轉(zhuǎn)捩雷諾數(shù)附近流場(chǎng)時(shí)頻域特性及模態(tài)分布.結(jié)果表明,剪切變稀流體和剪切增稠流體的第一轉(zhuǎn)捩雷諾數(shù)與牛頓流體差異顯著,且在轉(zhuǎn)捩臨界雷諾數(shù)附近監(jiān)控點(diǎn)處速度分量均呈現(xiàn)周期性變化趨勢(shì).通過(guò)對(duì)流場(chǎng)速度和渦量的本征正交分解發(fā)現(xiàn),不同類(lèi)型的流體在轉(zhuǎn)捩臨界雷諾數(shù)附近,前兩階模態(tài)均為流場(chǎng)的主模態(tài),能量占比超過(guò)95%,且同類(lèi)型流體不同雷諾數(shù)的主模態(tài)間具有相似的結(jié)構(gòu).

    1 引 言

    頂蓋驅(qū)動(dòng)流是計(jì)算流體力學(xué)中經(jīng)典物理模型之一,被廣泛應(yīng)用于層流流動(dòng)、流動(dòng)轉(zhuǎn)捩等問(wèn)題的研究[1?4].研究流動(dòng)轉(zhuǎn)捩問(wèn)題,可為界定流場(chǎng)狀態(tài)和調(diào)控流體動(dòng)力特性提供理論依據(jù),具有重要學(xué)術(shù)意義.轉(zhuǎn)捩問(wèn)題涉及多種流動(dòng)狀態(tài),屬?gòu)?fù)雜流體力學(xué)問(wèn)題,已有大量學(xué)者開(kāi)展過(guò)相關(guān)研究.例如,Peng 等[5]采用有限差分法研究了頂蓋驅(qū)動(dòng)流的轉(zhuǎn)捩雷諾數(shù),并給出臨界雷諾數(shù)附近流場(chǎng)的時(shí)頻域特性.Bruneau 和Saad[6]進(jìn)一步提高網(wǎng)格分辨率,改進(jìn)對(duì)流項(xiàng)差分格式,佐證了Peng 等[5]、Auteri 等[7]、Sahin 和Owens[8]關(guān)于第一轉(zhuǎn)捩雷諾數(shù)在8000 附近的結(jié)論.需要指出的是,上述研究均是以牛頓流體為研究對(duì)象.

    非牛頓流體的剪應(yīng)力與剪切應(yīng)變率之間為非線(xiàn)性關(guān)系,會(huì)表現(xiàn)出剪切增稠[9]、剪切變稀[10]等不同于牛頓流體的運(yùn)動(dòng)行為.非牛頓流體廣泛存在于自然界和工業(yè)生產(chǎn)中,如化工中的聚乙烯、聚氯乙烯,食品中的淀粉液,人體內(nèi)的血液等.非牛頓流體種類(lèi)眾多,其中具有代表性的是冪律流體.冪律流體按照冪律指數(shù)主要分為剪切變稀流體、剪切增稠流體.作為特例,牛頓流體可視為指數(shù)為1的冪律流體.學(xué)者們開(kāi)展了大量的冪律流體實(shí)驗(yàn)與仿真的研究,如Johnston 等[11]發(fā)現(xiàn)相比于牛頓流體,用冪律流體模擬血液流動(dòng),可獲更精準(zhǔn)的壁面切應(yīng)力;Zhao 等[12]模擬了冪律流體在微通道中的電滲流動(dòng);Hojjat 等[13]通過(guò)實(shí)驗(yàn)研究了含納米顆粒流體的剪切變稀行為.但是,關(guān)于冪律流體轉(zhuǎn)捩問(wèn)題的鮮有研究報(bào)道.

    與有限差分和有限體積等傳統(tǒng)方法不同,格子Boltzmann 方 法(lattice Boltzmann method,LBM)[14]通過(guò)對(duì)分布函數(shù)進(jìn)行碰撞和遷移處理,可實(shí)現(xiàn)宏觀(guān)流場(chǎng)時(shí)空演化過(guò)程的模擬,具有并行性高,邊界處理簡(jiǎn)單等優(yōu)勢(shì),已在牛頓流體仿真領(lǐng)域取得大量研究成果[15?17].近年來(lái),也有學(xué)者將LBM拓展至非牛頓流體模擬,如Boyd 等[18]在LBM 框架下實(shí)現(xiàn)了冪律流體中速度導(dǎo)數(shù)的快速求解;Wang和Ho[19]在LBM的框架下模擬了剪切變稀的血液流動(dòng);Chen 等[20]模擬了冪律流體通過(guò)多孔介質(zhì)的流動(dòng).然而,將LBM 應(yīng)用于冪律流體轉(zhuǎn)捩問(wèn)題研究仍有待探索.

    基于此,本文以二維頂蓋驅(qū)動(dòng)流為物理模型,以格子Boltzmann 方法為核心求解器,采用非牛頓流體冪律模型,開(kāi)展了非牛頓流體二維頂蓋驅(qū)動(dòng)流轉(zhuǎn)捩模擬,并借助快速傅里葉變換、本征正交分解(proper orthogonal decomposition,POD)等方法,分析了轉(zhuǎn)捩點(diǎn)附近的流場(chǎng)時(shí)頻域特性及模態(tài)分布.

    2 數(shù)值方法

    2.1 物理模型

    頂蓋驅(qū)動(dòng)流模型如圖1 所示,上頂板為運(yùn)動(dòng)壁面,速度為U0,其余3 個(gè)壁面均為靜止無(wú)滑移壁面.模型水平和豎直方向上長(zhǎng)度分別為L(zhǎng)x和Ly,且Lx=Ly=L0,L0為特征長(zhǎng)度.圖1 中,坐標(biāo)原點(diǎn)O位于底面與左壁面交點(diǎn)處.為定量表征流場(chǎng)狀態(tài)與雷諾數(shù)(Reynolds number,Re)之間的關(guān)系,選取監(jiān)控點(diǎn)坐標(biāo)為(0.6L0,0.6L0).

    圖1 算例模型示意圖Fig.1.Schematic diagram of calculation model.

    冪律流體運(yùn)動(dòng)黏性系數(shù)的計(jì)算公式如下:

    式中,S為流場(chǎng)中格點(diǎn)處應(yīng)變率張量,γ˙ 為剪切率,n為冪律指數(shù),ν為運(yùn)動(dòng)黏性系數(shù),m為冪律系數(shù)(與Re有關(guān)的常量).當(dāng)n=1 時(shí)為牛頓流體;n<1為剪切變稀流體;n>1 時(shí)為剪切增稠流體.

    雷諾數(shù)的定義為[21]

    2.2 數(shù)值方法

    采用LBM 多松弛時(shí)間模型[22],并使用基于統(tǒng)一計(jì)算設(shè)備構(gòu)架的顯卡并行加速技術(shù)[23]來(lái)提高計(jì)算速度.其中,多松弛時(shí)間模型表達(dá)式為

    式中,f(x,t) 為t時(shí)刻流場(chǎng)中坐標(biāo)位置為x處的分布函數(shù),ei為離散速度,δt為時(shí)間步長(zhǎng),M為從速度空間向矩空間進(jìn)行轉(zhuǎn)換的矩陣,為平衡態(tài)分布函數(shù).

    針對(duì)文中研究的二維頂蓋驅(qū)動(dòng)流,選用經(jīng)典D2Q9 模型[24],離散速度ei為

    LBM 平衡態(tài)分布函數(shù)采用不可壓模型[25],其表達(dá)式為

    式中,ρ0為不可壓流體的密度;p為壓強(qiáng);u為格點(diǎn)處的速度;權(quán)系數(shù)wi的取值為w0=4/9 ,w1?4=1/9,w5?8=1/36 ;cs為格子聲速.松弛參數(shù)對(duì)角矩陣選為Λf=[1.0,1.1,1.0,1.0,1.2,1.0,1.2,1/τ,1/τ],松弛時(shí)間τ與運(yùn)動(dòng)黏性關(guān)系為ν=(τ ?0.5).

    流場(chǎng)中宏觀(guān)速度與壓強(qiáng)的求解公式為

    為保證流場(chǎng)時(shí)空分辨率,本文最終選用的網(wǎng)格分辨率為512 × 512,計(jì)算時(shí)間總長(zhǎng)2000T0(T0為參考時(shí)間),即1600 萬(wàn)個(gè)格子時(shí)間步長(zhǎng).靜壁邊界采用半步長(zhǎng)反彈格式[26],動(dòng)壁采用動(dòng)量補(bǔ)償格式[27].為防止運(yùn)動(dòng)黏性系數(shù)過(guò)大或者過(guò)小造成計(jì)算結(jié)果發(fā)散,這里設(shè)置運(yùn)動(dòng)黏性系數(shù)上限和下限[28]分別為 2 0ν0,ν0/20 .不失一般性,文中取n=0.75,1.25 分別代表剪切變稀流體和剪切增稠流體.具體模擬流程如圖2 所示.

    圖2 LBM 與冪律模型耦合計(jì)算流程示意圖Fig.2.Diagram of coupling calculation process between LBM and power law model.

    2.3 POD 方法

    POD[29,30]是一種廣泛應(yīng)用于流場(chǎng)降階模型中的方法,可用于流場(chǎng)特征信息的提取.以速度POD為例,已知N個(gè)時(shí)刻流場(chǎng)速度信息u(x,ti),i=1—N,則有

    另 外,Σ=diag(σ1,σ2,···σi,···)中各σi模長(zhǎng)的平方可用來(lái)表征對(duì)應(yīng)模態(tài)的“能量”.按照模態(tài)能量的大小進(jìn)行排序即可得到流場(chǎng)的主要模態(tài).

    2.4 第一轉(zhuǎn)捩雷諾數(shù)計(jì)算

    在頂蓋驅(qū)動(dòng)流中,隨著Re增加,呈現(xiàn)出的流動(dòng)狀態(tài)往往是從層流逐步過(guò)渡到湍流,并存在Hopf 分岔現(xiàn)象.通過(guò)在流場(chǎng)的特定位置處設(shè)置監(jiān)控點(diǎn)的方法[5,31],可以研究Hopf 分岔現(xiàn)象.開(kāi)始產(chǎn)生分岔現(xiàn)象時(shí)的雷諾數(shù),即第一轉(zhuǎn)捩雷諾數(shù)可以通過(guò)如下的方式計(jì)算:

    其中 max[u(t)] 與 min[u(t)] 分別為流場(chǎng)充分發(fā)展后,監(jiān)控點(diǎn)處水平速度u的最大值與最小值.鑒于與Re成正比關(guān)系,因此可以通過(guò)計(jì)算不同雷諾數(shù)下的,使用最小二乘法線(xiàn)性擬合出Re-圖.擬合所得的直線(xiàn)與=0的交點(diǎn),即流場(chǎng)從穩(wěn)定向周期態(tài)轉(zhuǎn)變的第一轉(zhuǎn)捩臨界雷諾數(shù).

    3 計(jì)算結(jié)果與分析

    3.1 算例驗(yàn)證

    為了驗(yàn)證計(jì)算的準(zhǔn)確性,這里首先給出了牛頓流體頂蓋驅(qū)動(dòng)流模擬結(jié)果.圖3 為Re=5000 時(shí),方腔中心線(xiàn)上的速度分布與Peng 等[5]結(jié)果的對(duì)比,其中X-v曲線(xiàn)顯示的是豎直方向上方腔的中心線(xiàn)上速度分量分布,即縱坐標(biāo)y/L=0.5 上的X和v之間的對(duì)應(yīng)關(guān)系;Y-u曲線(xiàn)顯示的是水平方向上的中心線(xiàn)上速度分量分布,即橫坐標(biāo)x/L=0.5上的Y和u之間的關(guān)系.Peng 等[5]在牛頓流體頂蓋驅(qū)動(dòng)流轉(zhuǎn)捩問(wèn)題研究中,使用的是限差分法,本文的計(jì)算結(jié)果與Peng 等[5]對(duì)比良好.值得注意的是Peng 等[5]求解對(duì)流項(xiàng)使用的是七階迎風(fēng)格式,擴(kuò)散項(xiàng)使用的是六階中心差分格式.

    圖3 頂蓋驅(qū)動(dòng)流在中心線(xiàn)處的速度分布Fig.3.Velocity distribution at center line of lid-driven flow.

    3.2 第一轉(zhuǎn)捩雷諾數(shù)

    按照文中2.4 節(jié)的方法得到牛頓流體、剪切變稀流體和剪切增稠流體的第一轉(zhuǎn)捩雷諾數(shù)依次為7835,5496 和11546,結(jié)果如圖4 所示.

    圖4 轉(zhuǎn)捩臨界雷諾數(shù)Fig.4.Transition critical Reynolds number.

    將本文計(jì)算的牛頓流體轉(zhuǎn)捩雷諾數(shù)與文獻(xiàn)的結(jié)果進(jìn)行比對(duì),如表1 所列,可以看出,本文計(jì)算的牛頓流體轉(zhuǎn)捩臨界雷諾數(shù)相對(duì)于Peng 等[5]和Poliashenko 與Aidu[32]的計(jì)算結(jié)果而言較大,而相對(duì)于Cazemier 等[31]和Bruneau 與Saad[6]的計(jì)算結(jié)果而言較小,但是從總的相對(duì)誤差的值來(lái)看,誤差上限不超過(guò)5%,可以看出本文計(jì)算可靠性較高.

    表1 牛頓流體轉(zhuǎn)捩臨界雷諾數(shù)計(jì)算結(jié)果的對(duì)比Table 1. Comparison of calculation results of critical Reynolds number for Newtonian fluid transition.

    相比于牛頓流體,剪切變稀流體的轉(zhuǎn)捩雷諾數(shù)較低,即剪切變稀流體更容易從定常狀態(tài)向周期態(tài)轉(zhuǎn)變.而由剪切變稀流體本身的定義可知,剪切速率越大,其黏性會(huì)越低,流動(dòng)狀態(tài)也越容易發(fā)生轉(zhuǎn)變,這與計(jì)算結(jié)果表現(xiàn)出的性質(zhì)是相符合的.與牛頓流體相比,剪切增稠流體的轉(zhuǎn)捩臨界雷諾數(shù)較高,即剪切增稠流體更不易從定常狀態(tài)向周期態(tài)轉(zhuǎn)變.而根據(jù)剪切增稠流體本身的定義可知,剪切速率越大,其黏性會(huì)越大,流動(dòng)狀態(tài)也越不易發(fā)生轉(zhuǎn)變.

    3.3 速度監(jiān)控點(diǎn)時(shí)頻域特性

    選取三類(lèi)流體在轉(zhuǎn)捩雷諾數(shù)附近的流場(chǎng),進(jìn)行時(shí)頻域分析.圖5 給出了三類(lèi)流體在轉(zhuǎn)捩雷諾數(shù)附近,速度監(jiān)控點(diǎn)處的時(shí)頻域特性圖.其中圖5(a)—(c)給出了三類(lèi)流體在轉(zhuǎn)捩前后,水平方向速度隨時(shí)間的變化,可以看出,三類(lèi)流體在轉(zhuǎn)捩發(fā)生前的定常流動(dòng)狀態(tài)下,速度為一水平直線(xiàn),在轉(zhuǎn)捩雷諾數(shù)附近,速度均呈現(xiàn)周期性變化的特點(diǎn),在轉(zhuǎn)捩后湍流狀態(tài)下,速度有明顯的大幅度波動(dòng).對(duì)水平速度作快速傅里葉變換得圖5(d)—(f),在發(fā)生轉(zhuǎn)捩前,由于其直流部分即f=0 為其主要分量,在頻譜圖上顯示為一點(diǎn),在轉(zhuǎn)捩雷諾數(shù)附近,水平速度均存在一個(gè)主頻f1與諧頻f2,且滿(mǎn)足f2=2f1,但不同類(lèi)型流體的f1值并不相同,在轉(zhuǎn)捩后湍流狀態(tài)下,出現(xiàn)了明顯的多頻現(xiàn)象.圖5(g)—(i)為三類(lèi)流體在轉(zhuǎn)捩雷諾數(shù)附近,水平速度u與豎直方向速度v之間形成的速度相圖.可以發(fā)現(xiàn)不同類(lèi)型的流體的速度相圖的變化范圍并不相同,但均為單一光滑封閉曲線(xiàn).

    圖5 三類(lèi)流體在速度監(jiān)控點(diǎn)處的時(shí)頻域特性 (a),(d) n=1,Re=7500,8500,16000;(g) n=1,Re=8500;(b),(e) n=0.75,Re=5500,6500,10000;(h) n=0.75,Re=6500;(c),(f) n=1.25,Re=9500,13000,20000;(i) n=1.25,Re=13000Fig.5.Time-frequency spectrum characteristics at the velocity monitoring point for three types of fluids:(a),(d) n=1,Re=7500,8500,16000;(g) n=1,Re=8500;(b),(e) n=0.75,Re=5500,6500,10000;(h) n=0.75,Re=6500;(c),(f) n=1.25,Re=9500,13000,20000;(i) n=1.25,Re=13000.

    圖6 為監(jiān)控點(diǎn)處牛頓流體、剪切變稀與剪切增稠流體在不同雷諾數(shù)情況下的速度相圖,各封閉曲線(xiàn)中心點(diǎn)用虛線(xiàn)進(jìn)行連接.圖6(a)為牛頓流體在不同雷諾數(shù)下的速度相圖,其中Re=7500,8000時(shí),速度波動(dòng)極小,在相圖上呈現(xiàn)的是1 個(gè)點(diǎn),而在Re=8500—9500 時(shí),相圖為光滑封閉曲線(xiàn),曲線(xiàn)整體向右上方移動(dòng)且包圍的面積在不斷增加.圖6(b)為剪切變稀流體在不同雷諾數(shù)下的速度相圖,在Re=5500 時(shí),速度相圖為1 個(gè)點(diǎn),而Re=5700—6000 時(shí),速度相圖為封閉曲線(xiàn),隨著Re增加向上方移動(dòng),且包圍的面積在不斷增加.圖6(c)為剪切增稠流體在不同雷諾數(shù)下的速度相圖,其中Re=11500,12000,12500 時(shí),速度相圖為單一點(diǎn),而Re=13000—15000 時(shí),速度相圖為光滑封閉曲線(xiàn).速度相圖整體變化趨勢(shì)為先向右再向左上方移動(dòng),向右移動(dòng)的過(guò)程中,豎直方向上速度v幾乎沒(méi)有變化,而向左上方移動(dòng)的過(guò)程中,封閉曲線(xiàn)包圍的面積逐漸增大,有較為明顯的變化規(guī)律.

    圖6 在監(jiān)控點(diǎn)處的速度相圖 (a) 牛頓流體;(b) 剪切變稀流體;(c) 剪切增稠流體Fig.6.Velocity phase diagrams at the monitoring point:(a) Newtonian fluid;(b) shear-thinning fluid;(c) shear-thickening fluid.

    3.4 流場(chǎng)模態(tài)分析

    進(jìn)一步使用POD 對(duì)牛頓流體、剪切變稀和剪切增稠流體的流場(chǎng)速度、渦量和黏性進(jìn)行處理,結(jié)果如圖7 和圖8 所示.圖7(a)—(c)分別為牛頓流體、剪切變稀流體和剪切增稠流體的水平速度POD能量占比圖.從圖7 可以發(fā)現(xiàn),不同類(lèi)型的流體在第一臨界轉(zhuǎn)捩雷諾數(shù)附近且發(fā)生轉(zhuǎn)捩后,均呈現(xiàn)前兩階模態(tài)能量占主導(dǎo)地位的特點(diǎn),能量占比超95%,且前兩階模態(tài)能量占比接近.在圖7(a)Re=10500 時(shí),出現(xiàn)了多模態(tài)的結(jié)果,這里的流場(chǎng)可能是進(jìn)入Hopf 第二分岔第二模式中[5].在圖7(c)Re=11500,12500 時(shí),一階模態(tài)能量占比大幅度高于二階模態(tài)的能量占比,但結(jié)合圖6(c)可以發(fā)現(xiàn),此時(shí)流場(chǎng)速度波動(dòng)量較小幾乎為定常流場(chǎng),一階模態(tài)與二階模態(tài)相對(duì)于平均流場(chǎng)幾乎可以忽略.

    圖7 速度u的各階模態(tài)的能量占比 (a) 牛頓流體;(b) 剪切變稀流體;(c) 剪切增稠流體Fig.7.Energy share of each order of mode for velocity u:(a) Newtonian fluid;(b) shear thinning-fluid;(c) shear-thickening fluid.

    圖8(a)—(c)分別為牛頓流體、剪切變稀流體和剪切增稠流體的渦量POD 能量占比圖.渦量POD 能量占比與水平速度POD 能量占比結(jié)果相比,同種流體在同一雷諾數(shù)下,渦量POD 一階模態(tài)的占比略高于水平速度POD 能量占比,渦量POD 中前兩階模仍占主導(dǎo)地位.

    圖8 渦量的各階模態(tài)的能量占比 (a) 牛頓流體;(b) 剪切變稀流體;(c) 剪切增稠流體Fig.8.Vortex energy share of each order of mode:(a) Newtonian fluid;(b) shear-thinning fluid;(c) shear-thickening fluid.

    圖9 為牛頓流體在不同雷諾數(shù)下各階模態(tài)的流場(chǎng)模態(tài)的速度場(chǎng)云圖,圖9(a)和圖9(d)分別為Re=8500 與Re=9000 時(shí)平均水平速度場(chǎng)云圖,圖9(b)和圖9(e) 分別為Re=8500 與Re=9000時(shí)水平速度一階模態(tài)云圖,圖9(c)和圖9(f) 分別為Re=8500 與Re=9000 時(shí)水平速度二階模態(tài)云圖.從速度平均場(chǎng)來(lái)看,Re=8500 與Re=9000的結(jié)果相近;從一階模態(tài)的云圖來(lái)看,峰值區(qū)域均靠近壁面,且峰值區(qū)域的形狀相似;從二階模態(tài)的云圖來(lái)看,二階模態(tài)的結(jié)果相近,且峰值主要分布在壁面處.

    圖9 牛頓流體的水平速度的各階模態(tài)圖 (a),(b),(c) Re=8500 時(shí),平均場(chǎng)、一階模態(tài)與二階模態(tài);(d),(e),(f) Re=9000 時(shí),平均場(chǎng)、一階模態(tài)與二階模態(tài)Fig.9.Modal diagrams of horizontal velocity of Newtonian fluid:(a),(b) and(c) The mean field,the first and second modes when Re=8500;(d),(e),(f) mean field,first-order mode and second-order mode when Re=9000.

    圖10 為剪切變稀流體在不同Re下各階模態(tài)的流場(chǎng)模態(tài)的速度場(chǎng)云圖,圖10(a)和圖10(d)分別為Re=5700 與Re=6000 時(shí)平均速度場(chǎng)云圖,圖10(b)和圖10(e)分別為Re=5700 與Re=6000時(shí)速度一階模態(tài)云圖,圖10(c)和圖10(f)分別為Re=5700 與Re=6000 時(shí)速度二階模態(tài)云圖.從平均場(chǎng)的結(jié)果來(lái)看,Re=5700 與Re=6000 結(jié)果相近;從一階和二階模態(tài)的結(jié)果來(lái)看,模態(tài)峰值分布區(qū)域集中在壁面附近,峰值分布區(qū)域的形狀相似,但峰值分布區(qū)的模態(tài)值相反.

    圖11 為剪切增稠流體在不同Re下各階模態(tài)的流場(chǎng)模態(tài)的速度場(chǎng)云圖,圖11(a)和圖11(d)分別為Re=13000 與Re=13500 時(shí)平均速度場(chǎng)云圖,圖11(b)和圖11(e)分別為Re=13000 與Re=13500 時(shí)速度一階模態(tài)云圖,圖11(c)和圖11(f)分別為Re=13000 與Re=13500 時(shí)速度二階模態(tài)云圖.從平均場(chǎng)的結(jié)果來(lái)看,Re=13000 與Re=13500 結(jié)果相近;從一階和二階模態(tài)的結(jié)果來(lái)看,峰值分布區(qū)域集中在壁面附近,且峰值分布區(qū)域的形狀相似,這與牛頓流體和剪切變稀流體的結(jié)果是類(lèi)似的.剪切增稠流體一階和二階模態(tài)峰值分布區(qū)的模態(tài)值相反.

    圖11 剪切增稠流體的水平速度的各階模態(tài)圖 (a),(b),(c) Re=13000 時(shí),平均場(chǎng)、一階模態(tài)與二階模態(tài);(d),(e),(f) Re=13500 時(shí),平均場(chǎng)、一階模態(tài)與二階模態(tài)Fig.11.Modal diagrams of horizontal velocity of shear-thickening fluid:(a),(b),(c) The mean field,the first and second modes when Re=13000;(d),(e),(f) mean field,first-order mode and second-order mode when Re=13500.

    4 結(jié) 論

    論文針對(duì)典型的非牛頓流體—冪律流體,利用格子Boltzmann的計(jì)算框架,以頂蓋驅(qū)動(dòng)流為物理模型,研究得到冪律指數(shù)在0.75 和1.25 時(shí)的第一臨界轉(zhuǎn)捩Re分別為5496 和11546.在轉(zhuǎn)捩Re附近且發(fā)生轉(zhuǎn)捩后,流場(chǎng)監(jiān)控點(diǎn)的速度隨時(shí)間呈周期性變化,存在一個(gè)主頻與諧頻.與轉(zhuǎn)捩前速度相圖為單一點(diǎn)不同,轉(zhuǎn)捩后速度相圖為光滑封閉曲線(xiàn),且隨著Re增加,曲線(xiàn)中心點(diǎn)朝一定的方向移動(dòng)且包絡(luò)面積不斷增加,但不同流體中心點(diǎn)移動(dòng)的方向不相同.流場(chǎng)的速度、渦量POD 分解的結(jié)果均表明,流場(chǎng)在轉(zhuǎn)捩臨界雷諾數(shù)附近,其前兩階模態(tài)為其主要模態(tài),能量占比超95%.從模態(tài)云圖的結(jié)果看,一階模態(tài)與二階模態(tài)速度峰值分布區(qū)域主要在壁面附近,且同種流體不同Re下的平均場(chǎng)結(jié)果相近,一階模態(tài)與二階模態(tài)的峰值分布區(qū)域形狀相似.

    猜你喜歡
    牛頓流體雷諾數(shù)二階
    一類(lèi)二階迭代泛函微分方程的周期解
    非牛頓流體
    什么是非牛頓流體
    少兒科技(2019年3期)2019-09-10 07:22:44
    區(qū)別牛頓流體和非牛頓流體
    一類(lèi)二階中立隨機(jī)偏微分方程的吸引集和擬不變集
    二階線(xiàn)性微分方程的解法
    一類(lèi)二階中立隨機(jī)偏微分方程的吸引集和擬不變集
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    首款XGEL非牛頓流體“高樂(lè)高”系列水溶肥問(wèn)世
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    一边亲一边摸免费视频| 亚洲精品亚洲一区二区| 国产成人精品一,二区| 国产v大片淫在线免费观看| 一区二区三区乱码不卡18| 在线观看免费高清a一片| 久久久精品欧美日韩精品| 成人特级av手机在线观看| 亚洲自偷自拍三级| www.色视频.com| 日本wwww免费看| av播播在线观看一区| 国产亚洲5aaaaa淫片| 男的添女的下面高潮视频| 国内少妇人妻偷人精品xxx网站| 美女脱内裤让男人舔精品视频| 三级经典国产精品| 伦理电影大哥的女人| 免费大片黄手机在线观看| 亚洲欧美成人精品一区二区| 欧美高清成人免费视频www| 日韩成人伦理影院| 久久久久久久国产电影| 晚上一个人看的免费电影| 国产v大片淫在线免费观看| 一级av片app| 国产精品久久久久久久久免| 国产精品一二三区在线看| 久久人人爽人人片av| 在线免费观看不下载黄p国产| 别揉我奶头 嗯啊视频| 欧美一区二区亚洲| 一边亲一边摸免费视频| 国产欧美另类精品又又久久亚洲欧美| 99热国产这里只有精品6| 亚洲自拍偷在线| 中文天堂在线官网| 成年版毛片免费区| 久久鲁丝午夜福利片| 在线精品无人区一区二区三 | 直男gayav资源| 久久99热6这里只有精品| 97精品久久久久久久久久精品| 国产精品一区www在线观看| 亚洲国产精品成人久久小说| 亚洲色图综合在线观看| 偷拍熟女少妇极品色| 天美传媒精品一区二区| 久久精品久久精品一区二区三区| 国产免费视频播放在线视频| 毛片一级片免费看久久久久| 国产成人免费无遮挡视频| av在线天堂中文字幕| 亚洲,一卡二卡三卡| av国产精品久久久久影院| 久久6这里有精品| 一区二区三区四区激情视频| 舔av片在线| 国产视频首页在线观看| 爱豆传媒免费全集在线观看| 日日撸夜夜添| 国产精品三级大全| 狂野欧美激情性xxxx在线观看| 综合色丁香网| 久久韩国三级中文字幕| 性色avwww在线观看| 久久久色成人| 99热全是精品| 国产精品国产三级国产专区5o| 人妻少妇偷人精品九色| 亚洲精品日本国产第一区| 久久久色成人| 青春草国产在线视频| 欧美高清性xxxxhd video| 成年av动漫网址| 国产亚洲av片在线观看秒播厂| 日韩欧美一区视频在线观看 | 亚洲精品,欧美精品| 精品人妻视频免费看| 亚洲欧洲国产日韩| 七月丁香在线播放| 国国产精品蜜臀av免费| 狂野欧美激情性xxxx在线观看| 免费大片18禁| 午夜福利视频精品| 亚洲精品日本国产第一区| 欧美丝袜亚洲另类| 女的被弄到高潮叫床怎么办| 乱码一卡2卡4卡精品| 久久这里有精品视频免费| 尾随美女入室| 99精国产麻豆久久婷婷| 久久热精品热| 国产真实伦视频高清在线观看| 亚洲成人中文字幕在线播放| 高清日韩中文字幕在线| 亚洲自偷自拍三级| av一本久久久久| 欧美性猛交╳xxx乱大交人| 免费观看性生交大片5| 国产亚洲精品久久久com| 禁无遮挡网站| 免费观看的影片在线观看| 久久久久网色| 欧美3d第一页| 中文在线观看免费www的网站| 联通29元200g的流量卡| 观看免费一级毛片| 舔av片在线| 亚洲国产最新在线播放| 在现免费观看毛片| 简卡轻食公司| 亚洲美女搞黄在线观看| 国产伦精品一区二区三区四那| 大香蕉久久网| 精品少妇黑人巨大在线播放| 亚洲成人av在线免费| 只有这里有精品99| www.色视频.com| 综合色av麻豆| 亚洲va在线va天堂va国产| 狂野欧美激情性bbbbbb| 成年av动漫网址| av专区在线播放| 一个人看视频在线观看www免费| 国产精品一区二区性色av| av播播在线观看一区| 又大又黄又爽视频免费| 亚洲国产精品专区欧美| 噜噜噜噜噜久久久久久91| 激情 狠狠 欧美| 亚洲精品久久午夜乱码| 国产在视频线精品| 日韩一区二区三区影片| 丝瓜视频免费看黄片| 国产一区有黄有色的免费视频| 尾随美女入室| 久久女婷五月综合色啪小说 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 我的老师免费观看完整版| 亚洲精品久久久久久婷婷小说| 久久久久精品久久久久真实原创| 精品久久国产蜜桃| 午夜免费鲁丝| 99久国产av精品国产电影| 蜜桃亚洲精品一区二区三区| 国产乱人偷精品视频| 18+在线观看网站| 久久精品国产亚洲av涩爱| 99精国产麻豆久久婷婷| 蜜臀久久99精品久久宅男| 九九久久精品国产亚洲av麻豆| 国国产精品蜜臀av免费| 卡戴珊不雅视频在线播放| 在现免费观看毛片| 国产一级毛片在线| 国产免费又黄又爽又色| 国产精品.久久久| 日本黄大片高清| 免费播放大片免费观看视频在线观看| eeuss影院久久| 亚洲欧洲日产国产| 男男h啪啪无遮挡| 亚洲人与动物交配视频| 内地一区二区视频在线| 欧美97在线视频| 九草在线视频观看| 国产精品一区二区性色av| 婷婷色综合大香蕉| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 97热精品久久久久久| 亚洲av不卡在线观看| 一级爰片在线观看| 国产又色又爽无遮挡免| h日本视频在线播放| 日韩av不卡免费在线播放| 日韩中字成人| 日韩av免费高清视频| 国产成人a∨麻豆精品| 国产中年淑女户外野战色| 精品99又大又爽又粗少妇毛片| 人妻 亚洲 视频| 一个人看的www免费观看视频| 男女啪啪激烈高潮av片| 国产精品av视频在线免费观看| 国产精品精品国产色婷婷| 精品一区二区三区视频在线| 亚洲不卡免费看| 日韩在线高清观看一区二区三区| 亚洲电影在线观看av| 91aial.com中文字幕在线观看| 国产一区亚洲一区在线观看| 免费大片18禁| 久久人人爽人人片av| 人妻一区二区av| 欧美bdsm另类| 国产精品久久久久久av不卡| 特大巨黑吊av在线直播| 三级经典国产精品| 黑人高潮一二区| 在线观看免费高清a一片| 精品亚洲乱码少妇综合久久| 欧美高清成人免费视频www| 亚洲国产高清在线一区二区三| 精华霜和精华液先用哪个| 少妇裸体淫交视频免费看高清| 精品人妻熟女av久视频| 晚上一个人看的免费电影| 精品久久久噜噜| 伦理电影大哥的女人| 2022亚洲国产成人精品| 七月丁香在线播放| 国产高清有码在线观看视频| 亚洲国产精品成人久久小说| 新久久久久国产一级毛片| 日日摸夜夜添夜夜添av毛片| 日韩欧美 国产精品| 久久久久久伊人网av| 国产av国产精品国产| 下体分泌物呈黄色| 涩涩av久久男人的天堂| 久久久久国产网址| 亚洲av国产av综合av卡| 九九久久精品国产亚洲av麻豆| 在线观看国产h片| 热99国产精品久久久久久7| 新久久久久国产一级毛片| 91狼人影院| 国产免费一级a男人的天堂| 免费av毛片视频| 亚洲av一区综合| 乱系列少妇在线播放| tube8黄色片| 欧美激情久久久久久爽电影| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美在线一区| 大片免费播放器 马上看| 国产黄a三级三级三级人| 王馨瑶露胸无遮挡在线观看| 国产在线一区二区三区精| 免费大片黄手机在线观看| 国产精品国产三级专区第一集| 777米奇影视久久| 欧美xxxx黑人xx丫x性爽| 亚洲精品中文字幕在线视频 | 婷婷色麻豆天堂久久| 18禁裸乳无遮挡动漫免费视频 | 国产毛片在线视频| 22中文网久久字幕| 免费不卡的大黄色大毛片视频在线观看| 六月丁香七月| 伊人久久国产一区二区| .国产精品久久| 蜜臀久久99精品久久宅男| 久久人人爽人人爽人人片va| 2021少妇久久久久久久久久久| 精品亚洲乱码少妇综合久久| 少妇高潮的动态图| 免费观看a级毛片全部| 午夜亚洲福利在线播放| 在线免费观看不下载黄p国产| 另类亚洲欧美激情| freevideosex欧美| 国国产精品蜜臀av免费| 久久久久网色| 成人午夜精彩视频在线观看| 亚洲av中文av极速乱| a级毛片免费高清观看在线播放| 午夜爱爱视频在线播放| 国产免费又黄又爽又色| 91久久精品电影网| 国产亚洲午夜精品一区二区久久 | 久久女婷五月综合色啪小说 | 亚洲欧美清纯卡通| 涩涩av久久男人的天堂| .国产精品久久| 成人国产av品久久久| 97在线人人人人妻| 少妇高潮的动态图| 婷婷色综合大香蕉| 久久精品久久精品一区二区三区| 免费看av在线观看网站| 日本熟妇午夜| 日韩欧美精品v在线| 日韩成人av中文字幕在线观看| 亚洲成人久久爱视频| 美女被艹到高潮喷水动态| 亚洲精品视频女| 在线观看国产h片| 建设人人有责人人尽责人人享有的 | 国产伦精品一区二区三区四那| 男人舔奶头视频| 久久久精品94久久精品| 国产精品久久久久久精品电影小说 | 亚洲精品第二区| 18+在线观看网站| 黄片无遮挡物在线观看| 日本爱情动作片www.在线观看| 中文字幕制服av| 夫妻性生交免费视频一级片| 99久久中文字幕三级久久日本| 好男人在线观看高清免费视频| 黄色欧美视频在线观看| 中文天堂在线官网| 丰满乱子伦码专区| 久久精品夜色国产| 免费看不卡的av| 日本与韩国留学比较| 免费观看性生交大片5| 99久久人妻综合| 人妻夜夜爽99麻豆av| 亚洲内射少妇av| 国产精品99久久99久久久不卡 | 在现免费观看毛片| 黄色欧美视频在线观看| 亚洲精品aⅴ在线观看| 激情五月婷婷亚洲| 三级国产精品片| 伊人久久精品亚洲午夜| 婷婷色综合大香蕉| 性插视频无遮挡在线免费观看| 精品亚洲乱码少妇综合久久| 在线天堂最新版资源| 亚洲精品乱久久久久久| 一本一本综合久久| 日韩一区二区视频免费看| 永久免费av网站大全| 激情 狠狠 欧美| 国产精品蜜桃在线观看| 国产男女超爽视频在线观看| 亚洲av男天堂| 国产精品嫩草影院av在线观看| 久久久久久久亚洲中文字幕| 麻豆久久精品国产亚洲av| 国产欧美另类精品又又久久亚洲欧美| 在线免费十八禁| 国产又色又爽无遮挡免| 亚洲精品成人久久久久久| av国产久精品久网站免费入址| 成人综合一区亚洲| 天天躁夜夜躁狠狠久久av| 我要看日韩黄色一级片| 国产精品蜜桃在线观看| 人妻制服诱惑在线中文字幕| 99久久精品热视频| 日本欧美国产在线视频| 波多野结衣巨乳人妻| 婷婷色综合大香蕉| 最近中文字幕高清免费大全6| 国产久久久一区二区三区| 亚洲欧美一区二区三区黑人 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 老司机影院成人| 丝袜喷水一区| 亚洲精品aⅴ在线观看| 亚洲精品国产成人久久av| 国产69精品久久久久777片| 亚洲国产日韩一区二区| 国产精品一及| 久久久久网色| 亚洲无线观看免费| 精品国产三级普通话版| 国产片特级美女逼逼视频| 免费看不卡的av| 婷婷色综合大香蕉| 久久午夜福利片| 蜜桃亚洲精品一区二区三区| 欧美日韩综合久久久久久| 美女国产视频在线观看| 久久久成人免费电影| 一区二区三区精品91| av福利片在线观看| 九九爱精品视频在线观看| av在线天堂中文字幕| 99久久中文字幕三级久久日本| 亚洲av免费在线观看| 丰满少妇做爰视频| 中文字幕久久专区| 一区二区三区精品91| 嫩草影院新地址| 熟妇人妻不卡中文字幕| 在线看a的网站| 亚洲aⅴ乱码一区二区在线播放| 人人妻人人看人人澡| 18禁裸乳无遮挡动漫免费视频 | 亚洲激情五月婷婷啪啪| 午夜老司机福利剧场| 国产成人a区在线观看| 久久久亚洲精品成人影院| 国产亚洲最大av| 如何舔出高潮| 一级毛片黄色毛片免费观看视频| 国产成人91sexporn| 久久国内精品自在自线图片| 小蜜桃在线观看免费完整版高清| 亚洲真实伦在线观看| 国产免费视频播放在线视频| 久久久久久久久久久丰满| 人妻夜夜爽99麻豆av| 久久久久久久午夜电影| 亚洲久久久久久中文字幕| 国产熟女欧美一区二区| 亚洲精品国产成人久久av| 美女视频免费永久观看网站| 欧美三级亚洲精品| 欧美性猛交╳xxx乱大交人| 又爽又黄a免费视频| 日韩av在线免费看完整版不卡| 国产免费视频播放在线视频| 在线观看人妻少妇| 亚洲一区二区三区欧美精品 | 高清视频免费观看一区二区| 亚洲国产av新网站| 国产片特级美女逼逼视频| 另类亚洲欧美激情| 国产伦精品一区二区三区视频9| 少妇 在线观看| 国产精品一区二区在线观看99| 一级a做视频免费观看| 精品99又大又爽又粗少妇毛片| 欧美日韩一区二区视频在线观看视频在线 | 成人国产麻豆网| 国产精品嫩草影院av在线观看| 免费看不卡的av| 国产精品av视频在线免费观看| 免费观看的影片在线观看| 国产精品99久久久久久久久| 亚洲av男天堂| 人妻系列 视频| 波野结衣二区三区在线| 国产成人免费观看mmmm| 久久久久国产精品人妻一区二区| 国产精品.久久久| 欧美少妇被猛烈插入视频| 伊人久久国产一区二区| av网站免费在线观看视频| 女人十人毛片免费观看3o分钟| 真实男女啪啪啪动态图| 免费高清在线观看视频在线观看| 视频中文字幕在线观看| 乱码一卡2卡4卡精品| 免费看a级黄色片| 国产亚洲av嫩草精品影院| 亚洲国产精品999| 99热网站在线观看| 麻豆精品久久久久久蜜桃| 制服丝袜香蕉在线| 日日摸夜夜添夜夜爱| 欧美变态另类bdsm刘玥| 99久久精品热视频| 日本黄色片子视频| 一级毛片黄色毛片免费观看视频| 午夜免费鲁丝| 精品一区二区三卡| 精品国产三级普通话版| 中文天堂在线官网| 欧美xxⅹ黑人| 国产黄片视频在线免费观看| 免费少妇av软件| 日韩不卡一区二区三区视频在线| 在线免费观看不下载黄p国产| 三级男女做爰猛烈吃奶摸视频| 一级二级三级毛片免费看| 亚洲在线观看片| 国产精品无大码| 王馨瑶露胸无遮挡在线观看| 午夜免费男女啪啪视频观看| 国内精品美女久久久久久| 激情 狠狠 欧美| 亚洲天堂国产精品一区在线| 2021天堂中文幕一二区在线观| 国语对白做爰xxxⅹ性视频网站| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久亚洲| 在现免费观看毛片| 日韩国内少妇激情av| 在线观看三级黄色| 国产精品久久久久久久久免| 国产成人精品久久久久久| 色视频在线一区二区三区| 精品一区二区三卡| 99热国产这里只有精品6| 一本一本综合久久| 在线精品无人区一区二区三 | 在线观看av片永久免费下载| 观看美女的网站| 精品久久久噜噜| 在线观看一区二区三区| 18禁在线无遮挡免费观看视频| 亚洲三级黄色毛片| 日本色播在线视频| 国产淫片久久久久久久久| 99re6热这里在线精品视频| 五月玫瑰六月丁香| 中文字幕亚洲精品专区| 亚洲欧美一区二区三区黑人 | 亚洲自拍偷在线| 亚洲最大成人手机在线| 国产欧美亚洲国产| 91久久精品国产一区二区三区| 精品人妻一区二区三区麻豆| 久久国产乱子免费精品| 国产人妻一区二区三区在| 黄色一级大片看看| 我要看日韩黄色一级片| 又黄又爽又刺激的免费视频.| 97人妻精品一区二区三区麻豆| 日韩亚洲欧美综合| 久久久久国产网址| 在线观看国产h片| 亚洲人成网站在线播| 午夜福利在线在线| 国产在视频线精品| 日韩,欧美,国产一区二区三区| 夫妻午夜视频| 亚洲av在线观看美女高潮| 精华霜和精华液先用哪个| 日韩一区二区三区影片| av网站免费在线观看视频| 国产亚洲av嫩草精品影院| 伦理电影大哥的女人| 超碰av人人做人人爽久久| 人人妻人人看人人澡| 精品久久国产蜜桃| 国产日韩欧美在线精品| 亚洲国产av新网站| 亚洲久久久久久中文字幕| 欧美日韩视频精品一区| 久久久久精品性色| 身体一侧抽搐| 久久国内精品自在自线图片| 亚洲经典国产精华液单| 观看美女的网站| 成年女人在线观看亚洲视频 | 免费黄网站久久成人精品| 欧美亚洲 丝袜 人妻 在线| 国产精品嫩草影院av在线观看| 亚洲欧美精品自产自拍| 日本一二三区视频观看| 国产av码专区亚洲av| 97热精品久久久久久| 精品少妇黑人巨大在线播放| 久久ye,这里只有精品| 看黄色毛片网站| 色哟哟·www| 久热这里只有精品99| 日韩av免费高清视频| 男人添女人高潮全过程视频| 国产精品人妻久久久影院| 国产一区有黄有色的免费视频| 国产成年人精品一区二区| av在线老鸭窝| 一级毛片 在线播放| 精品人妻一区二区三区麻豆| 国产精品国产三级国产av玫瑰| 丝袜喷水一区| 中文精品一卡2卡3卡4更新| 亚洲av在线观看美女高潮| 99久久中文字幕三级久久日本| 久久久a久久爽久久v久久| 亚洲av男天堂| 日本一本二区三区精品| av又黄又爽大尺度在线免费看| 成人综合一区亚洲| 日韩av免费高清视频| 国产午夜精品一二区理论片| 亚洲欧洲日产国产| 熟女av电影| 久久综合国产亚洲精品| 天天一区二区日本电影三级| 丝袜美腿在线中文| 久久久久久久久久人人人人人人| 99久久精品国产国产毛片| 亚洲精品乱码久久久久久按摩| 欧美丝袜亚洲另类| 久久99蜜桃精品久久| 欧美成人精品欧美一级黄| 少妇被粗大猛烈的视频| 国产久久久一区二区三区| 日日撸夜夜添| 精品久久久噜噜| 久久久久久国产a免费观看| 在线播放无遮挡| 啦啦啦中文免费视频观看日本| 看免费成人av毛片| 毛片一级片免费看久久久久| 插逼视频在线观看| 在线 av 中文字幕| 亚洲av中文av极速乱| 啦啦啦啦在线视频资源| 草草在线视频免费看| 日本一本二区三区精品| 日韩av免费高清视频| 麻豆成人av视频| 国产精品久久久久久av不卡| 亚洲电影在线观看av| 麻豆精品久久久久久蜜桃| 免费观看a级毛片全部| 亚洲精品乱码久久久v下载方式| 中国三级夫妇交换| 成年女人在线观看亚洲视频 | 久热这里只有精品99| 99九九线精品视频在线观看视频| 精品少妇黑人巨大在线播放| 69人妻影院| 欧美最新免费一区二区三区| 国产精品嫩草影院av在线观看| 在线观看av片永久免费下载| 久久99热6这里只有精品| 国产一区二区在线观看日韩| av在线app专区| 国产av不卡久久| 少妇的逼水好多| 免费大片18禁| av福利片在线观看|