Sharp Convergence Rate of the Glimm Scheme for General Nonlinear Hyperbolic Systems

Consider a general strictly hyperbolic, quasilinear system, in one space dimension u(t) + A(u)u(x) = 0, (1) where u bar right arrow A(u), u is an element of Omega subset of R(N), is a smooth matrix-valued map. Given an initial datum u(0, .) with small total variation, let u(t, .) be the corresponding (unique) vanishing viscosity solution of (1) obtained as a limit of solutions to the viscous parabolic approximation u(t) + A(u)u(x) = mu u(xx), as mu --> 0. For every T >= 0, we prove the a-priori bound parallel to u(epsilon) (T, .) - u(T, .)parallel to(L1) = o(1) . root epsilon|log epsilon| (2) for an approximate solution u(epsilon) of (1) constructed by the Glimm scheme, with mesh size Delta x = Delta t = epsilon, and with a suitable choice of the sampling sequence. This result provides for general hyperbolic systems the same type of error estimates valid for Glimm approximate solutions of hyperbolic systems of conservation laws u(t) + F(u)(x) = 0 satisfying the classical Lax or Liu assumptions on the eigenvalues lambda(k) (u) and on the eigenvectors r(k) (u) of the Jacobian matrix A(u) = DF(u). The estimate (2) is obtained introducing a new wave interaction functional with a cubic term that controls the nonlinear coupling of waves of the same family and at the same time decreases at interactions by a quantity that is of the same order of the product of the wave strength times the change in the wave speeds. This is precisely the type of errors arising in a wave tracing analysis of the Glimm scheme, which is crucial to control in order to achieve an accurate estimate of the convergence rate as (2).