In this study, we consider a stabilizer free weak Galerkin (SFWG) finite element method to solve a second‐order Sobolev equation. The SFWG method has various assets, including the support for higher order of accuracy and fewer coefficients. We have analyzed error estimates for both semidiscrete and fully discrete schemes using a non‐standard projection operator. The convergence rate of numerical error estimates of 0.3emO0.3em()hk+2+τ2$$ \kern0.3em O\kern0.3em \left({h}^{k+2}+{\tau}^2\right) $$ in L∞()H1$$ {L}^{\infty}\left({H}^1\right) $$ norm and 0.3emO0.3em()hk+3+τ2$$ \kern0.3em O\kern0.3em \left({h}^{k+3}+{\tau}^2\right) $$ in L∞()L2$$ {L}^{\infty}\left({L}^2\right) $$ norm are obtained, which are two order higher than the optimal order associated with WG finite element space ()ℙk(K),ℙk+1(∂K),[]ℙk+1(K)2.$$ \left({\mathrm{\mathbb{P}}}_k(K),{\mathrm{\mathbb{P}}}_{k+1}\left(\partial K\right),{\left[{\mathrm{\mathbb{P}}}_{k+1}(K)\right]}^2\right). $$ At last, several numerical examples are reported to validate the supercloseness property and efficiency of the proposed SFWG method. These experiments confirm the accuracy of the theoretical findings.