There is a significant bug in the face integration term on the right-hand-side of the gradient system for the BR1 implementation. Specifically, it neglects the contribution of the boundary faces, which implicitly assumes that the interior state satisfies the boundary conditions.
There is a significant bug in the face integration term on the right-hand-side of the gradient system for the BR1 implementation. Specifically, it neglects the contribution of the boundary faces, which implicitly assumes that the interior state satisfies the boundary conditions.