【泵技術知識】射流泵內部構件的新型式摹擬剖析

2017-01-14 by:CAE仿真在線來源:互聯網

1控制方程

弱可壓縮流動控制方程組由可壓縮連續方程和動量方程導出。該模型能對非定常流動進行模擬且具有相當快的計算速度。

對于弱可壓縮流動,其連續方程和動量方程為

pt K u = 0(1)

ut uu p-2u = 0(2)

式中: K = a2.

在湍流運動中,各種流動的特征量均隨時間和空間坐標而呈現隨機的脈動,是三維的有渦流動而且伴隨著渦的強烈脈動。通過三維渦量場中旋渦的拉伸和變形,形成湍流中各種不同尺度的旋渦。

而這些不同尺度的旋渦在湍流運動中起著不同的作用。大尺度旋渦從時均流動中取得能量,能量由大尺度旋渦向小尺度旋渦逐級傳遞,并最后在小尺度旋渦中,通過流體的粘性將能量耗散。

大渦模擬克服了湍流模式理論的時均處理和普適性方面存在的缺陷。其基本思想是:首先,把包括脈動在內的湍流瞬時運動通過某種濾波方法分解成大尺度運動和小尺度運動兩部分,大尺度量要通過數值求解運動微分方程直接計算出來,小尺度運動對大尺度運動的影響通過建立亞格子模型來模擬(叫作次網格尺度模擬, Subgrid Scale,簡稱為SGS) ,這樣就大大減少了計算工作量和對內存的需求。

大渦模擬中對控制方程采取某種濾波處理,使得方程計算大渦時是精確的,而小渦的影響則以某種需要模擬的項出現。由于只需要模擬更小的渦,所以精度比雷諾平均方法要高。濾波的過程類似于對N S方程的時間平均,但濾波是對控制方程在空間上的處理,所以也可稱為空間平均。空間平均與時間平均在本質上是不同的,即便是統計定常的流動,濾波后的流動參數也不是常數,它仍然包含大渦的脈動。關于濾波方法目前有很多,其中高斯型濾波函數在實際應用中比較簡潔。

對于湍流流動,我們用湍流渦粘性系數t來考慮在旋渦擴散中無法求解的湍流流動的影響。計算t最簡單的模型是設定它為常數,該模型一般用于研究穩定流動。本文采用基于各向同性渦粘性假設的計算模型。

設計算域內小網格上的切應力為

ij= - 2tSij(3)

t=CS!22SijSij12(4)

Sij=12uixj ujxi(5)

式中: CS為子渦擴散系數; 為子網格的尺度; S ij為網格平均應變率。

上述模型是由Smagorinsky于1963年首先提出的。該模型將湍流粘性系數t與網格尺度及流體微團的應變率S ij聯系起來。

式中的常數系數CS目前還沒有適用于各種情況的通用值。對于二維湍流流動, Smagorinsky建議取CS= 0. 4.根據式( 4)可以分析,不同的C S值實際上反映了紊動場從時均場中提取能量的多少,顯然較大的CS值說明紊動場所提取的能量多,因而流體沿程流動的擴散就快些,反之則慢些。總之,CS值的選取需要反復試驗和比較。

!的計算公式為!= 2 Vi Vi 1 2 1 2( 6)式中: V i和V i 1為相鄰網格的體積。

大渦模擬法有其獨特的優點,但將其應用于實際三維湍流流動卻存在極大的困難。例如,通用的小渦模型需要十分密集的網格節點,進而對計算機的存儲能力要求很高;求解非線性偏微分方程要處理大量數據,需要計算機具有高速數值處理的能力等。若將三維大渦模擬簡化成實用的二維形式,就既利用了大渦模擬原有的優點,又可在通常的計算條件下(例如PC機)用于實際湍流運動的數值預測。文獻< 3>的研究成果表明大渦模擬的二維形式是可行的。

將控制方程在弱可壓縮流動計算域中過濾,表示相鄰網格的面積; C是比例常數(Smagorinsky常數) .

2數值模擬計算

2. 1生成計算網格

網格質量對CFD計算精度和計算效率有重要影響。對于復雜邊界的CFD問題,網格生成極為耗時,且極容易出錯,生成網格所需時間常常大于計算時間。

我們利用FLUENT軟件的網格劃分技術很容易地解決了本文的計算問題。具體做法是:首先在FLUENT的GAM BIT環境中,生成射流泵的幾何模型,然后利用GAMBIT的強大功能進行網格劃分,單劃分之前需要用戶輸入一些生成網格所需的相關參數。這里假定幾何模型已裝載到GAM BIT顯示區,從對話框可以看到,在劃分2D網格時,需要指定4組參數: facees, scheme, spacing和options.在scheme中我們選擇map,它使用指定的網格單元創造規則有序的結構網格;在spac ing我們用Interval Count指定邊界上分點時使用的間隔數。最后射流泵流場的計算網格(局部)如圖1所示。文中為了顯示清楚,選用較為稀疏的網格( 60 12 60 12) .顯示了射流泵的求解區域。

2. 2初始及邊界條件

當計算開始時,給定射流泵的工作水流、吸入水流和出口壓力均為同一數值,這顯然不是一個正常的運行工況。開始計算后,隨著時間的增加,工作水流壓力逐漸增大,吸入水流壓力逐漸減小,甚至可能成為負值,而出口壓力基本保持不變,這時的水流為非定常流動。經過一段時間后,工作水流壓力不再增大,并維持在某一數值。射流泵內的水流流動會逐步趨于定常,最終穩定在某一工況。這時計算將收斂并輸出結果。

具體說明如下:

(1)當t= 0時,上、下游斷面流動處于靜止狀態,令pS1= p01= p31uS= u0= u3= 0在初始時刻,射流泵的工作水流、吸入水流在進口處和出口處的壓力分別相等;各處的流速全部為0.

(2)當t 1 > t> 0時, p 0 = p t p t = p 01 p 0 C 1 C 2 t - t 0 t 1 C 3 t - t 0 2 t - 2 t - t 0 2 t p 0 x = 0,p S x = 0,u S x = 0開始計算后,射流泵工作壓力隨時間增加而逐漸增大,增幅符合上述的函數關系;射流泵工作水流和吸入水流在進口處沿x方向的壓力梯度為0,工作水流沿x方向的速度梯度為0.式中t 1, t 0為時間常數。除上述表達式外,還可以有其他定義,這取決于所假定的射流泵運行工況。

(3)當t> t 1時,u S x = 0,p S x = 0 p S = p 3 - !p 3式中: p S為吸入壓力; p 0為工作壓力; p 3為出口壓力; u S為吸入流速; u 0為工作流速; u 3為出口流速;!p 3為射流泵出口壓力差。

2. 3數值模擬結果及分析

理論上認為,射流泵屬于有限空間射流流動,其流場可以分為:流核區,射流核心速度保持不變; !基本流動區,射流邊界層迅速擴展到壁面,能量、質量交換加劇;回流區,這是一個可能存在的區域,如果射流在擴展到固壁之前,卷吸了所有的被引射流,則固壁邊界層會發生分離。在流動方向上產生回流,可以用無因次數Ct數來判定; 管流區。數值模擬的結果也如此。

顯示了射流泵在正常運行工況下數值試驗的可視化結果(流量比q= 0. 322;面積比m= 1. 59; L為無量綱化的軸方向長度; P s, P為無量綱化的吸入和工作流體壓力; u s, u為無量綱化的吸入和工作流體速度)。

從3射流泵內壓力變化曲線可以看出,在喉管入口段,由于工作水流本身能量還未較多地傳遞給吸入水流,再加上管道漸縮,故壓力是下降的。在喉管中,雖然管徑不變,但由于工作水流能量向吸入水流傳遞,管道內總動能下降而壓力是遞增的。

4射流泵內流速變化曲線表明,工作水流流速在噴嘴處達到最大值后便逐漸降低,而吸入水流流速卻逐漸增加,這進一步證明了能量傳遞的結果。 5則顯示了射流泵內流流場的速度分布。

由于受固壁的約束,射流泵在某些工況下會出現流線脫壁而形成回流漩渦區,回流區的產生是動態且不穩定的,因此數值模擬比較困難。有關文獻資料表明,在計算射流泵流場時一般采用拋物形流動假定,當流體發生回流或漩渦時就無法計算,還必須加入其他近似計算方法處理。而本文采用弱可壓縮流體控制方程和大渦模擬湍流模型較好地解決了這一問題。 6是比較典型的射流泵產生回流和旋渦的矢量圖。

從6( a)中可見,在該工況下射流泵吸入通道的最狹窄處出現了回流和漩渦。隨著運行工況的變化,回流和旋渦不斷增大,并向前后延伸,如6( b)、(c)所示,此時射流泵運行極不穩定。

大渦模擬中的亞格子模型系數的取值,是一個比較有爭議的問題。筆者對此也做了一些探討。在前面的計算中,統一取亞格子模型系數C S= 0. 4.為比較CS取值不同對計算結果的影響,現分別考察了另外2種不同的情況,即其他條件和參數不變,分別取C S = 0. 1和C S = 0. 5進行數值試驗。但從壓力沿程變化曲線、流速沿程變化曲線和流動矢量圖的比較來看,不存在明顯差異。也許正因為如此, C S的取值一直沒有定論。

一般來說,計算網格的疏密程度對計算結果會有影響。但在本文中,在一般精細網格(例如60 12 60 12)和非常致密網格(例如400 80 400 80)下,計算結果沒有本質的不同(非常致密網格的圖形在此很難顯示清楚,故從略)。計算中很少出現發散情況。可以認為,本文采用的數學模型和數值方法有良好的穩定性和收斂性,節省計算時間,適用于工程計算。

3數值模擬結果的驗證

由于本文數值模擬的是平面問題,而實際的實驗數據多為軸對稱情況,因此不能直接對比驗證,但可采用間接驗證的方法。在參考文獻< 4>中,陸宏圻教授建立了液體射流泵基本方程,該基本方程概括了國內外其他學者提出的各基本方程所表述的客觀規律,有大量實驗數據支持。該方程的簡化形式在工程中應用誤差總體上不超過5 .因此用基本性能曲線來間接驗證數值計算結果的正確性應當是可行的。

液體射流泵基本方程的函數關系為: h = f q, m, s,式中: h為壓力比; q為流量比; m為面積比;s為容重比。在射流泵內流場數值計算中, q值通過初始條件確定(即在初始條件中給定p c) ,面積比m和容重比s對特定的射流泵和工作介質是固定值,因此可算出相應的h值。不同的q值對應不同的h值,構成了相應的射流泵工況。計算出若干個工況后,就能得到一條q- h曲線,即射流泵基本性能曲線。而在q值相同的情況下,應用陸宏圻教授的理論公式,也可計算出相應的h值,并同樣能得到另一條q- h曲線。比較這2條曲線的吻合程度,即可間接驗證數值計算結果的正確性。

將采用理論公式計算和數值預報計算的2條q- h曲線共同繪制于7,從圖中可見,除極少部分外, 2條曲線基本上是重合或接近的。

4結論

湍流流動數值模擬是一項富有挑戰性的工作,本文僅僅探討了大渦模擬湍流模型在射流泵流場的應用。數值實驗的結果表明,本文采用的弱可壓縮流體控制方程和大渦模擬湍流模型特別適用于計算馬赫數較小的流體(例如水流)的非定常流動,且易于求解,計算效率高,是很好的射流泵優化設計和流動分析的輔助工具。

開放分享：優質有限元技術文章,助你自學成才