| 查看: 628 | 回复: 2 | ||
[求助]
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; } |
回复此楼» 猜你喜欢
301求调剂
已经有18人回复
-
化工学硕 285求调剂
已经有19人回复
-
328求调剂
已经有8人回复
-
求调剂
已经有28人回复
-
304求调剂（085602，过四级，一志愿985）
已经有21人回复
-
288环境专硕,求调材料方向
已经有33人回复
-
285求调剂
已经有16人回复
-
273求调剂
已经有41人回复
-
331求调剂
已经有4人回复
-
298求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
2楼2013-02-26 13:46:07
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢参与！ 2013-02-28 11:34:09
buwangchuxin: 金币+100 2013-03-04 19:33:30
感谢参与，应助指数 +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 } |
3楼2013-02-26 16:27:13
