| 查看: 1868 | 回復(fù): 12 | |||
TigersZhao新蟲(chóng) (初入文壇)
|
[求助]
UDF求助 已有3人參與
|
||
|
我用下面的UDF代碼模擬吸附過(guò)程,為什么編譯能通過(guò),在進(jìn)行初始化時(shí)fluent會(huì)報(bào)錯(cuò):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; } |
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +60 |
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +60 |
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +60 |
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +60 |
新蟲(chóng) (初入文壇)
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +60 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +22 | RichLi_ 2026-03-25 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 375求調(diào)劑 +4 | 雨夏整夜 2026-03-29 | 4/200 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專(zhuān)碩-298 +4 | 手機(jī)用戶(hù) 2026-03-26 | 4/200 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研]
|
平樂(lè)樂(lè)樂(lè) 2026-03-26 | 4/200 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|