Sharp-interface VOF method for phase-change simulations on unstructured meshes
多面體網格的相變模擬能徹底消除傳統結構化網格中高達四重的各向異性誤差。
- 結合代數與幾何 VOF，在多面體網格實現無經驗參數的相變模擬。
- 直接以局部溫度梯度計算相變速率，徹底摒棄傳統的經驗閉包模型。
- 多面體網格不規則性能抵消系統偏差，徹底消除四重各向異性誤差。
傳統 CFD 沸騰模擬多受限於結構化網格。瑞士保羅謝爾研究所近期發表新框架，首度將代數 VOF 與幾何重建應用於任意多面體網格。研究更揭示驚人結論：標準笛卡爾網格會引發高達四重的各向異性誤差使氣泡變形，反而不規則多面體網格能完美抵消偏差，實現各向同性生長。
突破複雜幾何限制的非結構化 VOF 相變模擬
沸騰是工業界最具效率的熱傳遞模式之一，被廣泛應用於核能反應爐、化學製造、發電以及高階電子設備的冷卻系統。為了預測氣泡的成核、生長與脫離，工程師高度仰賴計算流體力學中的介面捕捉技術。在眾多方法中，VOF（Volume-of-Fluid，流體體積法）憑藉其優異的質量守恆特性成為主流，主要透過計算單一網格內的液體體積分數來追蹤液氣交界面。
過去高精度的 VOF 相變模型幾乎都建立在結構化網格上，儘管能簡化介面對流計算，但在對應真實工業設備的非典型複雜幾何時卻捉襟見肘。將技術拓展至非結構化網格面臨多重挑戰，特別是在不規則單元中計算表面法向量極易產生誤差。同時，如何在相變發生時維持介面銳利度與溫度梯度準確性，也是傳統數值方法的重大瓶頸。
結合 CICSAM 傳輸與 isoap 重建的混合策略
為兼顧計算效率與幾何精度，研究團隊採用了結合代數與幾何的混合 VOF 策略。在流體體積分數的對流傳輸部分，他們選用了 CICSAM（任意網格高壓縮介面捕捉格式）代數演算法。這套方案混合了能保持介面銳利度的高壓縮格式（Hyper-C），以及防止介面產生不自然褶皺的高階迎風格式（Ultimate Quickest），在任意網格上皆可有效抑制數值擴散。
捕捉精確介面幾何特徵時，該求解器則導入了開源的 isoap 幾何處理庫。相較於僅依賴體積分數數值的代數方法，isoap 能夠在任意多面體單元中，以線性插值精準定位體積分數等於 0.5 的交會點，並將這些點組合成一塊塊多邊形介面。這個實體重建出來的幾何多邊形，能提供極為精確的面向量與法線方向，成為後續計算跨相界熱流的堅實基礎。
拋去經驗參數：基於局部溫度梯度的熱流驅動相變
這項研究最核心的技術突破，在於徹底拋棄了傳統 CFD 軟體中依賴人工調參的經驗質量轉移模型。在處理相變時，系統嚴格遵照熱力學的 Rankine-Hugoniot 跳躍條件，直接由液相與汽相兩側朝向幾何介面的熱通量來決定體積質量轉移率。這意味著蒸發與凝結的速度完全被真實的局部溫度梯度所驅動，無需任何額外的封閉模型。
執行這項計算時，非結構化網格的邊界溫度連續性成了一大考驗。由於液氣交界面的溫度固定為飽和溫度，溫度場在此處呈現不連續狀態，團隊為此引入了「介面修改最小平方法」。一旦系統偵測到兩個相鄰網格分屬不同相態，便會自動切斷傳統的梯度插值模板，改以網格中心到幾何重建介面的縮短向量來計算梯度。這種做法有效避免了跨相插值造成的溫度梯度污染，使得熱通量計算不僅具備物理一致性，且完全適用於不規則多面體網格。
克服非正交網格擴散與兩階段表面張力平滑化技術
不規則網格天然具備網格非正交性（non-orthogonality），流場的擴散項與壓力梯度很容易因此失去穩定。為了解決這個隱患，團隊在能量與動量方程式的擴散項中引入了延遲校正技術，將通量拆分為隱式求解的正交部分與顯式修正的非正交部分。藉由採用過度鬆弛（overrelaxed）修正因子，即便在嚴重扭曲的網格邊界上，計算過程依然能維持穩健的收斂能力。
針對 VOF 模擬中最難纏的表面張力與寄生電流問題，該求解器採用了連續表面力模型。為了避免不連續的體積分數導致曲率計算暴衝，團隊設計了兩階段的網格面積權重平滑化演算法來運算曲率。巧妙的是，在計算最終的表面張力體積力時，系統仍然保留銳利、未經平滑的體積分數梯度，藉此維持與壓力梯度之間的嚴格力平衡。搭配專屬的 Rhie-Chow 插值與非定常修正，徹底排除了網格不對稱造成的壓力與速度解耦。
笛卡爾網格的四重各向異性誤差與多面體網格優勢
團隊透過一維 Stefan 蒸發問題、Sucking 凝結問題以及複雜的三維 Scriven 氣泡生長測試，對該框架進行了全面驗證。模擬數據與理論解析解展現了高度的吻合，但在深入分析三維氣泡生長時，研究人員發現了一個嚴重違背傳統認知的現象。被業界公認為最工整、最理想的笛卡爾結構化六面體網格，反而產生了巨大的系統性偏差。
深入解析指出，笛卡爾網格的絕對對齊特性與修改後的最小平方法梯度模板產生了不良交互作用，系統性地高估了介面處的溫度梯度。這不僅導致氣泡半徑在增長過程中出現持續的過衝，更引發了明顯的四重各向異性（four-fold anisotropy）現象，使得原本應該是完美球體的氣泡順著網格對角線被拉伸。相反地，當相同的物理模型運行在非結構化多面體網格上時，多邊形網格面朝向的隨機性完美抵消了模板對稱性帶來的偏差，最終實現了真正的各向同性氣泡生長與單調收斂。
湍流環狀沸騰流實測與有限體積求解器的通用潛力
在確認基礎物理的準確性與多面體網格獨特優勢後，團隊進一步測試了該框架在複雜湍流場域中的表現。透過模擬向上並流環狀沸騰流，非結構化網格系統成功捕捉了波浪調製蒸發的早期瞬態變化。其運算結果不僅展現了穩定的數值表現，更與先前的 LES（大渦模擬）研究及實驗觀測在定性上保持高度一致。
更具指標意義的是，這套結合了 CICSAM 傳輸與 isoap 多邊形提取的純熱通量相變方法，完全建立在標準的單元中心有限體積法架構之上。這意味著該技術並不依附於單一的專有軟體環境，而是能無縫整合進其他支援非結構化網格的開源 CFD 程式碼中。這項突破為未來在工業散熱、核反應爐安全分析等複雜幾何環境中，執行高保真度的兩相流模擬開闢了一條深具擴展潛力的新路徑。
結合精確幾何重建與局部溫度梯度的 VOF 方法，不僅解放了複雜工業設備的相變模擬限制，更證明了不規則多面體網格在消除結構化梯度偏差上的絕對優勢。