在使用interFoam的时候，偶尔会需要计算液体质心的位置，以便于分析液体的载荷转移量，下面记录一种比较笨的实现方法。
在interFoam求解器源代码中添加以下代码，可以计算液体质心：
注：如果把代码中密度rho去掉，那么计算的是几何体模型的质心。
之后可以通过Info语句将质心数据输出于log文档中，最后用matlab提取log中每一个时间步中关于质心的信息，在matlab中可以使用下面代码来提取质心信息：
另外有一种比较好的方法，当然就是把它写成类似于forces，probs的后处理函数，用这种方法可以直接输出仅含质心信息的数据文档，有空在记录下。
来源：oschina
链接：https://my.oschina.net/LLLenhomepage/blog/4743584