This paper aims to construct a sixth-order weighted essentially nonoscillatory scheme for simulating the nonlinear degenerate parabolic equations in a finite difference framework.
To design this scheme, we approximate the second derivative in these equations in a different way, which of course is still in a conservative form. In this way, unlike the common practice of reconstruction, the approximation of the derivatives of odd order is needed to develop the numerical flux.
The results obtained by the new scheme produce less error compared to the results of other schemes in the literature that are recently developed for the nonlinear degenerate parabolic equations while requiring less computational times.
This research develops a new weighted essentially nonoscillatory scheme for solving the nonlinear degenerate parabolic equations in multidimensional space. Besides, any selection of the constants (sum equals one is the only requirement for them), named the linear weights, will obtain the desired accuracy.
