Commit e3baecec authored by Will Bainbridge's avatar Will Bainbridge Committed by Andrew Heather
Browse files

PairCollision: Removed unnecessary operations

Resolves bug report https://bugs.openfoam.org/view.php?id=2680
parent cf0d7a60
......@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
......@@ -229,8 +229,7 @@ void Foam::PairSpringSliderDashpot<CloudType>::evaluatePair
vector USlip_AB =
U_AB - (U_AB & rHat_AB)*rHat_AB
+ (pA.omega() ^ (dAEff/2*-rHat_AB))
- (pB.omega() ^ (dBEff/2*rHat_AB));
- ((dAEff/2*pA.omega() + dBEff/2*pB.omega()) ^ rHat_AB);
scalar deltaT = this->owner().mesh().time().deltaTValue();
......@@ -276,10 +275,7 @@ void Foam::PairSpringSliderDashpot<CloudType>::evaluatePair
}
else
{
fT_AB =
-kT*tangentialOverlapMag
*tangentialOverlap_AB/tangentialOverlapMag
- etaT*USlip_AB;
fT_AB = - kT*tangentialOverlap_AB - etaT*USlip_AB;
}
pA.f() += fT_AB;
......
......@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
......@@ -142,10 +142,7 @@ void Foam::WallSpringSliderDashpot<CloudType>::evaluateWall
}
else
{
fT_PW =
-kT*tangentialOverlapMag
*tangentialOverlap_PW/tangentialOverlapMag
- etaT*USlip_PW;
fT_PW = - kT*tangentialOverlap_PW - etaT*USlip_PW;
}
p.f() += fT_PW;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment