| 查看: 641 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者1294219827将赠送您 30 个金币 | ||
1294219827新虫 (初入文坛)
|
[求助]
VOF模型，管道流动，我想求截面气相所占面积，下面是UDF
|
|
|
VOF模型，我想求截面气相所占面积，下面是UDF，有错误，不知道哪错了，请各位同仁指正！ DEFINE_ON_DEMAND(VOID) { real air_area=0; //存储气相所占面积 Thread *tf; Domain *d; face_t f; cell_t c0; //存储与f相邻的单元 Thread *tc0 = NULL; //存储与f相邻的单元线指针 d = Get_Domain(3); //得到气相域指针 tf= Lookup_Thread(d, 12); //获取截面所在的线指针 begin_f_loop(f, tf) { tc0 = Lookup_Thread(d, 26); c0 = F_C0(f, tf); if (C_VOF(c0, tc0) <= 1 && C_VOF(c0, tc0) >= 0.5) { F_AREA(A, f, tf); //获得面向量 air_area+ = NV_MAG(A); //求得面的面积 } } end_f_loop(f, t) printf("air_area=%g,avg=%g\n", air_area); } |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
10