Numerical Preservation of Velocity Induced Invariant Regions for Reaction–Diffusion Systems on Evolving Surfaces