62 mesh.pointsInstance(),
68 auto& oDelta = *oDelta_;
75 mesh.pointsInstance(),
81 auto& nDelta = *nDelta_;
96 mag(
n[facei] & (
C[owner[facei]] - Cf[facei]));
98 mag(
n[facei] & (
C[neighbour[facei]] - Cf[facei]));
105 const fvPatch& currPatch =
mesh.boundary()[patchi];
108 const vectorField nPatch(currPatch.Sf()/currPatch.magSf());
111 if (currPatch.coupled())
118 const label own = pOwner[facei];
121 oDelta.boundaryFieldRef()[patchi][facei] =
122 mag(nPatch[facei] & (pCf[facei] -
C[own]));
126 nDelta.boundaryFieldRef()[patchi] =
127 currPatch.weights()*oDelta.boundaryFieldRef()[patchi]
128 /(scalar(1) - currPatch.weights());
137 const label own = pOwner[facei];
140 oDelta.boundaryFieldRef()[patchi][facei] =
141 mag(nPatch[facei] & (pCf[facei] -
C[own]));
143 nDelta.boundaryFieldRef()[patchi][facei] =
144 mag(nPatch[facei] & (pCf[facei] -
C[own]));
166 "weightedFlux::interpolate(" + vf.name() +
')',
167 mesh.time().timeName(),
173 auto& result = tresult.ref();
180 const scalar sigmaDeltaO = sigma_[owner[facei]]/oDelta[facei];
181 const scalar sigmaDeltaN = sigma_[neighbour[facei]]/nDelta[facei];
184 (vf[owner[facei]]*sigmaDeltaO + vf[neighbour[facei]]*sigmaDeltaN)
185 /(sigmaDeltaO + sigmaDeltaN);
189 auto& bfld = result.boundaryFieldRef();
207 sigma_.boundaryField()[patchi].patchNeighbourField()
215 const label own = pOwner[facei];
217 const scalar sigmaDeltaO =
220 const scalar sigmaDeltaN =
224 (vfO[facei]*sigmaDeltaO + vfN[facei]*sigmaDeltaN)
225 /(sigmaDeltaO + sigmaDeltaN);