A concise overview is given of various numerical methods that can be used to analyse localization and failure in engineering materials. The importance of the cohesive-zone approach is emphasized and various ways to incorporate the cohesive-zone methodology in discretization methods are discussed. Numerical representations of cohesive-zone models suffer from a certain mesh bias. For discrete representations this is caused by the initial mesh design, while for smeared representations it is rooted in the ill-posedness of the rate boundary value problem that arises upon the introduction of decohesion. A proper representation of the discrete character of cohesive-zone formulations which avoids any mesh bias can be obtained elegantly when exploit...