cfd理論過渡到編程的傻瓜入門教程-文庫吧資料
2025-07-04 08:13本頁面
【正文】 L[i]+)。i3。//Dealta+U(I) bL[i]=UI[i]Um1[i]。//Dealta+U(I+1) bR[i]=Up1[i]UI[i]。i3。 epsilon=*h。 } }}//重構(gòu)：MUSCL方法+Van Albada限制器void SolveReconstruction(double Um1[],double UI[],double Up1[],double Up2[],double UL[],double UR[]){ int i。cell[i].U[2]=。 } else if(cell[i].flag==2) { cell[i].U[0]=。cell[i].U[2]=。i++) { if(cell[i].flag==1) { cell[i].U[0]=。 //直接賦遠(yuǎn)場值 for(i=0。 }}//邊界條件：Dummy Cell方法void boundarycondition(){ int i。// printf(x=%f,d=%f,u=%f,p=%f\n,cell[i].xc,cell[i].U[0],cell[i].U[1],cell[i].U[2])。cell[i].U[2]=。 } else { cell[i].U[0]=。 cell[i].U[1]=。 x=cell[i].xc。i++) { cell[i].xc=(xi+h/2)+i*h。 for(i=0。 xi=XS12*h。 XS1=。//網(wǎng)格生成及流場初始化void initialsolve(){ int i。//entropy}。//jibenbianliang double R[3]。 double W[3],Wp[3]。//基本變量和守恒變量之間的轉(zhuǎn)換函數(shù)//前后各留兩個(gè)網(wǎng)格單元作為虛擬網(wǎng)格單元 計(jì)算網(wǎng)格單元從2到NOc3struct cell{ int flag。void SolveWtoU(double W[3],double U[3])。//時(shí)間步長double timesum。//限制器控制參數(shù)double XS1,XS2。通量計(jì)算方式: AUSM Scheme。通常，如果不是做基礎(chǔ)研究，你只需要知道它們的特點(diǎn)，會用它們就可以了，而不必深究它們怎么推導(dǎo)出來的?？傊?jì)算還是很不直接的。而、和采用所謂Roe平均值這才是Roe求解器關(guān)鍵的地方！總結(jié)一下，就是用Roe平均計(jì)算界面上的氣體狀態(tài)，然后計(jì)算得到，這樣就可以得到了。這個(gè)矩陣可以寫成特征矩陣和特征向量矩陣的形式而，和的具體表達(dá)式在許多書上都有，而且這里的矩陣表達(dá)有問題，所以就不寫了。這里和已經(jīng)用MUSCL求得了，的定義在第一講中已經(jīng)介紹了?！昂纳ⅰ边@個(gè)術(shù)語在激波捕捉格式中是最常見的。但是不要指望我把Roe求解器解釋清楚，因?yàn)檫@個(gè)不是很容易三言兩語說清的。在CFL3D軟件主頁（）上看手冊。關(guān)于這一點(diǎn)，我不打算做詳細(xì)介紹了，直接使用現(xiàn)有的近似黎曼解就可以了，都是通過和計(jì)算得到。OK，搞定了。此外還有基于特征變量的MUSCL，要復(fù)雜一點(diǎn)，但是被認(rèn)為適合更高精度的格式。密度和速度通過同樣的方法來搞定。加入斜率限制器后，上述公式就有點(diǎn)變化。因?yàn)橹苯佑枚味囗?xiàng)式去逼近一個(gè)間斷，會導(dǎo)致這樣的效果。但是這不是最終形式。所以。假設(shè)壓力p在控制體i內(nèi)部的分布是一個(gè)二次多項(xiàng)式，控制體i的中心處于處，左右兩個(gè)界面就是和。這里我只針對二次多項(xiàng)式來講解，你看完之后肯定能自己推導(dǎo)出一次多項(xiàng)式的結(jié)果（如果你搞不定，那我對你的智商表示懷疑）。MUSCL假設(shè)控制體內(nèi)原始變量（就是）的分布是一次或者二次多項(xiàng)式，如果得到了這個(gè)多項(xiàng)式，就可以求出控制體左右兩個(gè)界面的一側(cè)的值和。其兩個(gè)主要步驟就是1 怎么確定和2 怎么計(jì)算這里我們給出第一點(diǎn)一個(gè)具體的實(shí)現(xiàn)方法，就是基于原始變量的MUSCL格式（以下簡稱MUSCL）。而黎曼求解器對精度的影響不是那么大，但是對整個(gè)算法的物理適用性有影響，也就是說某種近似黎曼求解器可能對某些流動問題不合適，比如單純的Roe對于鈍頭體的脫體激波會算得亂七八糟，后來加了熵修正才算搞定。這個(gè)多項(xiàng)式的階數(shù)是和計(jì)算精度密切相關(guān)的，階數(shù)越高，誤差就越小。都是用多項(xiàng)式近似的，這是數(shù)值計(jì)算中的一個(gè)強(qiáng)大工具，你可以在很多地方看到這種近似。也不過是在這兩個(gè)問題上做文章，怎么確定，怎么計(jì)算。精確解也比較麻煩，要分析5種情況，用牛頓法迭代求解（牛頓法是什么？看數(shù)值計(jì)算的書去，哦，算了，現(xiàn)在暫時(shí)可以不必看）。什么是“精確黎曼解”，就是計(jì)算這個(gè)激波管問題的精確解。上述方法是非常重要的，是由一個(gè)蘇聯(lián)人Godunov在50年代首創(chuàng)的，后來被發(fā)展成為通用Godunov方法，著名的ENO/WENO就是其中的一種。好了，既然目的只是為了求，我在這里，只告訴你一種計(jì)算方法，也是非常重要、非常流行的一種方法。這里，初學(xué)者只需要記住這個(gè)結(jié)論：對于流動問題，不可以這
試題試卷相關(guān)推薦
鄂ICP備17016276號-1