fluent多孔介質(zhì)模型課件

多孔介質(zhì)模型多孔介質(zhì)是由多相物質(zhì)所占據(jù)的共同空間，也是多相物質(zhì)共存的一種組合體，沒有固體骨架的那部分空間叫做孔隙，由液體或氣體或氣液兩相共同占有，相對于其中一相來說，其他相都彌散在其中，并以固相為固體骨架，構(gòu)成空隙空間的某些空洞相互連通。多孔介質(zhì)模型可以應(yīng)用于很多問題，如通過充滿介質(zhì)的流動、通過過濾紙、穿孔圓盤、流量分配器以及管道堆的流動。1多孔介質(zhì)模型多孔介質(zhì)是由多相物質(zhì)所占據(jù)的共同空間，也動量方程能量方程的處理阻力系數(shù)的推導(dǎo)操作步驟（實例）后處理

多孔介質(zhì)模擬方法是將流動區(qū)域中固體結(jié)構(gòu)的作用看作是附加在流體上的分布阻力。多孔介質(zhì)模型應(yīng)用

基于多孔介質(zhì)模型的數(shù)值模擬方法最初被用于模擬換熱器和核反應(yīng)堆中流體的流動和傳熱問題。換熱器中存在大量的換熱管道和阻礙片,要模擬具有500根換熱管和10個阻礙片的換熱器中的流動,將需要1.5億個網(wǎng)格單元,這大大超出了目前計算機的計算能力。為此,PatankarSpalding提出了采用分布阻力的方法,也稱為多孔介質(zhì)模型的方法。之后sha等采用這種方法模擬了蒸汽發(fā)生器和核反應(yīng)堆堆芯中流體的流動,Karayannis等模擬了換熱器中的流動,Prithiviraj和Andrews模擬了三維換熱器中的流動。我們?yōu)槭裁匆獞?yīng)用多孔介質(zhì)模型？2動量方程多孔介質(zhì)模擬方法是將流動區(qū)域中固體結(jié)構(gòu)的作計算流體力學控制方程方程變數(shù)擴散系數(shù)源項連續(xù)性方程100X-動量方程uY-動量方程vZ-動量方程w能量方程T時間項對流項擴散項廣義源項3計算流體力學控制方程方程變數(shù)擴散系數(shù)源項連續(xù)性方程1動量方程通用形式：幫助文件：流體力學：運動方程：慣性力品質(zhì)力表面力

多孔介質(zhì)模擬方法是將流動區(qū)域中固體結(jié)構(gòu)的作用看作是附加在流體上的分布阻力。isexternalbodyforces.alsocontainsothermodel-dependentsourcetermssuchasporous-mediaanduser-definedsources.4動量方程通用形式：幫助文件：流體力學：運動方程：慣性力品質(zhì)力多孔介質(zhì)的源項

多孔介質(zhì)的作用是在動量方程中增加一個源項來模擬，源項由兩部分組成：一個粘性損失項和一個慣性損失項。其中Si是i向(x,y,z)動量源項，D和C是規(guī)定的矩陣。在多孔介質(zhì)單元中，動量損失對于壓力梯度有貢獻，壓降和流體速度（或速度平方）成比例。對于簡單的多孔介質(zhì)，只保留D和C矩陣的對角線元素。對于簡單的、各向同性的多孔介質(zhì)，各個方向的阻力特性一樣，對角線元素相等。5多孔介質(zhì)的源項多孔介質(zhì)的作用是在動量方程中增加一個源多孔介質(zhì)的源項在多孔介質(zhì)單元中，動量損失對于壓力梯度有貢獻，壓降和流體速度（或速度平方）成比例。

多孔介質(zhì)的Darcy定律：通過多孔介質(zhì)的層流流動中，壓降和速度成比例，常數(shù)C2

可以考慮為零。忽略對流加速以及擴散。

△Px,△Py,△Pz分別是x,y,z三個方向的壓力降?！鱪x,

△ny,

△nz分別是多孔介質(zhì)在x,y,z三個方向的真實厚度。6多孔介質(zhì)的源項在多孔介質(zhì)單元中，動量損失對于壓力梯度多孔介質(zhì)的源項多孔中的慣性損失：在高速流動中，多孔介質(zhì)動量源項中的常數(shù)C2

可以對慣性損失作出修正。C2可以看成沿著流動方向每一單位長度的損失系數(shù)，因此允許壓降指定為動壓頭的函數(shù)。

如果模擬的是穿孔板或者管道堆，有時可以消除滲透項而只是用內(nèi)部損失項，從而得到下面的多孔介質(zhì)簡化方程：

△Px,△Py,△Pz分別是x,y,z三個方向的壓力降?！鱪x,

△ny,

△nz分別是多孔介質(zhì)在x,y,z三個方向的真實厚度。7多孔介質(zhì)的源項多孔中的慣性損失：在高速流動中，多孔介能量方程的處理能量方程：多孔介質(zhì)對能量方程修正：對于多孔介質(zhì)流動，F(xiàn)LUENT仍然解標準能量輸運方程，只是修改了對流項和時間導(dǎo)數(shù)項。對對流項的計算采用了有效對流函數(shù)，時間導(dǎo)數(shù)項則計入了固體區(qū)域?qū)Χ嗫捉橘|(zhì)的熱慣性效應(yīng)。多孔區(qū)域的有效熱傳導(dǎo)率keff是由流體的熱傳導(dǎo)率和固體的熱傳導(dǎo)率的體積平均值計算得到：8能量方程的處理能量方程：多孔介質(zhì)對能量方程修正：對于阻力系數(shù)推導(dǎo)采用已知的壓強損失，基于名義速度推導(dǎo)多孔介質(zhì)參數(shù)采用Ergun方程計算充填床的多孔介質(zhì)參數(shù)用經(jīng)驗公式推導(dǎo)穿過多孔板的流動參數(shù)用表格數(shù)據(jù)計算流過纖維墊的多孔介質(zhì)參數(shù)用壓強和速度的實驗數(shù)據(jù)計算多孔介質(zhì)系數(shù)9阻力系數(shù)推導(dǎo)采用已知的壓強損失，基于名義速度推導(dǎo)多孔介質(zhì)參數(shù)采用Ergun方程計算充填床的多孔介質(zhì)參數(shù)在湍流中，充填床的數(shù)學模型是用穿透率和慣性損失系數(shù)來定義的。計算相關(guān)常數(shù)的一種辦法是使用半經(jīng)驗公式Ergun方程，這個方程適用的雷諾數(shù)范圍很廣，同時也使用于多種填充物：在計算層流時，上式中右端第二項可以被去掉，Ergun方程隨之簡化為Blake-Kozeny方程：在上述方程中μ為粘度，Dp為粒子的平均直徑，L為床的深度,ε為空腔比率,定義為空腔與填充床的體積比。對比多孔介質(zhì)的源項方程，各方向上的滲透率和慣性損失系數(shù)為：10采用Ergun方程計算充填床的多孔介質(zhì)參數(shù)在湍流中，用經(jīng)驗公式推導(dǎo)穿過多孔板的流動參數(shù)用VanWinkle方程計算帶方孔的多孔板上壓強的損失。方程的提出者認為該方程適用于呈三角形分布的等距方孔板的湍流計算，具體形式如下：式中m為通過板的質(zhì)量流量，fA為孔的總面積，pA板的總面積(固體與孔的和)，D/t孔直徑與板厚之比,C是隨雷諾數(shù)和D/t變化的系數(shù)，其值可以通過查表獲得。在t/D>1.6,且Re>4000時，C近似等于0.98，其中雷諾數(shù)是用孔的直徑做特征長，孔中流體的速度做特征速度求出的。將方程

代入上式，并除以板厚Δx=t可得：其中v是名義速度。與多孔介質(zhì)源項方程對比后可以發(fā)現(xiàn)：11用經(jīng)驗公式推導(dǎo)穿過多孔板的流動參數(shù)用VanWink用表格數(shù)據(jù)計算流過纖維墊的多孔介質(zhì)參數(shù)在計算流體流過纖維墊或過濾器這類問題時，除了Blake-Kozeny方程外，還可以用試驗數(shù)據(jù)進行計算。固體材料的體積比纖維的無量綱滲透率B0.2620.250.2580.260.2210.400.2180.410.1720.80表中，a為纖維直徑。α的含義與darcy定律的方程中相同，對于給定纖維直徑和體積比的多孔介質(zhì)是容易求出的。12用表格數(shù)據(jù)計算流過纖維墊的多孔介質(zhì)參數(shù)在計算流體流過纖維墊或用壓強和速度的實驗數(shù)據(jù)計算多孔介質(zhì)系數(shù)

在已知多孔介質(zhì)上的速度與壓強降的試驗數(shù)據(jù)時，可以通過插值求出多孔介質(zhì)上系數(shù)。速度（m/s）壓強降（Pa）20.078.050.0487.080.01432.0110.02964.0采用上表的數(shù)據(jù)可以擬合出一條“速度－壓強降”曲線，其方程為：對比上述兩式便可求出粘性阻力系數(shù)和慣性阻力系數(shù)。13用壓強和速度的實驗數(shù)據(jù)計算多孔介質(zhì)系數(shù)在已知多孔介質(zhì)實例計算上圖中的計算區(qū)域尺寸如下：總的計算域：長1m，寬0.1m；Poroustwo：長0.57m，寬0.02m；（處于正中間）Porousone：寬0.03m，高0.06m；Porousthree：寬0.03m，高0.06m；

邊界條件如上圖中所示，進口取velocityinlet，速度為0.01m/s；出口取pressureoutlet，壓力值為大氣壓。三個多孔介質(zhì)區(qū)中，porousone和porousthree的性質(zhì)一樣。Poroustwo的粘性阻力系數(shù)為1e+10，其余多孔介質(zhì)區(qū)為1e+13.由于是低速層流流動，不考慮慣性阻力的影響。進口出口WALLsymmetryPorousonePoroustwoPorousthree14實例計算上圖中的計算區(qū)域尺寸如下：進口出口WALLsymmeFluent中設(shè)置在GAMBIT中將多孔區(qū)單獨設(shè)置，但其性質(zhì)仍為fluid.在fluent的邊界條件設(shè)置多孔區(qū)的參數(shù)，方向設(shè)置如下圖。多孔區(qū)poroustwo的粘性阻力設(shè)為1e+10；其余多孔區(qū)粘性阻力設(shè)為1e+13，如右邊兩圖所示。15Fluent中設(shè)置在GAMBIT中將多孔區(qū)單獨設(shè)置，但多孔介質(zhì)的后處理在多孔介質(zhì)區(qū)域，由于粘性阻力的存在，流體在多孔區(qū)內(nèi)有較大的壓降如第一圖所示；poroustwo的粘性阻力系數(shù)是其他多孔區(qū)的千分之一，故流體幾乎不會通過porousone和porousthree，而全部由poroustwo通過，如第二圖和第三圖所示。16多孔介質(zhì)的后處理在多孔介質(zhì)區(qū)域，由于粘性阻力的存在，多孔介質(zhì)的后處理上圖是x=0.215、x=0.5、x=0.785三條線上的軸向速度分布：可以看出在多孔區(qū)中，軸向速度都為0.X=0.215X=0.5X=0.78517多孔介質(zhì)的后處理上圖是x=0.215、x=0.5、x=0.7多孔介質(zhì)的后處理Y=0.05上圖是y=0.05線上的徑向速度分布：可以看出，在porousone和porousthree中徑向速度都為0；由于porousone的阻擋，其前面有大量的流體往上流，而在porousthree后面則有一個回流。18多孔介質(zhì)的后處理Y=0.05上圖是y=0.05線上的徑向速度催化排氣凈化器汽車等的催化排氣凈化器（catalyticconverter）簡化模型如下圖：進口多孔區(qū)出口催化轉(zhuǎn)化器通常用于凈化天然氣和柴油發(fā)動機排除的廢氣，廢氣中包含一氧化碳，氧化氮以及未燃燒的碳氫燃料等有害物質(zhì)。排除的廢氣在催化凈化反應(yīng)器中通過一個涂有金屬催化劑的陶瓷結(jié)構(gòu)多孔區(qū)，所以在多孔區(qū)中的壓力梯度和速度分布對廢氣的凈化效率有很大的影響。因此CFD分析用于設(shè)計高效的催化轉(zhuǎn)化器。本文模擬的催化轉(zhuǎn)化器進出口兩端對稱，總長260mm.多孔區(qū)長80mm，直徑100mm，進出口直徑30mm。19催化排氣凈化器汽車等的催化排氣凈化器（catalyt網(wǎng)格劃分與計算設(shè)置在GAMBIT里面全部劃分為六面體網(wǎng)格，如下圖所示，將計算區(qū)域分為兩部分，中間黃色區(qū)域為多孔區(qū)。

在fluent中設(shè)置為湍流計算，在多孔區(qū)采用層流計算，并設(shè)Y,Z兩個方向的阻力系數(shù)較X方向大三個數(shù)量級，如右兩圖所示：20網(wǎng)格劃分與計算設(shè)置在GAMBIT里面全部劃分為六面體計算結(jié)果上圖Y=0截面上的速度矢量。由圖可以看出，氣流進入凈化器后在氣流兩邊形成了回流，而在多孔區(qū)中的流動變得非常均勻且只有X方向的流動。這樣金屬催化劑能起到更好的作用。21計算結(jié)果上圖Y=0截面上的速度矢量。由圖可以看出，氣計算結(jié)果上圖為x=95、130、165截面的X方向速度云圖?？梢钥闯觯瑲饬魉俣仍谕ㄟ^多孔區(qū)后變得更均勻。在通過多孔區(qū)前，中心區(qū)域(紅色區(qū))的速度值非常的大,而在通過多孔區(qū)后明顯減小。在整個界面上，通過多孔區(qū)之后速度變得更均勻。22計算結(jié)果上圖為x=95、130、165截面的X方向速計算結(jié)果上圖y=0截面的壓力云圖。可以看出，當氣流通過多孔區(qū)，速度變化大的中心區(qū)域，壓力變化也非常大。這是由于多孔區(qū)的粘性阻力和慣性阻力的影響，壓降很大。23計算結(jié)果上圖y=0截面的壓力云圖。可以看出，當氣流通計算結(jié)果上圖為在多孔區(qū)內(nèi)，沿中心線的壓強變化?？梢钥闯?，穿過多孔區(qū)的壓力降約為450Pa.24計算結(jié)果上圖為在多孔區(qū)內(nèi)，沿中心線的壓強變化。可以看謝謝觀賞！25謝謝觀賞！25多孔介質(zhì)模型多孔介質(zhì)是由多相物質(zhì)所占據(jù)的共同空間，也是多相物質(zhì)共存的一種組合體，沒有固體骨架的那部分空間叫做孔隙，由液體或氣體或氣液兩相共同占有，相對于其中一相來說，其他相都彌散在其中，并以固相為固體骨架，構(gòu)成空隙空間的某些空洞相互連通。多孔介質(zhì)模型可以應(yīng)用于很多問題，如通過充滿介質(zhì)的流動、通過過濾紙、穿孔圓盤、流量分配器以及管道堆的流動。26多孔介質(zhì)模型多孔介質(zhì)是由多相物質(zhì)所占據(jù)的共同空間，也動量方程能量方程的處理阻力系數(shù)的推導(dǎo)操作步驟（實例）后處理

多孔介質(zhì)模擬方法是將流動區(qū)域中固體結(jié)構(gòu)的作用看作是附加在流體上的分布阻力。多孔介質(zhì)模型應(yīng)用

基于多孔介質(zhì)模型的數(shù)值模擬方法最初被用于模擬換熱器和核反應(yīng)堆中流體的流動和傳熱問題。換熱器中存在大量的換熱管道和阻礙片,要模擬具有500根換熱管和10個阻礙片的換熱器中的流動,將需要1.5億個網(wǎng)格單元,這大大超出了目前計算機的計算能力。為此,PatankarSpalding提出了采用分布阻力的方法,也稱為多孔介質(zhì)模型的方法。之后sha等采用這種方法模擬了蒸汽發(fā)生器和核反應(yīng)堆堆芯中流體的流動,Karayannis等模擬了換熱器中的流動,Prithiviraj和Andrews模擬了三維換熱器中的流動。我們?yōu)槭裁匆獞?yīng)用多孔介質(zhì)模型？27動量方程多孔介質(zhì)模擬方法是將流動區(qū)域中固體結(jié)構(gòu)的作計算流體力學控制方程方程變數(shù)擴散系數(shù)源項連續(xù)性方程100X-動量方程uY-動量方程vZ-動量方程w能量方程T時間項對流項擴散項廣義源項28計算流體力學控制方程方程變數(shù)擴散系數(shù)源項連續(xù)性方程1動量方程通用形式：幫助文件：流體力學：運動方程：慣性力品質(zhì)力表面力

多孔介質(zhì)模擬方法是將流動區(qū)域中固體結(jié)構(gòu)的作用看作是附加在流體上的分布阻力。isexternalbodyforces.alsocontainsothermodel-dependentsourcetermssuchasporous-mediaanduser-definedsources.29動量方程通用形式：幫助文件：流體力學：運動方程：慣性力品質(zhì)力多孔介質(zhì)的源項

多孔介質(zhì)的作用是在動量方程中增加一個源項來模擬，源項由兩部分組成：一個粘性損失項和一個慣性損失項。其中Si是i向(x,y,z)動量源項，D和C是規(guī)定的矩陣。在多孔介質(zhì)單元中，動量損失對于壓力梯度有貢獻，壓降和流體速度（或速度平方）成比例。對于簡單的多孔介質(zhì)，只保留D和C矩陣的對角線元素。對于簡單的、各向同性的多孔介質(zhì)，各個方向的阻力特性一樣，對角線元素相等。30多孔介質(zhì)的源項多孔介質(zhì)的作用是在動量方程中增加一個源多孔介質(zhì)的源項在多孔介質(zhì)單元中，動量損失對于壓力梯度有貢獻，壓降和流體速度（或速度平方）成比例。

多孔介質(zhì)的Darcy定律：通過多孔介質(zhì)的層流流動中，壓降和速度成比例，常數(shù)C2

可以考慮為零。忽略對流加速以及擴散。

△Px,△Py,△Pz分別是x,y,z三個方向的壓力降?！鱪x,

△ny,

△nz分別是多孔介質(zhì)在x,y,z三個方向的真實厚度。31多孔介質(zhì)的源項在多孔介質(zhì)單元中，動量損失對于壓力梯度多孔介質(zhì)的源項多孔中的慣性損失：在高速流動中，多孔介質(zhì)動量源項中的常數(shù)C2

可以對慣性損失作出修正。C2可以看成沿著流動方向每一單位長度的損失系數(shù)，因此允許壓降指定為動壓頭的函數(shù)。

如果模擬的是穿孔板或者管道堆，有時可以消除滲透項而只是用內(nèi)部損失項，從而得到下面的多孔介質(zhì)簡化方程：

△Px,△Py,△Pz分別是x,y,z三個方向的壓力降?！鱪x,

△ny,

△nz分別是多孔介質(zhì)在x,y,z三個方向的真實厚度。32多孔介質(zhì)的源項多孔中的慣性損失：在高速流動中，多孔介能量方程的處理能量方程：多孔介質(zhì)對能量方程修正：對于多孔介質(zhì)流動，F(xiàn)LUENT仍然解標準能量輸運方程，只是修改了對流項和時間導(dǎo)數(shù)項。對對流項的計算采用了有效對流函數(shù)，時間導(dǎo)數(shù)項則計入了固體區(qū)域?qū)Χ嗫捉橘|(zhì)的熱慣性效應(yīng)。多孔區(qū)域的有效熱傳導(dǎo)率keff是由流體的熱傳導(dǎo)率和固體的熱傳導(dǎo)率的體積平均值計算得到：33能量方程的處理能量方程：多孔介質(zhì)對能量方程修正：對于阻力系數(shù)推導(dǎo)采用已知的壓強損失，基于名義速度推導(dǎo)多孔介質(zhì)參數(shù)采用Ergun方程計算充填床的多孔介質(zhì)參數(shù)用經(jīng)驗公式推導(dǎo)穿過多孔板的流動參數(shù)用表格數(shù)據(jù)計算流過纖維墊的多孔介質(zhì)參數(shù)用壓強和速度的實驗數(shù)據(jù)計算多孔介質(zhì)系數(shù)34阻力系數(shù)推導(dǎo)采用已知的壓強損失，基于名義速度推導(dǎo)多孔介質(zhì)參數(shù)采用Ergun方程計算充填床的多孔介質(zhì)參數(shù)在湍流中，充填床的數(shù)學模型是用穿透率和慣性損失系數(shù)來定義的。計算相關(guān)常數(shù)的一種辦法是使用半經(jīng)驗公式Ergun方程，這個方程適用的雷諾數(shù)范圍很廣，同時也使用于多種填充物：在計算層流時，上式中右端第二項可以被去掉，Ergun方程隨之簡化為Blake-Kozeny方程：在上述方程中μ為粘度，Dp為粒子的平均直徑，L為床的深度,ε為空腔比率,定義為空腔與填充床的體積比。對比多孔介質(zhì)的源項方程，各方向上的滲透率和慣性損失系數(shù)為：35采用Ergun方程計算充填床的多孔介質(zhì)參數(shù)在湍流中，用經(jīng)驗公式推導(dǎo)穿過多孔板的流動參數(shù)用VanWinkle方程計算帶方孔的多孔板上壓強的損失。方程的提出者認為該方程適用于呈三角形分布的等距方孔板的湍流計算，具體形式如下：式中m為通過板的質(zhì)量流量，fA為孔的總面積，pA板的總面積(固體與孔的和)，D/t孔直徑與板厚之比,C是隨雷諾數(shù)和D/t變化的系數(shù)，其值可以通過查表獲得。在t/D>1.6,且Re>4000時，C近似等于0.98，其中雷諾數(shù)是用孔的直徑做特征長，孔中流體的速度做特征速度求出的。將方程

代入上式，并除以板厚Δx=t可得：其中v是名義速度。與多孔介質(zhì)源項方程對比后可以發(fā)現(xiàn)：36用經(jīng)驗公式推導(dǎo)穿過多孔板的流動參數(shù)用VanWink用表格數(shù)據(jù)計算流過纖維墊的多孔介質(zhì)參數(shù)在計算流體流過纖維墊或過濾器這類問題時，除了Blake-Kozeny方程外，還可以用試驗數(shù)據(jù)進行計算。固體材料的體積比纖維的無量綱滲透率B0.2620.250.2580.260.2210.400.2180.410.1720.80表中，a為纖維直徑。α的含義與darcy定律的方程中相同，對于給定纖維直徑和體積比的多孔介質(zhì)是容易求出的。37用表格數(shù)據(jù)計算流過纖維墊的多孔介質(zhì)參數(shù)在計算流體流過纖維墊或用壓強和速度的實驗數(shù)據(jù)計算多孔介質(zhì)系數(shù)

在已知多孔介質(zhì)上的速度與壓強降的試驗數(shù)據(jù)時，可以通過插值求出多孔介質(zhì)上系數(shù)。速度（m/s）壓強降（Pa）20.078.050.0487.080.01432.0110.02964.0采用上表的數(shù)據(jù)可以擬合出一條“速度－壓強降”曲線，其方程為：對比上述兩式便可求出粘性阻力系數(shù)和慣性阻力系數(shù)。38用壓強和速度的實驗數(shù)據(jù)計算多孔介質(zhì)系數(shù)在已知多孔介質(zhì)實例計算上圖中的計算區(qū)域尺寸如下：總的計算域：長1m，寬0.1m；Poroustwo：長0.57m，寬0.02m；（處于正中間）Porousone：寬0.03m，高0.06m；Porousthree：寬0.03m，高0.06m；

邊界條件如上圖中所示，進口取velocityinlet，速度為0.01m/s；出口取pressureoutlet，壓力值為大氣壓。三個多孔介質(zhì)區(qū)中，porousone和porousthree的性質(zhì)一樣。Poroustwo的粘性阻力系數(shù)為1e+10，其余多孔介質(zhì)區(qū)為1e+13.由于是低速層流流動，不考慮慣性阻力的影響。進口出口WALLsymmetryPorousonePoroustwoPorousthree39實例計算上圖中的計算區(qū)域尺寸如下：進口出口WALLsymmeFluent中設(shè)置在GAMBIT中將多孔區(qū)單獨設(shè)置，但其性質(zhì)仍為fluid.在fluent的邊界條件設(shè)置多孔區(qū)的參數(shù)，方向設(shè)置如下圖。多孔區(qū)poroustwo的粘性阻力設(shè)為1e+10；其余多孔區(qū)粘性阻力設(shè)為1e+13，如右邊兩圖所示。40Fluent中設(shè)置在GAMBIT中將多孔區(qū)單獨設(shè)置，但多孔介質(zhì)的后處理在多孔介質(zhì)區(qū)域，由于粘性阻力的存在，流體在多孔區(qū)內(nèi)有較大的壓降如第一圖所示；poroustwo的粘性阻力系數(shù)是其他多孔區(qū)的千分之一，故流體幾乎不會通過porousone和porousthree，而全部由poroustwo通過，如第二圖和第三圖所示。41多孔介質(zhì)的后處理在多孔介質(zhì)區(qū)域，由于粘性阻力的存在，多孔介質(zhì)的后處理上圖是x=0.215、x=0.5、x=0.785三條線上的軸向速度分布：可以看出在多孔區(qū)中，軸向速度都為0.X=0.215X=0.5X=0.78542多孔介質(zhì)的后處理上圖是x=0.215、x=0.5、x=0.7多孔介質(zhì)的后處理Y=0.05上圖是y=0.05線上的徑向速度分布：可以看出，在porousone和porousthree中徑向速度都為0；由于porou