“…This radical approach has been shown to offer several advantages over conventional methods, such as accelerated convergence for steady state solution and higher order of accuracy for both primary and gradient variables, as demonstrated for diffusion [16], the incompressible/compressible Navier-Stokes equations [17,18], third-order dispersion equations [19], an incompressible magnetohydrodynamics model [20], an elliptic distance-function model [21], and so on. The original approach of Nishikawa has been further extended to a constant diffusion tensor by Lou et al [22,23] discretized with the reconstructed discontinuous Galerkin scheme (rDG). The hyperbolic approach has also been implemented for an anisotropic diffusion equation, based on highorder finite-volume schemes in Ref.…”