| 4 | 1/1 | 返回列表 |
| 查看: 1283 | 回復(fù): 3 | ||
白白的大蘿卜新蟲 (小有名氣)
|
[求助]
UFD調(diào)試問題! 已有1人參與
|
在Fluent 中加入源項(xiàng)的UDF 一初始化就開始報錯 Error: received a fatal signal(Segmentation fault).,UDF是照著別人的文獻(xiàn)寫得,沒有書寫錯誤,現(xiàn)在不知道問題出在哪里,本人對UDF和C也是剛?cè)腴T的新手,在這里 希望好心人能幫忙調(diào)試一下,有償也可以的。!![]() ![]() ![]() ![]() 救救孩子吧 萬分感謝!#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; } |
新蟲 (小有名氣)
新蟲 (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085602高分子方向求調(diào)劑 +7 | tlgudy 2026-03-04 | 7/350 |
|
|
[考研] 282求調(diào)劑 +7 | 夕~日 2026-03-05 | 8/400 |
|
|
[考研] 材料與化工,291,求調(diào)劑 +7 | 咕嚕咕嚕123123 2026-03-05 | 8/400 |
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
[考研] 材料類考研調(diào)劑 +7 | gemmgemm 2026-03-01 | 8/400 |
|
|
[考研] 22408-273求調(diào)劑-擔(dān)任3個項(xiàng)目的負(fù)責(zé)人-1篇國際期刊論文(一作)1篇核心期刊論文在投。 +3 | 沒想好取什么名 2026-03-03 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 曼殊2266 2026-02-28 | 8/400 |
|
|
[考研] 環(huán)境調(diào)劑 +7 | 柒槿levana 2026-03-01 | 7/350 |
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 一志愿西工大材料學(xué),英一數(shù)二總分321分,求調(diào)劑。 +3 | zz05zz 2026-03-02 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | 杰尼龜aaa 2026-03-03 | 4/200 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[基金申請] 請問大家,研究風(fēng)險與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
[考研] 272求調(diào)劑 +9 | 材紫有化 2026-02-28 | 9/450 |
|
|
[考研] 一志愿東北大學(xué)材料專碩328,求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
[考研] 1 +5 | 2026計(jì)算機(jī)_誠心 2026-03-01 | 8/400 |
|