Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix bug in LES coefficient calculation. #846

Merged
merged 1 commit into from
Sep 20, 2024
Merged

Fix bug in LES coefficient calculation. #846

merged 1 commit into from
Sep 20, 2024

Conversation

marchdf
Copy link
Contributor

@marchdf marchdf commented Sep 20, 2024

See 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

Thanks to @bronise777 for finding this. Discussed in #845

See 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

Thanks to @bronise777 for finding this.
@marchdf marchdf requested a review from baperry2 September 20, 2024 14:54
Copy link
Contributor

@baperry2 baperry2 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Digging up the old fortran version, it looks like we used to have - instead of *, so this was probably a transcription error when porting to C++

M(:) = betaij(:) - alphaij(i,j,k,:)

@marchdf
Copy link
Contributor Author

marchdf commented Sep 20, 2024

Wow nice sleuthing!! I wonder how this happened...

@marchdf marchdf merged commit 9c33a6f into development Sep 20, 2024
13 checks passed
@marchdf marchdf deleted the fix-les-bug branch September 20, 2024 18:46
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants