This paper deals with the non-conservative coupling of two one-dimensional barotropic Euler systems at an interface at x = 0. The closure pressure laws differ in the domains x < 0 and x > 0, and a Dirac source term concentrated at x = 0 models singular pressure losses. We propose two numerical methods. The first one relies on ghost state reconstructions at the interface while the second is based on a suitable relaxation framework. Both methods satisfy a well-balanced property for stationary solutions. In addition, the second method preserves mass conservation and exactly restores the prescribed singular pressure drops for both unsteady and steady solutions.
