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: