The BZ reaction exhibits chemical oscillations and spatial pattern formation. The Oregonator model: du/dt = (u − u² − fv(u−q)/(u+q))/ε + D_u∇²u, dv/dt = u − v. Spiral waves and target patterns emerge from the interplay of autocatalytic excitation and inhibition, diffusion-coupled across the 2D grid.