大渦模擬

大渦模擬(Large eddy simulation, LES) 是用於計算流體動力學的湍流數學模型。它最初由約瑟夫·斯瑪格林斯基於 1963 年提出，用於模擬大氣氣流， 並由迪爾多夫(1970)推廣。 LES 目前廣泛應用於多個研究領域，包括燃燒、 聲學、 和大氣邊界層模擬。透過數值求解Navier-Stokes 方程來模擬湍流需要解決非常廣泛的時間和長度尺度，所有這些都會影響流場。這樣的解析度可以透過直接數值模擬(DNS) 來實現，但 DNS 的計算成本很高，難以模擬具有複雜幾何形狀或流動配置的實際工程系統，例如湍流噴射、泵、車輛和起落架。LES 背後的主要思想是透過 Navier-Stokes 方程的低通濾波忽略最小長度尺度來降低計算成本。這種可以被視為時間和空間平均的低通濾波，有效地從數值解中去除了小尺度資訊。然而，這些資訊並不是無關緊要的，它對流場的影響必須被建模，由此衍生小尺度可以發揮重要作用的問題的活躍研究領域，例如近壁流， 反應流、 和多相流。過濾器定義和屬性.LES 濾波器可應用於空間和時間場formula_1並執行空間濾波操作、時間濾波操作或兩者一併進行。過濾後的流場加了上劃線，定義為：formula_2其中formula_3是濾波器卷積核。這也可以寫成：formula_4過濾器核心formula_3有一個相關的截止長度尺度formula_6和截止時間尺度formula_7 .小於這些的尺度將被消掉。使用上述過濾器定義，任何流場formula_8可以分為過濾和子過濾部分，如formula_9需要注意的是，大渦模擬濾波操作不滿足雷諾運算元的性質。過濾的控制方程.LES的控制方程是透過過濾控制流場的偏微分方程得到的formula_10 .不可壓縮和可壓縮 LES 控制方程之間存在差異，這導致了新濾波操作的定義。不可壓縮流動.對於不可壓縮流動，對連續性方程和 Navier-Stokes 方程進行濾波，得到濾波後的不可壓縮連續性方程，formula_11和過濾後的 Navier-Stokes 方程，formula_12其中formula_13是過濾後的壓力場和formula_14是使用過濾速度評估的應變率張量。非線性濾波平流項formula_15是LES建模困難的主要原因。它需要未過濾的速度場，這是未知的，因此必須對其進行建模。下面的分析說明了非線性帶來的困難，即它導致大小尺度之間的相互作用，防止尺度分離。過濾後的平流項可以按照萊昂納德(1975)拆分為：formula_16其中formula_17是殘餘應力張量，因此過濾後的 Navier-Stokes 方程變為formula_18與殘餘應力張量formula_17對所有未封閉的項進行分組。 Leonard 將這個應力張量分解為formula_20併為每個項提供物理解釋。 formula_21 是倫納德張量，代表大尺度之間的相互作用；formula_22 為類雷諾應力項，表示子濾波器尺度 (SFS) 之間的相互作用；formula_23 是克拉克張量， 表示大尺度和小尺度之間的跨尺度相互作用。 建模未閉合項formula_17是亞網格尺度模型需要解析的項。亞網格應力張量formula_17的存在，使這變得具有挑戰性必須考慮所有尺度之間的相互作用，包括過濾尺度與未過濾尺度。被動標量formula_8的濾波控制方程 ，例如混合分數或溫度，可以寫成formula_27其中formula_28是formula_8的擴散通量，formula_30是標量formula_8的子濾波器通量。過濾後的擴散通量formula_32是未封閉的，除非假設它具有特定的形式，例如梯度擴散模型formula_33 . formula_30的定義形式類似formula_17 ：formula_36並且可以類似地分解為各種尺度之間相互作用的貢獻。這種子過濾器通量也需要子過濾器模型。推導.使用愛因斯坦求和約定，笛卡爾座標中不可壓縮流體的 Navier-Stokes 方程為formula_37formula_38過濾動量方程導致formula_39如果我們假設過濾和微分對等，那麼formula_40這個方程模擬了過濾變數的時間變化formula_41 .由於未過濾的變數formula_42未知，無法直接計算formula_43 。然而，formula_44是已知的。進行了替換：formula_45令formula_46 。得到的方程組是 LES 方程：formula_47可壓縮控制方程.對於可壓縮流動的控制方程，每個方程都從質量守恆開始過濾。這給出了：formula_48這導致了一個額外的子過濾器項。然而，希望避免必須對質量守恆方程的子過濾器尺度進行建模。出於這個原因，Favre 提出了一種密度加權濾波操作，稱為 Favre 濾波，定義為任意量formula_8為：formula_50在不可壓縮的限制下，它變成了正常的過濾操作。這使得質量守恆方程化為：formula_51然後可以將此概念擴充套件到編寫可壓縮流的 Favre 過濾動量方程。根據弗雷曼：formula_52其中formula_53是剪應力張量，對於牛頓流體，由下式給出：formula_54formula_55這一項表示評估粘度的子過濾器粘性貢獻formula_56使用 Favre 過濾溫度formula_57 . Favre 濾波動量場的亞網格應力張量：formula_58以此類推，萊昂納德分解也可以寫成濾波三重積的殘餘應力張量formula_59 。三重乘積可以使用 Favre 過濾運算子重寫為formula_60 ，這是一個未封閉的項。它可以以類似於上面提到的formula_15的方式分解，這導致子濾波器應力張量formula_66 。這個子過濾器項可以分為三種類型的相互作用的貢獻：倫達德張量formula_67, 代表解析尺度之間的相互作用；克拉克張量formula_68 ，表示已解析和未解析的尺度之間的相互作用；和雷諾張量formula_69 ，它表示未解析的尺度之間的相互作用。濾波動能方程.除了過濾的質量和動量方程之外，過濾動能方程可以提供額外的見解。可以過濾動能場以產生總過濾動能：formula_70過濾後的總動能可以分解為兩項：過濾後的速度場的動能formula_71 ,formula_72和剩餘動能formula_73 ,formula_74這樣formula_75 .守恆方程為formula_71可以透過將濾波後的動量傳輸方程乘以formula_77產生：formula_78其中formula_79是粘性應力對過濾速度場的動能的耗散，並且formula_80表示動能的子過濾尺度 (SFS) 耗散。左邊的項代表運輸，右邊的項是耗散動能的匯項。formula_81 SFS 耗散項特別令人感興趣，因為它代表能量從大分辨尺度到小未分辨尺度的轉移。一般， formula_81將能量從大尺度轉移到小尺度。然而，瞬間formula_81可以是積極的"或"消極的，這意味著它也可以作為源項formula_71, 濾波後的速度場的動能。能量從未解析到已解析尺度的傳遞稱為反向散射。LES的數值方法.大渦模擬涉及使用計算流體動力學求解離散濾波控制方程。 LES 從域大小中解析尺度formula_85過濾器尺寸formula_6 ，因此必須解決很大一部分高波數湍流波動。這需要高階數值方案，或者如果使用低階數值方案，則需要精細的網格解析度。 波普 的第 13 章解決了網格解析度有多精細的問題formula_87需要解析過濾的速度場formula_88 . Ghosal 發現，對於低階離散化方案，例如在有限體積方法中使用的那些，截斷誤差可以與子濾波器尺度貢獻相同，除非濾波器寬度formula_6比網格間距formula_87大得多。雖然偶數階方案具有截斷誤差，但它們是非耗散的， 並且由於子濾波器尺度模型是耗散的，偶數階方案不會像耗散方案那樣強烈地影響子濾波器尺度模型的貢獻。過濾器實現.大渦模擬中的濾波操作可以是隱式的，也可以是顯式的。隱式過濾認識到子過濾器比例模型將以與許多數值方案相同的方式消散。透過這種方式，可以假設網格或數值離散化方案是 LES 低通濾波器。雖然這充分利用了網格解析度，並消除了計算子濾波器比例模型項的計算成本，但很難確定與一些數值問題相關的 LES 濾波器的形狀。此外，截斷誤差也可能成為問題。在顯式濾波中， LES 濾波器應用於離散的 Navier-Stokes 方程，提供明確定義的濾波器形狀並減少截斷誤差。然而，顯式過濾比隱式過濾需要更精細的網格，並且計算量與formula_91成正比。Sagaut (2006) 的第 8 章更詳細地介紹了 LES 的數值問題。大渦模擬的邊界條件.入口邊界條件對LES的精度影響很大，對於LES入口條件的處理是一個複雜的問題。理論上，一個良好的 LES 邊界條件應包含以下特徵：(1) 提供準確的流動特性資訊，即速度和湍流；(2) 滿足 Navier-Stokes 方程和其他物理；(3)易於實施和適應不同情況。目前，為 LES 生成入口條件的方法大致分為 Tabor 等人分類的兩類：產生湍流入口的第一種方法是根據具體情況合成它們，例如傅立葉技術、原理正交分解和渦流方法。合成技術試圖在入口處構建具有合適的類湍流特性的湍流場，並使其易於指定湍流引數，例如湍流動能和湍流耗散率。此外，使用隨機數生成的入口條件在計算上並不昂貴。然而，該方法存在一個嚴重的缺陷。合成的湍流不滿足由 Navier-Stokes 方程控制的流體流動的物理結構。第二種方法涉及一個單獨的前體計算，以生成一個湍流資料庫，該資料庫可以引入到入口處的主要計算中。資料庫可以透過多種方式生成，例如迴圈域、預先準備好的庫和內部對映。然而，前驅體模擬產生湍流流入的方法需要很大的計算能力。研究人員檢查了各種合成和前體計算的應用，發現入口湍流越真實，LES 預測結果就越準確。模型未解析的尺度.為了討論未解析尺度的建模，首先必須對未解析尺度進行分類。它們分為兩組：已解決的子過濾器尺度(SFS) 和子網格尺度(SGS)。解析的子濾波器尺度表示波數大於截止波數的尺度formula_92 ，但其影響被過濾器抑制。僅當使用波空間中的非區域性濾波器時，才存在已解析的子濾波器尺度。這些解析的子過濾器尺度必須使用過濾器重建來建模。子網格比例是小於截止濾波器寬度formula_6的任何比例。SGS 模型的形式取決於過濾器的實現。如LES 數值方法部分所述，如果考慮隱式 LES，則不實施 SGS 模型，並且假設離散化的數值效應模擬未解決的湍流運動的物理特性。子網格比例模型.如果沒有普遍有效的湍流描述，則在構建和應用 SGS 模型時必須利用經驗資訊，並輔以伽利略不變性等基本物理約束。存在兩類 SGS 模型；第一類是功能模型，第二類是結構模型。一些模型可能被歸類為兩者。功能模型.功能模型比結構模型更簡單，只關注以物理上正確的速率耗散能量。這些基於人工渦流粘度方法，其中湍流的影響集中在湍流粘度中。該方法將亞網格尺度上的動能耗散視為類似於分子擴散。在這種情況下，偏斜部分formula_17被建模為：formula_95其中formula_96是湍流渦流粘度和formula_97是應變率張量。根據量綱分析，渦流粘度的單位必須為formula_98 。大多數渦流粘度 SGS 模型將渦流粘度建模為特徵長度尺度和特徵速度尺度的乘積。Smagorinsky-Lilly 模型.第一個成功開發的 SGS 模型是 Smagorinsky-Lilly SGS 模型，它由Smagorinsky 開發並用於 Deardorff 的第一個 LES 模擬。 它將渦流粘度建模為：formula_99其中formula_6是網格大小，formula_101是一個常數。該方法假設小尺度的能量產生和耗散處於平衡狀態——即， formula_102 .動態模型.Germano 等 使用 Smagorinsky 模型確定了許多研究，每個研究都發現了不同的 Smagorinsky 常數值formula_101針對不同的流量配置。為了試圖為 SGS 模型制定一種更通用的方法，Germano 等提出了一個動態 Smagorinsky 模型，它使用了兩個過濾器：一個網格 LES 過濾器，表示為formula_104 ，以及一個測試 LES 濾波器formula_105，用於任何湍流的場formula_106 。測試過濾器的尺度大於網格過濾器，並在 LES 表示的已經平滑的場上增加了對湍流場的額外平滑。將測試濾波器應用於 LES 方程會產生一組新的方程，它們的形式相同但 SGS 應力項formula_107被formula_108所替代。 Germano等注意到即使兩者都由於存在未解析的尺度而無法精確計算，這兩個張量之間仍然存在一個精確的關係。這種關係被稱為 Germano 恆等式：formula_109這裡formula_110可以直接計算得出，因為它只涉及過濾的速度和測試過濾的操作。同一性的意義在於，如果假設湍流是自相似的，那麼網格和測試級別的 SGS 應力具有相同的形式formula_111和formula_112 ，則 Germano 恆等式提供了一個方程，從中可以得到 Smagorinsky 係數formula_101 可能得以推出。為了做到這一點，需要在原始推導中引入x兩個額外的步驟。首先，假設即使formula_101原則上是可變的，變化足夠慢，可以從過濾操作中移出：formula_115 。其次，既然formula_101是一個標量，Germano 恆等式與一個二階張量 formula_101值 。因此對formula_101的推算條件多餘了。他提議formula_101透過最小化殘差使用最小二乘擬合來確定。由此得：formula_122其中formula_123。為簡潔起見formula_124, formula_125在 LES 模擬中實施該模型的最初嘗試被證明是不成功的。首先，計算出的係數根本不像假設的那樣“緩慢變化”，而且變化與任何其他湍流場一樣多。其次，計算formula_101可以是積極的，也可以是消極的。後一個事實本身不應被視為缺點，因為使用過濾的 DNS 欄位的先驗測試表明，本地子網格耗散率formula_127即使在流體域上的積分始終為正，表示大尺度上的能量淨耗散，在湍流場中的積分幾乎與正數一樣可能為負數。與渦粘度的嚴格正相反，正值的輕微優勢導致觀察到的淨耗散。這種從小到大的能量“反向散射”確實對應於 Smagorinsky 模型中的負 C 值。然而，發現 Germano-Lilly 公式不能產生穩定的計算。透過在同質方向上平均分子和分母，採用了一種特別的方法formula_128當平均涉及足夠大的統計樣本時，計算formula_101是正的，無論有或沒有平均，也會提高計算穩定性。 Meneveau 提出對具有指數衰減“記憶”的拉格朗日流體軌跡進行平均。這可以應用於缺乏均勻方向的問題，並且如果進行平均的有效時間足夠長，但又不足以消除感興趣的空間不均勻性，那麼這可以是穩定的。Lilly 對 Germano 方法的修改，然後是統計平均或綜合去除負粘度區域似乎是臨時的，即使它可以“工作”。 Ghosal 等人提出了一種稱為“動態定位模型”的最小二乘最小化過程的替代公式。 在這種方法中，首先定義一個量formula_130張量formula_17和formula_132替換為適當的 SGS 模型。然後，該張量表示子網格模型未能在每個空間位置符 Germano 恆等式的數量。在Lilly的方法中，formula_101可以脫離帽符運算元：formula_134使得formula_135成為一個formula_101的代數函式，然後透過要求確定formula_137被認為是formula_101的函式具有最小的可能值。然而，由於formula_101如此獲得的結果與湍流中的任何其他波動量一樣可變，原始假設formula_101不能後天證明。在 DLM 方法中，透過不呼叫從測試過濾操作中刪除formula_101的步驟來避免這種不一致。相反，人們透過數量定義了整個流的全域性誤差：formula_142其中積分範圍在整個流體體積上。這個全域性誤差formula_143是空間變化函式formula_144 的函式，formula_101的結果需要讓誤差最小化。這個變分問題的解決方法是formula_101必須滿足第二類 Fredholm 積分方程formula_148其中函式formula_149和formula_150根據解析的項formula_151定義，因此在每個時間步長和整個流體域的積分範圍都是已知的。積分方程透過迭代程式進行數值求解，如果與預處理方案一起使用，則發現收斂通常很快。儘管這種變分方法消除了 Lilly 方法中固有的不一致性，但formula_144從積分方程獲得的仍然顯示出與負粘度相關的不穩定性。這可以透過堅持來解決formula_153受約束最小化formula_154 .這導致了一個formula_101的非線性方程：formula_156這裡字尾 + 表示“正數部分”，即， formula_157 。儘管這表面上看起來像“剪裁”，但它不是一個臨時方案，而是約束變分問題的真正解決方案。這種 DLM(+) 模型被發現是穩定的，並且對於強制和衰減的各向同性湍流、通道流動和各種其他更復雜的幾何形狀產生了出色的結果。如果流動恰好具有均勻的方向，那麼可以引入假定 formula_158 。然後，變分方法立即產生 Lilly 的結果，對均勻方向進行平均，無需對先前結果進行臨時修改。DLM(+) 模型的一個缺點是它沒有描述反向散射，而反向散射在分析 DNS 資料時被認為是真實的“事物”。開發了兩種方法來解決這個問題。由於 Carati 等人的一種方法。 類似於朗道的脈動流體力學理論，增加了一個由漲落耗散定理確定的振幅的脈動力。在第二種方法中，有人注意到任何“反向散射”能量出現在解析尺度中，只是以亞網格尺度中的能量為代價。 DLM 可以透過一種簡單的方式進行修改，以考慮到這一物理事實，從而在本質上穩定的同時允許反向散射。 在DLM 的這個 k 方程版本中，DLM(k) 替換在 Smagorinsky 渦流粘度模型中的formula_159一項，以formula_160作為適當的速度尺度。確定formula_101的流程保持與“無約束”版本相同，除了張量formula_162, formula_163。為了確定 k，我們現在使用傳輸方程formula_165其中formula_166是運動粘度，formula_167是分別代表動能耗散和擴散的正係數。這些可以按照 DLM(+) 中的約束最小化的動態過程來確定。這種方法雖然比 DLM(+) 實施起來更昂貴，但被發現是穩定的，並且與測試的各種流的實驗資料有很好的一致性。此外，DLM(k) 在數學上不可能導致計算不穩定，因為大尺度和 SGS 能量的總和不會因構造而增加。這兩種結合反向散射的方法都運作良好。與 DLM(+) 相比，它們產生的模型耗散稍小，效能有所提高。 DLM(k) 模型還產生子網格動能，這可能是一個感興趣的物理量。這些改進是在模型實現的成本有所增加的情況下實現的。動態模型起源於 1990 年斯坦福大學湍流研究中心(CTR) 的暑期專案 。