| 查看: 2753 | 回复: 6 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
lihongyang0金虫 (初入文坛)
|
[交流]
FLuent 空化 UDF能量源项
|
||||
|
我写一个能量方程的源项，但是计算结果显示，我的源项加没加计算结果都一样，证明说这里面有地方写错了。但是我实在看不出来哪里弄错了。 我加载代码的方式应该是正确的吧（先iinterpreted，然后Cell zone， source前打勾，最后把能量UDF源选进去） #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; } |
回复此楼» 收录本帖的淘帖专辑推荐
 软件及其学习资料 | fluent软件学习 |
» 猜你喜欢
今天务委会开完了，明天出结果吗
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有9人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有6人回复
-
放榜前的不淡定
已经有18人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
-
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一下大神们关于fluent中UDF的问题，谢谢！
已经有3人回复
- fluent UDF 源项含义
已经有4人回复
- 汽蚀（空化）UDF编写中遇到的问题
已经有13人回复
- fluent热源udf
已经有2人回复
- fluent在瞬态计算时候如何每次迭代都调用一次源项的udf
已经有6人回复
- fluent udf 热源
已经有9人回复
- FLUENT两相流空化模拟
已经有3人回复
- udf添加能量源项
已经有3人回复
- fluent 中二维模型中，材料属性随着温度变化的udf源代码？
已经有14人回复
- 一个fluent UDF问题
已经有4人回复
- fluent能量源项中温度T对时间的t的偏导数，即dT/dt如何用udf编写能量源项
已经有6人回复
- fluent物理模型及相关定义和使用
已经有90人回复
- Fluent中若干疑问
已经有5人回复
- FLUENT源项UDF非稳态迭代的一个疑问请教
已经有8人回复
- 【求助】通气空泡多相流（气、水蒸汽、水）空化模型的UDF如何入手啊？？？？
已经有4人回复
3楼2014-06-05 09:17:45
2楼2014-06-04 11:10:10
lihongyang0
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1057.2
- 散金: 18
- 红花: 1
- 帖子: 37
- 在线: 148.5小时
- 虫号: 2268108
- 注册: 2013-01-27
- 性别: GG
- 专业: 可再生与替代能源利用中的
4楼2014-06-05 09:34:03
独恋sophia
木虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 4044.5
- 帖子: 38
- 在线: 18.4小时
- 虫号: 1628317
- 注册: 2012-02-20
- 性别: MM
- 专业: 工程热物理相关交叉领域
6楼2014-11-26 11:01:34
