fix bug in adhesion (#308)

InteractiveComputerGraphics · Mar 6, 2024 · dc2f63f · dc2f63f
1 parent 9696086
commit dc2f63f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/SPlisHSPlasH/SurfaceTension/SurfaceTension_Jeske2023.cpp b/SPlisHSPlasH/SurfaceTension/SurfaceTension_Jeske2023.cpp
@@ -582,7 +582,7 @@ void SurfaceTension_Jeske2023::computeRHS(VectorXr &b, VectorXr &g) {
                                                             // + getDensityGrad(neighborIndex)
                                                             ) / (rhoij * rhoij)
                         );
-                        const Vector3r a = adhesion * xixj * density_i * W_firstOrder;
+                        const Vector3r a = adhesion * (xixj - dt*vj)  * density_i * W_firstOrder;
                         bi += a;
                         bm_neighbor->addForce(xi, -m_model->getMass(i) / density_i * a);
                 );
@@ -1159,7 +1159,7 @@ void SurfaceTension_Jeske2023::computeRHS(VectorXr &b, VectorXr &g) {
                                                             // + getDensityGrad(neighborIndex)
                                                             ) / (rhoij * rhoij)
                         );
-                        bi += adhesion * xixj * density_i * W_firstOrder;
+                        bi += adhesion * (xixj - dt * vj) * density_i * W_firstOrder;
 
                 );
             }