BUG: Curle - added missing division by r. See #2169

This commit is contained in:
Andrew Heather
2021-08-02 18:50:48 +01:00
parent 7f2cbd6cd4
commit 9645e1adc2

View File

@ -213,11 +213,11 @@ bool Foam::functionObjects::Curle::execute()
forAll(observerPositions_, pointi) forAll(observerPositions_, pointi)
{ {
const vectorField r(Cfp - observerPositions_[pointi]); const vectorField r(observerPositions_[pointi] - Cfp);
const scalarField invMagR(1/(mag(r) + ROOTVSMALL)); const scalarField invMagR(1/(mag(r) + ROOTVSMALL));
pDash[pointi] += pDash[pointi] +=
sum((pp*sqr(invMagR) + invMagR/c0_*dpdtp)*(Sfp & r)); sum((pp*sqr(invMagR) + invMagR/c0_*dpdtp)*(Sfp & (r*invMagR)));
} }
} }