Loading AI tools
来自维基百科,自由的百科全书
計算流體力學(英語:Computational Fluid Dynamics,簡稱CFD)是21世紀流體力學領域的重要技術之一,使用數值方法在計算機中對流體力學的控制方程進行求解,從而可預測流場的流動。目前有多種商業CFD軟件問世,比如 FLOW-3D、FLUENT、CFD-ACE+(CFDRC)、Phoenics、CFX、Star-cd等。
目前在工程領域CFD方法已經得到廣泛的應用。美國海空軍下一代F-35戰鬥機所使用的附面層分離進氣道是CFD的成果之一。附面層分離進氣道通過特殊設計形狀的突起分離流速較慢的附面層以改善渦輪風扇發動機的進氣流場。此設計比傳統的附面層隔板方法可以減輕數百公斤重量,同時在一定速度範圍內能夠維持很好的分離效率。
CFD最基本的考慮是如何把連續流體在計算機上用離散的方式處理。一個方法是把空間區域離散化成小胞腔,以形成一個立體網格或者格點,然後應用合適的算法來解運動方程(對於不粘滯流體用歐拉方程,對於粘滯流體用納維-斯托克斯方程)。另外,這樣的一個網格可以是不規則的(例如在二維由三角形組成,在三維由四面體組成)或者是規則的;前者的特徵是每個胞腔必須單獨存儲在內存中。最後,如果問題是高度動態的並且在尺度上跨越很大的範圍,網格本身應該可以動態隨時間調整,譬如在自適應網格細化方法中。
如果選擇不使用基於網格的方法,也有一些可選的替代,比較突出的有:
對於層流情況和對於所有相關的長度尺度都可以包含在格點中的湍流的情形,直接求解納維-斯托克斯方程是可能的(通過直接數值模擬)。但一般情況下,適合於問題的尺度的範圍甚至大於今天的大型並行計算機可以建模的範圍。
在這些情況下,湍流模擬需要引入湍流模型。大渦流模擬和雷諾平均納維-斯托克斯方程表述(Reynold-Averaged-Navier-Stokes, RANS)和k-ε模型或者雷諾應力模型一起,是處理這些尺度的兩種技術。
很多實例中,其他方程和納維-斯托克斯方程要同時被求解。這些其他的方程可能包括描述種類濃度、化學反應、熱傳導等。很多高級的模擬軟體允許更複雜的情形的模擬,涉及到多相流(例如,液/氣、固/氣、液/固)或者非牛頓流體(例如血液)。
所有這些方法都遵循同樣的基本的程序。
給定離散化的穩定性通常在數值上建立,而不是像在簡單的線性問題上那樣可以解析的建立。必須特別小心才能保證離散化能夠漂亮的處理不連續的解。歐拉方程和納維-斯托克斯方程兩者都可能有衝擊波和接觸表面。
在使用中的一些離散化方法包括:
直接數值模擬(DNS)捕捉了所有相關尺度的湍流運動,所以不需要對最小的尺度建立另外的模型。該途徑極其耗費資源,在現代的機器上對於複雜的問題甚至可能無法達成,因而對於流體運動最小尺度的模型還是有必要的。
雷諾平均納維-斯托克斯方程(RANS)是湍流最古老的方法。將流體運動的各物理量表達為一個統計平均值和一個脈動值的和,代入原始的N-S方程後再取統計平均,從而得到關於統計平均物理量的控制方程——RANS方程。在動量方程中,出現了一個新的應力項,它是脈動速度的二階相關,稱為雷諾應力。這個二階張量是新的未知量,必須通過一定的模型使整個方程組得到封閉(即使得方程個數與未知量個數相等),才能求解。對於這個張量有很多不同的封閉模型,常見的如零階模型,一階模型和二階模型。目前二階模型,如k-epsilon模型和k-omiga模型,是工程中最為常見的湍流模型。
下面是常見的誤解:RANS方程不適用於有一個時變平均流的流場,因為這些方程是『時間-平均』的。事實上,統計不穩定(在非靜態)流可以同等對待。這有時被稱為URANS。沒有任何雷諾平均方法中的先天缺陷使得這個情況不能被處理,但是用於閉合(封閉)方程的湍流模型必須要在均值中的變化發生的時間相對於包含大部分能量的湍流運動的時間尺度相對較大時才正確。
大渦流模擬(Large eddy simulation, LES)是一種技術,其中更小的渦流被濾掉並用亞格點尺度模型來建模,而更大的能量的渦流則被模擬出來。該方法通常比RANS模型需要更細化的網格,但是比DNS解所需要的網格粗的多。
分離渦流模擬(Detached eddy simulation,DES)是RANS模型的一個修改,其中模型在細到足夠使用LES計算的區域切換到亞格點尺度表述。接近固體邊界的區域和渦流長度尺度小於最大格點尺寸的區域被交給RANS模式的解。隨着渦流長度尺度超出格點尺寸,區域使用LES模式。所以格點解析度不像純LES要求那樣高,因此大大減少了計算開銷。雖然DES最初為Spalart-Allmaras模型而建立(Spalart等人, 1997年),它可以用其它RANS模型實現(Strelets, 2001年),通過適當的修改RANS模型中隱含或顯式地涉及到的長度尺度。所以雖然基於Spalart-Allamaras模型的DES像帶牆模型的LES一樣作用,基於其它模型的DES(譬如兩方程模型)行為就像混合RANS-LES模型。格點的生成比簡單的RANS或LES更為複雜因為要做RANS-LES切換。DES是一種非區域途徑並會產生一個跨越解的RANS和LES區域的單一的光滑速度場。
離散化完成後要進行的方程系統的基本求解由很多數值線性代數的為人熟知的許多算法完成。可以使用靜態的迭代方法,譬如對稱高斯-賽得爾(Gauss-Seidel)或者漸進鬆弛(successive overrelaxation),或者克雷洛夫子空間法(Krylov subspaces)。在後者中,解的餘數在非線性算子的一個子空間的正交基上最小化。克雷洛夫子空間法通常和一個預條件算子(preconditioner)和一個內牛頓迭代一起使用。對於非線性問題不幸的是,正交基不能通過短循環來構造(相在一般的共軛梯度法中那樣),而必須存儲整個的序列的向量。
過去幾十年間,計算流體力學進展快速,除了應用在空氣動力學、渦輪、航海工程、航空載具、氣象學模擬、環境工程學等重要科學領域;在過去的幾年中,計算流體力學在建築設計中越來越重要。計算流體力學(和熱傳學)是以電腦計算納維爾-斯托克斯方程等能量、動量方程,以求得流動系統或熱系統的控制微分方程 [1],可透過有限元素法或有限體積法求解。這些技術廣泛用於工程設計或者分析和流體交互作用的設備,例如車輛,飛行器,泵,化學裝置和通風系統等等。該技術也用於計算機圖形學,因為動畫家不能像控制固態人物那樣簡單的處理流體,因而必須使用CFD技術結合交互工具,達成動畫影片或遊戲中的流體模擬。
有許多求解納維-斯托克斯方程的商用軟體包括下面這些(按字母序排列): AVL/FIRE、CFX、Fluent、KIVA、NUMECA、Phoenics、和STAR-CD。其它軟件包作為附件或CFD工具的補充產品使用。這些包括用於後置處理的FieldView和求解細節化學動力學的KINetics。
計算流體力學可用於分析建築物廢氣對環境的影響,預測建築物中的煙霧和火災風險、量化室內環境品質以及設計自然通風系統等方面,這對於建築性能模擬有重要幫助[2]。建築性能模擬的目的是量化與建築物的設計、建造、營運和控制等各個方面[2],而計算流體力學正好提供具有科學根據的評估與計算。諸多建築能耗模擬軟體如ESP-r[3]、ENVI-met[4]、TRNSYS[5]都有求解建築計算流體力學問題之功能。
以前,大多數與建築有關的流體力學問題,例如通風分析,風荷載,風環境等,都是使用風洞測試進行檢查的,但如今,所有這些測試都可以使用計算流體力學有效地完成。計算流體力學可以在相對較短的時間內解決所有上述問題,並且與較舊的方法(實驗方法)相比,它更經濟,更有效。[6]計算流體力學可以用於預測建築物內部和周圍的空氣流動,熱傳遞和污染物輸送。計算流體力學可以輕鬆更改不同組件的類型和位置,送風條件和系統控制時間表,從而檢查各種供暖通空調(HVAC)系統的有效性和效率。此外,計算流體力學通過建模和最佳化建築場地平面圖和室內布局,幫助制定被動的加熱/冷卻/通風策略(例如自然通風)[7]。 在全球範圍內,建築部門約占總能耗的40%[8],加強建築節能設計對於總體節能有巨大效益。
由於建築部門占總消耗量的巨大比例,因此有必要研究建築物的最佳配置以減少建築物的能源消耗。為了實現這一目標,計算流體力學可以發揮重要作用。建築性能模擬和CFD程序是重要的建築設計工具,可用於評估建築性能,包括熱舒適度,室內空氣質量機械系統效率和能耗[9]。
建築物中的計算流體力學主要用於建築通風分析。
通過調節室內空氣參數(空氣溫度,相對濕度,空氣速度和空氣中的化學物質濃度)來進行建築物的通風研究,以找到室內溫度可接受的熱舒適環境。計算流體力學可有效計算室內空氣流動相關參數,用以預測建築物的通風性能[10]。
這些空氣參數對於設計舒適的室內以及室外環境中建築物的良好整合至關重要。這是因為適當的通風系統的設計和控制策略的開發需要有關以下參數的詳細信息;
上述資訊對於建築師設計建築結構與建築外殼也是有相當大的幫助。隨著電腦技術的進步,計算流體力學技術在建築物中得到了廣泛的成功應用[11]。
最近,通風及其相關領域已成為風能工程的重要組成部分。可以使用風洞調查(實驗性地)或通過計算流體力學建模(理論上)來進行通風研究。在某些應用中,自然通風系統可能比強制通風系統更可取,因為它消除或減少了機械通風系統,從而既可提供風扇能量又可節省第一筆費用。在當今時代,由於開發了許多計算流體力學軟件和其他建築物性能模擬軟件,因此評估建築物中自然/強制通風系統的可能性變得更加容易。計算流體力學分析比實驗方法非常有用,因為可以在後處理中找到變量之間的其他相關關係。通過實驗或數值獲得的數據可通過兩種方式使用:[12]
在選擇建築場地和位置時,風荷載起着重要作用。如果兩座建築物在某個位置並排存在並且有一定的間隙,那麼當通過該間隙吹動建築物兩端的風量為時,在第一瞬間分別繞過每座建築物的流量總和必須大於或等於圍繞單個建築物的末端以壓力損失為代價。
因此,將有壓力的建築物進入間隙,這將導致建築物側面的風荷載更高。當風吹過高層建築的表面時,前面的向下流動會產生渦流(如圖1所示)。靠近地面的反方向風速可能是參考風速的140%。因此,如果在該區域中存在任何建築物,則該建築物可能會受到損壞(尤其是建築物的屋頂可能會受到嚴重損壞)。如果在建築物的建造初期考慮到風荷載的影響,就可以成功地防止對建築物的這種破壞。在建築的早期,所有這些風荷載效應都是通過風洞試驗確定的,但是今天所有這些試驗都可以通過CFD分析成功地進行。為建築物提供宜人的環境的重要性正在增加。經常要求建築師和風力工程師在建築物的形成階段和施工的規劃階段仔細檢查設計(方向,位置,位置以及周圍建築物之間的間隙)。 [12]因此,通過計算流體力學分析,可以找到適合建築物方向,位置和局部最佳化問題(局部風速,對流係數等)。
Seamless Wikipedia browsing. On steroids.
Every time you click a link to Wikipedia, Wiktionary or Wikiquote in your browser's search results, it will show the modern Wikiwand interface.
Wikiwand extension is a five stars, simple, with minimum permission required to keep your browsing private, safe and transparent.