[fxzf]

Hello OF experters.

I want to ask a silly question, I am trying to understand how OF implement processor boundary condition.

Let's use the laplacian as example, so in build laplacian scheme, there is

fvmLaplacianUncorrected
{
....

  const fvsPatchScalarField& pDeltaCoeffs = deltaCoeffs.boundaryField()[patchi];

  if (pvf.coupled())
  {
    fvm.internalCoeffs()[patchi] = pGamma*pvf.gradientInternalCoeffs(pDeltaCoeffs)
    fvm.boundaryCoeffs()[patchi] = -pGamma*pvf.gradientBoundaryCoeffs(pDeltaCoeffs)

  }

}

First question, I believer, if it loops to processor type boundary patch, it will this pvf.coupled() is true? Correct me if I am wrong.

Second question, what is the pDeltaCoeffs here? I know deltaCoeffs here is tsnGradScheme_().deltaCoeffs(vf), so I believe it should be 1/distance(cell centre to neighbour cell centre) on the field, and on the processor boundary patch, it is 1/distance(cell centre to patch face centre) ???

However, for processor boundary patch, to build a matrix, I believe it should get the distance from cell centre to neighbour cell centre ???

Third Question, I checked processorFvPatchField.C, I didn't find gradientBoundaryCoeffs and gradientInternalCoeffs. So where and how do they implemented for processsor B.C.

Finally, for pGamma, I believe, it should be interpolated value from this processor cell to neighbour processor cell ? but I cannot find where does this interpolation implemented.

Thanks very much.