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

Fix implicit call for NEMO cases #1482

Merged
merged 1 commit into from
Jan 8, 2022
Merged

Fix implicit call for NEMO cases #1482

merged 1 commit into from
Jan 8, 2022

Conversation

WallyMaier
Copy link
Contributor

Proposed Changes

This changes a boolean check in CNEMONumerics.cpp which always seems to be true. This has been causing some segfaults in 10-11 species viscous cases. I suspect this is causes because the Kind_TimeNumScheme is not updated before the numerics constructor.

Some work that still needs to be done is confirm/test 10 and 11 species gases do work when the problem is supposed implicit.

Related Work

PR#1422

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.

@@ -53,7 +53,7 @@ CNEMONumerics::CNEMONumerics(unsigned short val_nDim, unsigned short val_nVar,
EDDY_VISC_INDEX = nSpecies+nDim+9;

/*--- Read from CConfig ---*/
implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
implicit = (config->GetKind_TimeIntScheme_Flow() == EULER_IMPLICIT);
Copy link
Contributor Author

Choose a reason for hiding this comment

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

This is always true for at least the first iteration.

Copy link
Member

Choose a reason for hiding this comment

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

Yep you need the physics specific versions (_Flow, _Turb, etc.) in constructors.

@WallyMaier
Copy link
Contributor Author

@CatarinaGarbacz please let me know if this fix resolves your issues.

@WallyMaier WallyMaier added the bug label Jan 7, 2022
@pr-triage pr-triage bot removed the PR: unreviewed label Jan 8, 2022
@WallyMaier WallyMaier merged commit 85b2b12 into develop Jan 8, 2022
@WallyMaier WallyMaier deleted the fix_implicit_call branch January 8, 2022 20:14
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.

3 participants