fluent中有两类face：

1. interior face
   

2. boundary face

fluent自带的face macro能对face进行一些简单的操作【link1】。

其中，针对boundary face的macro多一些，

而针对interior face的macro，只有F_P(f,t)F_P(f,t)P(f,t)  和  F_FLUX(f,t)F_P(f,t)

F_P(f,t)  和  F_FLUX(f,t)

但是跟cell macro相比【link2】，有点少。

可是，有的时候，

还真需要确定某个 interior face上的物理量。

比如根据速度判断是否倒流？

此时，需要借助cell macro。

基本思路是，

把interior face两侧的cell上的物理量做平均。

因为interior face是2个cell的结合面，

而boundary face只有1个cell。

先定义好两侧cell

cell_t  c0,  c1 = -1;
Thread *t0, *t1 = NULL;

锁定一侧

 c0 = F_C0(f,t);

 t0 = F_C0_THREAD(f,t);

锁定另一侧

c1 = F_C1(f,t);       /* Get cell on other side of face */

t1 = F_C1_THREAD(f,t);

对两侧做平均

vel_x=(C_U(c0,t0)+C_U(c1,t1))/2;

如果是多相流，

需要在上面乘以组分。

vel_x=(C_U(c0,t0)*C_VOF(c0,t0)+C_U(c1,t1)*C_VOF(c1,t1))/2;

下面的例子，

是对某个interior face计算质量流量

【link3】

如果是并行的话，

还需要加上并行保护。

参考如下例子，

【link4】