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

SA 2D axisymmetric source terms #2197

Merged
merged 23 commits into from
Aug 11, 2024
Merged

Conversation

Corentinthomee
Copy link
Contributor

Proposed Changes

Give a brief overview of your contribution here in a few sentences.
Additional source terms in the SA turbulence model have been added for the axisymmetric formulation.

In the baseline model, the only term that has an impact is:
image
Rewriting it in cylindrical coordinates, an additional source term is obtained:
image

In the Negative variant of the model, this source varies slightly due to the presence of $f_n$:
image

The compressibility correction also requires a new source. Rewriting
image
in cylindrical coordinates yields the new source
image

For the other model variants, it looks like the modification of strain magnitude (done for SST) takes care of the differences, but I might be mistaken.

Related Work

Resolve any issues (bug fix or feature request), note any related PRs, or mention interactions with the work of others, if any.
#1565

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.

@bigfooted
Copy link
Contributor

Great, we needed this! Do you have a testcase that we can use as a regression test to add to the PR?

@Corentinthomee
Copy link
Contributor Author

I'm thinking about reusing the SST test case (axisymmetric_rans/air_nozzle) so that the same mesh can be used. Would that be fine ?

@bigfooted
Copy link
Contributor

I'm thinking about reusing the SST test case (axisymmetric_rans/air_nozzle) so that the same mesh can be used. Would that be fine ?

That's a quite specific testcase. How did you test the implementation? Do you have some simple results for the turbulent axisymmetric jet? We know that the results for round jets are not super good for standard SA, but at least we know what the result should look like, here on p.32 - 35 they test SA for the round jet and compare to the Wygnanski & Fiedler measurements:
https://www.researchgate.net/publication/24296213_Turbulence_Modeling_Validation_Testing_and_Development

The domain is just a rectangle so pretty simple to set up.

You could also use the V&V test that we have, but it is a variable density jet:
https://su2code.github.io/vandv/SANDIA_jet/

Copy link
Member

@pcarruscag pcarruscag left a comment

Choose a reason for hiding this comment

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

Nice work 👍 clean code

@Corentinthomee
Copy link
Contributor Author

I'm thinking about reusing the SST test case (axisymmetric_rans/air_nozzle) so that the same mesh can be used. Would that be fine ?

That's a quite specific testcase. How did you test the implementation? Do you have some simple results for the turbulent axisymmetric jet? We know that the results for round jets are not super good for standard SA, but at least we know what the result should look like, here on p.32 - 35 they test SA for the round jet and compare to the Wygnanski & Fiedler measurements: https://www.researchgate.net/publication/24296213_Turbulence_Modeling_Validation_Testing_and_Development

The domain is just a rectangle so pretty simple to set up.

You could also use the V&V test that we have, but it is a variable density jet: https://su2code.github.io/vandv/SANDIA_jet/

I did some tests on a supersonic jet, since I already had the files, but it's a quite complex case so it's not that good to verify the implementation. I will look into the simpler cases you mentioned !

@bigfooted
Copy link
Contributor

Any update on a testcase to add to the PR?

@Corentinthomee
Copy link
Contributor Author

Sorry, I have been quite busy the past few weeks.
I will try to get it done in the next few days.

@pcarruscag
Copy link
Member

@Corentinthomee can you add a testcase? otherwise please update AUTHORS.md so we can merge

@bigfooted
Copy link
Contributor

I will add a testcase in a separate PR.

@bigfooted bigfooted mentioned this pull request Aug 7, 2024
6 tasks
@pcarruscag pcarruscag merged commit b843f3f into su2code:develop Aug 11, 2024
35 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants