Skip to content

Consistent calculation of dissipation terms for Roe - #299

Merged
economon merged 1 commit into
developfrom
fix_roe
Aug 17, 2016
Merged

economon merged 1 commit into
developfrom
fix_roe

Conversation

@economon

Copy link
Copy Markdown
Member

Consistent calculation of dissipation terms for the Roe scheme for both explicit and implicit time integration. The previous explicit implementation was unstable with the high-order DG solver.

@economon
economon merged commit 84aa0c7 into develop Aug 17, 2016
@economon
economon deleted the fix_roe branch August 17, 2016 08:19
rois1995 added a commit that referenced this pull request Sep 29, 2026
… explicit general-gas Roe

In 3D the first three columns of the P matrix combine the entropy wave, weighted by n_j, with the shear waves.
Their strengths are n_j (Delta rho - Delta p / c^2) + (Delta u x n)_j; the entropy jump was used unweighted in the
first one only, so a contact or shear layer on a face not aligned with x got a wrong flux. The standard Roe scheme
already projects the jump with P^-1 since #299. Unit test roe_contact_tests.cpp: a contact and a shear layer give
the upwind flux in 2D and 3D.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01LEL91DW5WPbPwgFtCvHga6
rois1995 added a commit that referenced this pull request Sep 29, 2026
… explicit general-gas Roe

In 3D the first three columns of the P matrix combine the entropy wave, weighted by n_j, with the shear waves.
Their strengths are n_j (Delta rho - Delta p / c^2) + (Delta u x n)_j; the entropy jump was used unweighted in the
first one only, so a contact or shear layer on a face not aligned with x got a wrong flux. The standard Roe scheme
already projects the jump with P^-1 since #299. Unit test roe_contact_tests.cpp: a contact and a shear layer give
the upwind flux in 2D and 3D.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01LEL91DW5WPbPwgFtCvHga6
rois1995 added a commit that referenced this pull request Sep 29, 2026
… explicit general-gas Roe

In 3D the first three columns of the P matrix combine the entropy wave, weighted by n_j, with the shear waves.
Their strengths are n_j (Delta rho - Delta p / c^2) + (Delta u x n)_j; the entropy jump was used unweighted in the
first one only, so a contact or shear layer on a face not aligned with x got a wrong flux. The standard Roe scheme
already projects the jump with P^-1 since #299. Unit test roe_contact_tests.cpp: a contact and a shear layer give
the upwind flux in 2D and 3D.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01LEL91DW5WPbPwgFtCvHga6
rois1995 added a commit that referenced this pull request Sep 29, 2026
… explicit general-gas Roe

In 3D the first three columns of the P matrix combine the entropy wave, weighted by n_j, with the shear waves.
Their strengths are n_j (Delta rho - Delta p / c^2) + (Delta u x n)_j; the entropy jump was used unweighted in the
first one only, so a contact or shear layer on a face not aligned with x got a wrong flux. The standard Roe scheme
already projects the jump with P^-1 since #299. Unit test roe_contact_tests.cpp: a contact and a shear layer give
the upwind flux in 2D and 3D.

Co-Authored-By: Claude Opus 5.5 (1M context) <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01LEL91DW5WPbPwgFtCvHga6
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant