突擴(kuò)管中高分子聚合物溶液湍流減阻的數(shù)值模擬

石若冉，張浩，2，*，張志，李茂林

（1.山東建筑大學(xué)熱能工程學(xué)院，山東 濟(jì)南 250101；2.山東省綠色建筑協(xié)同創(chuàng)新中心，山東 濟(jì)南 250101）

0 引言

供暖和制冷的實(shí)現(xiàn)建立在流體運(yùn)輸?shù)幕A(chǔ)上，流動帶來的阻力不僅會增加能耗，而且會使管道壁面受到磨損，因而減少管道輸運(yùn)的阻力以達(dá)到節(jié)能性、經(jīng)濟(jì)性的目的一直是重要的研究方向。添加劑減阻指的是在管道流體中加入高分子聚合物或者表面活性劑來降低管道阻力的技術(shù)。早在1948 年，TOMS［1］就發(fā)現(xiàn)高分子聚合物可使湍流流動阻力大幅減小，該現(xiàn)象后被稱為“TOMS”效應(yīng)。自此學(xué)者們便開始了針對添加劑減阻方向的研究，為降低能耗和提高經(jīng)濟(jì)效益做出貢獻(xiàn)。在實(shí)際工程應(yīng)用方面，1979年聚合物減阻在美國阿拉斯加的石油輸運(yùn)行業(yè)中實(shí)現(xiàn)商業(yè)應(yīng)用，在管道中添加減阻劑取得了良好的減阻效果［2］。為實(shí)現(xiàn)減阻節(jié)能，焦利芳等［3］開展減阻劑實(shí)際應(yīng)用研究，在青島某小區(qū)供熱系統(tǒng)中加入表面活性劑十六烷基三甲基氯化銨／水楊酸鈉（CTAC／NaSal），其效果理想［3］。這些應(yīng)用實(shí)驗(yàn)充分證明了添加劑減阻技術(shù)實(shí)際應(yīng)用方面的可行性，對其進(jìn)行深入研究勢必將在節(jié)能領(lǐng)域產(chǎn)生顯著的經(jīng)濟(jì)效益。

針對添加劑減阻方式，學(xué)者們從流動特性、減阻機(jī)理等方面開展了研究，且隨著計(jì)算機(jī)科學(xué)的發(fā)展，數(shù)值模擬也逐漸成為有力的工具。張紅霞［4］分別利用純黏性剪切稀化Carreau-Bird 模型與Giesekus黏彈性非牛頓流體本構(gòu)方程模擬了CTAC 水溶液，探究其流動與減阻作用。李昌烽等［5］利用Cross模型模擬了聚丙烯酰胺（Polyacrylamide，PAM）溶液在二維軸對稱圓管中的流場和溫度場，分析了其減阻效果、傳熱性能。朱涵［6］采用Cross 模型分別模擬了減阻劑溶液的水下彈體與魚形繞流減阻流動，分析了其減阻效果。李恩田［7］利用Carreau-Bird模型與冪律模型對CTAC減阻溶液、利用Cross流變模型對PAM溶液分別進(jìn)行數(shù)值模擬，并結(jié)合壁面微溝槽探究了聯(lián)合減阻。這些研究闡述了表面活性劑的減阻機(jī)理和規(guī)律，也充分證明了利用數(shù)值模擬研究減阻的可行性。

目前，關(guān)于添加劑減阻研究多數(shù)集中在直管，而實(shí)際管道系統(tǒng)中存在著很多種類的非直管管件，突擴(kuò)管作為其中之一，其減阻研究在理論認(rèn)識或?qū)嶋H工程應(yīng)用都具有很高的價值。PAK 等［8］分析了聚丙烯酰胺AP-273水溶液在不同突擴(kuò)比的突擴(kuò)管中的局部阻力損失特性。POOLE 等［9］利用實(shí)驗(yàn)研究了不同濃度的聚丙烯酸（Polyacrylic Acid，PAA）溶液在面積比1∶4 突擴(kuò)管中的湍流流動，探究了其流動的不對稱性和再附著長度的變化。國內(nèi)學(xué)者利用實(shí)驗(yàn)方法較多，且多數(shù)采用表面活性劑作為添加劑。焦利芳等［10］實(shí)驗(yàn)研究了CTAC／NaSal 減阻溶液在不同種類不規(guī)則管道（包括90°彎頭、三通以及變徑管）內(nèi)的流動和減阻特性，發(fā)現(xiàn)其減阻效果比起直管段有著減弱現(xiàn)象。段傳偉［11］利用實(shí)驗(yàn)法得到不同管徑比的突擴(kuò)管、突縮管中不同濃度的十六烷基三甲基溴化銨（Cetyltrimethylammonium Bromide，CTAB）溶液的流動，研究了壓力分布、阻力發(fā)展長度、減阻效果等。蔡書鵬等［12-13］研究了管徑比為1∶1.52及截面積比為1∶4 的突擴(kuò)管中兩種不同質(zhì)量分?jǐn)?shù)的CTAB水溶液的流動減阻行為，探究了局部阻力系數(shù)的改變及下游壓力恢復(fù)長度等問題。王侃宏等［14］在管徑比為1∶1.49 的突擴(kuò)管中使用了3 種不同濃度CTAB水溶液，研究了上下游管道的沿程阻力系數(shù)以及突擴(kuò)處的局部阻力系數(shù)的變化與減阻規(guī)律。

高分子聚合物和表面活性劑的主要不同之處在于，高分子聚合物起減阻作用的柔性長鏈結(jié)構(gòu)相較于表面活性劑的膠束結(jié)構(gòu)，前者在強(qiáng)剪切應(yīng)力作用下會被破壞，且無法自我修復(fù)，而后者需達(dá)到一定濃度才可實(shí)現(xiàn)減阻，并且表面活性劑的添加會對管壁造成不利影響［15］。文章選擇高分子聚合物PAM溶液作為研究對象，利用Cross模型描述其黏度變化，與純水工況作對比，在管徑比1∶1.5 的三維突擴(kuò)圓管模型中分析其流場分布和不同雷諾數(shù)下直管段及突擴(kuò)部分的減阻效果、回流區(qū)的變化。

1 數(shù)值模擬方法

1.1 幾何模型

建立的突擴(kuò)管模型如圖1所示，設(shè)置上、下游管管徑（d1、d2）分別為10、15 mm。為使得管內(nèi)流體充分發(fā)展，兩段管長分別設(shè)為20 倍自身管徑長度，即上、下游管長（l1、l2）分別為200、300 mm。

圖1 突擴(kuò)管幾何模型

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

文章的流體流動研究不考慮傳熱，所需控制方程張量形式，即連續(xù)性方程和N-S 方程分別由式（1）和（2）［7］表示為

式中μ為流體動力黏度，Pa·s；p為壓強(qiáng)，Pa；fi為流體質(zhì)點(diǎn)的單位質(zhì)量力，m／s2；ρ為流體密度，kg／m3。

湍流模型中的RNGk-ε模型適用于渦流的描述，對壁面的處理更準(zhǔn)確。考慮到突擴(kuò)管流動的特點(diǎn)，選其作為計(jì)算模型。RNGk-ε模型中湍動能（k）方程和湍動能耗散（ε）方程分別由式（3）和（4）［15-16］表示為

高分子聚合物減阻溶液屬于非牛頓流體，相較于純水而言，其溶液具有剪切稀化特性［7］，即黏度隨著剪切速率的增加而減小，因此需要建立本構(gòu)方程描述黏度。Cross 模型是其中一種用于描述非牛頓流體剪切稀化特征的流變模型，可由式（5）［18］表示為

式中λ為松弛時間，s；n為冪指數(shù)常數(shù)；μ0為零剪切黏度，Pa·s；為剪切速率，s-1。參考文獻(xiàn)［7］和［18］的研究，為描述質(zhì)量濃度為15 mg／L 聚丙烯酰胺溶液的動力黏度，參數(shù)μ0＝3.835 mPa·s、λ＝2.746 s、n＝0.7。

1.3 模擬設(shè)置

利用數(shù)值模擬法模擬突擴(kuò)管中純水和高分子聚合物PAM溶液的流動，二者物性參數(shù)設(shè)置如下：當(dāng)模型工況為純水時，動力黏度為0.001 kg／（m·s），密度為1 000 kg／m3；當(dāng)模型工況為PAM溶液時，黏度描述采用Cross模型，而密度設(shè)置同純水。

模擬采用三維穩(wěn)態(tài)的方式，忽略重力影響；湍流模型選擇RNGk-ε；邊界條件中入口采用velocityinlet，通過設(shè)置入口平均速度大小實(shí)現(xiàn)不同雷諾數(shù)的湍流工況，出口采用pressure-outlet，壁面為無滑移條件。

2 網(wǎng)格無關(guān)性驗(yàn)證與可行性分析

2.1 網(wǎng)格無關(guān)性驗(yàn)證

突擴(kuò)管幾何模型的網(wǎng)格劃分，因其結(jié)構(gòu)較簡單，采用適用圓柱體的o-block 結(jié)構(gòu)化網(wǎng)格劃分的方式。設(shè)置了網(wǎng)格數(shù)分別為42 萬、96 萬、121 萬的3種網(wǎng)格，在雷諾數(shù)為10 000 的純水工況下，對比突擴(kuò)管中軸線上的靜壓分布，從而驗(yàn)證網(wǎng)格無關(guān)性。不同網(wǎng)格數(shù)目對應(yīng)的靜壓分布如圖2 所示，靜壓分布趨勢一致，上游管在0 ～0.2 m 范圍內(nèi)下降至最小值，經(jīng)突擴(kuò)部分壓力先上升后又呈下降趨勢。在上游管段內(nèi)，96 萬與121 萬網(wǎng)格的靜壓分布高度重合，壓力恢復(fù)最大值也極為接近，可知42 萬網(wǎng)格計(jì)算所得結(jié)果并不精準(zhǔn)，網(wǎng)格數(shù)至96萬再增加網(wǎng)格數(shù)量壓力分布基本不再變化，因此選取96萬網(wǎng)格進(jìn)行計(jì)算，網(wǎng)格模型如圖3所示。

圖2 不同網(wǎng)格數(shù)目靜壓分布

圖3 突擴(kuò)管網(wǎng)格模型

2.2 模擬合理性驗(yàn)證

為了驗(yàn)證模擬計(jì)算的正確性，在純水工況下選取雷諾數(shù)在104～105范圍計(jì)算上游管沿程阻力系數(shù)，并與布拉休斯經(jīng)驗(yàn)公式所得結(jié)果進(jìn)行對比。設(shè)置不同的入口速度調(diào)整雷諾數(shù)Re，其計(jì)算公式由式（6）［19］表示為

式中μs為純水動力黏度，Pa·s；v1為上游管道平均流速，m／s。

文章相關(guān)分析皆以純水工況上游直管雷諾數(shù)為依據(jù)。沿程阻力系數(shù)f由伯努利方程推導(dǎo)得出，由式（7）［13］表示為

式中Δp為兩斷面靜壓差，Pa；l為兩斷面距離，m。為保證管流充分發(fā)展且不受突擴(kuò)處影響，選取距上游管入口處10d1、15d1兩處的測壓斷面。

布拉休斯經(jīng)驗(yàn)公式是計(jì)算牛頓流體沿程阻力系數(shù)的常用公式之一，適用于紊流光滑區(qū)、4 000＜Re＜100 000范圍內(nèi)，由式（8）［19］表示為

數(shù)值模擬結(jié)果與布拉休斯經(jīng)驗(yàn)值對比如圖4所示，兩種沿程阻力系數(shù)的計(jì)算結(jié)果趨勢一致，相對誤差最大為5.6%，由此可得該模擬方法合理。

圖4 沿程阻力系數(shù)對比

2.3 Cross模型準(zhǔn)確性驗(yàn)證

流動介質(zhì)為PAM溶液時，采用Cross模型描述其動力黏度。為驗(yàn)證Cross 模型準(zhǔn)確性，建立60 mm×100 mm×1 000 mm 的槽道幾何模型，參考文獻(xiàn)［7］和［18］研究的結(jié)果，設(shè)置各參數(shù)取值λ＝5.491 s、n＝0.7、μ0＝9.835 mPa·s，描述質(zhì)量濃度為60 mg／L的PAM胺溶液黏度，模擬計(jì)算的設(shè)置同1.3 節(jié)，得到范寧系數(shù)模擬值，與實(shí)驗(yàn)所得范寧系數(shù)［7］進(jìn)行對比。范寧系數(shù)Cf由式（9）和（10）［4］表示為

式中τw為壁面平均切應(yīng)力，Pa；H為槽道截面高度，取值為60 mm；W為寬度，取值為100 mm；Δpc為槽道充分發(fā)展流段所取兩截面的靜壓差，Pa；L為兩截面的間距，取值為200 mm；U為管段平均流速，m／s。

模擬中設(shè)置入口流速為0.553 6 m／s，通過模擬結(jié)果計(jì)算得到Cf＝0.002 99，與實(shí)驗(yàn)中相同流速下的范寧系數(shù)計(jì)算結(jié)果0.002 95 對比［7］，相對誤差約為1.4%，證明Cross模型具有一定的準(zhǔn)確性。

3 結(jié)果分析

3.1 突擴(kuò)管流場分析

3.1.1 速度分布

純水與PAM溶液在Re＝10 000 工況下的速度分布如圖5所示。在模擬采用的突擴(kuò)圓管模型中，兩者的速度分布皆呈軸對稱。由于流體與管壁之間產(chǎn)生的摩擦，上游直管段的速度呈現(xiàn)從中軸線位置沿徑向向壁面減小的規(guī)律，并且上游整體流動幾乎不受突擴(kuò)的影響。經(jīng)過突擴(kuò)面，流速整體減小，在臺階處形成速度死區(qū)［20］，隨著流動的發(fā)展，中軸線處速度下降至穩(wěn)定，且速度死區(qū)逐漸消失，再次形成充分發(fā)展流；與上游管相比，下游管速度更小，這是由于管徑增大導(dǎo)致的。速度死區(qū)的出現(xiàn)是因?yàn)樵诹黧w流動過程中，截面突然擴(kuò)大，使得一部分流體從主流當(dāng)中分離出去，產(chǎn)生回流現(xiàn)象，并且速度死區(qū)通常存在著渦旋［21］，為了消除渦旋，主流要提供能量，經(jīng)突擴(kuò)面后其流速就會降低。比較二者速度云圖可以明顯看出，當(dāng)流動介質(zhì)為PAM 溶液時，與純水相比下游管近壁處的速度出現(xiàn)略微增大，表明PAM使得摩擦阻力減小，因而速度有所提高。

圖5 Re＝10 000速度云圖

3.1.2 壓力分布

取管道中軸線上的靜壓值，無量綱化為壓力系數(shù)Cp，其分布如圖6所示。無量綱壓力用中軸線各點(diǎn)靜壓與突擴(kuò)面處的靜壓差與上游管動壓之比由式（11）表示為

圖6 不同工況下管道中軸線壓力分布

式中pm、p0分別為中軸線各點(diǎn)、突擴(kuò)面處的靜壓，Pa。

根據(jù)Borda-Canot 預(yù)測值計(jì)算式［9］，即2β2（1-β2），其中β為上下游管徑比，得到壓力恢復(fù)最大值約為0.494。純水工況下，Re分別為40 000和60 000時，模擬所得壓力恢復(fù)最大值分別為0.484、0.487，可見Re的變化對此值影響不大，與理論值相比誤差最大約為2.0%，可知模擬較為準(zhǔn)確。圖6 中，上游管在0～20內(nèi)的分布線斜率表示沿程阻力系數(shù)，下游管從壓力恢復(fù)最大值處至40 的分布線斜率表示沿程阻力系數(shù)與管徑比四次方的乘積，因管徑比為定值，通過對比斜率可比較沿程阻力系數(shù)。

比較兩種雷諾數(shù)的純水工況可以看出，無論上游直管還是下游直管，Re＝40 000 時沿程阻力系數(shù)更大；Re＝60 000 時的壓力恢復(fù)最大值出現(xiàn)的位置更靠后，說明回流區(qū)長度增加。比較相同雷諾數(shù)下純水與PAM溶液的壓力分布線，可以明顯觀察到加入PAM后直管中沿程阻力系數(shù)的減小，證明減阻效果的產(chǎn)生；同時PAM 溶液壓力恢復(fù)最大值更大，且此值越大，說明回流區(qū)的渦旋損失減小［13］，因此可得PAM溶液中突擴(kuò)處的局部阻力也表現(xiàn)出減阻效果，而相較于純水工況，其壓力恢復(fù)最大值的增加并不顯著，可以推斷局部阻力的減小幅度較小。另外，相同雷諾數(shù)下PAM溶液壓力恢復(fù)最大值所處位置比純水更靠后，說明減阻劑的加入使得回流區(qū)長度增加。

3.2 直管沿程阻力減阻分析

高分子聚合物溶液可實(shí)現(xiàn)管道湍流減阻，其直管中的減阻率可高達(dá)70%［7］。減阻率RD的計(jì)算公式由式（12）［12］表示為

式中fs、f′分別為純水、高分子聚合物PAM 溶液的沿程阻力系數(shù)。

利用式（7）得到雷諾數(shù)分別為5 000、10 000、20 000、40 000、60 000、80 000、100 000 工況下的PAM溶液上游直管的沿程阻力系數(shù)，通過式（12）得到不同雷諾數(shù)下的上游直管減阻率，結(jié)果如圖7所示。

圖7 上游直管減阻率與雷諾數(shù)的關(guān)系圖

隨著雷諾數(shù)的增加上游直管段的減阻率呈現(xiàn)先增大后減小的趨勢，其減阻率最大值約為26.2%，臨界雷諾數(shù)約為20 000。出現(xiàn)這種減阻趨勢的原因是，高分子聚合物的柔性鏈?zhǔn)芗羟辛Φ睦熳饔茫沟闷渑c湍流結(jié)構(gòu)相互作用，降低湍流強(qiáng)度，對湍流渦旋產(chǎn)生抑制，產(chǎn)生了減阻效果，并在一定雷諾數(shù)范圍內(nèi)，拉伸作用加強(qiáng)而減阻效果增強(qiáng)，但是隨著雷諾數(shù)繼續(xù)增大，剪切力進(jìn)一步的增強(qiáng)使得高分子聚合物結(jié)構(gòu)無法承受其強(qiáng)度，主鏈和側(cè)鏈開始斷裂，其減阻率下降［7］。

最大減阻率數(shù)值較小的原因在于減阻效果與濃度有關(guān)［22］。此次模擬所參考的PAM溶液濃度未達(dá)到或已超出最佳濃度；另外，減阻流體的流變特性包括剪切稀化特性和黏彈性［4］，數(shù)值模擬中黏彈性的描述通過相關(guān)的黏彈性本構(gòu)方程來實(shí)現(xiàn)，而模擬所使用到的Cross 模型用以描述其剪切稀化特性，未描述黏彈性，因此導(dǎo)致模擬得到的減阻率較小，同時也說明了使用Cross模型描述減阻溶液的局限性。

3.3 突擴(kuò)局部阻力減阻分析

純水中突擴(kuò)管的局部阻力系數(shù)理論值ξ由式（13）［19］表示為

式中A1、A2分別為上、下游管的斷面面積，m2。計(jì)算可得純水工況局部阻力系數(shù)理論值ξ＝0.309。

計(jì)算純水與PAM 溶液工況局部阻力系數(shù)模擬值時，采用伯努利方程［19］推導(dǎo)，得到式（14）和（15）

式中v2為下游管的平均流速，m／s；f1、f2分別為上、下游管的沿程阻力系數(shù)，計(jì)算時上游管兩斷面取距入口10d1與15d1處，下游管兩斷面取距突擴(kuò)面10d2與15d2處；p1、p2分別為上、下游充分發(fā)展段中所取兩斷面1、2的靜壓，Pa；l10、l02分別為斷面1 至突擴(kuò)面的長度和突擴(kuò)面至斷面2的長度，m。

為綜合分析突擴(kuò)部分對管內(nèi)上、下游整體減阻效果的影響，選定兩組不同距離的斷面，對比PAM溶液局部阻力系數(shù)。方案一的兩處斷面選定于突擴(kuò)面前5d1和突擴(kuò)面后10d2的位置；方案二的斷面選定在突擴(kuò)面前10d1和突擴(kuò)面后20d2的位置。在Re為5 000至100 000的湍流流動范圍內(nèi)，局部阻力系數(shù)的結(jié)果如圖8 所示。當(dāng)流動介質(zhì)為純水時，局部阻力系數(shù)模擬值隨著雷諾數(shù)的增大先減小后逐漸趨于穩(wěn)定，其值在0.315 ～0.288間變化。PAM 溶液工況下，兩種方案所得局部阻力系數(shù)的數(shù)值整體非常接近，相差最大不超過1.0%，可知即使上、下游斷面距離不同，但整體減阻效果基本一致。局部阻力系數(shù)模擬值變化整體趨勢表現(xiàn)為隨著雷諾數(shù)的增大而增加，且未超出純水工況的值，此與PAK 等［8］實(shí)驗(yàn)趨勢一致，但增加幅度略小。雷諾數(shù)為100 000工況與純水工況下的局部阻力系數(shù)已極為接近，PAM 溶液幾乎失去減阻作用。在此次模擬的雷諾數(shù)范圍內(nèi)，突擴(kuò)局部的減阻率最大約為10.8%，與直管中減阻率最大值（26.2%）相比，突擴(kuò)局部的減阻效果不甚明顯，分析其原因在于流體經(jīng)突擴(kuò)部分產(chǎn)生了渦旋，形成了回流區(qū)，剪切力隨之增大，高分子聚合物的柔性鏈抵抗不住剪切而發(fā)生斷裂，并且隨著雷諾數(shù)增加，剪切力影響越發(fā)明顯，局部阻力系數(shù)逐漸增大，導(dǎo)致減阻效果變差。

圖8 局部阻力系數(shù)與雷諾數(shù)的關(guān)系圖

3.4 回流區(qū)流動變化

回流區(qū)渦旋作為影響局部阻力的因素，其流動、能量變化的分析對了解減阻機(jī)理、實(shí)現(xiàn)突擴(kuò)管減阻有著重要意義。將純水、PAM溶液工況下所得回流區(qū)長度模擬結(jié)果通過Xr／h進(jìn)行無量綱化表示，其中，Xr表示回流區(qū)長度，h＝（d2-d1）／2，所得結(jié)果見表1。無論是純水還是PAM減阻溶液，回流區(qū)長度都是隨著雷諾數(shù)的增加而增加；而在相同雷諾數(shù)下，PAM溶液的回流區(qū)比純水中更長。減阻溶液中回流區(qū)長度增加，而依然表現(xiàn)出減阻效果。

表1 不同雷諾數(shù)下兩種流體的回流區(qū)長度

湍動能云圖如圖9 所示。Re＝10 000 時，PAM溶液經(jīng)突擴(kuò)面后，臺階處的湍動能等值線密度較純水明顯變得稀疏，其最大等值線面積也明顯變小，說明回流區(qū)內(nèi)的湍動能有所減小，體現(xiàn)出了減阻效果。而從等值線形狀可看出，湍動能最大處在回流區(qū)末端，對比相同數(shù)值的等值線，其在PAM 溶液中比純水中在軸向上的位置靠后，說明此處湍動能軸向的減小變緩，PAM 使得湍動能沿軸向增加，因此導(dǎo)致回流區(qū)長度有所增長。

圖9 不同雷諾數(shù)下的湍動能云圖

由圖9可知，雷諾數(shù)增大，流體的湍動能也隨之增大，因而使得回流區(qū)長度增大。比較相同雷諾數(shù)下不同流體，可以觀察到PAM 溶液相較于純水，最大湍動能等值線面積皆有所減小，上游直管的壁面處湍動能也有著較為明顯的下降，說明PAM對湍動能起到了抑制作用，且通過前面的分析已知兩種雷諾數(shù)工況下都有減阻效果，可以推斷這種抑制作用是實(shí)現(xiàn)減阻的原因。Re＝100 000 的工況下，PAM溶液與純水的局部阻力系數(shù)最為接近。二者的湍動能云圖分別如圖9（c）與（d）所示，雖然最大湍動能值仍有減小，但臺階后的回流區(qū)處的等值線密度變化卻并不明顯，因此表現(xiàn)出局部阻力減阻效果的下降。

4 結(jié)論

文章利用Cross 模型描述高分子聚合物PAM溶液的黏度，模擬其在管徑比為1∶1.5 的突擴(kuò)管中不同雷諾數(shù)下的湍流流動，并于相同條件下的純水工況做對比，分析了流場、減阻效果與回流區(qū)流動變化，得出結(jié)論如下：

（1）突擴(kuò)管中的流動無論介質(zhì)為純水還是PAM溶液，經(jīng)突擴(kuò)面后都產(chǎn)生了速度死區(qū)；PAM 溶液壓力恢復(fù)最大值大于純水的，且出現(xiàn)位置靠后，說明渦旋損失有所減小且回流區(qū)長度增加。

（2）上游直管中隨著雷諾數(shù)的增大，減阻率先增大后減小，其最大減阻率約為26.2%，臨界雷諾數(shù)約為20 000。

（3）隨著雷諾數(shù)的增加，純水工況下的局部阻力系數(shù)先減小后趨于穩(wěn)定，PAM溶液的局部阻力系數(shù)呈現(xiàn)增長的趨勢，在模擬所設(shè)雷諾數(shù)范圍內(nèi)，后者數(shù)值未超過前者，PAM 溶液表現(xiàn)出減阻效果，但最大減阻率約為10.8%，較直管段減阻有較大幅度的下降。

（4）純水與PAM 溶液的回流區(qū)皆隨雷諾數(shù)的增加而伸長，相同雷諾數(shù)下PAM溶液中回流區(qū)長度更大，原因在于其末端處湍動能的軸向增加；相對于純水，PAM溶液中湍動能最大值區(qū)域減小，PAM 對湍動能產(chǎn)生了抑制作用，表現(xiàn)出減阻效果。