| 7 | 1/1 | 返回列表 |
| 查看: 2494 | 回復(fù): 6 | |||||
lihongyang0金蟲 (初入文壇)
|
[交流]
FLuent 空化 UDF能量源項
|
|
我寫一個能量方程的源項,但是計算結(jié)果顯示,我的源項加沒加計算結(jié)果都一樣,證明說這里面有地方寫錯了。但是我實在看不出來哪里弄錯了。 我加載代碼的方式應(yīng)該是正確的吧(先iinterpreted,然后Cell zone, source前打勾,最后把能量UDF源選進(jìn)去) #include "udf.h" #include "sg_mphase.h" DEFINE_SOURCE(energy_source,c,m_t,dS,eqn) { Thread *l_t, *v_t; real m_dot, hc, source, fv, dT,m_dot2,source2; real cigma_m, cigma_m2,flag; real D_l, D_v,p_loc; real p_sat,p_sat2, T_loc, p_vapor,dp, dp2, T_loc2, p_vapor2; real p1, p2, p3, p4,m_dot_Re,m_dot_Re2,m_dot_Rc,m_dot_Rc2; real F_vap, F_cond, temp1, hc2; p_loc=C_P(c,m_t)*1000000; F_vap=0.02; F_cond=0.01; l_t=THREAD_SUB_THREAD(m_t,0); v_t=THREAD_SUB_THREAD(m_t,1); fv=C_VOF(c,v_t)*C_R(c,v_t)/(C_VOF(c,v_t)*C_R(c,v_t)+C_VOF(c,l_t)*C_R(c,l_t)); T_loc=C_T(c,m_t); p1=10.86; p2=-1980; p3=1.233*100000; p4=-2.617*1000000; p_sat=p1*T_loc*T_loc*T_loc+p2*T_loc*T_loc+p3*T_loc+p4; dT=0.01; T_loc2=T_loc+dT; p_sat2=p1*T_loc2*T_loc2*T_loc2+p2*T_loc2*T_loc2+p3*T_loc2+p4; p_vapor=p_sat+0.5*0.39*C_R(c,m_t)*C_K(c,m_t); p_vapor2=p_sat2+0.5*0.39*C_R(c,m_t)*C_K(c,m_t); p1=3.424e-09; p2=-1.867e-07; p3= -0.0002614; p4=0.02863; cigma_m=p1*T_loc*T_loc*T_loc+p2*T_loc*T_loc+p3*T_loc+p4; cigma_m2=p1*T_loc2*T_loc2*T_loc2+p2*T_loc2*T_loc2+p3*T_loc2+p4; p1=-0.0002051; p2=0.03527; p3=-3.064; p4=320.1; hc=p1*T_loc*T_loc*T_loc+p2*T_loc*T_loc+p3*T_loc+p4; hc2=p1*T_loc2*T_loc2*T_loc2+p2*T_loc2*T_loc2+p3*T_loc2+p4; dp=p_vapor-p_loc; dp2=p_vapor2-p_loc; temp1=sqrt(C_K(c,m_t)); D_l=C_R(c,l_t); D_v=C_R(c,v_t); m_dot_Re=0.0; m_dot_Rc=0.0; source=0.0; source2=0.0; flag=1000.0; if(p_loc<=p_vapor) { m_dot_Re=F_vap*temp1/cigma_m*D_l*D_v*sqrt(2/3*dp/D_l)*(1-fv); } if(p_loc>p_vapor) { m_dot_Rc=F_cond*temp1/cigma_m*D_l*D_l*sqrt(2/3*(-dp)/D_l)*fv; } source=(m_dot_Re-m_dot_Rc)*hc*flag; dS[eqn]=0.0; return source; } |
軟件及其學(xué)習(xí)資料 | fluent軟件學(xué)習(xí) |

金蟲 (初入文壇)
|
木蟲 (初入文壇)

鐵蟲 (初入文壇)

| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 +7 | kikiki7 2026-03-30 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +9 | Youth_ 2026-03-30 | 9/450 |
|
|
[考研] 材料專碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
[考研] 332求調(diào)劑 +13 | 032500 2026-03-25 | 13/650 |
|
|
[考研] 材料科學(xué)與工程調(diào)劑 +8 | 深V宿舍吧 2026-03-30 | 8/400 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7 | nhj_ 2026-03-25 | 7/350 |
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
[考研] 281求調(diào)劑 +4 | 亞克西good 2026-03-26 | 6/300 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|