We focus on optimal control problems governed by partial differential equations. The presence of constraints on the state provides numerical and analytical difficulties. We treat these obstacles introducing a Lavrentiev regularization. The key issue is addressed to the analytical investigation of the convergence order of solutions based on certain source representation arguments. A virtual control approach is investigated for problems with additional control constraints. Numerical experiments are presented.