huanghuang36

銅蟲 (小有名氣)

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

[求助] UDF 改為并行計算

我的UDF 在單機下能正常使用，但是并行不是需要改一下語句嘛，本人愚鈍，照著教程改了改但是一直不對，希望高人能幫忙指點一下，不勝感激啊
如下是我的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;
}

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

回復此樓

» 猜你喜歡

278求調劑已經有7人回復
26考研一志愿中國石油大學(華東)305分求調劑已經有5人回復
280求調劑已經有9人回復
279求調劑已經有5人回復
286求調劑已經有5人回復
材料與化工（0856）304求 B區(qū) 調劑已經有3人回復
268求調劑已經有8人回復
265求調劑已經有14人回復
22 350 本科985求調劑，求老登收留已經有3人回復
【考研調劑】化學專業(yè) 281分，一志愿四川大學，誠心求調劑已經有9人回復

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

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

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

Mr.Silence

新蟲 (初入文壇)

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

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

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

贊一下

回復此樓

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

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

huanghuang36

銅蟲 (小有名氣)

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

呵呵，目前我也不是很會啊，共同學習

贊一下

回復此樓

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

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

sunbird04

銀蟲 (小有名氣)

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

【答案】應助回帖

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

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

贊一下(1人)

回復此樓

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

huanghuang36

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

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

huanghuang36

銅蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

huanghuang36

銅蟲 (小有名氣)

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

送鮮花一朵

引用回帖:

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

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

贊一下(1人)

回復此樓

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

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

sunbird04

銀蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

watcher007

新蟲 (正式寫手)

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

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

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

贊一下

回復此樓

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

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

匿名

用戶注銷 (著名寫手)

Suhas V. Patankar

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

本帖僅樓主可見

贊一下

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

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

會飛的豬ss

銀蟲 (小有名氣)

應助: 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是一個源項，添加了你所說的語句之后，警告需要一個返回值，但是我明明寫了返回值，機器卻說需要一個返回值。有人指點說是因為語句的非正常使用導致函數(shù)體為空。請問，這個問題怎么解決？？？

贊一下

回復此樓

多使點勁，豬也能飛起來。。。

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

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

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

[考研] 【考研調劑】化學專業(yè) 281分，一志愿四川大學，誠心求調劑 +9	吃吃吃才有意義 2026-03-19	9/450	2026-03-21 13:28 by 鄒gv
[考研] 一志愿山大07化學 332分四六級已過本科山東雙非求調劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 265求調劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 307求調劑 +10	冷笙123 2026-03-17	10/500	2026-03-21 01:54 by JourneyLucky
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +4	葵梓衛(wèi)隊 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 一志愿蘇州大學材料求調劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 350求調劑 +5	weudhdk 2026-03-19	5/250	2026-03-20 22:04 by luoyongfeng
[考研] 0817 化學工程 299分求調劑有科研經歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 一志愿北京化工大學0703化學318分，有科研經歷，求調劑 +4	一瓶苯甲酸 2026-03-14	4/200	2026-03-20 20:36 by fen_rao
[考研] 一志愿吉林大學材料學碩321求調劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 0703化學調劑，六級已過，有科研經歷 +13	曦熙兮 2026-03-15	13/650	2026-03-20 19:35 by Dream007008
[考研] 0856調劑，是學校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 求調劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學校友
[考研] 生物學調劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 286求調劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 生物學071000 329分求調劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 一志愿南京大學，080500材料科學與工程，調劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 0856專碩279求調劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015