A numerical approach of friedrichs’ systems under constraints in bounded domains

Abstract

We present here an explicit finite volume scheme on unstructured meshes adapted to first-order hyperbolic systems under constraints in bounded domains. This scheme is based on the work (Coudière, Vila, Villedieu in C R Acad Sci Paris Sér I Math 331:95–100, 2000, [3]) in the unconstrained case and the splitting strategy of Després, Lagoutière, Seguin (Nonlinearity 24:3055–3081, 2011, [4]). We show that this scheme is stable under a Courant–Friedrichs–Lewy condition (and convergent for problems posed in the whole space), and we illustrate the solution constructed by this scheme on the example of the simplified model of perfect plasticity. From the theoretical point of view, the interaction between the constraint and the boundary of the domain in the model of perfect plasticity is encoded by a nonlinear boundary condition. With this numerical approach, we will show that, even if this scheme uses the underlying linear boundary condition, the results are consistent with the nonlinear model (and in particular with the nonlinear boundary condition).