In this paper, we report a way to condition the global stiffness matrix resulting from the Galerkin edge element formulation, in order to make the BICG solver converge for high Peclet number situations. The concept of artificial diffusion is introduced and extended to edge element cases. A general expression of the artificial diffusion scheme for 3D edge elements is derived and a numerical validation in 2D case is presented. It was observed that although nearly identical results were achieved both with or without introducing artificial parameters for a 2D test problem, the former provided a better conditioned global stiffness matrix to ensure the BICG's convergence.