From cccd62e45007290a62ebf951d9ae1dbb113e2a22 Mon Sep 17 00:00:00 2001 From: BehradEsg Date: Mon, 16 Oct 2023 15:40:21 +0200 Subject: [PATCH] getting the particle density from particle cloud --- .../subModels/forceModel/staticPressure/staticPressure.C | 3 +-- .../subModels/forceModel/staticPressure/staticPressure.H | 2 -- 2 files changed, 1 insertion(+), 4 deletions(-) diff --git a/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.C b/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.C index 02e8978c..c7b92f04 100644 --- a/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.C +++ b/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.C @@ -56,7 +56,6 @@ staticPressure::staticPressure g_(propsDict_.lookup ("g_DEM")), voidfractionFieldName_(propsDict_.lookup("voidfractionFieldName")), voidfraction_(sm.mesh().lookupObject (voidfractionFieldName_)), - rhoPart_(readScalar(propsDict_.lookup ("rhoPart"))), solidFraction_(readScalar(propsDict_.lookup ("DomainSolidVolumeFraction"))) { // init force sub model @@ -76,7 +75,7 @@ void staticPressure::setForce() const { #include "setupProbeModel.H" label cellI; - scalar rhoMix_ = solidFraction_*rhoPart_ + (1.0-solidFraction_)*rhoGas_; + scalar rhoMix_ = solidFraction_*particleCloud_.particleDensity(0) + (1.0-solidFraction_)*rhoGas_; vector force; for(int index = 0; index < particleCloud_.numberOfParticles(); ++index) diff --git a/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.H b/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.H index 4f1fba88..868619b8 100644 --- a/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.H +++ b/src/lagrangian/cfdemParticle/subModels/forceModel/staticPressure/staticPressure.H @@ -65,8 +65,6 @@ private: const volScalarField& voidfraction_; - scalar rhoPart_; - scalar solidFraction_; public: