-
Notifications
You must be signed in to change notification settings - Fork 113
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
Mean value in Zhang-Shu limiter on curved meshes #1945
Mean value in Zhang-Shu limiter on curved meshes #1945
Conversation
…n mesh types other than TreeMesh
Review checklistThis checklist is meant to assist creators of PRs (to let them know what reviewers will typically look for) and reviewers (to guide them in a structured review process). Items do not need to be checked explicitly for a PR to be eligible for merging. Purpose and scope
Code quality
Documentation
Testing
Performance
Verification
Created with ❤️ by the Trixi.jl community. |
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #1945 +/- ##
==========================================
- Coverage 96.15% 96.15% -0.00%
==========================================
Files 454 454
Lines 36496 36509 +13
==========================================
+ Hits 35092 35103 +11
- Misses 1404 1406 +2
Flags with carried forward coverage won't be shown. Click here to find out more. ☔ View full report in Codecov by Sentry. |
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.
LGTM. Thanks
This is a bug fix where all mesh types other than
TreeMesh
should compute the solution mean valueu_mean
including the element Jacobians. If this is not done, it is possible to lose conservation. Typically, the values fromextra_analysis_errors = (:conservation_error,)
in theAnalysisCallback
should stay between1e-13
to1e-15
. However, when one computes the incorrect value foru_mean
on a curved mesh these conservation errors might increase to something between1e-4
to1e-8
(depending on the setup, how often the limiter is called, and how "warped" the mesh elements are).