-
Notifications
You must be signed in to change notification settings - Fork 1
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
Covariant formulation of the shallow water equations #58
Conversation
…m/new_shallow_water
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #58 +/- ##
==========================================
+ Coverage 89.45% 90.39% +0.94%
==========================================
Files 20 22 +2
Lines 1801 2104 +303
==========================================
+ Hits 1611 1902 +291
- Misses 190 202 +12 ☔ 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.
Many thanks, @tristanmontoya! great work!
I have not checked the math thoroughly, but I leave a couple of comments below.
BTW, thanks for improving the documentation of my elixirs as well!
src/solvers/dgsem_p4est/containers_2d_manifold_in_3d_covariant.jl
Outdated
Show resolved
Hide resolved
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, @tristanmontoya!
Some minor comments below:
Co-authored-by: Andrés Rueda-Ramírez <aruedara@uni-koeln.de>
Great, thanks for your review! Going to make those changes and merge now. |
New version of #37 with entropy-stable split forms and auxiliary variables. Requires Trixi.jl v0.9.9 (see #59). The following new features are implemented:
CovariantShallowWaterEquations2D
which implements the split formulation of the covariant derivative using non-conservative terms, but does not use non-conservative terms for the standard formulation. The split formulation conserves total energy when no interface dissipation is added.func(u, aux_vars, equations)
rather than justfunc(u, equations)
, which is needed sinceentropy
is a function of auxiliary variables in this case.J
rather thansqrtG
for the area element.