Transport-equilibrium schemes for computing nonclassical shocks

This Note presents a very efficient numerical strategy for computing weak solutions of a scalar conservation law which fails to be genuinely nonlinear. In such a situation, the dynamics of shock solutions turns out to be mainly driven by a prescribed kinetic function that imposes the speed of propagation of the discontinuities. We show how to enforce the validity of the kinetic criterion at the discrete level. The resulting scheme provides, in addition, sharp profiles. Numerical evidence are included. To cite this article: C. Chalons, C. R. Acad. Sci. Paris, Ser. I 342 (2006).

RésuméCette Note présente un algorithme très efficace pour le calcul des solutions faibles d'une loi de conservation scalaire non vraiment nonlinéaire. Dans ce contexte, la dynamique des solutions choc repose principalement sur la donnée d'une fonction cinétique qui fixe la vitesse de propagation des discontinuités. Nous montrons comment forcer la validité du critère cinétique au niveau discret. Le schéma obtenu fournit par ailleurs des discontinuités sans diffusion numérique. Des résultats numériques sont présentés. Pour citer cet article : C. Chalons, C. R. Acad. Sci. Paris, Ser. I 342 (2006).