MD simulations are used to study the effects of defects on the H diffusion in NbC. Buckingham potential is more accurate for diffusion of H atoms than LJ potential. H diffusion coefficient (D) increases with carbon vacancy concentration. H diffusion coefficient for 6 Å pore (radius = 6 Å) is as high as that for 20 Å pore. For small pores, H diffusion coefficient drops notably at elevated temperatures.