如何在库文件中修正中心差分的表面格式

[吴琼]

你好。想问一下在openfoam中离散过程中对于面插值的中心差分格式的源文件代码openfoam/src/finiteVolume/interpolation/surfaceInterpolation/surfaceInterpolationScheme/surfaceInterpolationScheme.C
其中对于点乘面差值格式的代码形式主要由Foam::surfaceInterpolationScheme<Type>::dotInterpolate 函数决定的，我现在的主要想法是想通过修正插值格式在代码sfi[fi] = Sfi[fi] & (lambda[fi]*(vfi[P[fi]] - vfi[N[fi]]) + vfi[N[fi]] ) 后面加体心梯度项点乘偏移量向量，如图所示：

• 我该如何修改代码，或者说大致的思路，请各位指点，拜托了。