“…Then, a high order staggered DG method for general secondorder elliptic problems is developed in [51], and it is applied to various physical problems arising from practical applications [55,51,36,52,47,50,54,48]. For a pressure-robust method, Zhao et al [53] proposed a lowest-order pressure-robust staggered DG method for the Stokes equations. They used a reconstruction operator based on an H(div) conforming function space on polygonal meshes.…”