Fix bug in LES coefficient calculation. (#846)

AMReX-Combustion · Sep 20, 2024 · 9c33a6f · 9c33a6f
1 parent 400b84f
commit 9c33a6f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/Source/LES.H b/Source/LES.H
@@ -417,7 +417,7 @@ pc_dynamic_smagorinsky_coeffs(
 
   // "resolved turbulent stresses" and others
   for (int dir = 0; dir < AMREX_SPACEDIM * AMREX_SPACEDIM; dir++) {
-    M[dir] = betaij[dir] * alphaij(i, j, k, dir);
+    M[dir] = betaij[dir] - alphaij(i, j, k, dir);
   }
 
   L[i00] = Kij(i, j, k, 0) - q(i, j, k, QRHO) * q(i, j, k, QU) * q(i, j, k, QU);