Upwind difference approximations for degenerate parabolic convection-diffusion equations with a discontinuous coefficient

We analyse approximate solutions generated by an upwind differencescheme (of Engquist–Osher type) for nonlinear degenerateparabolic convection–diffusion equations where the nonlinearconvective flux function has a discontinuous coefficient (x)and the diffusion function A(u) is allowed to be strongly degenerate(the pure hyperbolic case is included in our setup). The mainproblem is obtaining a uniform bound on the total variationof the difference approximation u, which is a manifestationof resonance. To circumvent this analytical problem, we constructa singular mapping (, ·) such that the total variationof the transformed variable z = (, u) can be bounded uniformlyin . This establishes strong L¹ compactness of z and, since(, ·) is invertible, also u. Our singular mapping isnovel in that it incorporates a contribution from the diffusionfunction A(u). We then show that the limit of a converging sequenceof difference approximations is a weak solution as well as satisfyinga Krukov-type entropy inequality. We prove that the diffusionfunction A(u) is Hölder continuous, implying that the constructedweak solution u is continuous in those regions where the diffusionis nondegenerate. Finally, some numerical experiments are presentedand discussed.