A Monte Carlo Algorithm, capable of simulating numerically the time behavior of a spatially homogeneous reaction system, is expanded to include diffusion process. This expanded algorithm is applied to the Oregonator model of the Belousov-Zhabotinskii reaction to study the phenomenon of spatial oscillations. The results show that trigger waves travel at constant velocity which is uniquely determined by the given properties of the system and the velocity of trigger wave propagation increases linearly with the square root of both [H⁺] and [BrO₃^-].
The Monte Carlo algorithm for a homogeneous reaction system is also applied to the Oregonator model of the Belousov-Zhabotinskii reaction to study the phenomenon of temporal oscillations. The results show that the period of oscillation is 54 sec and the concentrations of intermediates exhibit the limit cycle behavior in phase plane. These results from the Monte Carlo algorithm are compared with the previous experimental results and with the simulation results obtained by a finite difference method.