wzhhgd3210

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 71.9
散金: 121
紅花: 2
帖子: 177
在線: 101小時
蟲號: 938276
注冊: 2010-01-06
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

[求助] UDS及其后處理求助！已有2人參與

樓主剛開始接觸CFD，現(xiàn)在模擬室內(nèi)空調(diào)通風(fēng)，需要利用每個節(jié)點上已計算出的溫度和風(fēng)速計算求解該點的PMV值并輸出云圖（該值可以看做室內(nèi)溫度和風(fēng)速的函數(shù)），看到paper上都說用UDF做，我在網(wǎng)上找了一個開源的修改了一下，可以interpret但是迭代計算就會報錯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)
}
}

請各位大俠幫忙看一下問題在哪，另外樓主是小白，還想請教一下UDS的具體操作：先interpret代碼，然后在Funtion Hooks里的Adjust選項里選中我的UDF激活，再把User-Difined Scalars里面Number添加1，設(shè)置邊界條件迭代。這樣的操作就可以了么？還是說Memory那里也要改？先謝謝各位大俠了~

回復(fù)此樓

» 猜你喜歡

一志愿東華大學(xué)化學(xué)070300，求調(diào)劑已經(jīng)有8人回復(fù)
277材料科學(xué)與工程080500求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有18人回復(fù)
293求調(diào)劑已經(jīng)有5人回復(fù)
280分求調(diào)劑一志愿085802 已經(jīng)有7人回復(fù)
0854電子信息求調(diào)劑已經(jīng)有3人回復(fù)
263求調(diào)劑已經(jīng)有4人回復(fù)
石河子大學(xué)（211、雙一流）碩博研究生長期招生公告已經(jīng)有3人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有13人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求助治療抑郁的欣百達(dá)(鹽酸度洛西汀腸溶膠囊的說明書，最好是PDF版本的！已經(jīng)有8人回復(fù)
【求助】與化學(xué)藥品相比，中藥及其制劑在分析前處理,鑒別檢查,含量測定的特點已經(jīng)有6人回復(fù)

1樓 2013-05-02 21:11:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

匿名

用戶注銷 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 989
散金: 666
紅花: 1
帖子: 322
在線: 65.6小時
蟲號: 0
注冊: 2012-05-10
專業(yè): 恒星與星際物質(zhì)

本帖僅樓主可見

贊一下

2樓2015-01-03 11:40:53

已閱申請仿真EPI 回復(fù)此樓編輯查看我的主頁

nhnt

銅蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 456.8
散金: 105
帖子: 87
在線: 60.1小時
蟲號: 3129029
注冊: 2014-04-11
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

引用回帖:

3樓: Originally posted by soriyoshi at 2015-01-03 14:26:15
剛剛鼓搗出了一種解決方法，為了方便以后看到此貼的有疑問的人，暫時記在這里�？梢杂肈EFINE_ADJUST，不過是用C_UDMI，然后在define-user define memory 里面把數(shù)量加成1，計算完以后再contour里面自然有一個udm， ...

我也要做關(guān)于PMV的udf，可是剛接觸，我按照這個代碼和你的方法，迭代的時候還是有錯誤，能不能給個詳細(xì)的操作��？

贊一下(1人)

回復(fù)此樓

5樓2015-06-27 16:20:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

匿名

用戶注銷 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 989
散金: 666
紅花: 1
帖子: 322
在線: 65.6小時
蟲號: 0
注冊: 2012-05-10
專業(yè): 恒星與星際物質(zhì)

本帖僅樓主可見

贊一下

3樓2015-01-03 14:26:15

已閱申請仿真EPI 回復(fù)此樓編輯查看我的主頁

nhnt

銅蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 456.8
散金: 105
帖子: 87
在線: 60.1小時
蟲號: 3129029
注冊: 2014-04-11
性別: GG
專業(yè): 傳熱傳質(zhì)學(xué)

樓主向你求助下啊，我按照你說的方法，還是會出現(xiàn)錯誤提示，能不能給說下詳細(xì)的過程？

贊一下

回復(fù)此樓

4樓2015-06-26 11:10:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyan2387

捐助貴賓 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 108
帖子: 105
在線: 67.9小時
蟲號: 2679233
注冊: 2013-09-25
性別: MM
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

3樓: Originally posted by soriyoshi at 2015-01-03 14:26:15
剛剛鼓搗出了一種解決方法，為了方便以后看到此貼的有疑問的人，暫時記在這里。可以用DEFINE_ADJUST，不過是用C_UDMI，然后在define-user define memory 里面把數(shù)量加成1，計算完以后再contour里面自然有一個udm， ...

親可否分享一下你的PMV的 UDF呢我也是想做熱環(huán)境舒適性想做PMV的云圖

贊一下

回復(fù)此樓

Wuha

6樓2016-08-31 17:12:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyan2387

捐助貴賓 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 108
帖子: 105
在線: 67.9小時
蟲號: 2679233
注冊: 2013-09-25
性別: MM
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

親我的UDF編寫的時候在循環(huán)迭代的地方總是提示參量多可否分享一下您的代碼謝謝啦

贊一下

回復(fù)此樓

Wuha

7樓2016-08-31 17:12:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

得得315

禁蟲 (正式寫手)

本帖內(nèi)容被屏蔽

8樓2016-09-09 16:50:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

靜靜長大123

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 149.8
散金: 40
帖子: 136
在線: 3.2小時
蟲號: 5052073
注冊: 2016-09-23
專業(yè): 計算機(jī)軟件

引用回帖:

求樓主寫一下具體的設(shè)置的項目和參數(shù)，寫好了UDF但是就是算不出來，很著急啊

贊一下

回復(fù)此樓

9樓2017-07-04 15:34:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

MoonLee777

新蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 2274.6
紅花: 17
帖子: 839
在線: 78小時
蟲號: 5947368
注冊: 2017-03-12
性別: GG
專業(yè): 流體力學(xué)

好像是udm要保存、不然每迭代一次上一次就消失了

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

10樓2017-07-06 00:10:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 求調(diào)劑院校信息 +6	CX 330 2026-03-21	6/300	2026-03-22 15:25 by 無懈可擊111
[考研] 269專碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 354求調(diào)劑 +7	Tyoumou 2026-03-18	10/500	2026-03-22 11:11 by 人來盛
[考研] 資源與環(huán)境調(diào)劑申請(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂