 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (3782)
- >導(dǎo)師招生 (500)
- >考研 (381)
- >文獻(xiàn)求助 (199)
- >蟲(chóng)友互識(shí) (131)
- >休閑灌水 (120)
- >考博 (99)
- >基金申請(qǐng) (81)
- >招聘信息布告欄 (46)
- >碩博家園 (40)
- >綠色求助(高懸賞) (31)
- >博后之家 (26)
- >論文投稿 (25)
- >SciFinder/Reaxys (17)
- >找工作 (9)
- >農(nóng)林 (8)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1166 | 回復(fù): 4 | |||
jiangge2585銅蟲(chóng) (小有名氣)
|
[求助]
udf 熱源 初始化
|
|
求大神幫我看看我的程序哪里出問(wèn)題了?這是一個(gè)激光高斯面熱源程序，我的意圖是：激光脈寬0.008s，頻率是10Hz，就是1秒激光打10次，每次持續(xù)0.008s。但是 模擬的結(jié)果卻是不能初始化，初始化時(shí)就卡主了，一直停在那。求大神救我！ 這是程序： DEFINE_PROFILE(heat_flux,t,i) { real xx[ND_ND]; real y,x,r,a,time_step; int n; face_t f; begin_f_loop(f,t) { F_CENTROID(xx,f,t); y=xx[1]; x=xx[0]; r=0.001; a=0; n=0; if(a<=1) { if(a>=0.108*n&&a<=0.108*n+0.08) { if(sqrt(x*x+y*y)<=0.001) { F_PROFILE(f,t,i)=1160*0.7/(3.1415926*r*r)*exp(-(x*x+y*y)/(r*r)); } else { F_PROFILE(f,t,i)=0; } RP_Get_Real("flow-time"  ;a=RP_Get_Real("flow-time" ;n=(int)a/0.108; } else { F_PROFILE(f,t,i)=0; } } else { Error("error reading file" ;} } end_f_loop(f,t) } |
回復(fù)此樓» 猜你喜歡
333求調(diào)劑
已經(jīng)有10人回復(fù)
-
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
考研化學(xué)308分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298求調(diào)劑
已經(jīng)有7人回復(fù)
-
263求調(diào)劑
已經(jīng)有5人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷
已經(jīng)有6人回復(fù)
-
求調(diào)劑材料學(xué)碩080500，總分289分
已經(jīng)有21人回復(fù)
-
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有8人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有9人回復(fù)
-
291 求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- FLUENT中編寫(xiě)UDF時(shí)指針變量是否要初始化
已經(jīng)有4人回復(fù)
- FLUENT UDF 編譯通過(guò)，初始化出錯(cuò)的原因？？？求助
已經(jīng)有4人回復(fù)
- UDF 熱源隨溫度變化的處理
已經(jīng)有4人回復(fù)
- 請(qǐng)教：模擬帶內(nèi)熱源的管道流動(dòng)
已經(jīng)有4人回復(fù)
- 使用UDF后，初始化錯(cuò)誤
已經(jīng)有11人回復(fù)
- UDF應(yīng)用之后的初始化錯(cuò)誤
已經(jīng)有6人回復(fù)
- 【求助】fluent計(jì)算中，udf編譯通過(guò)，初始化出錯(cuò)
已經(jīng)有11人回復(fù)
jiangge2585
銅蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 75.5
- 帖子: 70
- 在線: 55.4小時(shí)
- 蟲(chóng)號(hào): 1508657
- 注冊(cè): 2011-11-24
- 性別: GG
- 專業(yè): 傳動(dòng)機(jī)械學(xué)
|
這些天想了想，好像是begin_f_loop(f,t) 和循環(huán)程序沖突了，所以才會(huì)卡主，begin_f_loop(f,t) 掃描出一個(gè)面網(wǎng)格中心P,然后把符合(sqrt(x*x+y*y)<=0.001) 的P丟進(jìn)循環(huán)語(yǔ)句執(zhí)行，執(zhí)行過(guò)一次之后，a和n的值就已經(jīng)超出了循環(huán)條件所規(guī)定的范圍，不能再進(jìn)行循環(huán)了。 我就改了一下，但下面這個(gè)程序卻沒(méi)有發(fā)熱，不知道哪有問(wèn)題了？求大神救我！ #include "udf.h" DEFINE_PROFILE(heat_flux,t,i) { real xx[ND_ND]; real y,x,r,a; int n; face_t f; begin_f_loop(f,t) { F_CENTROID(xx,f,t); y=xx[1]; x=xx[0]; r=0.001; if(sqrt(x*x+y*y)<=0.001) for(a=n=0;a>=0.1*n&&a<=0.1*n+0.008&&a<=1;n=(int)a/0.1) { F_PROFILE(f,t,i)=1160*0.7/(3.1415926*r*r)*exp(-(x*x+y*y)/(r*r)); a=RP_Get_Real("flow-time" ;} } end_f_loop(f,t) } |
2樓2013-06-25 15:14:35
CFD項(xiàng)目承接 
新蟲(chóng) (初入文壇)
|
本帖內(nèi)容被屏蔽 |
3樓2013-06-25 15:24:10
CFD項(xiàng)目承接
新蟲(chóng) (初入文壇)
★
1592203609: 金幣+1, 謝謝回帖，建議給出更為詳細(xì)的答案 2013-06-25 17:58:02
1592203609: 金幣+1, 謝謝回帖，建議給出更為詳細(xì)的答案 2013-06-25 17:58:02
|
本帖內(nèi)容被屏蔽 |
4樓2013-06-25 15:27:00
perfectgem
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1138.9
- 散金: 6
- 帖子: 92
- 在線: 11.7小時(shí)
- 蟲(chóng)號(hào): 2101103
- 注冊(cè): 2012-10-29
- 性別: MM
- 專業(yè): 零件加工制造
5樓2014-01-06 15:49:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 jiangge2585 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
