  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.
