Different from long-lasting reconstruction processes which employ high order polynomials
enforced with some carefully designed limiting projections to seek stable solutions around
discontinuities, the current discretized scheme employs THINC (Tangent of Hyperbola for
INterface Capturing) functions with adaptive sharpness to solve both smooth and
discontinuous solutions. Due to the essentially monotone and bounded properties of THINC …