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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

huanghuang36

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 138.6
帖子: 92
在線: 43小時
蟲號: 2050668
注冊: 2012-10-09
性別: MM
專業(yè): 機(jī)械工程

[求助] UDF 改為并行計算

我的UDF 在單機(jī)下能正常使用，但是并行不是需要改一下語句嘛，本人愚鈍，照著教程改了改但是一直不對，希望高人能幫忙指點(diǎn)一下，不勝感激啊
如下是我的UDF：主要使用的就是CG_MOTION

#include
#include "udf.h"
#include "math.h"
real r_movefluid;
DEFINE_CG_MOTION(moveface, dt, cg_vel, cg_omega, time, dtime)
{
      Thread *t;
face_t f;
      real x0,y0,r;
real v0;
      real v;
real omega=M_PI;
      real theta=omega*time;
real x[ND_ND];
      t = DT_THREAD(dt);
begin_f_loop(f,t)
      {
            F_CENTROID(x,f,t);
            x0=x[0];
            y0=x[1];
            r=sqrt(pow(x0,2)+pow(y0,2));
      }
end_f_loop(f,t)
      v0=omega*r;
      r_movefluid=r;
if (time>=0 && time<=M_PI/(180*omega))
         v=0;
else if (time<=8*M_PI/(180*omega))
         v=12.96*(theta-(M_PI/180))*omega/(7*pow(M_PI,2));
else if (time<=10*M_PI/(180*omega))
         v=0.072*omega/M_PI;
else if (time<=29*M_PI/(180*omega))
         v=(0.072*omega/M_PI)-12.96*(theta-M_PI/18)*omega/(19*M_PI*M_PI);
else if (time<=31*M_PI/(180*omega))
         v=0;

      cg_vel[0]=v*cos(theta)-v0*sin(theta);
      cg_vel[1]=v*sin(theta)+v0*cos(theta);
      cg_vel[2]=0;
      cg_omega[0]=0;
      cg_omega[1]=0;
      cg_omega[2]=omega;
}
DEFINE_CG_MOTION(rotate, dt, cg_vel, cg_omega, time, dtime)
{
      real omega=M_PI;
      cg_omega[2]=omega;
}
DEFINE_CG_MOTION(movefluid, dt, cg_vel, cg_omega, time, dtime)
{
real v0;
      real v;
real omega=M_PI;
      real theta=omega*time;
      v0=omega*r_movefluid;
if (time>=0 && time<=M_PI/(180*omega))
         v=0;
else if (time<=8*M_PI/(180*omega))
         v=12.96*(theta-(M_PI/180))*omega/(7*pow(M_PI,2));
else if (time<=10*M_PI/(180*omega))
         v=0.072*omega/M_PI;
else if (time<=29*M_PI/(180*omega))
         v=(0.072*omega/M_PI)-12.96*(theta-M_PI/18)*omega/(19*M_PI*M_PI);
else if (time<=31*M_PI/(180*omega))
         v=0;

      cg_vel[0]=v*cos(theta)-v0*sin(theta);
      cg_vel[1]=v*sin(theta)+v0*cos(theta);
      cg_vel[2]=0;
      cg_omega[0]=0;
      cg_omega[1]=0;
      cg_omega[2]=omega;
}

剛注冊的，目前只有兩個金幣，呵呵，不好意思啊~~

回復(fù)此樓

» 猜你喜歡

【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有11人回復(fù)
一志愿深大，0703化學(xué)，總分302，求調(diào)劑已經(jīng)有4人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有3人回復(fù)
材料求調(diào)劑已經(jīng)有3人回復(fù)
333求調(diào)劑已經(jīng)有5人回復(fù)
求助已經(jīng)有6人回復(fù)
278求調(diào)劑已經(jīng)有9人回復(fù)
297求調(diào)劑已經(jīng)有12人回復(fù)
0703化學(xué)297求調(diào)劑已經(jīng)有3人回復(fù)
材料學(xué)碩301分求調(diào)劑已經(jīng)有6人回復(fù)

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

VASP并行計算出錯已經(jīng)有9人回復(fù)
100個金幣：udf改成并行的，求高人已經(jīng)有22人回復(fù)
曳力模型udf要并行計算，該怎么修改呢已經(jīng)有5人回復(fù)
并行計算如何實(shí)現(xiàn)k點(diǎn)平均分配？已經(jīng)有5人回復(fù)
【求助】用vof模型計算出的界面，有木有udf語句可以找得到已經(jīng)有6人回復(fù)
【求助】單機(jī)多核fluent并行udf 已經(jīng)有13人回復(fù)
【求助】單機(jī)多核UDF怎么并行計算已經(jīng)有15人回復(fù)
【求助】fluent計算中，udf編譯通過，初始化出錯已經(jīng)有11人回復(fù)
【求助】linux版fluent使用udf并行問題請教已經(jīng)有4人回復(fù)
【求助】castep并行計算問題已經(jīng)有15人回復(fù)

1樓 2012-10-09 10:43:48

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

Mr.Silence

新蟲 (初入文壇)

應(yīng)助: 4 (幼兒園)
金幣: 105.3
帖子: 34
在線: 5.3小時
蟲號: 2023898
注冊: 2012-09-23
專業(yè): 流體力學(xué)

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣-1, 應(yīng)助指數(shù)-1, 謝謝參與，非應(yīng)助請勿選擇“應(yīng)助回帖” 2012-10-10 21:19:51

同求并行計算的方法。我也在學(xué)習(xí)UDF，希望樓主不吝賜教。

贊一下

回復(fù)此樓

2樓2012-10-09 11:17:09

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

huanghuang36

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 138.6
帖子: 92
在線: 43小時
蟲號: 2050668
注冊: 2012-10-09
性別: MM
專業(yè): 機(jī)械工程

呵呵，目前我也不是很會啊，共同學(xué)習(xí)

贊一下

回復(fù)此樓

3樓2012-10-10 10:02:36

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

sunbird04

銀蟲 (小有名氣)

應(yīng)助: 36 (小學(xué)生)
金幣: 710.2
紅花: 4
帖子: 143
在線: 98.6小時
蟲號: 881246
注冊: 2009-10-23
性別: GG
專業(yè): 鋼鐵冶金

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

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝應(yīng)助 2012-10-10 21:20:01
huanghuang36: 金幣+2, ★★★很有幫助, 謝謝你的幫助哦~~ 2012-10-16 09:34:24

在每段程序最前頭加上 #if !RP_HOST 最末尾加上#endif 即可

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

huanghuang36

4樓2012-10-10 18:14:27

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

huanghuang36

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 138.6
帖子: 92
在線: 43小時
蟲號: 2050668
注冊: 2012-10-09
性別: MM
專業(yè): 機(jī)械工程

引用回帖:

4樓: Originally posted by sunbird04 at 2012-10-10 18:14:27
在每段程序最前頭加上 #if !RP_HOST 最末尾加上#endif 即可

太感謝你的回答了，我已經(jīng)照你說的改了程序，還有我想問一下，那個計算節(jié)點(diǎn)（也就是處理器）數(shù)量的設(shè)置是不是和我的子域的數(shù)量有關(guān)啊，我看書上說子域的數(shù)量應(yīng)該是計算節(jié)點(diǎn)的整數(shù)倍，是要這么設(shè)是吧~~

贊一下

回復(fù)此樓

5樓2012-10-10 21:40:46

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

huanghuang36

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 138.6
帖子: 92
在線: 43小時
蟲號: 2050668
注冊: 2012-10-09
性別: MM
專業(yè): 機(jī)械工程

送鮮花一朵

引用回帖:

4樓: Originally posted by sunbird04 at 2012-10-10 18:14:27
在每段程序最前頭加上 #if !RP_HOST 最末尾加上#endif 即可

而且我并行運(yùn)算的時候，它不是把一個模型分成幾塊分別給每個處理器進(jìn)行計算嘛，但是它的分割很不規(guī)則啊，不是按我子域的形狀分的，我運(yùn)算的時候在分割的邊界會有微度的變形產(chǎn)生，最后導(dǎo)致負(fù)體積，運(yùn)算錯誤，這是怎么回事？？呵呵，求指點(diǎn)~~！~~~~~

贊一下(1人)

回復(fù)此樓

6樓2012-10-11 09:45:41

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

sunbird04

銀蟲 (小有名氣)

應(yīng)助: 36 (小學(xué)生)
金幣: 710.2
紅花: 4
帖子: 143
在線: 98.6小時
蟲號: 881246
注冊: 2009-10-23
性別: GG
專業(yè): 鋼鐵冶金

引用回帖:

6樓: Originally posted by huanghuang36 at 2012-10-11 09:45:41
而且我并行運(yùn)算的時候，它不是把一個模型分成幾塊分別給每個處理器進(jìn)行計算嘛，但是它的分割很不規(guī)則啊，不是按我子域的形狀分的，我運(yùn)算的時候在分割的邊界會有微度的變形產(chǎn)生，最后導(dǎo)致負(fù)體積，運(yùn)算錯誤，這是怎 ...

這個，我沒做過動網(wǎng)格。。。所以不能給你具體得回答。
Fluent 并行時會自動劃分區(qū)域，當(dāng)然你也可以自己手動劃分區(qū)域：打開Fluet Parallel>Auto Partition 然后選擇你所想要的劃分方法（或者你選擇Partition，完全自己手動劃分區(qū)域），然后讀取case文件。

贊一下

回復(fù)此樓

7樓2012-10-17 17:18:35

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

watcher007

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6.9
帖子: 316
在線: 122.7小時
蟲號: 1808594
注冊: 2012-05-10
專業(yè): 機(jī)械動力學(xué)

★
xiegangmai: 金幣+1, 鼓勵討論交流 2012-10-21 17:24:27

樓主，看了你的UDF，感覺我們做的東西很像，有機(jī)會可以進(jìn)一步交流一下。我曾經(jīng)用在Fluent并行運(yùn)算，選擇使用自動分區(qū)，一般感覺是按照控制體進(jìn)行簡單的劃分，例如長方體的控制體，可能就是分成兩半，或四分之一，另外在并行使一般核數(shù)都是選偶數(shù)的。

贊一下

回復(fù)此樓

8樓2012-10-17 18:54:28

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

匿名

用戶注銷 (著名寫手)

Suhas V. Patankar

應(yīng)助: 165 (高中生)
金幣: 439.4
散金: 3001
紅花: 10
沙發(fā): 1
帖子: 1878
在線: 2074.9小時
蟲號: 0
注冊: 2012-05-12
性別: GG
專業(yè): 流體力學(xué)

本帖僅樓主可見

贊一下

9樓2013-07-13 18:54:30

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

會飛的豬ss

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 465.2
帖子: 79
在線: 30.9小時
蟲號: 3011673
注冊: 2014-03-03
性別: GG
專業(yè): 腫瘤物理治療

引用回帖:

4樓: Originally posted by sunbird04 at 2012-10-10 18:14:27
在每段程序最前頭加上 #if !RP_HOST 最末尾加上#endif 即可

試了這個辦法，我的UDF是一個源項(xiàng)，添加了你所說的語句之后，警告需要一個返回值，但是我明明寫了返回值，機(jī)器卻說需要一個返回值。有人指點(diǎn)說是因?yàn)檎Z句的非正常使用導(dǎo)致函數(shù)體為空。請問，這個問題怎么解決？？？

贊一下

回復(fù)此樓

多使點(diǎn)勁，豬也能飛起來。。。

10樓2014-07-23 11:32:45

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

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

返回列表

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

[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 材料求調(diào)劑 +3	@taotao 2026-03-21	3/150	2026-03-21 18:11 by WTUChen
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 302求調(diào)劑 +12	呼呼呼。。。。 2026-03-17	12/600	2026-03-21 17:29 by ColorlessPI
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	7/350	2026-03-21 16:09 by lemonzzn
[考研] 0856材料專碩353求調(diào)劑 +3	NIFFFfff 2026-03-20	3/150	2026-03-21 10:23 by luoyongfeng
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 310求調(diào)劑 +3	baibai1314 2026-03-16	3/150	2026-03-21 03:56 by JourneyLucky
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8	niko- 2026-03-19	9/450	2026-03-20 21:57 by luoyongfeng
[考研] 本人考085602 化學(xué)工程專碩 +19	不知道叫什么！ 2026-03-15	21/1050	2026-03-20 20:48 by zhukairuo
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 0703化學(xué) 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 304求調(diào)劑 +5	素年祭語 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海

24小時熱門版塊排行榜

[求助] UDF 改為并行計算

» 猜你喜歡

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

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

» 本帖已獲得的紅花（最新10朵）

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