dynamic_smagorinsky_coeffs

[bronise777]

In the pc_dynamic_smagorinsky_coeffs function, should M here be the difference of the stress tensor instead of the product?
// "resolved turbulent stresses" and others
for (int dir = 0; dir < AMREX_SPACEDIM * AMREX_SPACEDIM; dir++) {
M[dir] = betaij[dir] * alphaij(i, j, k, dir);
}

[marchdf]

Can you point to the line in a paper that is making you think this? You could be right.

[bronise777]

Formula 17 on page 365 in Martín, M. Pino, U. Piomelli, and G. V. Candler. "Subgrid-Scale Models for Compressible Large-Eddy
Simulations." Theoretical and Computational Fluid Dynamics 13, no. 5 (2000): 361–76. https://link.springer.com/article/10.1007/PL00020896

[marchdf]

Yeah it definitely looks like a bug... I opened a PR to fix. Thank you!