From 2f5da5af4b97ccfb8ed9643f802eae7b882c4302 Mon Sep 17 00:00:00 2001
From: Andrew Heather <>
Date: Mon, 2 Aug 2021 18:50:48 +0100
Subject: [PATCH] BUG: Curle - added missing division by r. See #2169
---
src/functionObjects/field/Curle/Curle.C | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/functionObjects/field/Curle/Curle.C b/src/functionObjects/field/Curle/Curle.C
index 39f38def4ec..1a8133a83de 100644
--- a/src/functionObjects/field/Curle/Curle.C
+++ b/src/functionObjects/field/Curle/Curle.C
@@ -213,11 +213,11 @@ bool Foam::functionObjects::Curle::execute()
forAll(observerPositions_, pointi)
{
- const vectorField r(Cfp - observerPositions_[pointi]);
+ const vectorField r(observerPositions_[pointi] - Cfp);
const scalarField invMagR(1/(mag(r) + ROOTVSMALL));
pDash[pointi] +=
- sum((pp*sqr(invMagR) + invMagR/c0_*dpdtp)*(Sfp & r));
+ sum((pp*sqr(invMagR) + invMagR/c0_*dpdtp)*(Sfp & (r*invMagR)));
}
}
--
GitLab