We propose a weak Galerkin finite element method for the Westervelt model of ultrasound
waves on polygonal meshes. More precisely, we study the spatial discretization of …