From 1cd4b292eb94446d0aae9d34b096e8a1da93630a Mon Sep 17 00:00:00 2001
From: Henry Weller <http://openfoam.org>
Date: Thu, 4 Jul 2019 14:25:21 +0100
Subject: [PATCH] fvOptions::damping: Corrected handling of multiple damping
 regions

Resolves bug-report https://bugs.openfoam.org/view.php?id=3302
---
 src/fvOptions/sources/derived/damping/damping/damping.C | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/fvOptions/sources/derived/damping/damping/damping.C b/src/fvOptions/sources/derived/damping/damping/damping.C
index 5c7c395140..395aebc127 100644
--- a/src/fvOptions/sources/derived/damping/damping/damping.C
+++ b/src/fvOptions/sources/derived/damping/damping/damping.C
@@ -58,15 +58,15 @@ Foam::tmp<Foam::volScalarField::Internal> Foam::fv::damping::forceCoeff() const
     );
     scalarField& forceCoeff = tforceCoeff.ref();
 
-    const scalar lambda = lambda_.value();
-
     forAll(origins_, i)
     {
         const vectorField& c = mesh_.cellCentres();
         const scalarField x((c - origins_[i]) & directions_[i]);
-        forceCoeff = lambda*max(forceCoeff, scale_->value(x));
+        forceCoeff = max(forceCoeff, scale_->value(x));
     }
 
+    forceCoeff *= lambda_.value();
+
     // Write out the force coefficient for debugging
     if (debug && mesh_.time().writeTime())
     {