Currently, the gradient at the physical boundary (soln_grad_ext) set by boundary_face_values is not being enforced; as it is being set to the interior value (soln_grad_int) in all cases inside evaluate_auxiliary_flux.
To fix it, we must pass on_boundary flag to PhysicsBase's dissipative_flux and modify the dissipative flux accordingly; the fix take effect in assemble_boundary_term
Currently, the gradient at the physical boundary (
soln_grad_ext
) set byboundary_face_values
is not being enforced; as it is being set to the interior value (soln_grad_int
) in all cases insideevaluate_auxiliary_flux
.To fix it, we must pass
on_boundary
flag to PhysicsBase'sdissipative_flux
and modify the dissipative flux accordingly; the fix take effect inassemble_boundary_term