56 template<
class Type,
class NVDweight>
99 if (edgeFlux_[
edge] > 0)
103 d /=
mag(d)/
mesh.edgeInterpolation::lPN().internalField()[
edge];
108 d -=
n[neighbour[
edge]]*(
n[neighbour[
edge]]&d);
109 d /=
mag(d)/
mesh.edgeInterpolation::lPN().internalField()[
edge];
120 gradc[neighbour[
edge]],
131 if (bWeights[patchI].coupled())
135 const scalarField& pEdgeFlux = edgeFlux_.boundaryField()[patchI];
153 mesh.areaCentres().boundaryField()[patchI].patchInternalField()
158 mesh.areaCentres().boundaryField()[patchI]
159 .patchNeighbourField()
164 mesh.faceAreaNormals().boundaryField()[patchI]
165 .patchInternalField()
170 mesh.faceAreaNormals().boundaryField()[patchI]
171 .patchNeighbourField()
176 mesh.edgeInterpolation::lPN().boundaryField()[patchI]
183 if (pEdgeFlux[edgeI] > 0)
185 d = CN[edgeI] - CP[edgeI];
186 d -= nP[edgeI]*(nP[edgeI]&d);
187 d /=
mag(d)/pLPN[edgeI];
191 d = CN[edgeI] - CP[edgeI];
192 d -= nN[edgeI]*(nN[edgeI]&d);
193 d /=
mag(d)/pLPN[edgeI];
211 return tWeightingFactors;