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

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

24小時(shí)熱門版塊排行榜

返回列表

上一頁(yè)

sfbook

銀蟲 (小有名氣)

應(yīng)助: 13 (小學(xué)生)
金幣: 432.5
紅花: 2
帖子: 66
在線: 12.7小時(shí)
蟲號(hào): 4555158
注冊(cè): 2016-03-30
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
歪歪_crazy: 金幣+100, ★★★★★最佳答案 2016-12-01 17:05:41

引用回帖:

10樓: Originally posted by 歪歪_crazy at 2016-12-01 10:27:57
首先，我試驗(yàn)了一下，確實(shí)如你所說的，DEFINE_SOURCE 是自帶循環(huán)的，但是它的循環(huán)應(yīng)該是針對(duì)每一個(gè) cell ,也就是說沒有一個(gè) cell ,它會(huì)循環(huán)一次，而嵌在其中的 Loop 循環(huán)也會(huì)循環(huán)一次，并且每次循環(huán)得到的也是我先 ...

同學(xué)，我感覺你需要系統(tǒng)的學(xué)習(xí)一下C語(yǔ)言和UDF的編程了。。。。。
默認(rèn)單相流domain只有1，多相流才需要用到更多的domain，如果你的Phases選項(xiàng)沒開就不用考慮Get_Domain了，有問題可以仔細(xì)看看udf幫助里的Special Macros這節(jié)。我給你修改了一下程序，有問題再跟我說。
#include "udf.h"
// #define zone_ID 0 貌似沒有編號(hào)為0的區(qū)域吧。。。
double P_out=0.;

DEFINE_ADJUST(TEG_adjust,domain) //定義里面有 domain，不用自己定義
{
      Thread *t;
      cell_t c;
   //  Domain *all; 沒用，你的函數(shù)就是在定義的domain里用的
      int  Zone_ID=2; //兩個(gè)ID有一個(gè)就夠啦，不管你是全局變量還是局部變量，你又不會(huì)給他重新賦值，除非你要算多個(gè)區(qū)域的值,這里給個(gè) 2是隨意定的，而且墻壁編號(hào)通常為2

      face_t f;
      real temp=0.,grad=0.,source=0.;
      real a=0.,b=0.;
      real u=0.,u_sum=0.,vol=0.,vol_sum=4.752*0.00001;
      real r=0.,r_sum=0.,s=1.76*0.001;
      real p=0.;
t=Lookup_Thread(domain,Zone_ID);
            begin_c_loop (c,t)
            {
               temp=C_T(c,t);
               vol=C_VOLUME(c,t);
               grad=-0.0003*C_T_G(c,t)[0];/*0 for x,1 for y, 2 for z*/

               a=-6.035*0.0001+6.270*0.000001*temp-1.571*0.00000001*temp*temp+1.260*0.00000000001*temp*temp*temp;
               u=a*grad*vol/vol_sum;
               u_sum+=u;

               b=-7.188*0.0000001+6.926*0.000000001*temp-3.573*0.000000000001*temp*temp;
               r=b*vol/s/s;
               r_sum+=r;
            }
            end_c_loop (c,t) //開始和結(jié)束至少一樣啊。。。。

      p=u_sum*u_sum/(2*r_sum);

      P_out=p;//直接輸出p不好么。。。。

Message("%g, %g %g\n", u_sum,r_sum,P_out);
}

贊一下

回復(fù)此樓

11樓2016-12-01 15:33:16

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

歪歪_crazy

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 207
紅花: 1
帖子: 16
在線: 17.9小時(shí)
蟲號(hào): 2762998
注冊(cè): 2013-10-29
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

11樓: Originally posted by sfbook at 2016-12-01 15:33:16
同學(xué)，我感覺你需要系統(tǒng)的學(xué)習(xí)一下C語(yǔ)言和UDF的編程了。。。。。
默認(rèn)單相流domain只有1，多相流才需要用到更多的domain，如果你的Phases選項(xiàng)沒開就不用考慮Get_Domain了，有問題可以仔細(xì)看看udf幫助里的Special ...

讓你見笑了…… 上面給的程序比較亂是因?yàn)槲乙恢痹诟�，在不斷地試，所以有些定義就沒有刪去。

贊一下

回復(fù)此樓

12樓2016-12-01 17:08:58

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

歪歪_crazy

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 207
紅花: 1
帖子: 16
在線: 17.9小時(shí)
蟲號(hào): 2762998
注冊(cè): 2013-10-29
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

首先，十分感謝你的回復(fù)，給了我很大幫助。因?yàn)槭莿偨佑|UDF，并且C語(yǔ)言也只是會(huì)點(diǎn)皮毛，所以錯(cuò)誤百出，讓您見笑了。
其次我想分享一下我這幾天編這個(gè)udf的一些心得，算是班門弄斧，并且給以后看到這個(gè)帖子的人一些幫助吧。
1.關(guān)于C_T_G(c,t) 這個(gè)宏的運(yùn)用，需要在fluent的控制臺(tái)中打開一些設(shè)置 solve/set>expert Keep temporary solver memory from being free? [yes]
通過打開這個(gè)在設(shè)置來(lái)保證 Fluent 不會(huì)在計(jì)算結(jié)束就清除數(shù)據(jù)而導(dǎo)致的不能計(jì)算梯度。同時(shí)根據(jù)我的試驗(yàn)，發(fā)現(xiàn)在使用到這個(gè)語(yǔ)句時(shí)，如果在最初讀入，F(xiàn)luent沒有儲(chǔ)存值的情況下，也會(huì)報(bào)錯(cuò)，不能計(jì)算。這時(shí)就需要先不加載含有這個(gè)語(yǔ)句的udf,運(yùn)行幾步之后，在加載udf就能計(jì)算了。
2.在使用 Lookup_Thread 語(yǔ)句之后，不能使用 thread_loop_c,thread_loop_f 因?yàn)檫@兩者都是針對(duì) domain 全局進(jìn)行的循環(huán)，跟在 Lookup_Thread 這個(gè)指定區(qū)域語(yǔ)句后應(yīng)該是因?yàn)楫a(chǎn)生矛盾而報(bào)錯(cuò)。
3.還是關(guān)于 Lookup_Thread(domain,Zone_ID) 的，如果在本語(yǔ)句之后進(jìn)行的循環(huán)是關(guān)于 cell 的 thread_loop_c，那么查找的 ID 應(yīng)為 Cell Zone Conditions 中你想要計(jì)算區(qū)域的 ID，不行選擇 Boundary Conditions之中的，不然也會(huì)報(bào)錯(cuò)。
4.上傳一張關(guān)于各個(gè)函數(shù)調(diào)用位置的彩圖，感覺挺有用的，有興趣的可以看看。

最后，再次感謝 @sfbook ，謝謝你的耐心指導(dǎo)！如果以上我說的有么不對(duì)，還請(qǐng)指出。以后如果有機(jī)會(huì)，希望還能向你請(qǐng)教！

3333.PNG

1111.PNG

222.PNG

贊一下(1人)

回復(fù)此樓

13樓2016-12-01 18:52:14

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

歐式發(fā)酵乳

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 366.5
帖子: 187
在線: 25.5小時(shí)
蟲號(hào): 6243283
注冊(cè): 2017-04-05
性別: GG
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

送紅花一朵

看到你這個(gè)讓我想起來(lái)我起前兩天遇到的問題，跟你一模一樣，不夠我是通過在循環(huán)外賦值給全局變量才避免被清零的。

贊一下

回復(fù)此樓

14樓2020-10-15 10:48:04

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主學(xué)員fIJLZY 的主題更新

返回列表

上一頁(yè)

普通表情龍兔虎貓

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

[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +8	SUICHILD 2026-03-12	8/400	2026-03-18 20:55 by winsuccess
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 學(xué)碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

sfbook

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

歪歪_crazy

歪歪_crazy

歐式發(fā)酵乳

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频