|
fluent中有两类face:
interior face
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】
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-19 04:18
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社