Numerical method for the time optimal control problem governed by the Benjamin-Bona-Mahony equation