It seems you have found a solution to this topic as it is closed, and the other post content is deleted.
If not, please refer to this example https://optics.ansys.com/hc/en-us/articles/360042800113-Group-delay-analysis
Once you get the function of phase vs wavelength you can get the group index/velocty by definition. This method works fine for many cases.