Summary: Numerical investigation of receptivity and flow transition in spatio-temporal framework have shown the central role of spatio-temporal wave-front (STWF) created by wall excitation for transition of a two-dimensional (2D) zero pressure gradient boundary layer (ZPGBL) in [the first author and S. Bhaumik, “Onset of turbulence from the receptivity stage of fluid flows”, Phys. Rev. Lett. 107, No. 15, Article ID 154501, 5 p. (2011; doi:10.1103/PhysRevLett.107.154501)]. Although the STWF is created by linear mechanism, it is the later nonlinear stage of evolution revealed by the solution of Navier-Stokes equation (NSE), which causes formation of turbulent spots merging together to create fully developed turbulent flow. Thus, computing STWF for ZPGBL from NSE is of prime importance, which has been reported by the present authors following earlier theoretical investigation. Similar computational efforts using NSE by other researchers do not report finding the STWF. In the present investigation we identify the main reason for other researchers to miss STWF, as due to taking a very short computational domain. Secondly, we show that even one takes a long enough domain and detect STWF, use of traditional low accuracy method will not produce the correct dynamics as reported by the first author and Bhaumik [loc. cit.]. The role of time integration plays a very strong role in the dynamics of transitional flows. We have shown here that implicit methods are more error prone, as compared to explicit time integration methods during flow transition. For the present problem, it is noted that the classical Crank-Nicolson method is unstable for 2D NSE. Same error-prone nature will also be noted for hybrid implicit-explicit time integration methods (known as the IMEX methods). One of the main feature of present analysis is to highlight the accuracy of computations by compact schemes used by the present investigators over a significantly longer domain and over unlimited time, as opposed to those reported earlier in the literature for the wall excitation problem. A consequence of taking long streamwise domain enables one to detect special properties of STWF and its nonlinear growth. The main focus of the present research is to highlight the importance of STWF, which is a new class of spatio-temporal solution obtained from the linear receptivity by solving Orr-Sommerfeld equation and nonlinear analysis of NSE.