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

Cleanup wall functions and add user controls #1360

Merged
merged 125 commits into from
Oct 20, 2021

Conversation

bigfooted
Copy link
Contributor

Proposed Changes

Give a brief overview of your contribution here in a few sentences.

  • fixed a small bug in the wall model
  • added config options (replace as many hard coded constants as possible)
  • aerodynamic heating (wall temperature imposed in compressible flow, still needs a high mach testcase)
  • added a test case for SA model
  • cleaned up a bit

Related Work

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

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 the '-Wall -Wextra -Wno-unused-parameter -Wno-empty-body' compiler flags, or simply --warnlevel=2 when using meson).
  • My contribution is commented and consistent with SU2 style.
  • 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 and others added 30 commits February 21, 2021 09:18
Co-authored-by: TobiKattmann <31306376+TobiKattmann@users.noreply.github.com>
Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
spaces

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
spaces

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
refvel

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
const bool

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
spaces

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
spaces

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
spaces

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
away with the comment blocks!

Co-authored-by: Pedro Gomes <38071223+pcarruscag@users.noreply.github.com>
@TobiKattmann
Copy link
Contributor

Changes that are necessary in a follow-up PR:

  • Move regression tests to parallel regression and hybrid reg. Maybe the 4 cases can be split over those.
  • Find similarities between IncNS and NS TauWallWF implementation and make it a common implementation.

Please add to this post other things (I am fine with others editing this)

@TobiKattmann
Copy link
Contributor

i made a little Test for the Non-dimensionalization: Turns out the inc Testcases did not specify INC_NONDIM where the default is INITIAL_VALUES. I will specify that for the reg tests.

I think the non-dim is currently correct: attached the incomp-SA case comparison INITIAL_VALUES against DIMENSIONAL.

image

Here the non-dim value is dimensionalized using what gets printed in the screen output
image

Comment on lines 2557 to 2560

if (!wallfunctions) {
YPlus[iMarker][iVertex] = WallDistMod * FrictionVel / (Viscosity / Density);
}
Copy link
Contributor

Choose a reason for hiding this comment

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

This change made the comp SST WF reg test fail. There are no changes for the inc SST case and the SA cases, though.

Comment on lines -499 to +502
if (Y_Plus < 5.0) continue;
if (Y_Plus < minYPlus) {
/* --- Use zero flux (Neumann) conditions, i.e. nothing has to be done. --- */
continue;
}
Copy link
Member

Choose a reason for hiding this comment

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

That might be a positive change @TobiKattmann, the recomputation of y+ could be interfering with the logic here, affecting the boundary conditions of the turb solver.

@@ -321,7 +321,7 @@ def main():
turb_wallfunction_flatplate_sst.cfg_dir = "wallfunctions/flatplate/compressible_SST"
turb_wallfunction_flatplate_sst.cfg_file = "turb_SST_flatplate.cfg"
turb_wallfunction_flatplate_sst.test_iter = 10
turb_wallfunction_flatplate_sst.test_vals = [-4.230000, -1.930543, -1.998684, 1.250334, -1.635534, 1.462491, 10, -2.151894, 0.072872, 0.002514] #last 10 columns
turb_wallfunction_flatplate_sst.test_vals = [-4.229955, -1.930560, -1.998477, 1.250383, -1.635663, 1.462396, 10.000000, -2.151959, 0.072873, 0.002514] #last 10 columns
Copy link
Contributor

Choose a reason for hiding this comment

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

Thanks for this "final" addition and all the help on this (where of course @bigfooted did the heavy lifting ;) )

Now this is good to go from my perspective, I have no further comments. Maybe we can merge this in todays dev meeting :)

Things tbd in the future are

  • consolidate compr and inc implementation into some unified entitiy
  • Distribute Testcases over parallel, hybrid regression test
  • explanation on the website, tutorial

Copy link
Contributor

Choose a reason for hiding this comment

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

Alright merging this in now ✈️

@TobiKattmann TobiKattmann merged commit 93f78c4 into su2code:develop Oct 20, 2021
@TobiKattmann
Copy link
Contributor

TobiKattmann commented Nov 10, 2021

See #1430 for some more things to fix

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