“…Since its introduction in the seminal paper by Bérenger [18], the Perfectly Matched Layer (PML) approach provides an alternative powerful tool to simulate the numerical solution of PDEs in unbounded domains [1,2,3,9,10,17,19,20,21,24,25,32,43]. Concerning nonlinear PDEs, the PML technique has been applied for example to the nonlinear Schrödinger equations [3,9,10,43], Euler and Navier-Stokes equations [31] or two-fluid plasma equations [40]. Concerning the nonlinear wave equations, Appelö and Kreiss [2] proposed and studied a first-order PML formulation with Dirichlet/Neumann boundary conditions related to Hagstrom techniques [1,29].…”