A numerical model was proposed to predict the hydrogen interstitial concentration in iron matrix both in the immediate surrounding of and far from an edge dislocation. The current formulation successfully incorporates the relaxation effect and, thus, correctly predicts the stress state around the dislocations.