View article

Numerical discretizations of the Generalized Porous Medium Equation (GPME) with discontinuous coefficients are analyzed with respect to the formation of numerical artifacts. In addition to the degeneracy and self-sharpening of the GPME with continuous coefficients, detailed in [1], increased numerical challenges occur in the discontinuous coefficients case. These numerical challenges manifest themselves in spurious temporal oscillations in second order finite volume discretizations with both arithmetic and harmonic averaging. The integral average, developed in [2], leads to improved solutions with monotone and reduced amplitude temporal oscillations. In this paper, we propose a new method called the Shock-Based Averaging Method (SAM) that incorporates the shock position into the numerical scheme. The shock position is numerically calculated by discretizing the theoretical speed of the front from the …