We investigate the dynamics of a two dimensional bi-axial next nearest neighbour Ising (BNNNI) model following a quench to zero temperature. The Hamiltonian is given by H = −J0 L i,j=1 [(Si,j Si+1,j + Si,jSi,j+1) − κ(Si,j Si+2,j + Si,jSi,j+2)] . For κ < 1, the system does not reach the equilibrium ground state and keep evolving in active states for ever. For κ ≥ 1, though the system reaches a final state, but it do not reach the ground state always and freezes to a striped state with a finite probability like two dimensional ferromagnetic Ising model and ANNNI model. The overall dynamical behaviour for κ > 1 and κ = 1 is quite different. The residual energy decays in a power law for both κ > 1 and κ = 1 from which the dynamical exponent z have been estimated. The persistence probability shows algebraic decay for κ > 1 with an exponent θ = 0.22 ± 0.002 while the dynamical exponent for ordering z = 2.33 ± 0.01. For κ = 1, the system belongs to a completely different dynamical class with θ = 0.332 ± 0.002 and z = 2.47 ± 0.04. We have computed the freezing probability for different values of κ. We have also studied the decay of autocorrelation function with time for different regime of κ values. The results have been compared with that of the two dimensional ANNNI model.