From b95b6a0afda796f3d9210c784a3202bc41fb90a3 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Tue, 21 Feb 2023 16:12:16 +0000
Subject: [PATCH] lagrangian: PatchCollisionDensity: Corrected number-density
 fields

---
 .../PatchCollisionDensity/PatchCollisionDensity.C            | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/lagrangian/parcel/submodels/CloudFunctionObjects/PatchCollisionDensity/PatchCollisionDensity.C b/src/lagrangian/parcel/submodels/CloudFunctionObjects/PatchCollisionDensity/PatchCollisionDensity.C
index 5f9d10dbfa..c64aa32829 100644
--- a/src/lagrangian/parcel/submodels/CloudFunctionObjects/PatchCollisionDensity/PatchCollisionDensity.C
+++ b/src/lagrangian/parcel/submodels/CloudFunctionObjects/PatchCollisionDensity/PatchCollisionDensity.C
@@ -301,9 +301,10 @@ void Foam::PatchCollisionDensity<CloudType>::postPatch
         const scalar magSf =
             this->owner().mesh().magSf().boundaryField()[patchi][patchFacei];
 
-        numberCollisionDensity_[patchi][patchFacei] += 1/magSf;
+        numberCollisionDensity_[patchi][patchFacei] +=
+            p.nParticle()/magSf;
         massCollisionDensity_[patchi][patchFacei] +=
-            p.mass()*p.nParticle()/magSf;
+            p.nParticle()*p.mass()/magSf;
     }
 }