[david]

Hi all

I encountered a problem with the forces functionObject. I have a simulation with different patches in the domain. Some of the patches have similar names, e.g.:

BLADE
HUB-BLADE
SHROUD-BLADE


I wanted to calculate the forces acting on the patch "BLADE". Unfortunatelly the forces functionObject adds the forces of all patches with the string "BLADE" in its name.

Additionally I have a question concerning the moments. As far as I understand the code, the moments are calculated with the cell centers:


vectorField Md = mesh.C().boundaryField()[patchi] - CofR_;

vectorField pf =
mesh.Sf().boundaryField()[patchi]*p.boundaryField()[patchi];

fm.first().first() += rho(p)*sum(pf);
fm.second().first() += rho(p)*sum(Md ^ pf);

Wouldn't it be better to use the face centers (mesh.Cf())?


Thanks a lot & best regards
David