大渦模擬 Large eddy simulation LES 是用於計算流體動力學的湍流數學模型。它最初由約瑟夫斯瑪格林斯基於 1963 年提出，用於模擬大氣氣流， 並由迪爾多夫 1970 推廣。 LES 目前廣泛應用於多個研究領域，包括燃燒、 聲學、 和大氣邊界層模擬。 透過數值求解 NavierStokes 方程來模擬湍流需要解決非常廣泛的時間和長度尺度，所有這些都會影響流場。這樣的解析度可以透過直接數值模擬 DNS 來實現，但 DNS 的計算成本很高，難以模擬具有複雜幾何形狀或流動配置的實際工程系統，例如湍流噴射、泵、車輛和起落架。 LES 背後的主要思想是透過 NavierStokes 方程的低通濾波忽略最小長度尺度來降低計算成本。這種可以被視為時間和空間平均的低通濾波，有效地從數值解中去除了小尺度資訊。然而，這些資訊並不是無關緊要的，它對流場的影響必須被建模，由此衍生小尺度可以發揮重要作用的問題的活躍研究領域，例如近壁流， 反應流、 和多相流。 過濾器定義和屬性 LES 濾波器可應用於空間和時間場 phi boldsymbol xt 並執行空間濾波操作、時間濾波操作或兩者一併進行。過濾後的流場加了上劃線，定義為： overline phi boldsymbol xt displaystyle int_inftyinfty int_inftyinfty phi boldsymbol rtau G boldsymbol xboldsymbol rt tau dtau d boldsymbol r 其中 G 是濾波器卷積核。這也可以寫成： overline phi G star phi 過濾器核心 G 有一個相關的截止長度尺度 Delta 和截止時間尺度 tau_c 小於這些的尺度將被消掉。使用上述過濾器定義，任何流場 phi 可以分為過濾和子過濾（用素數表示）部分，如 phi bar phi phiprime 需要注意的是，大渦模擬濾波操作不滿足雷諾運算元的性質。 過濾的控制方程 LES 的控制方程是透過過濾控制流場的偏微分方程得到的 rho boldsymbol uboldsymbol xt 不可壓縮和可壓縮 LES 控制方程之間存在差異，這導致了新濾波操作的定義。 不可壓縮流動 對於不可壓縮流動，對連續性方程和 NavierStokes 方程進行濾波，得到濾波後的不可壓縮連續性方程， frac partial bar u_i partial x_i 0 和過濾後的 NavierStokes 方程， frac partial bar u_i partial t frac partialpartial x_j left overline u_i u_j right frac 1rho frac partial overline p partial x_i nu frac partialpartial x_j left frac partial bar u_i partial x_j frac partial bar u_j partial x_i right frac 1rho frac partial overline p partial x_i 2 nu frac partialpartial x_j bar S_ij 其中 bar p 是過濾後的壓力場和 bar S_ij 是使用過濾速度評估的應變率張量。非線性濾波平流項 overline u_i u_j 是 LES 建模困難的主要原因。它需要未過濾的速度場，這是未知的，因此必須對其進行建模。下面的分析說明了非線性帶來的困難，即它導致大小尺度之間的相互作用，防止尺度分離。 過濾後的平流項可以按照萊昂納德 1975 拆分為： overline u_i u_j tau_ij overline u_i overline u_j 其中 tau_ij 是殘餘應力張量，因此過濾後的 NavierStokes 方程變為 frac partial bar u_i partial t frac partialpartial x_j left overline u_i overline u_j right frac 1rho frac partial overline p partial x_i 2 nu frac partialpartial x_j bar S_ij frac partial tau_ij partial x_j 與殘餘應力張量 tau_ij 對所有未封閉的項進行分組。 Leonard 將這個應力張量分解為 tau_ij L_ij C_ij R_ij 併為每個項提供物理解釋。 L_ij overline bar u_i bar u_j bar u_i bar u_j 是倫納德張量，代表大尺度之間的相互作用；R_ij overline uprime_i uprime_j 為類雷諾應力項，表示子濾波器尺度 SFS 之間的相互作用；C_ij overline bar u_i uprime_j overline bar u_j uprime_i 是克拉克張量， 表示大尺度和小尺度之間的跨尺度相互作用。 建模未閉合項 tau_ij 是亞網格尺度（SGS）模型需要解析的項。亞網格應力張量 tau_ij 的存在，使這變得具有挑戰性必須考慮所有尺度之間的相互作用，包括過濾尺度與未過濾尺度。 被動標量 phi 的濾波控制方程 ，例如混合分數或溫度，可以寫成 frac partial overline phi partial t frac partialpartial x_j left overline u_j overline phi right frac partial overline J_phi partial x_j frac partial q_jpartial x_j 其中 J_phi 是 phi 的擴散通量，q_j 是標量 phi 的子濾波器通量。過濾後的擴散通量 overline J_phi 是未封閉的，除非假設它具有特定的形式，例如梯度擴散模型 J_phi D_phi frac partial phipartial x_i q_j 的定義形式類似 tau_ij ： q_j bar phi overline u_j overline phi u_j 並且可以類似地分解為各種尺度之間相互作用的貢獻。這種子過濾器通量也需要子過濾器模型。 推導 使用愛因斯坦求和約定，笛卡爾座標中不可壓縮流體的 NavierStokes 方程為 frac partial u_ipartial x_i 0 frac partial u_ipartial t frac partial u_iu_jpartial x_j frac 1rho frac partial ppartial x_i nu frac partial2 u_ipartial x_j partial x_j 過濾動量方程導致 overline frac partial u_ipartial t overline frac partial u_iu_jpartial x_j overline frac 1rho frac partial ppartial x_i overline nu frac partial2 u_ipartial x_j partial x_j 如果我們假設過濾和微分對等，那麼 frac partial bar u_ipartial t overline frac partial u_iu_jpartial x_j frac 1rho frac partial bar ppartial x_i nu frac partial2 bar u_ipartial x_j partial x_j 這個方程模擬了過濾變數的時間變化 bar u_i 由於未過濾的變數 u_i 未知，無法直接計算 overline frac partial u_iu_jpartial x_j 。然而， frac partial bar u_ibar u_jpartial x_j 是已知的。進行了替換： frac partial bar u_ipartial t frac partial bar u_ibar u_jpartial x_j frac 1rho frac partial bar ppartial x_i nu frac partial2 bar u_ipartial x_j partial x_j left overline frac partial u_iu_jpartial x_j frac partial bar u_ibar u_jpartial x_jright 令 tau_ij overline u_i u_j bar u_i bar u_j 。得到的方程組是 LES 方程： frac partial bar u_ipartial t bar u_j frac partial bar u_ipartial x_j frac 1rho frac partial bar ppartial x_i nu frac partial2 bar u_ipartial x_j partial x_j frac partialtau_ijpartial x_j 可壓縮控制方程 對於可壓縮流動的控制方程，每個方程都從質量守恆開始過濾。這給出了： frac partial overline rhopartial t frac partial overline u_i rho partial x_i 0 這導致了一個額外的子過濾器項。然而，希望避免必須對質量守恆方程的子過濾器尺度進行建模。出於這個原因，Favre 提出了一種密度加權濾波操作，稱為 Favre 濾波，定義為任意量 phi 為： tilde phi frac overline rho phi overline rho 在不可壓縮的限制下，它變成了正常的過濾操作。這使得質量守恆方程化為： frac partial overline rhopartial t frac partial overline rho tilde u_i partial x_i 0 然後可以將此概念擴充套件到編寫可壓縮流的 Favre 過濾動量方程。根據弗雷曼： frac partial overline rho tilde u_i partial t frac partial overline rho tilde u_i tilde u_j partial x_j frac partial overline p partial x_i frac partial tilde sigma_ij partial x_j frac partial overline rho tau_ijr partial x_j frac partialpartial x_j left overline sigma_ij tilde sigma_ij right 其中 sigma_ij 是剪應力張量，對於牛頓流體，由下式給出： sigma_ij 2 mu T S_ij frac 23 mu T delta_ij S_kk frac partialpartial x_j left overline sigma_ij tilde sigma_ij right 這一項表示評估粘度的子過濾器粘性貢獻 mu T 使用 Favre 過濾溫度 tilde T Favre 濾波動量場的亞網格應力張量： tau_ijr widetilde u_i cdot u_j tilde u_i tilde u_j 以此類推，萊昂納德分解也可以寫成濾波三重積的殘餘應力張量 overline rho phi psi 。三重乘積可以使用 Favre 過濾運算子重寫為 overline rho widetilde phi psi ，這是一個未封閉的項（它需要場 phi 和 psi 的資訊，當只有 tilde phi 和 tilde psi 兩項是已知的時候）。它可以以類似於上面提到的 overline u_i u_j 的方式分解，這導致子濾波器應力張量 overline rho left widetilde phi psi tilde phi tilde psi right 。這個子過濾器項可以分為三種類型的相互作用的貢獻：倫達德張量 L_ij 代表解析尺度之間的相互作用；克拉克張量 C_ij ，表示已解析和未解析的尺度之間的相互作用；和雷諾張量 R_ij ，它表示未解析的尺度之間的相互作用。 濾波動能方程 除了過濾的質量和動量方程之外，過濾動能方程可以提供額外的見解。可以過濾動能場以產生總過濾動能： overline E frac 12 overline u_i u_i 過濾後的總動能可以分解為兩項：過濾後的速度場的動能 E_f E_f frac 12 overline u_i overline u_i 和剩餘動能 k_r k_r frac 12 overline u_i u_i frac 12 overline u_i overline u_i frac 12 tau_iir 這樣 overline E E_f k_r 守恆方程為 E_f 可以透過將濾波後的動量傳輸方程乘以 overline u_i 產生： frac partial E_fpartial t overline u_j frac partial E_fpartial x_j frac 1rho frac partial overline u_i bar p partial x_i frac partial overline u_i tau_ijrpartial x_j 2 nu frac partial overline u_i bar S_ij partial x_j epsilon_f Pi 其中 epsilon_f 2 nu bar S_ij bar S_ij 是粘性應力對過濾速度場的動能的耗散，並且 Pi tau_ijr bar S_ij 表示動能的子過濾尺度 SFS 耗散。 左邊的項代表運輸，右邊的項是耗散動能的匯項。 Pi SFS 耗散項特別令人感興趣，因為它代表能量從大分辨尺度到小未分辨尺度的轉移。一般， Pi 將能量從大尺度轉移到小尺度。然而，瞬間 Pi 可以是積極的或消極的，這意味著它也可以作為源項 E_f 濾波後的速度場的動能。能量從未解析到已解析尺度的傳遞稱為反向散射（同樣，能量從已解析到未解析尺度的傳遞稱為前向散射）。 LES 的數值方法 大渦模擬涉及使用計算流體動力學求解離散濾波控制方程。 LES 從域大小中解析尺度 L 過濾器尺寸 Delta ，因此必須解決很大一部分高波數湍流波動。這需要高階數值方案，或者如果使用低階數值方案，則需要精細的網格解析度。 波普 的第 13 章解決了網格解析度有多精細的問題 Delta x 需要解析過濾的速度場 overline uboldsymbol x Ghosal 發現，對於低階離散化方案，例如在有限體積方法中使用的那些，截斷誤差可以與子濾波器尺度貢獻相同，除非濾波器寬度 Delta 比網格間距 Delta x 大得多。雖然偶數階方案具有截斷誤差，但它們是非耗散的， 並且由於子濾波器尺度模型是耗散的，偶數階方案不會像耗散方案那樣強烈地影響子濾波器尺度模型的貢獻。 過濾器實現 大渦模擬中的濾波操作可以是隱式的，也可以是顯式的。隱式過濾認識到子過濾器比例模型將以與許多數值方案相同的方式消散。透過這種方式，可以假設網格或數值離散化方案是 LES 低通濾波器。雖然這充分利用了網格解析度，並消除了計算子濾波器比例模型項的計算成本，但很難確定與一些數值問題相關的 LES 濾波器的形狀。此外，截斷誤差也可能成為問題。 在顯式濾波中， LES 濾波器應用於離散的 NavierStokes 方程，提供明確定義的濾波器形狀並減少截斷誤差。然而，顯式過濾比隱式過濾需要更精細的網格，並且計算量與 Delta x4 成正比。Sagaut 2006 的第 8 章更詳細地介紹了 LES 的數值問題。 大渦模擬的邊界條件 入口邊界條件對 LES 的精度影響很大，對於 LES 入口條件的處理是一個複雜的問題。理論上，一個良好的 LES 邊界條件應包含以下特徵： 1 提供準確的流動特性資訊，即速度和湍流； 2 滿足 NavierStokes 方程和其他物理； 3 易於實施和適應不同情況。 目前，為 LES 生成入口條件的方法大致分為 Tabor 等人分類的兩類： 產生湍流入口的第一種方法是根據具體情況合成它們，例如傅立葉技術、原理正交分解（POD）和渦流方法。合成技術試圖在入口處構建具有合適的類湍流特性的湍流場，並使其易於指定湍流引數，例如湍流動能和湍流耗散率。此外，使用隨機數生成的入口條件在計算上並不昂貴。然而，該方法存在一個嚴重的缺陷。合成的湍流不滿足由 NavierStokes 方程控制的流體流動的物理結構。 第二種方法涉及一個單獨的前體計算，以生成一個湍流資料庫，該資料庫可以引入到入口處的主要計算中。資料庫（有時稱為 庫）可以透過多種方式生成，例如迴圈域、預先準備好的庫和內部對映。然而，前驅體模擬產生湍流流入的方法需要很大的計算能力。 研究人員檢查了各種合成和前體計算的應用，發現入口湍流越真實，LES 預測結果就越準確。 模型未解析的尺度 為了討論未解析尺度的建模，首先必須對未解析尺度進行分類。它們分為兩組：已解決的子過濾器尺度 SFS 和子網格尺度 SGS。 解析的子濾波器尺度表示波數大於截止波數的尺度 k_c ，但其影響被過濾器抑制。僅當使用波空間中的非區域性濾波器（例如盒式或高斯濾波器）時，才存在已解析的子濾波器尺度。這些解析的子過濾器尺度必須使用過濾器重建來建模。 子網格比例是小於截止濾波器寬度 Delta 的任何比例。SGS 模型的形式取決於過濾器的實現。如 LES 數值方法部分所述，如果考慮隱式 LES，則不實施 SGS 模型，並且假設離散化的數值效應模擬未解決的湍流運動的物理特性。 子網格比例模型 如果沒有普遍有效的湍流描述，則在構建和應用 SGS 模型時必須利用經驗資訊，並輔以伽利略不變性等基本物理約束。存在兩類 SGS 模型；第一類是功能模型，第二類是結構模型。一些模型可能被歸類為兩者。 功能（渦粘）模型 功能模型比結構模型更簡單，只關注以物理上正確的速率耗散能量。這些基於人工渦流粘度方法，其中湍流的影響集中在湍流粘度中。該方法將亞網格尺度上的動能耗散視為類似於分子擴散。在這種情況下，偏斜部分 tau_ij 被建模為： tau_ijr frac 13 tau_kk delta_ij 2 nu_mathrm t bar S_ij 其中 nu_mathrm t 是湍流渦流粘度和 bar S_ij frac 12 left frac partial bar u_i partial x_j frac partial bar u_jpartial x_i right 是應變率張量。 根據量綱分析，渦流粘度的單位必須為 left nu_mathrm t right frac mathrm m2mathrm s 。大多數渦流粘度 SGS 模型將渦流粘度建模為特徵長度尺度和特徵速度尺度的乘積。 SmagorinskyLilly 模型 第一個成功開發的 SGS 模型是 SmagorinskyLilly SGS 模型，它由 Smagorinsky 開發並用於 Deardorff 的第一個 LES 模擬。 它將渦流粘度建模為： nu_mathrm t C Delta2sqrt 2bar S_ijbar S_ij C Delta2 left bar S right 其中 Delta 是網格大小，C 是一個常數。 該方法假設小尺度的能量產生和耗散處於平衡狀態 即， epsilon Pi 動態模型（Germano 等） Germano 等 使用 Smagorinsky 模型確定了許多研究，每個研究都發現了不同的 Smagorinsky 常數值 C 針對不同的流量配置。為了試圖為 SGS 模型制定一種更通用的方法，Germano 等提出了一個動態 Smagorinsky 模型，它使用了兩個過濾器：一個網格 LES 過濾器，表示為 overline f ，以及一個測試 LES 濾波器 hat f，用於任何湍流的場 f 。測試過濾器的尺度大於網格過濾器，並在 LES 表示的已經平滑的場上增加了對湍流場的額外平滑。將測試濾波器應用於 LES 方程（透過將 網格 濾波器應用於 NavierStokes 方程獲得）會產生一組新的方程，它們的形式相同但 SGS 應力項 tau_ij overline u_i u_j bar u_i bar u_j 被 T_ij widehat overline u_i u_j hat bar u_i hat bar u_j 所替代。 Germano 等注意到即使兩者都由於存在未解析的尺度而無法精確計算，這兩個張量之間仍然存在一個精確的關係。這種關係被稱為 Germano 恆等式： L_ij T_ij hat tau_ij 這裡 L_ij widehat bar u_i bar u_j widehat bar u_i widehat bar u_j 可以直接計算得出，因為它只涉及過濾的速度和測試過濾的操作。同一性的意義在於，如果假設湍流是自相似的，那麼網格和測試級別的 SGS 應力具有相同的形式 tau_ij tau_kk3delta_ij 2 C Delta2 bar S_ij bar S_ij 和 T_ij T_kk3delta_ij 2 C hat Delta2 hat bar S_ij hat bar S_ij ，則 Germano 恆等式提供了一個方程，從中可以得到 Smagorinsky 係數 C （不再是 常數）可能得以推出。 為了做到這一點，需要在原始推導中引入 x 兩個額外的步驟。首先，假設即使 C 原則上是可變的，變化足夠慢，可以從過濾操作中移出： widehat C C widehat 。其次，既然 C 是一個標量，Germano 恆等式與一個二階張量（選擇應變張量的速率）聯絡起來，以將其轉換為一個標量方程，從而推出 C。 Lilly 發現了一種從張量恆等式中推出 C 的更好的方法。他指出，Germano 恆等式需要在空間中的每個點上滿足單個量的九個方程（其中只有五個是獨立的） C 值 。因此對 C 的推算條件多餘了。他提議 C 透過最小化殘差使用最小二乘擬合來確定。由此得： C frac L_ij m_ij m_kl m_kl 其中 m_ij alpha_ij widehat beta_ij 。為簡潔起見 alpha_ij 2 hat Delta2 hat bar S hat bar S_ij beta_ij 2 Delta2 bar S bar S_ij 在 LES 模擬中實施該模型的最初嘗試被證明是不成功的。首先，計算出的係數根本不像假設的那樣 緩慢變化，而且變化與任何其他湍流場一樣多。其次，計算 C 可以是積極的，也可以是消極的。後一個事實本身不應被視為缺點，因為使用過濾的 DNS 欄位的先驗測試表明，本地子網格耗散率 tau_ij bar S_ij 即使在流體域上的積分始終為正，表示大尺度上的能量淨耗散，在湍流場中的積分幾乎與正數一樣可能為負數。與渦粘度的嚴格正相反，正值的輕微優勢導致觀察到的淨耗散。這種從小到大的能量 反向散射 確實對應於 Smagorinsky 模型中的負 C 值。然而，發現 GermanoLilly 公式不能產生穩定的計算。透過在同質方向上平均分子和分母（在流動中存在這樣的方向），採用了一種特別的方法 C frac leftlangle L_ij m_ij rightrangle leftlangle m_kl m_kl rightrangle 當平均涉及足夠大的統計樣本時，計算 C 是正的（或至少極少負值）穩定的計算是可能的。簡單地將負值設定為零（一個稱為 剪裁 的過程），無論有或沒有平均，也會提高計算穩定性。 Meneveau 提出對具有指數衰減 記憶 的拉格朗日流體軌跡進行平均。這可以應用於缺乏均勻方向的問題，並且如果進行平均的有效時間足夠長，但又不足以消除感興趣的空間不均勻性，那麼這可以是穩定的。 Lilly 對 Germano 方法的修改，然後是統計平均或綜合去除負粘度區域似乎是臨時的，即使它可以 工作。 Ghosal 等人提出了一種稱為 動態定位模型（DLM）的最小二乘最小化過程的替代公式。 在這種方法中，首先定義一個量 E_ij L_ij T_ij hat tau_ij 張量 tau_ij 和 T_ij 替換為適當的 SGS 模型。然後，該張量表示子網格模型未能在每個空間位置符 Germano 恆等式的數量。在 Lilly 的方法中，C 可以脫離帽符運算元： widehat C C widehat 使得 E_ij 成為一個 C 的代數函式，然後透過要求確定 E_ij E_ij 被認為是 C 的函式具有最小的可能值。然而，由於 C 如此獲得的結果與湍流中的任何其他波動量一樣可變，原始假設 C 不能後天證明。在 DLM 方法中，透過不呼叫從測試過濾操作中刪除 C 的步驟來避免這種不一致。相反，人們透過數量定義了整個流的全域性誤差： E C int E_ij E_ij dV 其中積分範圍在整個流體體積上。這個全域性誤差 E C xyzt 是空間變化函式 C xyzt 的函式（這裡的瞬時 t 是固定的，只是作為一個引數出現），C 的結果需要讓誤差最小化。這個變分問題的解決方法是 C 必須滿足第二類 Fredholm 積分方程 C boldsymbol x f boldsymbol x int K boldsymbol x boldsymbol y C boldsymbol y dboldsymbol y 其中函式 K boldsymbol x boldsymbol y 和 f boldsymbol x 根據解析的項 L_ijalpha_ijbeta_ij 定義，因此在每個時間步長和整個流體域的積分範圍都是已知的。積分方程透過迭代程式進行數值求解，如果與預處理方案一起使用，則發現收斂通常很快。儘管這種變分方法消除了 Lilly 方法中固有的不一致性，但 C xyzt 從積分方程獲得的仍然顯示出與負粘度相關的不穩定性。這可以透過堅持來解決 E C 受約束最小化 C xyzt geq 0 這導致了一個 C 的非線性方程： C boldsymbol x left f boldsymbol x int K boldsymbol x boldsymbol y C boldsymbol y dboldsymbol y right_ 這裡字尾 表示 正數部分，即， x_ x x2 。儘管這表面上看起來像 剪裁，但它不是一個臨時方案，而是約束變分問題的真正解決方案。這種 DLM 模型被發現是穩定的，並且對於強制和衰減的各向同性湍流、通道流動和各種其他更復雜的幾何形狀產生了出色的結果。如果流動恰好具有均勻的方向（讓我們說方向 x 和 z），那麼可以引入假定 C C yt 。然後，變分方法立即產生 Lilly 的結果，對均勻方向進行平均，無需對先前結果進行臨時修改。 DLM 模型的一個缺點是它沒有描述反向散射，而反向散射在分析 DNS 資料時被認為是真實的 事物。開發了兩種方法來解決這個問題。由於 Carati 等人的一種方法。 類似於朗道的脈動流體力學理論，增加了一個由漲落耗散定理確定的振幅的脈動力。在第二種方法中，有人注意到任何 反向散射 能量出現在解析尺度中，只是以亞網格尺度中的能量為代價。 DLM 可以透過一種簡單的方式進行修改，以考慮到這一物理事實，從而在本質上穩定的同時允許反向散射。 在 DLM 的這個 k 方程版本中，DLM k 替換在 Smagorinsky 渦流粘度模型中的 Delta bar S 一項，以 sqrt k 作為適當的速度尺度。確定 C 的流程保持與 無約束 版本相同，除了張量 alpha_ij 2 hat Delta sqrt K hat bar S_ij beta_ij 2 hat Delta sqrt k bar S_ij （其中子測試尺度動能 K 與子網格尺度動能 k 的關係為 K k L_ii2 ）（接著追蹤 Germano 恆等式）。為了確定 k，我們現在使用傳輸方程 frac partial kpartial t u_j frac partial kpartial x_j tau_ij bar S_ij frac C_Delta k32 frac partialpartial x_j left D Delta sqrt k frac partial kpartial x_j right nu frac partial2 k partial x_j partial x_j 其中 nu 是運動粘度，C_D 是分別代表動能耗散和擴散的正係數。這些可以按照 DLM 中的約束最小化的動態過程來確定。這種方法雖然比 DLM 實施起來更昂貴，但被發現是穩定的，並且與測試的各種流的實驗資料有很好的一致性。此外，DLM k 在數學上不可能導致計算不穩定，因為大尺度和 SGS 能量的總和不會因構造而增加。這兩種結合反向散射的方法都運作良好。與 DLM 相比，它們產生的模型耗散稍小，效能有所提高。 DLM k 模型還產生子網格動能，這可能是一個感興趣的物理量。這些改進是在模型實現的成本有所增加的情況下實現的。 動態模型起源於 1990 年斯坦福大學湍流研究中心 CTR 的暑期專案 。一系列 CTRTea 研討會慶祝了湍流建模這一重要里程碑的 30 週年 。 另見 流體力學 伽利略不變性 某些型別濾波器的重要屬性 雷諾平均 NavierStokes 方程 湍流 延伸閱讀 Heus T van Heerwaarden C CJonker H J J Pier Siebesma A Axelsen S DOI 105194gmd34152010 ISSN 19919603 參考資料 Category 計算流體力學 Category 亂流 Category 流體力學 Category 流體動力學 Category 偏微分方程