bubble-dynamics

commit 568c00aa72a9cef6a69769409285fb88fcf922d4
parent faf316b78264cb00f2bb24ce055851267c0e8191
Author: Fabian Wermelinger <fabianw@mavt.ethz.ch>
Date:   Wed, 25 Jan 2017 17:07:18 +0100

Bug fix in rhs term involving derivative of external presure

Credit to U. Rasthofer

Diffstat:
M
src/kernels/KMCluster_TY.h
 | 
2
+-

1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/kernels/KMCluster_TY.h b/src/kernels/KMCluster_TY.h
@@ -121,7 +121,7 @@ public:
             // b(i) -= bd.rInv*(static_cast<Real>(1.0) + Rdc)*(pstat + bd.pext(t));
             b(i) += bd.rInv*(pbs - static_cast<Real>(2.0)*bd.sigma/Ui[0]) - static_cast<Real>(4.0)*bd.nuL*Ui[1]/Ui[0];
             b(i) += bd.rInv*static_cast<Real>(1.0 - 3.0*bd.gamma)*Rdc*pbs;
-            b(i) -= bd.rInv*Rdc*bd.dpext(t + Ui[0]*clInv);
+            b(i) -= bd.rInv*clInv*Ui[0]*bd.dpext(t + Ui[0]*clInv);
             // b(i) -= bd.rInv*Rdc*bd.dpext(t);
             b(i) -= static_cast<Real>(2.0)*bnbr;