diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C index 758e58743b1291dec9828bd73a8e33f0f02a656b..eed0c800d08d30ade1ba1f33d9f48f37f267adc2 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -95,7 +95,9 @@ Foam::forceSuSp Foam::SRFForce<CloudType>::calcNonCoupled const vector& r = p.position(); // Coriolis and centrifugal acceleration terms - value.Su() = mass*(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega))); + value.Su() = + mass*(1.0 - p.rhoc()/p.rho()) + *(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega))); return value; }