計算流體動力學從業者最頭疼的就是高質量分網問題，這通常占據整個仿真任務70%的時間。盡管在過去三十年算法的準確性和速度有了巨大進步，對復雜、真實世界模擬仿真來說，分網還是最耗時和CFD仿真過程可靠度最低的部分。無網格方法因此提供了一個可行的替代基于網格的流體計算方法，并且不需要傳統的網格結構，這樣就解決了很多分網相關的問題。下面，就來介紹當前CFD中主流的無網格方法。

什么是網格？

網格或格子定義為分析域或模型的離散單元格或單元，所有的流動變量和其他變量都在這些離散單元格中心求解。整個過程將物理域分解為更小的子域（單元/單元格）稱為分網，這些單元格分組形成邊界區域并且在這些區域施加邊界條件。不僅僅產生高質量網格并保持它是一項繁重任務，還可能有其他影響諸如：

• 收斂速率；

  
  

• 結果的準確性；

  
  

• 所需CPU時間。

盡管近年來有很多軟件有自動網格劃分功能，但大多數CFD從業者還是手動進行網格劃分。使用自動網格劃分方法用戶還是需要提供基本的輸入，如單元尺寸、需要劃分網格的區域、求解器來進行網格生成。可是這并非對所有案例都是可行的，對復雜幾何模型難于實現自動。因此，出現了CFD分析方法中的“無網格CFD”。

什么是無網格方法 ?

無網格方法用于建立整個問題域的代數方程系統，而不使用域離散的預定義網格。問題域內節點是分散的，節點在邊界上也是分散的節點組來代表（不離散）問題域及其邊界。無網格意味著無需節點間關系的信息，與傳統的有限體積或有限差分方法需要這種關系是不同的。

圖形顯示翼型劃分了網格（左側），右側相同的翼型則被點包圍。無網格方法在這些點處求解偏微分方程

為什么用無網格 ?

很多現有數值方法諸如有限體積方法，有限差分方法需要網格。在這樣的網格中，每個點有一個固定數量的預定義相鄰的點，與相鄰點的連接關系用于定義數學運算（諸如導數）并使用這個信息，求解整個域方程。

但對仿真材料可以移動或經歷大變形（例如移動網格問題），網格的連接不引入誤差則難以保持。盡管這種情況下仿真中可以重新劃分網格，最終還會導致進一步增加誤差。無網格方法可以避免這種問題。無網格方法的其他優勢：

• 節省分網的時間；

  
  

• 復雜幾何模型的仿真非常容易，而復雜幾何模型難于分網，可能需要花費幾個星期；

  
  

• 無需額外人員幫助劃分網格。

無網格方法有哪些？

近年來開發了很多無網格方法，我們來看一下發展的狀況。

光滑粒子流體動力學方法

SmoothedParticle Hydrodynamics(SPH)

SPH，用于天文學的最老的無網格方法之一，隨后越來越多用于流體流動研究。這種方法將節點視作物理顆粒，帶有質量和密度等信息，可以隨時間移動。這種方法中相鄰顆粒任何屬性值或其導數都是獨立的，顆粒可以按任何順序使用而沒有影響，顆粒可以來回移動甚至可交換位置。該方法的域離散、場函數近似和數值求解的基本步驟總結如下：

• 連續體分解為任意分布的一組沒有連接（無網格）的顆粒；

  
  

• 采用場函數近似的積分表示法；

  
  

• 引入顆粒近似轉換積分表示法為有限求和。

徑向基函數

Radial BasisFunctions(RBF)

RBF，為一種求解偏微分方程無網格方法，源于認識到徑向基函數插值對任何組節點、任意維度，都是光滑和準確的。RBF值依賴于到原點或任何其他指定點的距離的函數，來通過基于單一變量函數（徑向基函數）線性組合近似多變量函數。通常用來知道有限數量點（或太難于評估的情況）近似函數或數據。RBF一些常用類型有：

• 高斯函數；

  
  

• 多二次函數；

  
  

• 逆二次函數；

  
  

• 逆多二次函數。

有限點集方法

Finite PointsetMethod(FPM)

FPM，是使用拉格朗日法的顆粒方法，流體被有限數量的顆粒（點）代替，是非靜止顆粒。這些顆粒隨流體速度移動承載流體量，諸如密度、速度、壓力等。類似的邊界可以通過有限數量的邊界顆粒近似，邊界條件施加在上面。如同SPH方法，FPM也不使用剛性相鄰節點/顆粒（如有限體積方法FVM）列表。因此，允許所有的點/顆粒移動，每個時間步長相鄰點列表重新計算得到。這種方法適于復雜幾何模型，帶有自由表面，多相流的流動問題。

這種方法比廣泛使用的無網格方法SPH有一些優勢。SPH主要難點在于施加邊界條件，FPM方法通過使用移動最小二乘或最小二乘法非常自然地把顆粒置于邊界并施加邊界條件而解決這個困難。

其他商業無網格 CFD 軟件

讓我們了解一些商業無網格CFD軟件吧。

1

XFlow

這是商業無網格軟件之一，開發者是Next Limit Dynamics，現在屬于達索系統。它使用格子玻爾茲曼方程及基于無網格顆粒的運動求解器。XFlow能夠求解：移動邊界問題、多相流、流固耦合、瞬態分析、大渦模擬、聲學、非牛頓流體等。

風力渦輪機的CFD仿真（來源: XFlow）

XFlow仿真能夠評估渦輪效率，預測葉片載荷，尾跡湍流強度或風電場相互干擾。

2

NOGRID

NOGRID，始于2006年，用于CFD分析。無網格軟件NOGRID使用有限點集方法和N-S方程來求解CFD問題，這個軟件的特色功能包括：多相流、非牛頓流動、流固耦合、

NOGRID軟件基于網格方法，由于必須進行網格重劃分，受限的情況非常有優勢。使用快速穩健的NOGRID求解器，通常的建模和計算時間可以顯著縮短。

加速的船舶CFD仿真顯示速度動畫，使用NOGRID仿真得到的（來源：NOGRID）

（來源：NOGRID）

3

DualSPHysics

這個軟件基于光滑粒子流體動力學方法，開發用于研究歐拉方法遇到困難的自由表面流動現象，諸如水波或潰壩對海洋結構的影響。這個開源軟件產品由約翰霍普金斯大學（美國）、維戈大學（西班牙）和曼徹斯特大學（英國）研究者聯合開發，可以使用GPU來進行求解。下面案例顯示漂浮剛體相互作用問題，可以用這個軟件解決。

4

Algodoo

Algodoo是由Algoryx Simulation AB開發的，這是一個二維仿真框架使用SPH方法主要用于教學目的的軟件。易于使用，可以創建良好的教育目的和學習物理原理的可視化圖形。對學科教師和學生是一個好工具。

雖然不像高端CFD求解器能夠給出精確結果，但以非常少的時間給出可視化原理和進行很多不同仿真來說，還是一個很好的軟件。

除了流體仿真，Algodoo還支持結構、光學的仿真。

無網格技術節省了大量的網格劃分時間和精力。無網格CFD技術蓬勃發展未來前景光明，因為它解決了每個CFD工程師面對的最大困難之一—“分網”問題。

文章來源：CAE仿真學社