Scaled mass matrix needs to use kinematic viscosity (momentum diffusi…

…vity)
jedbrown · Jun 6, 2011 · cf03b85 · cf03b85
1 parent a8a9c14
commit cf03b85
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/fs/tests/vht.c b/src/fs/tests/vht.c
@@ -1844,7 +1844,7 @@ static dErr VHTJacobianAssemble_PressureEnergy(VHT vht,Mat Bpp,Mat Daux,Mat Bee,
         for (dInt j=0; j<Pp; j++) { /* trial functions */
           for (dInt i=0; i<Pp; i++) {
             /* Scaled mass matrx */
-            Kpp[i][j] += interpp[q][i] * jw[q] * (1./stash.eta.x) * interpp[q][j];
+            Kpp[i][j] += interpp[q][i] * jw[q] * (stash.rho.x/stash.eta.x) * interpp[q][j];
             /* Neumann Laplacian */
             Kppaux[i][j] += (+ derivp[q][i][0] * jw[q] * derivp[q][j][0]
                              + derivp[q][i][1] * jw[q] * derivp[q][j][1]