24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huanghuang36

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 138.6
帖子: 92
在线: 43小时
虫号: 2050668
注册: 2012-10-09
性别: MM
专业: 机械工程

[求助] 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;
}

刚注册的，目前只有两个金币，呵呵，不好意思啊~~

回复此楼

» 猜你喜欢

A期刊撤稿已经有3人回复
职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sunbird04

银虫 (小有名气)

应助: 36 (小学生)
金币: 710.2
红花: 4
帖子: 143
在线: 98.6小时
虫号: 881246
注册: 2009-10-23
性别: GG
专业: 钢铁冶金

引用回帖:

6楼: Originally posted by huanghuang36 at 2012-10-11 09:45:41
而且我并行运算的时候，它不是把一个模型分成几块分别给每个处理器进行计算嘛，但是它的分割很不规则啊，不是按我子域的形状分的，我运算的时候在分割的边界会有微度的变形产生，最后导致负体积，运算错误，这是怎 ...

这个，我没做过动网格。。。所以不能给你具体得回答。
Fluent 并行时会自动划分区域，当然你也可以自己手动划分区域：打开Fluet Parallel>Auto Partition 然后选择你所想要的划分方法（或者你选择Partition，完全自己手动划分区域），然后读取case文件。