CFD中文网

在做液滴撞击的模拟，因为要用到动态接触角所以需要获取三相点位置以及速度。目前是想用两步迭代的位移差来计算速度，但是现在不知道怎么去算这个三相点的位置。代码是这么写的

    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变成面积啊

具体工况是我按照我自己试验做的设置的，水滴是从10cm位置自然下落，结果正常来说铺展过程应该都是圆形的，收缩的时候可能会因为材料表面不均匀变成不规则形状。但是这模拟不应该存在因为平面不均匀的问题吧

@李东岳
哎这网格我也不知道为啥，之前设置整个流域是方形的得到结果是方的，然后我给改成圆柱形的结果还是方形的，也不知道到底是哪部分的问题。

@bestucan 感谢您的回复！
我仔细看了下您发的这个实验，实验是通过拉伸薄膜的改变液滴应力分布的。按照我设置的条件来看的话应该是不太一样，因为底面我设置的是60度接触角的壁面。但是不排除有您说的这种可能，所以我先试试把侧边和上边都改成wall或者其他边界条件试试。
感谢！

@bestucan 那我重新整理下问题
模型：VOF模型
相：空气、水（表面张力0.0728）
计算域以及边界条件：直径40mm高6mm圆柱，在中心处加密加密尺度约为0.03。上方压力进
口，侧方压力出口，底面接触角60度
离散格式：PISO其他默认
水滴设置：半径r=2.4mm/1.8mm 贴紧壁面-1.387速度下落
时间步长：10-6 / 5.10-7 / 10-7三种都试过每步迭代20次

迭代一万步后水滴接近完全展开的时候就会变成下边这种方块的形状，正常应该是圆形的

各位老师同学，我在模拟三维水滴的时候不知道为什么铺展过程水滴会变成方形的。网格尺度换过很多次，计算域边界形状方形圆形都用过，但都会出现这个情况，时间尺度-6 -7结果也有可能是这样，这是为啥啊