From 490c0e6efbe61ae97bb67b29d5bc5db32f56f5bb Mon Sep 17 00:00:00 2001
From: Thomas Lichtenegger <thomas.lichtenegger@jku.at>
Date: Thu, 4 Aug 2016 14:03:22 +0200
Subject: [PATCH] Add smoothing of dyn hold up.

---
 .../cfdemParticle/subModels/forceModel/Fines/FinesFields.C    | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/src/lagrangian/cfdemParticle/subModels/forceModel/Fines/FinesFields.C b/src/lagrangian/cfdemParticle/subModels/forceModel/Fines/FinesFields.C
index fb97fa50..4c483cab 100644
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/Fines/FinesFields.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/Fines/FinesFields.C
@@ -20,7 +20,7 @@ License
 #include "FinesFields.H"
 #include "addToRunTimeSelectionTable.H"
 #include "averagingModel.H"
-
+#include "smoothingModel.H"
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 namespace Foam
@@ -374,6 +374,8 @@ void FinesFields::integrateFields()
     alphaStEqn.solve();
     alphaDynEqn.solve();
     
+    particleCloud_.smoothingM().smoothen(alphaDyn_);
+    
     // limit hold-ups, should be done more elegantly
     
     scalar alphaStErr(0.0);