您说雷诺数然后说网格影响粘度计算，我以为说的是雷诺数里面的粘度，那个液相粘度是个定值，所以我才把问题描述了一遍
我做了三种气泡三种网格，第一个帖子是其中一种气泡，表面张力大粘度大所以气泡为球形上升速度慢，采用M3网格（梯度自适应，隔几个时间步刷新一次气液界面网格）后上浮速度明显减慢了，您说的网格引起的数值扩散启发了我，让我大概明白这其中原因了（我算法不是很懂，不一定对）：
因为VOF模型是根据cell的相体积率加权算粘度和密度的，CSF表面张力模型会引起相体积率的数值扩散，进而影响粘度，这种数值扩散与气液界面网格密度直接相关，所以采用动态网格自适应（M3）会造成每一次对气泡运动前方的网格更新后，相体积率数值扩散，由于网格加密了，所以界面更加尖锐，也使得液相相体积率为1的网格更偏下，这会阻碍气泡上浮，使气泡上浮速度变慢，这种情况在液相粘度大的时候更明显，我其他两种气泡由于液相粘度比球形气泡低一个量级，就没出现这种情况。我最近更改M3刷新频率，发现网格更新频率越高，上浮速度越低，这也在预料之内，更新越频繁阻碍作用越显著。
其实网格越密，气泡上浮速度应该是越高的，但由于这种情况（液相高粘+气液界面网格自适应）改变了上浮速度，我其他两种气泡（液相粘度低），随着网格加密，气泡上浮速度增加。
谢谢您的解答

@pajack 首先你定义了一个Thread *t，此时它只是一个空指针，然后你就用t去找mixture_thread，它都不知道t在哪，怎么找？thread_loop_c(t,d)只会遍历d下面的t，而相体积率和温度分属不同的domain，建议你用mp_thread_loop_c

@bestucan 您好！此处的气泡上浮雷诺数：

分子是液相密度、气泡上浮速度和等效直径的乘积，分母是液相粘度，我是用气泡质心位移除以对应时间间隔算上浮速度。
也就是说用了M3网格后，气泡上浮速度明显变慢了，这是球形气泡，我做的其他两种形状的气泡都没有出现这种问题。

模拟三维单气泡上浮，选择了三种网格M1、M2、M3，监测气泡上浮的雷诺数，发现M1和M2结果基本吻合，M3采用梯度自适应，测得的雷诺数明显偏小了，有大佬知道原因吗？
M1:

M2:

M3:

结果：