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 step-57 bilinear form #15442

Merged
merged 1 commit into from Jun 23, 2023
Merged

Conversation

tjhei
Copy link
Member

@tjhei tjhei commented Jun 22, 2023

This reverts part of #15127

fixes #15433

@tjhei tjhei added this to the Release 9.5 milestone Jun 22, 2023
Comment on lines -460 to +461
phi_u[i] * (present_velocity_values[q] *
present_velocity_gradients[q]) +
phi_u[i] * (present_velocity_gradients[q] *
present_velocity_values[q]) +
Copy link
Member

Choose a reason for hiding this comment

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

Does this mean that we accidentally multiplied by (grad u)^T instead of (grad u)?

@bangerth
Copy link
Member

Just to be extra sure: You tried this and it works again?

@tjhei
Copy link
Member Author

tjhei commented Jun 22, 2023

Does this mean that we accidentally multiplied by (grad u)^T instead of (grad u)?

Yes. When the terms got rearranged in #15127, we did not pay attention to the fact that grad_u is a 2-tensor.

Just to be extra sure: You tried this and it works again?

Yes. Without this change the program does not converge correctly (debug mode and release mode) and this is the minimal part of the last PR to revert to make it work again.

@kronbichler kronbichler merged commit c6b3d46 into dealii:master Jun 23, 2023
13 of 14 checks passed
@tjhei tjhei deleted the fix-step-57-convergence branch June 23, 2023 15:48
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.

step-57 stalls
4 participants