| 5 | 1/1 | 返回列表 |
| 查看: 2219 | 回復(fù): 12 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
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那里也要改?先謝謝各位大俠了~ |
銅蟲 (小有名氣)
銅蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +7 | pipiver 2026-03-30 | 7/350 |
|
|---|---|---|---|---|
|
[考研]
|
y7czhao 2026-03-26 | 12/600 |
|
|
[考研] 336求調(diào)劑 +5 | kiyy 2026-04-01 | 5/250 |
|
|
[考研] 279求調(diào)劑 +3 | 學(xué)而思兮知 2026-04-01 | 3/150 |
|
|
[考研] 265求調(diào)劑 +4 | 林深溫瀾 2026-04-01 | 6/300 |
|
|
[考研] 266分,一志愿電氣工程,本科材料,求材料專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-04-01 | 11/550 |
|
|
[考研] 085600,320分求調(diào)劑 +5 | 大饞小子 2026-04-01 | 6/300 |
|
|
[考研] 349求調(diào)劑 +6 | 吃的不少 2026-04-01 | 6/300 |
|
|
[基金申請(qǐng)] esi高被引論文是不是能對(duì)中標(biāo)有所加分和幫助呢 +3 | redcom 2026-04-01 | 3/150 |
|
|
[考研] 英一數(shù)二生物信息學(xué)287分,本科生物科學(xué),求調(diào)劑 +3 | 碧水xyz 2026-03-29 | 4/200 |
|
|
[考研] 0856初試324分求調(diào)劑 +5 | 想上學(xué)求調(diào) 2026-04-01 | 5/250 |
|
|
[考研] 生物學(xué)296求調(diào)劑 +10 | 湯圓包 2026-03-29 | 14/700 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +6 | 篤行致遠(yuǎn). 2026-03-28 | 6/300 |
|
|
[考研]
|
小羊36 2026-03-30 | 4/200 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|