| 1 | 1/1 | 返回列表 |
| 查看: 662 | 回復(fù): 0 | ||
shushanliu新蟲 (初入文壇)
|
[求助]
UDF實(shí)現(xiàn)周期性初始化錯(cuò)誤
|
|
我想要實(shí)現(xiàn)周期性變化,就是利用UDF將出口的參數(shù)賦值給進(jìn)口(下面是UDF程序),導(dǎo)進(jìn)fluent沒有問題,但是初始化時(shí)顯示Error: received a fatal signal (Segmentation fault).Error: received a fatal signal (Segmentation fault).Error Object: #f。全部的金幣,諸位大神看看是啥情況! #include "udf.h" #include "prop.h" #include "mem.h" int id_inlet=10;/*出口ID*/ int id_outlet=11;/*進(jìn)口ID*/ real avg_tempoutlet=293.15;/*出口初始溫度*/ real NV_VEC(A2); /*定義出口面積矢量*/ Thread *thread_inlet;/*定義進(jìn)口面索引*/ Thread *thread_outlet;/*定義出口面索引*/ face_t f; Domain *domain; real sum_Toutlet=0.0; real sum_Aoutlet=0.0; real tempoutlet=293.15; /*出口的實(shí)際溫度,節(jié)點(diǎn)溫度*/ real tw=323.15; /*壁面溫度*/ real theathoutlet=1.0;/*出口無量綱溫度*/ real uoutlet=0.25; real voutlet=0.0; real x[ND_ND]; DEFINE_ADJUST(myadjust,d) /*函數(shù)用于每一步迭代進(jìn)行調(diào)整*/ { domain=Get_Domain(1); /*指向區(qū)域的指針,可訪問mesh里面的所有單元和面索引*/ thread_inlet=Lookup_Thread(domain,id_inlet); /*獲取進(jìn)口的ID,線索*/ thread_outlet=Lookup_Thread(domain,id_outlet); /*獲取出口的ID,線索*/ begin_f_loop(f,thread_outlet) { F_AREA(A2,f,thread_outlet); tempoutlet=F_T(f,thread_outlet); sum_Aoutlet+=NV_MAG(A2); sum_Toutlet+=NV_MAG(A2)*tempoutlet; } end_f_loop(f,thread_outlet) /*計(jì)算出口平均溫度*/ //avg_tempinlet=sum_Tinlet/sum_Ainlet; /*進(jìn)口平均溫度*/ avg_tempoutlet=sum_Toutlet/sum_Aoutlet; /*出口平均溫度*/ begin_f_loop(f,thread_outlet) { F_AREA(A2,f,thread_outlet); F_CENTROID(x,f,thread_outlet); tempoutlet=F_T(f,thread_outlet); theathoutlet=(tempoutlet-tw)/(avg_tempoutlet-tw); F_UDMI(f,thread_outlet,1)=theathoutlet; F_UDMI(f,thread_outlet,2)=F_U(f,thread_outlet); F_UDMI(f,thread_outlet,3)=F_V(f,thread_outlet); } end_f_loop(f,thread_outlet)/*計(jì)算出口的無量綱溫度*/ } DEFINE_PROFILE(inlet_temp,thread_inlet,index) { real theathinlet=1.0; real tempinlet=293.15; real avg_tempinlet=293.15; begin_f_loop(f,thread_inlet) { avg_tempinlet=avg_tempoutlet; theathinlet=F_UDMI(f,thread_outlet,1); tempinlet=theathinlet*(avg_tempinlet-tw)+tw; F_PROFILE(f,thread_inlet,index)=tempinlet; /*進(jìn)口溫度邊界設(shè)置為出口的溫度*/ end_f_loop(f,thread_inlet) } } DEFINE_PROFILE(u_velocity,thread_inlet,index) { real uinlet=0.02;/*定義初始速度u*/ begin_f_loop(f,thread_inlet) { uinlet=F_UDMI(f,thread_outlet,2); F_PROFILE(f,thread_inlet,index)=uinlet; /*進(jìn)口u速度邊界設(shè)置為出口的速度*/ end_f_loop(f,thread_inlet) } } DEFINE_PROFILE(v_velocity,thread_inlet,index) { real vinlet=0.0;/*定義初始速度v*/ begin_f_loop(f,thread_inlet) { vinlet=F_UDMI(f,thread_outlet,3); F_PROFILE(f,thread_inlet,index)=vinlet; /*進(jìn)口v速度邊界設(shè)置為出口的速度*/ end_f_loop(f,thread_inlet) } }@wuming524 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083000,總分284,求調(diào)劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿華中科大0702,數(shù)一英一,293 +3 | 小羅露一二 2026-03-07 | 3/150 |
|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +4 | 鋼鐵大炮 2026-03-04 | 4/200 |
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +4 | 蟲友233 2026-03-07 | 4/200 |
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
[考研] 328求調(diào)劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 290求調(diào)劑 +5 | Bananaiy 2026-03-04 | 5/250 |
|
|
[考研] 一志愿東北大學(xué)080500,英一數(shù)二343分求調(diào)劑 +5 | Seele_v5 2026-03-04 | 5/250 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-04 | 3/150 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 成績(jī)276,專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|