| 查看: 410 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
fanerye2015新虫 (初入文坛)
|
[求助]
DEFINE_ADJUST编译总是报错
|
|
| 我的模型里面有气液两相，用DEFINE_ADJUST调整单元材料属性，将在一定范围内的单元赋予液体属性，代码如下 |
回复此楼» 猜你喜欢
求调剂
已经有3人回复
-
一志愿北京化工大学，初试成绩350求调剂
已经有14人回复
-
复试调剂
已经有4人回复
-
323分（计算机视觉和大模型项目）能直接上手
已经有3人回复
-
311分 22408 求调剂
已经有3人回复
-
320分人工智能调剂
已经有7人回复
-
一志愿郑大0705求调剂
已经有4人回复
-
0703化学
已经有10人回复
-
301求调剂
已经有10人回复
-
306分材料与化工求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
fanerye2015
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160.5
- 散金: 7
- 红花: 1
- 帖子: 42
- 在线: 24小时
- 虫号: 3716536
- 注册: 2015-03-07
- 专业: 金属材料的凝固与结晶学
4楼2015-12-16 09:55:25
fanerye2015
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160.5
- 散金: 7
- 红花: 1
- 帖子: 42
- 在线: 24小时
- 虫号: 3716536
- 注册: 2015-03-07
- 专业: 金属材料的凝固与结晶学
|
#include "udf.h" #include "sg_mphase.h" #include "math.h" #define v2 0.033333 #define r2 0.0004 DEFINE_ADJUST(powder,domain) { Thread *t; cell_t c; Thread *pri_th, *sec_th; real xx[ND_ND]; real x2,y2,z2,r; real Time2=RP_Get_Real("flow_time"  thread_loop_c(t,domain) { pri_th=THREAD_SUB_THREAD(t,0); sec_th=THREAD_SUB_THREAD(t,1); begin_c_loop_all(c,t) { C_CENTROID(xx,c,t); x2=xx[0]-0.003-v2*Time2; y2=xx[1]; z2=xx[2]-0.0005; r=sqrt(x2*x2+y2*y2+z2*z2); if(r<=r2) { C_VOF(c,sec_th)=1; C_T(c,t)=1928; } } end_c_loop_all(c,t); } } |
2楼2015-12-16 09:54:18
fanerye2015
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 160.5
- 散金: 7
- 红花: 1
- 帖子: 42
- 在线: 24小时
- 虫号: 3716536
- 注册: 2015-03-07
- 专业: 金属材料的凝固与结晶学
|
#include "udf.h" /*包括常规宏*/ #include "sg_mphase.h" /*包括体积分数宏*/ #include "math.h" #define v2 0.033333 #define r2 0.0004 DEFINE_ADJUST(powder,domain) { Thread *t; cell_t c; Thread *pri_th, *sec_th; real xx[ND_ND]; real x2,y2,z2,r; real Time2=RP_Get_Real("flow_time" thread_loop_c(t,domain) { pri_th=THREAD_SUB_THREAD(t,0); sec_th=THREAD_SUB_THREAD(t,1); begin_c_loop_all(c,t) { C_CENTROID(xx,c,t); x2=xx[0]-0.003-v2*Time2; y2=xx[1]; z2=xx[2]-0.0005; r=sqrt(x2*x2+y2*y2+z2*z2); if(r<=r2) { C_VOF(c,sec_th)=1; C_T(c,t)=1928; } } end_c_loop_all(c,t); } } |
3楼2015-12-16 09:54:54
10