| 查看: 1524 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
UFD调试问题！！ 已有1人参与
|
||
在Fluent 中加入源项的UDF 一初始化就开始报错 Error: received a fatal signal(Segmentation fault).，UDF是照着别人的文献写得，没有书写错误，现在不知道问题出在哪里，本人对UDF和C也是刚入门的新手，在这里 希望好心人能帮忙调试一下，有偿也可以的！！！ 救救孩子吧 万分感谢！#include "udf.h" /*Define which user-defined scalars to use.*/ enum { p,ads,d_ads,q_ads }; DEFINE_ADJUST(p_absolute,domain) { Thread*t; cell_t c; real p_operating; real p_static; /*solve the absolute pressure.*/ 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(adsorption,domain) { Thread*t; cell_t c; real n_o; real E_t; real x_p; real A; real x_b; real n; /*solve the adsorption equation.*/ 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(q_adsorption,domain) { Thread*t; cell_t c; /*solve the absolute pressure.*/ thread_loop_c(t,domain) { if(NULL!=THREAD_STORAGE(t,SV_UDS_I(q_ads))&&NULL!=T_STORAGE_R_NV(t,SV_UDSI_M1(q_ads))) { begin_c_loop(c,t) { real physical_dt; physical_dt=RP_Get_Real("physical-time-step" ;C_UDSI(c,t,q_ads)=C_UDSI_M1(c,t,q_ads)+0.001*physical_dt*(C_UDSI_M1(c,t,ads)-C_UDSI_M1(c,t,q_ads)); } end_c_loop(c,t) } } } DEFINE_ADJUST(d_adsorption,domain) { Thread*t; cell_t c; /*solve the energy source.*/ 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) { C_UDSI(c,t,d_ads)=-0.5422633*0.001*(C_UDSI(c,t,ads)-C_UDSI(c,t,q_ads)); } 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(enery_source,c,t,dS,eqn) { real e_source; e_source=-C_UDSI_M1(c,t,d_ads)*4*1.58E6; dS[eqn]=0.0; return e_source; } |
回复此楼» 猜你喜欢
有没有大神帮我看看基金代码
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有13人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
-
没有任何消息-是不是就凉了
已经有10人回复
4楼2021-01-23 10:02:02
白白的大萝卜2楼2019-12-02 13:54:29
3楼2019-12-02 15:10:55
