gambit里面網(wǎng)格質(zhì)量檢查方法及fluent計(jì)算udf的導(dǎo)入步驟.doc

第一部分：gambit中網(wǎng)格質(zhì)量的檢查方法如何檢查網(wǎng)格質(zhì)量，用什么指標(biāo)來說明網(wǎng)格好不好呢？怎么控制？一般是什么原因造成的? 一般也就是，網(wǎng)格的角度，網(wǎng)格變形的梯度等等吧 判斷網(wǎng)格質(zhì)量的方面有很多，不知你用的是什么軟件，下面總結(jié)的是針對Gambit幫助文件的簡單歸納，不同的軟件有不同的評價(jià)單元質(zhì)量的指標(biāo)，使用時(shí)最好仔細(xì)閱讀幫助文件。Area單元面積，適用于2D單元，較為基本的單元質(zhì)量特征。 Aspect Ratio長寬比，不同的網(wǎng)格單元有不同的計(jì)算方法，等于1是最好的單元，如正三角形，正四邊形，正四面體，正六面體等；一般情況下不要超過5：1. Diagonal Ratio對角線之比，僅適用于四邊形和六面體單元，默認(rèn)是大于或等于1的，該值越高，說明單元越不規(guī)則，最好等于1，也就是正四邊形或正六面體。 Edge Ratio長邊與最短邊長度之比，大于或等于1，最好等于1，解釋同上。 EquiAngle Skew通過單元夾角計(jì)算的歪斜度，在0到1之間，0為質(zhì)量最好，1為質(zhì)量最差。最好是要控制在0到0.4之間。 EquiSize Skew通過單元大小計(jì)算的歪斜度，在0到1之間，0為質(zhì)量最好，1為質(zhì)量最差。2D質(zhì)量好的單元該值最好在0.1以內(nèi)，3D單元在0.4以內(nèi)。 MidAngle Skew通過單元邊中點(diǎn)連線夾角計(jì)算的歪斜度，僅適用于四邊形和六面體單元，在0到1之間，0為質(zhì)量最好，1為質(zhì)量最差。 Size Change相鄰單元大小之比，僅適用于3D單元，最好控制在2以內(nèi)。 Stretch伸展度。通過單元的對角線長度與邊長計(jì)算出來的，僅適用于四邊形和六面體單元，在0到1之間，0為質(zhì)量最好，1為質(zhì)量最差。 Taper錐度。僅適用于四邊形和六面體單元，在0到1之間，0為質(zhì)量最好，1為質(zhì)量最差。 Volume單元體積，僅適用于3D單元，劃分網(wǎng)格時(shí)應(yīng)避免出現(xiàn)負(fù)體積。 Warpage翹曲。僅適用于四邊形和六面體單元，在0到1之間，0為質(zhì)量最好，1為質(zhì)量最差。 另外，在Fluent中的窗口鍵入：grid quality 然后回車，F(xiàn)luent能檢查網(wǎng)格的質(zhì)量，主要有以下三個(gè)指標(biāo)：1.Maxium cell squish: 如果該值等于1，表示得到了很壞的單元；2.Maxium cell skewness: 該值在0到1之間，0表示最好，1表示最壞；3.Maxium aspect-ratio: 1表示最好。第二部分：fluent中導(dǎo)入U(xiǎn)DF進(jìn)行計(jì)算的基本步驟：1. 導(dǎo)入mesh，F(xiàn)ilereadcase2. Grid-Scale 調(diào)整x， y的系數(shù)，達(dá)到需要尺寸3. 因?yàn)槭强蓧嚎s流動(dòng)，DefineModel-Solver選density based，其他保持不變4. 當(dāng)考慮溫度時(shí)，需要選中energy equation DefineModelEnergy5. Materials選擇需要的氣體 DefineMaterial，比如這里用氮?dú)猓c(diǎn)擊copy6. 選擇好氮?dú)庵?，在密度一欄設(shè)置成理想氣體。然后點(diǎn)擊Changeclose7. Define-Operating conditions里Operating pressure設(shè)為零8. 導(dǎo)入U(xiǎn)DF：Defineuser Defined-functions- compiled-，在跳出來的對話框里，點(diǎn)擊Add，選擇自己編寫的udf，如下圖。然后點(diǎn)擊Build，主窗口會(huì)出現(xiàn)以下信息：再在compiled- UDF 頁面點(diǎn)擊Load，主窗口有以下信息：表明編譯通過，已經(jīng)可以使用了9. DefineUser Defined-function hooks界面，點(diǎn)擊Edit，選擇自己編的UDF里的adjust程序，如下圖：10. DefineUser DefinedMemory界面，設(shè)定UDMI個(gè)數(shù)，根據(jù)自己的程序。11. 設(shè)置滑移邊界條件：DefineBoundary Conditions 選擇壁面，然后點(diǎn)擊Set，進(jìn)行設(shè)置： 壁面類型：Momentum 頁面，設(shè)置滑移速度，從下拉菜單中找到對應(yīng)的UDF，如下圖Thermal 頁面，設(shè)置溫度躍變，從下拉菜單中找到對應(yīng)的UDF，如下圖：12. 設(shè)置源項(xiàng)Boundary：DefineBoundary Conditions 選擇fluids 然后點(diǎn)擊Set進(jìn)入先選擇氣體種類，然后在source Terms前打勾，然后點(diǎn)擊X Momentum右邊的Edit，出現(xiàn)如下界面：把0變?yōu)?，然后從下拉菜單中找到合適的UDF另外兩個(gè)一樣處理13. 設(shè)置入口條件：DefineBoundary Conditions 選擇inlet 然后點(diǎn)擊Set進(jìn)入設(shè)置總壓和表壓，如果要設(shè)置溫度，點(diǎn)擊Thermal頁面，進(jìn)行設(shè)置。14. 設(shè)置出口條件：DefineBoundary Conditions 選擇outlet 然后點(diǎn)擊Set進(jìn)入根據(jù)出口壓強(qiáng)進(jìn)行設(shè)置15. 初始化流場: SolveInitializeinitionalize:根據(jù)經(jīng)驗(yàn)，稍微調(diào)整，速度不要太大，溫度不要太低。16. 設(shè)置殘差收斂準(zhǔn)則： SolveMonitorsResidual可以在左上角的Plot前打勾，這樣可以實(shí)時(shí)看到殘差曲線，同時(shí)設(shè)置Residual，一般設(shè)成10E-617. 設(shè)置計(jì)算時(shí)的Courant 數(shù)SolveControlsSolutionCourant數(shù)一開始小一點(diǎn)，有利于收斂，然后