 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 489 | 回復(fù): 2 | ||
steelbody新蟲 (初入文壇)
|
[求助]
載入U(xiǎn)DF之后無法初始化，求助高手
|
|
載入速度，湍動(dòng)能以及湍流耗散率的udf后無法初始化，求教，下面附上UDF，謝謝 #include "udf.h" #define YMIN 0.0 /* constants */ #define YMAX 0.17 #define UMEAN 1.0 #define B 1./7. #define DELOVRH 0.5 #define VISC 1.7894e-05 #define CMU 0.09 #define VKC 0.41 /* profile for x-velocity */ DEFINE_PROFILE(x_velocity, t, i) { real y, del, h, x[ND_ND], ufree; /* variable declarations */ face_t f; h = YMAX - YMIN; del = DELOVRH*h; ufree = 2.147; begin_f_loop(f, t) { F_CENTROID(x,f,t); y = x[1]; if (y <= del) F_PROFILE(f,t,i) = ufree*pow(y/del,B); else F_PROFILE(f,t,i) = ufree*pow((h-y)/del,B); } end_f_loop(f, t) } /* profile for kinetic energy */ DEFINE_PROFILE(k_profile, t, i) { real y, del, h, ufree, x[ND_ND]; real ff, utau, knw, kinf; face_t f; h = YMAX - YMIN; del = DELOVRH*h; ufree = 2.147; ff = 0.045/pow(ufree*del/VISC,0.25); utau=sqrt(ff*pow(ufree,2.)/2.0); knw=pow(utau,2.)/sqrt(CMU); kinf=0.002*pow(ufree,2.); begin_f_loop(f, t) { F_CENTROID(x,f,t); y=x[1]; if (y <= del) F_PROFILE(f,t,i)=knw+y/del*(kinf-knw); else F_PROFILE(f,t,i)=knw+(h-y)/del*(kinf-knw); } end_f_loop(f, t) } /* profile for dissipation rate */ DEFINE_PROFILE(dissip_profile, t, i) { real y, x[ND_ND], del, h, ufree; real ff, utau, knw, kinf; real mix, kay; face_t f; h = YMAX - YMIN; del = DELOVRH*h; ufree = 2.147; ff = 0.045/pow(ufree*del/VISC,0.25); utau=sqrt(ff*pow(ufree,2.)/2.0); knw=pow(utau,2.)/sqrt(CMU); kinf=0.002*pow(ufree,2.); begin_f_loop(f, t) { F_CENTROID(x,f,t); y=x[1]; if (y <= del) kay=knw+y/del*(kinf-knw); else kay=knw+(h-y)/del*(kinf-knw); if (VKC*y < 0.085*del) mix = VKC*y; else mix = 0.085*del; F_PROFILE(f,t,i)=pow(CMU,0.75)*pow(kay,1.5)/mix; } end_f_loop(f,t) } 出現(xiàn)問題如下 MBZ0BDJ)6NIO09USKTKQ[$9.png MBZ0BDJ)6NIO09USKTKQ[$9.png |
回復(fù)此樓» 猜你喜歡
085405 考的11408求各位老師帶走
已經(jīng)有3人回復(fù)
-
339求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有4人回復(fù)
-
085701環(huán)境工程，267求調(diào)劑
已經(jīng)有15人回復(fù)
-
材料求調(diào)劑一志愿哈工大324
已經(jīng)有9人回復(fù)
-
292求調(diào)劑
已經(jīng)有15人回復(fù)
-
304求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有17人回復(fù)
-
321求調(diào)劑
已經(jīng)有5人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
3_1415926
金蟲 (小有名氣)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1994.3
- 紅花: 6
- 帖子: 159
- 在線: 57.3小時(shí)
- 蟲號(hào): 1853296
- 注冊: 2012-06-09
- 性別: GG
- 專業(yè): 工程熱力學(xué)
2樓2017-04-06 12:54:22
steelbody
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 107.5
- 帖子: 20
- 在線: 7.2小時(shí)
- 蟲號(hào): 4641502
- 注冊: 2016-04-29
- 專業(yè): 流體力學(xué)
3樓2017-04-06 20:14:30
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 steelbody 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295材料工程專碩求調(diào)劑 +6 | 1428151015 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺_16 2026-03-24 | 4/200 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +8 | 北道巷 2026-03-24 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 316求調(diào)劑 +5 | 江辭666 2026-03-26 | 5/250 |
|
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
|
[考研] 321求調(diào)劑 +5 | 材料cailiao 2026-03-21 | 5/250 |
|
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|