CFD中的高階精度是不是扯淡

2017-08-12  by:CAE仿真在線  來(lái)源:互聯(lián)網(wǎng)

(本文對(duì)英文原文適當(dāng)概括)


先看高等數(shù)學(xué)對(duì)面積分的定義...,如圖:

CFD中的高階精度是不是扯淡fluent圖片1

(不要看文字...不要被各種專(zhuān)業(yè)名詞繞進(jìn)去,其實(shí)很好理解)

CFD中的高階精度是不是扯淡fluent圖片2

他的意思就是說(shuō),對(duì)f在面S上的面積分等于n個(gè)面微元的大小乘以對(duì)應(yīng)的n個(gè)f的值之和 (最后的公式)。


以上來(lái)自高等數(shù)學(xué)同濟(jì)版下冊(cè)



然后看CFD中的面積分(下面的公式和上面的一個(gè)道理):

CFD中的高階精度是不是扯淡fluent圖片3

為了準(zhǔn)確的計(jì)算上述方程表示的面積分,需要知道每個(gè)面上的f值。當(dāng)然我們不知道每個(gè)面的f值,因此我們引入了各種插值方法。這涉及到2種插值:


  1. 面積分的積分點(diǎn)可以從1到n,即k=1,2,3,...,n;

  2. 面上的f值需要從體心值插值而來(lái);


先看第一條:最簡(jiǎn)單的當(dāng)然是中點(diǎn)(midpoint)規(guī)則。它認(rèn)為積分等于f在面心的值和面積大小的乘積(即高等數(shù)學(xué)的公式中i=1,面微元就是整個(gè)面,f的值為f的面心值)。中點(diǎn)規(guī)則為二階精度的。


當(dāng)然我們可以達(dá)到更高的精度,比如如果采用辛普森(Simpson)規(guī)則,我們的面積分這樣計(jì)算:

CFD中的高階精度是不是扯淡fluent圖片4

就是說(shuō),k=4。我們使用了4個(gè)點(diǎn)來(lái)計(jì)算這個(gè)面積分(中點(diǎn)規(guī)則使用了1個(gè)點(diǎn))。


再看第二條:從體心插值到面心大家都很熟悉,就是迎風(fēng)、中心、QUICK等格式。具有不同的精度。


如果我們針對(duì)第一條使用了中點(diǎn)規(guī)則達(dá)到了二階精度,那么我們對(duì)面心f的值進(jìn)行插值的時(shí)候也要至少保留二階精度。


如果我們針對(duì)第一條使用了辛普森規(guī)則達(dá)到了四階精度,那么我們對(duì)面心f的值進(jìn)行插值的時(shí)候也要至少保留四階精度(比如Cubic格式)。


不難理解,如果在從“體心插值到面心”的這個(gè)過(guò)程調(diào)用了高階精度格式,比如調(diào)用了三階精度的QUICK格式,那么我們?cè)诘谝粭l就至少調(diào)用三階精度的面積分規(guī)則。舉例:


二階精度的中點(diǎn)規(guī)則+三階精度的QUICK格式=二階精度。

四階精度的辛普森規(guī)則+三階精度的QUICK格式=三階精度。


因此,只有在面積分調(diào)用高階精度的時(shí)候,插值的高階精度才有意義。


然而CFD中通常使用二階精度的中點(diǎn)規(guī)則,這是不是就意味著高階精度就是鬼扯...

CFD中的高階精度是不是扯淡fluent圖片5


下面我們請(qǐng)出Henry!看看他怎么說(shuō):


OpenFOAM使用的就是中點(diǎn)規(guī)則,因此不管你怎么插值,規(guī)則網(wǎng)格上的體積分都是二階精度。但是,對(duì)流項(xiàng)還是使用高階格式好一點(diǎn)。當(dāng)然了你可以使用高階精度的面積分,但我覺(jué)得沒(méi)必要?;蛘?可以使用FDM,如果守恒不是很重要的話。


CFD中的高階精度是不是扯淡fluent培訓(xùn)的效果圖片6

Computational Methods for Fluid dynamics 3rd Edition. Chapter 4.


開(kāi)放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:CFD中的高階精度是不是扯淡 Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做 Fluent、CFX流體分析 HFSS電磁分析 

編輯
在線報(bào)名:
  • 客服在線請(qǐng)直接聯(lián)系我們的客服,您也可以通過(guò)下面的方式進(jìn)行在線報(bào)名,我們會(huì)及時(shí)給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國(guó)服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號(hào)金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢(xún):點(diǎn)擊咨詢(xún) 點(diǎn)擊咨詢(xún)
項(xiàng)目QQ咨詢(xún):點(diǎn)擊咨詢(xún)
email:kf@1cae.com