一、流體網絡的數值算例
1.1流體網絡節點壓力法的簡介
在各種復雜的流體網絡系統中,有許多地方存在流體的匯流及分流現象,稱這種地方為節點。流體網絡的基本構成是支路、節點和邊界點。
支路:流體流過的一條通道。
節點:兩條以上支路的交匯點,通常為匯流及分流的地方。
邊界點:一種特殊的節點,其壓力是流體網絡方程組計算的已知條件。
1.2流體網絡的支路特性
在如圖1所示的系統中,設流體網絡中共有m條支路,其中第j條支路的特性方程為:
其中,△Pj——支路壓差(Pa)
Wj——支路流量(Kg/s)
Rj——管路的流阻(1/Kgm)
這里Rj為管路的流阻。流阻Rj*描述了一個支路的流道特性(流道面積、沿程和局部阻力特性)和流體特性(密度、溫度)對流體流動過程的影響。在不定常流體網絡中,各支路的流阻Rj是時間的函數。
下面以某一節點流體網絡為例,說明節點壓力法數學模塊的建立過程。
建立方程時,我們作以下假設:
?。?)P0周圍的容積都集中在P0處。
?。?)采用集總參數法,認為整個容積內的密度一致。
?。?)阻力集中在P0處。
將(1-1)式寫作:
式中Ci為第I條管路的導納,由管路結構特性決定。當管路一定時,Ci為定值。為了方便計算,可以將以上公式線性化,按下式計算:
式中Bi線性導納:
采用隱式歐拉公式將微分方程轉化為差分方程:
整理得:
式中,P0,P0'分別為前后兩個時刻的節點壓力
二、流體網絡的圖形組態及數值算例
圖2所示流體網絡包含四個節點,節點壓力分別為P2、P3、P5、P9,邊界點壓力分別為P1、P4、P6、P7、P8。設流體為不可壓縮流體,即K=0,同時假設沒有微小流量流入,即w=0。在某工況下,流體網絡的實驗測定值列于表1。
基于XDPS(XinHuaDistributedProcessingSystem)分散控制系統,可將組態系統模型的建立過程演變成為許多工程模塊的簡單連結關系。依照XDPS的建模規范,建立可壓縮流體網絡節點壓力法的組態模塊NodeP8,該模塊包括8個壓力輸入、8個流阻輸入、1個微小流量之和輸入,輸出包括1個節點壓力和8個支路流量,此外還有1個壓縮系數作為模塊參數。
將該組態頁下裝到XDPS分布式處理單元DPU,執行計算,將四個節點壓力進行輸出顯示,直到穩態。通過將某一路管系的導納發生改變,來觀察其它節點壓力的變化情況。通過管路導納的變化(如下表)對P2、P3、P5、P9這四個節點壓力的變化進行分析,我們不難發現當其中的規律:當某一個管系的綜合導納發生變化時,所有的節點壓力都會產生大小不一的變化。距離管系zui近的節點壓力變化zui大。如C1變化時,與C1直接相連的節點壓力P9變化zui大,次級相連的節點壓力P3、P2變化其次,再級相連的節點壓力P5變化zui小。
那么次級相連的壓力節點的變化是不是要比再級相連的節點壓力要大呢?其實不然,由表中可以看出,當C7變化時節點壓力P3變化zui大,次級相連的節點壓力P9變化較大,再級相連的節點壓力P2變化較小,次級相連的節點壓力P5變化zui小,而不是我們一般所認為的次級相連的節點壓力(P5、P9)變化一定比再相連的壓力節點(P2)的變化大。通過對管系綜合導納的計算,我們發現P3->P9、P3->P2、P3->P5的綜合導納分別為0.73、0.6、0.59。
由此,我們可以得出以下結論:
當某一個管系的導納發生變化時,將引起其它所有節點的壓力變化,其中與管系直接相連的節點壓力變化zui大;次級跟再級節點壓力的變化大小要根據其管系的綜合導納的大小來決定。其管系的綜合導納的大的節點壓力變化大,綜合導納的小的節點壓力變化小。
三、簡化流體網絡圖
圖2中,各個節點壓力跟導納都不一樣,無法確定導納跟壓力之間具體的關系。根據圖2,取消其中的P4邊界壓力,簡化成如下圖:
為了方便測試壓力與導納之間的關系,讓網絡中的所有的8個導納數值大小一樣,通過調整邊界壓力,使得P3、P5、P8的節點壓力分別為100、1、100。這時P3、P8與P5節點壓力相差很大,看看第四個節點的壓力變化會如何?如表3。
我們對C1進行改變時,得到如下數據:
當C1由0.1變化為1時:P2=2.73P5=1.22P8=99.94P3=99.89
當C1由0.1變化為0.001時:P2=49.99P5=1.34P8=100P3=99.94
由此發現:當C1發生變化時,
1、P2->P8之間的導納比P2->P3之間的導納要大一點,節點壓力P8的變化要比節點P3的壓力變化大一點。而由于節點壓力P8與節點壓力P3的壓力值很相近,所以壓力變化相對也很小。(根據線性綜合導納計算出P2->P8之間的導納為1577.35,P2->P3之間的導納為1414.21,P2->P3之間的導納為1577.35)
2、P2-P8之間的導納與P2-P5之間的導納數值大小一樣,通過比較發現:節點壓力P5的變化要比節點壓力P8的變化要大,可以說明當管系的導納一樣時,節點壓力越小的變化率要比節點壓力大的變化率大的多。
通過對上面兩個流體網絡節點壓力變化的規律,我們得出以下結論:
當P1=P2,C1>C2,P0發生變化時,則△P1>△P2.
當C1=C2,P1>P2,P0發生變化時,則△P1<△P2.
四、串連管路導納變化規律
上面我們討論了當管系比較復雜,我們把它再簡化成一個串連管路來看節點壓力與導納的變化規律。
設P1=100P2=50C1=100C2=1
通過組態計算出來P0=100。無論P2的值大小,P0保持不變。將C1值變大,P0值仍然不變;將C1值變小,當小到C=10時,P2開始有微弱的變化。
由此得出結論:如果某一管系的導納相對相鄰的管系導納可以忽略不計的化,那么導納大的管系兩端的節點壓力數值大小一樣。就是我們通常說的"大吃小"。
注:算例均采用單位制