在做液滴撞击的模拟，因为要用到动态接触角所以需要获取三相点位置以及速度。目前是想用两步迭代的位移差来计算速度，但是现在不知道怎么去算这个三相点的位置。代码是这么写的
Thread* thread = Lookup_Thread(domain, 9);
Thread** pt = THREAD_SUB_THREADS(thread);
cell_t cell;
face_t f;
real x[ND_ND], area[ND_ND];
double max_x = 0, f_Hoff_inverse, x_hoff, temp, Ca, volume;
int n;
sum = 0.0;
begin_c_loop_all(cell, pt[1])
{
if (C_VOF(cell, pt[1]) != 0)
{
C_CENTROID(x, cell, pt[1]);
if (x[0] > max_x)
max_x = x[0];
}
sum += C_VOF(cell, pt[1]);
}
end_c_loop_all(cell, pt[1])
我想问一下怎么把SUM变成面积啊
