Numerical method for the bone regeneration model, defined within the evolving 2D axisymmetric physical domain