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