While the deleterious effects of hydrogen on metals and alloys are well known, the precise role of hydrogen in the underlying microscopic mechanisms is still not well understood and as of yet, the modeling attempts on hydrogen embrittlement and hydrogen induced cracking have not led to a proper method for life-time prediction. This work aims at the development of a robust numerical strategy in order to solve the non-linear coupled problem presented in the work of Anand [1]. The numerical implementation is performed for finite element method and the analysis are done to address the issue of hydrogen transport and hydrogen-embrittlement-related failures in metals. Specifically, problems related to the mechanism of hydrogen enhanced localized plasticity (HELP) is studied and macroscale shear localization phenomenon resulting from hydrogen induced material softening is considered at the phenomenological level. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.