| 7 | 1/1 | 返回列表 |
| 查看: 2493 | 回復(fù): 6 | |||||
lihongyang0金蟲 (初入文壇)
|
[交流]
FLuent 空化 UDF能量源項(xiàng)
|
|
我寫一個能量方程的源項(xiàng),但是計(jì)算結(jié)果顯示,我的源項(xiàng)加沒加計(jì)算結(jié)果都一樣,證明說這里面有地方寫錯了。但是我實(shí)在看不出來哪里弄錯了。 我加載代碼的方式應(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 295分求調(diào)劑 +4 | W55j 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +11 | 匹克i 2026-03-23 | 11/550 |
|
|
[考研] 材料科學(xué)與工程調(diào)劑 +8 | 深V宿舍吧 2026-03-30 | 8/400 |
|
|
[考研] 277跪求調(diào)劑 +7 | 1915668 2026-03-27 | 11/550 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +25 | RichLi_ 2026-03-25 | 25/1250 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|