




下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領
文檔簡介
1、Fluent 蒸發(fā)相變 模擬 UDF經(jīng)過幾天的不懈折騰,終于找到一個較為完成的用于fluent蒸發(fā)相變模擬的udf的一個程序。而且注釋相對完整。#include udf.h /包括常規(guī)宏#include sg_mphase.h / 包括體積分數(shù)宏CVOF(C,T)#define T_SAT 373 /定義蒸發(fā)溫度100#define LAT_HT 1.e3 /定義蒸發(fā)潛熱J/KgDEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) /液相質(zhì)量源項UDFThread *mix_th, *sec_th; /定義計算區(qū)線指針real m_dot_l; /定義液相
2、質(zhì)量轉(zhuǎn)移 kg/(m2.s)mix_th = THREAD_SUPER_THREAD(pri_th); /指向混合區(qū)的主相即液相的指針sec_th = THREAD_SUB_THREAD(mix_th, 1); /指向單相控制區(qū)的氣相的指針,氣相為第二相if(C_T(cell, mix_th)=T_SAT) /如果液相單元的溫度高于蒸發(fā)溫度,液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移 m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dSeqn = -0.1*C_R(cell,
3、 pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; /定義源項對質(zhì)量轉(zhuǎn)移偏導 else m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)* fabs(T_SAT-C_T(cell,mix_th)/T_SAT; /如果指向混合區(qū)液相的單元溫度小于蒸發(fā)溫度,氣相向液相的質(zhì)量轉(zhuǎn)移,液相得 dSeqn = 0.;/由于是氣相向液相轉(zhuǎn)移,所以液相的質(zhì)量源項對質(zhì)量轉(zhuǎn)移的偏導為零 return m_dot_l;DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /氣相質(zhì)量
4、源項UDFThread * mix_th, *pri_th;real m_dot_v;mix_th = THREAD_SUPER_THREAD(sec_th); /指向混合區(qū)的第二相即氣相的指針pri_th = THREAD_SUB_THREAD(mix_th, 0); 指向單相控制區(qū)的液相的指針,液相為主相if(C_T(cell, mix_th)=T_SAT) /如果混合區(qū)單元的溫度高于蒸發(fā)溫度,液相向氣相的質(zhì)量質(zhì)量轉(zhuǎn)移 m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_
5、SAT; dSeqn = 0.;由于是液相向氣相轉(zhuǎn)移,所以氣相的質(zhì)量源項對來自液相的質(zhì)量轉(zhuǎn)移的偏導為零 else m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th)/T_SAT; /如果指向混合區(qū)的單元溫度小于蒸發(fā)溫度,氣相向液相的質(zhì)量轉(zhuǎn)移,氣相失 dSeqn = -0.1*C_R(cell, sec_th)* fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; /由于是氣相向液相轉(zhuǎn)移,所以氣相的質(zhì)量源項對自身的質(zhì)量轉(zhuǎn)移的偏導不為零 return m_dot_
6、v;DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) /混合模型能量源項UDFThread *pri_th, *sec_th;real m_dot;pri_th = THREAD_SUB_THREAD(mix_th, 0);/指向混合區(qū)的液相的指針sec_th = THREAD_SUB_THREAD(mix_th, 1); /指向混合區(qū)的氣相的指針if(C_T(cell, mix_th)=T_SAT) /如果混合區(qū)的單元溫度高于蒸發(fā)溫度。質(zhì)量轉(zhuǎn)移由液相向氣相轉(zhuǎn)移,吸熱,質(zhì)量轉(zhuǎn)移量前有負號 m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; dSeqn = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT; else /相反,氣相向液相轉(zhuǎn)移則放熱 m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*fabs(T_SAT-C_T(cell,mix_th)/T_SAT;
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025至2030中國白色家電行業(yè)市場運行分析及競爭格局與投資方向報告
- 2025至2030中國男士商務正裝行業(yè)深度研究及發(fā)展前景投資評估分析
- 2025至2030中國用于食品和飲料的金屬罐行業(yè)產(chǎn)業(yè)運行態(tài)勢及投資規(guī)劃深度研究報告
- 2025至2030中國玻璃門行業(yè)深度研究及發(fā)展前景投資評估分析
- 2025至2030中國玫瑰花露行業(yè)供需分析及發(fā)展前景報告
- 2025至2030中國物理治療軟件行業(yè)市場深度研究及發(fā)展前景投資可行性分析報告
- 商業(yè)培訓中激勵措施的心理機制研究
- 商業(yè)環(huán)境中殘疾人餐具使用的培訓與指導
- 招聘技巧培訓課件
- 智能教育設備應用中的隱私保護問題研究
- 干部履歷表(中共中央組織部2015年制)
- JCT1041-2007 混凝土裂縫用環(huán)氧樹脂灌漿材料
- SPA水療管理手冊
- 充電樁工程施工方案解決方案
- 7、煤礦安全管理二級質(zhì)量標準化驗收標準
- USSF-美國太空部隊數(shù)字服務遠景(英文)-2021.5-17正式版
- 靜配中心應急預案處理流程
- 江蘇省射陽中等專業(yè)學校工作人員招聘考試真題2022
- 廣東英語中考必背1600詞
- 醫(yī)療器械銷售代表工作計劃工作總結(jié)述職報告PPT模板下載
評論
0/150
提交評論