Enforcing bounds constraints through variational inequalities leads to accurate approximations in PDEs.