zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

COMP: avoid ambiguous construct from tmp - postProcessing

Browse Source
This commit is contained in:
Mark Olesen
2010-12-17 18:16:01 +01:00
parent 3d3c95a37a
commit d87a880ebc
3 changed files with 18 additions and 14 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
src/postProcessing/functionObjects/forces/forces/forces.C
View File

@ -408,22 +408,24 @@ Foam::forces::forcesMoments Foam::forces::calcForcesMoment() const
{
label patchi = iter.key();
vectorField Md = mesh.C().boundaryField()[patchi] - CofR_;
vectorField Md(mesh.C().boundaryField()[patchi] - CofR_);
scalarField sA = mag(Sfb[patchi]);
scalarField sA(mag(Sfb[patchi]));
// Normal force = surfaceUnitNormal * (surfaceNormal & forceDensity)
vectorField fN =
vectorField fN
(
Sfb[patchi]/sA
*(
Sfb[patchi] & fD.boundaryField()[patchi]
);
)
);
fm.first().first() += sum(fN);
fm.second().first() += sum(Md ^ fN);
// Tangential force (total force minus normal fN)
vectorField fT = sA*fD.boundaryField()[patchi] - fN;
vectorField fT(sA*fD.boundaryField()[patchi] - fN);
fm.first().second() += sum(fT);
fm.second().second() += sum(Md ^ fT);
@ -450,14 +452,14 @@ Foam::forces::forcesMoments Foam::forces::calcForcesMoment() const
{
label patchi = iter.key();
vectorField Md = mesh.C().boundaryField()[patchi] - CofR_;
vectorField Md(mesh.C().boundaryField()[patchi] - CofR_);
vectorField pf = Sfb[patchi]*(p.boundaryField()[patchi] - pRef);
vectorField pf(Sfb[patchi]*(p.boundaryField()[patchi] - pRef));
fm.first().first() += rho(p)*sum(pf);
fm.second().first() += rho(p)*sum(Md ^ pf);
vectorField vf = Sfb[patchi] & devRhoReffb[patchi];
vectorField vf(Sfb[patchi] & devRhoReffb[patchi]);
fm.first().second() += sum(vf);
fm.second().second() += sum(Md ^ vf);