| 查看: 2215 | 回復(fù): 12 | ||
wzhhgd3210銀蟲 (小有名氣)
|
[求助]
UDS及其后處理求助! 已有2人參與
|
|
|
樓主剛開始接觸CFD,現(xiàn)在模擬室內(nèi)空調(diào)通風(fēng),需要利用每個(gè)節(jié)點(diǎn)上已計(jì)算出的溫度和風(fēng)速計(jì)算求解該點(diǎn)的PMV值并輸出云圖(該值可以看做室內(nèi)溫度和風(fēng)速的函數(shù)),看到paper上都說用UDF做,我在網(wǎng)上找了一個(gè)開源的修改了一下,可以interpret但是迭代計(jì)算就會(huì)報(bào)錯(cuò)Fluent received fatal signal(ACCESS_VIOLATION)。 代碼如下: #include "udf.h" enum { pmv }; DEFINE_ADJUST(pmv_myudf,domain) { Thread *t; cell_t c; thread_loop_c(t,domain) { begin_c_loop(c,t) { double t_myudf=C_T(c,t)-273.15; double speed_u_myudf=C_U(c,t); double speed_v_myudf=C_V(c,t); double speed_w_myudf=C_W(c,t); double mh2o_myudf=C_YI(c,t,0); double clo=0.55; double icl_myudf=0.08525; double fcl_myudf; double tcl_myudf; double tcl1_myudf; double tcl2_myudf; double temperary1_myudf; double temperary2_myudf; double hc_myudf; double tr_myudf=t_myudf; double M_myudf=58.0; double W_myudf=0.0; double p_myudf=101325; double pa_myudf; double a_myudf; double b_myudf; double c_myudf; double d_myudf; double e_myudf; double f_myudf; double L_myudf; double speed_myudf=sqrt(pow(speed_u_myudf,2.0)+pow(speed_v_myudf,2.0)+pow(speed_w_myudf,2.0)); pa_myudf=29*mh2o_myudf/(18+11*mh2o_myudf)*p_myudf; if (icl_myudf<0.078) fcl_myudf=1.00+1.290*icl_myudf; else fcl_myudf=1.05+0.645*icl_myudf; tcl1_myudf=40; hc_myudf=1; tcl2_myudf=35.7-0.025*(M_myudf-W_myudf)-icl_myudf*(0.0000000396*fcl_myudf*(pow((tcl1_myudf+273),4)-pow((tr_myudf+273),4))+fcl_myudf*hc_myudf*(tcl1_myudf-t_myudf)); while (tcl1_myudf-tcl2_myudf>0.000001) { tcl1_myudf=tcl2_myudf; temperary1_myudf=2.38*pow((tcl1_myudf-t_myudf),0.25); temperary2_myudf=12.1*pow(speed_myudf,0.5); if (temperary1_myudf<temperary2_myudf) hc_myudf=temperary2_myudf; else hc_myudf=temperary1_myudf; tcl2_myudf=35.7-0.025*(M_myudf-W_myudf)-icl_myudf*(0.0000000396*fcl_myudf*(pow((tcl1_myudf+273),4)-pow((tr_myudf+273),4))+fcl_myudf*hc_myudf*(tcl1_myudf-t_myudf)); } tcl_myudf=tcl2_myudf; a_myudf=0.0000000396*fcl_myudf*(pow((tcl_myudf+273),4.0)-pow((tr_myudf+273),4.0)); b_myudf=fcl_myudf*hc_myudf*(tcl_myudf-t_myudf); c_myudf=0.00305*(5733-6.99*(M_myudf-W_myudf)-pa_myudf); d_myudf=0.42*(M_myudf-W_myudf-58.15); e_myudf=0.000017*M_myudf*(5867-pa_myudf); f_myudf=0.0014*M_myudf*(34-t_myudf); L_myudf=M_myudf-W_myudf-(a_myudf+b_myudf+c_myudf+d_myudf+e_myudf+f_myudf); C_UDSI(c,t,pmv)=(0.303*exp(-0.036*M_myudf)+0.028)*L_myudf; } end_c_loop(c,t) } } 請(qǐng)各位大俠幫忙看一下問題在哪,另外樓主是小白,還想請(qǐng)教一下UDS的具體操作:先interpret代碼,然后在Funtion Hooks里的Adjust選項(xiàng)里選中我的UDF激活,再把User-Difined Scalars里面Number添加1,設(shè)置邊界條件迭代。這樣的操作就可以了么?還是說Memory那里也要改?先謝謝各位大俠了~ |
銅蟲 (小有名氣)
銅蟲 (小有名氣)
捐助貴賓 (小有名氣)

捐助貴賓 (小有名氣)

|
本帖內(nèi)容被屏蔽 |
新蟲 (小有名氣)
新蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +6 | chuanzhu川燭 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|