Numerical Study of Phase Transition in An Exclusion Model with Parallel Dynamics