Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes

• Published in: Journal of computational physics 2008-09, Vol.227 (17), p.8107-8129
• Main Authors: Castro, Manuel J., LeFloch, Philippe G., Muñoz-Ruiz, María Luz, Parés, Carlos
• Format: Article
• Language: English
• Subjects: Computational techniques; Convergence error measure; Equivalent equation; Exact sciences and technology; Family of paths; Formally path-consistent scheme; Mathematical methods in physics; Nonconservative hyperbolic system; Physics; Shock wave
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2008.05.012

We are interested in nonlinear hyperbolic systems in nonconservative form arising in fluid dynamics, and, for solutions containing shock waves, we investigate the convergence of finite difference schemes applied to such systems. According to Dal Maso, LeFloch, and Murat’s theory, a shock wave theory for a given nonconservative system requires prescribing a priori a family of paths in the phase space. In the present paper, we consider schemes that are formally consistent with a given family of paths, and we investigate their limiting behavior as the mesh is refined. we first generalize to systems a property established earlier by Hou and LeFloch for scalar conservation laws, and we prove that nonconservative schemes generate, at the level of the limiting hyperbolic system, an convergence error source-term which, provided the total variation of the approximations remains uniformly bounded, is a locally bounded measure. This convergence error measure is supported on the shock trajectories and, as we demonstrate here, is usually “small”. In the special case that the scheme converges in the sense of graphs – a rather strong convergence property often violated in practice – then this measure source-term vanishes. We also discuss the role of the equivalent equation associated with a difference scheme; here, the distinction between scalar equations and systems appears most clearly since, for systems, the equivalent equation of a scheme that is formally path-consistent depends upon the prescribed family of paths. The core of this paper is devoted to investigate numerically the approximation of several (simplified or full) hyperbolic models arising in fluid dynamics. This leads us to the conclusion that for systems having nonconservative products associated with linearly degenerate characteristic fields, the convergence error vanishes. For more general models, this measure is evaluated very accurately, especially by plotting the shock curves associated with each scheme under consideration; as we demonstrate, plotting the shock curves provide a convenient approach for evaluating the range of validity of a given scheme.