trixi-framework / trixi-framework/Trixi.jl

Positivity violation in AMR

Open
#2,064 1 comment 2 reactions 1 assignee Claimed by @Arpit-Babbar View on GitHub
enhancement
Dominant language
Julia
Stars
731
Forks
166
Avg merge
2d 18h
Merged PRs (30d)
25

Description

Zhang and Shu's scaling limiter uses admissibility of the means to enforce admissibility of the solution at the solution points. Quoting the [remark from our preprint](https://arxiv.org/pdf/2402.11926),

> the procedure (Zhang and Shu's scaling) doesn’t ensure that the polynomial is admissible at points which are not the solution points. Adaptive mesh refinement introduces such points into the numerical method and can thus cause a failure of admissibility preservation in the following situations: (a) mortar solution values obtained by interpolation are not admissible, (b) mean values of the solution values obtained by interpolating from the larger element as in are not admissible.

I have not run into this issue in `Trixi.jl`, but I did when I was using Trixi's positivity limiter and AMR callbacks in `TrixiLW.jl`. Then, in our work, we did the natural fix as follows

> Since the scaling limiter can be used to enforce admissibility of solution at any desired points, the remedy to both the issues is further scaling; we simply perform scaling of mortar solution values, interpolated solution point values with the admissible mean value of the larger element. This will ensure that the mortar solution point values and the mean values are admissible.

To fix this, we need to begin with https://github.com/trixi-framework/Trixi.jl/issues/215.

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.