Fix sign on ueint pdivu

AMReX-Fluids · Jan 3, 2024 · 366c24c · 366c24c
1 parent c816b03
commit 366c24c
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/Source/Hydro/Hydro_utils_eb_K.H b/Source/Hydro/Hydro_utils_eb_K.H
@@ -104,7 +104,7 @@ eb_pdivu (int i, int j, int k,
                        (q1(i + 1, j, k, GDU   ) - q1(i, j, k, GDU   )) *dxinv[0]
                       +(q2(i, j + 1, k, GDPRES) + q2(i, j, k, GDPRES)) *
                        (q2(i, j + 1, k, GDV   ) - q2(i, j, k, GDV   )) * dxinv[1]);
-        divu(i,j,k,UEINT) += pdivu;
+        divu(i,j,k,UEINT) -= pdivu;
     }
     else if (!flag(i,j,k).isCovered())
     {
@@ -153,7 +153,7 @@ eb_pdivu (int i, int j, int k,
             fyp = amrex::Real(0.0);
         }
 
-        divu(i,j,k,UEINT) += (1.0/vfrc(i,j,k)) * q(i,j,k,QPRES)*
+        divu(i,j,k,UEINT) -= (1.0/vfrc(i,j,k)) * q(i,j,k,QPRES)*
             ( dxinv[0] * (apx(i+1,j,k)*fxp-apx(i,j,k)*fxm)
             + dxinv[1] * (apy(i,j+1,k)*fyp-apy(i,j,k)*fym) );
     }
@@ -169,7 +169,7 @@ eb_pdivu (int i, int j, int k,
                      +(q3(i,j,k+1, GDPRES) + q3(i,j,k, GDPRES)) *
                       (q3(i,j,k+1, GDW   ) - q3(i,j,k, GDW   )) * dxinv[2]);
 
-        divu(i,j,k,UEINT) += pdivu;
+        divu(i,j,k,UEINT) -= pdivu;
     }
     else if (!flag(i,j,k).isCovered())
     {
@@ -271,7 +271,7 @@ eb_pdivu (int i, int j, int k,
                 fzp = amrex::Real(0.0);
         }
 
-        divu(i,j,k,UEINT) += (1.0/vfrc(i,j,k)) * q(i,j,k,QPRES)*
+        divu(i,j,k,UEINT) -= (1.0/vfrc(i,j,k)) * q(i,j,k,QPRES)*
             ( dxinv[0] * (apx(i+1,j,k)*fxp-apx(i,j,k)*fxm)
             + dxinv[1] * (apy(i,j+1,k)*fyp-apy(i,j,k)*fym)
             + dxinv[2] * (apz(i,j,k+1)*fzp-apz(i,j,k)*fzm) );