蒸發(fā)過程UDF程序(fluent)_第1頁
蒸發(fā)過程UDF程序(fluent)_第2頁
蒸發(fā)過程UDF程序(fluent)_第3頁
蒸發(fā)過程UDF程序(fluent)_第4頁
蒸發(fā)過程UDF程序(fluent)_第5頁
全文預覽已結(jié)束

下載本文檔

版權(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. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論