Skip to content

MSW"+" and vectorization#2727

Open
pcarruscag wants to merge 5 commits intodevelopfrom
pedro/msw+
Open

MSW"+" and vectorization#2727
pcarruscag wants to merge 5 commits intodevelopfrom
pedro/msw+

Conversation

@pcarruscag
Copy link
Member

Proposed Changes

Version of MSW that uses a pressure sensor for switching between SW and MSW instead of just the i and j pressures (which are not very reliable with MUSCL). Also adding the vectorized version.
The paper I had as a reference for MSW mentioned this sensor, and I decided to give it a try.

PR Checklist

  • 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

what's the paper?

Comment on lines +91 to +98
/*--- For completeness, the original formulation for the switch to unmodified Steger-Warming is:
* dp = fabs(P_j - P_i) / fmin(P_j, P_i)
* w = 0.5 * (1 / (pow(alpha * dp, 2) + 1))
* where alpha is a tuning constant to make the scheme switch sooner (large alpha).
* We're using a version of the switch based on the maximum "dp" over neighbors, from
* "Development of an Unstructured Navier-Stokes Solver For Hypersonic Nonequilibrium
* Aerothermodynamics". The 0.06 constant is to match our default alpha (5) to the threshold
* of 0.3 in the paper. ---*/
Copy link
Member Author

Choose a reason for hiding this comment

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

Copy link
Contributor

Choose a reason for hiding this comment

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

Also nice to see their solution for the axisymmetry issue.

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.

2 participants