LES sfs flux #857
bronise777
started this conversation in
General
Replies: 1 comment 1 reply
-
hmmmm I would agree on the missing 1/2. And it seems like there is a sign issue. Is it on the Qj term or the Jj term? I am swamped with some things right now but if you have an argument for one or the other, I would be happy to look into it in a couple of days. |
Beta Was this translation helpful? Give feedback.
1 reply
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
-
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
This is the energy equation
But in the pc_smagorinsky_sfs_term function,
flx(iv, UEDEN) = -sigmadx * Uface[0] - sigmady * Uface[1] - sigmadz * Uface[2];
flx(iv, UEDEN) = flx(iv, UEDEN) - cp * Cs2 / PrT * flux_T;
I have two problems, following the above formula sigmadx * Uface[0] sign seems to be positive and should have a coefficient of 1/2.
Beta Was this translation helpful? Give feedback.
All reactions