Convergence of a numerical scheme for a mixed hyperbolic-parabolic system in two space dimensions (Q2806042)

The authors consider the numerical solution of the coupled system of partial differential equations NEWLINE\[NEWLINE\partial_t u + \mathrm{div} (f(u)v(w)) = (\alpha w - \beta) u,NEWLINE\]NEWLINE NEWLINE\[NEWLINE\partial_t w - \mu \Delta w = (\gamma - \delta u)w,\quad u(0,x,y) = u_0(x,y),\quad w(0,x,y) = w_0(x,y).NEWLINE\]NEWLINE For the discretization of the parabolic equation a finite-difference scheme is used and for the hyperbolic equation a Lax-Friedrichs type finite volume method. The basis of the convergence proof is an extension of Helly's theorem (see [\textit{C. M. Dafermos}, Hyperbolic conservation laws in continuum physics. 2nd ed. Berlin: Springer (2005; Zbl 1078.35001)]). For applying this theorem some properties and estimates for \(u\) and \(w\) are proved, namely the positivity of \(u\) and \(w\), \(L^1\) and \(L^\infty\) bounds on \(u\) and \(w\), a bound on the total variation in space, and the Lipschitz continuity in time of \(u\). Then it is proved that the approximations of \(u\) converge in \(L_{\mathrm{loc}}^1\) and the approximations of \(w\) weakly\(^\ast\) in \(L^\infty\). Finally, the convergence behaviour of the presented method is shown by numerical experiments.