We examine the Langevin diffusion confined to a closed, convex domain $D\subset\mathbb{R}^d$, represented as a reflected stochastic differential equation. We introduce a sequence of penalized stochastic differential equations and prove that their invariant measures converge, in Wasserstein-2 distance and with explicit polynomial rate, to the invariant measure of the reflected Langevin diffusion. We also analyze a time-discretization of the penalized process obtained via the Euler-Maruyama scheme and demonstrate the convergence to the original constrained measure. These results provide a rigorous approximation framework for reflected Langevin dynamics in both continuous and discrete time.