搜索CFD玩了這么多年，聽說過著名的“快”格式不？
2017-02-24 by:CAE仿真在線 來源:互聯網
“快”格式?沒聽說過。
QUICK格式~~
在數值理論領域,如果計算網格無限大,不管使用什么格式,對傳輸方程離散后的數值解會和精確解完全一樣。然而,實際上在CFD計算中我們僅僅使用有限的、非常少量的計算網格來進行計算,在這種情況下,如果要使得解的正確,格式需要具有以下特點:
-
守恒
-
有界
-
傳輸特性
我們考慮最簡單的情況:假定網格均勻分布且間距為1.....
滿足守恒的插值格式:
如中心插分,面值為左右節點值之和的一半,如圖
在這種情況下,不管從哪個控制體來計算,面值均為左右節點值之和的一半是不變的。
不滿足守恒的插值格式:
如圖:
此圖使用了一種二次型線,即:利用3個節點的值構造一個光滑的曲線,面值即為面處曲線的高度。在這種情況下:
1. 考慮控制體2的東面,其面的值為連接1,2,3點曲線上的點A;
2. 考慮控制體3的西面,其面的值為連接1,2,3點曲線上的點B;
可見,在這種情況下,從不同控制體計算的面值竟然不同了!
我們認為這種插值格式是不守恒的。
QUICK格式:
那么同樣采用二次型線的QUICK格式為什么是守恒的呢?
非常簡單。。
先看什么是QUICK格式:Quadratric Upstream Interpolation for the Convesction Kinetics scheme,其為一個使用三點迎風二次函數來計算面值的格式,由Leonard在1979年提出。
QUICK格式的面值由“這個面左右倆個節點附加上游節點”共三個節點的值構造的二次函數型線計算。
回到那個不守恒的二次型線插值圖:
使用QUICK格式計算控制體2的東面,則A點的值為面值。B點的值扔掉。其他節點以此類推,QUICK格式最后的型線即為這樣:
這樣,上文不守恒格式的問題迎刃而解。因為,不管從哪個控制體出發,面值都是相同的。
并且,QUICK格式在均一網格為一種二階格式。
注:QUICK格式調用了相鄰節點以及另外一個節點來構造面值,因此某些的矩陣求解器如1D問題下的三對角矩陣求解器不可用。可參考:
-
CFD控制方程離散后“亂七八糟”的矩陣
(CFD界:原文使用面梯度為范例對擴散方程舉例求解來表示這種非守恒特性,相關描述略)
其他CFD界出品有關的文章:
-
守恒
-
有界、守恒和穩態
-
控制方程的“守恒”和“非守恒”
-
一張圖弄懂“向前向后差分格式”!
-
插值格式
An Introduction to Computational Fluid Dynamics. H.K. Versteeg, W. Malalasekera. Page 143
相關標簽搜索：CFD玩了這么多年，聽說過著名的“快”格式不？ Fluent培訓 Fluent流體培訓 Fluent軟件培訓 fluent技術教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學反應 fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析
