derivative bug fix in divu calculation (#263)

divu in spherical is 1/r^2 d/dr(r^2 * u).
python-hydro · Sep 9, 2024 · d11ac15 · d11ac15
1 parent 73d8038
commit d11ac15
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/pyro/compressible/interface.py b/pyro/compressible/interface.py
@@ -329,8 +329,8 @@ def artificial_viscosity(ng, dx, dy, Lx, Ly,
                 vb = 0.5 * (v[i, j - 1] + v[i - 1, j - 1])
 
                 # Finite difference to get ux and vy
-                ux = (ur*rr - ul*rl) / (rc * dx)
-                vy = (sint*vt - sinb*vb) / (rc * sinc * dy)
+                ux = (ur * rr * rr - ul * rl * rl) / (rc * rc * dx)
+                vy = (sint * vt - sinb * vb) / (rc * sinc * dy)
 
                 # Find div(U)_{i-1/2, j-1/2}
                 divU[i, j] = ux + vy