From aa82dad8be8358e3985fd5e491d41c61aac9297c Mon Sep 17 00:00:00 2001
From: asanaz <31247086+asanaz@users.noreply.github.com>
Date: Wed, 16 Sep 2020 10:45:51 +0200
Subject: [PATCH] Update SchillerNaumannDrag.C

The changes are added.
---
 .../forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C      | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C b/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C
index 398d1590..c274e462 100644
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/SchillerNaumannDrag/SchillerNaumannDrag.C
@@ -87,7 +87,7 @@ SchillerNaumannDrag::SchillerNaumannDrag
     // read those switches defined above, if provided in dict
     forceSubM(0).readSwitches();
 
-    particleCloud_.checkCG(false);
+    particleCloud_.checkCG(true);
 
     if (propsDict_.found("scale"))
         scaleDia_ = scalar(readScalar(propsDict_.lookup("scale")));
@@ -151,7 +151,7 @@ void SchillerNaumannDrag::setForce() const
                     Cd = max(0.44,24.0/Rep*(1.0+0.15*pow(Rep,0.687)));
 
                     // calc particle's drag
-                    drag = 0.125*Cd*rho*M_PI*ds*ds_scaled*magUr*Ur;
+                    drag = 0.125*Cd*rho*M_PI*ds*ds*scaleDia_*magUr*Ur*scaleDrag_;
 
                     if (modelType_=="B")
                         drag /= voidfraction;