Hydrodynamic limit for a nongradient interacting particle system with stochastic reservoirs

We consider a nongradient interacting particle system whose macroscopic behavior is described by a $d$-dimensional nonlinear parabolic equation on a square with boundary conditions. Assuming that the diffusion coefficient is Lipschitz, we prove that the rescaled density field converges to a unique weak solution of the parabolic equation.