| 5 | 1/1 | 返回列表 |
| 查看: 2497 | 回復(fù): 6 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||
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í) |
鐵蟲 (初入文壇)


金蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 290分調(diào)劑求助 +4 | 吉祥止止陳 2026-03-25 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +5 | 里子木yy 2026-03-29 | 5/250 |
|
|
[考研] 362求調(diào)劑 +3 | 西南交材料專碩3 2026-03-31 | 3/150 |
|
|
[考研] 085600材料與化工329分求調(diào)劑 +5 | 這是有太陽哇 2026-03-25 | 5/250 |
|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑,不限專業(yè),有學(xué)上就行 +11 | lllllos 2026-03-30 | 11/550 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +7 | 蘿卜燉青菜 2026-03-28 | 8/400 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|