Proposed Changes

This fixes an issue introduced in #2525 where the under-relaxation factors for $k$ and $\omega$ were calculated incorrectly with the SST model.

Since the update for $k$ and $\omega$ is in terms of conserved transport quantities, their stored value must be multiplied by the local flow density when compared to the non-linear update. This was not the case previously, and the solver could overestimate the relative change and compute a much stricter relaxation factor.

This is particularly relevant to incompressible flow, such as the steam periodic flow tutorial case, where the density is several orders of magnitude. The consequence of this issue was a stalled convergence of the turbulence variables as discussed in #2594.

This issue did not affect problems using the SA model, which doesn't use a conserved quantity in its formulation, or where the normalized density is near unity.

Related Work

Related to #2525 and #2594

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.