43 patchMeanVelocityForce,
52 Foam::fv::patchMeanVelocityForce::patchMeanVelocityForce
54 const word& sourceName,
55 const word& modelType,
61 patch_(coeffs_.get<
word>(
"patch")),
62 patchi_(
mesh.boundaryMesh().findPatchID(patch_))
67 <<
"Cannot find patch " << patch_
84 (flowDir_ &
U.boundaryField()[patchi_])
85 *mesh_.boundary()[patchi_].magSf()
87 sum(mesh_.boundary()[patchi_].magSf())
96 if (Pstream::parRun() && isA<cyclicPolyPatch>(
patches[patchi_]))
100 processorCyclicPolyPatch::patchIDs(patch_,
patches)
103 forAll(processorCyclicPatches, pcpi)
105 const label patchi = processorCyclicPatches[pcpi];
110 (flowDir_ &
U.boundaryField()[patchi])
111 *mesh_.boundary()[patchi].magSf()
114 sumAmagUsumA.
y() +=
sum(mesh_.boundary()[patchi].magSf());
120 return sumAmagUsumA.
x()/sumAmagUsumA.
y();