并行边界上网格面心处数值通量(φ_f)的同步

[新喻庸]

最近在编写程序的时候遇到并行边界上网格面心处数值通量(φ_f)的同步问题，单核运行不存在并行边界，所以没有这个问题。以下图为例：

并行计算时，求解区域分为了4个区域，两个区域之间产生了并行边界面，对于内部网格面元，网格面心处的数值通量(φ_f)的计算跟单核运行时是一样的，但是对于并行边界处的网格面元，数值通量(φ_f)在两个相邻并行边界上应该是一样的，因为并行边界的存在不应该改变原来的解。因此，我用了如下的syncFaceList()函数来进行同步，alpha_f是定义在网格面心上的surfaceScalarField，这个函数在头文件"syncTools.H"中有定义，我用了以下三种同步方式，发现流体穿过并行边界时还是产生了误差，单核运行时误差(数值解与解析解之差)可以达到10E-12次方量级，但是并行运行时流体一穿过并行边界，误差就上升到了10E-4量级，我觉得这是因为并行边界上数据同步没有做好的原因，请问有小伙伴遇到过这个问题嘛？使用syncFaceList()函数的正确方式是什么？如果使用syncFaceList()函数进行同步网格面心处的数值通量(φ_f)还是有问题的话，打算自己写一下同步过程。

// syncTools::syncFaceList(mesh(), alpha_f, maxEqOp<scalar>(), false);
syncTools::syncFaceList(mesh(), alpha_f, plusEqOp<scalar>(), false);
// syncTools::syncFaceList(mesh(), alpha_f, eqOp<scalar>(), false);