



全文預(yù)覽已結(jié)束
下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認(rèn)領(lǐng)
文檔簡介
1 前言工程實際中遇到的各類流動問題,當(dāng)工質(zhì)熱物理性質(zhì)偏離理想氣體時,需要應(yīng)用實際氣體性質(zhì)進行求解,才能比較精確地反映流動特性。例如在蒸汽透平低壓部分,低壓水蒸汽的性質(zhì)與理想氣體性質(zhì)偏離較大,如采用理想氣體性質(zhì)進行求解,誤差較大。IAPWS-IF97標(biāo)準(zhǔn)中給出了計算水蒸汽氣體性質(zhì)的計算公式,但由于計算量很大,并不適合在大規(guī)模三維計算中直接采用。因此有必要發(fā)展一種高效高精度的實際氣體性質(zhì)計算方法。本文以水蒸汽實際氣體的求解為例,利用FLUENT 軟件所提供的實際氣體模型(UDRGM )實現(xiàn)了大規(guī)模復(fù)雜計算中水蒸汽實際氣體性質(zhì)的應(yīng)用。2 FLUENT 中的實際氣體模型(UDRGM )及實現(xiàn)FLUENT 軟件中提供了實際氣體性質(zhì)計算的接口UDRGM (User Defined Real Gas Model)。利用該接口,可以將水蒸汽性質(zhì)表導(dǎo)入FLUENT 中,實現(xiàn)對水蒸汽實際氣體性質(zhì)的計算。UDRGM 要求對實際氣體提供如表1 所示的12 個函數(shù)。對水蒸汽定義這12 個函數(shù)后,即可實現(xiàn)在FLUENT 中調(diào)用水蒸汽實際氣體性質(zhì)進行計算。表1 UDRGM 所定義的函數(shù)2.1 實現(xiàn)方法首先由IAPWS-IF97 計算公式生成表1 中12 個函數(shù)的數(shù)據(jù)表,在啟動FLUENT 時載入該數(shù)據(jù)表,在FLUENT 中進行水蒸汽實際氣體性質(zhì)計算時,直接利用該數(shù)據(jù)表進行插值得到所需的函數(shù)值。由于插值輸入變量的數(shù)目為2,因而插值可采用雙線性插值或雙三次插值。如圖1 所示,雙線性插值根據(jù)插值點(x0,y0)周圍4 個節(jié)點上的函數(shù)值進行插值,而雙三次插值根據(jù)插值點(x0,y0)周圍16 個節(jié)點上的函數(shù)值進行插值。因此雙線性插值的精度低于雙三次插值的精度,但節(jié)省計算時間。按照文獻(xiàn)3 的估計,采用雙線性插值比采用理想氣體狀態(tài)方程進行計算的時間增加13左右,采用雙三次插值比采用理想氣體狀態(tài)方程進行計算的時間增加19左右。這種方法省去了大量計算水蒸汽實際氣體性質(zhì)的時間,因而適用于大規(guī)模的三維CFD 計算分析。所生成的水蒸汽性質(zhì)表輸入變量等間距分布,如圖1 所示,可以得到插值公式為:圖1 插值方法io = int(xo-xmin)/(xmax-xmin)+1 (1)jo = int(yo-ymin)/(ymax-ymin)+1 (2)可以看到,插值所需要的時間與數(shù)據(jù)表中數(shù)據(jù)的密度無關(guān)。因此,只要計算機內(nèi)存條件允許,可以在數(shù)據(jù)表中生成足夠密度的數(shù)據(jù),從而提高插值的精度。2.2 誤差分析為了檢驗插值方法的精度,對插值誤差進行了分析,以壓力p 和溫度T 為輸入變量, 對表1 中函數(shù)1 的插值精度進行了檢驗。檢驗中插值函數(shù)用了雙線性插值。結(jié)果如圖2 和圖3 所示,可以看到與直接由IF-97 公式計算值相比,整個區(qū)域內(nèi)插值計算的最大誤差小于1,完全滿足工程計算的要求。圖2 由水蒸汽表插值和IF-97公式計算得到的密度 圖3 插值對數(shù)誤差2.3 算例分析為了對比采用水蒸汽實際氣體性質(zhì)和理想氣體性質(zhì)計算的差異,對某Laval 噴管中的流動進行了計算。計算采用H 型網(wǎng)格,網(wǎng)格數(shù)為19060。計算中分別采用了理想氣體狀態(tài)方程、FLUENT 內(nèi)置的維里型狀態(tài)方程以及本文發(fā)展的基于IAPWS-IF97 水蒸汽性質(zhì)表插值的方法。對該算例計算50 步所用的時間對比結(jié)果列在表2 中??梢钥吹?,與采用理想氣體狀態(tài)方程相比較,采用FLUENT 內(nèi)置的維里型狀態(tài)方程,計算時間增加了180,而采用本文發(fā)展的方法,計算時間增加了27。表2 計算50 步所用計算時間對比圖4 給出了三種方法計算得到的壓力、溫度、密度等值線分布,可以看到計算結(jié)果之間仍然存在明顯差異。采用維里型狀態(tài)方程與采用基于IF 97 的水蒸汽表插值計算結(jié)果基本是吻合的,二者與采用理想氣體狀態(tài)方程得到的結(jié)果有明顯差異。考慮到計算所費的時間, 在大規(guī)模三維計算中,采用基于IAPWS IF 97 的水蒸汽性質(zhì)表插值比采用FLUENT 內(nèi)置的維里型狀態(tài)方程可以節(jié)省大量計算時間。(a) 壓力(Pa)(b) 溫度(K)(c) 密度(kg/m3) 理想氣體狀態(tài)方程 - 維里型狀態(tài)方程 IF 97 水蒸汽表插值圖4 參數(shù)等值線分布3 結(jié)論基于FLUENT 軟件提供的UDRGM 接口,發(fā)展了一種在大規(guī)模數(shù)值計算中采用實際氣體性質(zhì)進
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年保險業(yè)數(shù)字化理賠服務(wù)成本降低與效率提升策略報告
- 2025年保險數(shù)字化理賠服務(wù)移動端用戶體驗研究報告
- 2025至2031年中國綠豐三落王可濕性粉劑行業(yè)投資前景及策略咨詢研究報告
- 2025至2031年中國立體壓紋鋁塑板行業(yè)投資前景及策略咨詢研究報告
- 2025至2031年中國煤氣柜密封油行業(yè)投資前景及策略咨詢研究報告
- 《Windows Server操作系統(tǒng)維護與管理項目教程》課件學(xué)習(xí)情境7
- 幼兒教育在線平臺行業(yè)跨境出海項目商業(yè)計劃書
- 放射工作人員證考試:放射防護(上)
- 長治市保安證考試試題及答案
- 怎么考交安c考試試題及答案
- 2025年離婚協(xié)議書內(nèi)容
- 西湖大學(xué)《土木工程CAD》2023-2024學(xué)年第二學(xué)期期末試卷
- 公司工傷報銷管理制度
- 拆遷款收款協(xié)議書
- 血液透析??茝?fù)習(xí)題及參考答案
- 安全生產(chǎn)月知識競賽試題(含答案)
- 2025年物流管理與工程考試試題及答案
- 風(fēng)險管理工具的應(yīng)用試題及答案
- 小學(xué)新聞宣傳培訓(xùn):創(chuàng)新思路與實效傳播
- 期貨營業(yè)AI應(yīng)用行業(yè)深度調(diào)研及發(fā)展項目商業(yè)計劃書
- 2025年4月自考00420物理(工)試題及答案含評分標(biāo)準(zhǔn)
評論
0/150
提交評論