24小时热门版块排行榜

返回列表

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

buwangchuxin

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 12.3
帖子: 74
在线: 43.4小时
虫号: 1734871
注册: 2012-04-03
专业: 传热传质学

[求助] 70个金币求帮忙讲一个udf改写为并行的udf

下面是原来的udf，描述的是对壁面添加一个质量源项，由于迭代所需时间太长需要改为并行，哪位大侠帮帮忙吧！！！！多谢了

#include"udf.h"
#include"mem.h"
#define D 1.82e-5
#define v1 360
#define r 1.0E-5
DEFINE_ADJUST(adsorption_adjust, domain)
{
cell_t c0;
Thread* t0;
Thread* wall_thread;
face_t wall =6;
real A[ND_ND] ;
real x[ND_ND] ;
real At;
real J,C;
wall_thread= Lookup_Thread(domain,wall);
begin_f_loop (wall, wall_thread)
{
c0  = F_C0(wall,wall_thread);
t0 = F_C0_THREAD(wall,wall_thread);
C_CENTROID(x,c0,t0);
C=C_YI(c0,t0,0)*1000;
J=(-r*v1/4)/(1+r*v1/4*x[1]/D)*C;
F_AREA(A,wall, wall_thread);
At = NV_MAG(A);
C_UDMI(c0,t0,0)=At;
C_UDMI(c0,t0,1)=C_VOLUME(c0,t0);
C_UDMI(c0, t0,2) = J*At/C_VOLUME(c0,t0);
}
end_f_loop (wall,wall_thread)
}
DEFINE_SOURCE(mass_source,c,t,dS,eqn)
{
real source ;
source= C_UDMI(c,t,2) ;
dS[eqn]=0 ;
return source;
}

回复此楼

» 猜你喜欢

天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有11人回复
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有9人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有10人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有10人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复

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

fluent UDF代码调用预定义程序的问题（如fortran程序）已经有15人回复
udf求边界上的点出错已经有3人回复
求助udf的输入输出已经有16人回复
UDF 改为并行计算已经有9人回复
如何编写UDF来读取txt格式的一组数据？已经有11人回复
边界条件UDF求助已经有4人回复
UDF添加动能源项，发现温度没什么变化已经有20人回复
求助！！！我想要在fluent中针对一个面进行积分，那么UDF应该用什么宏函数？怎么用？已经有5人回复
好心人帮帮忙吧，为什么这个UDF导入fluent以后，在边界条件处没有显示~~~ 已经有17人回复
新手求救啊！Fluent非稳态计算怎么样才能在每一个时间步上调用一次UDF啊？已经有15人回复
求一个液滴碰壁的UDF（或相似UDF）已经有13人回复
激光焊接udf，加载温度场不对已经有7人回复
请教指定浓度的法向梯度为0的UDF如何写已经有5人回复
UDF源项的一个疑问已经有14人回复
udf profile update interval 是什么意思？已经有5人回复
求助：UDF模型的编译问题已经有15人回复
使用UDF后，初始化错误已经有11人回复
重金！关于udf的问题！已经有5人回复
UDF应用之后的初始化错误已经有6人回复
udf 重力方向变化已经有17人回复
介个UDF有啥问题已经有10人回复
【求助】单机多核fluent并行udf 已经有13人回复
【求助】求UDF高手帮忙看看我的问题已经有18人回复
【求助】50金币求帮写一个并行的UDF 已经有4人回复
【求助】物料属性可以用UDF修改吗已经有8人回复

1楼 2013-02-26 10:13:50

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

znon

金虫 (小有名气)

应助: 40 (小学生)
金币: 739.2
红花: 8
帖子: 164
在线: 74.7小时
虫号: 325496
注册: 2007-03-17
性别: GG
专业: 能源化工

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2013-02-28 11:34:09
buwangchuxin: 金币+100 2013-03-04 19:33:30

单纯从程序来看，不需要大改动。你看看下面的程序能否运行。

#include"udf.h"
#include"mem.h"
#define D 1.82e-5
#define v1 360
#define r 1.0E-5
DEFINE_ADJUST(adsorption_adjust, domain)
{
  #if !RP_HOST
      cell_t c0;
      Thread* t0;
      Thread* wall_thread;
      face_t wall =6;
real A[ND_ND] ;
real x[ND_ND] ;
      real At;
      real J,C;
      wall_thread= Lookup_Thread(domain,wall);
      begin_f_loop (wall, wall_thread)
      {
c0  = F_C0(wall,wall_thread);
t0 = F_C0_THREAD(wall,wall_thread);
C_CENTROID(x,c0,t0);
C=C_YI(c0,t0,0)*1000;
J=(-r*v1/4)/(1+r*v1/4*x[1]/D)*C;
F_AREA(A,wall, wall_thread);
At = NV_MAG(A);
C_UDMI(c0,t0,0)=At;
C_UDMI(c0,t0,1)=C_VOLUME(c0,t0);
C_UDMI(c0, t0,2) = J*At/C_VOLUME(c0,t0);
}
      end_f_loop (wall,wall_thread)
  #endif
}
DEFINE_SOURCE(mass_source,c,t,dS,eqn)
{
#if !RP_HOST
      real source ;
      source= C_UDMI(c,t,2) ;
      dS[eqn]=0 ;
      return source;
  #endif
}

赞一下(1人)

回复此楼

3楼2013-02-26 16:27:13

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

查看全部 3 个回答

buwangchuxin

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 12.3
帖子: 74
在线: 43.4小时
虫号: 1734871
注册: 2012-04-03
专业: 传热传质学

高手指教一下吧！！！！

回复此楼

2楼2013-02-26 13:46:07

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

查看全部 3 个回答