Skip to content

feat:Implement Dynamic Smagorinsky SGS model for LES (Germano-Lilly, local approximation)#2850

Open
Vamshimadhava wants to merge 4 commits into
su2code:masterfrom
Vamshimadhava:implement-DynamicSmagorinsky
Open

feat:Implement Dynamic Smagorinsky SGS model for LES (Germano-Lilly, local approximation)#2850
Vamshimadhava wants to merge 4 commits into
su2code:masterfrom
Vamshimadhava:implement-DynamicSmagorinsky

Conversation

@Vamshimadhava

@Vamshimadhava Vamshimadhava commented Jul 22, 2026

Copy link
Copy Markdown

Proposed Changes

Implements the CDynamicSmagorinskyModel::ComputeEddyViscosity_3D method in sgs_model.inl.

The implementation follows the Germano–Lilly dynamic procedure using a single-point approximation of the test filter, making it compatible with SU2's unstructured DG FEM solver architecture, which does not expose neighbor stencils at the SGS model level.

The Lilly least-squares estimator is used to compute the local Smagorinsky coefficient (Cs²), which is clipped to zero to prevent backscatter.

Related Work
Germano, M., Piomelli, U., Moin, P., & Cabot, W. H. (1991). A dynamic subgrid-scale eddy viscosity model. Physics of Fluids A, 3(7), 1760–1765.
Lilly, D. K. (1992). A proposed modification of the Germano subgrid-scale closure method. Physics of Fluids A, 4(3), 633–635.
Related to the existing ones CSmagorinskyModel and CWALEModel implementations in the same file.

PR Checklist

Put an X by all that apply. You can fill this out after submitting the PR. If you have any questions, don't hesitate to ask! We want to help. These are a guide for you to know what the reviewers will be looking for in your contribution.

  • I am submitting my contribution to the develop branch.
  • My contribution generates no new compiler warnings (try with --warnlevel=3 when using meson).
  • My contribution is commented and consistent with SU2 style (https://su2code.github.io/docs_v7/Style-Guide/).
  • I used the pre-commit hook to prevent dirty commits and used pre-commit run --all to format old commits.
  • I have added a test case that demonstrates my contribution, if necessary.
  • I have updated appropriate documentation (Tutorials, Docs Page, config_template.cpp), if necessary.

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