• <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 建模方法研究
    伦理电影免费视频| 亚洲国产看品久久| 亚洲av电影在线进入| 成年人黄色毛片网站| 人妻久久中文字幕网| 99九九在线精品视频| 激情视频va一区二区三区| 国产不卡一卡二| 国内毛片毛片毛片毛片毛片| 欧美黄色淫秽网站| 一级毛片女人18水好多| 12—13女人毛片做爰片一| av视频免费观看在线观看| 国产免费福利视频在线观看| 久久这里只有精品19| 日本av手机在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 巨乳人妻的诱惑在线观看| 国产有黄有色有爽视频| 最新在线观看一区二区三区| 欧美日韩国产mv在线观看视频| 欧美亚洲日本最大视频资源| 丝瓜视频免费看黄片| 老熟女久久久| av电影中文网址| 下体分泌物呈黄色| 一区二区三区激情视频| 亚洲三区欧美一区| 精品国产亚洲在线| 丁香欧美五月| 黄色视频不卡| 母亲3免费完整高清在线观看| 久久九九热精品免费| 亚洲国产中文字幕在线视频| 窝窝影院91人妻| 亚洲午夜精品一区,二区,三区| 欧美日韩黄片免| avwww免费| 男女床上黄色一级片免费看| 国产午夜精品久久久久久| 狠狠婷婷综合久久久久久88av| 久久性视频一级片| 制服诱惑二区| 久久午夜亚洲精品久久| bbb黄色大片| 国产精品电影一区二区三区 | 老熟妇乱子伦视频在线观看| 国产淫语在线视频| 夜夜爽天天搞| 黑人欧美特级aaaaaa片| 老汉色∧v一级毛片| 亚洲男人天堂网一区| 最新的欧美精品一区二区| 夫妻午夜视频| 国产精品98久久久久久宅男小说| 精品乱码久久久久久99久播| 日本av免费视频播放| 黄片播放在线免费| 国产精品久久久久久精品古装| 欧美激情 高清一区二区三区| 久久精品成人免费网站| 国产99久久九九免费精品| 免费一级毛片在线播放高清视频 | 高清毛片免费观看视频网站 | 亚洲五月婷婷丁香| 国产成人免费无遮挡视频| 亚洲全国av大片| 日韩熟女老妇一区二区性免费视频| 99精品久久久久人妻精品| 丁香六月欧美| 搡老岳熟女国产| 国产有黄有色有爽视频| www.熟女人妻精品国产| 婷婷成人精品国产| 国产一区二区三区在线臀色熟女 | 精品熟女少妇八av免费久了| 丰满少妇做爰视频| 国产日韩欧美在线精品| 国产黄色免费在线视频| 丰满迷人的少妇在线观看| 中文亚洲av片在线观看爽 | 日韩欧美免费精品| 国产熟女午夜一区二区三区| 成人精品一区二区免费| 妹子高潮喷水视频| 亚洲精品美女久久av网站| 国产国语露脸激情在线看| 国产成人精品久久二区二区免费| 久久精品熟女亚洲av麻豆精品| 亚洲av欧美aⅴ国产| 国产精品 欧美亚洲| 精品少妇一区二区三区视频日本电影| 精品少妇黑人巨大在线播放| 欧美黄色淫秽网站| 中文亚洲av片在线观看爽 | 69精品国产乱码久久久| 我要看黄色一级片免费的| 99re在线观看精品视频| 欧美亚洲 丝袜 人妻 在线| 久久热在线av| 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 国产高清国产精品国产三级| 欧美+亚洲+日韩+国产| 日韩视频一区二区在线观看| 亚洲avbb在线观看| 飞空精品影院首页| 性色av乱码一区二区三区2| 人妻久久中文字幕网| 一二三四在线观看免费中文在| 18禁国产床啪视频网站| 精品一区二区三区四区五区乱码| 女性被躁到高潮视频| 人妻 亚洲 视频| 一区二区三区精品91| 国产99久久九九免费精品| 国产在线一区二区三区精| 午夜福利免费观看在线| 亚洲综合色网址| 少妇 在线观看| 一本一本久久a久久精品综合妖精| 久9热在线精品视频| 成人精品一区二区免费| 国产高清激情床上av| 黄片大片在线免费观看| 热re99久久国产66热| 国产精品麻豆人妻色哟哟久久| 纵有疾风起免费观看全集完整版| 久久99热这里只频精品6学生| 午夜精品国产一区二区电影| 啪啪无遮挡十八禁网站| 亚洲 欧美一区二区三区| 亚洲av成人不卡在线观看播放网| 久久免费观看电影| 脱女人内裤的视频| av在线播放免费不卡| 日韩大码丰满熟妇| 一区二区三区国产精品乱码| 少妇精品久久久久久久| 亚洲色图 男人天堂 中文字幕| 男男h啪啪无遮挡| 国产精品久久久av美女十八| 日韩大码丰满熟妇| 女人爽到高潮嗷嗷叫在线视频| 国产xxxxx性猛交| 丝袜人妻中文字幕| 亚洲午夜理论影院| av欧美777| 中文欧美无线码| 热re99久久国产66热| 日韩视频在线欧美| 在线观看免费午夜福利视频| 亚洲午夜精品一区,二区,三区| 亚洲av美国av| 国产单亲对白刺激| 欧美成人午夜精品| 丝袜美腿诱惑在线| 精品一区二区三区av网在线观看 | 一区在线观看完整版| 99久久99久久久精品蜜桃| 国产xxxxx性猛交| 国产精品成人在线| 免费日韩欧美在线观看| 十八禁网站网址无遮挡| 免费久久久久久久精品成人欧美视频| 大型黄色视频在线免费观看| a级毛片在线看网站| 窝窝影院91人妻| 在线观看一区二区三区激情| 久久香蕉激情| 757午夜福利合集在线观看| 国产精品免费一区二区三区在线 | 亚洲精品国产一区二区精华液| 另类精品久久| 高清黄色对白视频在线免费看| 青青草视频在线视频观看| 涩涩av久久男人的天堂| 制服诱惑二区| 97人妻天天添夜夜摸| 国产欧美日韩精品亚洲av| 操美女的视频在线观看| 妹子高潮喷水视频| 狠狠狠狠99中文字幕| 国产成人av激情在线播放| 国产精品影院久久| 久久午夜亚洲精品久久| www日本在线高清视频| 日本精品一区二区三区蜜桃| 国产97色在线日韩免费| 丝袜在线中文字幕| a级毛片在线看网站| 91精品国产国语对白视频| avwww免费| 欧美日韩一级在线毛片| 亚洲专区国产一区二区| 精品国产乱码久久久久久男人| 国产高清国产精品国产三级| 亚洲国产毛片av蜜桃av| 一进一出抽搐动态| 国产老妇伦熟女老妇高清| 十八禁网站网址无遮挡| 婷婷丁香在线五月| 飞空精品影院首页| 最新的欧美精品一区二区| 久久国产精品影院| 90打野战视频偷拍视频| 久久国产精品人妻蜜桃| 在线观看人妻少妇| 欧美日韩福利视频一区二区| 又紧又爽又黄一区二区| 欧美日韩亚洲综合一区二区三区_| 亚洲欧洲精品一区二区精品久久久| 嫩草影视91久久| 国产精品久久久久久人妻精品电影 | 亚洲专区字幕在线| 亚洲成av片中文字幕在线观看| 十分钟在线观看高清视频www| 黄网站色视频无遮挡免费观看| 叶爱在线成人免费视频播放| 亚洲情色 制服丝袜| tube8黄色片| 久久性视频一级片| 亚洲专区国产一区二区| 黄色视频不卡| 国产精品欧美亚洲77777| 久久人人97超碰香蕉20202| 啪啪无遮挡十八禁网站| 欧美日韩国产mv在线观看视频| 精品乱码久久久久久99久播| svipshipincom国产片| www.自偷自拍.com| 18在线观看网站| 久久精品国产综合久久久| a级片在线免费高清观看视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美日韩另类电影网站| 满18在线观看网站| 亚洲精品av麻豆狂野| 一区二区日韩欧美中文字幕| 夜夜夜夜夜久久久久| 欧美精品亚洲一区二区| 极品教师在线免费播放| 日日摸夜夜添夜夜添小说| 人妻 亚洲 视频| 桃花免费在线播放| 成人av一区二区三区在线看| 母亲3免费完整高清在线观看| 亚洲 欧美一区二区三区| 丁香欧美五月| 久久这里只有精品19| 动漫黄色视频在线观看| 丝袜人妻中文字幕| 黄色成人免费大全| 免费女性裸体啪啪无遮挡网站| 看免费av毛片| 在线观看免费视频网站a站| 国产精品秋霞免费鲁丝片| 日韩一区二区三区影片| 一级片免费观看大全| 精品久久久精品久久久| 免费高清在线观看日韩| 国产成人免费无遮挡视频| 国产一卡二卡三卡精品| 午夜久久久在线观看| 日韩欧美国产一区二区入口| 免费在线观看黄色视频的| 在线观看免费视频日本深夜| 人成视频在线观看免费观看| 我要看黄色一级片免费的| 国产三级黄色录像| 成人影院久久| 美女午夜性视频免费| 视频区图区小说| 丰满饥渴人妻一区二区三| 欧美精品一区二区大全| www日本在线高清视频| 免费高清在线观看日韩| 国产成人精品久久二区二区免费| 99国产精品免费福利视频| 国产激情久久老熟女| 男女午夜视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 欧美亚洲日本最大视频资源| 中文欧美无线码| 九色亚洲精品在线播放| 久久久久久久精品吃奶| 在线 av 中文字幕| 国产福利在线免费观看视频| 香蕉久久夜色| 国产精品一区二区在线观看99| 午夜免费鲁丝| 精品人妻熟女毛片av久久网站| 正在播放国产对白刺激| 考比视频在线观看| www.精华液| 亚洲成人免费av在线播放| 欧美日韩一级在线毛片| 妹子高潮喷水视频| 日本av免费视频播放| 女同久久另类99精品国产91| 中国美女看黄片| 成年人午夜在线观看视频| 免费高清在线观看日韩| 9191精品国产免费久久| 亚洲性夜色夜夜综合| 人人澡人人妻人| 免费黄频网站在线观看国产| 国产欧美亚洲国产| 十八禁网站免费在线| 在线观看免费日韩欧美大片| 伦理电影免费视频| 日本五十路高清| cao死你这个sao货| 色视频在线一区二区三区| 亚洲国产毛片av蜜桃av| 成年人黄色毛片网站| 青青草视频在线视频观看| 国产无遮挡羞羞视频在线观看| 国产欧美日韩精品亚洲av| 久久久久久亚洲精品国产蜜桃av| 女性被躁到高潮视频| 精品一品国产午夜福利视频| 中文欧美无线码| 黄色毛片三级朝国网站| 国产男靠女视频免费网站| 久久久久久久久久久久大奶| 中文字幕最新亚洲高清| 亚洲全国av大片| 99精国产麻豆久久婷婷| 中文字幕人妻熟女乱码| av电影中文网址| 欧美变态另类bdsm刘玥| tocl精华| 国精品久久久久久国模美| 亚洲欧美色中文字幕在线| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 午夜免费成人在线视频| 日本vs欧美在线观看视频| 欧美+亚洲+日韩+国产| 欧美日韩一级在线毛片| 伊人久久大香线蕉亚洲五| 国产黄频视频在线观看| 妹子高潮喷水视频| 美女主播在线视频| 久久久久久久大尺度免费视频| 日本av免费视频播放| 久久久精品区二区三区| 欧美日韩视频精品一区| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩福利视频一区二区| 韩国精品一区二区三区| e午夜精品久久久久久久| 十八禁高潮呻吟视频| 极品人妻少妇av视频| 国产精品影院久久| 首页视频小说图片口味搜索| 在线观看人妻少妇| 成人三级做爰电影| 两个人免费观看高清视频| svipshipincom国产片| 日韩欧美一区视频在线观看| 亚洲欧洲精品一区二区精品久久久| 久久精品国产亚洲av香蕉五月 | 真人做人爱边吃奶动态| 国产精品自产拍在线观看55亚洲 | 男女午夜视频在线观看| 可以免费在线观看a视频的电影网站| 色在线成人网| 一级片'在线观看视频| 国产日韩欧美在线精品| 国产亚洲午夜精品一区二区久久| 午夜福利一区二区在线看| av电影中文网址| 国产91精品成人一区二区三区 | 人人妻人人澡人人爽人人夜夜| 丁香六月欧美| 国产av精品麻豆| 丰满少妇做爰视频| 最近最新中文字幕大全免费视频| 九色亚洲精品在线播放| 亚洲五月婷婷丁香| 久久久久久久大尺度免费视频| 久久午夜亚洲精品久久| 18禁观看日本| 最新在线观看一区二区三区| 日本欧美视频一区| www.精华液| 国产精品国产av在线观看| 制服人妻中文乱码| 欧美日本中文国产一区发布| 王馨瑶露胸无遮挡在线观看| 中文字幕另类日韩欧美亚洲嫩草| 老司机亚洲免费影院| 99国产精品一区二区蜜桃av | 精品人妻在线不人妻| av网站在线播放免费| 国产在线观看jvid| 国产精品欧美亚洲77777| 国产不卡av网站在线观看| 91大片在线观看| 美女高潮喷水抽搐中文字幕| 日日夜夜操网爽| 91成年电影在线观看| 精品久久久精品久久久| 交换朋友夫妻互换小说| 一区二区三区激情视频| bbb黄色大片| 精品亚洲乱码少妇综合久久| 色婷婷av一区二区三区视频| 国产精品久久久av美女十八| 日韩一卡2卡3卡4卡2021年| 黄网站色视频无遮挡免费观看| 国产精品成人在线| 久久婷婷成人综合色麻豆| 亚洲黑人精品在线| 成年人午夜在线观看视频| 久久天躁狠狠躁夜夜2o2o| 建设人人有责人人尽责人人享有的| 91国产中文字幕| 久久久久久人人人人人| 国产区一区二久久| 精品一品国产午夜福利视频| 久久久久久久久免费视频了| 性高湖久久久久久久久免费观看| 欧美乱码精品一区二区三区| 51午夜福利影视在线观看| 777久久人妻少妇嫩草av网站| 亚洲精品自拍成人| 老司机靠b影院| 99热国产这里只有精品6| 老熟女久久久| 国产亚洲欧美精品永久| 蜜桃在线观看..| xxxhd国产人妻xxx| 91成年电影在线观看| 超碰成人久久| a在线观看视频网站| 999久久久精品免费观看国产| 久久香蕉激情| 一个人免费看片子| 多毛熟女@视频| 欧美黄色片欧美黄色片| 99国产精品99久久久久| 久久精品国产亚洲av高清一级| 久久国产精品影院| 精品国产乱码久久久久久小说| 国产日韩欧美亚洲二区| 黑人猛操日本美女一级片| 亚洲第一青青草原| 真人做人爱边吃奶动态| 久久99一区二区三区| 女同久久另类99精品国产91| 日韩有码中文字幕| 少妇精品久久久久久久| av天堂久久9| 啦啦啦 在线观看视频| 一区二区三区乱码不卡18| 满18在线观看网站| 999久久久精品免费观看国产| 热99re8久久精品国产| 国产成人系列免费观看| 国产精品一区二区在线不卡| 99精品欧美一区二区三区四区| 少妇猛男粗大的猛烈进出视频| 性色av乱码一区二区三区2| 少妇 在线观看| 久久久国产成人免费| 精品亚洲乱码少妇综合久久| 午夜成年电影在线免费观看| 国产激情久久老熟女| av视频免费观看在线观看| 一级毛片女人18水好多| 肉色欧美久久久久久久蜜桃| 变态另类成人亚洲欧美熟女 | videosex国产| 岛国在线观看网站| 欧美成狂野欧美在线观看| 伊人久久大香线蕉亚洲五| 精品欧美一区二区三区在线| 亚洲,欧美精品.| 真人做人爱边吃奶动态| 99国产精品一区二区蜜桃av | 啦啦啦 在线观看视频| 久久精品亚洲精品国产色婷小说| 汤姆久久久久久久影院中文字幕| av线在线观看网站| 国产精品免费大片| 岛国毛片在线播放| a级毛片黄视频| 91国产中文字幕| 午夜免费成人在线视频| 久久久精品94久久精品| 日韩三级视频一区二区三区| 欧美性长视频在线观看| 嫁个100分男人电影在线观看| 成人国产一区最新在线观看| 人妻久久中文字幕网| 99国产精品免费福利视频| 美国免费a级毛片| 久久久久国内视频| www.自偷自拍.com| svipshipincom国产片| 亚洲精品国产精品久久久不卡| 又黄又粗又硬又大视频| 国产黄色免费在线视频| 老熟女久久久| 午夜久久久在线观看| 99在线人妻在线中文字幕 | 午夜福利欧美成人| 亚洲色图综合在线观看| 亚洲精品中文字幕一二三四区 | 日韩大片免费观看网站| 国产伦理片在线播放av一区| 久久精品国产综合久久久| 1024视频免费在线观看| 日韩欧美国产一区二区入口| 丰满人妻熟妇乱又伦精品不卡| 久久久水蜜桃国产精品网| 免费在线观看视频国产中文字幕亚洲| 一区二区三区精品91| 女人高潮潮喷娇喘18禁视频| 国产aⅴ精品一区二区三区波| 人人妻人人澡人人爽人人夜夜| 国产精品98久久久久久宅男小说| 国产成人一区二区三区免费视频网站| 久久精品国产亚洲av香蕉五月 | 女人高潮潮喷娇喘18禁视频| 久久午夜综合久久蜜桃| 国内毛片毛片毛片毛片毛片| 动漫黄色视频在线观看| 少妇的丰满在线观看| 欧美激情 高清一区二区三区| 搡老熟女国产l中国老女人| 丝袜人妻中文字幕| 涩涩av久久男人的天堂| 两个人看的免费小视频| 久久人妻熟女aⅴ| 在线观看免费视频日本深夜| 女人精品久久久久毛片| 美女高潮喷水抽搐中文字幕| 不卡av一区二区三区| 欧美乱妇无乱码| 天天影视国产精品| 日日夜夜操网爽| 交换朋友夫妻互换小说| 久久人人97超碰香蕉20202| svipshipincom国产片| 午夜福利免费观看在线| 女同久久另类99精品国产91| 久久久久久亚洲精品国产蜜桃av| 99精品欧美一区二区三区四区| av不卡在线播放| 别揉我奶头~嗯~啊~动态视频| 日韩 欧美 亚洲 中文字幕| 亚洲精品av麻豆狂野| 国产精品 国内视频| 乱人伦中国视频| 久久久久精品国产欧美久久久| 国产av精品麻豆| 热re99久久精品国产66热6| 18禁裸乳无遮挡动漫免费视频| 免费在线观看完整版高清| 99国产精品一区二区三区| 久久久国产成人免费| 亚洲精品久久成人aⅴ小说| 99国产精品99久久久久| 汤姆久久久久久久影院中文字幕| 一个人免费看片子| 日日夜夜操网爽| 电影成人av| 99久久人妻综合| 国产一区二区 视频在线| 亚洲欧美一区二区三区黑人| 黄色 视频免费看| 免费久久久久久久精品成人欧美视频| 午夜福利免费观看在线| 建设人人有责人人尽责人人享有的| 成年女人毛片免费观看观看9 | 天堂中文最新版在线下载| 成人av一区二区三区在线看| 免费av中文字幕在线| 老司机影院毛片| 999久久久精品免费观看国产| 大片电影免费在线观看免费| 别揉我奶头~嗯~啊~动态视频| av欧美777| 久久精品国产综合久久久| 欧美精品一区二区免费开放| 精品视频人人做人人爽| 中文字幕高清在线视频| 免费人妻精品一区二区三区视频| 亚洲第一av免费看| 精品高清国产在线一区| 大香蕉久久成人网| 国产真人三级小视频在线观看| 在线观看免费高清a一片| 欧美亚洲 丝袜 人妻 在线| 啦啦啦免费观看视频1| 真人做人爱边吃奶动态| 嫁个100分男人电影在线观看| 久久国产精品男人的天堂亚洲| 国产免费视频播放在线视频| 啪啪无遮挡十八禁网站| 国产在线一区二区三区精| 淫妇啪啪啪对白视频| 一级片免费观看大全| 免费黄频网站在线观看国产| 国精品久久久久久国模美| 丝袜人妻中文字幕| 日韩视频在线欧美| netflix在线观看网站| 一区在线观看完整版| 国产精品国产av在线观看| 亚洲va日本ⅴa欧美va伊人久久|