-
Notifications
You must be signed in to change notification settings - Fork 911
Fix enthalpy related issues #2652
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
base: develop
Are you sure you want to change the base?
Conversation
Co-authored-by: Cristopher Morales <98025159+Cristopher-Morales@users.noreply.github.com>
…into fix_restart_scalars
…into fix_restart_scalars
| const su2double t_ref = config->GetStandard_RefTemperatureND(); | ||
| Enthalpy_Ref = 0.0; | ||
| su2double t_i = 1.0; | ||
| for (int i = 0; i < N; ++i) { | ||
| t_i *= t_ref; | ||
| t_i *= STD_REF_TEMP; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
How does non-dimensionalization work here?
| su2double t_i = 1.0; | ||
| for (int i = 0; i < N; ++i) { | ||
| t_i *= STD_REF_TEMP; | ||
| t_i *= STD_REF_TEMP / config->GetInc_Temperature_Ref(); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Having a std_ref_temp_nd in CConfig probably makes sense if this will be used a lot
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
we already have a temperature_ref and a temperature_init, I think we should re-use that.
| Note Cp = Cv, (gamma = 1).*/ | ||
| Temperature = t; | ||
| Enthalpy = Cp * Temperature; | ||
| Enthalpy = Cp * (Temperature - STD_REF_TEMP); // Sensible enthalpy relative to STD_REF_TEMP |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
What about here? Is it always dimensional?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
It looks like it. I'll introduce Temperature_Ref for enthalpy nondimensionalization upstream of this..
…into fix_restart_scalars
Proposed Changes
Related Work
small fixes for #2426
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.
pre-commit run --allto format old commits.