Carlos Parés

Manuel J. Castro, Philippe G. LeFloch, María Luz Muñoz-Ruiz et al.

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 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, a "convergence error" source-term which, provided the total variation of the approximations remains uniformly bounded, is a locally bounded measure. We 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 models arising in fluid dynamics. For systems having nonconservative products associated with linearly degenerate characteristic fields, the convergence error vanishes. For some other models, this measure is evaluated very accurately, especially by plotting the shock curves associated with each scheme under consideration.

Hugo Carrillo, Carlos Parés

A new family of high order methods for systems of conservation laws are introduced: the Compact Approximate Taylor (CAT) methods. These methods are based on centered (2p + 1)-point stencils where p is an arbitrary integer. We prove that the order of accuracy is 2p and that CAT methods are an extension of high-order Lax-Wendroff methods for linear problems. Due to this, they are linearly L2-stable under a CFL<1 condition. In order to prevent the spurious oscillations that appear close to discontinuities two shock-capturing techniques have been considered: a fux-limiter technique (FL-CAT methods) and WENO reconstruction for the frst time derivative (WENO-CAT methods). We follow the WENO-Lax Wendroff Approximate Taylor method of Zorio, Baeza and Mullet (2017) in the second approach. A number of test cases are considered to compare these methods with other WENO-based schemes: the linear transport equation, Burgers equation, and the 1D compressible Euler system are considered. Although CAT methods present an extra computational cost due to the local character, this extra cost is compensated by the fact that they still give good solutions with CFL values close to 1.