版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

gyctju

金蟲 (正式寫手)

應(yīng)助: 314 (大學(xué)生)
金幣: 942.6
散金: 227
紅花: 48
沙發(fā): 3
帖子: 996
在線: 427.6小時
蟲號: 1569491
注冊: 2012-01-08
專業(yè): 分離過程

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai(金幣+2): 謝謝應(yīng)助！ 2012-01-08 22:58:47

可能你想法錯了。邊界條件部分不應(yīng)該跟累加傳熱量的循環(huán)放在一起，這樣的實際效果跟你想的是不一樣的。
還有邊界條件設(shè)定時應(yīng)該用begin-f-loop
建議修改成如下形式試試
#include "udf.h"
real Q_tot=0;
DEFINE_PROFILE (unsteady_heatflux, thread, position)
{
Domain *d;
real cp=880;
real density=2180;
real tempn, templ, volume, Q;
Thread *t;
cell_t c;
d= Get_Domain(1); /*Get the domain using Fluent utility *//*Loop over all cell threads in the domain*/
real time= RP_Get_Real("flow-time";
real b=(int)(time/3600)+1;
int i=(int)((b-1)/24);
real Heat[3]={278.5, 280.7, 278.5};
thread_loop_c(t,d) /*Compute Q */ /*Loop over all cells */
begin_c_loop(c,t)
{

volume=C_VOLUME(c,t); /* get cell volume */
templ=C_T_M1(c,t); /*Get cell tempertuare of previous step*/
tempn=C_T(c,t); /*Get cell tempertuare*/
Q=cp*density*volume*(tempn-templ);
Q_tot +=Q;
}
end_c_loop(c,t)
}
begine_f_loop(c,t,)
{
if (Q_tot<= Heat )
F_PROFILE(c,t,position)=60;
else F_PROFILE(c,t,position)=0;
}
end_f_loop

贊一下(1人)

回復(fù)此樓

11樓2012-01-08 21:11:19

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

xiaoguai3

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 41
帖子: 23
在線: 8.4小時
蟲號: 1543104
注冊: 2011-12-19
專業(yè): 可再生與替代能源利用中的

引用回帖:

樓: Originally posted by gyctju at 2012-01-08 21:11:19:
可能你想法錯了。邊界條件部分不應(yīng)該跟累加傳熱量的循環(huán)放在一起，這樣的實際效果跟你想的是不一樣的。
還有邊界條件設(shè)定時應(yīng)該用begin-f-loop
建議修改成如下形式試試
#include "udf.h"
real Q_t ...

謝謝回復(fù)~請問累加傳熱量該怎么處理呢？

贊一下

回復(fù)此樓

12樓2012-01-09 09:29:51

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

gyctju

金蟲 (正式寫手)

應(yīng)助: 314 (大學(xué)生)
金幣: 942.6
散金: 227
紅花: 48
沙發(fā): 3
帖子: 996
在線: 427.6小時
蟲號: 1569491
注冊: 2012-01-08
專業(yè): 分離過程

★
xiegangmai(金幣+1): 謝謝參與！ 2012-01-09 21:27:48

第一個循環(huán)就是累加熱量的

贊一下(1人)

回復(fù)此樓

13樓2012-01-09 10:42:57

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

gyctju

金蟲 (正式寫手)

應(yīng)助: 314 (大學(xué)生)
金幣: 942.6
散金: 227
紅花: 48
沙發(fā): 3
帖子: 996
在線: 427.6小時
蟲號: 1569491
注冊: 2012-01-08
專業(yè): 分離過程

引用回帖:

12樓: Originally posted by xiaoguai3 at 2012-01-09 09:29:51:
謝謝回復(fù)~請問累加傳熱量該怎么處理呢？

我上邊回復(fù)的那個修改大括號那塊出錯了，把begine_f_loop前面的大括號挪到最后去

贊一下

回復(fù)此樓

14樓2012-01-09 10:45:24

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

xiaoguai3

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 41
帖子: 23
在線: 8.4小時
蟲號: 1543104
注冊: 2011-12-19
專業(yè): 可再生與替代能源利用中的

引用回帖:

樓: Originally posted by gyctju at 2012-01-09 10:45:24:
我上邊回復(fù)的那個修改大括號那塊出錯了，把begine_f_loop前面的大括號挪到最后去

你好，我累加放在邊界里面是因為需要加一個判斷語句，用到REAL_TIME等參數(shù)，累加到一定情況Q歸零重新開始加。

贊一下

回復(fù)此樓

15樓2012-01-09 11:26:19

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

gyctju

金蟲 (正式寫手)

應(yīng)助: 314 (大學(xué)生)
金幣: 942.6
散金: 227
紅花: 48
沙發(fā): 3
帖子: 996
在線: 427.6小時
蟲號: 1569491
注冊: 2012-01-08
專業(yè): 分離過程

★
xiegangmai(金幣+1): 謝謝參與！ 2012-01-09 21:28:09

again，邊界條件用begin_f_loop
估計你都沒看我修改后的語句。

贊一下(1人)

回復(fù)此樓

16樓2012-01-09 19:16:54

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

xiaoguai3

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 41
帖子: 23
在線: 8.4小時
蟲號: 1543104
注冊: 2011-12-19
專業(yè): 可再生與替代能源利用中的

引用回帖:

16樓: Originally posted by gyctju at 2012-01-09 19:16:54:
again，邊界條件用begin_f_loop
估計你都沒看我修改后的語句。

我在自己程序中改了。試圖用fprintf將計算的Q值存入文件，但運行后發(fā)現(xiàn)TXT文件中數(shù)字一直是零，不知道哪里出錯了。

CODE:

#include "udf.h"

#define cp 800.

#define density 2180.

FILE *fp;

DEFINE_PROFILE(unsteady_temp,t,index)

{

float tempn,templ,volume,Q,Q_tot,time,a,b,tempad;

int i;

float Heat[5]={2780000.,        280.7,        278.5,        280.5,        285.};  

Domain *d;

Thread *ct,*t0;

cell_t c,c0;

face_t f;

time = RP_Get_Real("flow-time");

a = time/3600;

b = a/24;

d = Get_Domain(1);

i=(int)(((int)a)/24);

if (b==(int)b)   

Q_tot=0;

else 

{

fp=fopen("Q_tot.txt","r");

fscanf(fp,"%f",Q_tot);

fclose(fp);

}

thread_loop_c(ct,d)  

{

        begin_c_loop(c,ct) /*Loop over all cells循環(huán)所有THREAD上所有CELL */

       {

volume=C_VOLUME(c,ct); /* get cell volume */

templ=C_T_M1(c,ct); /*Get cell tempertuare of previous step有問題�。。。。。。�！*/

tempn=C_T(c,ct);   /*Get cell tempertuare*/

Q=cp*density*volume*(tempn-templ);

Q_tot +=Q;

       }

     end_c_loop(c,ct)

}

fp=fopen("Q_tot.txt","w");

  fprintf(fp,"%f",Q_tot);    /*寫入累加后Q_tot值 */

fclose(fp);

if (Q_tot<= Heat[i])

  begin_f_loop(f,t)

{   

    F_PROFILE(f,t,index)=333;  /*BC:Temp=333K*/

}

end_f_loop(f,t) 

else

begin_f_loop(f,t)

{   

c0 = F_C0(f,t);

t0 = THREAD_T0(t);

tempad=C_T(c0,t0); /*temperature of adjacent cell*/

F_PROFILE(f,t,index)=tempad; 

}

end_f_loop(f,t)

}

贊一下

回復(fù)此樓

17樓2012-01-09 21:21:55

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 xiaoguai3 的主題更新

返回列表

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

[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	10/500	2026-03-18 20:41 by Wangjingyue
[考研] 331求調(diào)劑（0703有機化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 機械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky