搜索CFD中的高階精度是不是扯淡
2017-08-12 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
(本文對英文原文適當概括)
先看高等數(shù)學對面積分的定義...,如圖:
(不要看文字...不要被各種專業(yè)名詞繞進去,其實很好理解)
他的意思就是說,對f在面S上的面積分等于n個面微元的大小乘以對應的n個f的值之和 (最后的公式)。
以上來自高等數(shù)學同濟版下冊
然后看CFD中的面積分(下面的公式和上面的一個道理):
為了準確的計算上述方程表示的面積分,需要知道每個面上的f值。當然我們不知道每個面的f值,因此我們引入了各種插值方法。這涉及到2種插值:
-
面積分的積分點可以從1到n,即k=1,2,3,...,n;
-
面上的f值需要從體心值插值而來;
先看第一條:最簡單的當然是中點(midpoint)規(guī)則。它認為積分等于f在面心的值和面積大小的乘積(即高等數(shù)學的公式中i=1,面微元就是整個面,f的值為f的面心值)。中點規(guī)則為二階精度的。
當然我們可以達到更高的精度,比如如果采用辛普森(Simpson)規(guī)則,我們的面積分這樣計算:
就是說,k=4。我們使用了4個點來計算這個面積分(中點規(guī)則使用了1個點)。
再看第二條:從體心插值到面心大家都很熟悉,就是迎風、中心、QUICK等格式。具有不同的精度。
如果我們針對第一條使用了中點規(guī)則達到了二階精度,那么我們對面心f的值進行插值的時候也要至少保留二階精度。
如果我們針對第一條使用了辛普森規(guī)則達到了四階精度,那么我們對面心f的值進行插值的時候也要至少保留四階精度(比如Cubic格式)。
不難理解,如果在從“體心插值到面心”的這個過程調用了高階精度格式,比如調用了三階精度的QUICK格式,那么我們在第一條就至少調用三階精度的面積分規(guī)則。舉例:
二階精度的中點規(guī)則+三階精度的QUICK格式=二階精度。
四階精度的辛普森規(guī)則+三階精度的QUICK格式=三階精度。
因此,只有在面積分調用高階精度的時候,插值的高階精度才有意義。
然而CFD中通常使用二階精度的中點規(guī)則,這是不是就意味著高階精度就是鬼扯...
下面我們請出Henry!看看他怎么說:
OpenFOAM使用的就是中點規(guī)則,因此不管你怎么插值,規(guī)則網(wǎng)格上的體積分都是二階精度。但是,對流項還是使用高階格式好一點。當然了你可以使用高階精度的面積分,但我覺得沒必要。或者,可以使用FDM,如果守恒不是很重要的話。
Computational Methods for Fluid dynamics 3rd Edition. Chapter 4.
相關標簽搜索：CFD中的高階精度是不是扯淡 Fluent培訓 Fluent流體培訓 Fluent軟件培訓 fluent技術教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學反應 fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析
