Piecewise linear transformation in diffusive flux discretization

Abstract

To ensure the discrete maximum principle or solution positivity in finite volume schemes, diffusive flux is sometimes discretized as a conical combination of finite differences. Such a combination may be impossible to construct along material discontinuities using only cell concentration values. This is often resolved by introducing auxiliary node, edge, or face concentration values that are explicitly interpolated from the surrounding cell concentrations. We propose to discretize the diffusive flux after applying a local piecewise linear coordinate transformation that effectively removes the discontinuities. The resulting scheme does not need any auxiliary concentrations and is therefore remarkably simpler, while being second-order accurate under the assumption that the structure of the domain is locally layered.Comment: 11 pages, 1 figures, preprint submitted to Journal of Computational Physic