| 查看: 1969 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
UDF求助 已有3人参与
|
|||
|
我用下面的UDF代码模拟吸附过程，为什么编译能通过，在进行初始化时fluent会报错：received fatal signal (segmentation fault)!求大神帮忙解答一下！十分感谢！ #include "udf.h" enum { p, ads, d_ads, H_ad }; DEFINE_ADJUST(p_absolute,domain) /*calculate the absolute pressure*/ { Thread *t; cell_t c; real p_operating; real p_static; thread_loop_c(t,domain) { if (NULL!=THREAD_STORAGE(t,SV_UDS_I(p))) { begin_c_loop(c,t) { p_static=C_P(c,t); p_operating=RP_Get_Real("operating-pressure"  ; C_UDSI(c,t,p)=p_static+p_operating; } end_c_loop(c,t) } } } DEFINE_ADJUST(Heat_Adsorption,domain) /*calculate the released heat during adsorption*/ { Thread *t; cell_t c; real p_operating; real p_static; thread_loop_c(t,domain) { if(NULL!=THREAD_STORAGE(t,SV_UDS_I(H_ad))&&NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(H_ad))) { begin_c_loop(c,t) { C_UDSI(c,t,H_ad)=3080*sqrt(log(71.6/C_UDSI_M1(c,t,ads))); } end_c_loop(c,t) } } } DEFINE_ADJUST(adsorption,domain) /* calculate the absolute adsorption quantity*/ { Thread *t; cell_t c; real n_o; real E_t; real x_p; real A; real x_b; real n; thread_loop_c(t,domain) { if(NULL!=THREAD_STORAGE(t,SV_UDS_I(ads))&&NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(ads))) { begin_c_loop(c,t) { real tem=C_T(c,t); n_o=71.6; E_t=3080+18.9*tem; x_p=1.47E9/C_UDSI(c,t,p); A=8.31429*tem*log(x_p); x_b=-A*A/(E_t*E_t); C_UDSI(c,t,ads)=n_o*exp(x_b); } end_c_loop(c,t) } } } DEFINE_ADJUST(d_adsorption,domain) /*calculate the adsorption rate*/ { Thread *t; cell_t c; thread_loop_c(t,domain) { if(NULL != THREAD_STORAGE(t,SV_UDS_I(d_ads))&&NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(d_ads))) { begin_c_loop(c,t) { real tem=C_T(c,t); real physical_dt=RP_Get_Real("physical-time-step" ;C_UDSI(c,t,d_ads)=-0.532224*C_UDSI(c,t,ads)*(2*pow((8.314*tem/(3080+18.9*tem)), 2)*log(1470*1000000/C_UDSI(c,t,p))/C_UDSI(c,t,p)*(C_UDSI(c,t,p)-C_UDSI_M1(c,t,p))/physical_dt-2*pow(8.314,2)*3080*tem/pow((3080+18.9*tem),3)*pow(log( 1470*1000000/C_UDSI(c,t,p)),2)*(C_T(c,t)-C_T_M1(c,t))/physical_dt); } end_c_loop(c,t) } } } DEFINE_SOURCE(mass_source,c,t,dS,eqn) { real m_source; m_source=C_UDSI_M1(c,t,d_ads); dS[eqn]=0.0; return m_source; } DEFINE_SOURCE(energy_source,c,t,dS,eqn) { real e_source; e_source=-C_UDSI_M1(c,t,d_ads)*0.496056352e3*30803*sqrt(log(71.6/C_UDSI_M1(c,t,ads))); dS[eqn]=0.0; return e_source; } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
3楼2014-10-16 14:29:22
死天
专家顾问 (著名写手)
-
专家经验: +60 - 应助: 399 (硕士)
- 金币: 3648.9
- 红花: 89
- 帖子: 1336
- 在线: 443.8小时
- 虫号: 2643057
- 注册: 2013-09-09
- 专业: 流体力学
- 管辖: 仿真模拟
2楼2014-10-16 13:25:38
死天
专家顾问 (著名写手)
-
专家经验: +60 - 应助: 399 (硕士)
- 金币: 3648.9
- 红花: 89
- 帖子: 1336
- 在线: 443.8小时
- 虫号: 2643057
- 注册: 2013-09-09
- 专业: 流体力学
- 管辖: 仿真模拟
4楼2014-10-16 14:40:57
5楼2014-10-16 15:07:53
