湍流模型

湍流模型

确定湍流輸運項的代數
湍流模型,封閉方程組。通過這組方程,Reynolds方程得以封閉。它基于對湍流過程的假設,借助經驗常數或函數,建立高階湍輸運項與低階湍輸運項直至與平均流之間的某種關系。
  • 中文名:湍流模型
  • 外文名:turbulence model
  • 别名:
  • 表達式:
  • 提出者:
  • 适用領域:
  • 方法:平均N-S方程的求解
  • 類型:零方程模型

模型分類

常用的湍流模型可根據所采用的微分方程數進行分類為:零方程模型、一方程模型、兩方程模型、四方程模型、七方程模型等。

湍流的特點

湍流是不規則、多尺度、有結構的流動,一般是三維、非定常的,具有很強的擴散性和耗散性。從物理結構上看,湍流是由各種不同尺度的帶有旋轉結構的渦疊合而成的流動,這些渦的大小及旋轉軸的方向分布是随機的。大尺度的渦主要由流動的邊界條件決定,其尺寸可以與流場的大小相比拟,它主要受慣性影響而存在,是引起低頻脈動的原因;小尺度的渦主要是由粘性力決定,其尺寸可能隻有流場尺度的千分之一的量級,是引起高頻脈動的原因。大尺度的渦破裂後形成小尺度的渦,較小尺度的渦破裂後形成更小尺度的渦。在充分發展的湍流區域内,流體渦的尺寸可在相當寬的範圍内連續變化。大尺度的渦不斷地從主流獲得能量,通過渦間的相互作用,能量逐漸向小尺寸的渦傳遞。最後由于流體粘性的作用,小尺度的渦不斷消失,機械能就轉化為流體的熱能。同時由于邊界的作用、擾動及速度梯度的作用,新的渦旋又不斷産生,湍流運動得以發展和延續。

相比于一般湍流,旋轉湍流中的旋轉效應改變了近壁湍流脈動旋度,圓周方向湍流強度增強。在流體機械,由于強旋轉、大曲率和多壁面的共同影響,旋轉湍流的各向異性特性更加突出,更容易産生流動分離,在葉片表面存在更大範圍的強剪切流動,甚至是由層流到湍流的轉捩流動

湍流計算方法

無論湍流運動多麼複雜,非穩态的連續方程和Navier-Stokes方程對于湍流的瞬時運動仍然是适用的。但是,湍流所具有的強烈瞬态性和非線性使得與湍流三維時間相關的全部細節無法用解析的方法精确描述,況且湍流流動的全部細節對于工程實際來說意義不大,因為人們所關心的經常是湍流所引起的平均流場變化。這樣,就出現了對湍流進行不同簡化處理的數學計算方法。其中,最原始的方法是基于統計平均或其他平均方法建立起來的時均化模拟方法。但這種基于平均方程與湍流模型的研究方法隻适用于模拟小尺度的湍流運動,不能夠從根本上解決湍流計算問題。為了使湍流計算更能反映不同尺度的旋渦運動,研究人員後來又發展了大渦模拟、分離渦模拟與直接數值模拟等方法。總體來說,湍流的計算方法主要分為3類: 雷諾時均模拟、尺度解析模拟和直接數值模拟。其中,前2類方法可看成是非直接數值模拟方法。

雷諾時均模拟方法

雷諾時均模拟方法是指在時間域上對流場物理量進行雷諾平均化處理,然後求解所得到的時均化控制方程。比較常用的模型包括Spalart-Allmaras模型、k-ε模型、k-ω模型和雷諾應力模型等。雷諾時均模拟方法計算效率較高,解的精度也基本可以滿足工程實際需要,是流體機械領域使用最為廣泛的湍流數值模拟方法。

尺度解析模拟方法

尺度解析模拟方法是指對流場中一部分湍流進行直接求解,其餘部分通過數學模型來計算。比較常用的模型包括大渦模拟、尺度自适應模拟、分離渦模拟和嵌入式大渦模拟等。這種方法對流場計算網格要求較高,特别是近壁區的網格密度要遠大于雷諾時均法,因此所需要的計算機資源較大,但在求解瞬态性和分離性比較強的流動,特别是流體機械偏離設計工況的流動時具有優勢。

直接數值模拟方法

直接數值模拟方法(Direct numerical simulation,DNS) 是直接用瞬态Navier-Stokes方程對湍流進行計算,理論上可以得到準确的計算結果。但是,在高雷諾數的湍流中包含尺度為10 ~ 100 μm的渦,湍流脈動的頻率常大于10 k Hz,隻有在非常微小的空間網格長度和時間步長下,才能分辨出湍流中詳細的空間結構及變化劇烈的時間特性。對于這樣的計算要求,現有的計算機能力還是比較困難的,DNS目前還無法用于真正意義上的工程計算。但是,局部時均化模型為開展DNS模拟提供了一種間接方法。該模型是一種橋接模型,通過控制模型參數可以實現從雷諾時均模拟到接近DNS的數值計算,是一種有着發展潛力的計算模型。

三種方法

1、平均N-S方程的求解。

2、大渦模拟(LES)。

3、直接數值模拟(DNS)。

但是由于葉輪機械内部結構的複雜性以及計算機運算速度較慢,大渦模拟和直接數值模拟還很少用于葉輪機械内部湍流場的計算,更多的是通過求解平均N-S方程來進行數值模拟。因為平均N-S方程的不封閉性,人們引入了湍流模型來封閉方程組,所以模拟結果的好壞很大程度上取決于湍流模型的準确度。自70年代以來,湍流模型的研究發展迅速,建立了一系列的零方程、一方程、兩方程模型和二階矩模型,已經能夠十分成功的模拟邊界層和剪切層流動。但是,對于複雜的工業流動,比如航空發動機中的壓氣機動靜葉相互幹擾問題,大曲率繞流,激波與邊界層相互幹擾,流動分離,高速旋轉以及其他一些原因,常常會改變湍流的結構,使那些能夠預測簡單流動的湍流模型失效,所以完善現有湍流模型和尋找新的湍流模型在實際工作中顯得尤為重要。

模型理論

湍流模式理論或簡稱湍流模型。湍流運動物理上近乎無窮多尺度漩渦流動和數學上的強烈非線性,使得理論實驗和數值模拟都很難解決湍流問題。雖然N-S方程能夠準确地描述湍流運動的細節,但求解這樣一個複雜的方程會花費大量的精力和時間。實際上往往采用平均N-S方程來描述工程和物理學問題中遇到的湍流運動。當我們對三維非定常随機不規則的有旋湍流流動的N-S方程平均後,得到相應的平均方程,此時平均方程中增加了六個未知的雷諾應力項 ,從而形成了湍流基本方程的不封閉問題。根據湍流運動規律以尋找附加條件和關系式從而使方程封閉就促使了幾年來各種湍流模型的發展,而且在平均過程中失去了很多流動的細節信息,為了找回這些失去的流動信息,也必須引入湍流模型。雖然許多湍流模型已經取得了某些預報能力,但至今還沒有得到一個有效的統一的湍流模型。同樣,在葉輪機械内流研究中,如何找到一種更合适更準确的湍流模型也有待于進一步研究。

模型理論的思想可追溯到100多年前,為了求解雷諾應力使方程封閉,早期的處理方法是模仿粘性流體應力張量與變形率張量關聯表達式,直接将脈動特征速度與平均運動場中速度聯系起來。十九世紀後期,Boussinesq提出用渦粘性系數的方法來模拟湍流流動,通過渦粘度将雷諾應力和平均流場聯系起來,渦粘系數的數值用實驗方法确定。到二次世界大戰前,發展了一系列的所謂半經驗理論,其中包括得到廣泛應用的普朗特混合長度理論,以及G.I泰勒渦量傳遞理論和Karman相似理論。他們的基本思想都是建立在對雷諾應力的模型假設上,使雷諾平均運動方程組得以封閉。1940年,我國流體力學專家周培源教授在世界上首次推出了一般湍流的雷諾應力輸運微分方程;1951年在西德的Rotta又發展了周培源先生的工作,提出了完整的雷諾應力模型。他們的工作現在被認為是以二階封閉模型為主的現代湍流模型理論的最早奠基工作。但因為當時計算機水平的落後,方程組實際求解還不可能。70年代後期,由于計算機技術的飛速發展,周培源等人的理論重新獲得了生命力,湍流模型的研究得到迅速發展。建立的一系列的兩方程模型和二階矩模型,已經能十分成功地模拟邊界層和剪切層流動,但是對于複雜的工業流動,比如大曲率繞流,旋轉流動,透平葉栅動靜葉互相幹擾等,這些因素對湍流的影響還不清楚,這些複雜流動也構成了進入二十一世紀後學術上和應用上先進湍流模型的研究。

湍流模型可根據微分方程的個數分為零方程模型、一方程模型、二方程模型和多方程模型。這裡所說的微分方程是指除了時均N-S方程外,還要增加其他方程才能是方程封閉,增加多少個方程,則該模型就被成為多少個模型。下面分别介紹各種湍流模型的研究現狀和進展。

選擇原則

湍流模型選取的準則:流體是否可壓、建立特殊的可行的問題、精度的要求、計算機的能力、時間的限制。為了選擇最好的模型,你需要了解不同條件的适用範圍和限制。

FLUENT軟件中提供以下湍流模型:1Spalart-Allmaras 模型;2k-ε模型;3k-ω模型;4 雷諾應力模型(RSM);5 大渦模拟模型(LES)。

1 Spalart-Allmaras模型應用範圍:

Spalart-Allmaras模型是設計用于航空領域的,主要是牆壁束縛(wall-bounded)流動,而且已經顯示出很好的效果。在透平機械中的應用也愈加廣泛。在湍流模型中利用Boussinesq逼近,中心問題是怎樣計算漩渦粘度。這個模型被Spalart-Allmaras提出,用來解決因湍流動粘滞率而修改的數量方程。

模型評價:

Spalart-Allmaras模型是相對簡單的單方程模型,隻需求解湍流粘性的輸運方程,不需要求解當地剪切層厚度的長度尺度;由于沒有考慮長度尺度的變化,這對一些流動尺度變換比較大的流動問題不太适合;比如平闆射流問題,從有壁面影響流動突然變化到自由剪切流,流場尺度變化明顯等問題。

Spalart-Allmaras模型中的輸運變量在近壁處的梯度要比k-ε中的小,這使得該模型對網格粗糙帶來數值誤差不太敏感。

Spalart-Allmaras模型不能斷定它适用于所有的複雜的工程流體。例如不能依靠它去預測均勻衰退,各向同性湍流。

2 k-ε模型

① 标準的k-ε模型:

最簡單的完整湍流模型是兩個方程的模型,要解兩個變量,速度和長度尺度。在FLUENT中,标準k-ε模型自從被Launder and Spalding提出之後,就變成工程流場計算中主要的工具了。适用範圍廣、經濟、合理的精度。它是個半經驗的公式,是從實驗現象中總結出來的。

湍動能輸運方程是通過精确的方程推導得到,耗散率方程是通過物理推理,數學上模拟相似原型方程得到的。

應用範圍:該模型假設流動為完全湍流,分子粘性的影響可以忽略,此标準κ-ε模型隻适合完全湍流的流動過程模拟。

② RNG k-ε模型:

RNG k-ε模型來源于嚴格的統計技術。它和标準k-ε模型很相似,但是有以下改進:

a、RNG模型在ε方程中加了一個條件,有效的改善了精度。

b、考慮到了湍流漩渦,提高了在這方面的精度。

c、RNG理論為湍流Prandtl數提供了一個解析公式,然而标準k-ε模型使用的是用戶提供的常數。

d、标準k-ε模型是一種高雷諾數的模型,RNG理論提供了一個考慮低雷諾數流動粘性的解析公式。這些公式的作用取決于正确的對待近壁區域。

這些特點使得RNG k-ε模型比标準k-ε模型在更廣泛的流動中有更高的可信度和精度。

③ 可實現的k-ε模型:

可實現的k-ε模型是才出現的,比起标準k-ε模型來有兩個主要的不同點:

·可實現的k-ε模型為湍流粘性增加了一個公式。

·為耗散率增加了新的傳輸方程,這個方程來源于一個為層流速度波動而作的精确方程。

術語“realizable”,意味着模型要确保在雷諾壓力中要有數學約束,湍流的連續性。

應用範圍:

可實現的k-ε模型直接的好處是對于平闆和圓柱射流的發散比率的更精确的預測。而且它對于旋轉流動、強逆壓梯度的邊界層流動、流動分離和二次流有很好的表現。

可實現的k-ε模型和RNG k-ε模型都顯現出比标準k-ε模型在強流線彎曲、漩渦和旋轉有更好的表現。由于帶旋流修正的k-ε模型是新出現的模型,所以還沒有确鑿的證據表明它比RNG k-ε模型有更好的表現。但是最初的研究表明可實現的k-ε模型在所有k-ε模型中流動分離和複雜二次流有很好的作用。

該模型适合的流動類型比較廣泛,包括有旋均勻剪切流,自由流(射流和混合層),腔道流動和邊界層流動。對以上流動過程模拟結果都比标準k-ε模型的結果好,特别是可再現k-ε模型對圓口射流和平闆射流模拟中,能給出較好的射流擴張。

模型評價

可實現的k-ε模型的一個不足是在主要計算旋轉和靜态流動區域時不能提供自然的湍流粘度,這是因為可實現的k-ε模型在定義湍流粘度時考慮了平均旋度的影響。這種額外的旋轉影響已經在單一旋轉參考系中得到證實,而且表現要好于标準k-ε模型。由于這些修改,把它應用于多重參考系統中需要注意。

3 k-ω模型

① 标準的k-ω模型:

标準的k-ω模型是基于Wilcox k-ω模型,它是為考慮低雷諾數、可壓縮性和剪切流傳播而修改的。标準的k-ε模型的一個變形就是SST k-ω模型,它在FLUENT中也是可用的

應用範圍:

Wilcox k-ω模型預測了自由剪切流傳播速率,像尾流、混合流動、平闆繞流、圓柱繞流和放射狀噴射,因而可以應用于牆壁束縛流動和自由剪切流動。

② SST k-ω模型:

SST k-ω模型由Menter發展,以便使得在廣泛的領域中可以獨立于k-ε模型,使得在近壁自由流中k-ω模型有廣泛的應用範圍和精度。為了達到此目的,k-ε模型變成了k-ω公式。SST k-ω模型和标準的k-ω模型相似,但有以下改進:

·SST k-ω模型和k-ε模型的變形增長于混合功能和雙模型加在一起。混合功能是為近壁區域設計的,這個區域對标準的k-ω模型有效,還有自由表面,這對k-ε模型的變形有效。

·SST k-ω模型合并了來源于ω方程中的交叉擴散。

·湍流粘度考慮到了湍流剪應力的傳播。

·模型常量不同。

這些改進使得SST k-ω模型比标準k-ω模型在廣泛的流動領域中有更高的精度和可信度。

③ 兩個模型的對比

兩種模型有相似的形式,有方程k和ω。SST和标準模型的不同之處是:

·從邊界層内部的标準k-ω模型到邊界層外部的高雷諾數的k-e模型的逐漸轉變。

·考慮到湍流剪應力的影響修改了湍流粘性公式。

4 RSM模型

雷諾應力模型

①GLCraft

Gibson - Launder Reynolds Stress Model with Craft wall reflection terms.

②GLWR

Gibson - Launder Reynolds Stress Model with standard wall reflection terms.

③GLnoWR

Gibson - Launder Reynolds Stress Model without standard wall reflection terms.

④SSG

Speziale - Sarkar - Gatski Reynolds stress model.

相關詞條

相關搜索

其它詞條