Abstract: Nonlinear propagation of shock waves through periodic structures have the potential to exhibit interesting phenomena. Frequency content of the shock that lies within a bandgap of the periodic structure is strongly attenuated, but nonlinear frequency-frequency interactions pumps energy back into those bands. To investigate the relative importance of these propagation phenomena, numerical experiments using the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation are carried out. Two-dimensional propagation through a periodic array of rectangular waveguides is performed by iteratively using the output of one waveguide as the input for the next waveguide. Comparison of the evolution of the initial shock wave for both the linear and nonlinear cases is presented.